We discuss the production of a class of heavy sterile neutrinos $\nu_h$ in proto-neutron stars. The neutrinos, of mass around $50$ MeV, have a negligible mixing with the active species but relatively large dimension-5 electromagnetic couplings. In particular, a magnetic dipole moment $\mu\approx 10^{-6}$ GeV$^{-1}$ implies that they are thermally produced through $e^+ e^-\to \bar \nu_h \nu_h$ in the early phase of the core collapse, whereas a heavy--light transition moment $\mu_{\rm tr}\approx 10^{-8}$ GeV$^{-1}$ allows their decay $\nu_h\to \nu_i \gamma$ with a lifetime around $10^{-3}$ s. This type of electromagnetic couplings has been recently proposed to explain the excess of electron-like events in baseline experiments. We show that the production and decay of these heavy neutrinos would transport energy from the central regions of the star to distances $d\approx 400$ km, providing a very efficient mechanism to enhance the supernova shock front and heat the material behind it.
We present molecular line imaging observations of three massive molecular outflow sources, G333.6-0.2, G333.1-0.4, and G332.8-0.5, all of which also show evidence for infall, within the G333 giant molecular cloud (GMC). All three are within a beam size (36 arcseconds) of IRAS sources, 1.2-mm dust clumps, various masing species and radio continuum-detected HII regions and hence are associated with high-mass star formation. We present the molecular line data and derive the physical properties of the outflows including the mass, kinematics, and energetics and discuss the inferred characteristics of their driving sources. Outflow masses are of 10 to 40 solar masses in each lobe, with core masses of order 10^3 solar masses. outflow size scales are a few tenth of a parsec, timescales are of several x10^4 years, mass loss rates a few x10^-4 solar masses/year. We also find the cores are turbulent and highly supersonic.
Strong lensing time-delay systems constrain cosmological parameters via the so-called time-delay distance and the angular diameter distance to the lens. In previous studies, only the former information was used. In this paper, we show that the cosmological constraints improve significantly when the latter information is also included. Specifically, the angular diameter distance plays a crucial role in breaking the degeneracy between the curvature of the Universe and the time-varying equation of state of dark energy. Using a mock sample of 55 bright quadruple lens systems based on expectations for ongoing/future imaging surveys, we find that adding the angular diameter distance information to the time-delay distance information and the cosmic microwave background data of Planck improves the constraint on the constant equation of state by 30%, on the time variation in the equation of state by a factor of two, and on the Hubble constant in the flat $\Lambda$CDM model by a factor of two. Therefore, previous forecasts for the statistical power of time-delay systems were significantly underestimated, i.e., time-delay systems are more powerful than previously appreciated.
With the first phase of the Square Kilometre Array (SKA1) entering into its final pre-construction phase, we investigate how best to maximise its scientific return. Specifically, we focus on the statistical measurement of the 21 cm power spectrum (PS) from the epoch of reionization (EoR) using the low frequency array, SKA1-low. To facilitate this investigation we use the recently developed MCMC based EoR analysis tool 21CMMC (Greig & Mesinger). In light of the recent 50 per cent cost reduction, we consider several different SKA core baseline designs, changing: (i) the number of antenna stations; (ii) the number of dipoles per station; and also (iii) the distribution of baseline lengths. We find that a design with a reduced number of dipoles per core station (increased field of view and total number of core stations), together with shortened baselines, maximises the recovered EoR signal. With this optimal baseline design, we investigate three observing strategies, analysing the trade-off between lowering the instrumental thermal noise against increasing the field of view. SKA1-low intends to perform a three tiered observing approach, including a deep 100 deg$^{2}$ at 1000 h, a medium-deep 1000 deg$^{2}$ at 100 h and a shallow 10,000 deg$^{2}$ at 10 h survey. We find that the three observing strategies result in comparable ($\lesssim$ per cent) constraints on our EoR astrophysical parameters. This is contrary to naive predictions based purely on the total signal-to-noise, thus highlighting the need to use EoR parameter constraints as a figure of merit, in order to maximise scientific returns with next generation interferometers.
We use four years of Kepler photometry to study the non-eclipsing spectroscopic binary KIC 10080943. We find both components to be $\gamma$ Doradus/$\delta$ Scuti hybrids, which pulsate in both p and g modes. We present an analysis of the g modes, which is complicated by the fact that the two sets of $\ell=1$ modes partially overlap in the frequency spectrum. Nevertheless, it is possible to disentangle them by identifying rotationally split doublets from one component and triplets from the other. The identification is helped by the presence of additive combination frequencies in the spectrum that involve the doublets but not the triplets. The rotational splittings of the multiplets imply core rotation periods of about 11 d and 7 d in the two stars. One of the stars also shows evidence of $\ell=2$ modes.
Calibrating the photometric redshifts of >10^9 galaxies for upcoming weak lensing cosmology experiments is a major challenge for the astrophysics community. The path to obtaining the required spectroscopic redshifts for training and calibration is daunting, given the anticipated depths of the surveys and the difficulty in obtaining secure redshifts for some faint galaxy populations. Here we present an analysis of the problem based on the self-organizing map, a method of mapping the distribution of data in a high-dimensional space and projecting it onto a lower-dimensional representation. We apply this method to existing photometric data from the COSMOS survey selected to approximate the anticipated Euclid weak lensing sample, enabling us to robustly map the empirical distribution of galaxies in the multidimensional color space defined by the expected Euclid filters. Mapping this multicolor distribution lets us determine where - in galaxy color space - redshifts from current spectroscopic surveys exist and where they are systematically missing. Crucially, the method lets us determine whether a spectroscopic training sample is representative of the full photometric space occupied by the galaxies in a survey. We explore optimal sampling techniques and estimate the additional spectroscopy needed to map out the color-redshift relation, finding that sampling the galaxy distribution in color space in a systematic way can efficiently meet the calibration requirements. While the analysis presented here focuses on the Euclid survey, similar analysis can be applied to other surveys facing the same calibration challenge, such as DES, LSST, and WFIRST.
The ultimate goal of the observation of nonthermal emission from astrophysical sources is to understand the underlying particle acceleration and evolution processes, and few tools are publicly available to infer the particle distribution properties from the observed photon spectra from X-ray to VHE gamma rays. Here I present naima, an open source Python package that provides models for nonthermal radiative emission from homogeneous distribution of relativistic electrons and protons. Contributions from synchrotron, inverse Compton, nonthermal bremsstrahlung, and neutral-pion decay can be computed for a series of functional shapes of the particle energy distributions, with the possibility of using user-defined particle distribution functions. In addition, naima provides a set of functions that allow to use these models to fit observed nonthermal spectra through an MCMC procedure, obtaining probability distribution functions for the particle distribution parameters. Here I present the models and methods available in naima and an example of their application to the understanding of a galactic nonthermal source. naima's documentation, including how to install the package, is available at this http URL
We show that hydrodynamic turbulent cloud simulations naturally produce large filaments made up of a network of smaller and coherent sub-filaments. Such simulations resemble observations of filaments and fibres in nearby molecular clouds. The sub-filaments are dynamical features formed at the stagnation points of the turbulent velocity field where shocks dissipate the turbulent energy. They are a ubiquitous feature of the simulated clouds, which appear from the beginning of the simulation and are not formed by gradual fragmentation of larger filaments. Most of the sub-filaments are gravitationally sub-critical and do not fragment into cores, however, there is also a significant fraction of supercritical sub-filaments which break up into star-forming cores. The sub-filaments are coherent along their length, and the residual velocities along their spine show that they are subsonically contracting without any ordered rotation on scales of ~0.1 pc. Accretion flows along the sub-filaments can feed material into star forming cores embedded within the network. The overall mass in sub-filaments and the number of sub-filaments increases as the cloud evolves. We propose that the formation of filaments and sub-filaments is a natural consequence of the turbulent cascade in the complex multi-phase interstellar medium. Sub-filaments are formed by the high wavenumber, small scale modes in the turbulent velocity field. These are then stretched by local shear motions and gathered together by a combination of low wavenumber modes and gravitational contraction on larger scales, and by doing so build up the extended filaments observed in column density maps.
We present two NuSTAR observations of the local Seyfert 2 active galactic nucleus (AGN) and an ultraluminous X-ray source (ULX) candidate in NGC 5643. Together with archival data from Chandra, XMM-Newton and Swift-BAT, we perform a high-quality broadband spectral analysis of the AGN over two decades in energy ($\sim$0.5-100 keV). Previous X-ray observations suggested that the AGN is obscured by a Compton-thick (CT) column of obscuring gas along our line-of-sight. However, the lack of high-quality $\gtrsim$ 10 keV observations, together with the presence of a nearby X-ray luminous source, NGC 5643 X-1, had left significant uncertainties in the characterization of the nuclear spectrum. NuSTAR now enables the AGN and NGC 5643 X-1 to be separately resolved above 10 keV for the first time and allows a direct measurement of the absorbing column density toward the nucleus. The new data show that the nucleus is indeed obscured by a CT column of $N_{\rm{H}}$ $\gtrsim$ 5 $\times$ 10$^{24}$ cm$^{-2}$. The range of 2-10 keV absorption-corrected luminosity inferred from the best fitting models is $L_{2-10,\rm{int}} =$ (0.8-1.7) $\times$ 10$^{42}$ erg s$^{-1}$, consistent with that predicted from multiwavelength intrinsic luminosity indicators. We also study the NuSTAR data for NGC 5643 X-1, and show that it exhibits evidence for a spectral cut-off at energy, $E$ $\sim$ 10 keV, similar to that seen in other ULXs observed by NuSTAR. Along with the evidence for significant X-ray luminosity variations in the 3-8 keV band from 2003-2014, our results further strengthen the ULX classification of NGC 5643 X-1.
The usual treatments of the neutrino flavor-evolution, beyond a surface above the last scattering, assume identical angular distributions at this distance for the different initial (unmixed) flavors, and for particles and antiparticles. Taking into account differences in these distributions that must be present, as a result of the species-dependent scattering cross-sections lower down in the star, leads to a new set of non-linear equations. These equations are unstable, even at the initial surface, with respect to perturbations that have both tiny mixing from neutrino oscillations and that break all-over spherical-symmetry. There could be important consequences for: 1) the dynamics of the explosion, 2) R-process-nucleosynthesis, and 3) some future observable neutrino pulse.
The seek of high precision analyses in $\gamma$-ray astronomy leads to the implementation of multivariate combination, benefiting from several reconstruction methods. Such analysis, called $X_{eff}$, was developed for the H.E.S.S. data using three shower reconstruction methods. This paper presents the improvement granted to this analysis by refining the distribution calculation of discriminant variables, considering observation conditions, and adding new variables in the $X_{eff}$ combination. The efficiency of the analysis is presented using simulations and real data. A comparison with the standard analysis (model++), for a typical set of sources, shows a significant gain in sensitivity.
We consider the prospects for the indirect detection of low mass dark matter which couples dominantly to quarks. If the center of mass energy is below about 280 MeV, the kinematically allowed final states will be dominated by photons and neutral pions, producing striking signatures at gamma ray telescopes. In fact, an array of new instruments have been proposed, which would greatly improve sensitivity to photons in this energy range. We find that planned instruments can improve on current sensitivity to dark matter models of this type by up to a few orders of magnitude.
We present radio continuum light curves of the magnetar SGR J1745$-$2900 and Sgr A* obtained with multi-frequency, multi-epoch Very Large Array observations between 2012 and 2014. During this period, a powerful X-ray outburst from SGR J1745$-$2900 occurred on 2013-04-24. Enhanced radio emission is delayed with respect to the X-ray peak by about seven months. In addition, the flux density of the emission from the magnetar fluctuates by a factor of 2 to 4 at frequencies between 21 and 41 GHz and its spectral index varies erratically. Here we argue that the excess fluctuating emission from the magnetar arises from the interaction of a shock generated from the X-ray outburst with the orbiting ionized gas at the Galactic center. In this picture, variable synchrotron emission is produced by ram pressure variations due to inhomogeneities in the dense ionized medium of the Sgr A West bar. The pulsar with its high transverse velocity is moving through a highly blue-shifted ionized medium. This implies that the magnetar is at a projected distance of $\sim0.1$ pc from Sgr A* and that the orbiting ionized gas is partially or largely responsible for a large rotation measure detected toward the magnetar. Despite the variability of Sgr A* expected to be induced by the passage of the G2 cloud, monitoring data shows a constant flux density and spectral index during this period
We present our results of a combined analysis of radial velocity and light curves of a double-lined spectroscopic and eclipsing binary KIC 09246715, observed photometrically by the $Kepler$ satellite, and spectroscopically with the OAO-1.88m telescope with the HIgh-Dispertion Echelle Spectrograph (HIDES). The target was claimed to be composed of two red giants, one of which is showing solar-like oscillations. We have found that the mass and radius of the primary are $M_1=2.169\pm0.024$ M$_\odot$ and $R_1=8.47\pm0.13$ R$_\odot$, and of the secondary: $M_2=2.143\pm0.025$ M$_\odot$ and $R_2=8.18\pm0.09$ R$_\odot$, which confirms the double-giant status. Our secondary is the star to which the oscillations were attributed. Results of its previous asteroseismic analysis are in agreement with ours, only significantly less precise, but the subsequent light-curve-based study failed to derive correct mass and radius of our primary. KIC 09246715 is one of the rare cases where asteroseismic parameters of a solar-like oscillator were confirmed by an independent method, and only the third example of a Galactic double-giant eclipsing binary with masses and radii measured with precision below 2%.
We present a data-adaptive spectral method - Monte Carlo Singular Spectrum Analysis (MC-SSA) - and its modification to tackle astrophysical problems. Through numerical simulations we show the ability of the MC-SSA in dealing with $1/f^{\beta}$ power-law noise affected by photon counting statistics. Such noise process is simulated by a first-order autoregressive, AR(1) process corrupted by intrinsic Poisson noise. In doing so, we statistically estimate a basic stochastic variation of the source and the corresponding fluctuations due to the quantum nature of light. In addition, MC-SSA test retains its effectiveness even when a significant percentage of the signal falls below a certain level of detection, e.g., caused by the instrument sensitivity. The parsimonious approach presented here may be broadly applied, from the search for extrasolar planets to the extraction of low-intensity coherent phenomena probably hidden in high energy transients.
We present a dynamical classification system for galaxies based on the shapes of their circular velocity curves (CVCs). We derive the CVCs of 40 SAURON and 42 CALIFA galaxies across Hubble sequence via a full line-of-sight integration as provided by solutions of the axisymmetric Jeans equations. We use Principal Component Analysis (PCA) applied to the circular curve shapes to find characteristic features and use a k-means classifier to separate circular curves into classes. This objective classification method identifies four different classes, which we name Slow-Rising (SR), Flat (F), Sharp-Peaked (SP) and Round-Peaked (RP) circular curves. SR-CVCs are mostly represented by late-type spiral galaxies (Scd-Sd) with no prominent spheroids in the central parts and slowly rising velocities; F-CVCs span almost all morphological types (E,S0,Sab,Sb-Sbc) with flat velocity profiles at almost all radii; SP-CVCs are represented by early-type and early-type spiral galaxies (E,S0,Sb-Sbc) with prominent spheroids and sharp peaks in the central velocities. RP-CVCs are represented by only two morphological types (E,Sa-Sab) with prominent spheroids, but RP-CVCs have much rounder peaks in the central velocities than SP-CVCs. RP-CVCs are typical for high-mass galaxies, while SR-CVCs are found for low-mass galaxies. Intermediate-mass galaxies usually have F-CVCs and SP-CVCs. Circular curve classification presents an alternative to typical morphological classification and may be more tightly linked to galaxy evolution.
The near-Earth asteroid (196256) 2003 EH1 has been suggested to have a dynamical association with the Quadrantid meteoroid stream. We present photometric observations taken to investigate the physical character of this body and to explore its possible relation to the stream. We find no evidence for on-going mass-loss. A model fitted to the point-like surface brightness profile at 2.1 AU limits the fractional contribution to the integrated brightness by near-nucleus coma to $\leq$ 2.5 %. Assuming an albedo equal to those typical of cometary nuclei ($\it p_{\rm R}$=0.04), we find that the effective nucleus radius is $r_e$ = 2.0$\pm$0.2 km. Time-resolved ${\it R}$-band photometry can be fitted by a two-peaked lightcurve having a rotational period of 12.650$\pm$0.033 hr. The range of the lightcurve, $\Delta m_{\rm R}$= 0.44 $\pm$ 0 .01 mag, is indicative of an elongated shape having an axis ratio $\sim$1.5 projected into the plane of the sky. The asteroid shows colors slightly redder than the Sun, being comparable with those of C-type asteroids. The limit to the mass loss rate set by the absence of resolved coma is $\lesssim$ 2.5$\times$ 10$^{-2}$ kg ${\rm s^{-1}}$, corresponding to an upper limit on the fraction of the surface that could be sublimating water ice $f_A$ $\lesssim$ 10$^{-4}$. Even if sustained over the 200-500 yr dynamical age of the Quadrantid stream, the total mass loss from 2003 EH1 would be too small to supply the reported stream mass ($10^{13}$ kg), implying either that the stream has another parent or that mass loss from 2003 EH1 is episodic.
The CoRoT space mission was operating for almost 6 years, producing thousands of continuous photometric light curves. The temporal series of exposures are processed by the production pipeline, correcting the data for known instrumental effects. But even after these model-based corrections, some collective trends are still visible in the light curves. We propose here a simple exposure-based algorithm to remove instrumental effects. The effect of each exposure is a function of only two instrumental stellar parameters, position on the CCD and photometric aperture. The effect is not a function of the stellar flux, and therefore much more robust. As an example, we show that the $\sim2\%$ long-term variation of the early run LRc01 is nicely detrended on average. This systematics removal process is part of the CoRoT legacy data pipeline.
The purpose of the hardware presented in this contribution is to decrease the energy threshold of the MAGIC telescopes without significantly increasing the data acquisition rate. To achieve this purpose, we developed an additional level of trigger that relies on the location in both MAGIC cameras where the trigger is issued to rule out accidental events. This allows to decrease the Discriminator Threshold (DT), which results in a reduction of the energy threshold of the instrument. We simulated the Topo-trigger concept using the standard MAGIC Monte Carlo (MC) and tested it with real telescope data. In this paper we show the concept and results of these tests.
We present the status of the laue project devoted to develop a technology for building a 20 meter long focal length Laue lens for hard x-/soft gamma-ray astronomy (80 - 600 keV). The Laue lens is composed of bent crystals of Gallium Arsenide (GaAs, 220) and Germanium (Ge, 111), and, for the first time, the focusing property of bent crystals has been exploited for this field of applications. We show the preliminary results concerning the adhesive employed to fix the crystal tiles over the lens support, the positioning accuracy obtained and possible further improvements. The Laue lens petal that will be completed in a few months has a pass band of 80 - 300 keV and is a fraction of an entire Laue lens capable of focusing X-rays up to 600 keV, possibly extendable down to 20 - 30 keV with suitable low absorption crystal materials and focal length. The final goal is to develop a focusing optics that can improve the sensitivity over current telescopes in this energy band by 2 orders of magnitude.
The Galactic Centre has been studied with the High Energy Stereoscopic System (H.E.S.S.) for over 10 years, revealing a bright, complex gamma-ray morphology. Besides a strong point-like very-high-energy gamma-ray source coincident with the supermassive black hole Sgr A*, pre- vious analyses also revealed a diffuse ridge of gamma-ray emission, indicative of a powerful cosmic-ray accelerator in this region. The addition of a fifth telescope with 600 m 2 mirror area to the centre of the H.E.S.S. array has increased the energy range accessible, allowing observations to take place down to 100 GeV and potentially below. This wider energy range allows an important overlap in observations with satellite instruments such as the Fermi-LAT gamma-ray telescope. We will present the results of new H.E.S.S observations of the Galactic Centre region and show a detailed analysis of the central source, including comparisons to results at other wavelengths.
Molecular clouds, the birthplaces of stars in galaxies, form dynamically from the diffuse atomic gas of the interstellar medium (ISM). The ISM is also threaded by magnetic fields which have a large impact on its dynamics. In particular, star forming regions must be magnetically supercrit- ical in order to accomodate gas clumps which can collapse under their own weight. Based on a parameter study of three dimensional magneto-hydrodyamical (MHD) simulations, we show that the long-standing problem of how such supercritical regions are generated is still an open issue.
In~\citep{Shafieloo2009}, Shafieloo, Sanhi and Starobinsky firstly proposed the possibility that the current cosmic acceleration (CA) is slowing down. This is rather counterintuitive, because a slowing down CA cannot be accommodated in almost all the mainstream cosmological models. In this work, by exploring the evolutionary trajectories of dark energy equation of state $w(z)$ and deceleration parameter $q(z)$, we present a comprehensive investigation on the slowing down of CA from both the theoretical and the observational sides. For the theoretical side, we study the impacts of different $w(z)$ by using six parametrization models, and then discuss the effects of spatial curvature. For the observational side, we investigate the effects of different type Ia supernovae (SNe Ia), different baryon acoustic oscillation (BAO), and different cosmic microwave background (CMB) data, respectively. We find that the evolution of CA are insensitive to the specific form of $w(z)$; in contrast, a non-flat Universe more favors a slowing down CA than a flat Universe. Moreover, we find that SNLS3 SNe Ia datasets favor a slowing down CA at 1$\sigma$ confidence level, while JLA SNe Ia samples prefer an eternal CA; in contrast, the effects of different BAO data are negligible. In addition, compared with CMB distance prior data, full CMB data more favor a slowing down CA. Since the evolutionary behavior of CA depends on both the theoretical models and the observational data, the possibility of slowing down CA cannot be confirmed by current observations.
The 21cm-galaxy cross-power spectrum is expected to be one of the promising probes of the Epoch of Reionization (EoR), as it could offer information about the progress of reionization and the typical scale of ionized regions at different redshifts. With upcoming observations of 21cm emission from the EoR with the Low Frequency Array (LOFAR), and of high redshift Lyalpha emitters (LAEs) with Subaru's Hyper Suprime Cam (HSC), we investigate the observability of such cross-power spectrum with these two instruments, which are both planning to observe the ELAIS-N1 field at z=6.6. In this paper we use N-body + radiative transfer (both for continuum and Lyalpha photons) simulations at redshift 6.68, 7.06 and 7.3 to compute the 3D theoretical 21cm-galaxy cross-power spectrum, as well as to predict the 2D 21cm-galaxy cross-power spectrum expected to be observed by LOFAR and HSC. Once noise and projection effects are accounted for, our predictions of the 21cm-galaxy cross-power spectrum show clear anti-correlation on scales larger than ~ 60 h$^{-1}$ Mpc (corresponding to k ~ 0.1 h Mpc$^{-1}$), with levels of significance p=0.04 at z=6.6 and p=0.048 at z=7.3. On smaller scales, instead, the signal is completely contaminated.
A number of hard X-Ray (10 - 100 KeV) astronomical missions of near future will make use of multilayer-coated focusing mirrors. The technology based on Nickel electroformed replication is suitable for the multilayer optics realization, since multi-modular telescopes are foreseen. For example, for the Constellation-X mission there is the need of realizing up to 14 identical modules (12 flight modules plus two spares) which can be replicated by the same series of mandrels. The Ni replication approach is derived from the method already successfully used for making the Au coated soft X-ray mirrors with good imaging performances of the missions BeppoSAX, XMM-Newton and Swift. In the technological extension of the process to the multilayer optics fabrication, it would be convenient to overcoat the external surface of mandrels (normally in Kanigen) with a layer made of a very hard material. This would help to maintain the very low roughness level requested by the application (typically less than a couple of Angstroms for a 1 micrometer scan length with AFM) also after many replications and successive cleaning of the mandrel. Good material candidate are at this regard TiN and SiC, both characterized by a very high hardness. We have proven that flat prototypes with TiN and SiC overcoating can be superpolished at a level comparable to the traditional electroless Nickel coating. In this paper we will present a characterization by topographic measurement (AFM and WYKO) and by X-Ray scattering of two of these samples.
When very high-energy photons (VHE, E>100 GeV) propagate over cosmological distances, they interact with background light by pair production. Observations of spectral features in the VHE band of extragalactic sources related to this energy-dependent absorption process with the H.E.S.S. array of Cherenkov telescopes allow measuring the spectral energy distribution (SED) of the extragalactic background light (EBL), otherwise very difficult to determine. Preliminary results on the determination of the SED of the EBL will be presented, based on the measurements of the energy spectra of blazars with H.E.S.S. at redshifts up to z = 0.2. This model independent approach shows that the shape and overall normalization of the EBL SED is accessible.
The realization of X-ray telescopes with imaging capabilities in the hard (> 10 keV) X-ray band requires the adoption of optics with shallow (< 0.25 deg) grazing angles to enhance the reflectivity of reflective coatings. On the other hand, to obtain large collecting area, large mirror diameters (< 350 mm) are necessary. This implies that mirrors with focal lengths >10 m shall be produced and tested. Full-illumination tests of such mirrors are usually performed with on- ground X-ray facilities, aimed at measuring their effective area and the angular resolution; however, they in general suffer from effects of the finite distance of the X-ray source, e.g. a loss of effective area for double reflection. These effects increase with the focal length of the mirror under test; hence a "partial" full-illumination measurement might not be fully representative of the in-flight performances. Indeed, a pencil beam test can be adopted to overcome this shortcoming, because a sector at a time is exposed to the X-ray flux, and the compensation of the beam divergence is achieved by tilting the optic. In this work we present the result of a hard X-ray test campaign performed at the BL20B2 beamline of the SPring-8 synchrotron radiation facility, aimed at characterizing the Point Spread Function (PSF) of a multilayer-coated Wolter-I mirror shell manufactured by Nickel electroforming. The mirror shell is a demonstrator for the NHXM hard X-ray imaging telescope (0.3 - 80 keV), with a predicted HEW (Half Energy Width) close to 20 arcsec. We show some reconstructed PSFs at monochromatic X-ray energies of 15 to 63 keV, and compare them with the PSFs computed from post-campaign metrology data, self-consistently treating profile and roughness data by means of a method based on the Fresnel diffraction theory. The modeling matches the measured PSFs accurately.
Stellar limb darkening impacts a wide range of astronomical measurements. The accuracy to which it is modeled limits the accuracy in any covariant parameters of interest, such as the radius of a transiting planet. With the ever growing availability of precise observations and the importance of robust estimates of astrophysical parameters, an emerging trend has been to freely fit the limb darkening coefficients (LDCs) describing a limb darkening law of choice, in order to propagate our ignorance of the true intensity profile. In practice, this approach has been limited to two-parameter limb darkening laws, such as the quadratic law, due to the relative ease of sampling the physically allowed range of LDCs. Here, we provide a highly efficient method for sampling LDCs describing a more accurate three-parameter non-linear law. We first derive analytic criteria which can quickly test if a set of LDCs are physical, although naive sampling with these criteria leads to an acceptance rate less than 1%. We then show that the loci of allowed LDCs can be transformed into a cone-like volume, from which we are able to draw uniform samples. We show that samples drawn uniformly from the conal region are physically valid in 97.3% of realizations and encompass 94.4% of the volume of allowed parameter space. We provide Python and Fortran code (LDC3) to sample from this region (and perform the reverse calculation) at https://github.com/davidkipping/LDC3, which also includes a subroutine to efficiently test whether a sample is physically valid or not.
We present the results of 938 speckle measures of double stars and suspected double stars drawn mainly from the Hipparcos Catalogue, as well as 208 observations where no companion was noted. One hundred fourteen pairs have been resolved for the first time. The data were obtained during four observing runs in 2014 using the Differential Speckle Survey Instrument (DSSI) at Lowell Observatory's Discovery Channel Telescope. The measurement precision obtained when comparing to ephemeris positions of binaries with very well-known orbits is generally less than 2 mas in separation and 0.5 degrees in position angle. Differential photometry is found to have internal precision of approximately 0.1 magnitudes and to be in very good agreement with Hipparcos measures in cases where the comparison is most relevant. We also estimate the detection limit in the cases where no companion was found. Visual orbital elements are derived for 6 systems.
Many results from astrophysical observations point to a 27% contribution of non-baryonic dark matter to the mass-energy budget of the universe. Although still elusive, strongly motivated candidates in form of weakly interacting massive particles could explain the nature of dark matter, and their annihilation or decay would give rise to detectable signatures in gamma-rays. In 2012, the H.E.S.S. collaboration started taking data with the largest imaging atmospheric Cherenkov telescope in the world which significantly lowered the energy threshold of the already operational four-telescope system. In particular, due to its location and improved performance at low energies, the H.E.S.S. experiment is now in a position to extend the search for dark matter line signals down to the 100 GeV mass range. The sensitivity to line searches with a new full likelihood approach will be discussed and preliminary results from observations with the second phase of H.E.S.S. will be presented.
We explore modifications to the current scenario for the slow neutron capture process in asymptotic giant branch (AGB) stars to account for the Pb deficiency observed in post-AGB stars of low metallicity ([Fe/H] ~ -1.2) and low initial mass (~ 1 - 1.5 Msun) in the Large and Small Magellanic Clouds. We calculated the stellar evolution and nucleosynthesis for a 1.3 Msun star with [Fe/H]=-1.3 and tested different amounts and distributions of protons leading to the production of the main neutron source within the 13C-pocket and proton ingestion scenarios. No s-process models can fully reproduce the abundance patterns observed in the post-AGB stars. When the Pb production is lowered the abundances of the elements between Eu and Pb, such as Er, Yb, W, and Hf, are also lowered to below those observed. Neutron-capture processes with neutron densities intermediate between the s and the rapid neutron-capture processes may provide a solution to this problem and be a common occurrence in low-mass, low-metallicity AGB stars.
We show that the Galactic latitude distribution of IceCube astrophysical neutrino events with energies above 100 TeV is inconsistent with the isotropic model of the astrophysical neutrino flux. Namely, the Galactic latitude distribution of the events shows an excess at low latitudes |b|<10 degrees and a deficit at high Galactic latitude |b|> 50 degrees. We use Monte-Carlo simulations to show that the inconsistency of the isotropic signal model with the data is at >3 sigma level.
Relic gravitational waves (RGWs) generated in the early Universe form a stochastic GW background, which can be directly probed by measuring the timing residuals of millisecond pulsars. In this paper, we investigate the constraints on the RGWs and on the inflationary parameters by the observations of current and potential future pulsar timing arrays. In particular, we focus on effects of various cosmic phase transitions (e.g. $e^{+}e^{-}$ annihilation, QCD transition and SUSY breaking) and relativistic free-streaming gases (neutrinos and dark fluids) in the general scenario of the early Universe, which have been neglected in the previous works. We find that the phase transitions can significantly damp the RGWs in the sensitive frequency range of pulsar timing arrays, and the upper limits of tensor-to-scalar ratio $r$ increase by a factor $\sim 2$ for both current and future observations. However, the effects of free-steaming neutrinos and dark fluids are all too small to be detected. Meanwhile, we find that, if the effective equation of state $w$ in the early Universe is larger than $1/3$, i.e. deviating from the standard hot big bang scenario, the detection of RGWs by pulsar timing arrays becomes much more promising.
In 2012, the High Energy Stereoscopic System (H.E.S.S.) was expanded by a fifth telescope (CT5). With an effective mirror diameter of 28m, CT5 is able to detect the Cherenkov light of very faint gamma-ray air showers, thereby significantly lowering the energy threshold of this telescope compared to the other four telescopes. Extracting as much information as possible from the recorded shower image is crucial for background rejection and to reach an energy threshold of a few tens of GeV. The camera of CT5 is conceived to register the time of the charge pulse maximum with respect to the beginning of the 16 ns integration window of each pixel. This information can be utilised to improve the event reconstruction. It also helps to reduce the background contamination at low energies. We present new techniques for background rejection based on CT5 timing information and evaluate their performance.
Some approaches to Quantum Gravity (QG) predict a modification of photon dispersion relations due to a breaking of Lorentz invariance. The effect is expected to affect photons near an effective QG energy scale. This scale has been constrained by observing gamma rays emitted from variable astrophysical sources such as gamma-ray bursts and flaring active galactic nuclei. Pulsars exhibit a periodic emission of possibly ms time scale. In 2014, the H.E.S.S. experiment reported the detection down to 20 GeV of gamma rays from the Vela pulsar having a periodicity of 89 ms. Using a likelihood analysis, calibrated with a dedicated Monte-Carlo procedure, we obtain the first limit on QG energy scale with the Vela pulsar. In this paper, the method and calibration procedure in use will be described and the results will be discussed.
The study of eclipsing binaries is our primary source of measured properties of normal stars, achieved through analysis of light and radial velocity curves of eclipsing systems. The study of oscillations and pulsations is increasingly vital for determining the properties of single stars, and investigating the physical phenomena active in their interiors. Combining the two methods holds the promise of establishing stringent tests of stellar evolutionary theory, and of calibrating model-dependent asteroseismology with empirically measured stellar properties. I review recent advances and outline future work.
Stellar seismology reveals some interior properties of thousands of solar-type stars but the solar seismic sound speed stays puzzling since a decade as it disagrees with the Standard Solar Model (SSM) prediction. One of the explanations of this disagreement may be found in the treatment of the transport of radiation from the solar core to the surface. As the same framework is used for other stars, it is important to check precisely the reliability of the interacting cross sections of photons with each species in order to ensure the energy transport for temperature T > 2 - 10$^6$ K and density $\rho$ > 0.2 g/cm$^3$. In this paper, we propose a new technique to reach the domain of temperature and density found in the solar radiative interior. This technique called the Double Ablation Front (DAF) is based on a high conversion of the laser energy into X-rays thanks to moderated Z material irradiated by laser intensity between 1.5 $\times$ 10$^{15}$ W/cm$^2$ and 4 $\times$ 10$^{15}$ W/cm$^2$. This high conversion creates, in addition to the electronic front a second ablation front in the moderated Z material. Between the two fronts there is a plateau of density and temperature that we exploit to heat a sample of iron or of oxide. The first simulations realized with the hydrodynamic code CHIC show that this technique allows to reach conditions equivalent to half the radiative zone of the Sun with high stability both in time and space. We examine the possibility to measure both iron and oxygen absorption spectra.
In this contribution talk we summarize the results of our ongoing project of detailed analysis of the chemical content (chemical tagging) as a promising powerful method to provide clear constraints on the membership of FGK kinematic candidates to stellar kinematic groups of different ages that can be used as an alternative or complementary to the methods that use kinematics, photometry or age indicators. This membership information is very important to better understand the star formation history in the solar neighborhood discerning between field-like stars (associated with dynamical resonances (bar) or spiral structure) and real physical structures of coeval stars with a common origin (debris of star-forming aggregates in the disk). We have already applied the chemical tagging method to constrain the membership of FGK candidate stars to the Hyades supercluster and the Ursa Major moving group and in this contribution we present the preliminary results of our study of the Castor moving group.
We review the orbital stability of the planar circular restricted three-body
problem, in the case of massless particles initially located between both
massive bodies. We present new estimates of the resonance overlap criterion and
the Hill stability limit, and compare their predictions with detailed dynamical
maps constructed with N-body simulations. We show that the boundary between
(Hill) stable and unstable orbits is not smooth but characterized by a rich
structure generated by the superposition of different mean-motion resonances
which does not allow for a simple global expression for stability.
We propose that, for a given perturbing mass $m_1$ and initial eccentricity
$e$, there are actually two critical values of the semimajor axis. All values
$a < a_{\rm Hill}$ are Hill-stable, while all values $a > a_{\rm unstable}$ are
unstable in the Hill sense. The first limit is given by the Hill-stability
criterion and is a function of the eccentricity. The second limit is virtually
insensitive to the initial eccentricity, and closely resembles a new resonance
overlap condition (for circular orbits) developed in terms of the intersection
between first and second-order mean-motion resonances.
Over the duration of the Kepler mission, KIC 8462852 was observed to undergo irregularly shaped, aperiodic dips in flux down to below the 20% level. The dipping activity can last for between 5 and 80 days. We characterize the object with high-resolution spectroscopy, spectral energy distribution fitting, and Fourier analyses of the Kepler light curve. We determine thatKIC 8462852 is a main-sequence F3 V/IV star, with a rotation period ~0.88 d, that exhibits no significant IR excess. In this paper, we describe various scenarios to explain the mysterious events in the Kepler light curve, most of which have problems explaining the data in hand. By considering the observational constraints on dust clumps orbiting a normal main-sequence star, we conclude that the scenario most consistent with the data is the passage of a family of exocomet fragments, all of which are associated with a single previous breakup event. We discuss the necessity of future observations to help interpret the system.
The design and performance of the LUX-ZEPLIN (LZ) detector is described as of March 2015 in this Conceptual Design Report. LZ is a second-generation dark-matter detector with the potential for unprecedented sensitivity to weakly interacting massive particles (WIMPs) of masses from a few GeV/c2 to hundreds of TeV/c2. With total liquid xenon mass of about 10 tonnes, LZ will be the most sensitive experiment for WIMPs in this mass region by the end of the decade. This report describes in detail the design of the LZ technical systems. Expected backgrounds are quantified and the performance of the experiment is presented. The LZ detector will be located at the Sanford Underground Research Facility in South Dakota. The organization of the LZ Project and a summary of the expected cost and current schedule are given.
We introduce a fully-coherent method for searching for gravitational wave signals generated by the merger of black hole and/or neutron star binaries. This extends the coherent analysis previously developed and used for targeted gravitational wave searches to an all-sky, all-time search. We apply the search to one month of data taken during the fifth science run of the LIGO detectors. We demonstrate an increase in sensitivity of 25% over the coincidence search, which is commensurate with expectations. Finally, we discuss prospects for implementing and running a coherent search for gravitational wave signals from binary coalescence in the advanced gravitational wave detector data.
In this work I briefly review some of the effects of hyperons on the properties of neutron and proto-neutron stars. In particular, I revise the problem of the strong softening of the EoS, and the consequent reduction of the maximum mass, induced by the presence of hyperons, a puzzle which has become more intringuing and difficult to solve due the recent measurements of the unusually high masses of the millisecond pulsars PSR J1903+0327 ($1.667\pm 0.021 M_\odot$), PSR J1614-2230 ($1.97 \pm 0.04 M_\odot$), and PSR J0348+0432 ($2.01 \pm 0.04 M_\odot$). Some of the solutions proposed to tackle this problem are discussed. Finally, I re-examine also the role of hyperons on the cooling properties of newly born neutron stars and on the so-called r-mode instability.
Links to: arXiv, form interface, find, astro-ph, recent, 1509, contact, help (Access key information)
We examine the host morphologies of heavily obscured active galactic nuclei (AGN) at $z\sim1$ to test whether obscured supermassive black hole growth at this epoch is preferentially linked to galaxy mergers. Our sample consists of 154 obscured AGN with $N_{\rm H}>10^{23.5}$ cm$^{-2}$ and $z<1.5$. Using visual classifications, we compare the morphologies of these AGN to control samples of moderately obscured ($10^{22}$ cm$^{-2}$ $<N_{\rm H}< 10^{23.5}$ cm$^{-2}$) and unobscured ($N_{\rm H}<10^{22}$ cm$^{-2}$) AGN. These control AGN are matched in redshift and intrinsic X-ray luminosity to our heavily obscured AGN. We find that heavily obscured AGN at z~1 are twice as likely to be hosted by late-type galaxies relative to unobscured AGN ($65.3^{+4.1}_{-4.6}\%$ vs $34.5^{+2.9}_{-2.7}\%$) and three times as likely to exhibit merger or interaction signatures ($21.5^{+4.2}_{-3.3}\%$ vs $7.8^{+1.9}_{-1.3}\%$). The increased merger fraction is significant at the 3.8$\sigma$ level. We also find that the incidence of point-like morphologies is inversely proportional to obscuration. If we exclude all point sources and consider only extended hosts, we find the correlation between merger fraction and obscuration is still evident, however at a reduced statistical significance ($2.5\sigma$). The fact that we observe a different disk/spheroid fraction versus obscuration indicates that viewing angle cannot be the only thing differentiating our three AGN samples, as a simple unification model would suggest. The increased fraction of disturbed morphologies with obscuration supports an evolutionary scenario, in which Compton-thick AGN are a distinct phase of obscured SMBH growth following a merger/interaction event. Our findings also suggest that some of the merger-triggered SMBH growth predicted by recent AGN fueling models may be hidden among the heavily obscured, Compton-thick population.
The nonlinear evolution of the Kelvin-Helmholtz instability is a popular test for code verification. To date, most Kelvin-Helmholtz problems discussed in the literature are ill-posed: they do not converge to any single solution with increasing resolution. This precludes comparisons among different codes and severely limits the utility of the Kelvin-Helmholtz instability as a test problem. The lack of a reference solution has led various authors to assert the accuracy of their simulations based on ad-hoc proxies, e.g., the existence of small-scale structures. This paper proposes well-posed Kelvin-Helmholtz problems with smooth initial conditions and explicit diffusion. We show that in many cases numerical errors/noise can seed spurious small-scale structure in Kelvin-Helmholtz problems. We demonstrate convergence to a reference solution using both Athena, a Godunov code, and Dedalus, a pseudo-spectral code. Problems with constant initial density throughout the domain are relatively straightforward for both codes. However, problems with an initial density jump (which are the norm in astrophysical systems) exhibit rich behavior and are more computationally challenging. In the latter case, Athena simulations are prone to an instability of the inner rolled-up vortex; this instability is seeded by grid-scale errors introduced by the algorithm, and disappears as resolution increases. Both Athena and Dedalus exhibit late-time chaos. Inviscid simulations are riddled with extremely vigorous secondary instabilities which induce more mixing than simulations with explicit diffusion. Our results highlight the importance of running well-posed test problems with demonstrated convergence to a reference solution. To facilitate future comparisons, we include the resolved, converged solutions to the Kelvin-Helmholtz problems in this paper in machine-readable form.
The so called "star forming main sequence" of galaxies is the apparent tight relationship between the star formation rate and stellar mass of a galaxy. Previous studies exclude galaxies which are not strictly "star forming" from the main sequence, because they do not lie on the same tight relation. Using local galaxies in the Sloan Digital Sky Survey we have classified galaxies according to their emission line ratios, and studied their location on the star formation rate - stellar mass plane. We find that galaxies form a sequence from the "blue cloud" galaxies which are actively forming stars, through a combination of composite, Seyfert, and LINER (Low-ionization nuclear emission-line region) galaxies, ending as "red-and-dead" galaxies. The sequence supports an evolutionary pathway for galaxies in which star formation quenching by active galactic nuclei (AGN) plays a key role.
The progenitors of many CCSNe are expected to be in binary systems. After the SN explosion, the companion may suffer from mass stripping and be shock heated as a result of the impact of the SN ejecta. If the binary system is disrupted, the companion is ejected as a runaway or hypervelocity star. By performing a series of 3D hydrodynamical simulations of the collision of SN ejecta with the companion star, we investigate how CCSN explosions affect their companions. We use the BEC code to construct the detailed companion structure at the time of SN explosion. The impact of the SN blast wave on the companion is followed by means of 3D SPH simulations using the Stellar GADGET code. For main-sequence (MS) companions, we find that the amount of removed mass, impact velocity, and chemical contamination of the companion that results from the impact of the SN ejecta, strongly increases with decreasing binary separation and increasing explosion energy. Their relationship can be approximately fitted by power laws, which is consistent with the results obtained from impact simulations of SNe Ia. However, we find that the impact velocity is sensitive to the momentum profile of the outer SN ejecta and, in fact, may decrease with increasing ejecta mass, depending on the modelling of the ejecta. Because most companions to Ib/c CCSNe are in their MS phase at the moment of the explosion, combined with the strongly decaying impact effects with increasing binary separation, we argue that the majority of such SNe lead to inefficient mass stripping and shock heating of the companion star following the impact of the ejecta. Our simulations show that the impact effects of Ib/c SN ejecta on the structure of MS companions, and thus their long term post-explosion evolution, is in general not be dramatic. We find that at most 10% of their mass is lost, and their resulting impact velocities are less than 100 km/s.
The uncertain origin of the recently-discovered 'changing-looking' quasar phenomenon - in which a luminous quasar dims significantly to a quiescent state in repeat spectroscopy over ~10 year timescales - may present unexpected challenges to our understanding of quasar accretion. To better understand this phenomenon, we take a first step to building a statistical sample of changing-look quasars with a systematic but simple archival search for these objects in the Sloan Digital Sky Survey Data Release 12. By leveraging the >10 year baselines for objects with repeat spectroscopy, we uncover two new changing-look quasars. Decomposition of the multi-epoch spectra and analysis of the broad emission lines suggest that the quasar accretion disk emission dims due to rapidly decreasing accretion rates, while disfavoring changes in intrinsic dust extinction. Narrow emission line energetics also support intrinsic dimming of quasar emission as the origin for this phenomenon rather than transient tidal disruption events. Although our search criteria included quasars at all redshifts and quasar transitions from either quasar-like to galaxy-like states or the reverse, all the most confident changing-look quasars discovered thus far have been relatively low-redshift (z ~ 0.2 - 0.3) and only exhibit quasar-like to galaxy-like transitions.
We carry out three dimensional radiation hydrodynamical simulations of gravitationally unstable discs to explore the movement of mass in a disc following its initial fragmentation. We find that the radial velocity of the gas in some parts of the disc increases by up to a factor of approximately 10 after the disc fragments, compared to before. While the movement of mass occurs in both the inward and outward directions, the inwards movement can cause the inner spirals of a self-gravitating disc to become sufficiently dense such that they can potentially fragment. This suggests that the dynamical behaviour of fragmented discs may cause subsequent fragmentation to occur at smaller radii than initially expected, but only after an initial fragment has formed in the outer disc.
Using observations from the MOSFIRE Deep Evolution Field (MOSDEF) survey, we investigate the physical conditions of star-forming regions in $z\sim2.3$ galaxies, specifically the electron density and ionization state. From measurements of the [O II]$\lambda\lambda$3726,3729 and [S II]$\lambda\lambda$6716,6731 doublets, we find a median electron density of $\sim250$ cm$^{-3}$ at $z\sim2.3$, an increase of an order of magnitude compared to measurements of galaxies at $z\sim0$. While $z\sim2.3$ galaxies are offset towards significantly higher O$_{32}$ values relative to local galaxies at fixed stellar mass, we find that the high-redshift sample follows a similar distribution to the low-metallicity tail of the local distribution in the O$_{32}$ vs. R$_{23}$ and O3N2 diagrams. Based on these results, we propose that $z\sim2.3$ star-forming galaxies have the same ionization parameter as local galaxies at fixed metallicity. In combination with simple photoionization models, the position of local and $z\sim2.3$ galaxies in excitation diagrams suggests that there is no significant change in the hardness of the ionizing spectrum at fixed metallicity from $z\sim0$ to $z\sim2.3$. We find that $z\sim2.3$ galaxies show no offset compared to low-metallicity local galaxies in emission line ratio diagrams involving only lines of hydrogen, oxygen, and sulfur, but show a systematic offset in diagrams involving [N II]$\lambda$6584. We conclude that the offset of $z\sim2.3$ galaxies from the local star-forming sequence in the [N II] BPT diagram is primarily driven by elevated N/O at fixed O/H compared to local galaxies. These results suggest that the local gas-phase and stellar metallicity sets the ionization state of star-forming regions at $z\sim0$ and $z\sim2$.
Cosmological N-body hydrodynamic computations following atomic and molecular chemistry (e$^-$, H, H$^+$, H$^-$, He, He$^+$, He$^{++}$, D, D$^+$, H$_2$, H$_2^+$, HD, HeH$^+$), gas cooling, star formation and production of heavy elements (C, N, O, Ne, Mg, Si, S, Ca, Fe, etc.) from stars covering a range of mass and metallicity are used to explore the origin of several chemical abundance patterns and to study both the metal and molecular content during simulated galaxy assembly. The resulting trends show a remarkable similarity to up-to-date observations of the most metal-poor damped Lyman-$\alpha$ absorbers at redshift $z\gtrsim 2$. These exhibit a transient nature and represent collapsing gaseous structures captured while cooling is becoming effective in lowering the temperature below $\sim 10^4\,\rm K$, before they are disrupted by episodes of star formation or tidal effects. Our theoretical results agree with the available data for typical elemental ratios, such as [C/O], [Si/Fe], [O/Fe], [Si/O], [Fe/H], [O/H] at redshifts $z\sim 2-7$. Correlations between HI and H$_2$ abundances show temporal and local variations and large spreads as a result of the increasing cosmic star formation activity from $z\sim 6$ to $z\sim 3$. The scatter we find in the abundance ratios is compatible with the observational data and is explained by simultaneous enrichment by sources from different stellar phases or belonging to different stellar populations. Simulated synthetic spectra support the existence of metal-poor cold clumps with large optical depth at $z\sim 6$ that could be potential population~III sites at low or intermediate redshift. The expected dust content is in line with recent determinations.
We report the discovery of a new "changing-look" quasar, SDSS J101152.98+544206.4, through repeat spectroscopy from the Time Domain Spectroscopic Survey. This is an addition to a small but growing set of quasars whose blue continua and broad optical emission lines have been observed to decline by a large factor on a time scale of approximately a decade. The 5100 Angstrom monochromatic continuum luminosity of this quasar drops by a factor of > 9.8 in a rest-frame time interval of < 9.7 years, while the broad H-alpha luminosity drops by a factor of 55 in the same amount of time. The width of the broad H-alpha line increases in the dim state such that the black hole mass derived from the appropriate single-epoch scaling relation agrees between the two epochs within a factor of 3. The fluxes of the narrow emission lines do not appear to change between epochs. The light curve obtained by the Catalina Sky Survey suggests that the transition occurs within a rest-frame time interval of approximately 500 days. We examine three possible mechanisms for this transition suggested in the recent literature. An abrupt change in the reddening towards the central engine is disfavored by the substantial difference between the timescale to obscure the central engine and the observed timescale of the transition. A decaying tidal disruption flare is consistent with the decay rate of the light curve but not with the prolonged bright state preceding the decay, nor can this scenario provide the power required by the luminosities of the emission lines. An abrupt drop in the accretion rate onto the supermassive black hole appears to be the most plausible explanation for the rapid dimming.
We present a catalog of GALEX Near-UV (NUV) and Far-UV (FUV) photometry for the Palomar/MSU and SDSS DR7 spectroscopic M dwarf catalogs. The catalog contains NUV measurements matched to 577 spectroscopically confirmed M dwarfs and FUV measurements matched to 150 spectroscopically confirmed M dwarfs. Using these data, we find that NUV and FUV luminosities strongly correlate with H{\alpha} emission, a typical indicator of magnetic activity in M dwarfs. We also examine the fraction of M dwarfs with varying degrees of strong line emission at NUV wavelengths. Our results indicate that the frequency of M dwarf NUV emission peaks at intermediate spectral types, with at least ~30% of young M4-M5 dwarfs having some level of activity. For mid-type M dwarfs, we show that NUV emission decreases with distance from the Galactic plane, a proxy for stellar age. Our complete matched source catalog is available online.
In this paper, we implement a perturbative approach, first proposed by Bouchet & Gispert (1999), to estimate variation of spectral index of galactic polarized synchrotron emission, using linear combination of simulated Stokes Q polarization maps of selected frequency bands from WMAP and Planck observations on a region of sky dominated by the synchrotron Stokes Q signal. We find that, a first order perturbative analysis recovers input spectral index map well. Along with the spectral index variation map our method provides a fixed reference index, \hat \beta_{0s}, over the sky portion being analyzed. Using Monte Carlo simulations we find that, <\hat \beta_{0s}> = -2.84 \pm 0.01, which matches very closely with position of a peak at \beta_s(p) = -2.85, of empirical probability density function of input synchrotron indices, obtained from the same sky region. For thermal dust, mean recovered spectral index, <\hat \beta_d> = 2.00 \pm 0.004, from simulations, matches very well with spatially fixed input thermal dust spectral index \beta_d = 2.00. As accompanying results of the method we also reconstruct CMB, thermal dust and a synchrotron template component with fixed spectral indices over the {\it entire} sky region. We use full pixel-pixel noise covariance matrices of all frequency bands, estimated from the sky region being analyzed, to obtain reference spectral indices for synchrotron and thermal dust, spectral index variation map, CMB map, thermal dust and synchrotron template components. The perturbative technique as implemented in this work has the interesting property that it can build a model to describe the data with an arbitrary but enough degree of accuracy (and precession) as allowed by the data. We argue that, our method of reference spectral index determination, CMB map, thermal dust and synchrotron template component reconstruction is a maximum likelihood method.
The GAPS project is running since 2012 with the goal to optimize the science return of the HARPS-N instrument mounted at Telescopio Nazionale Galileo. A large number of astronomers is working together to allow the Italian community to gain an international position adequate to the HARPS-N capabilities in the exoplanetary researches. Relevant scientific results are being obtained on both the main guidelines of the collaboration, i.e., the discovery surveys and the characterization studies. The planetary system discovered around the southern component of the binary XO-2 and its characterization together with that of the system orbiting the northern component are a good example of the completeness of the topics matched by the GAPS project. The dynamics of some planetary systems are investigated by studying the Rossiter-McLaughlin effect, while host stars are characterized by means of asteroseismology and star-planet interaction.
We use a modified version of MOOG to generate large grids of synthetic spectra in an attempt to derive quantitative abundances for three CVs (GK Per, RU Peg, and SS Cyg) by comparing the models to moderate resolution (R $\sim$ 25,000) $K$-band spectra obtained with NIRSPEC on Keck. For each of the three systems we find solar, or slightly sub-solar values for [Fe/H], but significant deficits of carbon: for SS Cyg we find [C/Fe] = $-$0.50, for RU Peg [C/Fe] = $-$0.75, and for GK Per [C/Fe] = $-$1.00. We show that it is possible to use lower resolution (R $\sim$ 2,000) spectra to quantify carbon deficits. We examine realistic veiling scenarios and find that emission from H I or CO cannot reproduce the observations.
One of the most important parameters determining the sensitivity of X-ray telescopes is their effective area as a function of the X-ray energy. The computation of the effective area of a Wolter-I mirror, with either a single layer or multilayer coating, is a very simple task for a source on-axis at astronomical distance. Indeed, when the source moves off-axis the calculation is more complicated, in particular for new hard X-ray imaging telescopes (NuSTAR, ASTRO-H, NHXM, IXO) beyond 10 keV, that will make use of multilayer coatings to extend the reflectivity band in grazing incidence. Unlike traditional single-layer coatings (in Ir or Au), graded multilayer coatings exhibit an oscillating reflectivity as a function of the incidence angle, which makes the effective area not immediately predictable for a source placed off-axis within the field of view. For this reason, the computation of the off-axis effective area has been so far demanded to ray- tracing codes, able to sample the incidence of photons onto the mirror assembly. Even if this approach should not be disdained, it would be interesting to approach the same problem from an analytical viewpoint. This would speed up and simplify the computation of the effective area as a function of the off-axis angle, a considerable advantage especially whenever the mirror parameters are still to be optimized. In this work we present the application of a novel, analytical formalism to the computation of the off-axis effective area and the grasp of the NHXM optical modules, requiring only the standard routines for the multilayer reflectivity computation.
We report Karl G. Jansky Very Large Array (VLA) and Atacama Large Millimeter Array (ALMA) spectroscopy in the redshifted molecular oxygen (O$_2$) 56.265~GHz and 424.763~GHz transitions from the $z=0.88582$ gravitational lens towards PKS\,1830$-$21. The ALMA non-detection of O$_2$ 424.763~GHz absorption yields the $3\sigma$ upper limit $N({\rm O}_2) \leq 5.8 \times 10^{17}$~cm$^{-2}$ on the O$_2$ column density, assuming that the O$_2$ level populations are thermalized at the gas kinetic temperature of 80~K. The VLA spectrum shows absorption by the CH$_3$CHO 56.185~GHz and 56.265~GHz lines, with the latter strongly blended with the O$_2$ 56.265~GHz line. Since the two CH$_3$CHO lines have the same equilibrium strength, we used the known CH$_3$CHO 56.185~GHz line profile to subtract out the CH$_3$CHO 56.265~GHz feature from the VLA spectrum, and then carried out a search for O$_2$ 56.265~GHz absorption in the residual spectrum. The non-detection of redshifted O$_2$ 56.265~GHz absorption in the CH$_3$CHO-subtracted VLA spectrum yields $N({\rm O}_2) \leq 2.3 \times 10^{17}$~cm$^{-2}$. Our $3\sigma$ limits on the O$_2$ abundance relative to H$_2$ are then $X({\rm O}_2) \leq 9.1 \times 10^{-6}$ (VLA) and $X({\rm O}_2) \leq 2.3 \times 10^{-5}$ (ALMA). These are $5-15$ times lower than the best previous constraint on the O$_2$ abundance in an external galaxy. The low O$_2$ abundance in the $z= 0.88582$ absorber may arise due to its high neutral carbon abundance and the fact that its molecular clouds appear to be diffuse or translucent clouds with low number density and high kinetic temperature.
We investigate the spectral properties of the UV (2650-3050 A) and optical (4000-5500 A) FeII emission features in a sample of 293 type 1 AGNs from SDSS database. We explore different correlations between their emission line properties, as well as the correlations with the other emission lines from the spectral range. We find several interesting correlations and we can outline the most interesting results as: (i) there is a kinematical connection between the UV and optical FeII lines indicating that the UV and optical FeII lines originate from the outer part of the Broad Line Region, so-called Intermediate Line Region (ILR); (ii) the unexplained anticorrelations of the optical FeII (EW FeII_opt vs. EW [OIII]5007A and EW FeII_opt vs. FWHM Hbeta) have not been detected for the UV FeII lines; (iii) the significant averaged redshift in the UV FeII lines, that is not present in optical FeII, indicates an inflow in the UV FeII emitting clouds and probably their asymmetric distribution; (iv) we confirm the anticorrelation between the intensity ratio of optical and UV FeII lines vs. FWHM of Hbeta, and we find the anticorrelations of this ratio with the widths of MgII 2800 A, optical FeII and UV FeII. This indicate a very important role of the column density and microturbulence in the emitting gas. We discuss the starburst activity in high density regions of young AGNs, as a possible explanation of the detected optical FeII correlations and intensity line ratios of the UV and optical FeII lines.
The variability of broad absorption lines is investigated for a sample of 188 broad-absorption-line (BAL) quasars (QSOs) ($z > 1.7$) with at least two-epoch observations from the Sloan Digital Sky Survey Data Release 7 (SDSS DR7), covering a time-scale of about 0.001 -- 3 years in the rest frame. Considering only the longest time-scale between epochs for each QSO, 73 variable regions in the \civ BAL troughs are detected for 43 BAL QSOs. The proportion of BAL QSOs showing variable regions increases with longer time-interval than about 1 year in the rest frame. The velocity width of variable regions is narrow compared to the BAL-trough outflow velocity. For 43 BAL QSOs with variable regions, it is found that there is a medium strong correlation between the variation of the continuum luminosity at 1500 \AA\ and the variation of the spectral index. With respect to the total 188 QSOs, larger proportion of BAL QSOs with variable regions appears bluer during their brighter phases, which implies that the origin of BAL variable regions is related to the central accretion process. For 43 BAL QSOs with variable regions, it is possible that there is a negative medium strong correlation between the absolute variation of the equivalent width and the \mgii-based black hole mass, and a medium strong correlation between the maximum outflow velocity of variable regions and the Eddington ratio. These results imply the connection between the BAL-trough variation and the central accretion process.
The variability of the broad absorption lines is investigated for a broad absorption line (BAL) QSO, SDSS J022844.09+000217.0 (z = 2.719), with 18 SDSS/BOSS spectra covering 4128 days in the observed frame. With the ratio of the rms spectrum to the mean spectrum, the relative flux change of the BAL-trough is larger than that of the emission lines and the continuum. Fitting the power-law continuum and the emission line profiles of \civ $\lambda$1549 and \siiv$\lambda$1399, we calculate the equivalent width (EW) for different epochs, as well as the continuum luminosity and the spectral index. It is found that there is a strong correlation between the BAL-trough EW and the spectral index, and a weak negative correlation between the BAL-trough EW and the continuum luminosity. The strong correlation between the BAL-trough EW and the spectral index for this one QSO suggests that dust is intrinsic to outflows. The weak correlation between the BAL variability and the continuum luminosity for this one QSO implies that the BAL-trough variation is not dominated by photoionization.
We investigate spatially resolved specific star formation rate (SSFR) in the inner $\sim$ 40 pc for a nearby Seyfert 2 galaxy, M51 (NGC 5194) by analyzing spectra obtained with the \emph{Hubble Space Telescope (HST)} Space Telescope Imaging Spectrograph (STIS). We present 24 radial spectra measured along the STIS long slit in M51, extending $\sim 1\arcsec$ from the nucleus (i.e., -41.5 pc to 39.4 pc). By the simple stellar population synthesis, the stellar contributions in these radial optical spectra are modeled. Excluding some regions with zero young flux fraction near the center (from -6 pc to 2 pc), we find that the mean flux fraction of young stellar populations (younger than 24.5 Myr) is about 9 \%, the mean mass fraction is about 0.09\%. The young stellar populations are not required in the center inner $\sim$ 8 pc in M51, suggesting a possible SSFR suppression in the circumnuclear region ($\sim$ 10 pc) from the feedback of active galactic nuclei (AGNs). The radial distribution of SSFR in M51 is not symmetrical with respect to the long slit in STIS. This unsymmetrical SSFR distribution is possibly due to the unsymmetrical AGN feedback in M51, which is related to its jet.
G73.9+0.9 has been classified as a probable shell-type supernova remnant (SNR), although it has also been suggested that this object could be a pulsar wind nebula (PWN). Here, a broadband model of the non-thermal emission of G73.9+0.9 from radio to gamma-rays is presented. The model includes a new gamma-ray observation obtained with analysis of 7 years of data from the Fermi-LAT telescope. Above 200 MeV the source is detected with a significance of 13-sigma and the spectrum of the radiation is best described by a power law with an index of -2.5. The leptonic mechanisms are hard to reconcile with the measured radio and gamma-ray SED. A PWN origin for the high-energy emission is also not very likely, due to the lack of detection of pulsars and of X-ray emission in the region, as well as from the shape of the gamma-ray spectrum. Given the possibility that the object is interacting with molecular clouds, a hadronic origin of the high-energy emission is more likely, and the spectral properties of the cosmic rays responsible for this radiation are derived.
In the current era of large surveys and massive data sets, autoclassification of astrophysical sources using intelligent algorithms is becoming increasingly important. In this paper we present the catalog of variable sources in the Third XMM-Newton Serendipitous Source catalog (3XMM) autoclassified using the Random Forest machine learning algorithm. We used a sample of manually classified variable sources from the second data release of the XMM-Newton catalogs (2XMMi-DR2) to train the classifier, obtaining an accuracy of ~92%. We also evaluated the effectiveness of identifying spurious detections using a sample of spurious sources, achieving an accuracy of ~95%. Manual investigation of a random sample of classified sources confirmed these accuracy levels and showed that the Random Forest machine learning algorithm is highly effective at automatically classifying 3XMM sources. Here we present the catalog of classified 3XMM variable sources. We also present three previously unidentified unusual sources that were flagged as outlier sources by the algorithm: a new candidate supergiant fast X-ray transient, a 400 s X-ray pulsar, and an eclipsing 5 hr binary system coincident with a known Cepheid.
Contributions of the Pierre Auger Collaboration to the 34th International Cosmic Ray Conference, 30 July - 6 August 2015, The Hague, The Netherlands
We present optical and near infrared observations of the type IIb supernova (SN) 2011fu from a few days to $\sim300$ d after explosion. The SN presents a double-peaked light curve (LC) similar to that of SN 1993J, although more luminous and with a longer cooling phase after the primary peak. The spectral evolution is also similar to SN 1993J's, with hydrogen dominating the spectra to $\sim40$ d, then helium gaining strength, and nebular emission lines appearing from $\sim60$ d post-explosion. The velocities derived from the P-Cygni absorptions are overall similar to those of other type IIb SNe. We have found a strong similarity between the oxygen and magnesium line profiles at late times, which suggests that these lines are forming at the same location within the ejecta. The hydrodynamical modelling of the pseudo-bolometric LC and the observed photospheric velocities suggest that SN 2011fu was the explosion of an extended star ($\rm R\sim450$ R$_\odot$), in which 1.3 $\times 10^{51}$ erg of kinetic energy were released and 0.15 M$_{\rm \odot}$ of $^{56}$Ni were synthesised. In addition, a better reproduction of the observed early pseudo-bolometric LC is achieved if a more massive H-rich-envelope than for other type IIb SNe is considered (0.3 M$_{\rm \odot}$). The hydrodynamical modelling of the LC and the comparison of our late-time spectra with nebular spectral models for type IIb SNe, point to a progenitor for SN 2011fu with a ZAMS mass of 13-18 M$_{\rm \odot}$.
As part of the TANAMI multiwavelength progam, we discuss new X-ray observations of the $\gamma$-ray and radio-loud Narrow Line Seyfert galaxy ($\gamma$-NLS1) PKS 2004-447. The active galaxy is a member of a small sample of radio-loud NLS1s detected in $\gamma$-rays by the Fermi Large Area Telescope. It is the radio-loudest and only southern-hemisphere source in this sample. We present results from our X-ray monitoring program comprised of Swift snapshot observations from 2012 through 2014 and two new X-ray observations with XMM-Newton in 2012. We analyze the X-ray spectrum and variability of this peculiar source using supplementary archival data from 2004 and 2011. The (0.5-10) keV spectrum is well described by a flat power law, which can be interpreted as non-thermal emission from a relativistic jet. The source exhibits moderate flux variability on timescales of both months and years. Correlated brightness variations in the (0.5-2) keV and (2-10) keV bands are explained by a single variable spectral component, such as the jet. A possible soft excess seen in the data from 2004 cannot be confirmed by the new \xmm{} observations in 2012. Any contribution to the total flux in 2004 is less than $20\%$ of the power-law component. The (0.5-10) keV luminosities of PKS 2004-447 are in the range of (0.5--2.7)$\times10^{44}\,\mathrm{erg\,s}^{-1}$. A comparison of the X-ray properties among the known $\gamma$-NLS1 galaxies shows that X-ray spectrum is typically dominated by a flat power law without intrinsic absorption. These objects are moderately variable in their brightness, while spectral variability is observed in at least two sources. The major difference across the X-ray spectra of $\gamma$-NLS1s is the luminosity, which spans a range of almost two orders of magnitude, from $10^{44}\,\mathrm{erg\,s}^{-1}$ to $10^{46}\,\mathrm{erg\,s}^{-1}$ in the (0.5-10) keV band.
We investigate the dust and gas distribution in the disk around HD 142527 based on ALMA observations of dust continuum, 13CO(3-2), and C18O(3-2) emission. The disk shows strong azimuthal asymmetry in the dust continuum emission, while gas emission is more symmetric. In this paper, we investigate how gas and dust are distributed in the dust-bright northern part of the disk and in the dust-faint southern part. We construct two axisymmetric disk models. One reproduces the radial profiles of the continuum and the velocity moments 0 and 1 of CO lines in the north and the other reproduces those in the south. We have found that the dust is concentrated in a narrow ring having ~50AU width (in FWHM; w_d=30AU in our parameter definition) located at ~170-200AU from the central star. The dust particles are strongly concentrated in the north. We have found that the dust surface density contrast between the north and south amounts to ~70. Compared to the dust, the gas distribution is more extended in the radial direction. We find that the gas component extends at least from ~100AU to ~250AU from the central star, and there should also be tenuous gas remaining inside and outside of these radii. The azimuthal asymmetry of gas distribution is much smaller than dust. The gas surface density differs only by a factor of ~3-10 between the north and south. Hence, gas-to-dust ratio strongly depends on the location of the disk: ~30 at the location of the peak of dust distribution in the south and ~3 at the location of the peak of dust distribution in the north. Despite large uncertainties, the overall gas-to-dust ratio is inferred to be ~10-30, indicating that the gas depletion may have already been under way.
Two habitable planetary states are proposed: an aqua planet like the Earth and a land planet that has a small amount of water. Land planets keep liquid water under larger solar radiation compared to aqua planets. Water loss may change an aqua planet into a land planet, and the planet can remain habitable for a longer time than had it stayed an aqua planet. We calculate planetary evolution with hydrogen escape for different initial water inventories and different distances from the central star. We find that there are two conditions to evolve an aqua planet into a land planet: the critical amount of water on the surface M_ml consistent with a planet being a land planet, and the critical amount of water vapor in the atmosphere M_cv that defines the onset of the runaway greenhouse state. We find that Earth-size aqua planets with initial oceans < 10 % of the Earth's can evolve into land planets if M_cv = 3 m in precipitable water and M_ml = 5 % of the Earth's ocean mass. Such planets can keep liquid water on their surface for another 2 Gyrs. The initial amount of water and M_cv are shown to be important dividing parameters of the planetary evolution path. Our results indicate that massive hydrogen escape could give a fresh start as another kind of habitable planet rather than the end of its habitability.
The propagation of charged particles, including cosmic rays, in a partially ordered magnetic field is characterized by a diffusion tensor whose components depend on the particle's Larmor radius $R_L$ and the degree of order in the magnetic field. This prescription relies explicitly on the assumption of a scale separation between random and mean magnetic fields, which usually applies in laboratory plasmas, but not in most astrophysical environments such as the interstellar medium (ISM). Direct estimates of the cosmic-ray diffusion tensor from test particle simulations have explored the range of particle energies corresponding to $10^{-2} \lesssim R_L/l_c \lesssim 10^{3}$, where $l_c$ is the magnetic correlation length. Modern simulations of the ISM have numerical resolution of order 1 pc, so the Larmor radius of the cosmic ray particles that dominate in their energy density is at least $10^{6}$ times smaller than the numerically resolved scales of the random magnetic field. Large-scale simulations of cosmic ray propagation in the ISM thus rely on oversimplified forms of the diffusion tensor. In this work we make first steps towards a more realistic description of cosmic ray diffusion in such simulations using test particle simulations in a random magnetic field. We provide explicit expressions for the cosmic ray diffusion tensor for $R_L/l_c \ll 1$ that may be used in a sub-grid model of cosmic ray diffusion and presumably a variety of other applications. We find that the diffusion coefficients are closely connected with existing transport theories and the random walk of magnetic lines.
A deep search for the potential glycine precursor hydroxylamine (NH$_2$OH) using the Caltech Submillimeter Observatory (CSO) at $\lambda = 1.3$ mm and the Combined Array for Research in Millimeter-wave Astronomy (CARMA) at $\lambda = 3$ mm is presented toward the molecular outflow L1157, targeting the B1 and B2 shocked regions. We report non-detections of NH$_2$OH in both sources. We a perform non-LTE analysis of CH$_3$OH observed in our CSO spectra to derive kinetic temperatures and densities in the shocked regions. Using these parameters, we derive upper limit column densities of NH$_2$OH of $\leq1.4 \times 10^{13}$~cm$^{-2}$ and $\leq1.5 \times 10^{13}$~cm$^{-2}$ toward the B1 and B2 shocks, respectively, and upper limit relative abundances of $N_{NH_2OH}/N_{H_2} \leq1.4 \times 10^{-8}$ and $\leq1.5 \times 10^{-8}$, respectively.
Intensity mapping of the neutral hydrogen brightness temperature promises to provide a three-dimensional view of the universe on very large scales. Nonlinear effects are typically thought to alter only the small-scale power, but we show how they can bias the extraction of cosmological information contained in the power spectrum on ultra-large scales. For linear perturbations to remain valid on large scales, we need to renormalize perturbations at higher order. In the case of intensity mapping, the second-order contribution to clustering from weak lensing dominates nonlinear contribution at high redshift. Renormalization modifies the mean brightness temperature and therefore the evolution bias. It also introduces a term that mimics white noise. These effects can influence forecasting analysis on ultra-large scales.
We study ten galaxy groups and clusters suggested in the literature to be "fossil system (FS)" based on \chandra\ observations. According to the $M_{500}-T$ and $L_{\rm X}-T$ relations, the gas properties of FSs are not physically distinct from ordinary galaxy groups or clusters. We also first study the $f_{\rm gas,~2500}-T$ relation and find that the FS exhibits same as ordinary systems. The gas densities of FSs within $0.1r_{200}$, are $\sim 10^{-3}$ cm$^{-3}$, which is the same order as galaxy clusters. The entropies within $0.1r_{200}$ ($S_{0.1r_{200}}$) of FSs are systematically lower than those in ordinary galaxy groups which is consistent with previous report, but we find their $S_{0.1r_{200}}-T$ relation is more similar to galaxy clusters. The derived mass profiles of FSs are consistent with the Navarro, Frenk, \& White model in $(0.1-1)r_{200}$, and the relation between scale radius $r_{\rm s}$ and characteristic mass density $ta_{\rm c}$ indicates the self-similarity of dark matter halos of FSs. The range of $r_{\rm s}$ and $ta_{\rm c}$ of FSs are also close to those of galaxy clusters. Therefore, FSs share more common characteristics with galaxy clusters. The special birth place of the FS makes it as a distinct galaxy system type.
Band split of solar type II radio bursts, discovered several decades ago, is a fascinating phenomenon with the type-II lanes exhibiting two almost-parallel sub-bands with similar morphology. The underlying split mechanism remains elusive. One popular interpretation is that the splitting bands are emitted from the shock upstream and downstream, respectively, with their frequency ratio ({\gamma}) determined by the shock compression ratio. This interpretation has been taken as the physical basis for many published references. Here we report an observational analysis of type II events with nice split selected from the ground-based RSTN data from 2001 to 2014, in the metric-decametric wavelength. We investigate the temporal variation and distribution of {\gamma}, and conduct correlation analyses on the deduced spectral values. It is found that {\gamma} varies in a very narrow range with >80% of {\gamma} (one-minute averaged data) being between 1.15 to 1.25. For some well-observed and long-lasting events, {\gamma} does not show a systematic variation trend within observational uncertainties, from the onset to the termination of the splits. In addition, the parameters representing the propagation speed of the radio source (presumably the coronal shock) show a very weak or basically no correlation with {\gamma}. We suggest that these results do not favor the upstreamdownstream scenario of band splits.
The formation mechanisms of thick discs are under discussion. Thick discs
might have formed either at high redshift on a short time-scale or might have
been built slowly over time. They may have an internal or an external origin.
Here we study in detail the kinematics and the stellar populations of the thick
disc of ESO533-4. ESO533-4 is a nearby bulgeless galaxy.
We present the first ever IFU study of an edge-on galaxy with enough depth to
study the thick disc. We exposed ESO533-4 with VIMOS@VLT for 6.5hours. The FOV
covered an axial extent 0.1-0.7r_25 (1-7kpc). We used pPXF and the MILES
library to obtain velocity and stellar population maps. We compared our
kinematic data with simple GADGET-2 models.
The apparent rotational lag of the thick disc of ESO533-4 is compatible with
that expected from the combinations of two effects: differential asymmetric
drift and the projection effects arising from studying a disc a few degrees
(2-3) away from edge-on. Thus, ESO533-4 contains little or no retrograde
material. This is compatible with three formation scenarii: the secular heating
of an initially thin disc, the formation of the thick disc at high redshift in
a turbulent disc phase, and its creation in a major merger event. If happening
in all galaxies, this last mechanism would cause retrograde thick discs in half
of them. Retrograde discs have not been observed in the five massive disc
galaxies (v_c>120km s^-1) for which thick disc kinematics are known. The
populations of the thin and the thick discs are separated in the
Age-log(Z/Z_Sun) plane. Thus, the thin and thick discs are made of two distinct
stellar populations. Although the stellar population results are not conclusive
due to the high dust extinction in ESO533-4, they do not favour a secular
evolution origin for the thick disc. Hence, we suggest that the thick disc of
ESO533-4 formed in a relatively short event (Abridged).
We explore the possibility to explain the properties of the Be/X-ray pulsars observed in the Small Magellanic Cloud within the magnetic levitation accretion scenario. This implies that their X-ray emission is powered by a wind-fed accretion onto a neutron star (NS) which captures matter from a magnetized stellar wind. The NS in this case is accreting matter from a non-keplerian magnetically levitating disc (ML-disc) which is surrounding its magnetosphere. This allows us to explain the observed periods of the pulsars in terms of spin equilibrium without the need of invoking dipole magnetic fields outside the usual range ~ 10^11- 10^13 G inferred from cyclotron features of Galactic high mass X-ray binaries. We find that the equilibrium period of a NS, under certain conditions, depends strongly on the magnetization of the stellar wind of its massive companion and, correspondingly, on the magnetic field of the massive companion itself. This may help to explain why similar NSs in binaries with similar properties rotate with different periods yielding a large scatter of periods of the accretion-powered pulsar observed in SMC and our galaxy.
Amongst the population of TeV gamma-ray sources detected with the High Energy Stereoscopic System (H.E.S.S.) in the Galactic plane, clearly identified supernova remnant (SNR) shells constitute a small but precious source class. TeV-selected SNRs are prime candidates for sources of efficient cosmic-ray acceleration. In this work, we present new SNR candidates that have been identified in the entire H.E.S.S. phase I data set of the Galactic plane recorded over the past ten years. Identification with a known SNR shell candidate was successful for one new source, HESS J1534-571. In other cases, TeV-only shell candidates are challenging to firmly identify as SNRs due to their lack of detected non-thermal emission in lower energy bands. We will discuss how these objects may present an important link between young and evolved SNRs, since their shell emission may be dominated by hadronic processes.
In this paper, we apply a new statistical analysis technique, Mean Field approach to Bayesian Independent Component Analysis (MF-ICA), on galaxy spectral analysis. This algorithm can compress the stellar spectral library into a few Independent Components (ICs), and galaxy spectrum can be reconstructed by these ICs. Comparing to other algorithms which decompose a galaxy spectrum into a combination of several simple stellar populations, MF-ICA approach offers a large improvement in the efficiency. To check the reliability of this spectral analysis method, three different methods are used: (1) parameter-recover for simulated galaxies, (2) comparison with parameters estimated by other methods, and (3) consistency test of parameters from the Sloan Digital Sky Survey galaxies. We find that our MF-ICA method not only can fit the observed galaxy spectra efficiently, but also can recover the physical parameters of galaxies accurately. We also apply our spectral analysis method to the DEEP2 spectroscopic data, and find it can provide excellent fitting for those low signal-to-noise spectra.
In the paper we investigate three novel rising submillimeter (THz) bursts occurred sequentially in a super-Active Region NOAA 10486. The average rising rate of the flux density above 200 GHz is only 20 sfu/GHz (corresponding spectral index $\alpha$ of 1.6) for the THz spectral components of 2003 October 28 and November 4 bursts, while it can attain values of 235 sfu/GHz ($\alpha$=4.8) for 2003 November 2 burst. The steeply rising THz spectrum can be produced by a population of high relativistic electrons with a low-energy cutoff of 1 MeV , while it only requires a low-energy cutoff of 30 keV for the two slowly rising THz bursts, via gyrosynchrotron (GS) radiation based on our numerical simulations of burst spectra in the magnetic dipole field case. The electron density variation is much larger in the THz source than that in microwave (MW) one. It is interesting that the THz source radius decreased by 20--50$\%$ during the decay phase for the three events, but the MW one increased by 28$\%$ for the 2003 November 2 event. In the paper we will present a calculation formula of energy released by ultrarelativistic electrons, accounting the relativistic correction for the first time. We find that the energy released by energetic electrons in the THz source exceeds that in microwave one due to the strong GS radiation loss at THz range, although the modeled THz source area is 3--4 orders smaller than the modeled MW one. The total energies released by energetic electrons via the GS radiation in radio sources are estimated, respectively, to be $5.2\times10^{33}$, $3.9\times10^{33}$ and $3.7\times10^{32}$ erg for the October 28, November 2 and 4 bursts, which are 131, 76 and 4 times as large as the thermal energies of $2.9\times10^{31}$, $2.1\times10^{31}$ and $5.2\times10^{31}$ erg estimated from the soft x-ray GOES observations.
In this paper, we point out and study a generic type of signals existing in the primordial universe models, which can be used to model-independently distinguish the inflation scenario from alternatives. These signals are generated by massive fields that function as standard clocks. The role of massive fields as standard clocks has been realized in previous works. Although the existence of such massive fields is generic, the previous realizations require sharp features to classically excite the oscillations of the massive clock fields. Here, we point out that the quantum fluctuations of massive fields can actually serve the same purpose as the standard clocks. We show that they are also able to directly record the defining property of the scenario type, namely, the scale factor of the primordial universe as a function of time a(t), but through shape-dependent oscillatory features in non-Gaussianities. Since quantum fluctuating massive fields exist in any realistic primordial universe models, these quantum primordial standard clock signals are present in any inflation models, and should exist quite generally in alternative-to-inflation scenarios as well. However, the amplitude of such signals is very model-dependent.
The Evolutionary Map of the Universe (EMU) is a proposed radio continuum survey of the Southern Hemisphere up to declination +30 deg., with the Australian Square Kilometre Array Pathfinder (ASKAP). EMU will use an automated source identification and measurement approach that is demonstrably optimal, to maximise the reliability, utility and robustness of the resulting radio source catalogues. As part of the process of achieving this aim, a "Data Challenge" has been conducted, providing international teams the opportunity to test a variety of source finders on a set of simulated images. The aim is to quantify the accuracy of existing automated source finding and measurement approaches, and to identify potential limitations. The Challenge attracted nine independent teams, who tested eleven different source finding tools. In addition, the Challenge initiators also tested the current ASKAPsoft source-finding tool to establish how it could benefit from incorporating successful features of the other tools. Here we present the results of the Data Challenge, identifying the successes and limitations for this broad variety of the current generation of radio source finding tools. As expected, most finders demonstrate completeness levels close to 100% at 10sigma dropping to levels around 10% by 5sigma. The reliability is typically close to 100% at 10sigma, with performance to lower sensitivities varying greatly between finders. All finders demonstrate the usual trade-off between completeness and reliability, whereby maintaining a high completeness at low signal-to-noise comes at the expense of reduced reliability, and vice-versa. We conclude with a series of recommendations for improving the performance of the ASKAPsoft source-finding tool.
We investigate the dynamics of self-gravity wakes in dense planetary rings. In particular, we examine how the pitch angle of self-gravity wakes depend on ring parameters using N-body simulations. We calculate the pitch angles using the two-dimensional autocorrelation function of the ring surface density. We obtain the pitch angles for the inner and outer parts of the autocorrelation function separately. We confirm that the pitch angles are 15 to 30 degrees for reasonable ring parameters, which are consistent with previous studies. We find that the inner pitch angle increases with the Saturnicentric distance, while it barely depends on the optical depth and the restitution coefficient of ring particles. The increase of the inner pitch angle with the Saturnicentric distance is consistent with the observations of the A ring. The outer pitch angle does not have the clear dependence on any ring parameters and is about 10 - 15 degrees. This value is consistent with the pitch angle of spiral arms in collisionless systems.
Recent observations have discovered a number of extremely gas-rich very faint dwarf galaxies possibly embedded in low-mass dark matter halos. We investigate star formation histories of these gas-rich dwarf ("almost dark") galaxies both for isolated and interacting/merging cases. We find that although star formation rates (SFRs) are very low (<10^-5 M_sun/yr) in the simulated dwarfs in isolation for the total halo masses (M_h) of 10^8-10^9 M_sun, they can be dramatically increased to be ~ 10^{-4} M_sun/yr when they interact or merge with other dwarfs. These interacting faint dwarfs with central compact HII regions can be identified as isolated emission line dots ("ELdots") owing to their very low surface brightness envelopes of old stars. The remnant of these interacting and merging dwarfs can finally develop central compact stellar systems with very low metallicities (Z/Z_sun<0.1), which can be identified as extremely metal-deficient ("XMD") dwarfs. These results imply that although there would exist many faint dwarfs that can be hardly detected in the current optical observations, they can be detected as isolated ELdots or XMD dwarfs, when they interact with other galaxies and their host environments. We predict that nucleated ultra-faint dwarfs formed from merging of almost dark dwarfs can be identified as globular clusters owing to the very low surface brightness stellar envelopes, if they are distant objects.
We report the discovery of two methane dwarfs in the dark cloud L 1688 of the Rho Oph star-forming region. The two objects were among the T dwarf candidates with possible methane absorption and cool atmospheres, as diagnosed by infrared colors using deep WIRCam/CFHT HK plus CH4ON images, and IRAC/Spitzer c2d data. Follow-up spectroscopic observations with the FLAMINGOS-2/Gemini South confirmed the methane absorption at 1.6 micron. Compared with spectral templates of known T dwarfs in the field, i.e., of the old populations, Oph J162738-245240 (Oph-T3) is a T0/T1 type, whereas Oph J162645-241949 (Oph-T17) is consistent with a T3/T4 type in the H band but an L8/T1 in the K band. Compared with the BT-Settl model, both Oph-T3 and Oph-T17 are consistent with being cool, ~ 1000 K and ~ 900 K, respectively, and of low surface gravity, log(g) = 3.5. With an age no more than a couple Myr, these two methane dwarfs thereby represent the youngest T dwarfs ever confirmed. A young late L dwarf, OphJ162651-242110, was found serendipitously in our spectroscopic observations.
Blazars are the most abundant class of known extragalactic very-high-energy (VHE, E>100 GeV) gamma-ray sources. However, one of the biggest difficulties in investigating their VHE emission resides in their limited number, since less than 60 of them are known by now. In this contribution we report on H.E.S.S. observations of the BL Lac object PKS 1440-389. This source has been selected as target for H.E.S.S. based on its high-energy gamma-ray properties measured by Fermi-LAT. The extrapolation of this bright, hard-spectrum gamma-ray blazar into the VHE regime made a detection on a relatively short time scale very likely, despite its uncertain redshift. H.E.S.S. observations were carried out with the 4-telescope array from February to May 2012 and resulted in a clear detection of the source. Contemporaneous multi-wavelength data are used to construct the spectral energy distribution of PKS 1440-389 which can be described by a simple one-zone synchrotron-self Compton model.
We present preliminary results on Dark Matter searches from observations of the Perseus galaxy cluster with the MAGIC Telescopes. MAGIC is a system of two Imaging Atmospheric Cherenkov Telescopes located in the Canary island of La Palma, Spain. Galaxy clusters are the largest known gravitationally bound structures in the Universe, with masses of ~10^15 Solar masses. There is strong evidence that galaxy clusters are Dark Matter dominated objects, and therefore promising targets for Dark Matter searches, particularly for decay signals. MAGIC has taken almost 300 hours of data on the Perseus Cluster between 2009 and 2015, the deepest observational campaign on any galaxy cluster performed so far in the very high energy range of the electromagnetic spectrum. We analyze here a small sample of this data and search for signs of dark matter in the mass range between 100 GeV and 20 TeV. We apply a likelihood analysis optimized for the spectral and morphological features expected in the dark matter decay signals. This is the first time that a dedicated Dark Matter optimization is applied in a MAGIC analysis, taking into account the inferred Dark Matter distribution of the source. The results with the full dataset analysis will be published soon by the MAGIC Collaboration.
We reconsider the origin and processing of dust in elliptical galaxies. We theoretically formulate the evolution of grain size distribution, taking into account dust supply from asymptotic giant branch (AGB) stars and dust destruction by sputtering in the hot interstellar medium (ISM), whose temperature evolution is treated by including two cooling paths: gas emission and dust emission (i.e. gas cooling and dust cooling). With our new full treatment of grain size distribution, we confirm that dust destruction by sputtering is too efficient to explain the observed dust abundance even if AGB stars continue to supply dust grains, and that, except for the case where the initial dust-to-gas ratio in the hot gas is as high as $\sim 0.01$, dust cooling is negligible compared with gas cooling. However, we show that, contrary to previous expectations, cooling does not help to protect the dust; rather, the sputtering efficiency is raised by the gas compression as a result of cooling. We additionally consider grain growth after the gas cools down. Dust growth by the accretion of gas-phase metals in the cold medium increase the dust-to-gas ratio up to $\sim 10^{-3}$ if this process lasts >10/(n_H/10^3 cm^{-3}) Myr, where $n_\mathrm{H}$ is the number density of hydrogen nuclei. We show that the accretion of gas-phase metals is a viable mechanism of increasing the dust abundance in elliptical galaxies to a level consistent with observations, and that the steepness of observed extinction curves is better explained with grain growth by accretion.
Detection of the signature of Rayleigh scattering in the transmission spectrum of an exoplanet is increasingly becoming the target of observational campaigns because the spectral slope of the Rayleigh continuum enables one to determine the scaleheight of its atmosphere in the absence of hazes. However, this is only true when one ignores the refractive effects of the exoplanet's atmosphere. I illustrate with a suite of simple isothermal clear Jovian H2-He atmosphere models with various abundances of water that refraction can decrease significantly the spectral slope of the Rayleigh continuum and that it becomes flat in the infrared. This mimics a surface, or an optically thick cloud deck, at much smaller pressures than one can probe in the non-refractive case. Although the relative impact of refraction on an exoplanet's transmission spectrum increases with decreasing atmospheric temperatures as well as increasing stellar temperature, it is still quite important from a retrieval's perspective even for a Jovian-like planet with an atmospheric temperature as high as 1200 K. Indeed, the flat Rayleigh continuum in the infrared breaks in large part the retrieval degeneracy between abundances of chemical species and the planet's radius because the size of spectral features increases significantly with abundances, in stark contrast with the non-refractive case which simply shifts them to a larger or smaller effective radius. Abundances inferred assuming the atmosphere is cloud-free are lower limits. These results show how important it is to include refraction in retrieval algorithms to interpret transmission spectra of gas giants accurately.
We present the largest to date sample of hydrogen Lyman continuum (LyC) emitting galaxies at any redshift, with $18$ LyA Emitters (LAEs) and $7$ Lyman Break Galaxies (LBGs), obtained from the SSA22 field with Subaru/Suprime-Cam. The sample is based on the $136$ LBGs and $159$ LAEs observed in the field, all with spectroscopically confirmed redshifts, and they are selected as galaxies with counterpart in a narrow-band filter image which traces LyC at $z\geq 3.06$. Many LyC candidates show a spatial offset between the rest-frame non-ionizing ultraviolet (UV) detection and the LyC-emitting substructure or between the LyA emission and LyC. Statistically it is highly unlikely that all candidates in our sample are contaminants, and there should be $\sim9$ and $\sim2$ viable LyC candidates among the LAEs and LBGs. There is some evidence for a positive LyC/LyA correlation, suggesting that both LyC and LyA escape via a similar mechanism. "Standard" SED models cannot explain the observed LyC LAEs colors, instead requiring more exotic models like a top-heavy IMF with an average stellar mass of $\sim100M_{sun}$ and the "Lyman limit bump" model with an escape of nebular recombination LyC. The LyC LBGs colors are consistent with a Salpeter IMF with no nebular emission. Both types of galaxies seem to require extremely metal-poor or metal-free young (zero age) stellar populations. We have obtained estimates of LyC escape fraction of $\sim20\%$ for LyC LBGs and $\sim30\%$ for LyC LAEs, assuming an SMC dust law. We cannot simultaneously explain the observed $f_{LyC}/f_{UV}$ flux density ratios and the UV slopes if we assume a Calzetti attenuation curve. The LyC emission seems to be bimodal - stacking non-detections reveals no significant LyC signal strength. The $3\sigma$ median upper limits on the flux density ratios from non-detections are [...]
We present the first estimate of age, stellar metallicity and chemical abundance ratios, for an individual early-type galaxy at high-redshift (z = 1.426) in the COSMOS field. Our analysis is based on observations obtained with the X-Shooter instrument at the VLT, which cover the visual and near infrared spectrum at high (R >5000) spectral resolution. We measure the values of several spectral absorptions tracing chemical species, in particular Magnesium and Iron, besides determining the age-sensitive D4000 break. We compare the measured indices to stellar population models, finding good agreement. We find that our target is an old (t > 3 Gyr), high-metallicity ([Z/H] > 0.5) galaxy which formed its stars at z_{form} > 5 within a short time scale ~0.1 Gyr, as testified by the strong [\alpha/Fe] ratio ( > 0.4), and has passively evolved in the first > 3-4 Gyr of its life. We have verified that this result is robust against the choice and number of fitted spectral features, and stellar population model. The result of an old age and high-metallicity has important implications for galaxy formation and evolution confirming an early and rapid formation of the most massive galaxies in the Universe.
V745 Sco is a recurrent nova, with the most recent eruption occurring in February 2014. V745 Sco was first observed by Swift a mere 3.7 hr after the announcement of the optical discovery, with the super-soft X-ray emission being detected around four days later and lasting for only ~two days, making it both the fastest follow-up of a nova by Swift and the earliest switch-on of super-soft emission yet detected. Such an early switch-on time suggests a combination of a very high velocity outflow and low ejected mass and, together with the high effective temperature reached by the super-soft emission, a high mass white dwarf (>1.3 M_sun). The X-ray spectral evolution was followed from an early epoch where shocked emission was evident, through the entirety of the super-soft phase, showing evolving column density, emission lines, absorption edges and thermal continuum temperature. UV grism data were also obtained throughout the super-soft interval, with the spectra showing mainly emission lines from lower ionization transitions and the Balmer continuum in emission. V745 Sco is compared with both V2491 Cyg (another nova with a very short super-soft phase) and M31N 2008-12a (the most rapidly recurring nova yet discovered). The longer recurrence time compared to M31N 2008-12a could be due to a lower mass accretion rate, although inclination of the system may also play a part. Nova V745 Sco (2014) revealed the fastest evolving super-soft source phase yet discovered, providing a detailed and informative dataset for study.
We use an existing laboratory facility for space hardware calibration in vacuum to study the impact of energetic ions on water ice. The experiment is intended to simulate the conditions on the surface of Jupiter's icy moons. We present first results of ion sputtering in a sample of porous ice, including the first experimental results for sulphur ion sputtering of ice. The results confirm theoretical predictions and extrapolations from previous sputtering experiments obtained at different impact angles for non-porous water ice.
Asteroids are frequently colliding with small projectiles. Although each individual small collision is not very important, their cumulative effect can substantially change topography and also the overall shape of an asteroid. We run simulations of random collisions onto a single target asteroid represented by triaxial ellipsoid. We investigated asteroids of several hundred meters to about 18 km in diameter for which we assumed all material excavated by the collision to escape the asteroid. The cumulative effect of these collisions is an increasing elongation of the asteroid figure. However, the estimated timescale of this process is much longer than the collisional lifetime of asteroids. Therefore, we conclude that small collisions are probably not responsible for the overall shape of small asteroids.
The cross section of material in debris discs is thought to be dominated by the smallest grains that can still stay in bound orbits despite the repelling action of stellar radiation pressure. Thus the minimum (and typical) grain size $s_\text{min}$ is expected to be close to the radiation pressure blowout size $s_\text{blow}$. Yet a recent analysis of a sample of Herschel-resolved debris discs showed the ratio $s_\text{min}/s_\text{blow}$ to systematically decrease with the stellar luminosity from about ten for solar-type stars to nearly unity in the discs around the most luminous A-type stars. Here we explore this trend in more detail, checking how significant it is and seeking to find possible explanations. We show that the trend is robust to variation of the composition and porosity of dust particles. For any assumed grain properties and stellar parameters, we suggest a recipe of how to estimate the "true" radius of a spatially unresolved debris disc, based solely on its spectral energy distribution. The results of our collisional simulations are qualitatively consistent with the trend, although additional effects may also be at work. In particular, the lack of grains with small $s_\text{min}/s_\text{blow}$ for lower luminosity stars might be caused by the grain surface energy constraint that should limit the size of the smallest collisional fragments. Also, a better agreement between the data and the collisional simulations is achieved when assuming debris discs of more luminous stars to have higher dynamical excitation than those of less luminous primaries. This would imply that protoplanetary discs of more massive young stars are more efficient in forming big planetesimals or planets that act as stirrers in the debris discs at the subsequent evolutionary stage.
Radio interferometers suffer from the problem of missing information in their data, due to the gaps between the antennas. This results in artifacts, such as bright rings around sources, in the images obtained. Multiple deconvolution algorithms have been proposed to solve this problem and produce cleaner radio images. However, these algorithms are unable to correctly estimate uncertainties in derived scientific parameters or to always include the effects of instrumental errors. We propose an alternative technique called Bayesian Inference for Radio Observations (BIRO) which uses a Bayesian statistical framework to determine the scientific parameters and instrumental errors simultaneously directly from the raw data, without making an image. We use a simple simulation of Westerbork Synthesis Radio Telescope data including pointing errors and beam parameters as instrumental effects, to demonstrate the use of BIRO.
We have analysed the Fermi LAT data on the SNR G73.9+0.9. We have confirmed a previous detection of high-energy gamma-rays from this source at a high significance of $\simeq 12\sigma$. The observed spectrum shows a significant curvature, peaking in $E F_E$ at $\sim$1 GeV. We have also calculated the flux upper limits in the mm-wavelength and X-ray ranges from Planck and XMM-Newton, respectively. We have inspected the intensity of the CO (1$\rightarrow $0) emission line and found a large peak at a velocity range corresponding to the previously estimated source distance of $\sim$4 kpc, which may indicate an association between a molecular cloud and the SNR. The gamma-ray emission appears due to interaction of accelerated particles within the SNR with the matter of the cloud. The most likely radiative process responsible for the gamma-ray emission is decay of neutral pions produced in ion-ion collisions. While a dominant leptonic origin of this emission can be ruled out, the relativistic electron population related to the observed radio flux will necessarily lead to a certain level of bremsstrahlung gamma-ray emission. Based on this broad-band modelling, we have developed a method to estimate the magnetic field, yielding $B\geq 80\,\mu$G at our best estimate of the molecular cloud density (or less at a lower density). G73.9+0.9 appears similar, though somewhat weaker, to other SNRs interacting with a local dense medium detected by the LAT.
We measure the weak gravitational lensing shear power spectra and their cross-power in two photometric redshift bins from the Canada-France-Hawaii Telescope Lensing Survey (CFHTLenS). The measurements are performed directly in multipole space in terms of adjustable band powers. For the extraction of the band powers from the data we have implemented and extended a quadratic estimator, a maximum likelihood method that allows us to readily take into account irregular survey geometries, masks, and varying sampling densities. We find the 68 per cent credible intervals in the $\sigma_8$-$\Omega_{\rm m}$-plane to be marginally consistent with results from $Planck$ for a simple five parameter $\Lambda$CDM model. For the projected parameter $S_8 \equiv \sigma_8(\Omega_{\rm m}/0.3)^{0.5}$ we obtain a best-fitting value of $S_8 = 0.768_{-0.039}^{+0.045}$. This constraint is consistent with results from other CFHTLenS studies as well as the Dark Energy Survey. Our most conservative model, including modifications to the power spectrum due to baryon feedback and marginalization over photometric redshift errors, yields an upper limit on the total mass of three degenerate massive neutrinos of $\Sigma m_\nu < 4.53 \, {\rm eV}$ at 95 per cent credibility, while a Bayesian model comparison does not favour any model extension beyond a simple five parameter $\Lambda$CDM model. Combining the shear likelihood with $Planck$ breaks the $\sigma_8$-$\Omega_{\rm m}$-degeneracy and yields $\sigma_8=0.817_{-0.014}^{+0.013}$ and $\Omega_{\rm m} = 0.298 \pm 0.011$ which is fully consistent with results from $Planck$ alone.
Magnetohydronamic turbulence is believed to play a crucial role in heating the laboratorial, space, and astrophysical plasmas. However, the precise connection between the turbulent fluctuations and the particle kinetics has not yet been established. Here we present clear evidence of plasma turbulence heating based on diagnosed wave features and proton velocity distributions from solar wind measurements by the Wind spacecraft. For the first time, we can report the simultaneous observation of counter-propagating magnetohydrodynamic waves in the solar wind turbulence. Different from the traditional paradigm with counter-propagating Alfv\'en waves, anti-sunward Alfv\'en waves (AWs) are encountered by sunward slow magnetosonic waves (SMWs) in this new type of solar wind compressible turbulence. The counter-propagating AWs and SWs correspond respectively to the dominant and sub-dominant populations of the imbalanced Els\"asser variables. Nonlinear interactions between the AWs and SMWs are inferred from the non-orthogonality between the possible oscillation direction of one wave and the possible propagation direction of the other. The associated protons are revealed to exhibit bi-directional asymmetric beams in their velocity distributions: sunward beams appearing in short and narrow patterns and anti-sunward broad extended tails. It is suggested that multiple types of wave-particle interactions, i.e., cyclotron and Landau resonances with AWs and SMWs at kinetic scales, are taking place to jointly heat the protons perpendicularly and parallel.
The High Energy Stereoscopic System (H.E.S.S.) experiment is one of the largest observatories for gamma-ray astronomy. It consists of four telescopes with a reflecting dish diameter of 12m (CT1 to CT4) and a newer large telescope (CT5) with a reflecting dish diameter of 28m. On CT5 876 mirror facets are mounted, all of them equipped with a computerised system for their alignment. The design of the mirror alignment and control system and the performance of the hardware installed to the telescope are presented. Furthermore the achieved point spread function of the telescope over the full operational elevation range as well as the stability of the alignment over an extended period of time are shown.
The presence of dark matter in the Universe is nowadays widely supported by a large body of astronomical and cosmological observations. One of the best targets to look for dark matter particle self-annihilation into very high energy gamma-rays is the Galactic Center (GC) region. A search for annihilating dark matter in the central 300 parsecs around the GC is performed with the H.E.S.S. array of ground-based Cherenkov telescopes. Using the full H.E.S.S.- I dataset (2004-2014) for this region, new constraints are derived on the velocity-weighted annihilation cross section with a 2D likelihood method, taking advantage of differences in both spectral and spatial morphologies of signal and background. Higher statistics from the 10-years GC dataset of H.E.S.S. I together with a novel analysis technique, allow to improve the present constraints by a factor of 2 for a DM particle mass of 1 TeV. The expected improvement in sensitivity of the new analysis technique, applied to performances of the H.E.S.S.-II array, is also presented.
The NectarCAM is a camera proposed for the medium-sized telescopes in the framework of the Cherenkov Telescope Array (CTA), the next-generation observatory for very-high-energy gamma-ray astronomy. The cameras are designed to operate in an open environment and their mechanics must provide protection for all their components under the conditions defined for the CTA observatory. In order to operate in a stable environment and ensure the best physics performance, each NectarCAM will be enclosed in a slightly overpressurized, nearly air-tight, camera body, to prevent dust and water from entering. The total power dissipation will be ~7.7 kW for a 1855-pixel camera. The largest fraction is dissipated by the readout electronics in the modules. We present the design and implementation of the cooling system together with the test bench results obtained on the NectarCAM thermal demonstrator.
We report detections and constraints for the near infrared Ks band secondary eclipses of seven hot-Jupiters using the IRIS2 infrared camera on the Anglo-Australian Telescope. Eclipses in the Ks band for WASP-18b and WASP-36b have been measured for the first time. We also present new measurements for the eclipses of WASP-4b, WASP-5b, and WASP-46b, as well as upper limits for the eclipse depths of WASP-2b and WASP-76b. In particular, two full eclipses of WASP-46b were observed, allowing us to demonstrate the repeatability of our observations via independent analyses on each eclipse. Significant numbers of eclipse depths for hot-Jupiters have now been measured in both Ks and the four Spitzer IRAC bandpasses. We discuss these measurements in the context of the broadband colours and brightness temperatures of the hot-Jupiter atmosphere distribution. Specifically, we re-examine the proposed temperature dichotomy between the most irradiated, and mildly irradiated planets. We find no evidence for multiple clusters in the brightness temperature - equilibrium temperature distributions in any of these bandpasses, suggesting a continuous distribution of heat re-emission and circulation characteristics for these planets.
We describe the third telescope of the Wise Observatory, a 0.70-m Centurion
28 (C28IL) installed in 2013 and named the Jay Baum Rich telescope to enhance
significantly the wide-field imaging possibilities of the observatory. The
telescope operates from a 5.5-m diameter dome and is equipped with a
large-format red-sensitive CCD camera, offering a ~one square degree imaged
field sampled at 0".83/pixel. The telescope was acquired to provide an
alternative to the existing 1-m telescope for studies such as microlensing,
photometry of transiting exo-planets, the follow-up of supernovae and other
optical transients, and the detection of very low surface brightness extended
features around galaxies.
The operation of the C28IL is robotic, requiring only the creation of a night
observing plan that is loaded in the afternoon prior to the observations. The
entire facility was erected for a component and infrastructure cost of well
under 300k$ and a labor investment of about two person-year. The successful
implementation of the C28IL, at a reasonable cost, demonstrates the viability
of small telescopes in an age of huge light-collectors.
The polarization of cosmic microwave background (CMB) can be used to constrain cosmological birefringence, the rotation of the linear polarization of CMB photons potentially induced by parity violating physics beyond the standard model. This effect produces non-null CMB cross correlations between temperature and B mode-polarization, and between E- and B-mode polarization. Both cross-correlations are otherwise null in the standard cosmological model. We use the recently released 2015 Planck likelihood in combination with the Bicep2/Keck/Planck (BKP) likelihood to constrain the birefringence angle $\alpha$. Our findings, that are compatible with no detection, read $\alpha = 0.0^{\circ} \pm 1.3^{\circ} \mbox{ (stat)} \pm 1^{\circ} \mbox{ (sys)} $ for {\sc Planck} data and $\alpha = 0.30^{\circ} \pm 0.27^{\circ} \mbox{ (stat)} \pm 1^{\circ} \mbox{(sys)} $ for BKP data. We finally forecast the expected improvements over present constraints when the Planck BB, TB and EB spectra at high $\ell$ will be included in the analysis.
The Ku-band (13.8 GHz - 2.2 cm) RADAR instrument onboard the Cassini-Huygens spacecraft has revealed the richness of the surface of Titan, as numerous seas, lakes, rivers, cryo-volcanic flows and vast dune fields have been discovered. Linear dunes are a major geomorphological feature present on Titan, covering up to 17% of its surface, mainly in equatorial regions. However, the resolution of the RADAR instrument is not good enough to allow a detailed study of the morphology of these features. In addition, other linear wind-related landforms, such as mega-yardangs (linear wind-abraded ridges formed in cohesive rocks), are likely to present a comparable radar signature that could be confused with the one of dunes. We conducted a comparative study of the radar radiometry of both linear dunes and mega-yardangs, based on representative terrestrial analogues: the linear dunes located in the Great Sand Sea in western Egypt and in the Namib Desert in Namibia, and the mega-yardangs observed in the Lut Desert in eastern Iran and in the Borkou Desert in northern Chad. We analysed the radar scattering of both terrestrial linear dunes and mega-yardangs, using high-resolution radar images acquired by the X-band (9.6 GHz - 3.1 cm) sensor of the TerraSAR-X satellite. Variations seen in the radar response of dunes are the result of a contrast between the dune and interdune scattering, while for mega-yardangs these variations are the result of a contrast between ridges and erosion valleys. We tested a simple surface scattering model, with parameters derived from the local topography and surface roughness estimates, to accurately reproduce the radar signal variations for both landforms. It appears that we can discriminate between two types of dunes - bare interdunes as in Egypt and sand-covered interdunes as in Namibia, and between two types of mega-yardangs - young yardangs...
Observations from multiple spacecraft show that there are energy spectral "breaks" at 1-10MeV in some large CME-driven shocks. However, numerical models can hardly simulate this property due to high computational expense. The present paper focuses on analyzing these energy spectral "breaks" by Monte Carlo particle simulations of an isolated CME-driven shock. Taking the Dec 14 2006 CME-driven shock as an example, we investigate the formation of this energy spectral property. For this purpose, we apply different values for the scattering time in our isolated shock model to obtain the highest energy "tails", which can potentially exceed the "break" energy range. However, we have not found the highest energy "tails" beyond the "break" energy range, but instead find that the highest energy "tails" reach saturation near the range of energy at 5MeV. So, we believe that there exists an energy spectral "cut off" in an isolated shock. If there is no interaction with another shock, there would not be formation of the energy spectral "break" property.
We measure the binarity of detached M-dwarfs in the Kepler field with orbital periods in the range of 1-90 days. Kepler's photometric precision and nearly continuous monitoring of stellar targets over time baselines ranging from 3 months to 4 years make its detection efficiency for eclipsing binaries nearly complete over this period range and for all radius ratios. Our investigation employs a statistical framework akin to that used for inferring planetary occurrence rates from planetary transits. The obvious simplification is that eclipsing binaries have a vastly improved detection efficiency that is limited chiefly by their geometric probabilities to eclipse. For the M-dwarf sample observed by the Kepler Mission, the fractional incidence of eclipsing binaries implies that there are $0.11 ^{+0.02} _{-0.04}$ close stellar companions per apparently single M-dwarf. Our measured binarity is higher than previous inferences of the occurrence rate of close binaries via radial velocity techniques, at roughly the 2$\sigma$ level. This study represents the first use of eclipsing binary detections from a high quality transiting planet mission to infer binary statistics. Application of this statistical framework to the eclipsing binaries discovered by future transit surveys will establish better constraints on short-period M$+$M binary rate, as well as binarity measurements for stars of other spectral types.
Here we report our first spectral observations of Be/X-ray and gamma-ray binaries obtained with the new Echelle spectrograph of the National Astronomical Observatory Rozhen. For four objects (LSI+61303, gamma Cas, MWC 148, 4U 2206+54), we report the parameters and estimate the sizes of their circumstellar discs using different emission lines (H-alpha, H-beta, H-gamma, HeI and FeII). For MWC 148, we find that the compact object goes deeply through the disc. The flank inflections of H-alpha can be connected with inner ring formed at the periastron passage or radiation transfer effects. We point out an intriguing similarity between the optical emission lines of the $\gamma$-ray binary MWC 148 and the well known Be star $\gamma$ Cas.
We present a multi-wavelength analysis of the infrared dust bubble S24, and its environs, with the aim of investigating the characteristics of the molecular gas and the interstellar dust linked to them, and analyzing the evolutionary status of the young stellar objects (YSOs) identified there. Using APEX data, we mapped the molecular emission in the CO(2-1), $^{13}$CO(2-1), C$^{18}$O(2-1), and $^{13}$CO(3-2) lines in a region of about 5'x 5' in size around the bubble. The cold dust distribution was analyzed using ATLASGAL and Herschel images. Complementary IR and radio data were also used.The molecular gas linked to the S24 bubble, G341.220-0.213, and G341.217-0.237 has velocities between -48.0 km sec$^{-1}$ and -40.0 km sec$^{-1}$. The gas distribution reveals a shell-like molecular structure of $\sim$0.8 pc in radius bordering the bubble. A cold dust counterpart of the shell is detected in the LABOCA and Herschel images.The presence of extended emission at 24 $\mu$m and radio continuum emission inside the bubble indicates that the bubble is a compact HII region. Part of the molecular gas bordering S24 coincides with the extended infrared dust cloud SDC341.194-0.221. A cold molecular clump is present at the interface between S24 and G341.217-0.237. As regards G341.220-0.213, the presence of an arc-like molecular structure at the northern and eastern sections of this IR source indicates that G341.220-0.213 is interacting with the molecular gas. Several YSO candidates are found to be linked to the IR extended sources, thus confirming their nature as active star-forming regions. The total gas mass in the region and the H$_2$ ambient density amount to 10300 M$_{\odot}$ and 5900 cm$^{-3}$, indicating that G341.220-0.213, G341.217-0.237, and the S24 HII region are evolving in a high density medium. A triggering star formation scenario is also investigated.
We report on a comparison between space weather events that occurred around the two peaks in the sunspot number (SSN) during solar cycle 24. The two SSN peaks occurred in the years 2012 and 2014. Even though SSN was larger during the second peak, we find that there were more space weather events during the first peak. The space weather events we considered are large solar energetic particle (SEP) events and major geomagnetic storms associated with coronal mass ejections (CMEs). We also considered interplanetary type II radio bursts, which are indicative of energetic CMEs driving shocks. When we compared the CME properties between the two SSN peaks, we find that more energetic CMEs occurred during the 2012 peak. In particular, we find that CMEs accompanying IP type II bursts had an average speed of 1543 km/s during the 2012 peak compared to 1201 km/s during the 2014 peak. This result is consistent with the reduction in the average speed of the general population of CMEs during the second peak. All SEP events were associated with the interplanetary type II bursts, which are better than halo CMEs as indicators of space weather. The comparison between the two peaks also revealed the discordant behavior CME rate and SSN is more pronounced during the second peak. None of the 14 disk-center halo CMEs was associated with a major storm in 2014. The lone major storm in 2014 was due to the intensification of the (southward) magnetic field in the associated magnetic cloud by a shock that caught up and propagated into the magnetic cloud.
Aims. To develop a fully Bayesian least squares deconvolution (LSD) that can be applied to the reliable detection of magnetic signals in noise-limited stellar spectropolarimetric observations using multiline techniques. Methods. We consider LSD under the Bayesian framework and we introduce a flexible Gaussian Process (GP) prior for the LSD profile. This prior allows the result to automatically adapt to the presence of signal. We exploit several linear algebra identities to accelerate the calculations. The final algorithm can deal with thousands of spectral lines in a few seconds. Results. We demonstrate the reliability of the method with synthetic experiments and we apply it to real spectropolarimetric observations of magnetic stars. We are able to recover the magnetic signals using a small number of spectral lines, together with the uncertainty at each velocity bin. This allows the user to consider if the detected signal is reliable. The code to compute the Bayesian LSD profile is freely available.
Geminga was first detected as a gamma-ray point source by the SAS-2 gamma-ray satellite observatory and the COS-B X-ray satellite observatory. Subsequent observations have identified Geminga as a heavily obscured radio-quiet pulsar associated with a nearby (250 pc) late Sedov phase (300,000 year) supernova remnant. The Geminga pulsar is the second brightest source detected by the Large Area Telescope aboard the Fermi gamma-ray satellite (Fermi-LAT) and has been frequently advanced as a source of the anomalous excess of cosmic ray positrons reported by PAMELA, Fermi-LAT, and AMS-2. It is surrounded by a compact X-ray pulsar wind nebula. Observations above 10 TeV by the water Cherenkov observatory Milagro have also revealed a diffuse gamma-ray halo around Geminga extending over several square degrees. Since 2007 the VERITAS IACT observatory has performed observations of Geminga and the surrounding halo region. However, the standard methods of source detection in VERITAS data have insufficient sensitivity to angularly extended sources (>0.5 degrees) to reveal a source on the scale of the Milagro detection. In this talk, we describe two approaches being developed to search for angularly extended very high energy gamma-ray emission surrounding the Geminga pulsar.
Direct techniques for cosmic ray observations have reached an unprecedented level of precision, unveiling fine-details of the energy spectra. I will introduce the evidence for new spectral features which has been accumulated by new experiments over the past few years, and review the main ideas invoked in the theoretical explanations of the revealed spectral breaks and elemental spectra non-universality. I will also briefly comment on the complementary situation of antimatter observations.
We present results obtained from a set of cosmological hydrodynamic simulations of galaxy clusters, aimed at comparing predictions with observational data on the diversity between cool-core and non-cool-core clusters. Our simulations include the effects of stellar and AGN feedback and are based on an improved version of the Smoothed-Particle-Hydrodynamics code GADGET-3, which ameliorates gas mixing and better captures gas-dynamical instabilities by including a suitable artificial thermal diffusion. In this Letter, we focus our analysis on the entropy profiles, our primary diagnostic to classify the degree of cool-coreness of clusters, and on the iron profiles. In keeping with observations, our simulated clusters display a variety of behaviors in entropy profiles: they range from steadily decreasing profiles at small radii, characteristic of cool-core systems, to nearly flat core isentropic profiles, characteristic of non cool-core systems. Using observational criteria to distinguish between the two classes of objects, we find them to occur in similar proportions in simulations and in observations. Furthermore, we also find that simulated cool-core clusters have profiles of iron abundance that are steeper than those of non-cool-core clusters, also in agreement with observational results. We show that the capability of our simulations to generate a realistic cool-core structure in the cluster population is due to AGN feedback and artificial thermal diffusion: their combined action allows to naturally distribute the energy extracted from super-massive black holes and to compensate the radiative losses of low-entropy gas with short cooling time residing in the cluster core.
Muons produced in Extensive Air Showers (EAS) generate ring-like images in Imaging Atmospheric Cherenkov Telescopes when travelling near parallel to the optical axis. From geometrical parameters of these images, the absolute amount of light emitted may be calculated analytically. Comparing the amount of light recorded in these images to expectation is a well established technique for telescope optical efficiency calibration. However, this calculation is usually performed under the assumption of an approximately circular telescope mirror. The H.E.S.S. experiment entered its second phase in 2012, with the addition of a fifth telescope with a non-circular 600m$^2$ mirror. Due to the differing mirror shape of this telescope to the original four H.E.S.S. telescopes, adaptations to the standard muon calibration were required. We present a generalised muon calibration procedure, adaptable to telescopes of differing shapes and sizes, and demonstrate its performance on the H.E.S.S. II array.
Whether we are alone in the universe is one of the greatest mysteries facing humankind. Given the >100 billion stars in our galaxy, many have argued that it is statistically unlikely that life, including intelligent life, has not emerged anywhere else. The lack of any sign of extraterrestrial intelligence, even though on a cosmic timescale extraterrestrial civilizations would have enough time to cross the galaxy, is known as Fermi's Paradox. One possible explanation for Fermi's Paradox is the Zoo Hypothesis which states that one or more extraterrestrial civilizations know of our existence and can reach us, but have chosen not to disturb us or even make their existence known to us. I propose here a proactive test of the Zoo Hypothesis. Specifically, I propose to send a message using television and radio channels to any extraterrestrial civilization(s) that might be listening and inviting them to respond. Even though I accept this is unlikely to be successful in the sense of resulting in a response from extraterrestrial intelligences, the possibility that extraterrestrial civilizations are monitoring us cannot be dismissed and my proposal is consistent with current scientific knowledge. Besides, issuing an invitation is technically feasible, cheap and safe, and few would deny the profound importance of establishing contact with one or more extraterrestrial intelligences. A website has been set up (this http URL) to encourage discussion of this proposal and for drafting the invitation message.
We study the effects of dark-matter annihilation during the epoch of big-bang nucleosynthesis on the primordial abundances of light elements. We improve the calculation of the light-element abundances by taking into account the effects of anti-nucleons emitted by the annihilation of dark matter and the interconversion reactions of neutron and proton at inelastic scatterings of energetic nucleons. Comparing the theoretical prediction of the primordial light-element abundances with the latest observational constraints, we derive upper bounds on the dark-matter pair-annihilation cross section. Implication to some of particle-physics models are also discussed.
We investigate the gauged Nambu-Jona-Lasinio model in curved spacetime at the large $N_c$ limit and in slow-roll approximation. The model can be described by the renormalization group corrected gauge-Higgs-Yukawa theory with the corresponding compositeness conditions. Evaluating the renormalization group (RG) improved effective action, we show that such model can produce CMB fluctuations and find inflationary parameters: spectral index, tensor-to-scalar-ratio and running of the spectral index. We demonstrate that the model can naturally satisfy the Planck 2015 data and maybe considered as an alternative candidate for Higgs inflation.
In the framework of the concordance cosmological model the first-order scalar and vector perturbations of the homogeneous background are derived without any supplementary approximations in addition to the weak gravitational field limit. The sources of these perturbations (inhomogeneities) are presented in the discrete form of a system of separate point-like gravitating masses. The obtained expressions for the metric corrections are valid at all (sub-horizon and super-horizon) scales and converge in all points except the locations of the sources, and their average values are zero (thus, first-order backreaction effects are absent). Both the Minkowski background limit and the Newtonian cosmological approximation are reached under certain well-defined conditions. An important feature of the velocity-independent part of the scalar perturbation is revealed: up to an additive constant it represents a sum of Yukawa potentials produced by inhomogeneities with the same finite time-dependent Yukawa interaction range. The suggesting itself connection between this range and the homogeneity scale is briefly discussed along with other possible physical implications.
Astrophysical black hole candidates are thought to be the Kerr black holes of general relativity, but there is currently no direct observational evidence that the spacetime geometry around these objects is described by the Kerr solution. The study of the properties of the electromagnetic radiation emitted by gas or stars orbiting these objects can potentially test the Kerr black hole hypothesis. In this paper, I review the state of the art of this research field, describing the possible approaches to test the Kerr metric with current and future observational facilities and discussing current constraints.
We analyse a non-relativistic approximation of the Dirac equation for slow fermions, coupled to the chameleon field and torsion in the spacetime with the Schwarzschild metric, taken in the weak gravitational field of the Earth approximation. We follow the analysis of the Dirac equation in the curved spacetime with torsion, proposed by Kostelecky (Phys. Rev. D69, 105009 (2004)), and apply the Foldy--Wouthuysen transformations. We derive the effective low-energy gravitational potentials for slow fermions, coupled to the gravitational field of the Earth, the chameleon field and to torsion with minimal and non-minimal couplings.
The Hamamatsu R11410 and R8520 photomultipliers (PMTs) are designed for applications in xenon-based dark matter detectors, featuring sensitivity to VUV light and stable operability at cryogenic temperatures. For eleven R11410 and seven R8520 PMTs, the relative photo-detection sensitivity at 470 nm was measured as a function of the position of the incident light on the photocathode. Typical non-uniformity values, considering 80 % of the photocathode surfaces, are 5 % - 10 % and 25 % - 30 % for the R11410 and R8520 types, respectively. It was found that the non-uniformity in the inner region of the photocathode is dominated by light reflections on the inner parts of the photosensors and the border regions contribute most to the observed non-uniformity.
Folding uncertainty in theoretical models into Bayesian parameter estimation is necessary in order to make reliable inferences. A general means of achieving this is by marginalising over model uncertainty using a prior distribution constructed using Gaussian process regression (GPR). Here, we apply this technique to (simulated) gravitational-wave signals from binary black holes that could be observed using advanced-era gravitational-wave detectors. Unless properly accounted for, uncertainty in the gravitational-wave templates could be the dominant source of error in studies of these systems. We explain our approach in detail and provide proofs of various features of the method, including the limiting behaviour for high signal-to-noise, where systematic model uncertainties dominate over noise errors. We find that the marginalised likelihood constructed via GPR offers a significant improvement in parameter estimation over the standard, uncorrected likelihood. We also examine the dependence of the method on the size of training set used in the GPR; on the form of covariance function adopted for the GPR, and on changes to the detector noise power spectral density.
Starting with the idea of regularization of singularities due to the variability of the fundamental constants in cosmology we first study the cyclic universe models. We find two models of oscillating mass density and pressure regularized by varying gravitational constant $G$. Then, we extend this idea onto the multiverse containing cyclic individual universes with either growing or decreasing entropy though leaving the net entropy constant. In order to get the key idea, we consider the doubleverse with the same geometrical evolution of the two "parallel" universes with their physical evolution (physical coupling constants $c(t)$ and $G(t)$) being different. An interesting point is that there is a possibility to exchange the universes at the point of maximum expansion -- the fact which was already noticed in quantum cosmology. Similar scenario is also possible within the framework of Brans-Dicke theory.
We discuss the Mass -Radius diagram for static neutron star models obtained by the numerical solution of modified Tolman-Oppenheimer-Volkoff equations in $f(R)$ gravity where the Lagrangians $f(R)=R+\alpha R^2 (1+\gamma R)$ and $f(R)=R^{1+\epsilon }$ are adopted. Unlike the case of the perturbative approach previously reported, the solutions are constrained by the presence of an extra degree of freedom, coming from the trace of the field equations. In particular, the stiffness of the equation of state determines an upper limit on the central density $\rho_c$ above which the the positivity condition of energy-matter tensor trace $T^{\rm m}=\rho - 3 p$ holds. In the case of quadratic f(R)-gravity, we find higher masses and radii at lower central densities with an inversion of the behavior around a pivoting $\rho_c$ which depends on the choice of the equation of state. When considering the cubic corrections, we find solutions converging to the required asymptotic behavior of flat metric only for $\gamma < 0$. A similar analysis is performed for $f(R)=R^{1+\epsilon }$ considering $\epsilon$ as the leading parameter. We work strictly in the Jordan frame in order to consider matter minimally coupled with respect to geometry. This fact allows us to avoid ambiguities that could emerge in adopting the Einstein frame.
We consider a modification of GR with a special type of a non-local f(R). The structure of the non-local operators is motivated by the string field theory and p-adic string theory. We pay special account to the stability of the de Sitter solution in our model and formulate the conditions on the model parameters to have a stable configuration. Relevance of unstable configurations for the description of the de Sitter phase during inflation is discussed. Special physically interesting values of parameters are studied in details.
Links to: arXiv, form interface, find, astro-ph, recent, 1509, contact, help (Access key information)
The High Altitude Water Cherenkov (HAWC) high-energy gamma-ray observatory has recently been completed on the slopes of the Sierra Negra volcano in central Mexico. HAWC consists of 300 Water Cherenkov Detectors, each containing 180 m$^3$ of ultra-purified water, that cover a total surface area of 20,000 m$^2$. It detects and reconstructs cosmic- and gamma-ray showers in the energy range of 100 GeV to 100 TeV. The HAWC trigger for the highest energy gammas reaches an effective area of 10$^5$ m$^2$ but many of them are poorly reconstructed because the shower core falls outside the array. An upgrade that increases the present fraction of well reconstructed showers above 10 TeV by a factor of 3-4 can be done with a sparse outrigger array of small water Cherenkov detectors that pinpoint the core position and by that improve the angular resolution of the reconstructed showers. Such an outrigger array would be of the order of 200 small water Cherenkov detectors of 2.5 m$^3$ placed over an area four times larger than HAWC. Detailed simulations are being performed to optimize the layout.
We investigate the cores of fossil galaxy groups and clusters (`fossil
systems') using archival Chandra data for a sample of 17 fossil systems. We
determined the cool-core fraction for fossils via three observable diagnostics,
the central cooling time, cuspiness, and concentration parameter. We quantified
the dynamical state of the fossils by the X-ray peak/brightest cluster galaxy
(BCG), and the X-ray peak/emission weighted centre separations. We studied the
X-ray emission coincident with the BCG to detect the presence of potential
thermal coronae. A deprojection analysis was performed for z < 0.05 fossils to
obtain cooling time and entropy profiles, and to resolve subtle temperature
structures. We investigated the Lx-T relation for fossils from the 400d
catalogue to see if the scaling relation deviates from that of other groups.
Most fossils are identified as cool-core objects via at least two cool-core
diagnostics. All fossils have their dominant elliptical galaxy within 50 kpc of
the X-ray peak, and most also have the emission weighted centre within that
distance. We do not see clear indications of a X-ray corona associated with the
BCG unlike that has been observed for some other objects. Fossils do not have
universal temperature profiles, with some low-temperature objects lacking
features that are expected for ostensibly relaxed objects with a cool-core. The
entropy profiles of the z < 0.05 fossil systems can be well-described by a
power law model, albeit with indices smaller than 1. The 400d fossils Lx-T
relation shows indications of an elevated normalisation with respect to other
groups, which seems to persist even after factoring in selection effects.
Was PS1-10jh (Gezari et al. 2012), an optical/UV transient discovered by the Pan-STARRS Medium Deep Survey, the tidal disruption of a star by a massive black hole? We address two aspects of the problem: the composition of the putative disrupted object (using the spectroscopic data), and the energetics of the observed gas and radiation (using the photometric data). We perform photoionization calculations and compare with the observed lower limit of the line ratio L(He II 4686)/L(Halpha) > 5 to argue that this event was not the disruption of a solar-type star, and instead was likely the disruption of a helium core (as first proposed by Gezari et al. 2012). Disruption of such a dense object requires a relatively small central BH, M_BH <~ 2 x 10^5 M_sun. We use the photometric data to infer that PS1-10jh comprised an outflow of ~ 0.01 M_sun of gas, escaping from the BH at ~1000 km/s, and we propose that this outflow was driven primarily by radiation pressure trapped by Thomson and resonance line scattering. The large ratio of radiated energy to kinetic energy, E_rad/E_K ~ 10^4, together with the large value of E_rad ~ 2 x 10^(51) erg, suggests that the outflow was shocked at large radius (perhaps similar to super-luminous supernovae or the internal shock model for gamma-ray bursts). We describe puzzles in the physics of PS1-10jh, and discuss how this event may help us understand future tidal disruptions and super-Eddington accretion events as well.
Photo-evaporation and planet formation have both been proposed as mechanisms responsible for the creation of a transition disc. We have studied their combined effect through a suite of 2d simulations of protoplanetary discs undergoing X-ray photoevaporation with an embedded giant planet. In a previous work we explored how the formation of a giant planet triggers the dispersal of the inner disc by photo-evaporation at earlier times than what would have happened otherwise. This is particularly relevant for the observed transition discs with large holes and high mass accretion rates that cannot be explained by photo-evaporation alone. In this work we significantly expand the parameter space investigated by previous simulations. In addition, the updated model includes thermal sweeping, needed for studying the complete dispersal of the disc. After the removal of the inner disc the disc is a non accreting transition disc, an object that is rarely seen in observations. We assess the relative length of this phase, to understand if it is long lived enough to be found observationally. Depending on the parameters, especially on the X-ray luminosity of the star, we find that the fraction of time spent as a non-accretor greatly varies. We build a population synthesis model to compare with observations and find that in general thermal sweeping is not effective enough to destroy the outer disc, leaving many transition discs in a relatively long lived phase with a gas free hole, at odds with observations. We discuss the implications for transition disc evolution. In particular, we highlight the current lack of explanation for the missing non-accreting transition discs with large holes, which is a serious issue in the planet hypothesis.
We report the discovery of a stellar stream in the Dark Energy Survey (DES) Year 1 (Y1A1) data. The discovery was made through simple color-magnitude filters and visual inspection of the Y1A1 data. We refer to this new object as the Phoenix stream, after its residing constellation. Through the subtraction of the background stellar population we detect a clear signal of a simple stellar population. By fitting the ridge line of the stream in color-magnitude space, we find that a stellar population with age $\tau=11.5\pm0.5$ Gyr and ${\rm [Fe/H]}<-1.6$ located 17.5$\pm$0.9 kpc from the Sun gives an adequate description of the stream stellar population. The stream is detected over an extension of 8$^{\circ}$.1 (2.5 kpc) and has a width of $\sim$54 pc assuming a Gaussian profile, indicating that a globular cluster is a probable progenitor. There is no known globular cluster within 5 kpc compatible with being the progenitor of the stream, assuming that the stream traces its orbit. We examined overdensities along the stream, however no obvious counterpart bound stellar system is visible in the coadded images. We also find overdensities along the stream that appear to be symmetrically distributed - consistent with the epicyclic overdensity scenario for the formation of cold streams - as well as a misalignment between the Northern and Southern part of stream. We find evidence that this stream and the halo cluster NGC 1261 might have a common accretion origin linked to the recently found EriPhe overdensity (Li et al. in preparation).
We present the Rhapsody-G suite of cosmological hydrodynamic AMR zoom simulations of ten massive galaxy clusters at the $M_{\rm vir}\sim10^{15}\,{\rm M}_\odot$ scale. These simulations include cooling and sub-resolution models for star formation and stellar and supermassive black hole feedback. The sample is selected to capture the whole gamut of assembly histories that produce clusters of similar final mass. We present an overview of the successes and shortcomings of such simulations in reproducing both the stellar properties of galaxies as well as properties of the hot plasma in clusters. In our simulations, a long-lived cool-core/non-cool core dichotomy arises naturally, and the emergence of non-cool cores is related to low angular momentum major mergers. Nevertheless, the cool-core clusters exhibit a low central entropy compared to observations, which cannot be alleviated by thermal AGN feedback. For cluster scaling relations we find that the simulations match well the $M_{500}-Y_{500}$ scaling of Planck SZ clusters but deviate somewhat from the observed X-ray luminosity and temperature scaling relations in the sense of being slightly too bright and too cool at fixed mass, respectively. Stars are produced at an efficiency consistent with abundance matching constraints and central galaxies have star formation rates consistent with recent observations. While our simulations thus match various key properties remarkably well, we conclude that the shortcomings strongly suggest an important role for non-thermal processes (through feedback or otherwise) or thermal conduction in shaping the intra-cluster medium.
We present an investigation into the effects of survey systematics such as varying depth, point spread function (PSF) size, and extinction on the galaxy selection and correlation in photometric, multi-epoch, wide area surveys. We take the Canada-France-Hawaii Telescope Lensing Survey (CFHTLenS) as an example. Variations in galaxy selection due to systematics are found to cause density fluctuations of up to 10% for some small fraction of the area for most galaxy redshift slices and as much as 50% for some extreme cases of faint high-redshift samples. This results in correlations of galaxies against survey systematics of order $\sim$1% when averaged over the survey area. We present an empirical method for mitigating these systematic correlations from measurements of angular correlation functions using weighted random points. These weighted random catalogs are estimated from the observed galaxy over densities by mapping these to survey parameters. We are able to model and mitigate the effect of systematic correlations allowing for non-linear dependencies of density on systematics. Applied to CFHTLenS we find that the method reduces spurious correlations in the data by a factor two for most galaxy samples and as much as an order of magnitude in others. Such a treatment is particularly important for an unbiased estimation of very small correlation signals, as e.g. from weak gravitational lensing magnification bias. We impose a criterion for using a galaxy sample in a magnification measurement of the majority of the systematic correlations show improvement and are less than 10% of the expected magnification signal when combined in the galaxy cross correlation. After correction the galaxy samples in CFHTLenS satisfy this criterion for $z_{\rm phot}<0.9$ and will be used in a future analysis of magnification.
We have examined a subset of 18 active galactic nuclei (AGNs) drawn from a sample of 81 galaxies that possess double-peaked narrow optical emission line spectra in the Sloan Digital Sky Survey, have two optical AGN emission components separated by >0.2", and are detected in the Faint Images of the Radio Sky at Twenty-centimeters survey. Without follow-up observations, the sources of the double-peaked narrow emission lines are uncertain, and may be produced by kpc-scale separation dual active supermassive black holes, AGN outflows, or disk rotation. In this work, we propose a new methodology to characterize double-peaked narrow emission-line galaxies based on optical long-slit spectroscopy and high resolution multi-band Very Large Array observations. The nature of the radio emission in the sample galaxies is varied. Of the 18 galaxies, we detect two compact flat-spectrum radio cores with projected spatial separations on the sky between 0.6-1.6 kpc in three galaxies: J1023+3243, J1158+3231, and J1623+0808. The two radio sources are spatially coincident with the two optical components of ionized gas with AGN-like line ratios, which confirms the presence of dual AGNs in these three galaxies. Dual AGNs account for only ~15% (3/18) of the double-peaked AGNs in our sample. Gas kinematics produce ~75% (13/18) of the double-peaked narrow emission lines, distributed in the following way: 7 AGN wind-driven outflows, 5 radio-jet driven outflows, and one rotating narrow-line region. The remaining ~10% (2/18) are ambiguous cases. Our method demonstrates the power of spatially resolved spectroscopy and high resolution radio observations for the identification of AGN outflows and AGN pairs with angular separations as small as 0.18".
Measurements of the redshift-space galaxy clustering have been a prolific source of cosmological information in recent years. In the era of precision cosmology, accurate covariance estimates are an essential step for the validation of galaxy clustering models of the redshift-space two-point statistics. For cases where only a limited set of simulations is available, assessing the data covariance is not possible or only leads to a noisy estimate. Also, relying on simulated realisations of the survey data means that tests of the cosmology dependence of the covariance are expensive. With these two points in mind, this work aims at presenting a simple theoretical model for the linear covariance of anisotropic galaxy clustering observations with synthetic catalogues. Considering the Legendre moments (`multipoles') of the two-point statistics and projections into wide bins of the line-of-sight parameter (`clustering wedges'), we describe the modelling of the covariance for these anisotropic clustering measurements for galaxy samples with a trivial geometry in the case of a Gaussian approximation of the clustering likelihood. The explicit formulas are presented for Fourier space and for configuration space covariance matrices. To validate our model, we create synthetic HOD galaxy catalogues by populating the haloes of an ensemble of large-volume N-body simulations. Using linear and non-linear input power spectra, we find excellent agreement between the model predictions and the measurements on the synthetic catalogues.
We report the discovery of an excess of main sequence turn-off stars in the direction of the constellations of Eridanus and Phoenix from the first year data of the Dark Energy Survey (DES). The Eridanus-Phoenix (EriPhe) overdensity is centered around l~285 deg and b~-60 deg and spans at least 30 deg in longitude and 10 deg in latitude. The Poisson significance of the detection is at least 9 sigma. The stellar population in the overdense region is similar in brightness and color to that of the nearby globular cluster NGC 1261, indicating that the heliocentric distance of EriPhe is about d~16 kpc. The extent of EriPhe in projection is therefore at least ~4 kpc by ~3 kpc. On the sky, this overdensity is located between NGC 1261 and a new stellar stream discovered by DES at a similar heliocentric distance, the so-called Phoenix Stream. Given their similar distance and proximity to each other, it is possible that these three structures may be kinematically associated. Alternatively, the EriPhe overdensity is morphologically similar to the Virgo overdensity and the Hercules-Aquila cloud, which also lie at a similar Galactocentric distance. These three overdensities lie along a polar plane separated by ~120 deg and may share a common origin. Spectroscopic follow-up observations of the stars in EriPhe are required to fully understand the nature of this overdensity.
Using a sample of ~100 nearby line-emitting galaxy nuclei, we have built the currently definitive atlas of spectroscopic measurements of H_alpha and neighboring emission lines at subarcsecond scales. We employ these data in a quantitative comparison of the nebular emission in Hubble Space Telescope (HST) and ground-based apertures, which offer an order-of-magnitude difference in contrast, and provide new statistical constraints on the degree to which Transition Objects and low-ionization nuclear emission-line regions (LINERs) are powered by an accreting black hole at <10 pc. We show that while the small-aperture observations clearly resolve the nebular emission, the aperture dependence in the line ratios is generally weak, and this can be explained by gradients in the density of the line-emitting gas: the higher densities in the more nuclear regions potentially flatten the excitation gradients, suppressing the forbidden emission. The Transition Objects show a threefold increase in the incidence of broad H_alpha emission in the high-resolution data, as well as the strongest density gradients, supporting the composite model for these systems as accreting sources surrounded by star-forming activity. The narrow-line LINERs appear to be the weaker counterparts of the Type 1 LINERs, where the low accretion rates cause the disappearance of the broad-line component. The enhanced sensitivity of the HST observations reveals a 30% increase in the incidence of accretion-powered systems at z~0. A comparison of the strength of the broad-line emission detected at different epochs implies potential broad-line variability on a decade-long timescale, with at least a factor of three in amplitude.
Because most large galaxies contain a central black hole, and galaxies often merge, black-hole binaries are expected to be common in galactic nuclei. Although they cannot be imaged, periodicities in the light curves of quasars have been interpreted as evidence for binaries, most recently in PG~1302-102, with a short rest-frame optical period of 4 years. If the orbital period matches this value, then for the range of estimated black hole masses the components would be separated by 0.007-0.017 pc, implying relativistic orbital speeds. There has been much debate over whether black hole orbits could be smaller than 1 pc. Here we show that the amplitude and the sinusoid-like shape of the variability of PG~1302-102 can be fit by relativistic Doppler boosting of emission from a compact, steadily accreting, unequal-mass binary. We predict that brightness variations in the ultraviolet light curve track those in the optical, but with a 2-3 times larger amplitude. This prediction is relatively insensitive to the details of the emission process, and is consistent with archival UV data. Follow-up UV and optical observations in the next few years can test this prediction and confirm the existence of a binary black hole in the relativistic regime.
We present estimates for the size and the logarithmic slope of the disk temperature profile of the lensed quasar Q 2237+0305 independent of the component velocities. These estimates are based on 6 epochs of multi-wavelength narrow band images from the Nordic Optical Telescope. For each pair of lensed images and for each photometric band, we determine the microlensing amplitude and chromaticity using pre-existing mid-IR photometry to define the baseline for no microlensing magnification. A statistical comparison of the combined microlensing data (6 epochs $\times$ 5 narrow bands $\times$ 6 image pairs) with simulations based on microlensing magnification maps gives Bayesian estimates for the half-light radius of $R_{1/2}=8.3^{+11.8}_{-4.8}\sqrt{ \langle M \rangle/0.3\, M_\odot}$ light-days, and $p=0.7\pm0.3$ for the logarithmic temperature profile $T\propto R^{ -1/p}$ exponent. This size estimate is in good agreement with most recent studies. Other works based on the study of single microlensing events predict smaller sizes, but could be statistically biased by focussing on high magnification events.
We use cosmological simulations to identify dark matter subhalo host candidates of the Fornax dwarf spheroidal galaxy using the stellar kinematic properties of Fornax. We consider cold dark matter (CDM), warm dark matter (WDM), and decaying dark matter (DDM) simulations for our models of structure formation. The subhalo candidates in CDM typically have smaller mass and higher concentrations at z = 0 than the corresponding candidates in WDM and DDM. We examine the formation histories of the ~ 100 Fornax candidate subhalos identified in CDM simulations and, using approximate luminosity-mass relationships for subhalos, we find two of these subhalos that are consistent with both the Fornax luminosity and kinematics. These two subhalos have a peak mass over ten times larger than their z = 0 mass. We suggest that in CDM the dark matter halo hosting Fornax must have been severely stripped of mass and that it had an infall time into the Milky Way of ~ 9 Gyr ago. In WDM, we find that candidate subhalos consistent with the properties of Fornax have a similar infall time and a similar degree of mass loss, while in DDM we find a later infall time of ~ 3 - 4 Gyr ago and significantly less mass loss. We discuss these results in the context of the Fornax star formation history, and show that these predicted subhalo infall times can be linked to different star formation quenching mechanisms. This emphasizes the links between the properties of the dark matter and the mechanisms that drive galaxy evolution.
The young open cluster Dolidze 25, in the direction of the Galactic Anticentre, has been attributed a very low metallicity, with typical abundances between $-0.5$ and $-0.7$ dex below solar. We intend to derive accurate cluster parameters and accurate stellar abundances for some of its members. We have obtained a large sample of intermediate- and high-resolution spectra for stars in and around Dolidze 25. We used the FASTWIND code to generate stellar atmosphere models to fit the observed spectra. We derive stellar parameters for a large number of OB stars in the area, and abundances of oxygen and silicon for a number of stars with spectral types around B0. We measure low abundances in stars of Dolidze 25. For the three stars with spectral types around B0, we find $0.3$ dex (Si) and $0.5$ dex (O) below the values typical in the solar neighbourhood. These values, even though not as low as those given previously, confirm Dolidze 25 and the surrounding H II region Sh2-284 as the most metal-poor star-forming environment known in the Milky Way. We derive a distance $4.5\pm0.3\:$kpc to the cluster ($r_{\textrm{G}}\approx12.3\:$kpc). The cluster cannot be older than $\sim3\:$Myr, and likely is not much younger. One star in its immediate vicinity, sharing the same distance, has Si and O abundances at most $0.15\:$dex below solar. The low abundances measured in Dolidze 25 are compatible with currently accepted values for the slope of the Galactic metallicity gradient, if we take into account that variations of at least $\pm0.15\:$dex are observed at a given radius. The area traditionally identified as Dolidze 25 is only a small part of a much larger star-forming region that comprises the whole dust shell associated with Sh2-284 and very likely several other smaller H II regions in its vicinity.
Near-infrared (NIR) spectra that have an angular resolution of ~ 0.15 arcsec are used to examine the stellar content of the central regions of the nearby elliptical galaxy Maffei 1. The spectra were recorded at the Subaru Telescope, with wavefront distortions corrected by the RAVEN Multi-Object Adaptive Optics science demonstrator. The Ballick-Ramsey C_2 absorption bandhead near 1.76 microns is detected, and models in which 10 - 20% of the light near 1.8 microns originates from stars of spectral type C5 reproduce this feature. Archival NIR and mid-infrared images are also used to probe the structural and photometric properties of the galaxy. Comparisons with models suggest that an intermediate age population dominates the spectral energy distribution between 1 and 5 microns near the galaxy center. This is consistent not only with the presence of C stars, but also with the large HBeta index that has been measured previously for Maffei 1. The J-K color is more-or-less constant within 15 arcsec of the galaxy center, suggesting that the brightest red stars are well-mixed in this area.
The discrepancy between the amplitudes of matter fluctuations inferred from
Sunyaev-Zel'dovich (SZ) cluster number counts, the primary temperature, and the
polarization anisotropies of the cosmic microwave background (CMB) measured by
the Planck satellite can be reconciled if the local universe is embedded in an
under-dense region as shown by Lee, 2014. Here using a simple void model
assuming the open Friedmann-Robertson-Walker geometry and a Markov Chain Monte
Carlo technique, we investigate how deep the local under-dense region needs to
be to resolve this discrepancy. Such local void, if exists, predicts the local
Hubble parameter value that is different from the global Hubble constant. We
derive the posterior distribution of the local Hubble parameter from a joint
fitting of the Planck CMB data and SZ cluster number counts assuming the simple
void model. We show that the predicted local Hubble parameter value of $H_{\rm
loc}=70.1\pm0.34~{\rm km\,s^{-1}Mpc^{-1}}$ is in better agreement with direct
local Hubble parameter measurements, indicating that the local void model
provides a consistent solution to the cluster number counts and Hubble
parameter discrepancies.
We report the discovery of RG1M0150, a massive, recently quenched galaxy at z=2.636 that is multiply imaged by the cluster MACSJ0150.3-1005. We derive a stellar mass of log M_*=11.49+0.10-0.16 and a half-light radius of R_e,maj =1.8+-0.4 kpc. Taking advantage of the lensing magnification, we are able to spatially resolve a remarkably massive yet compact quiescent galaxy at z>2 in ground-based near-infrared spectroscopic observations using Magellan/FIRE and Keck/MOSFIRE. We find no gradient in the strength of the Balmer absorption lines over 0.6 R_e - 1.6 R_e, which are consistent with an age of 760 Myr. Gas emission in [NII] broadly traces the spatial distribution of the stars and is coupled with weak Halpha emission (log [NII]/Halpha = 0.6+-0.2), indicating that OB stars are not the primary ionizing source. The velocity dispersion within the effective radius is sigma_e = 271+-41 km/s. We detect rotation in the stellar absorption lines for the first time beyond z~1. Using a two-integral Jeans model that accounts for observational effects, we measure a dynamical mass of log M_dyn =11.25+-0.15 and V/sigma=0.70+-0.22. This is a high degree of rotation considering the modest observed ellipticity of 0.12+-0.08, but it is consistent with predictions from dissipational merger simulations that produce compact remnants. The mass of RG1M0150 implies that it is likely to become a slowly rotating elliptical. If it is typical, this suggests that the progenitors of massive ellipticals retain significant net angular momentum after quenching which later declines, perhaps through accretion of satellites.
Using time-series photometry from the Kilodegree Extremely Little Telescope (KELT) exoplanet survey, we are looking for eclipses of stars by their protoplanetary disks, specifically in young stellar associations. To date, we have discovered two previously unknown, large dimming events around the young stars RW Aurigae and V409 Tau. We attribute the dimming of RW Aurigae to an occultation by its tidally disrupted disk, with the disruption perhaps resulting from a recent flyby of its binary companion. Even with the dynamical environment of RW Aurigae, the distorted disk material remains very compact and presumably capable of forming planets. This system also shows that strong binary interactions with disks can also influence planet and core composition by stirring up and mixing materials during planet formation. We interpret the dimming of V409 Tau to be due to a feature, possibly a warp or perturbation, lying at least 10 AU from the host star in its nearly edge-on circumstellar disk.
We have searched for presence of current star formation in outer stellar rings of early-type disk (S0-Sb) galaxies by inspecting a representative sample of nearby galaxies with rings from the recent Spitzer catalog ARRAKIS (Comeron et al. 2014). We have found that regular rings (of R-type) reveal young stellar population with the age of less than 200~Myr in about half of all the cases, while in the pseudorings (open rings, R'), which inhabit only spiral galaxies, current star formation proceeds almost always.
Recent analyses of cosmic ray arrival directions have resulted in evidence for a positive correlation with active galactic nuclei source positions that has weak significance against an isotropic source distribution. In this paper, we explore the sample size needed to measure a highly statistically significant correlation to a parent source catalogue. We compare several scenarios for the directional scattering of ultra-high energy cosmic rays given our current knowledge of the galactic and intergalactic magnetic fields. We find significant correlations are possible for a sample of >1000 cosmic ray protons.
In this study, we investigate the dispersive properties of smoothed particle magnetohydrodynamics (SPM) in a strongly magnetized medium by using linear analysis. In modern SPM, a correction term proportional to the divergence of the magnetic fields is subtracted from the equation of motion to avoid a numerical instability arising in a strongly magnetized medium. From the linear analysis, it is found that SPM with the correction term suffer from significant dispersive errors, especially for slow waves propagating along magnetic fields. The phase velocity for all wave numbers is significantly larger than the exact solution and has a peak at a finite wavenumber. These excessively large dispersive errors occur because magnetic fields contribute an unphysical repulsive force along magnetic fields. The dispersive errors cannot be reduced, even with a larger smoothing length and smoother kernel functions such as the Gaussian or quintic spline kernels. We perform the linear analysis for this problem and find that the dispersive errors can be removed completely while keeping SPM stable if the correction term is reduced by half. These findings are confirmed by several simple numerical experiments.
In this article I will review the theory behind the gravitational wave driven r-mode instability in rapidly rotating neutron stars and discuss which constraints can be derived from observations of spins and temperatures in Low Mass X-ray Binaries. I will discuss how a standard, `minimal' neutron star model is not consistent with the data, and discuss some of the additional physical mechanisms that could reconcile theory with observations. In particular I will focus on additional forms of damping due to exotic cores and on strong mutual friction due to superfluid vortices cutting through superconducting flux tubes, and examine the repercussions these effects could have on the saturation amplitude of the mode. Finally I will also discuss the possibility that oscillations due to r-modes may have been recently observed in the X-ray light curves of two Low Mass X-ray Binaries.
We present mid-infrared (MIR) luminosity functions (LFs) of local star-forming (SF) galaxies in the AKARI NEP-Wide Survey field. In order to derive more accurate luminosity function, we used spectroscopic sample only. Based on the NEP-Wide point source catalogue containing a large number of infrared (IR) sources distributed over the wide (5.4 sq. deg.) field, we incorporated the spectroscopic redshift data for about 1790 selected targets obtained by optical follow-up surveys with MMT/Hectospec and WIYN/Hydra. The AKARI continuous 2 to 24 micron wavelength coverage as well as photometric data from optical u band to NIR H-band with the spectroscopic redshifts for our sample galaxies enable us to derive accurate spectral energy distributions (SEDs) in the mid-infrared. We carried out SED fit analysis and employed 1/Vmax method to derive the MIR (8, 12, and 15 micron rest-frame) luminosity functions. We fit our 8 micron LFs to the double power-law with the power index of alpha= 1.53 and beta= 2.85 at the break luminosity. We made extensive comparisons with various MIR LFs from several literatures. Our results for local galaxies from the NEP region are generally consistent with other works for different fields over wide luminosity ranges. The comparisons with the results from the NEP-Deep data as well as other LFs imply the luminosity evolution from higher redshifts towards the present epoch.
The X-ray Timing and Polarization (XTP) is a mission concept for a future space borne X-ray observatory and is currently selected for early phase study. We present a new design of X-ray polarimeter based on the time projection gas chamber. The polarimeter, placed above the focal plane, has an additional rear window that allows hard X-rays to penetrate (a transmission of nearly 80% at 6 keV) through it and reach the detector on the focal plane. Such a design is to compensate the low detection efficiency of gas detectors, at a low cost of sensitivity, and can maximize the science return of multilayer hard X-ray telescopes without the risk of moving focal plane instruments. The sensitivity in terms of minimum detectable polarization, based on current instrument configuration, is expected to be 3% for a 1mCrab source given an observing time of 10^5 s. We present preliminary test results, including photoelectron tracks and modulation curves, using a test chamber and polarized X-ray sources in the lab.
The Lightweight Asymmetry and Magnetism Probe (LAMP) is a micro-satellite mission concept dedicated for astronomical X-ray polarimetry and is currently under early phase study. It consists of segmented paraboloidal multilayer mirrors with a collecting area of about 1300 cm^2 to reflect and focus 250 eV X-rays, which will be detected by position sensitive detectors at the focal plane. The primary targets of LAMP include the thermal emission from the surface of pulsars and synchrotron emission produced by relativistic jets in blazars. With the expected sensitivity, it will allow us to detect polarization or place a tight upper limit for about 10 pulsars and 20 blazars. In addition to measuring magnetic structures in these objects, LAMP will also enable us to discover bare quark stars if they exist, whose thermal emission is expected to be zero polarized, while the thermal emission from neutron stars is believed to be highly polarized due to plasma polarization and the quantum electrodynamics (QED) effect. Here we present an overview of the mission concept, its science objectives and simulated observational results.
We present an analysis of the nuclear infrared (IR, 1.6 to 18 $\mu$m) emission of the ultraluminous IR galaxy UGC 5101 to derive the properties of its active galactic nucleus (AGN) and its obscuring material. We use new mid-IR high angular resolution ($0.3-0.5$ arcsec) imaging using the Si-2 filter ($\lambda_{C}=8.7\, \mu$m) and $7.5-13$ $\mu$m spectroscopy taken with CanariCam (CC) on the 10.4m Gran Telescopio CANARIAS. We also use archival HST/NICMOS and Subaru/COMICS imaging and Spitzer/IRS spectroscopy. We estimate the near- and mid-IR unresolved nuclear emission by modelling the imaging data with GALFIT. We decompose the Spitzer/IRS and CC spectra using a power-law component, which represents the emission due to dust heated by the AGN, and a starburst component, both affected by foreground extinction. We model the resulting unresolved near- and mid-IR, and the starburst subtracted CC spectrum with the CLUMPY torus models of Nenkova et al. The derived geometrical properties of the torus, including the large covering factor and the high foreground extinction needed to reproduce the deep $9.7\, \mu$m silicate feature, are consistent with the lack of strong AGN signatures in the optical. We derive an AGN bolometric luminosity $L_{bol}\sim1.9\times10^{45}\,$erg s$^{-1}$ that is in good agreement with other estimates in the literature.
In this paper, we present a study on persistent and gradual penumbral decay and correlated decline of the photospheric transverse field component during 10-20 hours before a major flare (X1.8) eruption on 2011 September 7. This long-term pre-eruption behavior is corroborated with the well-imaged pre-flare filament rising, the consistent expansion of coronal arcades overlying the filament, as well as the NLFFF modelling results in the literature. We suggest that both the long-term pre-flare penumbral decay and the transverse field decline are the photospheric manifestation of the gradual rise of the coronal filament-flux rope system. We also suggest that a C3 flare and subsequent reconnection process preceding the X1.8 flare play an important role in triggering the later major eruption.
We present a study of spectral properties of galaxies in underdense large-scale structures, voids. Our void galaxy sample (75,939 galaxies) is selected from the Sloan Digital Sky Survey (SDSS) Data Release 7 (DR7) with $\rm z < 0.107$. We find that there are no significant differences in the luminosities, stellar masses, stellar populations, and specific star formation rates between void galaxies of specific spectral types and their wall counterparts. However, the fraction of star-forming galaxies in voids is significantly higher ($\ge 9\%$) than that in walls. Void galaxies, when considering all spectral types, are slightly fainter, less massive, have younger stellar populations and of higher specific star formation rates than wall galaxies. These minor differences are totally caused by the higher fraction of star-forming galaxies in voids. We confirm that AGNs exist in voids, already found by \cite{co08}, with similar abundance as in walls. Type I AGNs contribute $\sim$ 1\%-2\% of void galaxies, similar to their fraction in walls. The intrinsic [O III] luminosities , spanning over $10^6\ L_{\sun} \sim 10^9\ L_{\sun}$, and Eddington ratios are similar comparing our void AGNs versus wall AGNs. Small scale statistics show that all spectral types of void galaxies are less clustered than their counterparts in walls. Major merger may not be the dominant trigger of black hole accretion in the luminosity range we probe. Our study implies that the growth of black holes relies weakly on large scale structures.
The material used to form the CEMP-no stars presents signatures of material processed by the CNO cycle and by He-burning from a previous stellar generation called the source stars. In order to reproduce the relative abundance ratios like for instance C/N or $^{12}$C/$^{13}$C, some mixing between the two burning regions must have occured in the source stars and only the outer layers of the stars, with modest amount coming from the CO core, must have been expelled either through stellar winds or at the time of the (faint) supernova event. With new models at low metallicity including rotational mixing, we shall discuss how the variety of abundances observed for CEMP-no stars can be reproduced.
Observations at Very High Energies (VHE, E > 100 GeV) of the BL Lac object Mrk 501 taken with the High Energy Stereoscopic System (H.E.S.S.) in four distinct periods between 2004 and 2014 are presented, with focus on the 2012 and 2014 flaring states. The source is detected with high significance above $\sim$ 2 TeV in $\sim$ 13.1 h livetime. The observations comprise low flux states and strong flaring events, which in 2014 show a flux level comparable to the 1997 historical maximum. Such high flux states enable spectral variability and flux variability studies down to a timescale of four minutes in the 2-20 TeV energy range. During the 2014 flare, the source is clearly detected in each of these bins. The intrinsic spectrum is well described by a power law of index $\Gamma=2.15\pm0.06$ and does not show curvature in this energy range. Flux dependent spectral analyses show a clear harder-when-brighter behaviour. The high flux levels and the high sensitivity of H.E.S.S. allow studies in the unprecedented combination of short timescales and an energy coverage that extends significantly above 20 TeV. The high energies allow us to probe the effect of EBL absorption at low redshifts, jet physics and LIV. The multiwavelength context of these VHE observations is presented as well.
We analyse simulations of turbulent, magnetised molecular cloud cores focussing on the formation of Class 0 stage protostellar discs and the physical conditions in their surroundings. We show that for a wide range of initial conditions Keplerian discs are formed in the Class 0 stage already. Furthermore, we show that the accretion of mass and angular momentum in the surroundings of protostellar discs occurs in a highly anisotropic manner, by means of a few narrow accretion channels. The magnetic field structure in the vicinity of the discs is highly disordered, revealing field reversals up to distances of 1000 AU. These findings demonstrate that as soon as even mild turbulent motions are included, the classical disc formation scenario of a coherently rotating environment and a well-ordered magnetic field breaks down.
The high-frequency peaked BL Lac object 1ES 0229+200 (z = 0.14) was first detected in very high energy (VHE, E > 100 GeV) $\gamma$-rays by the H.E.S.S. (High Energy Stereoscopic System) collaboration in 2006. No flux variability was reported in the initial study and its spectral characteristics have been used to derive constraints on the extragalactic background light (EBL) and on the intergalactic magnetic field (IGMF). 1ES 0229+200 has been observed with H.E.S.S. for $\sim$130 hours from 2004 to 2013: the full dataset analysed with a more sensitive method will be presented here. The results indicate that the source is not constant and displays flux variability on yearly and monthly timescales. The existence of flux variability affects the derivation of the constraints on the IGMF. The H.E.S.S. observations cover several simultaneous multi-frequency campaigns and the VHE variations are compared with those reported in different bands.
The spectral shape of the prompt emissions of gamma-ray bursts (GRBs) is typically expressed by the Band function: smooth joining of two power-law functions for high-energy and low-energy regions. To reveal the origin of the Band function, we revisit the stochastic acceleration model, in which electrons are accelerated via scattering with turbulent waves in the GRB outflow. The balance between the acceleration and synchrotron cooling yields a narrow energy-distribution similar to the Maxwellian distribution. The synchrotron spectrum becomes consistent with the observed hard photon index for the low-energy region. On the other hand, the narrow electron energy distribution contradicts the power-law spectrum for the high-energy region. We consider an evolution of the electron energy distribution to solve this problem. The turbulence and magnetic field induced by a certain hydrodynamical instability gradually decay. According to this evolution, the typical synchrotron photon energy also decreases with time. The time-integrated spectrum forms the power-law shape for the high-energy region. We discuss the required evolutions of the turbulence and magnetic field to produce a typical Band function. Although the decay of the turbulence is highly uncertain, recent numerical simulations for decaying turbulence seem comparatively positive for the stochastic acceleration model. Another condition required to reconcile observations is a much shorter duration of the stochastic acceleration than the dynamical timescale.
We present the updated orbital solution for a transient Be X-ray binary V 0332+53 currently undergoing a major outburst using the data from gamma-ray burst monitor onboard the Fermi mission. We model the observed changes in the spin-frequency of the pulsar and deduce the orbital parameters of the system significantly improving existing constrains. The intrinsic spin-up of the neutron star is found to be comparable to values observed in previous giant outburst.
We describe the observing strategy, data reduction tools and early results of
a supernova (SN) search project, named SUDARE, conducted with the ESO VST
telescope aimed at measuring the rate of the different types of SNe in the
redshift range 0.2<z<0.8. The search was performed in two of the best-studied
extragalactic fields, CDFS and COSMOS, for which a wealth of ancillary data are
available in the literature or public archives.
(abridged)
We obtained a final sample of 117 SNe, most of which are SNIa (57%) and the
remaining core collapse events of which 44% type II, 22% type IIn and 34% type
Ib/c. In order to link the transients, we built a catalog of ~1.3x10^5 galaxies
in the redshift range 0<z<1 with a limiting magnitude K_AB=23.5 mag. We
measured the SN rate per unit volume for SN Ia and core collapse SNe in
different bin of redshifts. The values are consistent with other measurements
from the literature. The dispersion of the rate measurements for SNe Ia is
comparable with the scatter of the theoretical tracks for single (SD) and
double degenerate (DD) binary systems models, therefore the data do not allow
to disentangle among the two different progenitor scenarios. However, we may
notice that among the three tested models, SD and two flavours of DD, either
with a steep (DDC) or a wide (DDW) delay time distribution, the SD gives a
better fit across the whole redshift range whereas the DDC better matches the
steep rise up to redshift ~1.2. The DDW appears instead less favoured. The core
collapse SN rate is fully consistent, unlike recent claims, with the prediction
based on recent estimates of the star formation history, and standard
progenitor mass range.
The thin-shell instability has been named as one process, which can generate entangled structures in astrophysical plasma on collisional (fluid) scales. It is driven by a spatially varying imbalance between the ram pressure of the inflowing upstream plasma and the downstream's thermal pressure at a non-planar shock. Here we show by means of a particle-in-cell (PIC) simulation that an analogue process can destabilize a thin shell formed by two interpenetrating, unmagnetized and collisionless plasma clouds. The amplitude of the shell's spatial modulation grows and saturates after about ten inverse proton plasma frequencies, when the shell consists of connected piecewise linear patches.
Context: One of the most prominent processes suggested to heat the corona to
well above 10^6 K builds on nanoflares, short bursts of energy dissipation.
Aims: We compare observations to model predictions to test the validity of
the nanoflare process.
Methods: Using extreme UV data from AIA/SDO and HMI/SDO line-of-sight
magnetograms we study the spatial and temporal evolution of a set of loops in
active region AR 11850.
Results: We find a transient brightening of loops in emission from Fe xviii
forming at about 7.2 MK while at the same time these loops dim in emission from
lower temperatures. This points to a fast heating of the loop that goes along
with evaporation of material that we observe as apparent upward motions in the
image sequence. After this initial phases lasting for some 10 min, the loops
brighten in a sequence of AIA channels showing cooler and cooler plasma,
indicating the cooling of the loops over a time scale of about one hour. A
comparison to the predictions from a 1D loop model shows that this observation
supports the nanoflare process in (almost) all aspects. In addition, our
observations show that the loops get broader while getting brighter, which
cannot be understood in a 1D model.
Radio synchrotron emission is a powerful tool to study the strength and
structure of magnetic fields in galaxies. Unpolarized synchrotron emission
traces isotropic turbulent fields which are strongest in spiral arms and bars
(20-30\mu G) and in central starburst regions (50-100\mu G). Such fields are
dynamically important; they affect gas flows and drive gas inflows in central
regions.
Polarized emission traces ordered fields, which can be regular or anisotropic
turbulent, where the latter originates from isotropic turbulent fields by the
action of compression or shear. The strongest ordered fields (10-15\mu G) are
generally found in interarm regions. In galaxies with strong density waves,
ordered fields are also observed at the inner edges of spiral arms. Ordered
fields with spiral patterns exist in grand-design, barred and flocculent
galaxies, and in central regions. Ordered fields in interacting galaxies have
asymmetric distributions and are a tracer of past interactions between galaxies
or with the intergalactic medium.
Faraday rotation measures of the diffuse polarized radio emission from galaxy
disks reveal large-scale spiral patterns that can be described by the
superposition of azimuthal modes; these are signatures of regular fields
generated by mean-field dynamos. "Magnetic arms" between gaseous spiral arms
may also be products of dynamo action, but need a stable spiral pattern to
develop. Helically twisted field loops winding around spiral arms were found in
two galaxies so far. Large-scale field reversals, like the one found in the
Milky Way, could not yet be detected in external galaxies.
In radio halos around edge-on galaxies, ordered magnetic fields with X-shaped
patterns are observed. Halo fields of even symmetry probably dominate.
The origin and evolution of cosmic magnetic fields will be studied with
forthcoming radio telescopes like the Square Kilometre Array.
A joint analysis is done of the radio and X-ray observations of SN 1993J. It is argued that neither synchrotron cooling behind the forward shock nor thermal cooling behind the reverse shock is supported by observations. In order for adiabatic models to be consistent, a reinterpretation of the radius of the spatially resolved VLBI-source is needed during the first few hundred days. Instead of reflecting the position of the forward shock, it is then associated with the expansion of the Rayleigh-Taylor unstable region emanating from the contact discontinuity. Although observations imply a constant ratio between the energy densities in magnetic fields and relativistic electrons, they do not appear to scale individually with the thermal energy density behind the forward shock; rather, in adiabatic models, the evolution of the magnetic field strength is best understood as scaling inversely with the supernova radius.
Quasi-periodic propagating intensity disturbances have been observed in large coronal loops in EUV images over a decade, and are widely accepted to be slow magnetosonic waves. However, spectroscopic observations from Hinode/EIS revealed their association with persistent coronal upflows, making this interpretation debatable. We perform a 2.5D magnetohydrodynamic simulation to imitate the chromospheric evaporation and the following reflected patterns in a flare loop. Our model encompasses the corona, transition region, and chromosphere. We demonstrate that the quasi periodic propagating intensity variations captured by the synthesized \textit{Solar Dynamics Observatory}/Atmospheric Imaging Assembly (AIA) 131, 94~\AA~emission images match the previous observations well. With particle tracers in the simulation, we confirm that these quasi periodic propagating intensity variations consist of reflected slow mode waves and mass flows with an average speed of 310 km/s in an 80 Mm length loop with an average temperature of 9 MK. With the synthesized Doppler shift velocity and intensity maps of the \textit{Solar and Heliospheric Observatory}/Solar Ultraviolet Measurement of Emitted Radiation (SUMER) Fe XIX line emission, we confirm that these reflected slow mode waves are propagating waves.
We present a novel technique that uses the autocorrelation of the spectrum of a star to measure the line broadening caused by the modulus of its average surface magnetic field. The advantage of the autocorrelation comes from the fact that it can detect very small spectral line broadening effects because it averages over many spectral lines and therefore gives an average with a very high signal to noise ratio. We validate the technique with the spectra of known magnetic stars and obtain autocorrelation curves that are in full agreement with published magnetic curves obtained with Zeeman splitting. The autocorrelation also gives less noisy curves so that it can be used to obtain very accurate curves. We degrade the resolution of spectra of these magnetic stars to lower spectral resolutions where the Zeeman splitting is undetectable. At these resolutions, the autocorrelation still gives good quality curves, thereby showing that it can be used to measure magnetic fields in spectra where the Zeeman splitting is significantly smaller than the width of the spectral line. This would therefore allow observing magnetic fields in very faint Ap stars with low-resolution spectrographs, thereby greatly increasing the number of known magnetic stars. It also demonstrates that the autocorrelation can measure magnetic fields in rapidly rotating stars as well as weak magnetic fields that give a Zeeman splitting smaller than the intrinsic width of the spectral lines. Finally, it shows that the autocorrelation can be used to find unknown magnetic stars in low resolution spectroscopic surveys.
We investigate the role of the H_2^+ channel on H_2 molecule formation during the collapse of primordial gas clouds immersed in strong radiation fields which are assumed to have the shape of a diluted black-body spectra with temperature T_rad. Since the photodissociation rate of H_2^+ depends on its level population, we take full account of the vibrationally-resolved H_2^+ kinetics. We find that in clouds under soft but intense radiation fields with spectral temperature T_rad < 7000 K, the H_2^+ channel is the dominant H_2 formation process. On the other hand, for harder spectra with T_rad > 7000 K, the H^- channel takes over H_2^+ in the production of molecular hydrogen. We calculate the critical radiation intensity needed for supermassive star formation by direct collapse and examine its dependence on the H_2^+ level population. Under the assumption of local thermodynamic equilibrium (LTE) level population, the critical intensity is underestimated by a factor of a few for soft spectra with T_rad < 7000 K. For harder spectra, the value of the critical intensity is not affected by the level population of H_2^+. This result justifies previous estimates of the critical intensity assuming LTE populations since radiation sources like young and/or metal-poor galaxies are predicted to have rather hard spectra.
For cylindrical differentially rotating plasmas threaded with a uniform vertical magnetic field, we study large-scale magnetic field generation from finite amplitude perturbations using analytic theory and direct numerical simulations. Analytically, we impose helical fluctuations, a seed field, and a background flow and use quasi-linear theory for a single mode. The predicted large-scale field growth agrees with numerical simulations in which the magnetorotational instability (MRI) arises naturally. The vertically and azimuthally averaged toroidal field is generated by a fluctuation-induced EMF that depends on differential rotation. Given fluctuations, the method also predicts large-scale field growth for MRI-stable rotation profiles and flows with no rotation but shear.
We have carried out an ALMA survey of 15 confirmed or candidate low-mass (<0.2M$_\odot$) members of the TW Hya Association (TWA) with the goal of detecting molecular gas in the form of CO emission, as well as providing constraints on continuum emission due to cold dust. Our targets have spectral types of M4-L0 and hence represent the extreme low end of the TWA's mass function. Our ALMA survey has yielded detections of 1.3mm continuum emission around 4 systems (TWA 30B, 32, 33, & 34), suggesting the presence of cold dust grains. All continuum sources are unresolved. TWA 34 further shows 12CO(2-1) emission whose velocity structure is indicative of Keplerian rotation. Among the sample of known ~7-10 Myr-old star/disk systems, TWA 34, which lies just ~50 pc from Earth, is the lowest mass star thus far identified as harboring cold molecular gas in an orbiting disk.
Most exoplanetary systems in binary stars are of S--type, and consist of one or more planets orbiting a primary star with a wide binary stellar companion. Gravitational forcing of a single planet by a sufficiently inclined binary orbit can induce large amplitude oscillations of the planet's eccentricity and inclination through the Kozai-Lidov (KL) instability. KL cycling was invoked to explain: the large eccentricities of planetary orbits; the family of close--in hot Jupiters; and the retrograde planetary orbits in eccentric binary systems. However, several kinds of perturbations can quench the KL instability, by inducing fast periapse precessions which stabilize circular orbits of all inclinations: these could be a Jupiter--mass planet, a massive remnant disc or general relativistic precession. Indeed, mutual gravitational perturbations in multiplanet S--type systems can be strong enough to lend a certain dynamical rigidity to their orbital planes. Here we present a new and faster process that is driven by this very agent inhibiting KL cycling. Planetary perturbations enable secular oscillations of planetary eccentricities and inclinations, also called Laplace--Lagrange (LL) eigenmodes. Interactions with a remnant disc of planetesimals can make planets migrate, causing a drift of LL mode periods which can bring one or more LL modes into resonance with binary orbital motion. The results can be dramatic, ranging from excitation of large eccentricities and mutual inclinations to total disruption. Not requiring special physical or initial conditions, binary resonant driving is generic and could have profoundly altered the architecture of many S--type multiplanet systems. It can also weaken the multiplanet occurrence rate in wide binaries, and affect planet formation in close binaries.
Polarimetric surveys of the microwave sky at large angular scales are crucial in testing cosmic inflation, as inflation predicts a divergence-free $B$-mode angular power spectrum that extends to the largest scales on the sky. A promising technique for realizing such large surveys is through the use of rapid polarization modulation to mitigate variations in the atmosphere, coupling to the environment, and drifts in instrumental response. VPMs change the state of polarization by introducing a controlled, adjustable delay between orthogonal linear polarizations resulting in transformations between linear and circular polarization states. VPMs are currently being implemented in experiments designed to measure the polarization of the cosmic microwave background on large angular scales because of their capability for providing rapid, front-end polarization modulation and control over systematic errors. Despite the advantages provided by the VPM, it is important to identify and mitigate any time-varying effects that leak into the synchronously modulated component of the signal.In this paper, we consider and address the emission from a 300 K VPM on the system performance. Though instrument alignment can greatly reduce the influence of modulated VPM emission, some residual modulated signal is expected. We consider VPM emission in the presence of system misalignments and temperature variation. We use simulations of TOD to evaluate the effect of these residual errors on the power spectrum. The analysis and modeling in this paper guides experimentalists on the critical aspects of observations using VPMs as front-end modulators. By implementing the characterizations and controls as described front-end VPM modulation can be very powerful. None of the systematic errors studied fundamentally limit the detection and characterization of B-modes on large scales for a tensor-to-scalar ratio of $r=0.01$.
The known FU Ori-type young eruptive stars are exceedingly rare (a dozen or so confirmed objects) and 2MASS J06593158-0405277, with its 2014 outburst, is likely the latest addition to the family. All members have displayed just one such eruption in their recorded history, an event lasting for decades. To test the FU Ori nature of 2MASS J06593158-0405277, we have reconstructed its photometric history by measuring its brightness on Harvard photographic plates spanning the time interval 1899-1989. No previous large amplitude eruption similar to that initiated in 2014 has been found, as in bona fide FU Ori-type objects. The median value of the brightness in quiescence of 2MASS J06593158-0405277 is B=15.5, with the time interval 1935-1950 characterized by a large variability (~1 mag amplitude) that contrasts with the remarkable photometric stability displayed at later epochs. The variability during 1935-1950 can either be ascribed to some T Tau like activity of 2MASS J06593158-0405277 itself or to the also young and fainter star 2MASS J06593168-0405224 that lies 5 arcsec to the north and forms an unresolved pair at the astrometric scale of Harvard photographic plates.
Imaging and spectroscopy at (sub-)millimeter wavelengths are key frontiers in astronomy and cosmology. Large area spectral surveys with moderate spectral resolution (R=50-200) will be used to characterize large scale structure and star formation through intensity mapping surveys in emission lines such as the CO rotational transitions. Such surveys will also be used to study the SZ effect, and will detect the emission lines and continuum spectrum of individual objects. WSPEC is an instrument proposed to target these science goals. It is a channelizing spectrometer realized in rectangular waveguide, fabricated using conventional high-precision metal machining. Each spectrometer is coupled to free space with a machined feed horn, and the devices are tiled into a 2D array to fill the focal plane of the telescope. The detectors will be aluminum Lumped-Element Kinetic Inductance Detectors (LEKIDs). To target the CO lines and SZ effect, we will have bands at 135-175 GHz and 190-250 GHz, each Nyquist-sampled at R~200 resolution. Here we discuss the instrument concept and design, and successful initial testing of a WR10 (i.e. 90 GHz) prototype spectrometer. We recently tested a WR5 (180 GHz) prototype to verify that the concept works at higher frequencies, and also designed a resonant backshort structure that may further increase the optical efficiency. We are making progress towards integrating a spectrometer with a LEKID array and deploying a prototype device to a telescope for first light.
The NASA New Horizons spacecraft flies past Pluto on July 14, 2015, carrying two instruments that detect charged particles. Pluto has a tenuous, extended atmosphere that is escaping the weak gravity of the planet. The interaction of the solar wind with the escaping atmosphere of Pluto depends on solar wind conditions as well as the vertical structure of the atmosphere. We have analyzed Voyager 2 particles and fields measurements between 25 and 39 AU and present their statistical variations. We have adjusted these predictions to allow for the declining activity of the Sun and solar wind output. We summarize the range of SW conditions that can be expected at 33 AU and survey the range of scales of interaction that New Horizons might experience. Model estimates for the solar wind stand-off distance vary from approximately 7 to 1000 RP with our best estimate being around 40 RP (where we take the radius of Pluto to be RP=1184 km).
Turbulence is ubiquitous in molecular clouds (MCs), but its origin is still unclear because MCs are usually assumed to live longer than the turbulence dissipation time. It has been shown that interstellar medium (ISM) turbulence is likely driven by SN explosions, but it has never been demonstrated that SN explosions can establish and maintain a turbulent cascade inside MCs consistent with the observations. In this work, we carry out a simulation of SN-driven turbulence in a volume of (250 pc)^3, specifically designed to test if SN driving alone can be responsible for the observed turbulence inside MCs. We find that SN driving establishes velocity scaling consistent with the usual scaling laws of supersonic turbulence. This also means that previous idealized simulations of MC turbulence, driven with a random, large-scale volume force, were correctly adopted as appropriate models for MC turbulence, despite the artificial driving. We also find the same scaling laws extend to the interior of MCs, and their normalization is consistent with the observations, as shown by the velocity-size relation of the MCs selected from our simulation. Although selected in three dimensional space and without accounting for chemistry and radiative transfer effects, MCs from our simulations have properties very similar to those of observed MCs, such as velocity-size and mass-size relations, mass and size probability distributions, and magnetic field-density relation, besides the velocity scaling. Finally, we find that MC turbulence is super-Alfv\'{e}nic with respect to both the mean and rms cloud magnetic field strength. We conclude that MCs and their internal structure and dynamics are the natural result of SN-driven turbulence.
The abundance anomalies in lambda Boo stars are popularly explained by element-specific mass inflows at rates that are much greater than empirically-inferred bounds for interstellar accretion. Therefore, a lambda Boo star's thin outer envelope must derive from a companion star, planet, analogs to Kuiper Belt Objects or a circumstellar disk. Because radiation pressure on gas-phase ions might selectively allow the accretion of carbon, nitrogen, and oxygen and inhibit the inflow of elements such as iron, the source of the acquired matter need not contain dust. We propose that at least some lambda Boo stars accrete from the winds of hot Jupiters.
We have conducted ALMA cycle 2 observations in the 1.3 mm continuum and in the C18O (2-1) and SO (5_6-4_5) lines at a resolution of ~0.3" toward the Class 0 protostar B335. The 1.3 mm continuum, C18O, and SO emission all show central compact components with sizes of ~40-180 AU within more extended components. The C18O component shows signs of infalling and rotational motion. By fitting simple kinematic models to the C18O data, the protostellar mass is estimated to be 0.05 Msun. The specific angular momentum, on a 100 AU scale, is ~4.3E-5 km/s*pc. A similar specific angular momentum, ~3E-5 to 5E-5 km/s*pc, is measured on a 10 AU scale from the velocity gradient observed in the central SO component, and there is no clear sign of an infalling motion in the SO emission. By comparing the infalling and rotational motion, our ALMA results suggest that the observed rotational motion has not yet reached Keplerian velocity neither on a 100 AU nor even on a 10 AU scale. Consequently, the radius of the Keplerian disk in B335 (if present) is expected to be 1-3 AU. The expected disk radius in B335 is one to two orders of magnitude smaller than those of observed Keplerian disks around other Class 0 protostars. Based on the observed infalling and rotational motion from 0.1 pc to inner 100 AU scales, there are two possible scenarios to explain the presence of such a small Keplerian disk in B335: magnetic braking and young age. If our finding is the consequence of magnetic braking, ~50% of the angular momentum of the infalling material within a 1000 AU scale might have been removed, and the magnetic field strength on a 1000 AU scale is estimated to be ~200 uG. If it is young age, the infalling radius in B335 is estimated to be ~2700 AU, corresponding to a collapsing time scale of ~5E4 yr.
Next generation galaxy surveys demand the development of massive ensembles of galaxy mocks to model the observables and their covariances, what is computationally prohibitive using $N$-body simulations. COLA is a novel method designed to make this feasible by following an approximate dynamics but with up to 3 orders of magnitude speed-ups when compared to an exact $N$-body. In this paper we investigate the optimization of the code parameters in the compromise between computational cost and recovered accuracy in observables such as two-point clustering and halo abundance. We benchmark those observables with a state-of-the-art $N$-body run, the MICE Grand Challenge simulation (MICE-GC). We find that using 40 time steps linearly spaced since $z_i \sim 20$, and a force mesh resolution three times finer than that of the number of particles, yields a matter power spectrum within $1\%$ for $k \lesssim 1\,h {\rm Mpc}^{-1}$ and a halo mass function within $5\%$ of those in the $N$-body. In turn the halo bias is accurate within $2\%$ for $k \lesssim 0.7\,h {\rm Mpc}^{-1}$ whereas, in redshift space, the halo monopole and quadrupole are within $4\%$ for $k \lesssim 0.4\,h {\rm Mpc}^{-1}$. These results hold for a broad range in redshift ($0 < z < 1$) and for all halo mass bins investigated ($M > 10^{12.5} \, h^{-1} \, {\rm M_{\odot}}$). To bring accuracy in clustering to one percent level we study various methods that re-calibrate halo masses and/or velocities. We thus propose an optimized choice of COLA code parameters as a powerful tool to optimally exploit future galaxy surveys.
We present new radio measurements of the expansion rate of the Crab nebula's synchrotron nebula over a ~30-yr period. We find a convergence date for the radio synchrotron nebula of CE 1255 +- 27. We also re-evaluated the expansion rate of the optical line emitting filaments, and we show that the traditional estimates of their convergence dates are slightly biased. Using an un-biased Bayesian analysis, we find a convergence date for the filaments of CE 1091 +- 34 (~40 yr earlier than previous estimates). Our results show that both the synchrotron nebula and the optical line-emitting filaments have been accelerated since the explosion in CE 1054, but that the synchrotron nebula has been relatively strongly accelerated, while the optical filaments have been only slightly accelerated. The finding that the synchrotron emission expands more rapidly than the filaments supports the picture that the latter are the result of the Rayleigh-Taylor instability at the interface between the pulsar-wind nebula and the surrounding freely-expanding supernova ejecta, and rules out models where the pulsar wind bubble is interacting directly with the pre-supernova wind of the Crab's progenitor.
We present new spectroscopy and astrometry to characterize the nearby brown dwarf WISEP J180026.60+013453.1. The optical spectral type, L7.5, is in agreement with the previously reported near-infrared spectral type. The preliminary trigonometric parallax places it at a distance of $8.01 \pm 0.21$ pc, confirming that it is the fourth closest known late-L (L7-L9) dwarf. The measured luminosity, our detection of lithium, and the lack of low surface gravity indicators indicates that WISEP J180026.60+013453.1 has a mass $0.03 < M < 0.06 M_\odot$ and an age between 300 million and 1.5 billion years according to theoretical substellar evolution models. The low space motion is consistent with this young age. We have measured the rotational broadening ($v \sin i = 13.5 \pm 0.5$ km/s), and use it to estimate a maximum rotation period of 9.3 hr.
It is an intriguing possibility that dark matter (DM) could have flavor quantum numbers like the quarks. We propose and investigate a class of UV-complete models of this kind, in which the dark matter is in a scalar triplet of an SU(3) flavor symmetry, and interacts with quarks via a colored flavor-singlet fermionic mediator. Such mediators could be discovered at the LHC if their masses are $\sim 1$ TeV. We constrain the DM-mediator couplings using relic abundance, direct detection, and flavor-changing neutral-current considerations. We find that, for reasonable values of its couplings, scalar flavored DM can contribute significantly to the real and imaginary parts of the $B_s$-$\bar B_s$ mixing amplitude. We further assess the potential for such models to explain the galactic center GeV gamma-ray excess.
We study periods of elevated energetic particle intensities observed by STEREO-A when the partial pressure exerted by energetic ($\geq$83 keV) protons ($P_{EP}$) is larger than the pressure exerted by the interplanetary magnetic field ($P_{B}$). In the majority of cases, these periods are associated with the passage of interplanetary shocks. Periods when $P_{EP}$ exceeds $P_{B}$ by more than one order of magnitude are observed in the upstream region of fast interplanetary shocks where depressed magnetic field regions coincide with increases of the energetic particle intensities. When solar wind parameters are available, $P_{EP}$ also exceeds the pressure exerted by the solar wind thermal population ($P_{TH}$). Prolonged periods ($>$12 h) with both $P_{EP}$$>$$P_{B}$ and $P_{EP}$$>$$P_{TH}$ may also occur when energetic particles accelerated by an approaching shock encounter a region well-upstream of the shock characterized by low magnetic field magnitude and tenuous solar wind density. Quasi-exponential increases of the sum $P_{SUM}$=$P_{B}$+$P_{TH}$+$P_{EP}$ are observed in the immediate upstream region of the shocks regardless of individual changes in $P_{EP}$, $P_{B}$ and $P_{TH}$, indicating a coupling between $P_{EP}$ and the pressure of the background medium characterized by $P_{B}$ and $P_{TH}$. The quasi-exponential increase of $P_{SUM}$ implies a convected exponential radial gradient $\partial{P_{SUM}}/\partial{r}$$>$0 that results in an outward force applied to the plasma upstream of the shock. This force can be maintained by the mobile energetic particles streaming upstream of the shocks that, in the most intense events, drive electric currents able to generate diamagnetic cavities and depressed solar wind density regions.
In the context of the Higgs model involving gauge and Yukawa interactions with the spontaneous gauge symmetry breaking, we consider $\lambda \phi^4$ inflation with non-minimal gravitational coupling, where the Higgs field is identified as inflaton. Since the inflaton quartic coupling is very small, once quantum corrections through the gauge and Yukawa interactions are taken into account, the inflaton effective potential most likely becomes unstable. In order to avoid this problem, we need to impose stability conditions on the effective inflation potential, which lead to not only non-trivial relations amongst the particle mass spectrum of the model, but also correlations between the inflationary predictions and the mass spectrum. For concrete discussion, we investigate the minimal $B-L$ extension of the Standard Model with identification of the $B-L$ Higgs field as inflaton. The stability conditions for the inflaton effective potential fix the mass ratio amongst the $B-L$ gauge boson, the right-handed neutrinos and the inflaton. This mass ratio also correlates with the inflationary predictions. In other words, if the $B-L$ gauge boson and the right-handed neutrinos are discovered in future, their observed mass ratio provides constraints on the inflationary predictions.
This year marks the hundredth anniversary of Einstein's 1915 landmark paper "Die Feldgleichungen der Gravitation" in which the field equations of general relativity were correctly formulated for the first time, thus rendering general relativity a complete theory. Over the subsequent hundred years physicists and astronomers have struggled with uncovering the consequences and applications of these equations. This contribution, which was written as an introduction to six chapters dealing with the connection between general relativity and cosmology that will appear in the two-volume book "One Hundred Years of General Relativity: From Genesis and Empirical Foundations to Gravitational Waves, Cosmology and Quantum Gravity," endeavors to provide a historical overview of the connection between general relativity and cosmology, two areas whose development has been closely intertwined.
An ultra-sensitive opto-mechanical force sensor has been built and tested in the optics laboratory at INFN Trieste. Its application to experiments in the Dark Energy sector, such as those for Chameleon-type WISPs, is particularly attractive, as it enables a search for their direct coupling to matter. We present here the main characteristics and the absolute force calibration of the KWISP (Kinetic WISP detection) sensor. It is based on a thin Si3N4 micro-membrane placed inside a Fabry-Perot optical cavity. By monitoring the cavity characteristic frequencies it is possible to detect the tiny membrane displacements caused by an applied force. Far from the mechanical resonant frequency of the membrane, the measured force sensitivity is 5.0e-14 N/sqrt(Hz), corresponding to a displacement sensitivity of 2.5e-15 m/sqrt(Hz), while near resonance the sensitivity is 1.5e-14 N/sqrt(Hz), reaching the estimated thermal limit, or, in terms of displacement, 7.5e-16 N/sqrt(Hz). These displacement sensitivities are comparable to those that can be achieved by large interferometric gravitational wave detectors.
We study the inflation scenarios, in the framework of superstring theory, where the inflaton is an axion producing the adiabatic curvature perturbations while there exists another light axion producing the isocurvature perturbations. We discuss how the non-trivial couplings among string axions can generically arise, and calculate the consequent cross-correlations between the adiabatic and isocurvature modes through concrete examples. Based on the Planck analysis on the generally correlated isocurvature perturbations, we show that there is a preference for the existence of the correlated isocurvature modes for the axion monodromy inflation while the natural inflation disfavors such isocurvature modes.
Gravity theories beyond General Relativity typically predict dipolar gravitational emission by compact-star binaries. This emission is sourced by "sensitivity" parameters depending on the stellar compactness. We introduce a general formalism to calculate these parameters, and show that in shift-symmetric Horndeski theories stellar sensitivities and dipolar radiation vanish, provided that the binary's dynamics is perturbative (i.e. the post-Newtonian formalism is applicable) and cosmological-expansion effects can be neglected. This allows reproducing the binary-pulsar observed orbital decay.
We provide a model-independent argument indicating that for a black hole of entropy N the non-thermal deviations from Hawking radiation, per each emission time, are of order 1/N, as opposed to exp(-N). This fact abolishes the standard a priory basis for the information paradox.
We seek wormholes among rotating cylindrically symmetric configurations in general relativity. Exact wormhole solutions are presented with such sources of gravity as a massless scalar field, a cosmological constant, and a scalar field with an exponential potential. However, none of these solutions are asymptotically flat, which excludes the existence of wormhole entrances as local objects in our Universe. To overcome this difficulty, we try to build configurations with flat asymptotic regions using the cut-and-paste procedure: on both sides of the throat, a wormhole solution is matched to a properly chosen region of flat space-time at some surfaces $\Sigma_-$ and $\Sigma_+$. It is shown, however, that if the source of gravity in the throat region is a scalar field with an arbitrary potential, then one or both thin shells appearing on $\Sigma_-$ and $\Sigma_+$ inevitably violate the null energy condition. Thus, although rotating wormhole solutions are easily found without exotic matter, such matter is still necessary for obtaining asymptotic flatness.
Links to: arXiv, form interface, find, astro-ph, recent, 1509, contact, help (Access key information)
The Galactic centre hosts a crowded, dense nuclear star cluster with a half-light radius of 4 pc. Most of the stars in the Galactic centre are cool late-type stars, but there are also >100 hot early-type stars in the central parsec of the Milky Way. These stars are only 3-8 Myr old. Our knowledge of the number and distribution of early-type stars in the Galactic centre is incomplete. Only a few spectroscopic observations have been made beyond a projected distance of 0.5 pc of the Galactic centre. The distribution and kinematics of early-type stars are essential to understand the formation and growth of the nuclear star cluster. We cover the central >4pc^2 of the Galactic centre using the integral-field spectrograph KMOS. We extracted more than 1,000 spectra from individual stars and identified early-type stars based on their spectra. Our data set contains 114 bright early-type stars: 6 have narrow emission lines, 23 are Wolf-Rayet stars, 9 stars have featureless spectra, and 76 are O/B type stars. Our wide-field spectroscopic data confirm that the distribution of young stars is compact, with 90% of the young stars identified within 0.5 pc of the nucleus. We identify 24 new O/B stars primarily at large radii. We estimate photometric masses of the O/B stars and show that the total mass in the young population is >12,000M_sun. The O/B stars all appear to be bound to the Milky Way nuclear star cluster, while less than 30% belong to the clockwise rotating disk. The central concentration of the early-type stars is a strong argument that they have formed in situ. A large part of the young O/B stars is not on the disk, which either means that the early-type stars did not all form on the same disk or that the disk is dissolving rapidly. [abridged]
We revisit the effect of peculiar velocities on low-redshift type Ia supernovae. Velocities introduce an additional guaranteed source of correlations between supernova magnitudes that should be considered in all analyses of nearby supernova samples but has largely been neglected in the past. Applying a likelihood analysis to the latest compilation of nearby supernovae, we find no evidence for the presence of these correlations, although, given the significant noise, the data is also consistent with the correlations predicted for the standard LCDM model. We then consider the dipolar component of the velocity correlations - the frequently studied "bulk velocity" - and explicitly demonstrate that including the velocity correlations in the data covariance matrix is crucial for drawing correct and unambiguous conclusions about the bulk flow. In particular, current supernova data is consistent with no excess bulk flow on top of what is expected in LCDM and effectively captured by the covariance. We further clarify the nature of the apparent bulk flow that is inferred when the velocity covariance is ignored. We show that a significant fraction of this quantity is expected to be noise bias due to uncertainties in supernova magnitudes and not any physical peculiar motion.
Using Subaru/Suprime-Cam wide-field imaging and both Keck/ESI and LBT/MODS spectroscopy, we identify and characterize a compact star cluster, which we term NGC 3628-UCD1, embedded in a stellar stream around the spiral galaxy NGC 3628. The size and luminosity of UCD1 are similar to $\omega$ Cen, the most luminous Milky Way globular cluster, which has long been suspected to be the stripped remnant of an accreted dwarf galaxy. The object has a magnitude of $i=19.3$ mag (${\rm L}_{\rm i}=1.4\times10^{6}~{\rm L}_{\odot}$). UCD1 is marginally resolved in our ground-based imaging, with a half-light radius of $\sim10$ pc. We measure an integrated brightness for the stellar stream of $i=13.1$ mag, with $(g-i)=1.0$. This would correspond to an accreted dwarf galaxy with an approximate luminosity of ${\rm L}_i\sim4.1\times10^{8}~{\rm L}_{\odot}$. Spectral analysis reveals that UCD1 has an age of $6.6$ Gyr , $[\rm{Z}/\rm{H}]=-0.75$, an $[{\alpha}/\rm{Fe}]=-0.10$. We propose that UCD1 is an example of an $\omega$ Cen-like star cluster possibly forming from the nucleus of an infalling dwarf galaxy, demonstrating that at least some of the massive star cluster population may be created through tidal stripping.
The NASA Kepler mission has revolutionised time-domain astronomy and has massively expanded the number of known extrasolar planets. However, the effect of wide multiplicity on exoplanet occurrence has not been tested with this dataset. We present a sample of 401 wide multiple systems containing at least one Kepler target star. Our method uses Pan-STARRS1 and archival data to produce an accurate proper motion catalogue of the Kepler field. Combined with Pan-STARRS1 SED fits and archival proper motions for bright stars, we use a newly developed probabilistic algorithm to identify likely wide binary pairs which are not chance associations. As by-products of this we present stellar SED templates in the Pan-STARRS1 photometric system and conversions from this system to Kepler magnitudes. We find that Kepler target stars in our binary sample with separations above 6 arcseconds are no more or less likely to be identified as confirmed or candidate planet hosts than a weighted comparison sample of Kepler targets of similar brightness and spectral type. Therefore we find no evidence that binaries with projected separations greater than 3,000AU affect the occurrence rate of planets with P<300days around FGK stars.
Measurements of large-scale B-mode polarization in the cosmic microwave background (CMB) are a fundamental goal of current and future CMB experiments. However, because of the much higher instrumental sensitivity, CMB experiments will be more sensitive to any imperfect modelling of the Galactic foreground polarization in the estimation of the primordial B-mode signal. We compare the sensitivity to B-modes for different concepts of CMB satellite missions (LiteBIRD, COrE, COrE+, PRISM, EPIC, PIXIE) in the presence of Galactic foregrounds that are either correctly or incorrectly modelled. We quantify the impact on the tensor-to-scalar parameter of imperfect foreground modelling in the component separation process. Using Bayesian parametric fitting and Gibbs sampling, we perform the separation of the CMB and the Galactic foreground B-mode polarization. The resulting CMB B-mode power spectrum is used to compute the likelihood distribution of the tensor-to-scalar ratio. We focus the analysis to the very large angular scales (l<12) that can be probed only by CMB space missions, where primordial CMB B-modes dominate over spurious B-modes induced by gravitational lensing while Galactic foregrounds still show the most significant polarization intensity. We find that fitting a single modified blackbody component for the thermal dust where the real sky should account for two dust components may strongly bias the estimation of the tensor-to-scalar ratio by more than 5{\sigma}, at least for the most sensitive experiments. Neglecting in the parametric fitting model a positive curvature of the synchrotron spectral index may bias the estimate of the tensor-to-scalar ratio by 1{\sigma} to 2{\sigma}. For sensitive CMB experiments, omitting in the foreground modelling a 1% polarized spinning dust component can induce a non-negligible bias in the reconstructed tensor-to-scalar ratio.
At a projected distance of ~26 pc from Sgr A*, the Arches cluster provides insight to star formation in the extreme Galactic Center (GC) environment. Despite its importance, many key properties such as the cluster's internal structure and orbital history are not well known. We present an astrometric and photometric study of the outer region of the Arches cluster (R > 6.25") using HST WFC3IR. Using proper motions we calculate membership probabilities for stars down to F153M = 20 mag (~2.5 M_sun) over a 120" x 120" field of view, an area 144 times larger than previous astrometric studies of the cluster. We construct the radial profile of the Arches to a radius of 75" (~3 pc at 8 kpc), which can be well described by a single power law. From this profile we place a 3-sigma lower limit of 2.8 pc on the observed tidal radius, which is larger than the predicted tidal radius (1 - 2.5 pc). Evidence of mass segregation is observed throughout the cluster and no tidal tail structures are apparent along the orbital path. The absence of breaks in the profile suggests that the Arches has not likely experienced its closest approach to the GC between ~0.2 - 1 Myr ago. If accurate, this constraint indicates that the cluster is on a prograde orbit and is located front of the sky plane that intersects Sgr A*. However, further simulations of clusters in the GC potential are required to interpret the observed profile with more confidence.
Simple assumptions made regarding electron thermodynamics often limit the extent to which general relativistic magnetohydrodynamic (GRMHD) simulations can be applied to observations of low-luminosity accreting black holes. We present, implement, and test a model that self-consistently evolves an electron entropy equation and takes into account the effects of spatially varying electron heating and relativistic anisotropic thermal conduction along magnetic field lines. We neglect the back-reaction of electron pressure on the dynamics of the accretion flow. Our model is appropriate for systems accreting at $\ll 10^{-5}$ of the Eddington rate, so radiative cooling by electrons can be neglected. It can be extended to higher accretion rates in the future by including electron cooling and proton-electron Coulomb collisions. We present a suite of tests showing that our method recovers the correct solution for electron heating under a range of circumstances, including strong shocks and driven turbulence. Our initial applications to axisymmetric simulations of accreting black holes show that (1)~physically-motivated electron heating rates yield electron temperature distributions significantly different from the constant electron to proton temperature ratios assumed in previous work, with higher electron temperatures concentrated in the coronal region between the disc and the jet; (2)~electron thermal conduction significantly modifies the electron temperature in the inner regions of black hole accretion flows if the effective electron mean free path is larger than the local scale-height of the disc (at least for the initial conditions and magnetic field configurations we study). The methods developed in this work are important for producing more realistic predictions for the emission from accreting black holes such as Sagittarius A* and M87; these applications will be explored in future work.
We derive the total cold gas, atomic hydrogen, and molecular gas masses of approximately 24 000 galaxies covering four decades in stellar mass at redshifts 0.5 < z < 3.0, taken from the CANDELS survey. Our inferences are based on the inversion of a molecular hydrogen based star formation law, coupled with a prescription to separate atomic and molecular gas. We find that: 1) there is an increasing trend between the inferred cold gas (HI and H2), HI, and H2 mass and the stellar mass of galaxies down to stellar masses of 10^8 Msun already in place at z = 3; 2) the molecular fractions of cold gas increase with increasing stellar mass and look-back time; 3) there is hardly any evolution in the mean HI content of galaxies at fixed stellar mass; 4) the cold gas fraction and relative amount of molecular hydrogen in galaxies decrease at a relatively constant rate with time, independent of stellar mass; 5) there is a large population of low-stellar mass galaxies dominated by atomic gas. These galaxies are very gas rich, but only a minor fraction of their gas is molecular; 6) the ratio between star-formation rate (SFR) and inferred total cold gas mass (HI + H2) of galaxies (i.e., star-formation efficiency; SFE) increases with star-formation at fixed stellar masses. Due to its simplicity, the presented approach is valuable to assess the impact of selection biases on small samples of directly-observed gas masses and to extend scaling relations down to stellar mass ranges and redshifts that are currently difficult to probe with direct measurements of gas content.
Images obtained by the Cassini spacecraft between 2012 and 2015 reveal a periodic brightness variation in a region of Saturn's D ring that previously appeared to be rather featureless. Furthermore, the intensity and radial wavenumber of this pattern have decreased steadily with time since it was first observed. Based on analogies with similar structures elsewhere in the D ring, we propose that this structure was created by some event that disturbed the orbital motions of the ring particles, giving them finite orbital eccentricities and initially aligned pericenters. Differential orbital precession then transformed this structure into a spiral pattern in the ring's optical depth that became increasingly tightly wound over time. The observed trends in the pattern's radial wavenumber are roughly consistent with this basic model, and also indicate that the ring-disturbing event occurred in early December 2011. Similar events in 1979 may have generated the periodic patterns seen in this same region by the Voyager spacecraft. The 2011 event could have been caused by debris striking the rings, or by a disturbance in the planet's electromagnetic environment. The rapid reduction in the intensity of the brightness variations over the course of just a few years indicates that some process is either damping orbital eccentricities in this region or causing the orbital pericenters of particles with the same semi-major axis to become misaligned.
Sublimating gas molecules scatter off of the surface of an icy body in the same manner as photons. This means that for every photon-driven body force, there should be a sublimation-driven analogue that affects icy bodies. Thermal photons emitted from the surfaces of asymmetrically shaped bodies in the Solar System generate net torques that change the spin rates of these bodies over time. The long-term averaging of this torque is called the YORP effect. Here we propose a sublimation-driven analogue to the YORP effect (SYORP), in which sublimating gas molecules emitted from the surfaces of icy bodies also generate net torques on the bodies. However, sublimating molecules carry momentum away from the body at a rate ~10^4-10^5 greater than thermal photons, resulting in much greater body torques. While previous studies of sublimative torques focused on emissions from highly localized sources on the surfaces of Jupiter Family Comet nuclei, SYORP applies to non-localized emissions across the entire body, which likely dominates sublimation-drive torques on small icy chunks and Dynamically Young Comets, and can therefore be applied without high-resolution spacecraft observations of their surfaces. Instead, we repurpose the equations of the YORP effect to account for sublimation-driven torques. We show how an SYORP-driven mechanism best matches observations of the Sun-oriented lineaments (striae) of comet tails, whose formation mechanism has remained enigmatic for decades. The SYORP effect explains why striae are observed between near-perihelion and ~1 AU from the Sun for comets with perihelia less than 0.6 AU, and solves longstanding problems with moving enough material into the cometary tail to form visible striae. We show that SYORP can form striae that match those of Comet West, and produce a power-law size frequency distribution of their parent chunks with an index of -1.4 (-1.1- 2).
Future hard X-ray telescopes (e.g. SIMBOL-X and Constellation-X) will make use of hard X-ray optics with multilayer coatings, with angular resolutions comparable to the achieved ones in the soft X-rays. One of the crucial points in X-ray optics, indeed, is multilayer interfacial microroughness that causes effective area reduction and X-Ray Scattering (XRS). The latter, in particular, is responsible for image quality degradation. Interfacial smoothness deterioration in multilayer deposition processes is commonly observed as a result of substrate profile replication and intrinsic random deposition noise. For this reason, roughness growth should be carefully investigated by surface topographic analysis, X-ray reflectivity and XRS measurements. It is convenient to express the roughness evolution in terms of interface Power Spectral Densities (PSD), that are directly related to XRS and, in turn, in affecting the optic HEW (Half Energy Width). In order to interpret roughness amplification and to help us to predict the imaging performance of hard X-ray optics, we have implemented a well known kinetic continuum equation model in a IDL language program (MPES, Multilayer PSDs Evolution Simulator), allowing us the determination of characteristic growth parameters in multilayer coatings. In this paper we present some results from analysis we performed on several samples coated with hard X-ray multilayers (W/Si, Pt/C, Mo/Si) using different deposition techniques. We show also the XRS predictions resulting from the obtained modelizations, in comparison to the experimental XRS measurements performed at the energy of 8.05 keV.
Future large X-ray observatories like ATHENA will be equipped with very large optics, obtained by assembling modular optical elements, named X-ray Optical Units (XOU) based on the technology of either Silicon Pore Optics or Slumped Glass Optics. In both cases, the final quality of the modular optic (a 5 arcsec HEW requirement for ATHENA) is determined by the accuracy alignment of the XOUs within the assembly, but also by the angular resolution of the individual XOU. This is affected by the mirror shape accuracy, its surface roughness, and the mutual alignment of the mirrors within the XOU itself. Because of the large number of XOUs to be produced, quality tests need to be routinely done to select the most performing stacked blocks, to be integrated into the final optic. In addition to the usual metrology based on profile and roughness measurements, a direct measurement with a broad, parallel, collimated and uniform X- ray beam would be the most reliable test, without the need of a focal spot reconstruction as usually done in synchrotron light. To this end, we designed the BEaTriX (Beam Expander Testing X-ray facility) to be realized at INAF-OAB, devoted to the functional tests of the XOUs. A grazing incidence parabolic mirror and an asymmetrically cut crystal will produce a parallel X-ray beam broad enough to illuminate the entire aperture of the focusing elements. An X-ray camera at the focal distance from the mirrors will directly record the image. The selection of different crystals will enable to test the XOUs in the 1 - 5 keV range, included in the X-ray energy band of ATHENA (0.2-12 keV). In this paper we discuss a possible BEaTriX facility implementation. We also show a preliminary performance simulation of the optical system.
We aim to unveil the most massive central cluster black holes in the universe. We present a new search strategy which is based on a black hole mass gain sensitive 'calorimeter' and which links the innermost stellar density profile of a galaxy to the adiabatic growth of its central SMBH. In a first step we convert observationally inferred feedback powers into SMBH growth rates by using reasonable energy conversion efficiency parameters, $\epsilon$. In the main part of this paper we use these black hole growth rates, sorted in logarithmically increasing steps encompassing our whole parameter space, to conduct $N$-Body computations of brightest cluster galaxies with the newly developed MUESLI software. For the initial setup of galaxies we use core-Sersic models in order to account for SMBH scouring. We find that adiabatically driven core re-growth is significant at the highest accretion rates. As a result, the most massive black holes should be located in BCGs with less pronounced cores when compared to the predictions of empirical scaling relations which are usually calibrated in less extreme environments. For efficiency parameters $\epsilon<0.1$, BCGs in the most massive, relaxed and X-ray luminous galaxy clusters might even develop steeply rising density cusps. Finally, we discuss several promising candidates for follow up investigations, among them the nuclear black hole in the Phoenix cluster. Based on our results, it might have a mass of the order of $10^{11} M_\odot$.
We present a sample of 11 M31 Cepheids in stellar clusters, derived from the overlap of the Panchromatic Hubble Andromeda Treasury (PHAT) cluster catalog and the Pan-STARRS1 (PS1) disk Cepheid catalog. After identifying the PS1 Cepheids in the HST catalog, we calibrate the PS1 mean magnitudes using the higher resolution HST photometry, revealing up to 1 magnitude offsets due to crowding effects in the ground-based catalog. We measure ages of the clusters by performing single stellar population fits to their color-magnitude diagrams (CMDs) excluding their Cepheids. From these cluster age measurements, we derive an empirical period-age relation which agrees well with the existing literature values. By confirming this relation for M31 Cepheids, we justify its application in high-precision pointwise age estimation across M31.
We examine energy conversion from accreting pair plasma to outgoing Poynting flux by black hole rotation. Our approach is based on a two-fluid model consisting of collisionless pair plasma. The electric potential is not constant along magnetic field lines, unlike an ideal magnetohydrodynamics approximation. We show how and where longitudinal electric fields and toroidal magnetic fields are generated by the rotation, whereas they vanish everywhere for radial flow in a split monopole magnetic field in a Schwarzschild black hole. Outgoing electromagnetic power in a steady state is calculated by applying the WKB method to the perturbation equations for a small spin parameter. In our model, the luminosity has a peak in the vicinity of the black hole, but is damped toward the event horizon and infinity. The power at the peak is of the same order as that in the Blandford--Znajek process, although the physical mechanism is different.
Under standard assumptions including stationary and serially uncorrelated Gaussian gravitational wave stochastic background signal and noise distributions, as well as homogenous detector sensitivities, the standard cross-correlation detection statistic is known to be optimal in the sense of minimizing the probability of a false dismissal at a fixed value of the probability of a false alarm. The focus of this paper is to analyze the comparative efficiency of this statistic, versus a simple alternative statistic obtained by cross-correlating the \textit{squared} measurements, in situations that deviate from such standard assumptions. We find that differences in detector sensitivities have a large impact on the comparative efficiency of the cross-correlation detection statistic, which is dominated by the alternative statistic when these differences reach one order of magnitude. This effect holds even when both the signal and noise distributions are Gaussian. While the presence of non-Gaussian signals has no material impact for reasonable parameter values, the relative inefficiency of the cross-correlation statistic is less prominent for fat-tailed noise distributions but it is magnified in case noise distributions have skewness parameters of opposite signs. Our results suggest that introducing an alternative detection statistic can lead to noticeable sensitivity gains when noise distributions are possibly non-Gaussian and/or when detector sensitivities exhibit substantial differences, a situation that is expected to hold in joint detections from Advanced LIGO and Advanced Virgo, in particular in the early phases of development of the detectors, or in joint detections from Advanced LIGO and Einstein Telescope.
A recent Chandra observation of the nearby galaxy cluster Abell 585 has led to the discovery of an extended X-ray jet associated with the high-redshift background quasar B3 0727+409, a luminous radio source at redshift z=2.5. This is one of only few examples of high-redshift X-ray jets known to date. It has a clear extension of about 10-12", corresponding to a projected length of 80-100 kpc, with a possible hot spot as far as 35" from the quasar. The archival high resolution VLA maps surprisingly reveal no extended jet emission, except for one knot about 1.4" away from the quasar. The high X-ray to radio luminosity ratio for this source appears consistent with the $\propto (1+z)^{4}$ amplification expected from the inverse Compton radiative model. This serendipitous discovery may signal the existence of an entire population of similar systems with bright X-ray and faint radio jets at high redshift, a selection bias which must be accounted for when drawing any conclusions about the redshift evolution of jet properties and indeed about the cosmological evolution of supermassive black holes and active galactic nuclei in general.
We numerically investigate the dynamics of a 2D non-magnetised protoplanetary disc surrounded by an inflow coming from an external envelope. We find that the accretion shock between the disc and the inflow is unstable, leading to the generation of large-amplitude spiral density waves. These spiral waves propagate over long distances, down to radii at least ten times smaller than the accretion shock radius. We measure spiral-driven outward angular momentum transport with 1e-4 < alpha < 1e-2 for an inflow accretion rate Mout>1e-8 Msun/yr. We conclude that the interaction of the disc with its envelope leads to long-lived spiral density waves and radial angular momentum transport with rates that cannot be neglected in young non-magnetised protostellar discs.
This Letter investigates the spectral characteristics of the interplanetary magnetic field fluctuations at proton scales during several time intervals chosen along the speed profile of a fast stream. The character of the fluctuations within the first frequency decade, beyond the high frequency break located between the fluid and the kinetic regime, strongly depends on the type of wind. While the fast wind shows a clear signature of both right handed and left handed polarized fluctuations, possibly associated with KAW and Ion-Cyclotron waves, respectively, the rarefaction region, where the wind speed and the Alfv\'{e}nicity of low frequency fluctuations decrease, shows a rapid disappearance of the ion-cyclotron signature followed by a more gradual disappearance of the KAWs. Moreover, also the power associated to perpendicular and parallel fluctuations experiences a rapid depletion, keeping, however, the power anisotropy in favour of the perpendicular spectrum.
Tracing molecular hydrogen content with carbon monoxide in low-metallicity galaxies has been exceedingly difficult. Here we present a new effort, with IRAM 30-m observations of 12CO(1-0) of a sample of 8 dwarf galaxies having oxygen abundances ranging from 12+logO/H=7.7 to 8.4. CO emission is detected in all galaxies, including the most metal-poor galaxy of our sample (0.1 Zsun); to our knowledge this is the largest number of 12CO(1-0) detections ever reported for galaxies with 12+logO/H<=8 (0.2 Zsun) outside the Local Group. We calculate stellar masses (Mstar) and star-formation rates (SFRs), and analyze our results by combining our observations with galaxy samples from the literature. Extending previous results for a correlation of the molecular gas depletion time, tau(dep), with Mstar and specific SFR (sSFR), we find a variation in tau(dep) of a factor of 200 or more (from <50 Myr to 10 Gyr) over a spread of 1000 in sSFR and Mstar. We exploit the variation of tau(dep) to constrain the CO-to-H2 mass conversion factor alpha(CO) at low metallicity, and assuming a power-law variation find alpha(CO) \propto (Z/Zsun)^1.9, similar to results based on dust continuum measurements compared with gas mass. By including HI measurements, we show that the fraction of total gas mass relative to the baryonic mass is higher in galaxies that are metal poor, of low mass, and of high sSFR. Finally, comparisons of the data with star-formation models of the molecular gas phases suggest that, at metallicities Z/Zsun<=0.2, there are some discrepancies with model predictions.
When emission in a conical relativistic jet ceases abruptly, the observed decay light curve is controlled by the high-latitude "curvature effect". If the zero time is defined properly, the decay slope and the spectral index has a simple relation \alpha=2+\beta if the relativistic jet moves with a constant Lorentz factor. Uhm & Zhang recently found that the decay is steeper than this standard value if the jet is undergoing bulk acceleration when the emission ceases. By applying this theory to the flare data of GRBs, they found that the decay properties of flares demand that the emission region is undergoing significant bulk acceleration. This suggests that the jet is PFD, and that emission is powered by significant dissipation of Poynting-flux energy within the jet. Uhm & Zhang presented three X-ray flares as the examples. In this paper, we systematically analyze the flare data released by Swift to investigate whether bulk acceleration is common among flares. We select a sample of 85 bright flares detected in 63 GRBs, and investigate the relationship between \alpha and \beta during the steep decay phase of these flares. The decay index depends on the choice of the zero time point t0. We adopt two methods. "Method I" takes t0I as the first rising data point of each flare, and is the most conservative approach. We find that 71/85 flares have decay slopes steeper than the simplest curvature effect prediction, and therefore, are in the acceleration regime. "Method II" extrapolates the rising light curve of each flare backwards until the flux density is three orders of magnitude lower than the peak flux density, and defines the corresponding time as the time zero point (t0II). We find that 81/85 flares fall into the acceleration regime. This suggests that bulk acceleration is ubiquitous among X-ray flares, pointing towards a ubiquitous Poynting-flux-dominated jet composition for these events.
The observed velocities of the gas in barred galaxies are a combination of the azimuthally-averaged circular velocity and non-circular motions, primarily caused by gas streaming along the bar. These non-circular flows must be accounted for before the observed velocities can be used in mass modeling. In this work, we examine the performance of the tilted-ring method and the DiskFit algorithm for transforming velocity maps of barred spiral galaxies into rotation curves (RCs) using simulated data. We find that the tilted-ring method, which does not account for streaming motions, under/over-estimates the circular motions when the bar is parallel/perpendicular to the projected major axis. DiskFit, which does include streaming motions, is limited to orientations where the bar is not-aligned with either the major or minor axis of the image. Therefore, we propose a method of correcting RCs based on numerical simulations of galaxies. We correct the RC derived from the tilted-ring method based on a numerical simulation of a galaxy with similar properties and projections as the observed galaxy. Using observations of NGC 3319, which has a bar aligned with the major axis, as a test case, we show that the inferred mass models from the uncorrected and corrected RCs are significantly different. These results show the importance of correcting for the non-circular motions and demonstrate that new methods of accounting for these motions are necessary as current methods fail for specific bar alignments.
A major flaring state of the BL Lac object Mrk 501 was observed by the High Energy Stereoscopic System (H.E.S.S.) in June, 2014. Flux levels higher than one Crab unit were recorded and rapid variability at very high energies ($\sim$2-20 TeV) was revealed. The high statistics afforded by the flares allowed us to probe the presence of minutes timescale variability and study its statistical characteristics exclusively at TeV energies owing to the high energy threshold of approximately 2 TeV. Doubling times of a few minutes are estimated for fluxes greater than 2 TeV. Statistical tests on the light curves show interesting temporal structure in the variations including deviations from a normal flux distribution similar to those found in the PKS 2155-304 flare of July 2006, at nearly an order of magnitude higher threshold energy. Rapid variations at such high energies put strong constraints on the physical mechanisms in the blazar jet.
Both gamma-ray bursts (GRBs) and blazars have relativistic jets pointing at a small angle from our line of sight. Several recent studies suggested that these two kinds of sources may share similar jet physics. In this work, we explore the variability properties for GRBs and blazars as a whole. We find that the correlation between minimum variability timescale (MTS) and Lorentz factor, $\Gamma$, as found only in GRBs by Sonbas et al. can be extended to blazars with a joint correlation of $\rm MTS\propto\Gamma^{-4.7\pm0.3}$. The same applies to the $\rm MTS\propto \it L_{\gamma}^{\rm -1.0\pm0.1}$ correlation as found in GRBs, which can be well extended into blazars as well. These results provide further evidence that the jets in these two kinds of sources are similar despite of the very different mass scale of their central engines. Further investigations of the physical origin of these correlations are needed, which can shed light on the nature of the jet physics.
Most stellar remnants so far have been found in binary systems, where they
interact with matter from their companions. Isolated neutron stars and black
holes are difficult to find as they are dark, yet they are predicted to exist
in our Galaxy in vast numbers.
We explored the OGLE-III database of 150 million objects observed in years
2001-2009 and found 59 microlensing events exhibiting a parallax effect due to
the Earth's motion around the Sun. Combining parallax and brightness
measurements from microlensing light curves with expected proper motions in the
Milky Way, we identified 15 microlensing events which are consistent with
having a white dwarf, neutron star or a black hole lens and we estimated their
masses and distances. The most massive of our black hole candidates has 8.3
M_Sun and is at a distance of 2.4 kpc. The distribution of masses of our
candidates indicates a continuum in mass distribution with no mass gap between
neutron stars and black holes. We also present predictions on how such events
will be observed by the astrometric Gaia mission.
ESPRESSO is an extremely stable high-resolution spectrograph which is currently being developed for the ESO VLT. With its groundbreaking characteristics it is aimed to be a "science machine", i.e., a fully-integrated instrument to directly extract science information from the observations. In particular, ESPRESSO will be the first ESO instrument to be equipped with a dedicated tool for the analysis of data, the Data Analysis Software (DAS), consisting in a number of recipes to analyze both stellar and quasar spectra. Through the new ESO Reflex GUI, the DAS (which will implement new algorithms to analyze quasar spectra) is aimed to get over the shortcomings of the existing software providing multiple iteration modes and full interactivity with the data.
Deep spectral-line surveys in the mm and sub-mm range can detect thousands of lines per band uncovering the rich chemistry of molecular clouds, star forming regions and circumstellar envelopes, among others objects. The ability to study the faintest features of spectroscopic observation is, nevertheless, limited by a number of factors. The most important are the source complexity (line density), limited spectral resolution and insufficient sideband (image) rejection (SRR). Dual Sideband (2SB) millimeter receivers separate upper and lower sideband rejecting the unwanted image by about 15 dB, but they are difficult to build and, until now, only feasible up to about 500 GHz (equivalent to ALMA Band 8). For example ALMA Bands 9 (602-720 GHz) and 10 (787-950 GHz) are currently DSB receivers. Aims: This article reports the implementation of an ALMA Band 9 2SB prototype receiver that makes use of a new technique called calibrated digital sideband separation. The new method promises to ease the manufacturing of 2SB receivers, dramatically increase sideband rejection and allow 2SB instruments at the high frequencies currently covered only by Double Sideband (DSB) or bolometric detectors. Methods: We made use of a Field Programmable Gate Array (FPGA) and fast Analog to Digital Converters (ADCs) to measure and calibrate the receiver's front end phase and amplitude imbalances to achieve sideband separation beyond the possibilities of purely analog receivers. The technique could in principle allow the operation of 2SB receivers even when only imbalanced front ends can be built, particularly at very high frequencies. Results: This digital 2SB receiver shows an average sideband rejection of 45.9 dB while small portions of the band drop below 40 dB. The performance is 27 dB (a factor of 500) better than the average performance of the proof-of-concept Band 9 purely-analog 2SB prototype receiver.
We study the intermediate evolution model and show that, compared with the recent study of a power-law evolution, the intermediate evolution is a better description of the low-redshift regime supported by observations from type Ia supernovae and BAO. We found also that recent data suggest that the intermediate evolution is as good a fit to this redshift range as the $\Lambda$CDM model.
We present a detailed analysis of the Horizontal Branch of the Carina Dwarf Spheroidal Galaxy by means of synthetic modelling techniques, taking consistently into account the star formation history and metallicity evolution as determined from main sequence and red giant branch spectroscopic observations. We found that a range of integrated red giant branch mass loss values of 0.1-0.14 M, increasing with metallicity, is able to reproduce the colour extension of the old Horizontal Branch. However, leaving the mass loss as the only free parameter is not enough to match the detailed morphology of Carina Horizontal Branch. We explored the role played by the star formation history on the discrepancies between synthetic and observed Horizontal Branches. We derived a toy bursty star formation history that reproduces the horizontal branch star counts, and also matches qualitatively the red giant and the turn off regions. This star formation history is made of a subset of age and [M/H] components of the star formation history based on turn off and red giants only, and entails four separate bursts of star formation of different strenghts, centred at 2, 5, 8.6 and 11.5 Gyr, with mean [M/H] decreasing from \sim -1.7 to \sim -2.2 for increasing ages, and a Gaussian spread of 0.1 dex. The comparison between the metallicity distribution function of our star formation history and the one measured from the infrared CaT feature using a CaT-[Fe/H] calibration shows a qualitative agreement, once taken into account the range of [Ca/Fe] abundances measured in a sample of Carina stars, that biases the derived [Fe/H] distribution toward too low values. In conclusion, we have shown how the information contained within the horizontal branch of Carina (and dwarf galaxies in general) can be extracted and interpreted to refine the star formation history derived from red giants and turn off stars only. Abridged
Recent studies of the distribution and kinematics of the Milky Way and Andromeda satellite galaxy systems have confirmed the existence of coplanar, corotating structures of galaxies. In addition to the 'missing satellite problem', these structures pose a major challenge to the standard $\Lambda$CDM scenario of structure formation. We complement the efforts made by the dwarf galaxy community to extend these studies to other nearby galaxy groups by systematically searching for faint, unresolved dwarf members with a low surface brightness in the Southern Centaurus group of galaxies. The aim is to determine whether these coplanar, corotating structures are a universal phenomenon. We imaged an area of 60 square degrees (0.3 Mpc$^2$) around the M83 subgroup with the wide-field Dark Energy Camera (DECam) at the CTIO 4 m Blanco telescope in $g$ and $r$ down to a limiting surface brightness of $\mu_r\approx 30$ mag arcsec$^{-2}$. Various image-filtering techniques were applied to the DECam data to enhance the visibility of extremely low-surface brightness objects. We report the discovery of 16 new dwarf galaxy candidates in the direction of the M83 subgroup, roughly doubling the number of known dwarfs in that region. The photometric properties of the candidates, when compared to those of the Local Group, suggest membership in the M83 subgroup. The faintest objects have a central star density of $\approx1.3 L_\odot$ pc$^{-2}$ and a total magnitude of $g = 20.25$, corresponding to $M_g = -9.55$ at the nominal distance of 4.9 Mpc. The sky distribution of the new objects is significantly prolonged toward Cen A, suggesting that many of them belong to the Cen A subgroup or a common halo. We also provide updated surface photometry for the brighter, known dwarf members in the surveyed area.
We report on simultaneous broadband observations of the TeV-emitting blazar Markarian 501 between 1 April and 10 August 2013, including the first detailed characterization of the synchrotron peak with Swift and NuSTAR. During the campaign, the nearby BL Lac object was observed in both a quiescent and an elevated state. The broadband campaign includes observations with NuSTAR, MAGIC, VERITAS, the Fermi Large Area Telescope (LAT), Swift X-ray Telescope and UV Optical Telescope, various ground-based optical instruments, including the GASP-WEBT program, as well as radio observations by OVRO, Mets\"ahovi and the F-Gamma consortium. Some of the MAGIC observations were affected by a sand layer from the Saharan desert, and had to be corrected using event-by-event corrections derived with a LIDAR (LIght Detection And Ranging) facility. This is the first time that LIDAR information is used to produce a physics result with Cherenkov Telescope data taken during adverse atmospheric conditions, and hence sets a precedent for the current and future ground-based gamma-ray instruments. The NuSTAR instrument provides unprecedented sensitivity in hard X-rays, showing the source to display a spectral energy distribution between 3 and 79 keV consistent with a log-parabolic spectrum and hard X-ray variability on hour timescales. None (of the four extended NuSTAR observations) shows evidence of the onset of inverse-Compton emission at hard X-ray energies. We apply a single-zone equilibrium synchrotron self-Compton model to five simultaneous broadband spectral energy distributions. We find that the synchrotron self-Compton model can reproduce the observed broadband states through a decrease in the magnetic field strength coinciding with an increase in the luminosity and hardness of the relativistic leptons responsible for the high-energy emission.
We present a direct determination of the stellar metallicity in the close pair galaxy NGC~4038 (D = 20 Mpc) based on the quantitative analysis of moderate resolution KMOS/VLT spectra of three super star clusters (SSCs). The method adopted in our analysis has been developed and optimised to measure accurate metallicities from atomic lines in the $J$-band of single red supergiant (RSG) or RSG-dominated star clusters. Hence, our metallicity measurements are not affected by the biases and poorly understood systematics inherent to strong line HII methods which are routinely applied to massive data sets of galaxies. We find [Z]= +0.07 $\pm$ 0.03 and compare our measurements to HII strong line calibrations. Our abundances and literature data suggest the presence of a flat metallicity gradient, which can be explained as redistribution of metal-rich gas following the strong interaction.
Precise spectra of 3C 279 in the 0.5-70 keV range, obtained during two epochs
of
\emph{Swift} and \emph{NuSTAR} observations, are analyzed using a
near-equipartition model. We apply a one-zone leptonic model with a
three-parameter log-parabola electron energy distribution (EED) to fit the
\emph{Swift} and \emph{NuSTAR} X-ray data, as well as simultaneous optical and
\emph{Fermi}-LAT $\gamma$-ray data. The Markov Chain Monte Carlo (MCMC)
technique is used to search the high-dimensional parameter space and evaluate
the uncertainties on model parameters. We show that the two spectra can be
successfully fit in near-equipartition conditions, defined by the ratio of the
energy density of relativistic electrons to magnetic field $\zeta_{\rm e}$
being close to unity. In both spectra, the observed X-rays are dominated by
synchrotron-self Compton photons, and the observed $\gamma$ rays are dominated
by Compton scattering of external infrared photons from a surrounding dusty
torus.
Model parameters are well constrained. From the low state to the high state,
both the curvature of the log-parabola width parameter and the synchrotron peak
frequency significantly increase. The derived magnetic fields in the two states
are nearly identical ($\sim1$\ G), but the Doppler factor in the high state is
larger than that in the low state ($\sim$28 versus $\sim$18). We derive that
the gamma-ray emission site takes place outside the broad-line region, at
$\gtrsim$ 0.1 pc from the black hole, but within the dusty torus. Implications
for 3C 279 as a source of high-energy cosmic-rays are discussed.
Understanding the trajectory of a coronal mass ejection (CME), including any deflection from a radial path, and the orientation of its magnetic field is essential for space weather predictions. Kay et al. (2015b) developed a model, Forecasting a CME's Altered Trajectory (ForeCAT), of CME deflections and rotation due to magnetic forces, not including the effects of reconnection. ForeCAT is able to reproduce the deflection of observed CMEs (Kay et al. 2015a). The deflecting CMEs tend to show a rapid increase of their angular momentum close to the Sun, followed by little to no increase at farther distances. Here we quantify the distance at which the CME deflection is "determined," which we define as the distance after which the background solar wind has negligible influence on the total deflection. We consider a wide range in CME masses and radial speeds and determine that the deflection and rotation of these CMEs can be well-described by assuming they propagate with constant angular momentum beyond 10 Rs. The assumption of constant angular momentum beyond 10 Rs yields underestimates of the total deflection at 1 AU of only 1% to 5% and underestimates of the rotation of 10%. Since the deflection from magnetic forces is determined by 10 Rs, non-magnetic forces must be responsible for any observed interplanetary deflections or rotations where the CME has increasing angular momentum.
High-energy cosmic rays impinging onto the atmosphere of the Earth initiate extensive air showers. With the LOFAR radio telescope and the Auger Engineering Radio Array (AERA) at the Pierre Auger Observatory radio emission from air showers is detected. Recent results are presented from both experiments. The measured properties of the radio emission are described. The measurements are used to derive the properties of high-energy cosmic rays: their arrival direction, energy, and particle type (mass).
The results of our speckle interferometric observations of young binaries VY Tau and XZ Tau are presented. For the first time, we found a relative displacement of VY Tau components as well as a preliminary orbit for XZ Tau. It appeared that the orbit is appreciably non-circular and is inclined by $i \lesssim 47^o$ from the plane of the sky. It means that the rotation axis of XZ Tau A and the axis of its jet are significantly non-perpendicular to the orbital plane. We found that the average brightness of XZ Tau had been increasing from the beginning of the last century up to the mid-thirties and then it decreased by $\Delta B > 2$ mag. The maximal brightness has been reached significantly later on the time of periastron passage. The total brightness of XZ Tau's components varied in a non-regular way from 1970 to 1985 when eruptions of hot gas from XZ Tau A presumably had occurred. In the early nineties the variations became regular following which a chaotic variability had renewed. We also report that a flare activity of VY Tau has resumed after 40 yr pause, parameters of the previous and new flares are similar, and the flares are related with the A component.
Experimental efforts to measure neutrinos by radio-frequency (RF) signals resulting from neutrino interactions in-ice have intensified over the last decade. Recent calculations indicate that one may dramatically improve the sensitivity of ultra-high energy ("UHE"; >EeV) neutrino experiments via detection of radio waves trapped along the air-ice surface. Detectors designed to observe the "Askaryan effect" currently search for RF electromagnetic pulses propagating through bulk ice, and could therefore gain sensitivity if signals are confined to the ice-air boundary. To test the feasibilty of this scenario, measurements of the complex radio-frequency properties of several air-dielectric interfaces were performed for a variety of materials. Two-dimensional surfaces of granulated fused silica (sand), both in the lab as well as occurring naturally, water doped with varying concentrations of salt, natural rock salt formations, granulated salt and ice itself were studied, both in North America and also Antarctica. In no experiment do we observe unambiguous surface wave propagation, as would be evidenced by signals travelling with reduced signal loss and/or superluminal velocities, compared to conventional EM wave propagation.
AA Tau is a well-studied, nearby classical T Tauri star, which is viewed almost edge-on. A warp in its inner disk periodically eclipses the central star, causing a clear modulation of its optical light curve. The system underwent a major dimming event beginning in 2011 caused by an extra absorber, which is most likely associated with additional disk material in the line of sight toward the central source. We present new XMM-Newton X-ray, Hubble Space Telescope FUV, and ground based optical and near-infrared data of the system obtained in 2013 during the long-lasting dim phase. The line width decrease of the fluorescent H$_2$ disk emission shows that the extra absorber is located at $r>1\,$au. Comparison of X-ray absorption ($N_H$) with dust extinction ($A_V$), as derived from measurements obtained one inner disk orbit (eight days) after the X-ray measurement, indicates that the gas-to-dust ratio as probed by the $N_H$ to $A_V$ ratio of the extra absorber is compatible with the ISM ratio. Combining both results suggests that the extra absorber, i.e., material at $r>1\,$au, has no significant gas excess in contrast to the elevated gas-to-dust ratio previously derived for material in the inner region ($\lesssim0.1\,$au).
The O IV intercombination lines observed by the Interface Region Imaging Spectrograph (IRIS) between 1397 and 1407 A provide useful density diagnostics. This document presents data that address two issues related to these lines: (1) the contribution of S IV to the O IV 1404.8 line; and (2) the range of sensitivity of the O IV 1399.8/1401.2 ratio.
Binary fraction and orbital characteristics provide indications on the
conditions of star formation, as they shed light on the environment they were
born in. Multiple systems are more common in low density environments rather
than in higher density ones. In the current debate about the formation of
Globular Clusters and their multiple populations, studying the binary incidence
in the populations they host offers a crucial piece of information on the
environment of their birth and their subsequent dynamical evolution.
Through a multi-year observational campaign using FLAMES at VLT, we monitored
the radial velocity of 968 Red-Giant Branch stars located around the half-light
radii in a sample of 10 Galactic Globular Clusters. We found a total of 21
radial velocity variables identified as {\it bona fide} binary stars, for a
binary fraction of 2.2%$\pm$0.5%. When separating the sample into first
generation and second generation stars, we find a binary fraction of
4.9%$\pm$1.3% and 1.2%$\pm$0.4% respectively. Through simulations that take
into account possible sources of bias in detecting radial velocity variations
in the two populations, we show that the difference is significant and only
marginally affected by such effects.
Such a different binary fraction strongly suggests different conditions in
the environment of formation and evolution of first and second generations
stars, with the latter being born in a much denser environment. Our result
hence strongly supports the idea that the second generation forms in a dense
subsystem at the center of the loosely distributed first generation, where
(loose) binaries are efficiently destroyed.
New-generation X-ray polarimeters currently under development promise to open a new window in the study of high-energy astrophysical sources. Among them, neutron stars appear particularly suited for polarization measurements. Radiation from the (cooling) surface of a neutron star is expected to exhibit a large intrinsic polarization degree due to the star strong magnetic field ($\approx 10^{12}-10^{15}$ G), which influences the plasma opacity in the outermost stellar layers. The polarization fraction and polarization angle as measured by an instrument, however, do not necessary coincide with the intrinsic ones derived from models of surface emission. This is due to the effects of quantum electrodynamics in the highly magnetized vacuum around the star (the vacuum polarization) coupled with the rotation of the Stokes parameters in the plane perpendicular to the line of sight induced by the non-uniform magnetic field. Here we revisit the problem and present an efficient method for computing the observed polarization fraction and polarization angle in the case of radiation coming from the entire surface of a neutron star, accounting for both vacuum polarization and geometrical effects due to the extended emitting region. Our approach is fairly general and is illustrated in the case of blackbody emission from a neutron star with either a dipolar or a (globally) twisted magnetic field.
We revisit the simplest model for dark matter. In this context the dark matter candidate is a real scalar field which interacts with the Standard Model particles through the Higgs portal. We discuss the relic density constraints as well as the predictions for direct and indirect detection. The final state radiation processes are investigated in order to understand the visibility of the gamma lines from dark matter annihilation. We find two regions where one could observe the gamma lines at gamma-ray telescopes. We point out that the region where the dark matter mass is between 100 and 300 GeV can be tested in the near future at direct and indirect detection experiments.
Newtonian cosmological perturbation equations valid to full nonlinear order are well known in the literature. Assuming the absence of the transverse-tracefree part of the metric, we present the general relativistic counterpart valid to full nonlinear order. The relativistic equations are presented without taking the slicing (temporal gauge) condition. The equations do have the proper Newtonian and first post-Newtonian limits. We also present the relativistic pressure correction terms in the Newtonian hydrodynamic equations.
One of the problems often encountered in X-ray mirror manufacturing is setting proper manufacturing tolerances to guarantee an angular resolution - often expressed in terms of Point Spread Function (PSF) - as needed by the specific science goal. To do this, we need an accurate metrological apparatus, covering a very broad range of spatial frequencies, and an affordable method to compute the PSF from the metrology dataset. [...] However, the separation between these spectral ranges is difficult do define exactly, and it is also unclear how to affordably combine the PSFs, computed with different methods in different spectral ranges, into a PSF expectation at a given X-ray energy. For this reason, we have proposed a method entirely based on the Huygens-Fresnel principle to compute the diffracted field of real Wolter-I optics, including measured defects over a wide range of spatial frequencies. Owing to the shallow angles at play, the computation can be simplified limiting the computation to the longitudinal profiles, neglecting completely the effect of roundness errors. Other authors had already proposed similar approaches in the past, but only in far-field approximation, therefore they could not be applied to the case of Wolter-I optics, in which two reflections occur in sequence within a short range. The method we suggest is versatile, as it can be applied to multiple reflection systems, at any X-ray energy, and regardless of the nominal shape of the mirrors in the optical system. The method has been implemented in the WISE code, successfully used to explain the measured PSFs of multilayer-coated optics for astronomic use, and of a K-B optical system in use at the FERMI free electron laser.
In this paper we present a model for dark energy in minimal supergravity with flat K\"ahler metric and a power-law superpotential. These choices of $K$ and $W$ can lead to a spontaneous supersymmetry breaking, with the minimum of the potential at $V(\varphi_0)=0$. We assume that the massive gravitino can decay into a scalar field with same potential than before, but expanded around $\Phi\equiv\varphi-\varphi_0=0$. This expanded potential $V(\Phi)$ leads to an accelerated expansion of the universe with density parameter equals $\Omega_\Phi=0.7$ today.
Links to: arXiv, form interface, find, astro-ph, recent, 1509, contact, help (Access key information)
As the evolutionary link between the radio millisecond pulsars (MSPs) and the low mass X-ray binaries or intermediate mass X-ray binaries, the millisecond X-ray pulsars (MSXPs) are important objects in testing theories of pulsar formation and evolution. In general, neutron stars in MSXPs can form via core collapse supernova (CC channel) of massive stars or accretion induced collapse (AIC channel) of an accreting ONeMg WD whose mass reaches the Chandrasekhar limit. Here, in addition to CC and AIC channels we also consider another channel, i.e., evolution induced collapse (EIC channel) of a helium star with mass between $1.4$ and $2.5 M_{\odot}$. Using a population synthesis code, we have studied MSXPs arising from three different evolutionary channels. We find that the Galactic birthrates of transient MSXPs and persistent MSXPs are about 0.7---$1.4\times 10^{-4}$ yr$^{-1}$. Our population synthesis calculations have shown that about 50\%---90\% of the MSXPs have undergone CC channel, about 10\%---40\% of them have undergone EIC channel, and the MSXPs via AIC channel are the least.
We present the results of a search for astrophysical sources of brief transient neutrino emission using IceCube and DeepCore data acquired between May 15th 2012 and April 30th 2013. While the search methods employed in this analysis are similar to those used in previous IceCube point source searches, the data set being examined consists of a sample of predominantly sub-TeV muon neu- trinos from the Northern Sky (-5$^{\circ}$ < {\delta} < 90$^{\circ}$ ) obtained through a novel event selection method. This search represents a first attempt by IceCube to identify astrophysical neutrino sources in this relatively unexplored energy range. The reconstructed direction and time of arrival of neutrino events is used to search for any significant self-correlation in the dataset. The data revealed no significant source of transient neutrino emission. This result has been used to construct limits at timescales ranging from roughly 1$\,$s to 10 days for generic soft-spectra transients. We also present limits on a specific model of neutrino emission from soft jets in core-collapse supernovae.
We use measurements from the Planck satellite mission and galaxy redshift surveys over the last decade to test three of the basic assumptions of the standard model of cosmology, $\Lambda$CDM: the spatial curvature of the universe, the nature of dark energy and the laws of gravity on large scales. We obtain improved constraints on several scenarios that violate one or more of these assumptions. We measure $w_0=-0.94\pm0.17$ (18\% measurement) and $1+w_a=1.16\pm0.36$ (31\% measurement) for models with a time-dependent equation of state, which is an improvement over current best constraints \citep{Aubourg2014}. In the context of modified gravity, we consider popular scalar tensor models as well as a parametrization of the growth factor. In the case of one-parameter $f(R)$ gravity models with a $\Lambda$CDM background, we constrain $B_0 < 1.36 \times 10^{-5} $ (1$\sigma$ C.L.), which is an improvement by a factor of 4 on the current best \citep{XU2015}. We provide the very first constraint on the coupling parameters of general scalar-tensor theory and stringent constraint on the only free coupling parameter of Chameleon models. We also derive constraints on extended Chameleon models, improving the constraint on the coupling by a factor of 6 on the current best \citep{Hojjati2011} . We also measure $\gamma = 0.612 \pm 0.072$ (11.7\% measurement) for growth index parametrization which is an improvement over the current best measurement of $\gamma = 0.699\pm0.110$ (16\%) \citep{Samushia14}. We improve all the current constraints by combining results from various galaxy redshift surveys in a coherent way, which includes a careful treatment of scale-dependence introduced by modified gravity.
Using a sample of dwarf galaxies observed using the VIMOS IFU on the VLT, we investigate the mass-metallicity relation (MZR) as a function of star formation rate (FMR$_{\text{SFR}}$) as well as HI-gas mass (FMR$_{\text{HI}}$). We combine our IFU data with a subsample of galaxies from the ALFALFA HI survey crossmatched to the Sloan Digital Sky Survey to study the FMR$_{\text{SFR}}$ and FMR$_{\text{HI}}$ across the stellar mass range 10$^{6.6}$ to 10$^{8.8}$ M$_\odot$, with metallicities as low as 12+log(O/H) = 7.67. We find the 1$\sigma$ mean scatter in the MZR to be 0.05 dex. The 1$\sigma$ mean scatter in the FMR$_{\text{SFR}}$ (0.02 dex) is significantly lower than that of the MZR. The FMR$_{\text{SFR}}$ is not consistent between the IFU observed galaxies and the ALFALFA/SDSS galaxies for SFRs lower than 10$^{-2.4}$ M$_\odot$ yr$^{-1}$, however this could be the result of limitations of our measurements in that regime. The lowest mean scatter (0.01 dex) is found in the FMR$_{\text{HI}}$. We also find that the FMR$_{\text{HI}}$ is consistent between the IFU observed dwarf galaxies and the ALFALFA/SDSS crossmatched sample. We introduce the fundamental metallicity luminosity counterpart to the FMR, again characterized in terms of SFR (FML$_{\text{SFR}}$) and HI-gas mass (FML$_{\text{HI}}$). We find that the FML$_{\text{HI}}$ relation is consistent between the IFU observed dwarf galaxy sample and the larger ALFALFA/SDSS sample. However the 1$\sigma$ scatter for the FML$_{\text{HI}}$ relation is not improved over the FMR$_{\text{HI}}$ scenario. This leads us to conclude that the FMR$_{\text{HI}}$ is the best candidate for a physically motivated fundamental metallicity relation.
The influence of both major and minor mergers is expected to significantly affect gradients of stellar ages and metallicities in the outskirts of galaxies. Measurements of observed gradients are beginning to reach large radii in galaxies, but a theoretical framework for connecting the findings to a picture of galactic build-up is still in its infancy. We analyze stellar populations of a statistically representative sample of quiescent galaxies over a wide mass range from the Illustris simulation. We measure metallicity and age profiles in the stellar halos of quiescent Illustris galaxies ranging in stellar mass from $10^{10}$ to $10^{12} M_\odot$, accounting for observational projection and luminosity-weighting effects. We find wide variance in stellar population gradients between galaxies of similar mass, with typical gradients agreeing with observed galaxies. We show that, at fixed mass, the fraction of stars born in-situ within galaxies is correlated with the metallicity gradient in the halo, confirming that stellar halos contain unique information about the build-up and merger histories of galaxies.
The interaction of ionizing and far-ultraviolet radiation with the interstellar medium is of great importance. It results in the formation of regions in which the gas is ionized, beyond which are photodissociation regions (PDRs) in which the gas transitions to its atomic and molecular form. Several numerical codes have been implemented to study these two main phases of the interstellar medium either dynamically or chemically. In this paper we present TORUS-3DPDR, a new self-consistent code for treating the chemistry of three-dimensional photoionization and photodissociation regions. It is an integrated code coupling the two codes TORUS, a hydrodynamics and Monte Carlo radiation transport code, and 3D-PDR, a photodissociation regions code. The new code uses a Monte Carlo radiative transfer scheme to account for the propagation of the ionizing radiation including the diffusive component as well as a ray-tracing scheme based on the HEALPix package in order to account for the escape probability and column density calculations. Here, we present the numerical techniques we followed and we show the capabilities of the new code in modelling three-dimensional objects including single or multiple sources. We discuss the effects introduced by the diffusive component of the UV field in determining the thermal balance of PDRs as well as the effects introduced by a multiple sources treatment of the radiation field. We find that diffuse radiation can positively contribute to the formation of CO. With this new code, three-dimensional synthetic observations for the major cooling lines are possible, for making feasible a detailed comparison between hydrodynamical simulations and observations.
The connection between dark matter halos and galactic baryons is often not well-constrained nor well-resolved in cosmological hydrodynamical simulations. Thus, Halo Occupation Distribution (HOD) models that assign galaxies to halos based on halo mass are frequently used to interpret clustering observations, even though it is well-known that the assembly history of dark matter halos is related to their clustering. In this paper we use high-resolution hydrodynamical cosmological simulations to compare the halo and stellar mass growth of galaxies in a large-scale overdensity to those in a large-scale underdensity (on scales of about 20 Mpc). The simulation reproduces assembly bias, that halos have earlier formation times in overdense environments than in underdense regions. We find that the stellar mass to halo mass ratio is larger in overdense regions in central galaxies residing in halos with masses between 10$^{11}$-10$^{12.9}$ M$_{\odot}$. When we force the local density (within 2 Mpc) at z=0 to be the same for galaxies in the large-scale over- and underdensities, we find the same results. We posit that this difference can be explained by a combination of earlier formation times, more interactions at early times with neighbors, and more filaments feeding galaxies in overdense regions. This result puts the standard practice of assigning stellar mass to halos based only on their mass, rather than considering their larger environment, into question.
This work describes a full Bayesian analysis of the Nearby Universe as traced by galaxies of the 2M++ survey. The analysis is run in two sequential steps. The first step self-consistently derives the luminosity dependent galaxy biases, the power-spectrum of matter fluctuations and matter density fields within a Gaussian statistic approximation. The second step makes a detailed analysis of the three dimensional Large Scale Structures, assuming a fixed bias model and a fixed cosmology. This second step allows for the reconstruction of both the final density field and the initial conditions at z=1000 assuming a fixed bias model. From these, we derive fields that self-consistently extrapolate the observed large scale structures. We give two examples of these extrapolation and their utility for the detection of structures: the visibility of the Sloan Great Wall, and the detection and characterization of the Local Void using DIVA, a Lagrangian based technique to classify structures.
I present a novel method for measuring lens masses for microlensing events. By combining a measured lens flux with the microlens parallax vector pi_E, it is possible to derive the mass of the lens system without knowing the angular size of the Einstein ring, theta_E. This enables mass and distance measurements for single, luminous lenses, as well as binary and planetary lenses without caustic crossings. I discuss applications of this method in the contexts of the Spitzer, Kepler, and WFIRST microlensing missions.
The Sloan Digital Sky Survey IV extended Baryonic Oscillation Spectroscopic Survey (SDSS-IV/eBOSS) will observe approximately 270,000 emission-line galaxies (ELGs) to measure the Baryonic Acoustic Oscillation standard ruler (BAO) at redshift 0.9. To test different ELG selection algorithms, based on data from several imaging surveys, 9,000 spectra were observed with the SDSS spectrograph as a pilot survey. First, we provide a detailed description of each target selection algorithm tested. Then, using visual inspection and redshift quality flags, we find that the automated spectroscopic redshifts assigned by the pipeline meet the quality requirements for a robust BAO measurement. Also, we show the correlations between sky emission, signal-to-noise ratio in the emission lines and redshift error. As a result, we provide robust redshift distributions for the different target selection schemes tested. Finally, we infer two optimal target selection algorithms to be applied on DECam photometry that fulfill the eBOSS survey efficiency requirements.
The origin of the stellar initial mass function (IMF) is a fundamental issue in the theory of star formation. It is generally fit with a composite power law. Some clues on the progenitors can be found in dense starless cores that have a core mass function (CMF) with a similar shape. In the low-mass end, these mass functions increase with mass, albeit the sample may be somewhat incomplete; in the high-mass end, the mass functions decrease with mass. There is an offset in the turn-over mass between the two mass distributions. The stellar mass for the IMF peak is lower than the corresponding core mass for the CMF peak in the Pipe Nebula by about a factor of three. Smaller offsets are found between the IMF and the CMFs in other nebulae. We suggest that the offset is likely induced during a starburst episode of global star formation which is triggered by the formation of a few O/B stars in the multi-phase media, which naturally emerged through the onset of thermal instability in the cloud-core formation process. We consider the scenario that the ignition of a few massive stars photoionizes the warm medium between the cores, increases the external pressure, reduces their Bonnor?Ebert mass, and triggers the collapse of some previously stable cores. We quantitatively reproduce the IMF in the low-mass end with the assumption of additional rotational fragmentation.
We present a detailed study of the internal kinematics of the Galactic Globular Cluster M 4 (NGC 6121), by deriving the radial velocities from 7250 spectra for 2771 stars distributed from the upper part of the Red Giant Branch down to the Main Sequence. We describe new approaches to determine the wavelength solution from day-time calibrations and to determine the radial velocity drifts that can occur between calibration and science observations when observing with the GIRAFFE spectrograph at VLT. Two techniques to determine the radial velocity are compared, after a qualitative description of their advantages with respect to other commonly used algorithm, and a new approach to remove the sky contribution from the spectra obtained with fibre-fed spectrograph and further improve the radial velocity precision is presented. The average radial velocity of the cluster is $\langle v \rangle = 71.08 \pm 0.08$ km s$^{-1}$ with an average dispersion of $\mu_{v_c} = 3.97$ km s$^{-1}$. Using the same dataset and the same statistical approach of previous analyses, 20 additional binary candidates are found, for a total of 87 candidates. A new determination of the internal radial velocity dispersion as a function of cluster distance is presented, resulting in a dispersion of $4.5$ km s$^{-1}$ within 2$^{\prime}$ from the center of cluster and steadily decreasing outward. We statistically confirm the small amplitude of the cluster rotation, as suggested in the past by several authors. This new analysis represents a significant improvement with respect to previous results in literature and provides a fundamental observational input for the modeling of the cluster dynamics.
We use a sample of 149 spectroscopically confirmed UV-selected galaxies at $z\sim 2$ to investigate the relative dust attenuation of the stellar continuum and the nebular emission lines. For each galaxy in the sample, at least one rest-frame optical emission line (H$\alpha$/[NII]$\lambda6583$ or [OIII]$\lambda5007$) measurement has been taken from the litterature, and 41 galaxies have additional Spitzer/MIPS 24$\mu$m observations that are used to infer infrared luminosities. We use a spectral energy distribution (SED) fitting code that predicts nebular line strengths when fitting the stellar populations of galaxies in our sample, and we perform comparisons between the predictions of our models and the observed/derived physical quantities. We find that on average our code is able to reproduce all the physical quantities (e.g., UV $\beta$ slopes, infrared luminosities, emission line fluxes), but we need to apply a higher dust correction to the nebular emission compared to the stellar emission for the largest SFR ($\log\mathrm{(SFR/M}_\odot\mathrm{yr}^{-1})>1.87$, Salpeter IMF). We find a correlation between SFR and the difference in nebular and stellar color excesses, which could resolve the discrepant results regarding nebular dust correction at $z\sim2$ from previous results.
We present results from a subset of simulations from the "Evolution and Assembly of GaLaxies and their Environments" (EAGLE) suite in which the formulation of the hydrodynamics scheme is varied. We compare simulations that use the same subgrid models without re-calibration of the parameters but employing the standard GADGET flavour of smoothed particle hydrodynamics (SPH) instead of the more recent state-of-the-art ANARCHY formulation of SPH that was used in the fiducial EAGLE runs. We find that the properties of most galaxies, including their masses and sizes, are not significantly affected by the details of the hydrodynamics solver. However, the star formation rates of the most massive objects are affected by the lack of phase mixing due to spurious surface tension in the simulation using standard SPH. This affects the efficiency with which AGN activity can quench star formation in these galaxies and it also leads to differences in the intragroup medium that affect the X-ray emission from these objects. The differences that can be attributed to the hydrodynamics solver are, however, likely to be less important at lower resolution. We also find that the use of a time step limiter is important for achieving the feedback efficiency required to match observations of the low-mass end of the galaxy stellar mass function.
One of the most powerful techniques to study the dark sector of the Universe is weak gravitational lensing. In practice, to infer the reduced shear, weak lensing measures galaxy shapes, which are the consequence of both the intrinsic ellipticity of the sources and of the integrated gravitational lensing effect along the line of sight. Hence, a very large number of galaxies is required in order to average over their individual properties and to isolate the weak lensing cosmic shear signal. If this `shape noise' can be reduced, significant advances in the power of a weak lensing surveys can be expected. This paper describes a general method for extracting the probability distributions of parameters from catalogues of data using Voronoi cells, which has several applications, and has synergies with Bayesian hierarchical modelling approaches. This allows us to construct a probability distribution for the variance of the intrinsic ellipticity as a function of galaxy property using only photometric data, allowing a reduction of shape noise. As a proof of concept the method is applied to the CFHTLenS survey data. We use this approach to investigate trends of galaxy properties in the data and apply this to the case of weak lensing power spectra.
We conduct a spectroscopic search of quasars observed by the Sloan Digital Sky Survey (SDSS) with broad absorption line (BAL) troughs due to Mg II and troughs due to Fe II that simultaneously exhibit strong Balmer narrow emission lines (NELs). We find that in a redshift range of 0.4 less than or equal to z less than or equal to 0.9 approximately 23 of the 70 Mg II BALs and 4 of a subset of 15 Fe II BALs exhibit strong Balmer emission. We also find significant fractions of Mg II BALs (approximately 23%) and those Mg II BALs with Fe II troughs (approximately 27%) have strong continuum reddening, E(B - V) greater than or equal to 0.1. From measurements of the Balmer decrement in three objects, we find similarly significant reddening of the NEL region in three of the four objects; the NELs in the fourth object are not measurable. We also include one object in this study not taken from the SDSS sample that shows Fe II absorption and strong narrow emission, but due to measurement uncertainty and low continuum reddening the comparison is consistent but inconclusive. We find a trend in both the Mg II and Fe II BAL samples between the NEL reddening and continuum reddening. Because the narrow line reddening is consistent with the continuum reddening in every object in the two SDSS samples, it suggests that the reddening sources in these objects likely exist at larger radial distances than the narrow line regions from the central nucleus.
Voyager 2 observations revealed that the internal luminosity of Neptune is an order of magnitude higher than that of Uranus. If the two planets have similar interior structures and cooling histories, the luminosity of Neptune can only be explained by invoking some energy source beyond gravitational contraction. This paper investigates whether Centaur impacts could provide the energy necessary to produce the luminosity of Neptune. The major findings are (1) that impacts on both Uranus and Neptune are too infrequent to provide luminosities of order the observed value for Neptune, even for optimistic impact-rate estimates, and (2) that Uranus and Neptune rarely have significantly different impact-generated luminosities at any given time. Uranus and Neptune most likely have structural differences that force them to cool and contract at different rates.
We investigate a novel Bayesian analysis method, based on the Stochastically Lighting Up Galaxies (slug) code, to derive the masses, ages, and extinctions of star clusters from integrated light photometry. Unlike many analysis methods, slug correctly accounts for incomplete IMF sampling, and returns full posterior probability distributions rather than simply probability maxima. We apply our technique to 621 visually-confirmed clusters in two nearby galaxies, NGC 628 and NGC 7793, that are part of the Legacy Extragalactic UV Survey (LEGUS). LEGUS provides Hubble Space Telescope photometry in the NUV, U, B, V, and I bands. We analyze the sensitivity of the derived cluster properties to choices of prior probability distribution, evolutionary tracks, IMF, metallicity, treatment of nebular emission, and extinction curve. We find that slug's results for individual clusters are insensitive to most of these choices, but that the posterior probability distributions we derive are often quite broad, and sometimes multi-peaked and quite sensitive to the choice of priors. In contrast, the properties of the cluster population as a whole are relatively robust against all of these choices. We also compare our results from slug to those derived with a conventional non-stochastic fitting code, Yggdrasil. We show that slug's stochastic models are generally a better fit to the observations than the deterministic ones used by Yggdrasil. However, the overall properties of the cluster populations recovered by both codes are qualitatively similar.
Hierarchical triple-star systems are expected to form frequently via close binary-binary encounters in the dense cores of globular clusters. In a sufficiently inclined triple, gravitational interactions between the inner and outer binary can cause large-amplitude oscillations in the eccentricity of the inner orbit ("Lidov-Kozai cycles"), which can lead to a collision and merger of the two inner components. In this paper we use Monte Carlo models of dense star clusters to identify all triple systems formed dynamically and we compute their evolution using a highly accurate three-body integrator which incorporates relativistic and tidal effects. We find that a large fraction of these triples evolve through a non-secular dynamical phase which can drive the inner binary to higher eccentricities than predicted by the standard secular perturbation theory (even including octupole-order terms). We place constraints on the importance of Lidov-Kozai-induced mergers for producing: (i) gravitational wave sources detectable by Advanced LIGO (aLIGO), for triples with an inner pair of stellar black holes; and (ii) blue straggler stars, for triples with main-sequence-star components. We find a realistic aLIGO detection rate of black hole mergers due to the Lidov-Kozai mechanism of 2yr^-1, with about 20% of these having a finite eccentricity when they first chirp into the aLIGO frequency band. While rare, these events are likely to dominate among eccentric compact object inspirals that are potentially detectable by aLIGO. For blue stragglers, we find that the Lidov-Kozai mechanism can contribute up to ~10% of their total numbers in globular clusters. In clusters with low central densities, ~10^{3}-10^{4} M_Sun pc^-3, up to ~40% of binary blue stragglers could have formed in dynamically assembled triples.
The WFIRST-AFTA 2.4 m telescope will provide in the next decade the opportunity to host a coronagraph for the imaging and spectroscopy of planets and disks. The telescope, however, is not ideal, given its obscured aperture. Only recently have coronagraph designs been thoroughly investigated that can efficiently work with this configuration. Three coronagraph designs, the hybrid Lyot, the shaped pupil, and the phase-induced amplitude-apodization complex mask coronagraph (PIAA-CMC) have been selected for further development by the AFTA project. Real-world testbed demonstrations of these have just begun, so for now the most reliable means of evaluating their potential performance comes from numerical modeling incorporating diffraction propagation, realistic system models, and simulated wavefront sensing and control. Here we present the methods of performance evaluation and results for the current coronagraph designs.
It is very controversial whether radio--X-ray correlation as defined in LH state of XRBs can extend to quiescent state (e.g., X-ray luminosity less than a critical value of $L_{\rm X,c} \sim10^{-5.5}L_{\rm Edd}$) or not. In this work, we collect a sample of XRBs and low luminosity active galactic nuclei (LLAGNs) with wide distribution of Eddington ratios to reexplore the fundamental plane between 5 GHz radio luminosity, $L_{\rm R}$, 2-10 keV X-ray luminosity, $L_{\rm X}$, and black hole (BH) mass, $M_{\rm BH}$, namely $\log L_{\rm R}=\xi_{\rm X} \log L_{\rm X}+\xi_{\rm M}\log M_{\rm BH}+\rm constant$. For the whole sample, we confirm the former fundamental plane of Merloni et al. and Falcke et al. that $\xi_{\rm X}\sim 0.6$ and $\xi_{\rm M}\sim 0.8$ even after including more quiescent BHs. The quiescent BHs follow the fundamental plane very well, and, however, FR I radio galaxies follow a steeper track comparing other BH sources. After excluding FR Is, we investigate the fundamental plane for BHs in quiescent state with $L_{\rm X}< L_{\rm X,c}$ and sub-Eddington BHs with $L_{\rm X}> L_{\rm X,c}$ respectively, and both subsamples have a similar slope, $\xi_{\rm X}\sim0.6$, which support that quiescent BHs may behave similar to those in low-hard state. We further select two subsamples of AGNs with BH mass in a narrow range (FR Is with $M_{\rm BH}=10^{8.8\pm0.4}$ and other LLAGNs with $M_{\rm BH}=10^{8.0\pm0.4}$) to simulate the behavior of a single supermassive BH evolving from sub-Eddington to quiescent state. We find that the highly sub-Eddington sources with $L_{\rm X}/L_{\rm Edd}\sim10^{-6}-10^{-9}$ still roughly stay on the extension of radio--X-ray correlation as defined by other sub-Eddington BHs. Our results are consistent with several recent observations in XRBs that the radio--X-ray correlation as defined in low-hard state can extend to highly sub-Eddington quiescent state.
Existing planet-finding spectrometers are limited by systematic errors that result from their seeing-limited design. Of particular concern is the use of multi-mode fibers (MMFs), which introduce modal noise and accept significant amounts of background radiation from the sky. We present the design of a single-mode fiber-based acquisition camera for a diffraction-limited spectrometer named "iLocater." By using the "extreme" adaptive optics (AO) system of the Large Binocular Telescope (LBT), iLocater will overcome the limitations that prevent Doppler instruments from reaching their full potential, allowing precise radial velocity (RV) measurements of terrestrial planets around nearby bright stars. The instrument presented in this paper, which we refer to as the acquisition camera "demonstration system," will measure on-sky single-mode fiber (SMF) coupling efficiency using one of the 8.4m primaries of the LBT in fall 2015.
Cosmic rays are a fundamental source of ionization for molecular and diffuse clouds, influencing their chemical, thermal, and dynamical evolution. The amount of cosmic rays inside a cloud also determines the $\gamma$-ray flux produced by hadronic collisions between cosmic rays and cloud material. We study the spectrum of cosmic rays inside and outside of a diffuse cloud, by solving the stationary transport equation for cosmic rays including diffusion, advection and energy losses due to ionization of neutral hydrogen atoms. We found that the cosmic ray spectrum inside a diffuse cloud differs from the one in the interstellar medium (ISM) for energies smaller than $E_{br}\approx 100$ MeV, irrespective of the model details. Below $E_{br}$, the spectrum is harder (softer) than that in the ISM if the latter is a power law $\propto p^{-s}$ with $s$ larger (smaller) than $\sim0.42$. As a consequence also the ionization rate due to CRs is strongly affected. Assuming an average Galactic spectrum similar to the one inferred from AMS-2 and Voyager 1 data, we discuss the resulting ionization rate in a typical diffuse cloud.
We calculate the one-point probability density distribution functions (PDF) and the power spectra of the thermal and kinetic Sunyaev-Zeldovich (tSZ and kSZ) effects and the mean Compton Y parameter using the Magneticum Pathfinder simulations, state-of-the-art cosmological hydrodynamical simulations of a large cosmological volume of (896 Mpc/h)^3. These simulations follow in detail the thermal and chemical evolution of the intracluster medium as well as the evolution of super-massive black holes and their associated feedback processes. We construct full-sky maps of tSZ and kSZ from the light-cones out to z=0.17, and one realization of 8.8x8.8 degree wide, deep light-cone out to z=5.2. The local universe at z<0.027 is simulated by a constrained realisation. The tail of the one-point PDF of tSZ from the deep light-cone follows a power-law shape with an index of -3.2. Once convolved with the effective beam of Planck, it agrees with the PDF measured by Planck. The predicted tSZ power spectrum agrees with that of the Planck data at all multipoles up to l~1000, once the calculations are scaled to the Planck 2015 cosmological parameters with \Omega_m=0.308 and \sigma_8=0.8149. Consistent with the results in the literature, however, we continue to find the tSZ power spectrum at l=3000 that is significantly larger than that estimated from the high-resolution ground-based data. The simulation predicts the mean fluctuating Compton Y value of <Y>=1.18x10^{-6} for \Omega_m=0.272 and \sigma_8=0.809. Nearly half (~ 5x10^{-7}) of the signal comes from halos below a virial mass of 10^{13}M_\odot/h. Scaling this to the Planck 2015 parameters, we find <Y>=1.57x10^{-6}. The PDF and the power spectrum of kSZ from our simulation agree broadly with the previous work.
Gamma rays and microwave observations of the Galactic Center and surrounding areas indicate the presence of anomalous emission, whose origin remains ambiguous. The possibility of dark matter (DM) annihilation explaining both signals through prompt emission at gamma-rays and secondary emission at microwave frequencies from interactions of high-energy electrons produced in annihilation with the Galactic magnetic fields has attracted much interest in recent years. We investigate the DM interpretation of the Galactic Center gamma-ray excess by searching for the associated synchrotron in the WMAP-Planck data. Considering various magnetic field and cosmic-ray propagation models, we predict the synchrotron emission due to DM annihilation in our Galaxy, and compare it with the WMAP-Planck data at 23-70GHz. In addition to standard microwave foregrounds, we separately model the microwave counterpart to the Fermi Bubbles and the signal due to DM, and use component separation techniques to extract the signal associated with each template from the total emission. We confirm the presence of the Haze at the level of 7% of the total sky intensity at 23GHz in our chosen region of interest, with a harder spectrum $I \sim \nu^{-0.8}$ than the synchrotron from regular cosmic-ray electrons. The data do not show a strong preference towards fitting the Haze by either the Bubbles or DM emission only. Inclusion of both components provides a better fit with a DM contribution to the Haze emission of 20% at 23GHz, however, due to significant uncertainties in foreground modeling, we do not consider this a clear detection of a DM signal. We set robust upper limits on the annihilation cross section by ignoring foregrounds, and also report best-fit DM annihilation parameters obtained from a complete template analysis. We conclude that the WMAP-Planck data are consistent with a DM interpretation of the gamma-ray excess.
We comment on the paper "Dark Matter collisions with the Human Body" by K.~Freese and C.~Savage (Phys.\ Lett.\ B {\bf 717}, 25 (2012) [arXiv:1204.1339]) and describe a dark matter model for which the results of the previous paper do not apply. Within this mirror dark matter model, potentially hazardous objects, mirror micrometeorites, can exist potentially leading to diseases triggered by multiple mutations, such as cancer.
We report the detection of extended Ly alpha emission around individual star-forming galaxies at redshifts z = 3-6 in an ultradeep exposure of the Hubble Deep Field South obtained with MUSE on the ESO-VLT. The data reach a limiting surface brightness (1sigma) of ~1 x 10^-19 erg s^-1 cm^-2 arcsec^-2 in azimuthally averaged radial profiles, an order of magnitude improvement over previous narrowband imaging. Our sample consists of 26 spectroscopically confirmed Ly alpha-emitting, but mostly continuum-faint (m_AB >~ 27) galaxies. In most objects the Ly alpha emission is considerably more extended than the UV continuum light. While 5 of the faintest galaxies in the sample show no significantly detected Ly alpha haloes, the derived upper limits suggest that this is just due to insufficient S/N. Ly alpha haloes therefore appear to be (nearly) ubiquitous even for low-mass (~10^8-10^9 M_sun) star-forming galaxies at z>3. We decompose the Ly alpha emission of each object into a compact `continuum-like' and an extended halo component, and infer sizes and luminosities of the haloes. The extended Ly alpha emission approximately follows an exponential surface brightness distribution with a scale length of a few kpc. While these haloes are thus quite modest in terms of their absolute sizes, they are larger by a factor of 5-15 than the corresponding rest-frame UV continuum sources as seen by HST. They are also much more extended, by a factor ~5, than Ly alpha haloes around low-redshift star-forming galaxies. Between ~40% and >90% of the observed Ly alpha flux comes from the extended halo component, with no obvious correlation of this fraction with either the absolute or the relative size of the Ly alpha halo. Our observations provide direct insights into the spatial distribution of at least partly neutral gas residing in the circumgalactic medium of low to intermediate mass galaxies at z > 3.
Stars on the asymptotic giant branch (AGB) lose substantial amounts of matter, to the extent that they are important for the chemical evolution of, and dust production in, the universe. The mass loss is believed to increase gradually with age on the AGB, but it may also occur in the form of bursts, possibly related to the thermal pulsing phenomenon. Detached, geometrically thin, CO shells around carbon stars are good signposts of brief and intense mass ejection. We aim to put further constraints on the physical properties of detached CO shells around AGB stars. The photodissociation of CO and other carbon-bearing species in the shells leads to the possibility of detecting lines from neutral carbon. We have therefore searched for the CI($^3P_1-\,^3P_0$) line at 492 GHz towards two carbon stars, S Sct and R Scl, with detached CO shells of different ages, about 8000 and 2300 years, respectively. The CI($^3P_1-\,^3P_0$) line was detected towards R Scl. The line intensity is dominated by emission from the detached shell. The detection is at a level consistent with the neutral carbon coming from the full photodissociation of all species except CO, and with only limited photoionisation of carbon. The best fit to the observed $^{12}$CO and $^{13}$CO line intensities, assuming a homogeneous shell, is obtained for a shell mass of about 0.002 $M_\odot$, a temperature of about 100 K, and a CO abundance with respect to H$_2$ of 10$^{-3}$. The estimated CI/CO abundance ratio is about 0.3 for the best-fit model. However, a number of arguments point in the direction of a clumpy medium, and a viable interpretation of the data within such a context is provided.
We investigate the influence of interactions on the star formation by studying a sample of almost 1500 of the nearest galaxies, all within a distance of ~45 Mpc. We define the massive star formation rate (SFR), as measured from far-IR emission, and the specific star formation rate (SSFR), which is the former quantity normalised by the stellar mass of the galaxy, and explore their distribution with morphological type and with stellar mass. We then calculate the relative enhancement of these quantities for each galaxy by normalising them by the median SFR and SSFR values of individual control populations of similar non-interacting galaxies. We find that both SFR and SSFR are enhanced in interacting galaxies, and more so as the degree of interaction is higher. The increase is, however, moderate, reaching a maximum of a factor of 1.9 for the highest degree of interaction (mergers). The SFR and SSFR are enhanced statistically in the population, but in many individual interacting galaxies they are not enhanced at all. We discuss how those galaxies with the largest SFR and/or SSFR enhancement can be defined as starbursts. This study is based on a representative sample of nearby galaxies, including many low-mass and dwarf/irregular galaxies, and we argue that it should be used to place constraints on studies based on samples of galaxies at larger distances, beyond the local Universe.
We use the GALEX data of the GUViCS survey to construct the NUV luminosity function of the Virgo cluster over ~ 300 deg.2, an area covering the cluster and its surrounding regions up to ~ 1.8 virial radii. The NUV luminosity function is also determined for galaxies of different morphological type and NUV-i colour, and for the different substructures within the cluster. These luminosity functions are robust vs. statistical corrections since based on a sample of 833 galaxies mainly identified as cluster members with spectroscopic redshift (808) or high-quality optical scaling relations (10). We fit these luminosity functions with a Schechter function, and compare the fitted parameters with those determined for other nearby clusters and for the field. The faint end slope of the Virgo NUV luminosity function (alpha = -1.19), here sampled down to ~ NUV = -11.5 mag, is significantly flatter than the one measured in other nearby clusters and similar to the field one. Similarly M* = -17.56 is one-to-two magnitudes fainter than measured in Coma, A1367, the Shapley supercluster, and the field. These differences seem due to the quite uncertain statistical corrections and the small range in absolute magnitude sampled in these clusters. We do not observe strong systematic differences in the overall NUV luminosity function of the core of the cluster with respect to that of its periphery. We notice, however, that the relative contribution of red-to-blue galaxies at the faint end is inverted, with red quiescent objects dominating the core of the cluster and star forming galaxies dominating beyond one virial radius. These observational evidences are discussed in the framework of galaxy evolution in dense environments.
Gamma-ray bursts (GRBs) are some of the most energetic and exotic events in the Universe, however their behaviour at the highest energies (>10 GeV) is largely unknown. Although the Fermi-LAT space telescope has detected several GRBs in this energy range, it is limited by the relatively small collection area of the instrument. The H.E.S.S. experiment has now entered its second phase by adding a fifth telescope of 600 m$^{2}$ mirror area to the centre of the array. This new telescope increases the energy range of the array, allowing it to probe the sub-100 GeV range while maintaining the large collection area of ground based gamma-ray observatories, essential to probing short-term variability at these energies. We will present a description of the GRB observation scheme used by the H.E.S.S. experiment, summarising the behaviour and performance of the rapid GRB repointing system, the conditions under which potential GRB repointings are made and the data analysis scheme used for these observations.
We present unpublished data from a tidal disruption candidate in NGC 3599 which show that the galaxy was already X-ray bright 18 months before the measurement which led to its classification. This removes the possibility that the flare was caused by a classical, fast-rising, short-peaked, tidal disruption event. Recent relativistic simulations indicate that the majority of disruptions will actually take months or years to rise to a peak, which will then be maintained for longer than previously thought. NGC 3599 could be one of the first identified examples of such an event. The optical spectra of NGC 3599 indicate that it is a low-luminosity Seyfert/LINER with L_bol~10^40 ergs/s The flare may alternatively be explained by a thermal instability in the accretion disc, which propagates through the inner region at the sound speed, causing an increase of the disc scale height and local accretion rate. This can explain the <9 years rise time of the flare. If this mechanism is correct then the flare may repeat on a timescale of several decades as the inner disc is emptied and refilled.
Soft X-ray solar bursts are studied. Weak bursts with powers up to 10-8 W/m2 were detected. All the events were confirmed by GOES observations. Parameters of these microflares are determined. A physical mechanism for the low-intensity solar events is discussed.
Planck data towards the galaxy M82 are analyzed in the 70, 100 and 143 GHz bands. A substantial north-south and East-West temperature asymmetry is found, extending up to 1 degree from the galactic center. Being almost frequency-independent, these temperature asymmetries are indicative of a Doppler-induced effect regarding the line-of-sight dynamics on the halo scale, the ejections from the galactic center and, possibly, even the tidal interaction with M81 galaxy. The temperature asymmetry thus acts as a model-independent tool to reveal the bulk dynamics in nearby edge-on spiral galaxies, like the Sunyaev-Zeldovich effect for clusters of galaxies.
We conduct a multiplicity survey of members of the rho Ophiuchus cloud complex with high resolution imaging to characterize the multiple star population of this nearby star forming region and investigate the relation between stellar multiplicity and star and planet formation. Our aperture masking survey reveals the presence of 5 new stellar companions beyond the reach of previous studies, but does not result in the detection of any new substellar companions. We find that 43+/-6% of the 114 stars in our survey have stellar mass companions between 1.3-780AU, while 7 (+8 -5)% host brown dwarf companions in the same interval. By combining this information with knowledge of disk-hosting stars, we show that the presence of a close binary companion (separation < 40 AU) significantly influences the lifetime of protoplanetary disks, a phenomenon previously seen in older star forming regions. At the ~1-2Myr age of our Ophiuchus members ~2/3 of close binary systems have lost their disks, compared to only ~30% of single stars and wide binaries. This has significant impact on the formation of giant planets, which are expected to require much longer than 1 Myr to form via core accretion and thus planets formed via this pathway should be rare in close binary systems.
We study the HI K-band Tully-Fisher relation and the baryonic Tully-Fisher relation for a sample of 16 early-type galaxies, taken from the ATLAS3D sample, which all have very regular HI disks extending well beyond the optical body (> 5 R_eff). We use the kinematics of these disks to estimate the circular velocity at large radii for these galaxies. We find that the Tully-Fisher relation for our early-type galaxies is offset by about 0.5-0.7 magnitudes from the relation for spiral galaxies. The residuals with respect to the spiral Tully-Fisher relation correlate with estimates of the stellar mass-to-light ratio, suggesting that the offset between the relations is mainly driven by differences in stellar populations. We also observe a small offset between our Tully-Fisher relation with the relation derived for the ATLAS3D sample based on CO data representing the galaxies' inner regions (< 1 R_eff). This indicates that the circular velocities at large radii are systematically 10% lower than those near 0.5-1 R_eff, in line with recent determinations of the shape of the mass profile of early-type galaxies. The baryonic Tully-Fisher relation of our sample is distinctly tighter than the standard one, in particular when using mass-to-light ratios based on dynamical models of the stellar kinematics. We find that the early-type galaxies fall on the spiral baryonic Tully-Fisher relation if one assumes M/L_K = 0.54 M_sun/L_sun for the stellar populations of the spirals, a value similar to that found by recent studies of the dynamics of spiral galaxies. Such a mass-to-light ratio for spiral galaxies would imply that their disks are 60-70% of maximal. Our analysis increases the range of galaxy morphologies for which the baryonic Tully-Fisher relations holds, strengthening previous claims that it is a more fundamental scaling relation than the classical Tully-Fisher relation.
The lunar Askaryan technique, which involves searching for Askaryan radio
pulses from particle cascades in the outer layers of the Moon, is a method for
using the lunar surface as an extremely large detector of ultra-high-energy
particles. The high time resolution required to detect these pulses, which have
a duration of around a nanosecond, puts this technique in a regime quite
different from other forms of radio astronomy, with a unique set of associated
technical challenges which have been addressed in a series of experiments by
various groups. Implementing the methods and techniques developed by these
groups for detecting lunar Askaryan pulses will be important for a future
experiment with the Square Kilometre Array (SKA), which is expected to have
sufficient sensitivity to allow the first positive detection using this
technique.
Key issues include correction for ionospheric dispersion, beamforming,
efficient triggering, and the exclusion of spurious events from radio-frequency
interference. We review the progress in each of these areas, and consider the
further progress expected for future application with the SKA.
This paper is devoted to discuss the difference in the thermodynamic entropy budget {\it per baryon} in each type of stellar object found in Universe. We track and discuss the actual {\it decrease} of the stored baryonic thermodynamic entropy from the most primitive molecular cloud up to the final fate of matter in the black holes, passing through evolved states of matter as found in white dwarfs and neutron stars. We then discuss the case of actual stars of different masses throughout their {\it evolution}, clarifying the role of virial equilibrium condition for the decrease of the entropy and related issues. Finally, we discuss how gravity ultimately drives composition, hence structural changes along the stellar evolution all the way until the ultimate collapse to black holes, which may increase dramatically their entropy because of the gravitational contribution itself.
EUV (Extreme-Ultraviolet) waves are globally propagating disturbances that have been observed since the era of the SoHO/EIT instrument. Although the kinematics of the wave front and secondary wave components have been widely studied, there is not much known about the generation and plasma properties of the wave. In this paper we discuss the effect of an EUV wave on the local plasma as it passes through the corona. We studied the EUV wave, generated during the 2011 February 15 X-class flare/CME event, using Differential Emission Measure diagnostics. We analyzed regions on the path of the EUV wave and investigated the local density and temperature changes. From our study we have quantitatively confirmed previous results that during wave passage the plasma visible in the Atmospheric Imaging Assembly (AIA) 171A channel is getting heated to higher temperatures corresponding to AIA 193A and 211A channels. We have calculated an increase of 6 - 9% in density and 5 - 6% in temperature during the passage of the EUV wave. We have compared the variation in temperature with the adiabatic relationship and have quantitatively demonstrated the phenomenon of heating due to adiabatic compression at the wave front. However, the cooling phase does not follow adiabatic relaxation but shows slow decay indicating slow energy release being triggered by the wave passage. We have also identified that heating is taking place at the front of the wave pulse rather than at the rear. Our results provide support for the case that the event under study here is a compressive fast-mode wave or a shock.
Polarised light from astronomical targets can yield a wealth of information about their source radiation mechanisms, and about the geometry of the scattered light regions. Optical observations, of both the linear and circular polarisation components, have been impeded due to non-optimised instrumentation. The need for suitable observing conditions and the availability of luminous targets are also limiting factors. GASP uses division of amplitude polarimeter (DOAP) (Compain and Drevillon) to measure the four components of the Stokes vector simultaneously, which eliminates the constraints placed upon the need for moving parts during observation, and offers a real-time complete measurement of polarisation. Results from the GASP calibration are presented in this work for both a 1D detector system, and a pixel-by-pixel analysis on a 2D detector system. Following Compain et al. we use the Eigenvalue Calibration Method (ECM) to measure the polarimetric limitations of the instrument for each of the two systems. Consequently, the ECM is able to compensate for systematic errors introduced by the calibration optics, and it also accounts for all optical elements of the polarimeter in the output. Initial laboratory results of the ECM are presented, using APD detectors, where errors of 0.2% and 0.1{\deg} were measured for the degree of linear polarisation and polarisation angle respectively. Channel-to-channel image registration is an important aspect of 2-D polarimetry. We present our calibration results of the measured Mueller matrix of each sample, used by the ECM. A set of Zenith flat-field images were recorded during an observing campaign at the Palomar 200 inch telescope in November 2012. From these we show the polarimetric errors from the spatial polarimetry indicating both the stability and absolute accuracy of GASP.
Understanding the history and the evolution of the Milky Way disc is one of the main goals of modern astrophysics. We study the velocity dispersion behaviour of Galactic disc stars as a function of the [Mg/Fe] ratio, which can be used as a proxy of relative age. This key relation is essential to constrain the formation mechanisms of the disc stellar populations as well as the cooling processes. We used the recommended parameters and chemical abundances of 7800 FGK Milky Way field stars from the second internal data release of the Gaia-ESO Survey. These stars were observed with the GIRAFFE spectrograph, and cover a large spatial volume (6<R<10kpc and |Z|<2kpc). Based on the [Mg/Fe] and [Fe/H] ratios, we separated the thin- from the thick-disc sequence. From analysing the Galactocentric velocity of the stars for the thin disc, we find a weak positive correlation between Vphi and [Fe/H], due to a slowly rotating Fe-poor tail. For the thick disc, a strong correlation with [Fe/H] and [Mg/Fe] is established. We have detected an inversion of the radial velocity dispersion with [Mg/Fe] for thick-disc stars with [Fe/H]<-0.1dex and [Mg/Fe]>+0.2dex. First, the velocity dispersion increases with [Mg/Fe] at all [Fe/H] ratios for the thin-disc stars, and then it decreases for the thick-disc at the highest [Mg/Fe] abundances. Similar trends are observed within the errors for the azimuthal velocity dispersion, while a continuous increase with [Mg/Fe] is observed for the vertical velocity dispersion. The velocity dispersion decrease agrees with previous measurements of the RAVE survey, although it is observed here for a greater metallicity interval and a larger spatial volume. We confirm the existence of [Mg/Fe]-rich thick-disc stars with cool kinematics in the generally turbulent context of the primitive Galactic disc. This is discussed in the framework of the different disc formation scenarios.
The detection of trends or gradients in the transmission spectrum of extrasolar planets is possible with observations at very low spectral resolution. Transit measurements of sufficient accuracy using selected broad-band filters allow for an initial characterization of the atmosphere of the planet. We obtained time series photometry of 20 transit events and analyzed them homogeneously, along with eight light curves obtained from the literature. In total, the light curves span a range from 0.35 to 1.25 microns. During two observing seasons over four months each, we monitored the host star to constrain the potential influence of starspots on the derived transit parameters. We rule out the presence of a Rayleigh slope extending over the entire optical wavelength range, a flat spectrum is favored for HAT-P-12b with respect to a cloud-free atmosphere model spectrum. A potential cause of such gray absorption is the presence of a cloud layer at the probed latitudes. Furthermore, in this work we refine the transit parameters, the ephemeris and perform a TTV analysis in which we found no indication for an unseen companion. The host star showed a mild non-periodic variability of up to 1%. However, no stellar rotation period could be detected to high confidence.
Measurements of microwave attenuation at room temperature and 4.2 K have been performed on some conductors commonly used in receiver input circuits. The reduction in loss on cooling is substantial, particularly for copper and plated gold, both of which showed a factor of 3 loss reduction. Copper passivated with benzotriazole shows the same loss as without passivation. The residual resistivity ratio between room temperature and 4.2 K, deduced from the measurements using the classical skin effect formula, was smaller than the measured DC value to a degree consistent with conduction in the extreme anomalous skin effect regime at cryogenic temperatures. The measurements were made in the 5-10 GHz range. The materials tested were: aluminum alloys 1100-T6 and 6061-O, C101 copper, benzotriazole treated C101 copper, and brass plated with electroformed copper, Pur-A-Gold 125-Au soft gold, and BDT200 bright gold.
Since M4.5 - M9 dwarfs exhibit equatorial rotation velocities of order 10 km/s on average, radial velocity surveys targeting this stellar population will likely need to find methods to effectively remove starspot jitter. We present the first high resolution Doppler images of the M4.5 dwarf, GJ 791.2A, and the M9 dwarf, LP 944-20. The time series spectra of both objects reveal numerous line profile distortions over the rotation period of each star which we interpret as starspots. The transient distortions are modelled with spot/photosphere contrast ratios that correspond to model atmosphere temperature differences of Tphot-Tspot = 200 and 300 K. GJ 791.2A is a fully convective star with vsini = 35.1 km/s. Although we find more starspot structure at high latitudes, we reconstruct spots at a range of phases and latitudes with a mean spot filling of ~3%. LP 944-20 is one of the brightest known late-M dwarfs, with spectral type M9V and vsini = 30.8 km/s. Its spectral time series exhibits two dominant transient line distortions that are reconstructed as high latitude spots, while a mean spot filling factor of only 1.5% is found. The occurrence of low-contrast spots at predominantly high latitudes, which we see in both targets here, is in general likely to be responsible for the low amplitude photometric variability seen in late-M dwarfs. For GJ 791.2A, the radial velocities induced by the starspot features yield an rms velocity variability of 138 m/s, which can be reduced by a factor of 1.9 using our reconstructed surface brightness distributions.
We present the results of SPH simulations in which two clouds, each having mass $M_{_{\rm{o}}}\!=\!500\,{\rm M}_{_\odot}$ and radius $R_{_{\rm{o}}}\!=\!2\,{\rm pc}$, collide head-on at relative velocities of $\Delta v_{_{\rm{o}}} =2.4,\;2.8,\;3.2,\;3.6\;{\rm and}\;4.0\,{\rm km}\,{\rm s}^{-1}$. There is a clear trend with increasing $\Delta v_{_{\rm{o}}}$. At low $\Delta v_{_{\rm{o}}}$, star formation starts later, and the shock-compressed layer breaks up into an array of predominantly radial filaments; stars condense out of these filaments and fall, together with residual gas, towards the centre of the layer, to form a single large-$N$ cluster, which then evolves by competitive accretion, producing one or two very massive protostars and a diaspora of ejected (mainly low-mass) protostars; the pattern of filaments is reminiscent of the hub and spokes systems identified recently by observers. At high $\Delta v_{_{\rm{o}}}$, star formation occurs sooner and the shock-compressed layer breaks up into a network of filaments; the pattern of filaments here is more like a spider's web, with several small-$N$ clusters forming independently of one another, in cores at the intersections of filaments, and since each core only spawns a small number of protostars, there are fewer ejections of protostars. As the relative velocity is increased, the {\it mean} protostellar mass increases, but the {\it maximum} protostellar mass and the width of the mass function both decrease. We use a Minimal Spanning Tree to analyse the spatial distributions of protostars formed at different relative velocities.
We present new constraints on the high-temperature emission measure of a non-flaring solar active region using observations from the recently flown Focusing Optics X-ray Solar Imager sounding rocket payload. FOXSI has performed the first focused hard X-ray (HXR) observation of the Sun in its first successful flight on 2012 November 2. Focusing optics, combined with small strip detectors, enable high-sensitivity observations with respect to previous indirect imagers. This capability, along with the sensitivity of the HXR regime to high-temperature emission, offers the potential to better characterize high-temperature plasma in the corona as predicted by nanoflare heating models. We present a joint analysis of the differential emission measure (DEM) of active region 11602 using coordinated observations by FOXSI, Hinode/XRT and Hinode/EIS. The Hinode-derived DEM predicts significant emission measure between 1 MK and 3 MK, with a peak in the DEM predicted at 2.0-2.5 MK. The combined XRT and EIS DEM also shows emission from a smaller population of plasma above 8 MK. This is contradicted by FOXSI observations that significantly constrain emission above 8 MK. This suggests that the Hinode DEM analysis has larger uncertainties at higher temperatures and that >8 MK plasma above an emission measure of 3x10^44 cm^-3 is excluded in this active region.
We present results for Vela C obtained during the 2012 flight of the Balloon-borne Large Aperture Submillimeter Telescope for Polarimetry (BLASTPol). We mapped polarized intensity across almost the entire extent of this giant molecular cloud, in bands centered at 250, 350, and 500 {\mu}m. In this initial paper, we show our 500 {\mu}m data smoothed to a resolution of 2.5 arcminutes (approximately 0.5 pc). We show that the mean level of the fractional polarization p and most of its spatial variations can be accounted for using an empirical three-parameter power-law fit, p = p_0 N^(-0.4) S^(-0.6), where N is the hydrogen column density and S is the polarization-angle dispersion on 0.5 pc scales. The decrease of p with increasing S is expected because changes in the magnetic field direction within the cloud volume sampled by each measurement will lead to cancellation of polarization signals. The decrease of p with increasing N might be caused by the same effect, if magnetic field disorder increases for high column density sightlines. Alternatively, the intrinsic polarization efficiency of the dust grain population might be lower for material along higher density sightlines. We find no significant correlation between N and S. Comparison of observed submillimeter polarization maps with synthetic polarization maps derived from numerical simulations provides a promising method for testing star formation theories. Realistic simulations should allow for the possibility of variable intrinsic polarization efficiency. The measured levels of correlation among p, N, and S provide points of comparison between observations and simulations.
The Fermi-LAT data revealed giant bubbles of emission above and below the Galactic plane with an energy spectrum significantly harder than seen from other directions. How the bubbles connect to the Galactic plane is unclear. We find that bubble-like emission is not only found in the halo, but it is strongly present in the Galactic plane as well with a morphology close to the spatial distribution of the 1.8 MeV gamma-ray line from 26Al, a radioactive nucleus synthesized in SNRs. In addition, the spectral shape of this hard component coincides with the predicted spectrum from cosmic rays trapped in sources. Hence, we propose that the bubble-like emission in the plane has a hadronic origin, which arises from SCRs. The bubbles in the halo have the same energy spectrum, which suggests that they are outflows from the plane with the gamma-rays arising from hadronic interactions of protons trapped in a plasma of advected gas. Towards the Galactic center (GC) we observe the 1-3 GeV excess, but observe this excess in addition in all regions where there is strong 26Al production. This excludes the dark matter annihilation interpretation. Instead, we propose that the excess is caused by a deficit of the CRs at low rigidities in the regions of dense molecular clouds, which are characterized by higher energy losses and stronger stellar winds. The 26Al line is a tracer of such regions. If we introduce a break in the cosmic ray injection spectra in the regions of strong 26Al production to effectively take the deficit of low rigidity CRs into account, the excess in the Galactic plane and towards the GC both disappear. We thus conclude that the excess in the Galactic center is largely an artefact from the excess in the Galactic plane. The correlation between the two originates simply from the lines-of-sight towards the GC crossing the Galactic plane.
Preliminary results of the ongoing search for symbiotic binary stars in the Local Group of Galaxies are presented and discussed.
We present the first X-ray observations of the Taffy galaxies (UGC 12914/5) with the Chandra observatory, and detect soft X-ray emission in the region of the gas-rich, radio-continuum-emitting Taffy bridge. The results are compared to Herschel observations of dust and diffuse [CII] line-emitting gas. The diffuse component of the Taffy bridge has an X-ray luminosity of L(X) (0.5-8keV) =5.4 x 10^39 erg s^-1, which accounts for 19% of the luminosity of the sum for the two galaxies. The total mass in hot gas is (0.8--1.3) x 10^8 M_sun, which is approximately 1% of the total (HI~+~H2) gas mass in the bridge, and ~11% of the mass of warm molecular hydrogen discovered by Spitzer. The soft X-ray and dense CO-emitting gas in the bridge have offset distributions, with the X-rays peaking along the north-western side of the bridge in the region where stronger far-IR dust and diffuse [CII] gas is observed by Herschel. We detect nine Ultra Luminous X-ray sources (ULXs) in the system, the brightest of which is found in the bridge, associated with an extragalactic HII region. We suggest that the X-ray--emitting gas has been shocked--heated to high temperatures and "splashed" into the bridge by the collision. The large amount of gas stripped from the galaxies into the bridge and its very long gas depletion timescale (>10 Gyr) may explain why this system, unlike most major mergers, is not a powerful IR emitter.
Diffuse interstellar bands (DIBs) trace warm neutral and weakly-ionized diffuse interstellar medium (ISM). Here we present a dedicated, high signal-to-noise spectroscopic study of two of the strongest DIBs, at 5780 and 5797 \AA, in optical spectra of 666 early-type stars in the Small and Large Magellanic Clouds, along with measurements of the atomic Na\,{\sc i}\,D and Ca\,{\sc ii}\,K lines. The resulting maps show for the first time the distribution of DIB carriers across large swathes of galaxies, as well as the foreground Milky Way ISM. We confirm the association of the 5797 \AA\ DIB with neutral gas, and the 5780 \AA\ DIB with more translucent gas, generally tracing the star-forming regions within the Magellanic Clouds. Likewise, the Na\,{\sc i}\,D line traces the denser ISM whereas the Ca\,{\sc ii}\,K line traces the more diffuse, warmer gas. The Ca\,{\sc ii}\,K line has an additional component at $\sim200$--220 km s$^{-1}$ seen towards both Magellanic Clouds; this may be associated with a pan-Magellanic halo. Both the atomic lines and DIBs show sub-pc-scale structure in the Galactic foreground absorption; the 5780 and 5797 \AA\ DIBs show very little correlation on these small scales, as do the Ca\,{\sc ii}\,K and Na\,{\sc i}\,D lines. This suggests that good correlations between the 5780 and 5797 \AA\ DIBs, or between Ca\,{\sc ii}\,K and Na\,{\sc i}\,D, arise from the superposition of multiple interstellar structures. Similarity in behaviour between DIBs and Na\,{\sc i} in the SMC, LMC and Milky Way suggests the abundance of DIB carriers scales in proportion to metallicity.
A catalogue of the morphological features for the complete Spitzer Survey of Stellar Structure in Galaxies (S$^4$G), including 2352 nearby galaxies, is presented. The measurements are made using 3.6 $\mu$m images, largely tracing the old stellar population; at this wavelength the effects of dust are also minimal. The measured features are the sizes, ellipticities, and orientations of bars, rings, ringlenses, and lenses. Measured in a similar manner are also barlenses (lens-like structures embedded in the bars), which are not lenses in the usual sense, being rather the more face-on counterparts of the boxy/peanut structures in the edge-on view. In addition, pitch angles of spiral arm segments are measured for those galaxies where they can be reliably traced. More than one pitch angle may appear for a single galaxy. All measurements are made in a human-supervised manner so that attention is paid to each galaxy. We used isophotal analysis, unsharp masking, and fitting ellipses to measured structures. We find that the sizes of the inner rings and lenses normalized to barlength correlate with the galaxy mass: the normalized sizes increase toward the less massive galaxies; it has been suggested that this is related to the larger dark matter content in the bar region in these systems. Bars in the low mass galaxies are also less concentrated, likely to be connected to the mass cut-off in the appearance of the nuclear rings and lenses. We also show observational evidence that barlenses indeed form part of the bar, and that a large fraction of the inner lenses in the non-barred galaxies could be former barlenses in which the thin outer bar component has dissolved.
The surface density and vertical distribution of stars, stellar remnants, and gas in the solar vicinity form important ingredients for understanding the star formation history of the Galaxy as well as for inferring the local density of dark matter by using stellar kinematics to probe the gravitational potential. In this paper we review the literature for these baryonic components, reanalyze data, and provide tables of the surface densities and exponential scale heights of main sequence stars, giants, brown dwarfs, and stellar remnants. We also review three components of gas (H2, HI, and HII), give their surface densities at the solar circle, and discuss their vertical distribution. We find a local total surface density of M dwarfs of 17.3 pm 2.3 Mo/pc^2. Our result for the total local surface density of visible stars, 27.0 pm 2.7 Mo/pc^2, is close to previous estimates due to a cancellation of opposing effects: more mass in M dwarfs, less mass in the others. The total local surface density in white dwarfs is 4.9 pm 0.6 Mo/pc^2; in brown dwarfs, it is ~1.2 Mo/pc^2. We find that the total local surface density of stars and stellar remnants is 33.4 pm 3 Mo/pc^2, somewhat less than previous estimates. We analyze data on 21 cm emission and absorption and obtain good agreement with recent results on the local amount of neutral atomic hydrogen obtained with the Planck satellite. The local surface density of gas is 13.7 pm 1.6 Mo/pc^2. The total baryonic mass surface density that we derive for the solar neighborhood is 47.1 pm 3.4 Mo/pc^2. Combining these results with others' measurements of the total surface density of matter within 1-1.1 kpc of the plane, we find that the local density of dark matter is 0.013 pm 0.003Mo/pc^3.The local density of all matter is 0.097 pm 0.013 Mo/pc^3. We discuss limitations on the properties of a possible thin disk of dark matter.
We present a comprehensive and detailed study of the stellar populations of the Fornax dwarf spheroidal galaxy. We analyse their spatial distributions along the main body of the galaxy, obtaining their surface density maps, together with their radial density profiles. Results are based on the largest and most complete catalogue of stars in Fornax, with more than $3.5\times10^{5}$ stars covering the main body of the galaxy up to $\rm V \sim 24$. We find a differentiated structure in Fornax depending on the stellar ages. Old stars ($\gtrsim 10$ Gyr) follow an elliptical distribution well fitted by King profiles with relatively large core radius ($r_c = 760\pm60$ pc). On another hand, young populations ($\lesssim 3$ Gyr) concentrate in the central region of the galaxy ($r_c = 210\pm10$ pc), and are better fitted by Sersic profiles with $0.8 < n < 1.2$, indicating some disky shape. These stars show strong asymmetries and substructures not aligned with the main optical axes of Fornax. This together with the observed differences between metallicity and age distribution maps, strongly suggests accretion of material with different angular momentum. These results lead us to propose a scenario in which Fornax has suffered a major merger at $z\sim1$.
We have detected the Rossiter-Mclaughlin effect during a transit of WASP-47b, the only known hot Jupiter with close planetary companions. By combining our spectroscopic observations with Kepler photometry, we show that the projected stellar obliquity is $\lambda = 0^\circ \pm 24^\circ$. We can firmly exclude a retrograde orbit for WASP-47b, and rule out strongly misaligned prograde orbits. Low obliquities have also been found for most of the other compact multiplanet systems that have been investigated. The Kepler-56 system, with two close-in gas giants transiting their subgiant host star with an obliquity of at least 45$^\circ$, remains the only clear counterexample.
In this work we present and discuss the observations of the Mn abundances for 247 FGK dwarfs, located in the Galactic disc with metallicity -1<Fe/H]<+0.3. The observed stars belong to the substructures of the Galaxy thick and thin discs, and to the Hercules stream. The observations were conducted using the 1.93 m telescope at Observatoire de Haute-Provence (OHP, France) equipped with the echelle type spectrographs ELODIE and SOPHIE. The abundances were derived under the LTE approximation, with an average error for the [Mn/Fe] ratio of 0.10 dex. For most of the stars in the sample Mn abundances are not available in the literature. We obtain an evolution of [Mn/Fe] ratio with the metallicity [Fe/H] consistent with previous data compilations. In particular, within the metallicity range covered by our stellar sample the [Mn/Fe] ratio is increasing with the increase of metallicity. This due to the contribution to the Galactic chemical evolution of Mn and Fe from thermonuclear supernovae. We confirm the baseline scenario where most of the Mn in the Galactic disc and in the Sun is made by thermonuclear supernovae. In particular, the effective contribution from core-collapse supernovae to the Mn in the Solar system is about 10-20%. However, present uncertainties affecting the production of Mn and Fe in thermonuclear supernovae are limiting the constraining power of the observed [Mn/Fe] trend in the Galactic discs on, e.g., the frequency of different thermonuclear supernovae populations. The different production of these two elements in different types of thermonuclear supernovae needs to be disentangled by the dependence of their relative production on the metallicity of the supernova progenitor.
The detailed chemical composition of most metal-poor halo stars has been found to be highly uniform, but a minority of stars exhibit dramatic enhancements in their abundances of heavy neutron-capture elements and/or of carbon. The key question for Galactic chemical evolution models is whether these peculiarities reflect the composition of the natal clouds, or if they are due to later mass transfer of processed material from a binary companion. If the former case applies, the observed excess of certain elements was implanted within selected clouds in the early ISM from a production site at interstellar distances. Our aim is to determine the frequency and orbital properties of binaries among these chemically peculiar stars. This information provides the basis for deciding whether mass transfer from a binary companion is necessary and sufficient to explain their unusual compositions. This paper discusses our study of a sample of 17 moderately (r-I) and highly (r-II) r-process-element enhanced VMP and EMP stars. High-resolution, low signal-to-noise spectra of the stars were obtained at roughly monthly intervals over 8 years with the FIES spectrograph at the Nordic Optical Telescope. From these spectra, radial velocities with an accuracy of ~100 m/s were determined by cross-correlation against an optimized template. 14 of the programme stars exhibit no significant RV variation over this period, while 3 are binaries with orbits of typical eccentricity for their periods, resulting in a normal binary frequency of ~18+-6% for the sample. Our results confirm our preliminary conclusion from 2011, based on partial data, that the chemical peculiarity of the r-I and r-II stars is not caused by any putative binary companions. Instead, it was imprinted on the natal molecular clouds of these stars by an external, distant source. Models of the ISM in early galaxies should account for such mechanisms.
Though the small-scale, low-ionization knots, filaments and jets (LISs) of planetary nebulae (PNe) are known for ~30yr, some of their observational properties are not well established. In consequence our ability to include them in the wider context of the formation and evolution of PNe is directly affected. Why most structures have lower densities than the PN shells hosting them? Is their intense emission in low-ionization lines the key to their main excitation mechanism? Therefore, if considered altogether, can LISs line ratios, chemical abundances and kinematics enlighten the interplay between the different excitation and formation processes? Here we present a spectroscopic analysis of five PNe that possess LISs confirming that all nebular components have comparable electron temperatures, whereas the electron density is systematically lower in LISs than in the surrounding nebula. Chemical abundances of LISs versus other PN components do not show significant differences as well. By using diagnostic diagrams from shock models, we demonstrate that LISs' main excitation is due to shocks, whereas the other components are mainly photo-ionized. We also propose new diagnostic diagrams involving a few emission lines ([NII], [OIII], [SII]) and $\rm{log}$(f$_{shocks}$/f$_{\star}$), where f$_{shocks}$ and f$_{\star}$ are the ionization photon fluxes due to the shocks and the central star ionizing continuum, respectively. A robust relation differentiating the structures is found, with the shock-excited clearly having $\rm{log}$(f$_{shocks}$/f$_{\star}$)$>$-1; while the photo-ionized show $\rm{log}$(f$_{shocks}$/f$_{\star}$)$<$-2. A transition zone, with -2$<\rm{log}$(f$_{shocks}$/f$_{\star}$)$<$-1 where both mechanisms are equally important, is also defined.
Our goal is to relate the photometric and spectroscopic variability of classical T Tauri stars, of the star-forming cluster NGC 2264, to the physical processes acting in the stellar and circumstellar environment, within a few stellar radii from the star. NGC 2264 was the target of a multiwavelength observational campaign with CoRoT, MOST, Spitzer, and Chandra satellites and observations from the ground. We classified the CoRoT light curves of accreting systems according to their morphology and compared our classification to several accretion diagnostics and disk parameters. The morphology of the CoRoT light curve reflects the evolution of the accretion process and of the inner disk region. Accretion burst stars present high mass-accretion rates and optically thick inner disks. AA Tau-like systems, whose light curves are dominated by circumstellar dust obscuration, show intermediate mass-accretion rates and are located in the transition of thick to anemic disks. Classical T Tauri stars with spot-like light curves correspond mostly to systems with a low mass-accretion rate and low mid-IR excess. About 30% of the classical T Tauri stars observed in the 2008 and 2011 CoRoT runs changed their light-curve morphology. Transitions from AA Tau-like and spot-like to aperiodic light curves and vice versa were common. The analysis of the $H\alpha$ emission line variability of 58 accreting stars showed that 8 presented a periodicity that in a few cases was coincident with the photometric period. The blue and red wings of the $H\alpha$ line profiles often do not correlate with each other, indicating that they are strongly influenced by different physical processes. Accreting stars have a dynamic stellar and circumstellar environment that can be explained by magnetospheric accretion and outflow models, including variations from stable to unstable accretion regimes on timescales of a few years
Upcoming or future deep galaxy samples with wide sky coverage can provide independent measurement of the kinematic dipole - our motion relative to the rest frame defined by the large-scale structure. Such a measurement would present an important test of the standard cosmological model, as the standard model predicts the galaxy measurement should precisely agree with the existing precise measurements made using the CMB. However, the required statistical precision to measure the kinematic dipole typically makes the measurement susceptible to bias from the presence of the local-structure-induced dipole contamination. In order to minimize the latter, a sufficiently deep survey is required. We forecast both the statistical error and the systematic bias in the kinematic dipole measurements. We find that a survey covering $\sim 75\%$ of the sky in both hemispheres and having $\sim 30$ million galaxies can detect the kinematic dipole at $5\sigma$, while its median redshift should be at least $z_{med} \sim 0.75$ for negligible bias from the local structure.
We review results from our monitoring observations of several lensed quasars performed in the optical, UV, and X-ray bands. Modeling of the multi-wavelength light curves provides constraints on the extent of the optical, UV, and X-ray emission regions. One of the important results of our analysis is that the optical sizes as inferred from the microlensing analysis are significantly larger than those predicted by the theoretical-thin-disk estimate. In a few cases we also constrain the slope of the size-wavelength relation. Our size constraints of the soft and hard X-ray emission regions of quasars indicate that in some objects of our sample the hard X-ray emission region is more compact than the soft and in others the soft emission region is smaller. This difference may be the result of the relative strengths of the disk-reflected (harder and extended) versus corona-direct (softer and compact) components in the quasars of our sample. Finally, we present the analysis of several strong microlensing events where we detect an evolution of the relativistic Fe line profile as the magnification caustic traverses the accretion disk. These caustic crossings are used to provide constraints on the innermost stable circular orbit (ISCO) radius and the accretion disk inclination angle of the black hole in quasar RX J1131-1231.
We present photometric and spectroscopic observations of SN 2013fc, a bright type II supernova (SN) in a circumnuclear star-forming ring in the luminous infrared galaxy ESO 154-G010, observed as part of the Public ESO Spectroscopic Survey of Transient Objects (PESSTO). SN 2013fc is both photometrically and spectroscopically similar to the well-studied type IIn SN 1998S and to the bright type II-L SN 1979C. It exhibits an initial linear decline after maximum, followed by a short plateau phase and a tail phase with a decline too fast for $^{56}$Co decay with full gamma-ray trapping. Initially the spectrum was blue and featureless. Later on, a strong broad ($\sim 8000$ km s$^{-1}$) H$\alpha$ emission profile became prominent. We apply a Starlight stellar population model fit to the SN location (observed when the SN had faded) to estimate both a high extinction of $A_V = 2.9 \pm 0.2$ mag and an age of $10_{-2}^{+3}$ Myr for the underlying cluster. We compare the SN to SNe 1998S and 1979C and discuss its possible progenitor star in light of the similarities to these events. With a peak brightness of $B = -20.46 \pm 0.21$ mag, SN 2013fc is 0.9 mag brighter than SN 1998S and of comparable brightness to SN 1979C. We suggest that SN 2013fc, like SNe 1998S and 1979C, was consistent with a massive red supergiant (RSG) progenitor. Recent mass loss probably due to a strong RSG wind created the circumstellar matter illuminated through its interaction with the SN ejecta. We also observe a near-infrared excess, possibly due to newly condensed dust.
The B-mode Foreground Experiment (BFORE) is a proposed NASA balloon project designed to make optimal use of the sub-orbital platform by concentrating on three dust foreground bands (270, 350, and 600 GHz) that complement ground-based cosmic microwave background (CMB) programs. BFORE will survey ~1/4 of the sky with 1.7 - 3.7 arcminute resolution, enabling precise characterization of the Galactic dust that now limits constraints on inflation from CMB B-mode polarization measurements. In addition, BFORE's combination of frequency coverage, large survey area, and angular resolution enables science far beyond the critical goal of measuring foregrounds. BFORE will constrain the velocities of thousands of galaxy clusters, provide a new window on the cosmic infrared background, and probe magnetic fields in the interstellar medium. We review the BFORE science case, timeline, and instrument design, which is based on a compact off-axis telescope coupled to >10,000 superconducting detectors.
Models of the dynamical evolution of the early solar system following the dispersal of the gaseous protoplanetary disk have been widely successful in reconstructing the current orbital configuration of the giant planets. Statistically, some of the most successful dynamical evolution simulations have initially included a hypothetical fifth giant planet, of ice giant mass, which gets ejected by a gas giant during the early solar system's proposed instability phase. We investigate the likelihood of an ice giant ejection event by either Jupiter or Saturn through constraints imposed by the current orbits of their wide-separation regular satellites Callisto and Iapetus respectively. We show that planetary encounters that are sufficient to eject an ice giant, often provide excessive perturbations to the orbits of Callisto and Iapetus making it difficult to reconcile a planet ejection event with the current orbit of either satellite. Quantitatively, we compute the likelihood of reconciling a regular Jovian satellite orbit with the current orbit of Callisto following an ice giant ejection by Jupiter of ~ 42% and conclude that such a large likelihood supports the hypothesis of a fifth giant planet's existence. A similar calculation for Iapetus reveals that it is much more difficult for Saturn to have ejected an ice giant and reconcile a Kronian satellite orbit with that of Iapetus (likelihood ~ 1%), although uncertainties regarding the formation of Iapetus, on its unusual orbit, complicates the interpretation of this result.
We present the first radiation magnetohydrodynamics simulations of the atmosphere of white dwarf stars. We demonstrate that convective energy transfer is seriously impeded by magnetic fields when the plasma-beta parameter, the thermal to magnetic pressure ratio, becomes smaller than unity. The critical field strength that inhibits convection in the photosphere of white dwarfs is in the range B = 1-50 kG, which is much smaller than the typical 1-1000 MG field strengths observed in magnetic white dwarfs, implying that these objects have radiative atmospheres. We have then employed evolutionary models to study the cooling process of high-field magnetic white dwarfs, where convection is entirely suppressed during the full evolution (B > 10 MG). We find that the inhibition of convection has no effect on cooling rates until the effective temperature (Teff) reaches a value of around 5500 K. In this regime, the standard convective sequences start to deviate from the ones without convection owing to the convective coupling between the outer layers and the degenerate reservoir of thermal energy. Since no magnetic white dwarfs are currently known at the low temperatures where this coupling significantly changes the evolution, effects of magnetism on cooling rates are not expected to be observed. This result contrasts with a recent suggestion that magnetic white dwarfs with Teff < 10,000 K cool significantly slower than non-magnetic degenerates.
Based on the relationship between thermodynamics and gravity, and with the aid of Verlinde's formalism, we propose an alternative interpretation of the dynamical evolution of the Friedmann-Robertson-Walker Universe, which takes into account the entropy and temperature intrinsic to the horizon of the universe due to the information holographically stored there through non-gaussian statistical theories proposed by Tsallis and Kaniadakis. We use the most recent data of type Ia supernovae, baryon acoustic oscillations, and the Hubble expansion rate function to constrain the free parameters on the $\Lambda$CDM and $w$CDM models modified by the non-gaussian statistics. We evaluate the problem of age and we note that such modifications solve the problem at 1$\sigma$ level confidence. Also we analyze the effects on the linear growth of matter density perturbations.
The pseudoscalar resonance or "A-funnel" in the Minimal Supersymmetric Standard Model~(MSSM) is a widely studied framework for explaining dark matter that can yield interesting indirect detection and collider signals. The well-known Galactic Center excess (GCE) at GeV energies in the gamma ray spectrum, consistent with annihilation of a $\lesssim 40$ GeV dark matter particle, has more recently been shown to be compatible with significantly heavier masses following reanalysis of the background. In this paper, we explore the LHC and direct detection implications of interpreting the GCE in this extended mass window within the MSSM A-funnel framework. We find that compatibility with relic density, signal strength, collider constraints, and Higgs data can be simultaneously achieved with appropriate parameter choices. The compatible regions give very sharp predictions of 200-600 GeV CP-odd/even Higgs bosons at low tan$\beta$ at the LHC and spin-independent cross sections $\approx 10^{-11}$ pb at direct detection experiments. Regardless of consistency with the GCE, this study serves as a useful template of the strong correlations between indirect, direct, and LHC signatures of the MSSM A-funnel region.
Recently, there have been claims in the literature that the cosmological constant problem can be dynamically solved by specific compactifications of gravity from higher-dimensional toy models. These models have the novel feature that in the four-dimensional theory, the cosmological constant $\Lambda$ is much smaller than the Planck density and in fact accumulates at $\Lambda=0$. Here we show that while these are very interesting models, they do not properly address the real cosmological constant problem. As we explain, the real problem is not simply to obtain $\Lambda$ that is small in Planck units in a toy model, but to explain why $\Lambda$ is much smaller than other mass scales (and combinations of scales) in the theory. Instead, in these toy models, all other particle mass scales have been either removed or sent to zero, thus ignoring the real problem. To this end, we provide a general argument that the included moduli masses are generically of order Hubble, so sending them to zero trivially sends the cosmological constant to zero. We also show that the fundamental Planck mass is being sent to zero, and so the central problem is trivially avoided by removing high energy physics altogether. On the other hand, by including various large mass scales from particle physics with a high fundamental Planck mass, one is faced with a real problem, whose only known solution involves accidental cancellations in a landscape.
Combining the covariant coalescence model and a blast-wave-like analytical parametrization for (anti-)nucleon phase-space freezeout configuration, we explore light (anti-)nucleus production in central Au+Au collisions at $\sqrt{s_{NN}} = 200$ GeV. Using the nucleon freezeout configuration (denoted by FO1) determined from the measured spectra of protons (p), deutrons (d) and $^{3}$He, we find the predicted yield of $^{4}$He is significantly smaller than the experimental data. We show this disagreement can be removed by using a nucleon freezeout configuration (denoted by FO2) in which the nucleons are assumed to freeze out earlier than those in FO1 to effectively consider the effect of large binding energy value of $^{4}$He. Assuming the binding energy effect also exists for the production of $^5\text{Li}$, $^5\overline{\text{Li}}$, $^6\text{Li}$ and $^6\overline{\text{Li}}$ due to their similar binding energy values as $^{4}$He, we find the yields of these heavier (anti-)nuclei can be enhanced by a factor of about one order, implying that although the stable (anti-)$^6$Li nucleus is unlikely to be observed, the unstable (anti-)$^5$Li nucleus could be produced in observable abundance in Au+Au collisions at $\sqrt{s_{NN}} = 200$ GeV where it may be identified through the p-$^4\text{He}$ ($\overline{\text{p}}$-$^4\overline{\text{He}}$) invariant mass spectrum. The future experimental measurement on (anti-)$^5\text{Li}$ would be very useful to understand the production mechanism of heavier antimatter.
Collisionless shocks occur in various fields of physics. In the context of space and astrophysics they have been investigated for many decades. However, a thorough understanding of shock formation and particle acceleration is still missing. Collisionless shocks can be distinguished into electromagnetic and electrostatic shocks. Electromagnetic shocks are of importance mainly in astrophysical environments and they are mediated by the Weibel or filamentation instability. In such shocks, charged particles gain energy by diffusive shock acceleration. Electrostatic shocks are characterized by a strong electrostatic field, which leads to electron trapping. Ions are accelerated by reflection from the electrostatic potential. Shock formation and particle acceleration will be discussed in theory and simulations.
Inspired by the recent diboson excess observed at the LHC and possible interpretation within a TeV-scale Left-Right symmetric framework, we explore its implications for low-energy experiments searching for lepton number and flavor violation. Assuming a simple Type-II seesaw mechanism for neutrino masses, we show that for the right-handed (RH) gauge boson mass and coupling values required to explain the LHC anomalies, the RH contribution to the lepton number violating process of neutrinoless double beta decay ($0\nu\beta\beta$) is already constrained by current experiments for relatively low-mass (MeV-GeV) RH neutrinos. The future ton-scale $0\nu\beta\beta$ experiments could probe most of the remaining parameter space, irrespective of the neutrino mass hierarchy and uncertainties in the oscillation parameters and nuclear matrix elements. On the other hand, the RH contribution to the lepton flavor violating process of $\mu\to e\gamma$ is constrained for relatively heavier (TeV) RH neutrinos, thus providing a complementary probe of the model. Finally, a measurement of the absolute light neutrino mass scale from future precision cosmology could make this scenario completely testable.
Links to: arXiv, form interface, find, astro-ph, recent, 1509, contact, help (Access key information)