Links to: arXiv, form interface, find, astro-ph, recent, 1711, contact, help (Access key information)
We propose a novel method utilizing stellar kinematic data to detect low-mass substructure in the Milky Way's dark matter halo. By probing characteristic wakes that a passing dark matter subhalo leaves in the phase space distribution of ambient halo stars, we estimate sensitivities down to subhalo masses $\sim 10^7\,M_\odot$ or below. The detection of such subhalos would have implications for dark-matter and cosmological models that predict modifications to the halo-mass function at low halo masses. We develop an analytic formalism for describing the perturbed stellar phase-space distributions, and we demonstrate through simulations the ability to detect subhalos using the phase-space model and a likelihood framework. Our method complements existing methods for low-mass subhalo searches, such as searches for gaps in stellar streams, in that we can localize the positions and velocities of the subhalos today.
We use KMOS Deep Survey (KDS) galaxies, combined with results from a range of spectroscopic studies in the literature, to investigate the evolution of the stellar-mass Tully-Fisher relation since z ~ 4. We determine the slope and normalisation of the local rotation-velocity -- stellar-mass (Vc - $M_{\star}$) relationship using a reference sample of local spiral galaxies; thereafter we fix the slope, and focus on the evolution of velocity normalisation with redshift. The rotation-dominated KDS galaxies at z ~ 3.5 have rotation velocities ~ -0.1 dex lower than local reference galaxies at fixed stellar mass. By fitting 16 distant comparison samples spanning 0 < z < 3 (containing ~ 1200 galaxies), we show that the size and sign of the inferred Vc offset depends sensitively on the fraction of the parent samples used in the Tully-Fisher analysis, and how strictly the criterion of 'rotation dominated' is enforced. Confining attention to subsamples of galaxies that are especially 'disky' results in a consistent positive offset in Vc of ~ +0.1 dex, however these galaxies are not representative of the evolving-disk population at z > 1. We investigate the addition of pressure support, traced by intrinsic velocity dispersion ($\sigma_{int}$) to the KDS dynamical mass budget by adopting a 'total' effective velocity of form $V_{tot} = (Vc^{2} + 4.0\sigma_{int}^{2})^{0.5}$. The rotation-dominated and dispersion-dominated KDS galaxies fall on the same locus in the total-velocity versus stellar-mass plane, removing the need for debate over the precise selection threshold for rotation-dominated galaxies. The comparison sample offsets are in the range +0.08 to +0.15 dex in total-velocity zero-point (-0.30 to -0.55 dex in stellar-mass zero-point) from the local Tully-Fisher relation at z > 1, consistent with steady evolution of the ratio of dynamical to stellar mass with cosmic time.
Using data from MALT90 (Millimetre Astronomy Legacy Team Survey at 90 GHz), we present molecular line study of a sample of ATLASGAL (APEX Telescope Large Area Survey of the Galaxy) clumps. Twelve emission lines have been detected in all. We found that in most sources, emissions of HC$_3$N, HN$^{13}$C, CH$_3$CN, HNCO and SiO show more compact distributions than those of HCO$^+$, HNC, HCN and N$_2$H$^+$. By comparing with other molecular lines, we found that the abundance of HNCO ($\chi$(HNCO)) correlates well with other species such as HC3N, HNC, C2H, H13CO+ and N2H+. Previous studies indicate the HNCO abundance could be enhanced by shocks. However, in this study, we found the abundance of HNCO does not correlate well with that of SiO, which is also a good tracer of shocks. We suggest this may be because HNCO and SiO trace different parts of shocks. Our analysis indicates that the velocity of shock traced by HNCO tends to be lower than that traced by SiO. In the low-velocity shocks traced by HNCO, the HNCO abundance increases faster than that of SiO. While in the relatively high-velocity shocks traced by SiO, the SiO abundance increases faster than that of HNCO. We suggest that in the infrared dark cloud (IRDC) of MSXDC G331.71+00.59, high-velocity shocks are destroying the molecule of HNCO.
In this Letter, we present an analysis of the relation between the variability of broad absorption lines (BALs) and that of the continuum. Our sample is multi-epoch observations of 483 quasars by the Sloan Digital Sky Survey-I/II/III (SDSS-I/II/III). We derive the fractional flux variations of the continuum and fractional equivalent width (EW) variations for C IV and Si IV BALs, and explore the correlations between the three. Our results reveal moderate anticorrelations with high significance level between the fractional flux variations of the continuum and fractional EW variations for both C IV and Si IV BALs. We also prove a significant positive correlation between the fractional EW variations for C IV and Si IV BALs, which is in agreement with several previous studies. Our discoveries can serve as evidence for the idea: Change of an ionizing continuum is the primary driver of BAL variability.
We investigate emission lines of Lya, [OIII]5007, Ha, and [CII]158um from 1125 galaxies at z=4.9-7.0. Our sample is composed of 1098 Lya emitters (LAEs) at z=4.9, 5.7, 6.6, and 7.0 identified by Subaru/Hyper Suprime-Cam (HSC) narrowband surveys in the 4deg2 area covered by Spitzer large area survey with Subaru/HSC (SPLASH) and 29 galaxies at z=5.148-7.508 with deep [CII] ALMA or PdBI data in the literature. Fluxes of strong rest-frame optical emission lines of [OIII] and Ha (Hb) are constrained by significant excesses found in the SPLASH 3.6 and 4.5um photometry. At z=4.9, the rest-frame Ha equivalent width, EW^0_Ha, and the Lya escape fraction, f_Lya, positively correlate with the rest-frame Lya equivalent width, EW_Lya. The f_Lya-EW_Lya correlation is similarly found at z~0-2, suggesting no evolution of the correlation over z~0-5. The typical ionizing photon production efficiency of LAEs is logxi_ion/[Hz erg^-1]~25.5, significantly (60-100%) higher than those of LBGs at a given UV magnitude. At z=5.7-7.0, there exists an interesting turn-over trend that the [OIII]/Ha flux ratio increases in EW_Lya~0-30A, where it reaches a maximum, and then decreases out to EW_Lya~130A. We also identify a negative correlation between a [CII] luminosity to star-formation rate ratio (L_[CII]/SFR) and EW_Lya at the 4.5sigma confidence level. We carefully investigate the physical origins of the trend and the correlations with stellar-synthesis and photoionization models covering the vast parameter space of metallicity, ionization parameter, and stellar age, and find that the [OIII]/Ha turn-over trend is explained by a simple anti-correlation between EW_Lya and metallicity, indicative of the detections of the very metal-poor (~0.04Zo) galaxies with EW_Lya~200A. Moreover, this simple anti-correlation explains self-consistently all of the correlations of Lya, Ha, [OIII]/Ha, and [CII] identified in our study.
Recently a theory about the formation of over-densities of stars along tidal tails of globular clusters has been presented, this theory predicts the position and time of formation of such over-densities and was successfully tested with N-body simulations of globular clusters in a point mass galactic potential. In this work we present a comparison between this theory and our simulations using a dwarf galaxy orbiting two differently shaped dark matter halos to study the effects of a cored and a cuspy halo on the formation and evolution of tidal tails. We find no difference using a cuspy or a cored halo, however, we find an intriguing asymmetry between the leading and trailing arm of the tidal tails. The trailing arm grows faster than the leading arm. This asymmetry is seen in the distance to first over-density and its size as well. We establish a relation between the distance to the first over-density and the size of this over-density.
We present an analysis of the average spectral properties of 12,000 SDSS quasars as a function of accretion disc inclination, as measured from the equivalent width of the [O III] 5007{\AA} line. The use of this indicator on a large sample of quasars from the SDSS DR7 has proven the presence of orientation effects on the features of UV/optical spectra, confirming the presence of outflows in the NLR gas and that the geometry of the BLR is disc-like. Relying on the goodness of this indicator, we are now using it to investigate other bands/components of AGN. Specifically, the study of the UV/optical/IR SED of the same sample provides information on the obscuring "torus". The SED shows a decrease of the IR fraction moving from face-on to edge-on sources, in agreement with models where the torus is co-axial with the accretion disc. Moreover, the fact we are able to observe the broad emission lines also in sources in an edge-on position, suggests that the torus is rather clumpy than smooth as in the Unified Model. The behaviour of the SED as a function of EW[OIII] is in agreement with the predictions of the clumpy torus models as well.
We study the spiral arm morphology of a sample of the local spiral galaxies in the Illustris simulation and explore the supermassive black hole - galaxy connection beyond the bulge (e.g., spiral arm pitch angle, total stellar mass, dark matter mass, and total halo mass), finding good agreement with other theoretical studies and observational constraints. It is important to study the properties of supermassive black holes and their host galaxies through both observations and simulations and compare their results in order to understand their physics and formative histories.We find that Illustris prediction for supermassive black hole mass relative to pitch angle is in rather good agreement with observations and that barred and non-barred galaxies follow similar scaling relations. Our work shows that Illustris presents very tight correlations between supermassive black hole mass and large-scale properties of the host galaxy, not only for early-type galaxies but also low-mass, blue and star-forming galaxies. These tight relations beyond the bulge suggest that halo properties determine those of a disc galaxy and its supermassive black hole.
A more fine-tuned method for probing planet-forming regions, such as protoplanetary discs, could be rovibrational molecular spectroscopy observation of particular premineral molecules instead of more common but ultimately less related volatile organic compounds. Planets are created when grains aggregate, but how molecules form grains is an ongoing topic of discussion in astrophysics and planetary science. Using the spectroscopic data of molecules specifically involved in mineral formation could help to map regions where planet formation is believed to be occurring in order to examine the interplay between gas and dust. Four atoms are frequently associated with planetary formation: Fe, Si, Mg, and O. Magnesium, in particular, has been shown to be in higher relative abundance in planet-hosting stars. Magnesium oxide crystals comprise the mineral periclase making it the chemically simplest magnesium-bearing mineral and a natural choice for analysis. The monomer, dimer, and trimer forms of (MgO)_n with n = 1 - 3 are analyzed in this work using high-level quantum chemical computations known to produce accurate results. Strong vibrational transitions at 12.5 {\mu}m, 15.0 {\mu}m, and 16.5 {\mu}m are indicative of magnesium oxide monomer, dimer, and trimer making these wavelengths of particular interest for the observation of protoplanetary discs and even potentially planet-forming regions around stars. If such transitions are observed in emission from the accretion discs or absorptions from stellar spectra, the beginning stages of mineral and, subsequently, rocky body formation could be indicated.
We present a novel strategy to uncover the Galactic population of quiescent black holes (BHs). This is based on a new concept, the photometric mass function (PMF), which opens up the possibility of an efficient identification of dynamical BHs in large fields-of-view. This exploits the width of the disc Halpha emission line, combined with orbital period information. We here show that Halpha widths can be recovered using a combination of customized Halpha filters. By setting a width cut-off at 2200 km/s we are able to cleanly remove other Galactic populations of Halpha emitters, including ~99.9% of cataclysmic variables (CVs). Only short period (Porb<2.1 h) eclipsing CVs and AGNs will contaminate the sample but these can be easily flagged through photometric variability and, in the latter case, also mid-IR colours. We also describe the strategy of a deep (r=22) Galactic plane survey based on the concept of PMFs: HAWKs, the HAlpha-Width Kilo-deg Survey. We estimate that ~800 sqr deg are required to unveil ~50 new dynamical BHs, a three-fold improvement over the known population. For comparison, a century would be needed to produce an enlarged sample of 50 dynamical BHs from X-ray transients at the current discovery rate.
The first four gravitational wave events detected by LIGO were all interpreted as merging black hole binaries (BHBs), opening a new perspective on the study of such systems. Here we use our new population-synthesis code MOBSE, an upgraded version of BSE (Hurley et al. 2002), to investigate the demography of merging BHBs. MOBSE includes metallicity-dependent prescriptions for mass loss of massive hot stars. It also accounts for the impact of the electron-scattering Eddington factor on mass loss. We perform >10^8 simulations of isolated massive binaries, with 12 different metallicities, to study the impact of mass loss, core-collapse supernovae and common envelope on merging BHBs. Accounting for the dependence of stellar winds on the Eddington factor leads to the formation of black holes (BHs) with mass up to 65 Msun at metallicity Z~0.0002. However, most BHs in merging BHBs have masses <40 Msun. We find merging BHBs with mass ratios in the 0.1 - 1.0 range, even if mass ratios >0.6 are more likely. We predict that systems like GW150914, GW170814 and GW170104 can form only from progenitors with metallicity Z<=0.006, Z<=0.008 and Z<=0.012, respectively. Most merging BHBs have gone through a common envelope phase, but up to ~17 per cent merging BHBs at low metallicity did not undergo any common envelope phase. We find a much higher number of mergers from metal-poor progenitors than from metal-rich ones: the number of BHB mergers per unit mass is ~10^-4 Msun^-1 at low metallicity (Z = 0.0002 - 0.002) and drops to ~10^-7 Msun^-1 at high metallicity (Z ~ 0.02).
Important insights into the formation and evolution of the Galactic disc(s) are contained in the chemical compositions of stars. We analysed high-resolution and high signal to noise HARPS spectra of 79 solar twin stars in order to obtain precise determinations of their atmospheric parameters, ages ($\sigma$$\sim$0.4 Gyr) and chemical abundances ($\sigma$$<$0.01~dex) of 12 neutron-capture elements (Sr, Y, Zr, Ba, La, Ce, Pr, Nd, Sm, Eu, Gd, and Dy). This valuable dataset allows us to study the [X/Fe]-age relations over a time interval of $\sim$10 Gyr and among stars belonging to the thin and thick discs. These relations show that i) the $s$-process has been the main channel of nucleosynthesis of $n$-capture elements during the evolution of the thin disc; ii) the thick disc is rich in $r$-process elements which suggests that its formation has been rapid and intensive. %; iii) a chemical continuity between the thin and thick discs is evident in the abundances of Ba. In addition, the heavy (Ba, La, Ce) and light (Sr, Y, Zr) $s$-process elements revealed details on the dependence between the yields of AGB stars and the stellar mass or metallicity. Finally, we confirmed that both [Y/Mg] and [Y/Al] ratios can be employed as stellar clocks, allowing ages of solar twin stars to be estimated with an average precision of $\sim$0.5~Gyr.
We use the first data release of the Gaia mission to explore the three dimensional arrangement and the age ordering of the many stellar groups towards the Orion OB association, aiming at a new classification and characterization of the stellar population. We make use of the parallaxes and proper motions provided in the Tycho Gaia Astrometric Solution (TGAS) sub-set of the Gaia catalogue, and of the combination of Gaia and 2MASS photometry. In TGAS we find evidence for the presence of a young population, at a parallax $\varpi \sim 2.65 \, \mathrm{mas}$, loosely distributed around some known clusters: 25 Ori, $\epsilon$ Ori and $\sigma$ Ori, and NGC 1980 ($\iota$ Ori). The low mass counterpart of this population is visible in the color-magnitude diagrams constructed by combining Gaia and 2MASS photometry. We study the density distribution of the young sources in the sky. We find the same groups as in TGAS, and also some other density enhancements that might be related to the recently discovered Orion X group, the Orion dust ring, and to the $\lambda$ Ori complex. We estimate the ages of this population and we infer the presence of an age gradient going from 25 Ori (13-15 Myr) to the ONC (1-2 Myr). We confirm this age ordering by repeating the Bayesian fit using the Pan-STARRS1 data. The estimated ages towards the NGC 1980 cluster span a broad range of values. This can either be due to the presence of two populations coming from two different episodes of star formation or to a large spread along the line of sight of the same population. Our results form the first step towards using the Gaia data to unravel the complex star formation history of the Orion region in terms of the different star formation episodes, their duration, and their effects on the surrounding interstellar medium.
In light of the observed Galactic center gamma-ray excess, we investigate a simplified model, for which the scalar dark matter interacts with quarks through a pseudoscalar mediator. The viable regions of the parameter space, that can also account for the relic density and evade the current searches, are identified, if the low-velocity dark matter annihilates through an $s$-channel off-shell mediator mostly into $\bar{b} b$, and/or annihilates directly into two ${\it hidden}$ on-shell mediators, which subsequently decay into the quark pairs. These two kinds of annihilations are $s$-wave. The projected monojet limit set by the high luminosity LHC sensitivity could constrain the favored parameter space, where the mediator's mass is larger than the dark matter mass by a factor of 2. We show that the projected sensitivity of 15-year Fermi-LAT observations of dwarf spheroidal galaxies can provide a stringent constraint on the most parameter space allowed in this model. If the on-shell mediator channel contributes to the dark matter annihilation cross sections over 50$\%$, this model with a lighter mediator can be probed in the projected PICO-500L experiment.
Links to: arXiv, form interface, find, astro-ph, recent, 1711, contact, help (Access key information)
We propose a novel method utilizing stellar kinematic data to detect low-mass substructure in the Milky Way's dark matter halo. By probing characteristic wakes that a passing dark matter subhalo leaves in the phase space distribution of ambient halo stars, we estimate sensitivities down to subhalo masses $\sim 10^7\,M_\odot$ or below. The detection of such subhalos would have implications for dark-matter and cosmological models that predict modifications to the halo-mass function at low halo masses. We develop an analytic formalism for describing the perturbed stellar phase-space distributions, and we demonstrate through simulations the ability to detect subhalos using the phase-space model and a likelihood framework. Our method complements existing methods for low-mass subhalo searches, such as searches for gaps in stellar streams, in that we can localize the positions and velocities of the subhalos today.
We use KMOS Deep Survey (KDS) galaxies, combined with results from a range of spectroscopic studies in the literature, to investigate the evolution of the stellar-mass Tully-Fisher relation since z ~ 4. We determine the slope and normalisation of the local rotation-velocity -- stellar-mass (Vc - $M_{\star}$) relationship using a reference sample of local spiral galaxies; thereafter we fix the slope, and focus on the evolution of velocity normalisation with redshift. The rotation-dominated KDS galaxies at z ~ 3.5 have rotation velocities ~ -0.1 dex lower than local reference galaxies at fixed stellar mass. By fitting 16 distant comparison samples spanning 0 < z < 3 (containing ~ 1200 galaxies), we show that the size and sign of the inferred Vc offset depends sensitively on the fraction of the parent samples used in the Tully-Fisher analysis, and how strictly the criterion of 'rotation dominated' is enforced. Confining attention to subsamples of galaxies that are especially 'disky' results in a consistent positive offset in Vc of ~ +0.1 dex, however these galaxies are not representative of the evolving-disk population at z > 1. We investigate the addition of pressure support, traced by intrinsic velocity dispersion ($\sigma_{int}$) to the KDS dynamical mass budget by adopting a 'total' effective velocity of form $V_{tot} = (Vc^{2} + 4.0\sigma_{int}^{2})^{0.5}$. The rotation-dominated and dispersion-dominated KDS galaxies fall on the same locus in the total-velocity versus stellar-mass plane, removing the need for debate over the precise selection threshold for rotation-dominated galaxies. The comparison sample offsets are in the range +0.08 to +0.15 dex in total-velocity zero-point (-0.30 to -0.55 dex in stellar-mass zero-point) from the local Tully-Fisher relation at z > 1, consistent with steady evolution of the ratio of dynamical to stellar mass with cosmic time.
Using data from MALT90 (Millimetre Astronomy Legacy Team Survey at 90 GHz), we present molecular line study of a sample of ATLASGAL (APEX Telescope Large Area Survey of the Galaxy) clumps. Twelve emission lines have been detected in all. We found that in most sources, emissions of HC$_3$N, HN$^{13}$C, CH$_3$CN, HNCO and SiO show more compact distributions than those of HCO$^+$, HNC, HCN and N$_2$H$^+$. By comparing with other molecular lines, we found that the abundance of HNCO ($\chi$(HNCO)) correlates well with other species such as HC3N, HNC, C2H, H13CO+ and N2H+. Previous studies indicate the HNCO abundance could be enhanced by shocks. However, in this study, we found the abundance of HNCO does not correlate well with that of SiO, which is also a good tracer of shocks. We suggest this may be because HNCO and SiO trace different parts of shocks. Our analysis indicates that the velocity of shock traced by HNCO tends to be lower than that traced by SiO. In the low-velocity shocks traced by HNCO, the HNCO abundance increases faster than that of SiO. While in the relatively high-velocity shocks traced by SiO, the SiO abundance increases faster than that of HNCO. We suggest that in the infrared dark cloud (IRDC) of MSXDC G331.71+00.59, high-velocity shocks are destroying the molecule of HNCO.
In this Letter, we present an analysis of the relation between the variability of broad absorption lines (BALs) and that of the continuum. Our sample is multi-epoch observations of 483 quasars by the Sloan Digital Sky Survey-I/II/III (SDSS-I/II/III). We derive the fractional flux variations of the continuum and fractional equivalent width (EW) variations for C IV and Si IV BALs, and explore the correlations between the three. Our results reveal moderate anticorrelations with high significance level between the fractional flux variations of the continuum and fractional EW variations for both C IV and Si IV BALs. We also prove a significant positive correlation between the fractional EW variations for C IV and Si IV BALs, which is in agreement with several previous studies. Our discoveries can serve as evidence for the idea: Change of an ionizing continuum is the primary driver of BAL variability.
We investigate emission lines of Lya, [OIII]5007, Ha, and [CII]158um from 1125 galaxies at z=4.9-7.0. Our sample is composed of 1098 Lya emitters (LAEs) at z=4.9, 5.7, 6.6, and 7.0 identified by Subaru/Hyper Suprime-Cam (HSC) narrowband surveys in the 4deg2 area covered by Spitzer large area survey with Subaru/HSC (SPLASH) and 29 galaxies at z=5.148-7.508 with deep [CII] ALMA or PdBI data in the literature. Fluxes of strong rest-frame optical emission lines of [OIII] and Ha (Hb) are constrained by significant excesses found in the SPLASH 3.6 and 4.5um photometry. At z=4.9, the rest-frame Ha equivalent width, EW^0_Ha, and the Lya escape fraction, f_Lya, positively correlate with the rest-frame Lya equivalent width, EW_Lya. The f_Lya-EW_Lya correlation is similarly found at z~0-2, suggesting no evolution of the correlation over z~0-5. The typical ionizing photon production efficiency of LAEs is logxi_ion/[Hz erg^-1]~25.5, significantly (60-100%) higher than those of LBGs at a given UV magnitude. At z=5.7-7.0, there exists an interesting turn-over trend that the [OIII]/Ha flux ratio increases in EW_Lya~0-30A, where it reaches a maximum, and then decreases out to EW_Lya~130A. We also identify a negative correlation between a [CII] luminosity to star-formation rate ratio (L_[CII]/SFR) and EW_Lya at the 4.5sigma confidence level. We carefully investigate the physical origins of the trend and the correlations with stellar-synthesis and photoionization models covering the vast parameter space of metallicity, ionization parameter, and stellar age, and find that the [OIII]/Ha turn-over trend is explained by a simple anti-correlation between EW_Lya and metallicity, indicative of the detections of the very metal-poor (~0.04Zo) galaxies with EW_Lya~200A. Moreover, this simple anti-correlation explains self-consistently all of the correlations of Lya, Ha, [OIII]/Ha, and [CII] identified in our study.
Recently a theory about the formation of over-densities of stars along tidal tails of globular clusters has been presented, this theory predicts the position and time of formation of such over-densities and was successfully tested with N-body simulations of globular clusters in a point mass galactic potential. In this work we present a comparison between this theory and our simulations using a dwarf galaxy orbiting two differently shaped dark matter halos to study the effects of a cored and a cuspy halo on the formation and evolution of tidal tails. We find no difference using a cuspy or a cored halo, however, we find an intriguing asymmetry between the leading and trailing arm of the tidal tails. The trailing arm grows faster than the leading arm. This asymmetry is seen in the distance to first over-density and its size as well. We establish a relation between the distance to the first over-density and the size of this over-density.
We present an analysis of the average spectral properties of 12,000 SDSS quasars as a function of accretion disc inclination, as measured from the equivalent width of the [O III] 5007{\AA} line. The use of this indicator on a large sample of quasars from the SDSS DR7 has proven the presence of orientation effects on the features of UV/optical spectra, confirming the presence of outflows in the NLR gas and that the geometry of the BLR is disc-like. Relying on the goodness of this indicator, we are now using it to investigate other bands/components of AGN. Specifically, the study of the UV/optical/IR SED of the same sample provides information on the obscuring "torus". The SED shows a decrease of the IR fraction moving from face-on to edge-on sources, in agreement with models where the torus is co-axial with the accretion disc. Moreover, the fact we are able to observe the broad emission lines also in sources in an edge-on position, suggests that the torus is rather clumpy than smooth as in the Unified Model. The behaviour of the SED as a function of EW[OIII] is in agreement with the predictions of the clumpy torus models as well.
We study the spiral arm morphology of a sample of the local spiral galaxies in the Illustris simulation and explore the supermassive black hole - galaxy connection beyond the bulge (e.g., spiral arm pitch angle, total stellar mass, dark matter mass, and total halo mass), finding good agreement with other theoretical studies and observational constraints. It is important to study the properties of supermassive black holes and their host galaxies through both observations and simulations and compare their results in order to understand their physics and formative histories.We find that Illustris prediction for supermassive black hole mass relative to pitch angle is in rather good agreement with observations and that barred and non-barred galaxies follow similar scaling relations. Our work shows that Illustris presents very tight correlations between supermassive black hole mass and large-scale properties of the host galaxy, not only for early-type galaxies but also low-mass, blue and star-forming galaxies. These tight relations beyond the bulge suggest that halo properties determine those of a disc galaxy and its supermassive black hole.
A more fine-tuned method for probing planet-forming regions, such as protoplanetary discs, could be rovibrational molecular spectroscopy observation of particular premineral molecules instead of more common but ultimately less related volatile organic compounds. Planets are created when grains aggregate, but how molecules form grains is an ongoing topic of discussion in astrophysics and planetary science. Using the spectroscopic data of molecules specifically involved in mineral formation could help to map regions where planet formation is believed to be occurring in order to examine the interplay between gas and dust. Four atoms are frequently associated with planetary formation: Fe, Si, Mg, and O. Magnesium, in particular, has been shown to be in higher relative abundance in planet-hosting stars. Magnesium oxide crystals comprise the mineral periclase making it the chemically simplest magnesium-bearing mineral and a natural choice for analysis. The monomer, dimer, and trimer forms of (MgO)_n with n = 1 - 3 are analyzed in this work using high-level quantum chemical computations known to produce accurate results. Strong vibrational transitions at 12.5 {\mu}m, 15.0 {\mu}m, and 16.5 {\mu}m are indicative of magnesium oxide monomer, dimer, and trimer making these wavelengths of particular interest for the observation of protoplanetary discs and even potentially planet-forming regions around stars. If such transitions are observed in emission from the accretion discs or absorptions from stellar spectra, the beginning stages of mineral and, subsequently, rocky body formation could be indicated.
We present a novel strategy to uncover the Galactic population of quiescent black holes (BHs). This is based on a new concept, the photometric mass function (PMF), which opens up the possibility of an efficient identification of dynamical BHs in large fields-of-view. This exploits the width of the disc Halpha emission line, combined with orbital period information. We here show that Halpha widths can be recovered using a combination of customized Halpha filters. By setting a width cut-off at 2200 km/s we are able to cleanly remove other Galactic populations of Halpha emitters, including ~99.9% of cataclysmic variables (CVs). Only short period (Porb<2.1 h) eclipsing CVs and AGNs will contaminate the sample but these can be easily flagged through photometric variability and, in the latter case, also mid-IR colours. We also describe the strategy of a deep (r=22) Galactic plane survey based on the concept of PMFs: HAWKs, the HAlpha-Width Kilo-deg Survey. We estimate that ~800 sqr deg are required to unveil ~50 new dynamical BHs, a three-fold improvement over the known population. For comparison, a century would be needed to produce an enlarged sample of 50 dynamical BHs from X-ray transients at the current discovery rate.
The first four gravitational wave events detected by LIGO were all interpreted as merging black hole binaries (BHBs), opening a new perspective on the study of such systems. Here we use our new population-synthesis code MOBSE, an upgraded version of BSE (Hurley et al. 2002), to investigate the demography of merging BHBs. MOBSE includes metallicity-dependent prescriptions for mass loss of massive hot stars. It also accounts for the impact of the electron-scattering Eddington factor on mass loss. We perform >10^8 simulations of isolated massive binaries, with 12 different metallicities, to study the impact of mass loss, core-collapse supernovae and common envelope on merging BHBs. Accounting for the dependence of stellar winds on the Eddington factor leads to the formation of black holes (BHs) with mass up to 65 Msun at metallicity Z~0.0002. However, most BHs in merging BHBs have masses <40 Msun. We find merging BHBs with mass ratios in the 0.1 - 1.0 range, even if mass ratios >0.6 are more likely. We predict that systems like GW150914, GW170814 and GW170104 can form only from progenitors with metallicity Z<=0.006, Z<=0.008 and Z<=0.012, respectively. Most merging BHBs have gone through a common envelope phase, but up to ~17 per cent merging BHBs at low metallicity did not undergo any common envelope phase. We find a much higher number of mergers from metal-poor progenitors than from metal-rich ones: the number of BHB mergers per unit mass is ~10^-4 Msun^-1 at low metallicity (Z = 0.0002 - 0.002) and drops to ~10^-7 Msun^-1 at high metallicity (Z ~ 0.02).
Important insights into the formation and evolution of the Galactic disc(s) are contained in the chemical compositions of stars. We analysed high-resolution and high signal to noise HARPS spectra of 79 solar twin stars in order to obtain precise determinations of their atmospheric parameters, ages ($\sigma$$\sim$0.4 Gyr) and chemical abundances ($\sigma$$<$0.01~dex) of 12 neutron-capture elements (Sr, Y, Zr, Ba, La, Ce, Pr, Nd, Sm, Eu, Gd, and Dy). This valuable dataset allows us to study the [X/Fe]-age relations over a time interval of $\sim$10 Gyr and among stars belonging to the thin and thick discs. These relations show that i) the $s$-process has been the main channel of nucleosynthesis of $n$-capture elements during the evolution of the thin disc; ii) the thick disc is rich in $r$-process elements which suggests that its formation has been rapid and intensive. %; iii) a chemical continuity between the thin and thick discs is evident in the abundances of Ba. In addition, the heavy (Ba, La, Ce) and light (Sr, Y, Zr) $s$-process elements revealed details on the dependence between the yields of AGB stars and the stellar mass or metallicity. Finally, we confirmed that both [Y/Mg] and [Y/Al] ratios can be employed as stellar clocks, allowing ages of solar twin stars to be estimated with an average precision of $\sim$0.5~Gyr.
We use the first data release of the Gaia mission to explore the three dimensional arrangement and the age ordering of the many stellar groups towards the Orion OB association, aiming at a new classification and characterization of the stellar population. We make use of the parallaxes and proper motions provided in the Tycho Gaia Astrometric Solution (TGAS) sub-set of the Gaia catalogue, and of the combination of Gaia and 2MASS photometry. In TGAS we find evidence for the presence of a young population, at a parallax $\varpi \sim 2.65 \, \mathrm{mas}$, loosely distributed around some known clusters: 25 Ori, $\epsilon$ Ori and $\sigma$ Ori, and NGC 1980 ($\iota$ Ori). The low mass counterpart of this population is visible in the color-magnitude diagrams constructed by combining Gaia and 2MASS photometry. We study the density distribution of the young sources in the sky. We find the same groups as in TGAS, and also some other density enhancements that might be related to the recently discovered Orion X group, the Orion dust ring, and to the $\lambda$ Ori complex. We estimate the ages of this population and we infer the presence of an age gradient going from 25 Ori (13-15 Myr) to the ONC (1-2 Myr). We confirm this age ordering by repeating the Bayesian fit using the Pan-STARRS1 data. The estimated ages towards the NGC 1980 cluster span a broad range of values. This can either be due to the presence of two populations coming from two different episodes of star formation or to a large spread along the line of sight of the same population. Our results form the first step towards using the Gaia data to unravel the complex star formation history of the Orion region in terms of the different star formation episodes, their duration, and their effects on the surrounding interstellar medium.
In light of the observed Galactic center gamma-ray excess, we investigate a simplified model, for which the scalar dark matter interacts with quarks through a pseudoscalar mediator. The viable regions of the parameter space, that can also account for the relic density and evade the current searches, are identified, if the low-velocity dark matter annihilates through an $s$-channel off-shell mediator mostly into $\bar{b} b$, and/or annihilates directly into two ${\it hidden}$ on-shell mediators, which subsequently decay into the quark pairs. These two kinds of annihilations are $s$-wave. The projected monojet limit set by the high luminosity LHC sensitivity could constrain the favored parameter space, where the mediator's mass is larger than the dark matter mass by a factor of 2. We show that the projected sensitivity of 15-year Fermi-LAT observations of dwarf spheroidal galaxies can provide a stringent constraint on the most parameter space allowed in this model. If the on-shell mediator channel contributes to the dark matter annihilation cross sections over 50$\%$, this model with a lighter mediator can be probed in the projected PICO-500L experiment.
Links to: arXiv, form interface, find, astro-ph, recent, 1711, contact, help (Access key information)
We propose a novel method utilizing stellar kinematic data to detect low-mass substructure in the Milky Way's dark matter halo. By probing characteristic wakes that a passing dark matter subhalo leaves in the phase space distribution of ambient halo stars, we estimate sensitivities down to subhalo masses $\sim 10^7\,M_\odot$ or below. The detection of such subhalos would have implications for dark-matter and cosmological models that predict modifications to the halo-mass function at low halo masses. We develop an analytic formalism for describing the perturbed stellar phase-space distributions, and we demonstrate through simulations the ability to detect subhalos using the phase-space model and a likelihood framework. Our method complements existing methods for low-mass subhalo searches, such as searches for gaps in stellar streams, in that we can localize the positions and velocities of the subhalos today.
We use KMOS Deep Survey (KDS) galaxies, combined with results from a range of spectroscopic studies in the literature, to investigate the evolution of the stellar-mass Tully-Fisher relation since z ~ 4. We determine the slope and normalisation of the local rotation-velocity -- stellar-mass (Vc - $M_{\star}$) relationship using a reference sample of local spiral galaxies; thereafter we fix the slope, and focus on the evolution of velocity normalisation with redshift. The rotation-dominated KDS galaxies at z ~ 3.5 have rotation velocities ~ -0.1 dex lower than local reference galaxies at fixed stellar mass. By fitting 16 distant comparison samples spanning 0 < z < 3 (containing ~ 1200 galaxies), we show that the size and sign of the inferred Vc offset depends sensitively on the fraction of the parent samples used in the Tully-Fisher analysis, and how strictly the criterion of 'rotation dominated' is enforced. Confining attention to subsamples of galaxies that are especially 'disky' results in a consistent positive offset in Vc of ~ +0.1 dex, however these galaxies are not representative of the evolving-disk population at z > 1. We investigate the addition of pressure support, traced by intrinsic velocity dispersion ($\sigma_{int}$) to the KDS dynamical mass budget by adopting a 'total' effective velocity of form $V_{tot} = (Vc^{2} + 4.0\sigma_{int}^{2})^{0.5}$. The rotation-dominated and dispersion-dominated KDS galaxies fall on the same locus in the total-velocity versus stellar-mass plane, removing the need for debate over the precise selection threshold for rotation-dominated galaxies. The comparison sample offsets are in the range +0.08 to +0.15 dex in total-velocity zero-point (-0.30 to -0.55 dex in stellar-mass zero-point) from the local Tully-Fisher relation at z > 1, consistent with steady evolution of the ratio of dynamical to stellar mass with cosmic time.
Using data from MALT90 (Millimetre Astronomy Legacy Team Survey at 90 GHz), we present molecular line study of a sample of ATLASGAL (APEX Telescope Large Area Survey of the Galaxy) clumps. Twelve emission lines have been detected in all. We found that in most sources, emissions of HC$_3$N, HN$^{13}$C, CH$_3$CN, HNCO and SiO show more compact distributions than those of HCO$^+$, HNC, HCN and N$_2$H$^+$. By comparing with other molecular lines, we found that the abundance of HNCO ($\chi$(HNCO)) correlates well with other species such as HC3N, HNC, C2H, H13CO+ and N2H+. Previous studies indicate the HNCO abundance could be enhanced by shocks. However, in this study, we found the abundance of HNCO does not correlate well with that of SiO, which is also a good tracer of shocks. We suggest this may be because HNCO and SiO trace different parts of shocks. Our analysis indicates that the velocity of shock traced by HNCO tends to be lower than that traced by SiO. In the low-velocity shocks traced by HNCO, the HNCO abundance increases faster than that of SiO. While in the relatively high-velocity shocks traced by SiO, the SiO abundance increases faster than that of HNCO. We suggest that in the infrared dark cloud (IRDC) of MSXDC G331.71+00.59, high-velocity shocks are destroying the molecule of HNCO.
In this Letter, we present an analysis of the relation between the variability of broad absorption lines (BALs) and that of the continuum. Our sample is multi-epoch observations of 483 quasars by the Sloan Digital Sky Survey-I/II/III (SDSS-I/II/III). We derive the fractional flux variations of the continuum and fractional equivalent width (EW) variations for C IV and Si IV BALs, and explore the correlations between the three. Our results reveal moderate anticorrelations with high significance level between the fractional flux variations of the continuum and fractional EW variations for both C IV and Si IV BALs. We also prove a significant positive correlation between the fractional EW variations for C IV and Si IV BALs, which is in agreement with several previous studies. Our discoveries can serve as evidence for the idea: Change of an ionizing continuum is the primary driver of BAL variability.
We investigate emission lines of Lya, [OIII]5007, Ha, and [CII]158um from 1125 galaxies at z=4.9-7.0. Our sample is composed of 1098 Lya emitters (LAEs) at z=4.9, 5.7, 6.6, and 7.0 identified by Subaru/Hyper Suprime-Cam (HSC) narrowband surveys in the 4deg2 area covered by Spitzer large area survey with Subaru/HSC (SPLASH) and 29 galaxies at z=5.148-7.508 with deep [CII] ALMA or PdBI data in the literature. Fluxes of strong rest-frame optical emission lines of [OIII] and Ha (Hb) are constrained by significant excesses found in the SPLASH 3.6 and 4.5um photometry. At z=4.9, the rest-frame Ha equivalent width, EW^0_Ha, and the Lya escape fraction, f_Lya, positively correlate with the rest-frame Lya equivalent width, EW_Lya. The f_Lya-EW_Lya correlation is similarly found at z~0-2, suggesting no evolution of the correlation over z~0-5. The typical ionizing photon production efficiency of LAEs is logxi_ion/[Hz erg^-1]~25.5, significantly (60-100%) higher than those of LBGs at a given UV magnitude. At z=5.7-7.0, there exists an interesting turn-over trend that the [OIII]/Ha flux ratio increases in EW_Lya~0-30A, where it reaches a maximum, and then decreases out to EW_Lya~130A. We also identify a negative correlation between a [CII] luminosity to star-formation rate ratio (L_[CII]/SFR) and EW_Lya at the 4.5sigma confidence level. We carefully investigate the physical origins of the trend and the correlations with stellar-synthesis and photoionization models covering the vast parameter space of metallicity, ionization parameter, and stellar age, and find that the [OIII]/Ha turn-over trend is explained by a simple anti-correlation between EW_Lya and metallicity, indicative of the detections of the very metal-poor (~0.04Zo) galaxies with EW_Lya~200A. Moreover, this simple anti-correlation explains self-consistently all of the correlations of Lya, Ha, [OIII]/Ha, and [CII] identified in our study.
Recently a theory about the formation of over-densities of stars along tidal tails of globular clusters has been presented, this theory predicts the position and time of formation of such over-densities and was successfully tested with N-body simulations of globular clusters in a point mass galactic potential. In this work we present a comparison between this theory and our simulations using a dwarf galaxy orbiting two differently shaped dark matter halos to study the effects of a cored and a cuspy halo on the formation and evolution of tidal tails. We find no difference using a cuspy or a cored halo, however, we find an intriguing asymmetry between the leading and trailing arm of the tidal tails. The trailing arm grows faster than the leading arm. This asymmetry is seen in the distance to first over-density and its size as well. We establish a relation between the distance to the first over-density and the size of this over-density.
We present an analysis of the average spectral properties of 12,000 SDSS quasars as a function of accretion disc inclination, as measured from the equivalent width of the [O III] 5007{\AA} line. The use of this indicator on a large sample of quasars from the SDSS DR7 has proven the presence of orientation effects on the features of UV/optical spectra, confirming the presence of outflows in the NLR gas and that the geometry of the BLR is disc-like. Relying on the goodness of this indicator, we are now using it to investigate other bands/components of AGN. Specifically, the study of the UV/optical/IR SED of the same sample provides information on the obscuring "torus". The SED shows a decrease of the IR fraction moving from face-on to edge-on sources, in agreement with models where the torus is co-axial with the accretion disc. Moreover, the fact we are able to observe the broad emission lines also in sources in an edge-on position, suggests that the torus is rather clumpy than smooth as in the Unified Model. The behaviour of the SED as a function of EW[OIII] is in agreement with the predictions of the clumpy torus models as well.
We study the spiral arm morphology of a sample of the local spiral galaxies in the Illustris simulation and explore the supermassive black hole - galaxy connection beyond the bulge (e.g., spiral arm pitch angle, total stellar mass, dark matter mass, and total halo mass), finding good agreement with other theoretical studies and observational constraints. It is important to study the properties of supermassive black holes and their host galaxies through both observations and simulations and compare their results in order to understand their physics and formative histories.We find that Illustris prediction for supermassive black hole mass relative to pitch angle is in rather good agreement with observations and that barred and non-barred galaxies follow similar scaling relations. Our work shows that Illustris presents very tight correlations between supermassive black hole mass and large-scale properties of the host galaxy, not only for early-type galaxies but also low-mass, blue and star-forming galaxies. These tight relations beyond the bulge suggest that halo properties determine those of a disc galaxy and its supermassive black hole.
A more fine-tuned method for probing planet-forming regions, such as protoplanetary discs, could be rovibrational molecular spectroscopy observation of particular premineral molecules instead of more common but ultimately less related volatile organic compounds. Planets are created when grains aggregate, but how molecules form grains is an ongoing topic of discussion in astrophysics and planetary science. Using the spectroscopic data of molecules specifically involved in mineral formation could help to map regions where planet formation is believed to be occurring in order to examine the interplay between gas and dust. Four atoms are frequently associated with planetary formation: Fe, Si, Mg, and O. Magnesium, in particular, has been shown to be in higher relative abundance in planet-hosting stars. Magnesium oxide crystals comprise the mineral periclase making it the chemically simplest magnesium-bearing mineral and a natural choice for analysis. The monomer, dimer, and trimer forms of (MgO)_n with n = 1 - 3 are analyzed in this work using high-level quantum chemical computations known to produce accurate results. Strong vibrational transitions at 12.5 {\mu}m, 15.0 {\mu}m, and 16.5 {\mu}m are indicative of magnesium oxide monomer, dimer, and trimer making these wavelengths of particular interest for the observation of protoplanetary discs and even potentially planet-forming regions around stars. If such transitions are observed in emission from the accretion discs or absorptions from stellar spectra, the beginning stages of mineral and, subsequently, rocky body formation could be indicated.
We present a novel strategy to uncover the Galactic population of quiescent black holes (BHs). This is based on a new concept, the photometric mass function (PMF), which opens up the possibility of an efficient identification of dynamical BHs in large fields-of-view. This exploits the width of the disc Halpha emission line, combined with orbital period information. We here show that Halpha widths can be recovered using a combination of customized Halpha filters. By setting a width cut-off at 2200 km/s we are able to cleanly remove other Galactic populations of Halpha emitters, including ~99.9% of cataclysmic variables (CVs). Only short period (Porb<2.1 h) eclipsing CVs and AGNs will contaminate the sample but these can be easily flagged through photometric variability and, in the latter case, also mid-IR colours. We also describe the strategy of a deep (r=22) Galactic plane survey based on the concept of PMFs: HAWKs, the HAlpha-Width Kilo-deg Survey. We estimate that ~800 sqr deg are required to unveil ~50 new dynamical BHs, a three-fold improvement over the known population. For comparison, a century would be needed to produce an enlarged sample of 50 dynamical BHs from X-ray transients at the current discovery rate.
The first four gravitational wave events detected by LIGO were all interpreted as merging black hole binaries (BHBs), opening a new perspective on the study of such systems. Here we use our new population-synthesis code MOBSE, an upgraded version of BSE (Hurley et al. 2002), to investigate the demography of merging BHBs. MOBSE includes metallicity-dependent prescriptions for mass loss of massive hot stars. It also accounts for the impact of the electron-scattering Eddington factor on mass loss. We perform >10^8 simulations of isolated massive binaries, with 12 different metallicities, to study the impact of mass loss, core-collapse supernovae and common envelope on merging BHBs. Accounting for the dependence of stellar winds on the Eddington factor leads to the formation of black holes (BHs) with mass up to 65 Msun at metallicity Z~0.0002. However, most BHs in merging BHBs have masses <40 Msun. We find merging BHBs with mass ratios in the 0.1 - 1.0 range, even if mass ratios >0.6 are more likely. We predict that systems like GW150914, GW170814 and GW170104 can form only from progenitors with metallicity Z<=0.006, Z<=0.008 and Z<=0.012, respectively. Most merging BHBs have gone through a common envelope phase, but up to ~17 per cent merging BHBs at low metallicity did not undergo any common envelope phase. We find a much higher number of mergers from metal-poor progenitors than from metal-rich ones: the number of BHB mergers per unit mass is ~10^-4 Msun^-1 at low metallicity (Z = 0.0002 - 0.002) and drops to ~10^-7 Msun^-1 at high metallicity (Z ~ 0.02).
Important insights into the formation and evolution of the Galactic disc(s) are contained in the chemical compositions of stars. We analysed high-resolution and high signal to noise HARPS spectra of 79 solar twin stars in order to obtain precise determinations of their atmospheric parameters, ages ($\sigma$$\sim$0.4 Gyr) and chemical abundances ($\sigma$$<$0.01~dex) of 12 neutron-capture elements (Sr, Y, Zr, Ba, La, Ce, Pr, Nd, Sm, Eu, Gd, and Dy). This valuable dataset allows us to study the [X/Fe]-age relations over a time interval of $\sim$10 Gyr and among stars belonging to the thin and thick discs. These relations show that i) the $s$-process has been the main channel of nucleosynthesis of $n$-capture elements during the evolution of the thin disc; ii) the thick disc is rich in $r$-process elements which suggests that its formation has been rapid and intensive. %; iii) a chemical continuity between the thin and thick discs is evident in the abundances of Ba. In addition, the heavy (Ba, La, Ce) and light (Sr, Y, Zr) $s$-process elements revealed details on the dependence between the yields of AGB stars and the stellar mass or metallicity. Finally, we confirmed that both [Y/Mg] and [Y/Al] ratios can be employed as stellar clocks, allowing ages of solar twin stars to be estimated with an average precision of $\sim$0.5~Gyr.
We use the first data release of the Gaia mission to explore the three dimensional arrangement and the age ordering of the many stellar groups towards the Orion OB association, aiming at a new classification and characterization of the stellar population. We make use of the parallaxes and proper motions provided in the Tycho Gaia Astrometric Solution (TGAS) sub-set of the Gaia catalogue, and of the combination of Gaia and 2MASS photometry. In TGAS we find evidence for the presence of a young population, at a parallax $\varpi \sim 2.65 \, \mathrm{mas}$, loosely distributed around some known clusters: 25 Ori, $\epsilon$ Ori and $\sigma$ Ori, and NGC 1980 ($\iota$ Ori). The low mass counterpart of this population is visible in the color-magnitude diagrams constructed by combining Gaia and 2MASS photometry. We study the density distribution of the young sources in the sky. We find the same groups as in TGAS, and also some other density enhancements that might be related to the recently discovered Orion X group, the Orion dust ring, and to the $\lambda$ Ori complex. We estimate the ages of this population and we infer the presence of an age gradient going from 25 Ori (13-15 Myr) to the ONC (1-2 Myr). We confirm this age ordering by repeating the Bayesian fit using the Pan-STARRS1 data. The estimated ages towards the NGC 1980 cluster span a broad range of values. This can either be due to the presence of two populations coming from two different episodes of star formation or to a large spread along the line of sight of the same population. Our results form the first step towards using the Gaia data to unravel the complex star formation history of the Orion region in terms of the different star formation episodes, their duration, and their effects on the surrounding interstellar medium.
In light of the observed Galactic center gamma-ray excess, we investigate a simplified model, for which the scalar dark matter interacts with quarks through a pseudoscalar mediator. The viable regions of the parameter space, that can also account for the relic density and evade the current searches, are identified, if the low-velocity dark matter annihilates through an $s$-channel off-shell mediator mostly into $\bar{b} b$, and/or annihilates directly into two ${\it hidden}$ on-shell mediators, which subsequently decay into the quark pairs. These two kinds of annihilations are $s$-wave. The projected monojet limit set by the high luminosity LHC sensitivity could constrain the favored parameter space, where the mediator's mass is larger than the dark matter mass by a factor of 2. We show that the projected sensitivity of 15-year Fermi-LAT observations of dwarf spheroidal galaxies can provide a stringent constraint on the most parameter space allowed in this model. If the on-shell mediator channel contributes to the dark matter annihilation cross sections over 50$\%$, this model with a lighter mediator can be probed in the projected PICO-500L experiment.
Links to: arXiv, form interface, find, astro-ph, recent, 1711, contact, help (Access key information)
The coevolution of galaxies and their metal content serves as an important test for galaxy feedback models. We analyze the distribution and evolution of metals within the IllustrisTNG simulation suite with a focus on the gas-phase mass-metallicity relation (MZR). We find that the IllustrisTNG model broadly reproduces the slope and normalization evolution of the MZR across the redshift range $0<z<2$ and mass range $10^9 < M_*/\mathrm{M}_\odot < 10^{10.5}$. We make predictions for the high redshift ($2<z<10$) metal content of galaxies which is described by a gradual decline in the normalization of the metallicity with an average high redshift ($z>2$) evolution fit by $\mathrm{d\;log(Z)}/\mathrm{dz} \approx - 0.064$. Our simulations indicate that the metal retention efficiency of the interstellar medium (ISM) is low: a majority of gas-phase metals ($\sim$ 85 per cent at $z=0$) live outside of the ISM, either in an extended gas disk, the circumgalactic medium, or outside the halo. Nevertheless, the redshift evolution in the simulated MZR normalization is driven by the higher gas fractions of high redshift galaxies, not by changes to the metal retention efficiency. The scatter in the simulated MZR contains a clear correlation with the gas-mass or star formation rate of the system, in agreement with the observed fundamental metallicity relation. The scatter in the MZR is driven by a competition between periods of enrichment- and accretion-dominated metallicity evolution. We expect that while the normalization of the MZR declines with redshift, the strength of the correlation between metallicity and gas-mass at fixed stellar mass is not a strong function of redshift. Our results indicate that the "regulator" style models are best suited for simultaneously explaining the shape, redshift evolution, and existence of correlated scatter with gas fraction about the MZR.
Formation of supermassive stars (SMSs) with mass ~10^4 Msun is a promising pathway to seed the formation of supermassive black holes in the early universe. The so-called direct-collapse (DC) model postulates that such an SMS forms in a hot gas cloud irradiated by a nearby star-forming galaxy. We study the DC SMS formation in a fully cosmological context using three-dimensional radiation hydrodynamics simulations. We initialize our simulations using the outputs of the cosmological simulation of Chon et al. (2016), where two DC gas clouds are identified. The long-term evolution over a hundred thousand years is followed from the formation of embryo protostars through their growth to SMSs. We show that the strength of the tidal force by a nearby galaxy determines the multiplicity of the formed stars and affects the protostellar growth. In one case, where a collapsing cloud is significantly stretched by strong tidal force, multiple star-disk systems are formed via filament fragmentation. Small-scale fragmentation occurs in each circumstellar disk, and more than 10 stars with masses of a few times 10^3 Msun are finally formed. Interestingly, about a half of them are found as massive binary stars. In the other case, the gas cloud collapses nearly spherically under a relatively weak tidal field, and a single star-disk system is formed. Only a few SMSs with masses ~ 10^4 Msun are found already after evolution of a hundred thousand years, and the SMSs are expected to grow further by gas accretion and to leave massive blackholes at the end of their lives.
Galaxy clusters play a crucial role in constraining cosmological parameters. Cluster detection and characterization have therefore become an important field of modern cosmology. In this work, we aim at constructing reasonably representative composite luminosity function of cluster members in the WISE survey that can be used in the future as an input in matching algorithms for cluster detection. We use a sample of massive clusters from the redMaPPer catalog to match the positions of the WISE sources. We build the composite luminosity function of the WISE members in different redshift bins. We find that galaxy cluster members have a characteristic composite luminosity function, with a clear change in the slope at a given magnitude apparent $M^{*}$, which becomes fainter with increasing redshift. The best-fit bright-end slope $\beta$ is compatible with a constant value of $3.54$ with no clear trend with redshift, while the faint-end slope $\alpha$ is remarkably different, ranging from $2$ to $3$. We present the first characterization of the composite luminosity function of the WISE counterparts of redMaPPer cluster members and in this way we provide an element for building cluster detection techniques based on matching algorithms.
We present observational constraints on the stellar populations of two ultra-diffuse galaxies (UDGs) using optical through near-infrared (NIR) spectral energy distribution (SED) fitting. Our analysis is enabled by new $Spitzer$-IRAC 3.6 $\mu$m and 4.5 $\mu$m imaging, archival optical imaging, and the \texttt{prospector} fully Bayesian SED fitting framework. Our sample contains one field UDG (DGSAT I), one Virgo cluster UDG (VCC 1287), and one Virgo cluster dwarf elliptical for comparison (VCC 1122). Independently of SED fitting, we find that the optical--NIR colors of the three galaxies are significantly different from each other. We infer that VCC 1287 has an old ($\gtrsim7.7$ Gyr) and surprisingly metal-poor ($[Z/Z_{\odot}]\lesssim-1.0$) stellar population, even after marginalizing over uncertainties on diffuse interstellar dust. In contrast, the field UDG DGSAT I shows evidence of being younger than the Virgo UDG, with an extended star formation history and an age posterior extending down to $\sim3$ Gyr. The stellar metallicity of DGSAT I is sub-solar but higher than that of the Virgo UDG, with $[Z/Z_{\odot}]=-0.63^{+0.35}_{-0.62}$; in the case of exactly zero diffuse interstellar dust, DGSAT I may even have solar metallicity. The spectroscopically confirmed globular clusters of VCC 1287 have similar optical--NIR colors as the UDG itself, with empirical color relations suggesting sub-solar metallicities and supporting the metal-poor nature of VCC 1287. With VCC 1287 and several Coma UDGs, a general picture is emerging where cluster UDGs may be "failed" galaxies, but the field UDG DGSAT I seems more consistent with a stellar feedback-induced expansion scenario. [abridged]
Synthetic observations are playing an increasingly important role across astrophysics, both for interpreting real observations and also for making meaningful predictions from models. In this review, we provide an overview of methods and tools used for generating, manipulating and analysing synthetic observations and their application to problems involving star formation and the interstellar medium. We also discuss some possible directions for future research using synthetic observations.
Accurately predicting the demographics of dark matter (DM) substructure is of paramount importance for many fields of astrophysics, including gravitational lensing, galaxy evolution, halo occupation modeling, and constraining the nature of dark matter. Because of its strongly non-linear nature, DM substructure is typically modeled using N-body simulations, which reveal that large fractions of DM subhaloes undergo complete disruption. In this paper we use both analytical estimates and idealized numerical simulations to investigate whether this disruption is mainly physical, due to tidal heating and stripping, or numerical (i.e., artificial). We show that, contrary to naive expectation, subhaloes that experience a tidal shock $\Delta E$ that exceeds the subhalo's binding energy, $|E_{\rm b}|$, do not undergo disruption, even when $\Delta E/|E_{\rm b}|$ is as large as 100. Along the same line, and contrary to existing claims in the literature, instantaneously stripping matter from the outskirts of a DM subhalo also does not result in its complete disruption, even when the instantaneous remnant has positive binding energy. In addition, we show that tidal heating due to high-speed (impulsive) encounters with other subhaloes (`harassment'), is negligible compared to the tidal effects due to the host halo. Hence, we conclude that, in the absence of baryonic processes, the complete, physical disruption of CDM substructure is extremely rare, and that most disruption in numerical simulations therefore must be artificial. We discuss various processes that have been associated with numerical overmerging, and conclude that inadequate force-softening is the most likely culprit.
We present a multi-wavelength catalog in the Subaru-XMM Deep Field (SXDF) as part of the Spitzer Large Area Survey with Hyper-Suprime-Cam (SPLASH). We include the newly acquired optical data from the Hyper-Suprime Cam Subaru Strategic Program, accompanied by IRAC coverage from the SPLASH survey. All available optical and near-infrared data is homogenized and resampled on a common astrometric reference frame. Source detection is done using a multi-wavelength detection image including the $u$-band to recover the bluest objects. We measure multi-wavelength photometry and compute photometric redshifts as well as stellar masses for $\sim$800,000 objects over $\sim$4.3 deg$^2$. Using the available spectroscopic redshifts from various surveys over the range of $0<z<6$, we verify the performance of the photometric redshifts and we find a normalized median absolute deviation of 0.023 and outlier fraction of 3.4%. The SPLASH-SXDF catalog is a valuable, publicly available resource that is perfectly suited for studying galaxies in the early universe and tracing their evolution through cosmic time. The catalog is available for download from https://z.umn.edu/SXDF .
Energetic feedback from active galactic nuclei (AGNs) is difficult to observe, but often used in simulations to resolve several outstanding issues in galaxy formation. Here we derive new constraints on AGN feedback by comparing observations and simulations of the thermal Sunyaev-Zel'dovich (tSZ) effect, which causes spectral distortions in the cosmic microwave background with signal strengths that are proportional to the total energy deposited into the medium surrounding large galaxies. We draw on observational results presented in Spacek et al. (2016, 2017) who used data from the South Pole Telescope (SPT) and Atacama Cosmology Telescope (ACT) to measure the tSZ signal from >= 10^11 M_Sun and >= 1 Gyr galaxies in the redshift ranges of z=0.5-1.0 (low-z) and z=1.0-1.5 (high-z). Using two large-scale cosmological hydrodynamical simulations, one with AGN feedback (Horizon-AGN) and one without (Horizon-NoAGN), we extract simulated tSZ measurements around a population of galaxies equivalent to those observed. We find that the Horizon-AGN results only differ from the SPT measurements at levels of 0.4 sigma at low-z and 0.5 sigma at high-z, but they differ from the ACT measurements at levels of 6.9 sigma at low-z and 14.6 sigma at high-z. The Horizon-NoAGN results provide a slightly worse fit to the SPT measurements at levels of 1.8 sigma at low-z and 0.6 sigma at high-z, but a drastically better match to the ACT measurements at levels of 0.8 sigma at low-z and 1.9 sigma at high-z. We conclude that, while the lower-mass (<~ 5 x 10^11 M_Sun) SPT results show a hint of AGN feedback energy, the higher-mass (>~ 5 x 10^11 M_Sun) ACT results show significantly less energy than predicted in the simulation including feedback, indicating that AGN feedback may be milder than often predicted.
Feedback in the form of mass outflows driven by star formation or active galactic nuclei is a key component of galaxy evolution. The luminous infrared galaxy Zw 049.057 harbours a compact obscured nucleus with a possible far-IR signature of outflowing molecular gas. Due to the high optical depths at far-IR wavelengths, the interpretation of the outflow signature is uncertain. At mm and radio wavelengths, the radiation is better able to penetrate the large columns of gas and dust. We used high resolution observations from the SMA, ALMA, and the VLA to image the CO 2-1 and 6-5 emission, the 690 GHz continuum, the radio cm continuum, and absorptions by rotationally excited OH. The CO line profiles exhibit wings extending 300 km/s beyond the systemic velocity. At cm wavelengths, we find a compact (40 pc) continuum component in the nucleus, with weaker emission extending several 100 pc approximately along the major and minor axes of the galaxy. In the OH absorption lines toward the compact continuum, wings extending to a similar velocity as for the CO are seen on the blue side of the profile. The weak cm continuum emission along the minor axis is aligned with a highly collimated, jet-like dust feature previously seen in near-IR images of the galaxy. Comparison of the apparent optical depths in the OH lines indicate that the excitation conditions in Zw 049.057 differ from those in other OH megamaser galaxies. We interpret the wings in the spectral lines as signatures of a molecular outflow. A relation between this outflow and the minor axis radio feature is possible, although further studies are required to investigate this possible association and understand the connection between the outflow and the nuclear activity. Finally, we suggest that the differing OH excitation conditions are further evidence that Zw 049.057 is in a transition phase between megamaser and kilomaser activity.
This overview provides a historical perspective highlighting the pioneering
role which the fairly modest observational facilities of ARIES have played
since the 1990s in systematically characterizing the optical variability on
hour-like time scale (intra-night optical variability, or INOV) of several
major types of high-luminosity Active Galactic Nuclei (AGN). Such information
was previously available only for blazars. Similar studies have since been
initiated in at least a dozen countries, giving a boost to AGN variability
research.
Our work has, in particular, provided strong indication that mild INOV occurs
in radio-quiet QSOs (amplitude up to $\sim 3-5\%$ and duty cycle $\sim 10\%$)
and, moreover, has demonstrated that similarly mild INOV is exhibited even by
the vast majority of radio-loud quasars which possess powerful relativistic
jets (even including many that are beamed towards us). The solitary outliers
are blazars, the tiny strongly polarized subset of powerful AGN, which
frequently exhibit a pronounced INOV. Among the blazars, BL Lac objects often
show a bluer-when-brighter chromatic behavior, while the flat spectrum radio
quasars seem not to. Quantifying any differences of INOV among the major
subclasses of non-blazar type AGNs will require dedicated monitoring programs
using $2-3$ metre class telescopes.
We present high sensitivity ($\sigma_P \simeq 0.6\,$mJy) polarimetric observations in seven bands, from $2.1$ to $38\,$GHz, of a complete sample of $104$ compact extragalactic radio sources brighter than $200\,$mJy at $20\,$GHz. Polarization measurements in six bands, in the range $5.5-38\,$GHz, for $53$ of these objects were reported by \citet{Galluzzi2017}. We have added new measurements in the same six bands for another 51 sources and measurements at $2.1\,$GHz for the full sample of $104$ sources. Also, the previous measurements at $18$, $24$, $33$ and $38\,$GHz were re-calibrated using the updated model for the flux density absolute calibrator, PKS1934-638, not available for the earlier analysis. The observations, carried out with the Australia Telescope Compact Array (ATCA), achieved a $90\%$ detection rate (at $5\sigma$) in polarization. $89$ of our sources have a counterpart in the $72$ to $231\,$MHz GLEAM survey \citep{HurleyWalker2017}, providing an unparalleled spectral coverage of $2.7$ decades of frequency for these sources. While the total intensity data from $5.5$ to $38\,$GHz could be interpreted in terms of single component emission, a joint analysis of more extended total intensity spectra presented here, and of the polarization spectra, reveals that over $90\%$ of our sources show clear indications of at least two emission components. We interpret this as an evidence of recurrent activity. Our high sensitivity polarimetry has allowed a $5\,\sigma$ detection of the weak circular polarization for $\sim 38\%$ of the dataset, and a deeper estimate of $20\,$GHz polarization source counts than has been possible so far.
We present optical spectroscopy secured at the 10m Gran Telescopio Canarias of the counterparts of 20 extragalactic gamma-ray sources detected by the Fermi satellite. The observations allow us to investigate the nature of these sources and to determine their redshift. We find that all optical counterparts have a spectrum that is consistent with a BL Lac object nature. We are able to determine the redshift for 11 objects and set spectroscopic redshift limits for five targets. Only for four sources the optical spectrum is found featureless. In the latter cases we can set lower limits on the redshift based on the assumption that they are hosted by a typical massive elliptical galaxy whose spectrum is diluted by the non thermal continuum. The observations allow us to unveil the nature of these gamma-ray sources and provide a sanity check of a tool to discover the counterparts of gamma-ray emitters/blazars based on their multiwavelength emission.
Many stars, active galactic nuclei, accretion discs etc. are affected by the stochastic variations of temperature, turbulent gas motions, magnetic fields, number densities of atoms and dust grains. These stochastic variations influence on the extinction factors, Doppler widths of lines and so on. The presence of many reasons for fluctuations gives rise to Gaussian distribution of fluctuations. The usual models leave out of account the fluctuations. In many cases the consideration of fluctuations improves the coincidence of theoretical values with the observed data. The objective of this paper is the investigation of the influence of the number density fluctuations on the form of radiative transfer equations. We consider non-magnetized atmosphere in continuum.
First proposed by Paczynski in 1986, microlensing has been instrumental in the search of compact dark matter as well as discovery and characterization of exoplanets. In this article, we provide a brief history of microlensing, especially on the discoveries of compact objects and exoplanets. We then review the basics of microlensing and how astrometry can help break the degeneracy, providing a more robust determination of the nature of the microlensing events. We also outline prospects that will be made by on-going and forth-coming experiments/observatories
Links to: arXiv, form interface, find, astro-ph, recent, 1711, contact, help (Access key information)
[abridged] The enhanced degrees of deuterium fractionation observed in envelopes around protostars demonstrate the importance of chemistry at low temperatures, relevant in pre- and protostellar cores. formaldehyde is an important species in the formation of methanol and more complex molecules. Here, we present the first study of formaldehyde deuteration on small scales around the prototypical low-mass protostar IRAS 16293-2422 using high spatial and spectral resolution Atacama Large Millimeter/submillimeter Array (ALMA) observations. Numerous isotopologues of formaldehyde are detected, among them H$_2$C$^{17}$O, and D$_2^{13}$CO for the first time in the ISM. The large range of upper energy levels covered by the HDCO lines help constrain the excitation temperature to 106$\pm$13 K. Using the derived column densities, formaldehyde shows a deuterium fractionation of HDCO/H$_2$CO=6.5$\pm$1%, D$_2$CO/HDCO=12.8$^{+3.3}_{-4.1}$%, and D$_2$CO/H$_2$CO=0.6(4)$\pm$0.1%. The isotopic ratios derived are $^{16}$O/$^{18}$O=805, $^{18}$O/$^{17}$O=3.2 and $^{12}$C/$^{13}$C=56. The HDCO/H$_2$CO ratio is lower than found in previous studies, highlighting the uncertainties involved in interpreting single dish observations of the inner warm regions. The D$_2$CO/HDCO ratio is only slightly larger than the HDCO/H$_2$CO ratio. This is consistent with formaldehyde forming in the ice as soon as CO has frozen onto the grains, with most of the deuteration happening towards the end of the prestellar core phase. A comparison with available time-dependent chemical models indicates that the source is in the early Class 0 stage.
Pair-instability supernovae (PISNe) are very luminous explosions of massive, low metallicity stars. They can potentially be observed out to high redshifts due to their high explosion energies, thus providing a probe of the Universe prior to reionisation. The near-infrared camera, NIRCam, on board the James Webb Space Telescope is ideally suited for detecting their redshifted ultraviolet emission. We calculate the photometric signature of high-redshift PISNe and derive the optimal survey strategy for identifying their prompt emission and possible afterglow. We differentiate between PISNe and other sources that could have a similar photometric signature, such as active galactic nuclei or high redshift galaxies. We demonstrate that the optimal survey strategy, which maximises the visibility time of the PISN lightcurve per invested exposure time, consists of the two wide-band filters F200W and F356W with an exposure time of 600s. The PISN afterglow, caused by nebular emission and reverberation, is very faint and requires unfeasibly long exposure times to be uniquely identified. However, this afterglow would be visible for several hundred years, about two orders of magnitude longer than the prompt emission, rendering PISNe promising targets for future, even more powerful telescopes.
We present a morphological catalogue for $\sim$ 670,000 galaxies in the Sloan
Digital Sky Survey in two flavours: T-Type, related to the Hubble sequence, and
Galaxy Zoo 2 (GZ2 hereafter) classification scheme. By combining accurate
existing visual classification catalogues with machine learning, we provide the
largest and most accurate morphological catalogue up to date. The
classifications are obtained with Deep Learning algorithms using Convolutional
Neural Networks (CNNs).
We use two visual classification catalogues, the GZ2 and the Nair et al.
2010, for training CNNs with color images in order to obtain T-Types and a
series of GZ2 type questions (disk/features, edge-on galaxies, bar signature,
bulge prominence, roundness and mergers). We also provide an additional
probability enabling a separation between pure elliptical (E) from S0, where
the T-Type model is not so efficient. For the T-Type, our results show smaller
offset and scatter than previous models trained with support vector machines.
For the GZ2 type questions, our models have large accuracy (> 97\%), precision
and recall values (> 90\%) when applied to a test sample with the same
characteristics as the one used for training, i.e., with small uncertainties in
the GZ2 classification. When we apply our models to a subset of galaxies with
larger uncertainties in GZ2 (low agreement between classifiers), the resulting
probability distributions show a clear bimodality. This allows us to recover a
significant fraction of galaxies with a robust classification, for which the
GZ2 was classification was uncertain. The catalogue is publicly released with
the paper.
We attempt to determine the relative fraction of stars that have undergone significant radial migration by studying the orbital properties of metal-rich ([Fe/H]$>0.1$) stars within 2 kpc of the Sun using a sample of more than 3,000 stars selected from iDR4 of the Gaia-ESO Survey. We investigate the kinematic properties, such as velocity dispersion and orbital parameters, of stellar populations near the sun as a function of [Mg/Fe] and [Fe/H], which could show evidence of a major merger in the past history of the Milky Way. This was done using the stellar parameters from the Gaia-ESO Survey along with proper motions from PPMXL to determine distances, kinematics, and orbital properties for these stars to analyze the chemodynamic properties of stellar populations near the Sun. Analyzing the kinematics of the most metal-rich stars ([Fe/H]$>0.1$), we find that more than half have small eccentricities ($e<0.2$) or are on nearly circular orbits. Slightly more than 20\% of the metal-rich stars have perigalacticons $R_p>7$ kpc. We find that the highest [Mg/Fe], metal-poor populations have lower vertical and radial velocity dispersions compared to lower [Mg/Fe] populations of similar metallicity by $\sim10$ km s$^{-1}$. The median eccentricity increases linearly with [Mg/Fe] across all metallicities, while the perigalacticon decreases with increasing [Mg/Fe] for all metallicities. Finally, the most [Mg/Fe]-rich stars are found to have significant asymmetric drift and rotate more than 40 km s$^{-1}$ slower than stars with lower [Mg/Fe] ratios. While our results cannot constrain how far stars have migrated, we propose that migration processes are likely to have played an important role in the evolution of the Milky Way, with metal-rich stars migrating from the inner disk toward to solar neighborhood and past mergers potentially driving enhanced migration of older stellar populations in the disk.
We find two chemically distinct populations separated relatively cleanly in the [Fe/H] - [Mg/Fe] plane, but also distinguished in other chemical planes, among metal-poor stars (primarily with metallicities [Fe/H] $< -0.9$) observed by the Apache Point Observatory Galactic Evolution Experiment (APOGEE) and analyzed for Data Release 13 (DR13) of the Sloan Digital Sky Survey. These two stellar populations show the most significant differences in their [X/Fe] ratios for the $\alpha$-elements, C+N, Al, and Ni. In addition to these populations having differing chemistry, the low metallicity high-Mg population (which we denote the HMg population) exhibits a significant net Galactic rotation, whereas the low-Mg population (or LMg population) has halo-like kinematics with little to no net rotation. Based on its properties, the origin of the LMg population is likely as an accreted population of stars. The HMg population shows chemistry (and to an extent kinematics) similar to the thick disk, and is likely associated with $\it in$ $\it situ$ formation. The distinction between the LMg and HMg populations mimics the differences between the populations of low- and high-$\alpha$ halo stars found in previous studies, suggesting that these are samples of the same two populations.
We present the first constraints on stellar mass-to-light ratio gradients in an early-type galaxy (ETG) using multiple dynamical tracer populations to model the dark and luminous mass structure simultaneously. We combine the kinematics of the central starlight, two globular cluster populations and satellite galaxies in a Jeans analysis to obtain new constraints on M87's mass structure, employing a flexible mass model which allows for radial gradients in the stellar mass-to-light ratio. We find that, in the context of our model, a radially declining stellar-mass-to-light ratio is strongly favoured. Modelling the stellar mass-to-light ratio as following a power law, $\Upsilon_{\star} \sim R^{-\mu}$, we infer a power-law slope $\mu = -0.54 \pm 0.05$; equally, parameterising the stellar-mass-to-light ratio via a central mismatch parameter relative to a Salpeter IMF, $\alpha$, and scale radius $R_M$, we find $\alpha > 1.48$ at $95\%$ confidence and $R_M = 0.35 \pm 0.04$ kpc. We use stellar population modelling of high-resolution 11-band HST photometry to show that such a steep gradient cannot be achieved by variations in only the metallicity, age, dust extinction and star formation history if the stellar initial mass function (IMF) remains spatially constant. On the other hand, the stellar mass-to-light ratio gradient that we find is consistent with an IMF whose inner slope changes such that it is Salpeter-like in the central $\sim 0.5$ kpc and becomes Chabrier-like within the stellar effective radius. This adds to recent evidence that the non-universality of the IMF in ETGs may be confined to their core regions, and points towards a picture in which the stars in these central regions may have formed in fundamentally different physical conditions.
We present preliminary results of the high resolution $ (0.10^{\prime\prime} \times 0.15^{\prime\prime}) $ observations of the high mass star forming region S255IR with ALMA in several spectral windows from ~ 335 GHz to ~ 350 GHz. The main target lines were C$^{34}$S(7-6), CH$_3$CN($19_K-18_K$), CO(3-2) and SiO(8-7), however many other lines of various molecules have been detected, too. We present sample spectra and maps, discuss briefly the source structure and kinematics. A new, never predicted methanol maser line has been discovered.
We present a study of the inexplicit connection between radio jet activity and gamma-ray emission of BL Lacertae (BL Lac; 2200+420). We analyze the long-term millimeter activity of BL Lac via interferometric observations with the Korean VLBI Network (KVN) obtained at 22, 43, 86, and 129 GHz simultaneously over three years (from January 2013 to March 2016); during this time, two gamma-ray outbursts (in November 2013 and March 2015) can be seen in gamma-ray light curves obtained from Fermi observations. The KVN radio core is optically thick at least up to 86 GHz; there is indication that it might be optically thin at higher frequencies. To first order, the radio light curves decay exponentially over the time span covered by our observations, with decay timescales of 411+/-85 days, 352+/-79 days, 310+/-57 days, and 283+/-55 days at 22, 43, 86, and 129 GHz, respectively. Assuming synchrotron cooling, a cooling time of around one year is consistent with magnetic field strengths B~2microT and electron Lorentz factors gamma~10,000. Taking into account that our formal measurement errors include intrinsic variability and thus over-estimate the statistical uncertainties, we find that the decay timescale tau scales with frequency nu like tau~nu^{-0.2}. This relation is much shallower than the one expected from opacity effects (core shift), but in agreement with the (sub-)mm radio core being a standing recollimation shock. We do not find convincing radio flux counterparts to the gamma$ray outbursts. The spectral evolution is consistent with the `generalized shock model' of Valtaoja et al. (1992). A temporary increase in the core opacity and the emergence of a knot around the time of the second gamma-ray event indicate that this gamma-ray outburst might be an `orphan' flare powered by the `ring of fire' mechanism.
Out of the $\sim$ 3000 exoplanets detected so far, only fourteen planets are members of open clusters: among them an exoplanet system around Pr 0211 in the cluster M44 which consists of at least two planets with the outer planet moving on a highly eccentric orbit at 5.5 AU. One hypothesis is that a close fly-by of a neighbouring star was responsible for the eccentric orbit. We test this hypothesis. First we determine the type of fly-by that would lead to the observed parameters and then use this result to determine the history of such fly-bys in simulations of the early dynamics in an M44-like environment. We find that although very close fly-bys are required to obtain the observed properties of Pr 0211c, such fly-bys are relatively common due to the high stellar density and longevity of the cluster. Such close fly-bys are most frequent during the first 1-2 Myr after cluster formation, corresponding to a cluster age $\leq$ 3 Myr. During the first 2 to 3 Myr about 6.5% of stars actually experience a fly-by that would lead to such a small system-size as observed for Pr0211 or even smaller. It is unclear whether planets generally form on such short timescales. However, afterwards the close fly-by rate is still 0.2-0.5 Myr$^{-1}$, which means extrapolating this to the age of M44 12%-20% of stars would experience such close fly-bys over this timespan. Our simulations show that the fly-by scenario is a realistic option for the formation of eccentricity orbits of the planets in M44. The occurrence of such events is relatively high leading to the expectation that similar systems are likely common in open clusters in general.
We present properties of the low-surface-brightness galaxy KDG 218 observed with the HST/ACS. The galaxy has a half-light (effective) diameter of $a_e =47^{\prime\prime}$ and a central surface brightness of $SB_V(0) = 24.4^m/sq^{\prime\prime}$. The galaxy remains unresolved with the HST/ACS, which implies its distance of {$D > 13.1$ Mpc} and linear effective diameter of $A_e >3.0$ kpc. We notice that KDG 218 is most likely associated with a galaxy group around the massive lenticular NGC 4958 galaxy at approximately $22$ Mpc, or with the Virgo Southern Extension filament at approximately $16.5$ Mpc. At these distances, the galaxy is classified as an ultra-diffuse galaxy (UDG) similar to those found in the Virgo, Fornax, and Coma clusters. We also present a sample of 15 UDG candidates in the Local Volume. These sample galaxies have the following mean parameters: $\langle D\rangle = 5.1$ Mpc, $\langle A_e\rangle = 4.8$ kpc, and $\langle SB_B (e) \rangle = 27.4^m/sq^{\prime\prime}$. All the local UDG candidates reside near massive galaxies located in the regions with the mean stellar mass density (within 1 Mpc) about 50 times greater than the average cosmic density. The local fraction of UDGs does not exceed 1.5% of the Local Volume population. We notice that the presented sample of local UDGs is a heterogeneous one containing irregular, transition, and tidal types, as well as objects consisting of an old stellar population.
On the basis of the SDSS survey and spectral observations with the 6-m telescope of SAO RAS, we have peformed a detailed study of SDSS J170745+302056. By combination of its characteristics -- exponential surface brightness distribution, central surface brightness of stellar disk $\mu_0(B) = 23.25\,^m/\square"$, blue colors, low metallicity and low star formation rate -- the galaxy is a typical low surface brightness spiral galaxy. Exponential scalelength of the galaxy is $\approx$3 kpc, its optical diameter exceeds 20 kpc. SDSS J170745+302056 is a member of a group of five galaxies and probably it is in interaction with UGC 10716. The existence of a large low surface brightness galaxy in such a dense environment is very unusual.
We apply four statistical learning methods to a sample of $7941$ galaxies ($z<0.06$) from the Galaxy and Mass Assembly (GAMA) survey to test the feasibility of using automated algorithms to classify galaxies. Using $10$ features measured for each galaxy (sizes, colours, shape parameters \& stellar mass) we apply the techniques of Support Vector Machines (SVM), Classification Trees (CT), Classification Trees with Random Forest (CTRF) and Neural Networks (NN), returning True Prediction Ratios (TPRs) of $75.8\%$, $69.0\%$, $76.2\%$ and $76.0\%$ respectively. Those occasions whereby all four algorithms agree with each other yet disagree with the visual classification (`unanimous disagreement') serves as a potential indicator of human error in classification, occurring in $\sim9\%$ of ellipticals, $\sim9\%$ of Little Blue Spheroids, $\sim14\%$ of early-type spirals, $\sim21\%$ of intermediate-type spirals and $\sim4\%$ of late-type spirals \& irregulars. We observe that the choice of parameters rather than that of algorithms is more crucial in determining classification accuracy. Due to its simplicity in formulation and implementation, we recommend the CTRF algorithm for classifying future galaxy datasets. Adopting the CTRF algorithm, the TPRs of the 5 galaxy types are : E, $70.1\%$; LBS, $75.6\%$; S0-Sa, $63.6\%$; Sab-Scd, $56.4\%$ and Sd-Irr, $88.9\%$. Further, we train a binary classifier using this CTRF algorithm that divides galaxies into spheroid-dominated (E, LBS \& S0-Sa) and disk-dominated (Sab-Scd \& Sd-Irr), achieving an overall accuracy of $89.8\%$. This translates into an accuracy of $84.9\%$ for spheroid-dominated systems and $92.5\%$ for disk-dominated systems.
Fossil galaxy systems are classically thought to be the end result of galaxy group/cluster evolution, as galaxies experiencing dynamical friction sink to the center of the group potential and merge into a single, giant elliptical that dominates the rest of the members in both mass and luminosity. Most fossil systems discovered lie within $z < 0.2$, which leads to the question: what were these systems' progenitors? Such progenitors are expected to have imminent or ongoing major merging near the brightest group galaxy (BGG) that, when concluded, will meet the fossil criteria within the look back time. Since strong gravitational lensing preferentially selects groups merging along the line of sight, or systems with a high mass concentration like fossil systems, we searched the CASSOWARY survey of strong lensing events with the goal of determining if lensing systems have any predisposition to being fossil systems or progenitors. We find that $\sim$13% of lensing groups are identified as traditional fossils while only $\sim$3% of non-lensing control groups are. We also find that $\sim$23% of lensing systems are traditional fossil progenitors compared to $\sim$17% for the control sample. Our findings show that strong lensing systems are more likely to be fossil/pre-fossil systems than comparable non-lensing systems. Cumulative galaxy luminosity functions of the lensing and non-lensing groups also indicate a possible, fundamental difference between strong lensing and non-lensing systems' galaxy populations with lensing systems housing a greater number of bright galaxies even in the outskirts of groups.
The formation processes that led to the current Galactic stellar halo are still under debate. Previous studies have provided evidence for different stellar populations in terms of elemental abundances and kinematics, pointing to different chemical and star-formation histories. In the present work we explore, over a broader range in metallicity (-2.2 < [Fe/H] < -0.5), the two stellar populations detected in the first paper of this series from metal-poor stars in DR13 of the Apache Point Observatory Galactic Evolution Experiment (APOGEE). We aim to infer signatures of the initial mass function (IMF) and the most APOGEE-reliable alpha-elements (O, Mg, Si and Ca). Using simple chemical-evolution models, for each population. Compared with the low-alpha population, we obtain a more intense and longer-lived SFH, and a top-heavier IMF for the high-alpha population.
The Outer Scutum-Centaurus arm (OSC) is the most distant molecular spiral arm known in the Milky Way. The OSC may be the very distant end of the well-known Scutum-Centaurus arm, which stretches from the end of the Galactic bar to the outer Galaxy. At this distance the OSC is seen in the first Galactic quadrant. The population of star formation tracers in the OSC remains largely uncharacterized. Extragalactic studies show a strong correlation between molecular gas and star formation, and carbon monoxide (CO) emission was recently discovered in the OSC. Here we use the Arizona Radio Observatory (ARO) 12-m telescope to observe the $^{12}$CO J = 1-0 and $^{13}$CO J = 1-0 transitions toward 78 HII region candidates chosen from the WISE Catalog of Galactic HII Regions. These targets are spatially coincident with the Galactic longitude-latitude ($\ell, b$) OSC locus as defined by HI emission. We detect CO emission in $\sim 80$% of our targets. In total, we detect 117 $^{12}$CO and 40 $^{13}$CO emission lines. About 2/3 of our targets have at least one emission line originating beyond the Solar orbit. Most of the detections beyond the Solar orbit are associated with the Outer Arm, but there are 17 $^{12}$CO emission lines and 8 $^{13}$CO emission lines with LSR velocities that are consistent with the velocities of the OSC. There is no apparent difference between the physical properties (e.g., molecular column density) of these OSC molecular clouds and non--OSC molecular clouds within our sample.
Recently, a new form of dark matter has been suggested to naturally reproduce the empirically successful aspects of Milgrom's law in galaxies. The dark matter particle candidates are axion-like, with masses of order eV and strong self-interactions. They Bose-Einstein condense into a superfluid phase in the central regions of galaxy halos. The superfluid phonon excitations in turn couple to baryons and mediate an additional long-range force. For a suitable choice of the superfluid equation of state, this force can mimic Milgrom's law. In this paper we develop in detail some of the main phenomenological consequences of such a formalism, by revisiting the expected dark matter halo profile in the presence of an extended baryon distribution. In particular, we show how rotation curves of both high and low surface brightness galaxies can be reproduced, with a slightly rising rotation curve at large radii in massive high surface brightness galaxies, thus subtly different from Milgrom's law. We finally point out other expected differences with Milgrom's law, in particular in dwarf spheroidal satellite galaxies, tidal dwarf galaxies, and globular clusters, whose Milgromian or Newtonian behavior depends on the position with respect to the superfluid core of the host galaxy. We also expect ultra-diffuse galaxies within galaxy clusters to have velocities slightly above the baryonic Tully-Fisher relation. Finally, we note that, in this framework, photons and gravitons follow the same geodesics, and that galaxy-galaxy lensing, probing larger distances within galaxy halos than rotation curves, should follow predictions closer to the standard cosmological model than those of Milgrom's law.
We present new Chandra observations of the outer halo of the giant elliptical galaxy NGC 4472 (M49) in the Virgo Cluster. The data extend to 130 kpc (28'), and have a combined exposure time of 150 ks. After identifying optical counterparts using the Next Generation Virgo Cluster Survey to remove background active galactic nuclei and globular cluster (GC) sources, and correcting for completeness, we find that the number of field low-mass X-ray binaries (LMXBs) per unit stellar V-band light increases significantly with galactocentric radius. Because the flux limit of the complete sample corresponds to the Eddington limit for neutron stars in NGC 4472, many of the ~90 field LMXBs in this sample could host black holes. The excess of field LMXBs at large galactocentric radii may be partially caused by natal kicks on black holes and neutron stars in binary systems in the inner part of the galaxy. Furthermore, since the metallicity in the halo of NGC 4472 strongly decreases towards larger galactocentric radii, the number of field LMXBs per unit stellar mass is anti-correlated with metallicity, opposite to what is observed in GCs. Another way to explain the spatial distribution of field LMXBs is therefore a reversed metallicity effect, although we have not identified a mechanism to explain this in terms of stellar and binary evolution.
The University of Michigan 26-m paraboloid was dedicated to obtaining linear polarization and total flux density observations of blazars from the mid-1960s until June 2012 providing an unprecedented record tracking centimeter-band variability over decades at 14.5, 8.0, and 4.8 GHz for both targeted objects and members of flux-limited samples. In the mid-1970s through the mid-1980s, and during the last decade of the program, observations were additionally obtained of circular polarization for a small sample of radio-bright (S>5Jy), active sources. Key program results include evidence supporting class-dependent differences in the magnetic field geometry of BL Lac and QSO jets, identification of linear polarization changes temporally associated with flux outbursts supporting a shock-in-jet scenario, and determination of the spectral evolution of the Stokes V amplitude and polarity for testing proposed models. Recent radiative transfer modeling during large flares supports a jet scenario with a kinetically-dominated, relativistic flow at parsec scales with embedded turbulent magnetic fields and dynamically-weak ordered components which may be helical; the circular polarization observations are consistent with linear-to-circular mode conversion within this turbulent jet environment.
Observationally, supernovae (SNe) are divided into subclasses pertaining to their distinct characteristics. This diversity reflects the diversity in the progenitor stars. It is not entirely clear how different evolutionary paths leading massive stars to become a SN are governed by fundamental parameters such as progenitor initial mass and metallicity. This paper places constraints on progenitor initial mass and metallicity in distinct core-collapse SN subclasses, through a study of the parent stellar populations at the explosion sites. Integral field spectroscopy (IFS) of 83 nearby SN explosion sites with a median distance of 18 Mpc has been collected and analysed, enabling detection and spectral extraction of the parent stellar population of SN progenitors. From the parent stellar population spectrum, the initial mass and metallicity of the coeval progenitor are derived by means of comparison to simple stellar population models and strong-line methods. Additionally, near-infrared IFS was employed to characterise the star formation history at the explosion sites. No significant metallicity differences are observed among distinct SN types. The typical progenitor mass is found to be highest for SN Ic, followed by type Ib, then types IIb and II. SN IIn is the least associated with young stellar populations and thus massive progenitors. However, statistically significant differences in progenitor initial mass are observed only when comparing SNe IIn with other subclasses. Stripped-envelope SN progenitors with initial mass estimate lower than 25~$M_\odot$ are found; these are thought to be the result of binary progenitors. Confirming previous studies, these results support the notion that core-collapse SN progenitors cannot arise from single-star channel only, and both single and binary channels are at play in the production of core-collapse SNe. [ABRIDGED]
The driving mechanism of massive outflows observed in high-mass star-forming regions is investigated using three-dimensional magnetohydrodynamics (MHD) and protostellar evolution calculations. In our previous paper, we showed that the mass outflow rate depends strongly on the mass accretion rate onto the circumstellar disk around a high-mass protostar, and massive outflows may be driven by the magnetic effect in high-mass star-forming cores. In the present study, in order to verify that the MHD disk wind is the primary driving mechanism of massive outflows, we quantitatively compare outflow properties obtained through simulations and observations. Since the outflows obtained through simulations are slightly younger than those obtained through observations, the time-integrated quantities of outflow mass, momentum, and kinetic energy are slightly smaller than those obtained through observations. On the other hand, time-derivative quantities of mass ejection rate, outflow momentum flux, and kinetic luminosity obtained through simulations are in very good agreement with those obtained through observations. This indicates that the MHD disk wind greatly contributes to the massive outflow driving from high-mass protostars, and the magnetic field might significantly control the high-mass star formation process.
This is the first of a series of papers presenting the THROES (A caTalogue of HeRschel Observations of Evolved Stars) project, intended to provide a comprehensive overview of the spectroscopic results obtained in the far-infrared (55-670 {\mu}m) with the Her- schel space observatory on low-to-intermediate mass evolved stars in our Galaxy. Here we introduce the catalogue of interac- tively reprocessed PACS (Photoconductor Array Camera and Spectrometer) spectra covering the 55-200 {\mu}m range for 114 stars in this category for which PACS range spectroscopic data is available in the Herschel Science Archive (HSA). Our sample in- cludes objects spanning a range of evolutionary stages, from the asymptotic giant branch to the planetary nebula phase, displaying a wide variety of chemical and physical properties. The THROES/PACS catalogue is accessible via a dedicated web-based inter- face (https://throes.cab.inta-csic.es/) and includes not only the science-ready Herschel spectroscopic data for each source, but also complementary photometric and spectroscopic data from other infrared observatories, namely IRAS (Infrared Astronomical Satellite), ISO (Infrared Space Observatory) or AKARI, at overlapping wavelengths. Our goal is to create a legacy-value Herschel dataset that can be used by the scientific community in the future to deepen our knowledge and understanding of these latest stages of the evolution of low-to-intermediate mass stars.
Around the turn of the last century, star clusters of all kinds were considered "simple" stellar populations. Over the past decade, this situation has changed dramatically. At the same time, star clusters are among the brightest stellar population components and, as such, they are visible out to much greater distances than individual stars, even the brightest, so that understanding the intricacies of star cluster composition and their evolution is imperative for understanding stellar populations and the evolution of galaxies as a whole. In this review of where the field has moved to in recent years, we place particular emphasis on the properties and importance of binary systems, the effects of rapid stellar rotation, and the presence of multiple populations in Magellanic Cloud star clusters across the full age range. Our most recent results imply a reverse paradigm shift, back to the old simple stellar population picture for at least some intermediate-age (~1--3 Gyr-old) star clusters, opening up exciting avenues for future research efforts.
In late 2016 and early 2017 the flat spectrum radio quasar CTA 102 exhibited a very strong and long-lasting outburst. The event can be described by a roughly 2 months long increase of the baseline flux in the monitored energy bands (optical to $\gamma$ rays) by a factor 8, and a subsequent decrease over another 2 months back to pre-flare levels. The long-term trend was superseded by short but very strong flares, resulting in a peak flux that was a factor 50 above pre-flare levels in the $\gamma$-ray domain and almost a factor 100 above pre-flare levels in the optical domain. In this paper we explain the long-term evolution of the outburst by the ablation of a gas cloud penetrating the relativistic jet. The slice-by-slice ablation results in a gradual increase of the particle injection until the center of the cloud is reached, after which the injected number of particles decreases again. With reasonable cloud parameters we obtain excellent fits of the long-term trend.
Links to: arXiv, form interface, find, astro-ph, recent, 1711, contact, help (Access key information)