We present BOND, a Bayesian code to simultaneously derive oxygen and nitrogen abundances in giant H II regions. It compares observed emission lines to a grid of photoionization models without assuming any relation between O/H and N/O. Our grid spans a wide range in O/H, N/O and ionization parameter U, and covers different starburst ages and nebular geometries. Varying starburst ages accounts for variations in the ionizing radiation field hardness, which arise due to the ageing of H II regions or the stochastic sampling of the initial mass function. All previous approaches assume a strict relation between the ionizing field and metallicity. The other novelty is extracting information on the nebular physics from semi-strong emission lines. While strong lines ratios alone ([O III]/Hbeta, [O II]/Hbeta and [N II]/Hbeta) lead to multiple O/H solutions, the simultaneous use of [Ar III]/[Ne III] allows one to decide whether an H II region is of high or low metallicity. Adding He I/Hbeta pins down the hardness of the radiation field. We apply our method to H II regions and blue compact dwarf galaxies, and find that the resulting N/O vs O/H relation is as scattered as the one obtained from the temperature-based method. As in previous strong-line methods calibrated on photoionization models, the BOND O/H values are generally higher than temperature-based ones, which might indicate the presence of temperature fluctuations or kappa distributions in real nebulae, or a too soft ionizing radiation field in the models.
To investigate what drives the reversal of the morphology-density relation at intermediate/high redshift, we present a multi-wavelength analysis of 27 dusty starburst galaxies in the massive cluster Cl 0024+17 at z = 0.4. We combine H-alpha dynamical maps from the VLT/FLAMES multi-IFU system with far-infrared imaging using Herschel SPIRE and millimetre spectroscopy from IRAM/NOEMA, in order to measure the dynamics, star formation rates and gas masses of this sample. Most galaxies appear to be rotationally supported, with a median ratio of rotational support to line-of-sight velocity dispersion v/sigma ~ 5 +/- 2, and specific angular momentum lambda_R = 0.83 +/- 0.06 - comparable to field spirals of a similar mass at this redshift. The star formation rates of 3 - 26 M_solar/yr and average 12 CO derived gas mass of 1 x 10^10 M_solar suggest gas depletion timescales of ~ 1Gyr (~ 0.25 of the cluster crossing time). We derive characteristic dust temperatures (mean T_dust = 26 +/- 1 K) consistent with local galaxies of similar far-infrared luminosity, suggesting that the low density gas is yet to be stripped. Taken together, these results suggest that these starbursts have only recently accreted from the field, with star formation rates likely enhanced due to the effects of ram pressure. In order to make the transition to cluster S0s these galaxies must lose ~ 40% of their specific angular momentum. We suggest this must occur > 1 Gyr later, after the molecular gas has been depleted and/or stripped, via multiple tidal interactions with other cluster members.
$\Lambda$-Warm Dark Matter (WDM) has been proposed as alternative scenario to $\Lambda$ cold dark matter (CDM), motivated by discrepancies at the scale of dwarf galaxies, with less small-scale power and realized by collisionless particles with energies in the range $1-3$ keV. We present a new approach to constrain the viability of such WDM models using star formation histories of the dwarf spheroidal galaxies (dSphs) in the Local Group. We compare their high time-resolution star formation histories (SFHs) obtained with HST-based color magnitude diagrams with the range of possible collapse redshifts of their dark matter halos expected in CDM and in different WDM scenarios. The collapse redshift is inferred after determining a plausible infall mass of the subhalo. This is based on the current mass of individual dwarf inferred from stellar kinematics combined with results of cosmological simulations providing information on the subhalo evolution. Since WDM subhalos close to the filtering mass scale form significantly later than CDM, we show that they are in the first place difficult to reconcile with a truncation of star formation occurring as early as $z\geq 3$. The Ultra-Faint Dwarfs (UFDs) provide the most stringent constraints. Using 6 UFDs with the best determination of the SFHs, we show that we can exclude a 1 keV warm particle to a 2-$\sigma$ confidence interval consistently with other methods reported in the literature. For some objects the $2$ keV model is also excluded. We discuss the various caveats of the method, most notably the low number of dwarfs with accurately determined star formation histories and the uncertainties in the determination of the infall mass of the subhalos. Our preliminary analysis serves as a pathfinder for future investigations that will combine upcoming accurate SFHs for more local dSphs with direct analysis of WDM cosmological simulations with baryons.
The shapes of cluster central galaxies are not randomly oriented, but rather exhibit coherent alignments with the shapes of their parent clusters as well as with large-scale structure. In this work, we undertake a comprehensive study of the alignments of central galaxies at low redshift. Based on a sample of 8237 clusters and 94817 members in the redMaPPer cluster catalog with 0.1 < z < 0.35, we first quantify the alignment between the projected central galaxy shapes and the distribution of member satellites, to understand what central galaxy and cluster properties most strongly correlate with these alignments. Next, we investigate the angular segregation of satellites with respect to their central galaxy major axis directions, to identify the satellite properties that most strongly predict their angular segregation. We find that central galaxies are more aligned with their member galaxy distributions in clusters that are more elongated and have higher richness, and for central galaxies with larger physical size, higher luminosity and centering probability, and redder color. Satellites with redder color, higher luminosity, located closer to the central galaxy, and with smaller ellipticity show a stronger angular segregation toward their central galaxy major axes. Finally, we provide physical explanations for some of the identified correlations, and discuss the connection to theories of central galaxy alignments, the impact of primordial alignments with tidal fields, and the importance of anisotropic accretion.
To link the physical and star formation properties of structures ranging from
Giant Molecular Clouds (GMCs), to Molecular Cloud Complexes (MCCs), and to
Galaxies, we compare the mutual relations between their masses $M$, mass
surface densities $\Sigma_{M_{\rm gas}}$, radii $R$, velocity dispersions
$\sigma$, star formation rates $SFR$, and SFR densities $\Sigma_{\rm SFR}$
using data from the $^{12}$CO 1-0 CfA survey and from the literature. We derive
universal scaling relations for a comprehensive compilation of molecular cloud
structures, spanning 8 orders of magnitudes in size and 13 orders of magnitudes
in mass:
$\sigma\sim R^{0.47}$, $M\sim R^{1.96}$, $\Sigma_{\rm SFR}\sim \Sigma_{M_{\rm
gas}}^{1.37}$, ${SFR}\sim {M}^{0.87}$, and ${SFR}\sim {\sigma}^{2.66}$.
We also find that the slopes and the coefficients are different for
individual scales. Additionally, there is a break at the MCC scale in the
$\sigma-R$ relation and breaks between the starburst and the normal
star-forming objects in the $SFR-M$ and $\Sigma_{\rm SFR}$-$\Sigma_{\rm M_{\rm
gas}}$ relations.
The Schmidt-Kennicutt diagram is used to distinguish starburst from normal
star-forming structures by applying a $\Sigma_{M_{\rm gas}}$ threshold of
$100\,M_\odot$ pc$^{-2}$ and a $\Sigma_{\rm SFR}$ threshold of $1\,M_\odot$
yr$^{-1}$ kpc$^{-2}$. Mini-starburst complexes are MCCs that have enhanced
$\Sigma_{\rm SFR}$ ($>1\,M_\odot$ yr$^{-1}$ kpc$^{-2}$), probably caused by
dynamic events such as radiation pressure, colliding flows, or spiral arm
gravitational instability which compress material within the MCCs. Because of
the dynamical evolution, gravitational boundedness does not play a significant
role in characterizing the star formation activity of MCC, especially the
mini-starburst complex, which leads to the conclusion that the formation of
massive stars and clusters is dynamic.
There is evidence in 21cm HI emission for voids several kpc in size centered approximately on the Galactic centre, both above and below the Galactic plane. These appear to map the boundaries of the Galactic nuclear wind. An analysis of HI at the tangent points, where the distance to the gas can be estimated with reasonable accuracy, shows a sharp transition at Galactic radii $R\lesssim 2.4$ kpc from the extended neutral gas layer characteristic of much of the Galactic disk, to a thin Gaussian layer with FWHM $\sim 125$ pc. An anti-correlation between HI and $\gamma$-ray emission at latitudes $10^{\circ} \leq |b| \leq 20^{\circ}$ suggests that the boundary of the extended HI layer marks the walls of the Fermi Bubbles. With HI we are able to trace the edges of the voids from $|z| > 2$ kpc down to $z\approx0$, where they have a radius $\sim 2$ kpc. The extended HI layer likely results from star formation in the disk, which is limited largely to $R \gtrsim 3$ kpc, so the wind may be expanding into an area of relatively little HI. Because the HI kinematics can discriminate between gas in the Galactic center and foreground material, 21cm HI emission may be the best probe of the extent of the nuclear wind near the Galactic plane.
Radio continuum observations using the Australia telescope compact array at 5.5, 9.0, 17.0 and 22.8 GHz have detected free-free emission associated with 45 of 49 massive young stellar objects and HII regions. Of these, 26 sources are classified as ionized jets (12 of which are candidates), 2 as ambiguous jets or disc winds, 1 as a disc-wind, 14 as HII regions and 2 were unable to be categorised. Classification as ionized jets is based upon morphology, radio flux and spectral index, in conjunction with previous observational results at other wavelengths. Radio-luminosity and momentum are found to scale with bolometric luminosity in the same way as low-mass jets, indicating a common mechanism for jet production across all masses. In 13 of the jets, we see associated non-thermal/optically-thin lobes resulting from shocks either internal to the jet and/or at working surfaces. Ten jets display non-thermal (synchrotron emission) spectra in their lobes, with an average spectral index of -0.55 consistent with Fermi acceleration in shocks. This shows that magnetic fields are present, in agreement with models of jet formation incorporating magnetic fields. Since the production of collimated radio jets is associated with accretion processes, the results presented in this paper support the picture of disc-mediated accretion for the formation of massive stars with an upper-limit on the jet phase lasting approximately $6.5 \times 10^4 yr$. Typical mass loss rates in the jet are found to be $1.4 \times 10^{-5} M_\odot yr^{-1}$ with associated momentum rates of the order $(1-2) \times 10^{-2} M_\odot km s^{-1} yr^{-1}$.
How did galaxies form and evolve? This is one of the most challenging questions in astronomy to- day. Answering it requires a careful combination of observational and theoretical work to reliably determine the observed properties of cosmic bodies over large portions of the distant Universe on the one hand, and accurately model the physical processes driving their evolution on the other. Most importantly, it requires bringing together disparate multi-wavelength and multi-resolution spectro-photometric datasets in an homogeneous and well-characterized manner so that they are suitable for a rigorous statistical analysis. The Herschel Extragalactic Legacy Project (HELP) funded by the EC FP7 SPACE program aims to achieve this goal by combining the expertise of optical, infrared and radio astronomers to provide a multi-wavelength database for the dis- tant Universe as an accessible value-added resource for the astronomical community. It will do so by bringing together multi-wavelength datasets covering the 1000 deg2 mapped by Herschel extragalactic surveys and thus creating a joint lasting legacy from several ambitious sky surveys.
We present Herschel, ALMA Compact Array (ACA), and Caltech Submillimeter Observatory (CSO) observations of the prestellar core in L1689N, which has been suggested to be interacting with a molecular outflow driven by the nearby solar type protostar IRAS 16293-2422. This source is characterized by some of the highest deuteration levels seen in the interstellar medium. The change in the NH2D line velocity and width across the core provides clear evidence of an interaction with the outflow, traced by the high-velocity water emission. Quiescent, cold gas, characterized by narrow line widths is seen in the NE part of the core, while broader, more disturbed line profiles are seen in the W/SW part. Strong N2D+ and ND3 emission is detected with the ACA, extending S/SW from the peak of the single-dish NH2D emission. The ACA data also reveal the presence a compact dust continuum source, with a mean size of ~1100 au, a central density of (1-2) 10^7 cm-3, and a mass of 0.2-0.4 Msun. The dust emission peak is displaced ~5" to the south with respect to the N2D+ and ND3 emission, as well as the single-dish dust continuum peak, suggesting that the northern, quiescent part of the core is characterized by spatially extended continuum emission, which is resolved out by the interferometer. We see no clear evidence of fragmentation in this quiescent part of the core, which could lead to a second generation of star formation, although a weak dust continuum source is detected in this region in the ACA data.
We present a molecular line study towards 31 extended green object (EGO) clumps in the southern sky using data from MALT90 (Millimetre Astronomy Legacy Team 90 GHz). According to previous multiwavelength observations, we divide our sample into two groups: massive young stellar objects (MYSOs) and HII regions. Our results seem to support that N2H+ and C2H emissions mainly come from the gas inside quiescent clumps. In addition, we also find that the [N2H+]/[H13CO+] and [C2H]/[H13CO+] relative abundance ratios decrease from MYSOs to HII regions. These results suggest depletion of N2H+ and C2H in the late stages of massive-star formation, probably caused by the formation of HII regions inside. N2H+ and C2H might be used as chemical clocks for massive-star formation by comparing with other molecules such as H13CO+ and HC3N.
We search for hints to the origin and nature of compact stellar systems in the magnitude range of ultracompact dwarf galaxies in deep wide-field imaging data of the Fornax cluster core. We visually investigate a large sample of 355 spectroscopically confirmed cluster members with V-band equivalent magnitudes brighter than -10 mag for faint extended structures. Our data reveal peculiar compact stellar systems, which appear asymmetric or elongated from their outer light distribution. We characterize the structure of our objects by quantifying their core concentration, as well as their outer asymmetry and ellipticity. For the brighter objects of our sample we also investigate their spatial and phase-space distribution within the cluster. We argue that the distorted outer structure alone that is seen for some of our objects, is not sufficient to decide whether these systems have a star cluster or a galaxy origin. However, we find that objects with low core concentration and high asymmetry (or high ellipticity) are primarily located at larger cluster-centric distances as compared to the entire sample. This supports the hypothesis that at least some of these objects may originate from tidally stripped galaxies.
Whether the progenitors of Type-Ia Supernovae, single degenerate or double-degenerate white dwarf (WD) systems, is a highly debated topic. To address the origin of the Type Ia Tycho's supernova remnant (SNR), SN 1572, we have carried out a 12CO J=1-0 mapping and a 3-mm line survey towards the remnant using the IRAM 30 m telescope. We show that Tycho is surrounded by a clumpy molecular bubble at the local standard of rest velocity $\sim 61$ km s$^{-1}$ which expands at a speed $\sim 4.5$ km s$^{-1}$ and has a mass of $\sim 220$ $M_\odot$ (at the distance of 2.5 kpc). Enhanced 12CO J=2-1 line emission relative to 12CO J=1-0 emission and possible line broadenings (in velocity range -64-- -60 km s$^{-1}$) are found at the northeastern boundary of the SNR where the shell is deformed and decelerated. These features, combined with the morphological correspondence between the expanding molecular bubble and Tycho, suggest that the SNR is associated with the bubble at velocity range -66-- -57 km s$^{-1}$. The most plausible origin for the expanding bubble is the fast outflow (with velocity $> 100$ km s$^{-1}$) driven from the vicinity of a WD as it accreted matter from a non-degenerate companion star. The SNR has been expanding in the low-density wind-blown bubble and the shock wave has just reached the molecular cavity wall. This is the first unambiguous detection of the expanding bubble driven by the progenitor of the Type-Ia SNR, which constitutes evidence for a single degenerate progenitor for this Type-Ia supernova.
Neutral Hydrogen (HI) provides a very important fuel for star formation, but is difficult to detect at high redshift due to weak emission, limited sensitivity of modern instruments, and terrestrial radio frequency interference (RFI) at low frequencies. We the first attempt to use gravitational lensing to detect HI line emission from three gravitationally lensed galaxies behind the cluster Abell 773, two at redshift of 0.398 and one at z=0.487, using the Green Bank Telescope. We find a 3 sigma upper limit for a galaxy with a rotation velocity of 200 km/s is M_HI=6.58x10^9 and 1.5x10^10 M_solar at z=0.398 and z=0.487. The estimated HI masses of the sources at z=0.398 and z=0.487 are a factor of 3.7 and ~30 times lower than our detection limits at the respective redshifts. To facilitate these observations we have used sigma clipping to remove both narrow- and wide-band RFI but retain the signal from the source. We are able to reduce the noise of the spectrum by ~25% using our routine instead of discarding observations with too much RFI. The routine is most effective when ~10 of the integrations or fewer contains RFI. These techniques can be used to study HI in highly magnified distant galaxies that are otherwise too faint to detect.
We present Atacama Large Millimeter/submillimeter Array (ALMA) Cycle 3 observations of CO(2-1) emission from the circumnuclear disk in the E/S0 galaxy NGC 1332 at 0.044" resolution. The disk exhibits regular rotational kinematics and central high-velocity emission (+/-500 km/s) consistent with the presence of a compact central mass. We construct models for a thin, dynamically cold disk in the gravitational potential of the host galaxy and black hole, and fit the beam-smeared model line profiles directly to the ALMA data cube. Model fits successfully reproduce the disk kinematics out to r=200 pc. Fitting models just to spatial pixels within projected r=50 pc of the nucleus (two times larger than the black hole's gravitational radius of influence), we find M_BH=6.64(-0.63,+0.65)*10^8 solar masses. This observation demonstrates ALMA's powerful capability to determine the masses of supermassive black holes by resolving gas kinematics on small angular scales in galaxy nuclei.
We present the numerical code PRECESSION: a new open-source python module to study the dynamics of precessing black-hole binaries in the post-Newtonian regime. The code provides a comprehensive toolbox to (i) study the evolution of the black-hole spins along their precession cycles, (ii) perform gravitational-wave driven binary inspirals using both orbit-averaged and precession-averaged integrations, and (iii) predict the properties of the merger remnant through fitting formulae obtained from numerical-relativity simulations. PRECESSION is a ready-to-use tool to add the black-hole spin dynamics to larger-scale numerical studies such as gravitational-wave parameter estimation codes, population synthesis models to predict gravitational-wave event rates, galaxy merger trees and cosmological simulations of structure formation. PRECESSION provides fast and reliable integration methods to propagate statistical samples of black-hole binaries from/to large separations where they form to/from small separations where they become detectable, thus linking gravitational-wave observations of spinning black-hole binaries to their astrophysical formation history. The code is also a useful tool to compute initial parameters for numerical-relativity simulations targeting specific precessing systems. PRECESSION can be installed from the Python Package Index and it is freely distributed under version control on Github, where further documentation is provided.
We have previously calculated the intergalactic background light (IBL) as a function of redshift in the far ultraviolet to near infrared range, based purely on data from deep galaxy surveys. Here we utilize similar methods to determine the mid- and far infrared IBL out to a wavelength of 850 microns. Our approach enables us to constrain the range of photon densities, based on the uncertainties from observationally determined luminosity densities and colors. By also including the effect of the 2.7 K cosmic background photons, we determine 68% confidence upper and lower limits on the opacity of the universe to gamma-rays up to PeV energies. Our direct results on the IBL are consistent with those from complimentary gamma-ray analyses using observations from the Fermi $\gamma$-ray space telescope and the H.E.S.S. air Cherenkov telescope. Thus, we find no evidence of previously suggested processes for the modification of gamma-ray spectra other than that of absorption by pair production alone.
Links to: arXiv, form interface, find, astro-ph, recent, 1605, contact, help (Access key information)
We present BOND, a Bayesian code to simultaneously derive oxygen and nitrogen abundances in giant H II regions. It compares observed emission lines to a grid of photoionization models without assuming any relation between O/H and N/O. Our grid spans a wide range in O/H, N/O and ionization parameter U, and covers different starburst ages and nebular geometries. Varying starburst ages accounts for variations in the ionizing radiation field hardness, which arise due to the ageing of H II regions or the stochastic sampling of the initial mass function. All previous approaches assume a strict relation between the ionizing field and metallicity. The other novelty is extracting information on the nebular physics from semi-strong emission lines. While strong lines ratios alone ([O III]/Hbeta, [O II]/Hbeta and [N II]/Hbeta) lead to multiple O/H solutions, the simultaneous use of [Ar III]/[Ne III] allows one to decide whether an H II region is of high or low metallicity. Adding He I/Hbeta pins down the hardness of the radiation field. We apply our method to H II regions and blue compact dwarf galaxies, and find that the resulting N/O vs O/H relation is as scattered as the one obtained from the temperature-based method. As in previous strong-line methods calibrated on photoionization models, the BOND O/H values are generally higher than temperature-based ones, which might indicate the presence of temperature fluctuations or kappa distributions in real nebulae, or a too soft ionizing radiation field in the models.
To investigate what drives the reversal of the morphology-density relation at intermediate/high redshift, we present a multi-wavelength analysis of 27 dusty starburst galaxies in the massive cluster Cl 0024+17 at z = 0.4. We combine H-alpha dynamical maps from the VLT/FLAMES multi-IFU system with far-infrared imaging using Herschel SPIRE and millimetre spectroscopy from IRAM/NOEMA, in order to measure the dynamics, star formation rates and gas masses of this sample. Most galaxies appear to be rotationally supported, with a median ratio of rotational support to line-of-sight velocity dispersion v/sigma ~ 5 +/- 2, and specific angular momentum lambda_R = 0.83 +/- 0.06 - comparable to field spirals of a similar mass at this redshift. The star formation rates of 3 - 26 M_solar/yr and average 12 CO derived gas mass of 1 x 10^10 M_solar suggest gas depletion timescales of ~ 1Gyr (~ 0.25 of the cluster crossing time). We derive characteristic dust temperatures (mean T_dust = 26 +/- 1 K) consistent with local galaxies of similar far-infrared luminosity, suggesting that the low density gas is yet to be stripped. Taken together, these results suggest that these starbursts have only recently accreted from the field, with star formation rates likely enhanced due to the effects of ram pressure. In order to make the transition to cluster S0s these galaxies must lose ~ 40% of their specific angular momentum. We suggest this must occur > 1 Gyr later, after the molecular gas has been depleted and/or stripped, via multiple tidal interactions with other cluster members.
$\Lambda$-Warm Dark Matter (WDM) has been proposed as alternative scenario to $\Lambda$ cold dark matter (CDM), motivated by discrepancies at the scale of dwarf galaxies, with less small-scale power and realized by collisionless particles with energies in the range $1-3$ keV. We present a new approach to constrain the viability of such WDM models using star formation histories of the dwarf spheroidal galaxies (dSphs) in the Local Group. We compare their high time-resolution star formation histories (SFHs) obtained with HST-based color magnitude diagrams with the range of possible collapse redshifts of their dark matter halos expected in CDM and in different WDM scenarios. The collapse redshift is inferred after determining a plausible infall mass of the subhalo. This is based on the current mass of individual dwarf inferred from stellar kinematics combined with results of cosmological simulations providing information on the subhalo evolution. Since WDM subhalos close to the filtering mass scale form significantly later than CDM, we show that they are in the first place difficult to reconcile with a truncation of star formation occurring as early as $z\geq 3$. The Ultra-Faint Dwarfs (UFDs) provide the most stringent constraints. Using 6 UFDs with the best determination of the SFHs, we show that we can exclude a 1 keV warm particle to a 2-$\sigma$ confidence interval consistently with other methods reported in the literature. For some objects the $2$ keV model is also excluded. We discuss the various caveats of the method, most notably the low number of dwarfs with accurately determined star formation histories and the uncertainties in the determination of the infall mass of the subhalos. Our preliminary analysis serves as a pathfinder for future investigations that will combine upcoming accurate SFHs for more local dSphs with direct analysis of WDM cosmological simulations with baryons.
The shapes of cluster central galaxies are not randomly oriented, but rather exhibit coherent alignments with the shapes of their parent clusters as well as with large-scale structure. In this work, we undertake a comprehensive study of the alignments of central galaxies at low redshift. Based on a sample of 8237 clusters and 94817 members in the redMaPPer cluster catalog with 0.1 < z < 0.35, we first quantify the alignment between the projected central galaxy shapes and the distribution of member satellites, to understand what central galaxy and cluster properties most strongly correlate with these alignments. Next, we investigate the angular segregation of satellites with respect to their central galaxy major axis directions, to identify the satellite properties that most strongly predict their angular segregation. We find that central galaxies are more aligned with their member galaxy distributions in clusters that are more elongated and have higher richness, and for central galaxies with larger physical size, higher luminosity and centering probability, and redder color. Satellites with redder color, higher luminosity, located closer to the central galaxy, and with smaller ellipticity show a stronger angular segregation toward their central galaxy major axes. Finally, we provide physical explanations for some of the identified correlations, and discuss the connection to theories of central galaxy alignments, the impact of primordial alignments with tidal fields, and the importance of anisotropic accretion.
To link the physical and star formation properties of structures ranging from
Giant Molecular Clouds (GMCs), to Molecular Cloud Complexes (MCCs), and to
Galaxies, we compare the mutual relations between their masses $M$, mass
surface densities $\Sigma_{M_{\rm gas}}$, radii $R$, velocity dispersions
$\sigma$, star formation rates $SFR$, and SFR densities $\Sigma_{\rm SFR}$
using data from the $^{12}$CO 1-0 CfA survey and from the literature. We derive
universal scaling relations for a comprehensive compilation of molecular cloud
structures, spanning 8 orders of magnitudes in size and 13 orders of magnitudes
in mass:
$\sigma\sim R^{0.47}$, $M\sim R^{1.96}$, $\Sigma_{\rm SFR}\sim \Sigma_{M_{\rm
gas}}^{1.37}$, ${SFR}\sim {M}^{0.87}$, and ${SFR}\sim {\sigma}^{2.66}$.
We also find that the slopes and the coefficients are different for
individual scales. Additionally, there is a break at the MCC scale in the
$\sigma-R$ relation and breaks between the starburst and the normal
star-forming objects in the $SFR-M$ and $\Sigma_{\rm SFR}$-$\Sigma_{\rm M_{\rm
gas}}$ relations.
The Schmidt-Kennicutt diagram is used to distinguish starburst from normal
star-forming structures by applying a $\Sigma_{M_{\rm gas}}$ threshold of
$100\,M_\odot$ pc$^{-2}$ and a $\Sigma_{\rm SFR}$ threshold of $1\,M_\odot$
yr$^{-1}$ kpc$^{-2}$. Mini-starburst complexes are MCCs that have enhanced
$\Sigma_{\rm SFR}$ ($>1\,M_\odot$ yr$^{-1}$ kpc$^{-2}$), probably caused by
dynamic events such as radiation pressure, colliding flows, or spiral arm
gravitational instability which compress material within the MCCs. Because of
the dynamical evolution, gravitational boundedness does not play a significant
role in characterizing the star formation activity of MCC, especially the
mini-starburst complex, which leads to the conclusion that the formation of
massive stars and clusters is dynamic.
There is evidence in 21cm HI emission for voids several kpc in size centered approximately on the Galactic centre, both above and below the Galactic plane. These appear to map the boundaries of the Galactic nuclear wind. An analysis of HI at the tangent points, where the distance to the gas can be estimated with reasonable accuracy, shows a sharp transition at Galactic radii $R\lesssim 2.4$ kpc from the extended neutral gas layer characteristic of much of the Galactic disk, to a thin Gaussian layer with FWHM $\sim 125$ pc. An anti-correlation between HI and $\gamma$-ray emission at latitudes $10^{\circ} \leq |b| \leq 20^{\circ}$ suggests that the boundary of the extended HI layer marks the walls of the Fermi Bubbles. With HI we are able to trace the edges of the voids from $|z| > 2$ kpc down to $z\approx0$, where they have a radius $\sim 2$ kpc. The extended HI layer likely results from star formation in the disk, which is limited largely to $R \gtrsim 3$ kpc, so the wind may be expanding into an area of relatively little HI. Because the HI kinematics can discriminate between gas in the Galactic center and foreground material, 21cm HI emission may be the best probe of the extent of the nuclear wind near the Galactic plane.
Radio continuum observations using the Australia telescope compact array at 5.5, 9.0, 17.0 and 22.8 GHz have detected free-free emission associated with 45 of 49 massive young stellar objects and HII regions. Of these, 26 sources are classified as ionized jets (12 of which are candidates), 2 as ambiguous jets or disc winds, 1 as a disc-wind, 14 as HII regions and 2 were unable to be categorised. Classification as ionized jets is based upon morphology, radio flux and spectral index, in conjunction with previous observational results at other wavelengths. Radio-luminosity and momentum are found to scale with bolometric luminosity in the same way as low-mass jets, indicating a common mechanism for jet production across all masses. In 13 of the jets, we see associated non-thermal/optically-thin lobes resulting from shocks either internal to the jet and/or at working surfaces. Ten jets display non-thermal (synchrotron emission) spectra in their lobes, with an average spectral index of -0.55 consistent with Fermi acceleration in shocks. This shows that magnetic fields are present, in agreement with models of jet formation incorporating magnetic fields. Since the production of collimated radio jets is associated with accretion processes, the results presented in this paper support the picture of disc-mediated accretion for the formation of massive stars with an upper-limit on the jet phase lasting approximately $6.5 \times 10^4 yr$. Typical mass loss rates in the jet are found to be $1.4 \times 10^{-5} M_\odot yr^{-1}$ with associated momentum rates of the order $(1-2) \times 10^{-2} M_\odot km s^{-1} yr^{-1}$.
How did galaxies form and evolve? This is one of the most challenging questions in astronomy to- day. Answering it requires a careful combination of observational and theoretical work to reliably determine the observed properties of cosmic bodies over large portions of the distant Universe on the one hand, and accurately model the physical processes driving their evolution on the other. Most importantly, it requires bringing together disparate multi-wavelength and multi-resolution spectro-photometric datasets in an homogeneous and well-characterized manner so that they are suitable for a rigorous statistical analysis. The Herschel Extragalactic Legacy Project (HELP) funded by the EC FP7 SPACE program aims to achieve this goal by combining the expertise of optical, infrared and radio astronomers to provide a multi-wavelength database for the dis- tant Universe as an accessible value-added resource for the astronomical community. It will do so by bringing together multi-wavelength datasets covering the 1000 deg2 mapped by Herschel extragalactic surveys and thus creating a joint lasting legacy from several ambitious sky surveys.
We present Herschel, ALMA Compact Array (ACA), and Caltech Submillimeter Observatory (CSO) observations of the prestellar core in L1689N, which has been suggested to be interacting with a molecular outflow driven by the nearby solar type protostar IRAS 16293-2422. This source is characterized by some of the highest deuteration levels seen in the interstellar medium. The change in the NH2D line velocity and width across the core provides clear evidence of an interaction with the outflow, traced by the high-velocity water emission. Quiescent, cold gas, characterized by narrow line widths is seen in the NE part of the core, while broader, more disturbed line profiles are seen in the W/SW part. Strong N2D+ and ND3 emission is detected with the ACA, extending S/SW from the peak of the single-dish NH2D emission. The ACA data also reveal the presence a compact dust continuum source, with a mean size of ~1100 au, a central density of (1-2) 10^7 cm-3, and a mass of 0.2-0.4 Msun. The dust emission peak is displaced ~5" to the south with respect to the N2D+ and ND3 emission, as well as the single-dish dust continuum peak, suggesting that the northern, quiescent part of the core is characterized by spatially extended continuum emission, which is resolved out by the interferometer. We see no clear evidence of fragmentation in this quiescent part of the core, which could lead to a second generation of star formation, although a weak dust continuum source is detected in this region in the ACA data.
We present a molecular line study towards 31 extended green object (EGO) clumps in the southern sky using data from MALT90 (Millimetre Astronomy Legacy Team 90 GHz). According to previous multiwavelength observations, we divide our sample into two groups: massive young stellar objects (MYSOs) and HII regions. Our results seem to support that N2H+ and C2H emissions mainly come from the gas inside quiescent clumps. In addition, we also find that the [N2H+]/[H13CO+] and [C2H]/[H13CO+] relative abundance ratios decrease from MYSOs to HII regions. These results suggest depletion of N2H+ and C2H in the late stages of massive-star formation, probably caused by the formation of HII regions inside. N2H+ and C2H might be used as chemical clocks for massive-star formation by comparing with other molecules such as H13CO+ and HC3N.
We search for hints to the origin and nature of compact stellar systems in the magnitude range of ultracompact dwarf galaxies in deep wide-field imaging data of the Fornax cluster core. We visually investigate a large sample of 355 spectroscopically confirmed cluster members with V-band equivalent magnitudes brighter than -10 mag for faint extended structures. Our data reveal peculiar compact stellar systems, which appear asymmetric or elongated from their outer light distribution. We characterize the structure of our objects by quantifying their core concentration, as well as their outer asymmetry and ellipticity. For the brighter objects of our sample we also investigate their spatial and phase-space distribution within the cluster. We argue that the distorted outer structure alone that is seen for some of our objects, is not sufficient to decide whether these systems have a star cluster or a galaxy origin. However, we find that objects with low core concentration and high asymmetry (or high ellipticity) are primarily located at larger cluster-centric distances as compared to the entire sample. This supports the hypothesis that at least some of these objects may originate from tidally stripped galaxies.
Whether the progenitors of Type-Ia Supernovae, single degenerate or double-degenerate white dwarf (WD) systems, is a highly debated topic. To address the origin of the Type Ia Tycho's supernova remnant (SNR), SN 1572, we have carried out a 12CO J=1-0 mapping and a 3-mm line survey towards the remnant using the IRAM 30 m telescope. We show that Tycho is surrounded by a clumpy molecular bubble at the local standard of rest velocity $\sim 61$ km s$^{-1}$ which expands at a speed $\sim 4.5$ km s$^{-1}$ and has a mass of $\sim 220$ $M_\odot$ (at the distance of 2.5 kpc). Enhanced 12CO J=2-1 line emission relative to 12CO J=1-0 emission and possible line broadenings (in velocity range -64-- -60 km s$^{-1}$) are found at the northeastern boundary of the SNR where the shell is deformed and decelerated. These features, combined with the morphological correspondence between the expanding molecular bubble and Tycho, suggest that the SNR is associated with the bubble at velocity range -66-- -57 km s$^{-1}$. The most plausible origin for the expanding bubble is the fast outflow (with velocity $> 100$ km s$^{-1}$) driven from the vicinity of a WD as it accreted matter from a non-degenerate companion star. The SNR has been expanding in the low-density wind-blown bubble and the shock wave has just reached the molecular cavity wall. This is the first unambiguous detection of the expanding bubble driven by the progenitor of the Type-Ia SNR, which constitutes evidence for a single degenerate progenitor for this Type-Ia supernova.
Neutral Hydrogen (HI) provides a very important fuel for star formation, but is difficult to detect at high redshift due to weak emission, limited sensitivity of modern instruments, and terrestrial radio frequency interference (RFI) at low frequencies. We the first attempt to use gravitational lensing to detect HI line emission from three gravitationally lensed galaxies behind the cluster Abell 773, two at redshift of 0.398 and one at z=0.487, using the Green Bank Telescope. We find a 3 sigma upper limit for a galaxy with a rotation velocity of 200 km/s is M_HI=6.58x10^9 and 1.5x10^10 M_solar at z=0.398 and z=0.487. The estimated HI masses of the sources at z=0.398 and z=0.487 are a factor of 3.7 and ~30 times lower than our detection limits at the respective redshifts. To facilitate these observations we have used sigma clipping to remove both narrow- and wide-band RFI but retain the signal from the source. We are able to reduce the noise of the spectrum by ~25% using our routine instead of discarding observations with too much RFI. The routine is most effective when ~10 of the integrations or fewer contains RFI. These techniques can be used to study HI in highly magnified distant galaxies that are otherwise too faint to detect.
We present Atacama Large Millimeter/submillimeter Array (ALMA) Cycle 3 observations of CO(2-1) emission from the circumnuclear disk in the E/S0 galaxy NGC 1332 at 0.044" resolution. The disk exhibits regular rotational kinematics and central high-velocity emission (+/-500 km/s) consistent with the presence of a compact central mass. We construct models for a thin, dynamically cold disk in the gravitational potential of the host galaxy and black hole, and fit the beam-smeared model line profiles directly to the ALMA data cube. Model fits successfully reproduce the disk kinematics out to r=200 pc. Fitting models just to spatial pixels within projected r=50 pc of the nucleus (two times larger than the black hole's gravitational radius of influence), we find M_BH=6.64(-0.63,+0.65)*10^8 solar masses. This observation demonstrates ALMA's powerful capability to determine the masses of supermassive black holes by resolving gas kinematics on small angular scales in galaxy nuclei.
We present the numerical code PRECESSION: a new open-source python module to study the dynamics of precessing black-hole binaries in the post-Newtonian regime. The code provides a comprehensive toolbox to (i) study the evolution of the black-hole spins along their precession cycles, (ii) perform gravitational-wave driven binary inspirals using both orbit-averaged and precession-averaged integrations, and (iii) predict the properties of the merger remnant through fitting formulae obtained from numerical-relativity simulations. PRECESSION is a ready-to-use tool to add the black-hole spin dynamics to larger-scale numerical studies such as gravitational-wave parameter estimation codes, population synthesis models to predict gravitational-wave event rates, galaxy merger trees and cosmological simulations of structure formation. PRECESSION provides fast and reliable integration methods to propagate statistical samples of black-hole binaries from/to large separations where they form to/from small separations where they become detectable, thus linking gravitational-wave observations of spinning black-hole binaries to their astrophysical formation history. The code is also a useful tool to compute initial parameters for numerical-relativity simulations targeting specific precessing systems. PRECESSION can be installed from the Python Package Index and it is freely distributed under version control on Github, where further documentation is provided.
We have previously calculated the intergalactic background light (IBL) as a function of redshift in the far ultraviolet to near infrared range, based purely on data from deep galaxy surveys. Here we utilize similar methods to determine the mid- and far infrared IBL out to a wavelength of 850 microns. Our approach enables us to constrain the range of photon densities, based on the uncertainties from observationally determined luminosity densities and colors. By also including the effect of the 2.7 K cosmic background photons, we determine 68% confidence upper and lower limits on the opacity of the universe to gamma-rays up to PeV energies. Our direct results on the IBL are consistent with those from complimentary gamma-ray analyses using observations from the Fermi $\gamma$-ray space telescope and the H.E.S.S. air Cherenkov telescope. Thus, we find no evidence of previously suggested processes for the modification of gamma-ray spectra other than that of absorption by pair production alone.
Links to: arXiv, form interface, find, astro-ph, recent, 1605, contact, help (Access key information)
In the list of young stellar objects compiled by Megeath et al. (2012) for the Orion A molecular cloud, only 44 out of 1208 sources found projected onto low extinction (Ak<0.8 mag) gas are identified as protostars. These objects are puzzling because protostars are not typically expected to be associated with extended low extinction material. Here, we use high resolution extinction maps generated from Herschel data, optical/infrared and Spitzer Space Telescope photometry and spectroscopy of the low extinction protostellar candidate sources to determine if they are likely true protostellar sources or contaminants. Out of 44 candidate objects, we determine that 10 sources are likely protostars, with the rest being more evolved young stellar objects (18), galaxies (4), false detections of nebulosity and cloud edges (9), or real sources for which more data are required to ascertain their nature (3). We find none of the confirmed protostars to be associated with recognizable dense cores and we briefly discuss possible origins for these orphaned objects.
We have investigated the impact of photoionization and radiation pressure on a dusty star-forming cloud by one-dimensional radiation hydrodynamic simulations, which include absorption and re-emission of photons by dust. We find that even in a moderately dusty cloud with the infrared optical depth of 0.15, radiation pressure has strong impact on driving an outflow, while the effect of radiation pressure is negligible in a dustless cloud. The radiation pressure on dust creates an HII region whose density is much lower than that in a dustless cloud where an outflow is driven by thermal pressure of ionized gas. Due to the radiation pressure, a shocked shell expands with high velocity, > 100 km s^-1. Absorption of re-emitted photons by dust plays a significant role in driving an outflow when the infrared optical depth becomes unity and it increases the importance of radiation pressure. The column density of clouds decreases with very short timescale owing to the shell expansion. Because of the decline of the infrared optical depth, the radiation pressure becomes less important as the shell expands. This makes it difficult to incorporate the effects of radiation feedback as subgrid physics in simulations which do not solve radiative transfer directly.
Using a sample of 98 galaxy clusters recently imaged in the near infra-red with the ESO NTT, WIYN and WHT telescopes, supplemented with 33 clusters from the ESO archive, we measure how the stellar mass of the most massive galaxies in the universe, namely Brightest Cluster Galaxies (BCG), increases with time. Most of the BCGs in this new sample lie in the redshift range $0.2<z<0.6$, which has been noted in recent works to mark an epoch over which the growth in the stellar mass of BCGs stalls. From this sample of 132 clusters, we create a subsample of 102 systems that includes only those clusters that have estimates of the cluster mass. We combine the BCGs in this subsample with BCGs from the literature, and find that the growth in stellar mass of BCGs from 10 billion years ago to the present epoch is broadly consistent with recent semi-analytic and semi-empirical models. As in other recent studies, tentative evidence indicates that the stellar mass growth rate of BCGs may be slowing in the past 3.5 billion years. Further work in collecting larger samples, and in better comparing observations with theory using mock images is required if a more detailed comparison between the models and the data is to be made.
The evolution of dwarf satellites of the Milky Way is affected by the combination of ram pressure and tidal stripping, and internal feedback from massive stars. We investigate gas loss processes in the smallest satellites of the Milky Way using three-dimensional, high resolution, idealized wind tunnel simulations, accounting for gas loss through both ram pressure stripping and expulsion by supernova feedback. Using initial conditions appropriate for a dwarf galaxy like Leo T, we investigate whether or not environmental gas stripping and internal feedback can quench these low mass galaxies on the expected timescales, shorter than 2 Gyr. We find that supernova feedback contributes negligibly to the stripping rate for these low star formation rate galaxies. However, we also find that ram pressure stripping is less efficient than expected in the stripping scenarios we consider. Our work suggests that, although ram pressure stripping can eventually completely strip these galaxies, other physics is likely at play to reconcile our computed stripping times with the rapid quenching timescales deduced from observations of low mass Milky Way dwarf galaxies. We discuss the roles additional physics may play in this scenario, including host-satellite tidal interactions, cored vs. cuspy dark matter profiles, reionization, and satellite pre-processing. We conclude that a proper accounting of these physics together is necessary to understand the quenching of low mass Milky Way satellites.
We present the metallicity distribution of a sample of 471 RR Lyrae (RRL) stars in the Sculptor dSph, obtained from the $I$-band Period-Luminosity relation. It is the first time that the early chemical evolution of a dwarf galaxy is characterized in such a detailed and quantitative way, using photometric data alone. We find a broad metallicity distribution (FWHM=0.8 dex) that is peaked at [Fe/H]$\simeq$-1.90 dex, in excellent agreement with literature values obtained from spectroscopic data. Moreover, we are able to directly trace the metallicity gradient out to a radius of $\sim$55 arcmin. We find that in the outer regions (r$>\sim$32 arcmin) the slope of the metallicity gradient from the RRLs (-0.025 dex arcmin$^{-1}$) is comparable to the literature values based on red giant (RG) stars. However, in the central part of Sculptor we do not observe the latter gradients. This suggests that there is a more metal-rich and/or younger population in Sculptor that does not produce RRLs. This scenario is strengthened by the observation of a metal-rich peak in the metallicity distribution of RG stars by other authors, which is not present in the metallicity distribution of the RRLs within the same central area.
We investigate the presence and importance of dark matter discs in a sample of 24 simulated Milky Way galaxies in the APOSTLE project, part of the EAGLE programme of hydrodynamic simulations in Lambda-CDM cosmology. It has been suggested that a dark disc in the Milky Way may boost the dark matter density and modify the velocity modulus relative to a smooth halo at the position of the Sun, with ramifications for direct detection experiments. From a kinematic decomposition of the dark matter and a real space analysis of all 24 halos, we find that only one of the simulated Milky Way analogues has a detectable dark disc component. This unique event was caused by a merger at late time with an LMC-mass satellite at very low grazing angle. Considering that even this rare scenario only enhances the dark matter density at the solar radius by 35% and affects the high energy tail of the dark matter velocity distribution by less than 1%, we conclude that the presence of a dark disc in the Milky Way is unlikely, and is very unlikely to have a significant effect on direct detection experiments.
In this second paper on the entire virial region of the relaxed fossil cluster RXJ1159+5531, we present a hydrostatic analysis of the hot intracluster medium (ICM). For a model consisting of ICM, stellar mass from the central galaxy (BCG), and an NFW dark matter (DM) halo, we obtain good descriptions of the projected radial profiles of ICM emissivity and temperature. The BCG stellar mass is clearly detected with M_star/L_K = 0.61 +/- 0.11 solar, consistent with stellar population synthesis models for a Milky-Way IMF. We obtain a halo concentration, c_200 =8.4 +/- 1.0, and virial mass, M_200 = 7.9 +/- 0.6 x 10^13 M_sun. For its mass, the inferred concentration is larger than most relaxed halos produced in cosmological simulations with Planck parameters, consistent with RXJ1159+5531 forming earlier than the general halo population. The baryon fraction at r_200, f_b,200 = 0.134 +/- 0.007, is slightly below the Planck value (0.155) for the universe. When we account for the stellar baryons associated with non-central galaxies and the uncertain intracluster light, f_b,200 increases by ~0.015, consistent with the cosmic value. Performing our analysis in the context of MOND still requires a large DM fraction (85.0% +/- 2.5% at r=100 kpc) similar to that obtained using the standard Newtonian approach. The detection of a plausible stellar BCG mass component distinct from the NFW DM halo in the total gravitational potential suggests that ~10^14 M_sun represents the mass scale above which dissipation is unimportant in the formation of the central regions of galaxy clusters. (Abridged)
By adopting the empirical constraints related to the estimates of Helium enhancement ($\Delta Y$), present mass ratio between first and second stellar generations ($M_{1G}/M_{2G}$) and the actual mass of Galactic globular clusters ($M_{GC}$), we envisage a possible scenario for the formation of these stellar systems. Our approach allows for the possible loss of stars through evaporation or tidal interactions and different star formation efficiencies. In our approach the star formation efficiency of the first generation ($\epsilon_{1G}$) is the central factor that links the stellar generations as it not only defines both the mass in stars of the first generation and the remaining mass available for further star formation, but it also fixes the amount of matter required to contaminate the second stellar generation. In this way, $\epsilon_{1G}$ is fully defined by the He enhancement between successive generations in a GC. We also show that globular clusters fit well within a $\Delta Y$ {\it vs} $M_{1G}/M_{2G}$ diagram which indicates three different evolutionary paths. The central one is for clusters that have not loss stars, through tidal interactions, from either of their stellar generations, and thus their present $M_{GC}$ value is identical to the amount of low mass stars ($M_* \le$ 1 M$_\odot$) that resulted from both stellar generations. Other possible evolutions imply either the loss of first generation stars or the combination of a low star formation efficiency in the second stellar generation and/or a loss of stars from the second generation. From these considerations we derive a lower limit to the mass ($M_{tot}$) of the individual primordial clouds that gave origin to globular clusters.
The tadpole galaxy Kiso 5639 has a slowly rotating disk with a drop in metallicity at its star-forming head, suggesting that star formation was triggered by the accretion of metal-poor gas. We present multi-wavelength HST WFC3 images of UV through I band plus Halpha to search for peripheral emission and determine the properties of various regions. The head has a mass in young stars of ~10^6 Mo and an ionization rate of 6.4x10^51 s^{-1}, equivalent to ~2100 O9-type stars. There are four older star-forming regions in the tail, and an underlying disk with a photometric age of ~1 Gyr. The mass distribution function of 61 star clusters is a power law with a slope of -1.73+-0.51. Fourteen young clusters in the head are more massive than 10^4 Mo, suggesting a clustering fraction of 30%-45%. Wispy filaments of Halpha emission and young stars extend away from the galaxy. Shells and holes in the head HII region could be from winds and supernovae. Gravity from the disk should limit the expansion of the HII region, although hot gas might escape through the holes. The star formation surface density determined from Halpha in the head is compared to that expected from likely pre-existing and accreted gas. Unless the surface density of the accreted gas is a factor of ~3 or more larger than what was in the galaxy before, the star formation rate has to exceed the usual Kennicutt-Schmidt rate by a factor of >5.
The radiative torque (RAT) alignment of interstellar grains with ordinary paramagnetic susceptibilities has been supported by a number of earlier studies. The alignment of such grains depends on the so-called RAT parameter $q^{\max}$ that is determined by the grain shape. For interstellar grains with a broad range of $q^{\max}$, a significant fraction of grains is expected to get aligned with low angular momentum at the so-called low-J attractor points, which entail degrees of alignment between 20 or 30 percent, irrespectively of the strength of RATs. The latter value may not be sufficient for explaining the observed interstellar alignment in the diffuse medium. In this paper, we elaborate our model of radiative alignment for grains with enhanced magnetic susceptibility due to magnetic inclusions, such that both Magnetic torque and RAdiative Torque (MRAT) play a role in grain alignment. Such grains can get aligned with high angular momentum at the so-called high-J attractor points, which achieve a high degree of alignment. Using our analytical model (AMO) of RATs we derive the critical value of the magnetic relaxation parameter $\delta_{m}$ to produce high-J attractor points as functions of $q^{\max}$ and the anisotropic radiation angle relative to the magnetic field $\psi$. To calculate degree of grain alignment, we carry out numerical simulations of MRAT alignment by including stochastic excitations by gas collisions and magnetic fluctuations. We numerically demonstrate that the degree of MRAT alignment varies with $\psi$ and increases with increasing grain size for different values of $\delta_{m}$. Our obtained results pave the way for physical modeling of polarization spectrum of thermal dust emission as well as magnetic dipole emission from grains of arbitrary magnetic susceptibilities.
We present a kinematic catalogue for 21 M51-type galaxies. It consists of radial velocity distributions both from main and satellite components, along different position angles, which we obtained from long-slit spectroscopy. We detect deviations from circular motion in most of the main galaxies of each pair, due to the gravitational perturbation produced by the satellite galaxy. However somesystems do not show significant distortions in their radial velocity curves. We found some differences between the directions of photometric and kinematic major axes in main galaxies with a bar subsystem. The Tully-Fisher relation in the B-band and Ks-band for the present sample of M51-type systems is flatter when compared with isolated galaxies. Using the radial velocity data set, we built a synthetic normalized radial velocity distribution, as a reference for future modeling of these peculiar systems. The relative position angles between main galaxy major axis and companion location, as well as the velocity difference amplitude, indicate that the orbital motion of the satellite has a large projection on the main galaxy equatorial plane and that the satellite orbital motion is within the range of amplitudes of the main galaxy rotation curve. All the M51-type systems studied here except for one, are gravitationally bounded.
We have analyzed the frequency and properties of the nuclear activity in a sample of galaxies with circumnuclear rings and spirals (CNRs). This sample was compared with a control sample of galaxies with very similar global properties but without circumnuclear rings. We discuss the relevance of the results in regard to the AGN feeding processes and present the following results: (i) bright companion galaxies seem not to be important for the appearance of CNRs, which appear to be more related to intrinsic properties of the host galaxies or to minor merger processes; (ii) the proportion of strong bars in galaxies with an AGN and a CNR is somewhat higher than the expected ratio of strongly barred AGN galaxies from the results of Ho and co-workers; (iii) the incidence of Seyfert activity coeval with CNRs is clearly larger than the rate expected from the morphological distribution of the host galaxies; (iv) the rate of Sy 2 to Sy 1 type galaxies with CNRs is about three times larger than the expected ratio for galaxies without CNRs and is opposite to that predicted by the geometric paradigm of the classical unified model for AGNs, although it does support the hy-pothesis that Sy 2 activity is linked to circumnuclear star formation. The possible selection effects of the sample are discussed, and we conclude that the detected trends are strong enough to justify high quality observations of as large as possible sets of galaxies with circumnuclear rings and their matched control samples.
We present a dynamical model for gas transport, star formation, and winds in the nuclear regions of galaxies, focusing on the Milky Way's Central Molecular Zone (CMZ). In our model angular momentum and mass are transported by a combination of gravitational and bar-driven acoustic instabilities. In gravitationally-unstable regions the gas can form stars, and the resulting feedback drives both turbulence and a wind that ejects mass from the CMZ. We show that the CMZ is in a quasi-steady state where mass deposited at large radii by the bar is transported inward to a star-forming, ring-shaped region at $\sim 100$ pc from the Galactic Centre, where the shear reaches a minimum. This ring undergoes episodic starbursts, with bursts lasting $\sim 5-10$ Myr occurring at $\sim 20-40$ Myr intervals. During quiescence the gas in the ring is not fully cleared, but is driven out of a self-gravitating state by the momentum injected by expanding supernova remnants. Starbursts also drive a wind off the star-forming ring, with a time-averaged mass flux comparable to the star formation rate. We show that our model agrees well with the observed properties of the CMZ, and places it near a star formation minimum within the evolutionary cycle. We argue that such cycles of bursty star formation and winds should be ubiquitous in the nuclei of barred spiral galaxies, and show that the resulting distribution of galactic nuclei on the Kennicutt-Schmidt relation is in good agreement with that observed in nearby galaxies.
We have compiled a catalogue of HII regions detected with the Murchison Widefield Array (MWA) between 72 and 231MHz. The multiple frequency bands provided by the MWA allow us identify the characteristic spectrum generated by the thermal Bremsstrahlung process in HII regions. We detect 306 HII regions between 260 < l < 340 and report on the positions, sizes, peak, integrated flux density, and spectral indices of these HII regions. By identifying the point at which HII regions transition from the optically thin to thick regime we derive the physical properties including the electron density, ionised gas mass and ionising photon flux, towards 61 HII regions. This catalogue of HII regions represents the most extensive and uniform low frequency survey of HII regions in the Galaxy to date.
Emission-line galaxies (ELGs) are one of the main tracers of the large-scale structure to be targeted by the next-generation dark energy surveys. To provide a better understanding of the properties and statistics of these galaxies, we have collected spectroscopic data from the VVDS and DEEP2 deep surveys and estimated the galaxy luminosity functions (LFs) of three distinct emission lines, [OII], H$\beta$ and [OIII] at redshifts ($0.2 < z < 1.3$). Our measurements are based on the largest sample so far. We present the first measurement of the \Hb LF at these redshifts. We have also compiled LFs from the literature that were based on independent data or covered different redshift ranges, and we fit the entire set over the whole redshift range with analytic Schechter and Saunders models, assuming a natural redshift dependence of the parameters. We find that the characteristic luminosity ($L_*$) and density ($\phi_*$) of all LFs increase with redshift. Using the Schechter model, we find that $L^*$ of [OII] emitters increase by a factor 12 between redshift 0.2 and 1.3, that $L^*$ of [OIII] emitters increase by a factor 16 between redshift 0.3 and 1 and that $L^*$ of H$\beta$ emitters increase by a factor 13 between redshift 0.3 and 0.8.
Globular clusters which exhibit chemical and dynamical complexity have been suggested to be the stripped nuclei of dwarf galaxies (e.g., M54, $\omega$ Cen). We use $N$-body simulations of nuclear star clusters forming via the mergers of star clusters to explore the persistence of substructure in the phase space. We find that the observed level of differentiation is difficult to reconcile with the observed if nuclear clusters form wholly out of the mergers of star clusters. Only the star clusters that merged most recently retain sufficiently distinct kinematics to be distinguishable from the rest of the nuclear cluster though the critical factor is the number of merger events not the elapsed time. In situ star formation must therefore be included to explain the observed properties of nuclear star clusters, in good agreement with previous results.
Orientation, together with accretion and evolution, is one of the three main drivers in the Grand Unification of Active Galactic Nuclei (AGN). Being unresolved, determining the true inclination of those powerful sources is always difficult and indirect, yet it remains a vital clue to apprehend the numerous, panchromatic and complex spectroscopic features we detect. There are only a hundred inclinations derived so far; in this context, can we be sure that we measure the true orientation of AGN? To answer this question, four methods to estimate the nuclear inclination of AGN are investigated and compared to inclination-dependent observables (hydrogen column density, Balmer linewidth, optical polarization, and flux ratios within the IR and relative to X-rays). Among these orientation indicators, the method developed by Fisher, Crenshaw, Kraemer et al., mapping and modeling the radial velocities of the [O iii] emission region in AGN, is the most successful. The [O iii]-mapping technique shows highly statistically significant correlations at >95% confidence level for rejecting null hypothesis for all the test cases. Such results confirm that the Unified Model is correct at a scale ranging from kiloparsec to a fraction of a parsec. However, at a radial distance less than 0.01 pc from the central black hole, warps and misalignments may change this picture.
Could the velocity spread, increasing with time, in the Galactic disk be
explained as a result of gravitational interactions of stars with giant
molecular clouds (GMCs) and spiral arms? Do the old open clusters high above
the Galactic plane provide clues to this question? We explore the effects on
stellar orbits of scattering by inhomogeneities in the Galactic potential due
to GMCs, spiral arms and the Galactic bar, and whether high-altitude clusters
could have formed in orbits closer to the Galactic plane and later been
scattered.
Simulations of test-particle motions are performed in a realistic Galactic
potential. The effects of the internal structure of GMCs are explored. The
destruction of clusters in GMC collisions is treated in detail with N-body
simulations of the clusters.
The observed velocity dispersions of stars as a function of time are well
reproduced. The GMC structure is found to be significant, but adequate models
produce considerable scattering effects. The fraction of simulated massive old
open clusters, scattered into orbits with |z| > 400 pc, is typically 0:5%, in
agreement with the observed number of high-altitude clusters and consistent
with the present formation rate of massive open clusters.
The heating of the thin Galactic disk is well explained by gravitational
scattering by GMCs and spiral arms, if the local correlation between the GMC
mass and the corresponding voids in the gas is not very strong. Our results
suggest that the high-altitude metal-rich clusters were formed in orbits close
to the Galactic plane and later scattered to higher orbits. It is possible,
though not very probable, that the Sun formed in such a cluster before
scattering occurred.
In this paper we investigate the mass assembly history of ultramassive (Mstar > 10^11Msun) dense (Sigma = Mstar/(2*pi*Re^2) > 2500 Msun/pc^2) early-type galaxies (ETGs) over the last 9 Gyr. We have traced the evolution of the number density rho of ultramassive dense ETGs and have compared their structural (effective radius Re and stellar mass Mstar) and dynamical (velocity dispersion sigma_e) parameters over the redshift range 0 < z < 1.6. We have derived the number density at 1.6 < z < 1 from the MUNICS and GOODS-South surveys, while we have used the COSMOS and SDSS spectroscopic surveys to probe the intermediate and local redshift range. For the comparison of the dynamical and structural parameters, we have collected the ultramassive dense ETGs at 1.2 < z < 1.6 for which velocity dispersion measurements are available (11 ETGs). For 4 of them we present unpublished estimates of sigma_e. We probe the intermediate redshift range, and the local universe using the samples of ETGs by Saglia et al. (2010), Zahid et al. (2015), and by Thomas et al. (2010). We find that the number density of ultramassive dense ETGs evolves as rho(z) = K*(1 + z)^(0.3\pm0.8) implying a decrease of ~ 25% of the population since z = 1.6. By comparing the values of Re, Mstar, and sigma_e of ultramassive dense ETGs over the range 0 < z < 1.6 we find that all the high-z ETGs have a counterpart in the local universe. This implies either that the majority (~70%) of ultramassive dense ETGs has already completed its assembly and its shaping at <z> = 1.4, or that, if a significant fraction of them evolves in size, new ultramassive dense ETGs must form at z < 1.5 to maintain their number density almost constant. The difficulty into identify good progenitors for these new dense ETGs at z < 1.5, and the stellar populations properties of local ultramassive dense ETGs point toward the first hypothesis.
To explore the role of titanium- and calcium-dust depletion in gas in and around galaxies we systematically study Ti/Ca abundance ratios in intervening absorption-line systems at low and high redshift. We investigate high-resolution optical spectra obtained by the UVES instrument at the Very Large Telescope and spectroscopically analyze 34 absorption-line systems at z<=0.5 to measure column densities (or limits) for CaII and TiII. We complement our UVES data set with previously published absorption-line data on Ti/Ca for redshifts up to z~3.8. Our absorber sample contains 110 absorbers (DLAs, sub-DLAs & LLSs). We compare our Ti/Ca findings with results from the Milky Way and the Magellanic Clouds and discuss the properties of Ti/Ca absorbers in the general context of quasar absorption-line systems.Our analysis indicates that there are two distinct populations of absorbers with either high or low Ti/Ca ratios with a separation at [Ti/Ca}]~1. While the calcium dust depletion in most of the absorbers appears to be severe, the titanium depletions are mild in systems with high Ti/Ca ratios. The derived trend indicates that absorbers with high Ti/Ca ratios have dust-to-gas ratios that are substantially lower than in the Milky Way. We characterize the overall nature of the absorbers by correlating Ti/Ca with other observables (e.g., metallicity, velocity-component structure) and by modeling the ionization properties of singly-ionized Ca and Ti in different environments. We conclude that CaII and TiII bearing absorption-line systems trace predominantly neutral gas in the disks and inner halo regions of galaxies, where the abundance of Ca and Ti reflects the local metal and dust content of the gas. Our study suggests that the Ti/Ca ratio represents a useful measure for the gas-to-dust ratio and overall metallicity in intervening absorption-line systems.
The formation mechanism of brown dwarfs (BDs) is one of the long-standing problems in star formation because the typical Jeans mass in molecular clouds is too large to form these substellar objects. To answer this question, it is crucial to study a BD at the embedded phase. IRAS 16253-2429 is classified as a very low luminosity object (VeLLO) with internal luminosity 0.1 Lsun. VeLLOs are believed to be very low-mass protostars or even proto-BDs. We observed the jet/outflow driven by IRAS 16253-2429 in CO (2-1), (6-5), and (7-6) using the IRAM 30 m and APEX telescopes and the SMA in order to study its dynamical features and physical properties. Our SMA map reveals two protostellar jets, indicating the existence of a proto-binary system as implied by the precessing jet detected in H2 emission. We detect a wiggling pattern in the position-velocity diagrams along the jet axes, which is likely due to the binary orbital motion. Based on this, we derive the current mass of the binary as ~0.032 Msun. Given the low envelope mass, IRAS 16253-2429 will form a binary that probably consist of one or two BDs. Furthermore, we found that the outflow force as well as the mass accretion rate are very low based on the multi-transition CO observations, which suggests that the final masses of the binary components are at the stellar/substellar boundary. Since IRAS 16253 is located in an isolated environment, we suggest that BDs can form through fragmentation and collapse like low-mass stars.
Context: In the last years, there have been many studies on the omnipresence and structures of filaments in star-forming regions, as well as their role in the process of star formation. Those filaments are normally identified as elongated fibres across the plane of the sky. But how would we detect filaments that are inclined? Aims: We aim to learn more about whether, and how, total column density or dust temperature change with respect to the line of sight. Such variations would enable observers to use dust observations to identify and study filaments at any inclination and gain more insight on the distribution and orientations of filaments within the Galactic plane. Methods: As a first step, we perform numerical calculations on simple cylindrical models to evaluate the influence of filament geometry on the average flux density. After that, we apply our three-dimensional Monte Carlo dust radiative transfer code on two models of star-forming regions and derive maps of effective total column density and dust temperature at different viewing angles. Results: We see only slight changes of average flux density for all cylinders we study. For our more complex models, we find that the effective dust temperature is not sensitive to viewing angle, while the total column density is strongly influenced, with differences exceeding an order of magnitude. The variations are not injective with the viewing angle and depend on the structure of the object. Conclusions: We conclude that there is no single quantity in our analysis that can uniquely trace the inclination and three-dimensional structure of a filament based on dust observations alone. However, observing wide variations in total column density at a given effective dust temperature is indicative of inclined filaments.
We here report the properties of Wolf-Rayet (W-R) stars in 14 locations in the nearby spiral galaxy M81. These locations were found serendipitously while analysing the slit spectra of a sample of ~150 star-forming complexes, taken using the long-slit and Multi-Object spectroscopic modes of the OSIRIS instrument at the 10.4-m Gran Telescopio Canarias. Colours and magnitudes of the identified point sources in the Hubble Space Telescope images compare well with those of individual W-R stars in the Milky Way. Using templates of individual W-R stars, we infer that the objects responsible for the observed W-R features are single stars in 12 locations, comprising of 3 WNLs, 3 WNEs, 2 WCEs and 4 transitional WN/C types. In diagrams involving bump luminosities and the width of the bumps, the W-R stars of the same sub-class group together, with the transitional stars occupying locations intermediate between the WNE and WCE groups, as expected from the evolutionary models. However, the observed number of 4 transitional stars out of our sample of 14 is statistically high as compared to the 4% expected in stellar evolutionary models.
We study the cluster environment for a sample of 21 radio loud AGN from the 3CR catalog at z>1, 12 radio galaxies and 9 quasars with HST images in the optical and IR. We use two different approaches to determine cluster candidates. We identify the early type galaxies (ETGs) in every field by modeling each of the sources within a 40" radius of the targets with a Sersic profile. Using a simple passive evolution model, we derive the expected location of the ETGs on the red sequence (RS) in the color-magnitude diagram for each of the fields of our sources. For seven targets, the model coincides with the position of the ETGs. A second approach involves a search for over densities. We compare the object densities of the sample as a whole and individually against control fields taken from the GOODS-S region of 3D-HST survey. With this method we determine the fields of 10 targets to be cluster candidates. Four cluster candidates are found by both methods. The two methods disagree in some cases, depending on the specific properties of each field. For the most distant radio galaxy in the 3CR catalog (3C257 at z = 2.47), we identify a population of bluer ETGs that lie on the expected location of the RS model for that redshift. This appears to be the general behavior of ETGs in our fields and it is possibly a signature of the evolution of such galaxies. Our results are consistent with half of the z > 1 radio galaxies being located in dense, rapidly evolving environments.
The effects of many physical processes in the intracluster medium of galaxy clusters imprint themselves in X-ray surface brightness images. It is therefore important to choose optimal methods for extracting information from and enhancing the interpretability of such images. We describe in detail a gradient filtering edge detection method that we previously applied to images of the Centaurus cluster of galaxies. The Gaussian gradient filter measures the gradient in the surface brightness distribution on particular spatial scales. We apply this filter on different scales to Chandra X-ray observatory images of two clusters with AGN feedback, the Perseus cluster and M87, and a merging system, A3667. By combining filtered images on different scales using radial filters spectacular images of the edges in a cluster are produced. We describe how to assess the significance of features in filtered images. We find the gradient filtering technique to have significant advantages for detecting many kinds of features compared to other analysis techniques, such as unsharp-masking. Filtering cluster images in this way in a hard energy band allows shocks to be detected.
We have conducted a feasibility study to determine the effectiveness of using USNO-B1.0 data to preferentially detect objects with displaced red and blue components. A procedure was developed to search catalogue entries for such objects, which include M dwarfs paired with white dwarfs or with earlier main-sequence stars, and galaxies with asymmetric colour distributions. Residual differences between red and blue and infrared and blue scanned emulsion images define vectors, which, when appropriately aligned and of sufficient length, signal potential candidates. Test sample sets were analysed to evaluate the effective discrimination of the technique. Over 91,000 USNO-B1.0 catalogue entries at points throughout the celestial sphere were then filtered for acceptable combinations of entry observations and magnitudes and the resulting total of about 17,000 entries was winnowed down to a little more than 200 objects of interest. These were screened by visual examination of photo images to a final total of 146 candidates. About one quarter of these candidates coincide with SDSS data. Those constituents fall into two groups, single and paired objects. SDSS identified several galaxies in the first group. Regarding the second group, at least half of its members were tentatively identified as main-sequence pairs, the greater portion being of widely separated spectral types. Two white dwarf--main-sequence pairs were also identified. Most importantly, the vectors formed from USNO-B1.0 residuals were in alignment with corresponding SDSS pair position angles, thereby supporting this work`s central thesis.
Carbon, nitrogen, and oxygen (CNO) are key elements in stellar formation and evolution, and their abundances should also have a significant impact on planetary formation and evolution. We present a detailed spectroscopic analysis of 74 solar-type stars, 42 of which are known to harbour planets. We determine the nitrogen abundances of these stars and investigate a possible connection between N and the presence of planetary companions. We used VLT/UVES to obtain high-resolution near-UV spectra of our targets. Spectral synthesis of the NH band at 3360A was performed with the spectral synthesis codes MOOG and FITTING. We identify several spectral windows from which accurate N abundance can be obtained. Nitrogen distributions for stars with and without planets show that planet hosts are nitrogen-rich when compared to single stars. However, given the linear trend between [N/Fe] vs [Fe/H], this fact can be explained as being due to the metal-rich nature of planet hosts. We conclude that reliable N abundances can be derived for metal-rich solar type stars from the near UV molecular band at 3360A. We confirm a linear trend between [N/Fe] and metallicity expected from standard models of Galactic chemical evolution.
Links to: arXiv, form interface, find, astro-ph, recent, 1605, contact, help (Access key information)
The [NII] 122 and 205 \mu m transitions are powerful tracers of the ionized gas in the interstellar medium. By combining data from 21 galaxies selected from the Herschel KINGFISH and Beyond the Peak surveys, we have compiled 141 spatially resolved regions with a typical size of ~1 kiloparsec, with observations of both [NII] far-infrared lines. We measure [NII] 122/205 line ratios in the ~0.6-6 range, which corresponds to electron gas densities $n_e$~1-300 cm$^{-3}$, with a median value of $n_e$=30 cm$^{-3}$. Variations in the electron density within individual galaxies can be as a high as a factor of ~50, frequently with strong radial gradients. We find that $n_e$ increases as a function of infrared color, dust-weighted mean starlight intensity, and star formation rate surface density ($\Sigma_{SFR}$). As the intensity of the [NII] transitions is related to the ionizing photon flux, we investigate their reliability as tracers of the star formation rate (SFR). We derive relations between the [NII] emission and SFR in the low-density limit and in the case of a log-normal distribution of densities. The scatter in the correlation between [NII] surface brightness and $\Sigma_{SFR}$ can be understood as a property of the $n_e$ distribution. For regions with $n_e$ close to or higher than the [NII] line critical densities, the low-density limit [NII]-based SFR calibration systematically underestimates the SFR since [NII] emission is collisionally quenched. Finally, we investigate the relation between [NII] emission, SFR, and $n_e$ by comparing our observations to predictions from the MAPPINGS-III code.
We study the effects of local environment and stellar mass on galaxy properties using a mass complete sample of quiescent and star-forming systems in the COSMOS field at $z\lesssim$ 3. We show that at $z\lesssim$ 1, the median star-formation rate (SFR) and specific SFR (sSFR) of all galaxies depend on environment, but they become independent of environment at $z\gtrsim$ 1. However, we find that only for \textit{star-forming} galaxies, the median SFR and sSFR are similar in different environments, regardless of redshift and stellar mass. We find that the quiescent fraction depends on environment at $z\lesssim$ 1, and on stellar mass out to $z\sim$ 3. We show that at $z\lesssim$ 1, galaxies become quiescent faster in denser environments and that the overall environmental quenching efficiency increases with cosmic time. Environmental and mass quenching processes depend on each other. At $z\lesssim$ 1, denser environments more efficiently quench galaxies with higher masses (log($M/M_{\odot}$)$\gtrsim$ 10.7), possibly due to a higher merger rate of massive galaxies in denser environments, and that mass quenching is more efficient in denser regions. We show that the overall mass quenching efficiency ($\epsilon_{mass}$) for more massive galaxies (log($M/M_{\odot}$)$\gtrsim$ 10.2) rises with cosmic time until $z\sim$ 1 and flattens out since then. However, for less massive galaxies, the rise in $\epsilon_{mass}$ continues to the present time. Our results suggest that environmental quenching is only relevant at $z\lesssim$ 1, likely a fast process, whereas mass quenching is the dominant mechanism at $z\gtrsim$ 1, with a possible stellar feedback physics.
The UV/optical variability of active galactic nuclei and quasars is useful for understanding the physics of the accretion disk and is gradually attributed to the stochastic fluctuations over the accretion disk. Quasars generally appear bluer when they brighten in the UV/optical, the nature of which remains controversial. Recently \citeauthor{Sun2014} discovered that the color variation of quasars is timescale dependent, in the way that faster variations are even bluer than longer term ones. While this discovery can directly rule out models that simply attribute the color variation to contamination from the host galaxies, or to changes in the global accretion rates, it favors the stochastic disk fluctuation model as fluctuations in the innermost hotter disk could dominate the short-term variations. In this work, we show that a revised inhomogeneous disk model, where the characteristic timescales of thermal fluctuations in the disk are radius-dependent (i.e., $\tau \sim r$; based on the one originally proposed by \citeauthor{DexterAgol2011}), can well reproduce a timescale dependent color variation pattern, similar to the observed one and unaffected by the un-even sampling and photometric error. This demonstrates that one may statistically use variation emission at different timescales to spatially resolve the accretion disk in quasars, thus opens a new window to probe and test the accretion disk physics in the era of time domain astronomy. Caveats of the current model, which ought to be addressed in future simulations, are discussed.
Differentiating between active galactic nuclei (AGN) activity and star formation in z ~ 2 galaxies is difficult because traditional methods, such as line ratio diagnostics, change with redshift while multi-wavelength methods (X-ray, radio, IR) are sensitive to only the brightest AGN. We have developed a new method for spatially resolving emission lines in HST/WFC3 G141 grism spectra and quantifying AGN activity through the spatial gradient of the [O III]/H$\beta$ line ratio. Through detailed simulations, we show that our novel line-ratio gradient approach identifies ~ sim 40% more low-mass and obscured AGN than obtained by classical methods. Based on our simulations, we developed a relationship that maps stellar mass, star formation rate, and measured [O III]/H$\beta$ gradient to AGN Eddington ratio. We apply our technique to previously studied stacked samples of galaxies at z ~2 and find that our results are consistent with these studies. Using this gradient method will also be able to inform other galaxy evolution science, such as inside-out quenching and metallicity gradients, and will be widely applicable to future spatially resolved JWST data.
We investigate the shallow increase in globular cluster half-light radii with projected galactocentric distance $R_{gc}$ observed in the giant galaxies M87, NGC 1399, and NGC 5128. To model the trend in each galaxy, we explore the effects of orbital anisotropy and tidally under-filling clusters. While a strong degeneracy exists between the two parameters, we use kinematic studies to help constrain the distance $R_\beta$ beyond which cluster orbits become anisotropic, as well as the distance $R_{f\alpha}$ beyond which clusters are tidally under-filling. For M87 we find $R_\beta > 27$ kpc and $20 < R_{f\alpha} < 40$ kpc and for NGC 1399 $R_\beta > 13$ kpc and $10 < R_{f\alpha} < 30$ kpc. The connection of $R_{f\alpha}$ with each galaxy's mass profile indicates the relationship between size and $R_{gc}$ may be imposed at formation, with only inner clusters being tidally affected. The best fitted models suggest the dynamical histories of brightest cluster galaxies yield similar present-day distributions of cluster properties. For NGC 5128, the central giant in a small galaxy group, we find $R_\beta > 5$ kpc and $R_{f\alpha} > 30$ kpc. While we cannot rule out a dependence on $R_{gc}$, NGC 5128 is well fitted by a tidally filling cluster population with an isotropic distribution of orbits, suggesting it may have formed via an initial fast accretion phase. Perturbations from the surrounding environment may also affect a galaxy's orbital anisotropy profile, as outer clusters in M87 and NGC 1399 have primarily radial orbits while outer NGC 5128 clusters remain isotropic.
We present the discovery of a giant $\gtrsim$100 kpc Ly$\alpha$ nebula detected in the core of the X-ray emitting cluster CL J1449+0856 at $z=1.99$ through Keck/LRIS narrow-band imaging. This detection extends the known relation between Ly$\alpha$ nebulae and overdense regions of the Universe to the dense core of a $5-7\times10^{13}$ M$_{\odot}$ cluster. The most plausible candidates to power the nebula are two Chandra-detected AGN host cluster members. Given the physical conditions of the Ly$\alpha$-emitting gas and the possible interplay with the X-ray phase, we argue that the Ly$\alpha$ nebula would be short-lived ($\lesssim10$ Myr) if not continuously replenished with cold gas at a rate of $\gtrsim1000$ Myr. Cooling from the X-ray phase is disfavored as the replenishing mechanism, primarily because of the high Ly$\alpha$ to X-ray luminosity ratio ($L_{\mathrm{Ly\alpha}}/L_{\mathrm{X}} \approx0.3$), $\gtrsim10-1000\times$ higher than in local cool-core clusters. Cosmological cold flows are disfavored by current modeling. Thus, the cold gas is most plausibly supplied by cluster galaxies through massive outflows. An independent estimate of the total mass outflow rate of core members, based on the observed star formation and black hole accretion rates, matches the required replenishment to sustain the nebula. This scenario directly implies the extraction of energy from galaxies and its deposition in the surrounding intracluster medium, as required to explain the thermodynamic properties of local clusters. We estimate an energy injection of the order of $\thickapprox2$ keV per particle in the intracluster medium over a $2$ Gyr interval. AGN provide $75-85$% of the injected energy and $\approx66$% of the mass, while the rest is supplied by supernovae-driven winds.
We present new rest frame UV and visible observations of 22 high-redshift (1 < z < 2.5) 3C radio galaxies and QSOs obtained with the Hubble Space Telescope's (HST) Wide Field Camera 3 (WFC3) instrument. Using a custom data reduction strategy in order to assure the removal of cosmic rays, persistence signal, and other data artifacts, we have produced high-quality science-ready images of the targets and their local environments. We observe targets with regions of UV emission suggestive of active star formation. In addition, several targets exhibit highly distorted host galaxy morphologies in the rest frame visible images. Photometric analyses reveals that brighter QSOs tend to be generally redder than their dimmer counterparts. Using emission line fluxes from the literature, we estimate that emission line contamination is relatively small in the rest frame UV images for the QSOs. Using archival VLA data, we have also created radio map overlays for each of our targets, allowing for analysis of the optical and radio axes alignment.
We show, for the first time, that ${\rm H_2}$ formation on dust grains can be enhanced in disk galaxies under strong ram-pressure (RP). We numerically investigate how the time evolution, of ${\rm H}$ {\sc i} and ${\rm H_2}$ components in disk galaxies orbiting a group/cluster of galaxies, can be influenced by hydrodynamical interaction between the gaseous components of the galaxies and the hot intra-cluster medium (ICM). We find that compression of ${\rm H}$ {\sc i} caused by RP increases ${\rm H_2}$ formation in disk galaxies, before RP rapidly strips ${\rm H}$ {\sc i}, cutting off the fuel supply and causing a drop in ${\rm H_2}$ density. We also find that the level of this ${\rm H_2}$ formation enhancement in a disk galaxy under RP depends on the mass of its host cluster dark matter (DM) halo, initial positions and velocities of the disk galaxy, and disk inclination angle with respect to the orbital plane. We demonstrate that dust growth is a key factor in the evolution of the ${\rm H}$ {\sc i} and ${\rm H_2}$ mass in disk galaxies under strong RP. We discuss how the correlation between ${\rm H_2}$ fractions and surface gas densities of disk galaxies evolves with time in the galaxies under RP. We also discuss whether or not galaxy-wide star formation rates (SFRs) in cluster disk galaxies can be enhanced by RP if the SFRs depend on ${\rm H_2}$ densities.
As part of the Z-PAndAS Keck II DEIMOS survey of resolved stars in our neighboring galaxy, Andromeda (M31), we have built up a unique data set of measured velocities and chemistries for thousands of stars in the Andromeda stellar halo, particularly probing its rich and complex substructure. In this contribution, we will discuss the structural, dynamical and chemical properties of Andromeda's dwarf spheroidal galaxies, and how there is no observational evidence for a difference in the evolutionary histories of those found on and off M31's vast plane of satellites. We will also discuss a possible extension to the most significant merger event in M31 - the Giant Southern Stream - and how we can use this feature to refine our understanding of M31's mass profile, and its complex evolution.
The secular thickening of a self-gravitating stellar galactic disc embedded in a fluctuating potential is investigated. The thick WKB limit for the diffusion coefficient of the corresponding dressed Fokker-Planck equation is found using the epicyclic approximation, while assuming that only radially tightly wound transient spirals are sustained by the disc. This yields a simple quadrature, providing a clear understanding of the positions of maximum vertical orbital diffusion within the disc. This thick limit also offers a consistent derivation of a thick disc Toomre parameter, which is shown to be exponentially boosted by the ratio of the vertical to radial scale heights. When applied to a tepid stable tapered disc perturbed by shot noise, this formalism predicts the formation of ridges of resonant orbits towards larger vertical actions, as found in simulations. Potential fluctuations within the disc statistically induce a vertical bending of a subset of resonant orbits, triggering the corresponding increase in vertical velocity dispersion. A simple ad hoc model of fading sequences of slowing bars explains how such a diffusion process also provides a possible mechanism allowing for the resonant thickening of galactic discs over secular timescales: the faster the fading the broader the effect.
The secular thickening of a discrete self-gravitating galactic disc is investigated using the inhomogeneous multi-component Balescu-Lenard equation. The thick WKB limit for the diffusion and drift coefficients is found using the epicyclic approximation, while assuming that only radially tightly wound transient spirals are sustained by the disc. This yields a simple double quadrature for the drift and diffusion coefficients, providing a clear understanding of the positions of maximum vertical orbital diffusion within the disc induced by the effects of a finite number of particles. When applied to a tepid stable tapered disc, the Balescu-Lenard formalism predicts the formation of ridges of resonant orbits towards larger vertical actions, as found in direct numerical simulations, but over-estimates the timescale involved in their appearance. Swing amplication is likely needed to resolve this discrepancy, as demonstrated in the case of razor-thin discs. The joint evolution of a population of giant molecular clouds within the disc may accelerate the secular disc's thickening induced by finite${-N}$ effects, but the observed number of clouds in the Milky Way does not seem to be sufficient to explain its thick disc.
By cross-matching the currently largest optical catalog of galaxy clusters and the NVSS radio survey database, we obtain the largest complete sample of brightest cluster galaxies (BCGs) in the redshift range of 0.05<z<0.45, which have radio emission and redshift information. We confirm that more powerful radio BCGs tend to be these optically very bight galaxies located in more relaxed clusters. We derived the radio luminosity functions of BCGs from the largest complete sample of BCGs, and find that the functions depend on the optical luminosity of BCGs and the dynamical state of galaxy clusters. However, the radio luminosity function does not show significant evolution with redshift.
The physical properties of galactic winds are one of the keys to understand galaxy formation and evolution. These properties can be constrained thanks to background quasar lines of sight (LOS) passing near star-forming galaxies (SFGs). We present the first results of the MusE GAs FLOw and Wind (MEGAFLOW) survey obtained of 2 quasar fields which have 8 MgII absorbers of which 3 have rest-equivalent width greater than 0.8 \AA. With the new Multi Unit Spectroscopic Explorer (MUSE) spectrograph on the Very Large Telescope (VLT), we detect 6 (75$\%$) MgII host galaxy candidates withing a radius of 30 arcsec from the quasar LOS. Out of these 6 galaxy--quasar pairs, from geometrical arguments, one is likely probing galactic outflows, two are classified as "ambiguous", two are likely probing extended gaseous disks and one pair seems to be a merger. We focus on the wind$-$pair and constrain the outflow using a high resolution quasar spectra from Ultraviolet and Visual Echelle Spectrograph (UVES). Assuming the metal absorption to be due to gas flowing out of the detected galaxy through a cone along the minor axis, we find outflow velocities of the order of $\approx$ 150 km/s (i.e. smaller than the escape velocity) with a loading factor, $\eta =\dot M_{\rm out}/$SFR, of $\approx$ 0.7. We see evidence for an open conical flow, with a low-density inner core. In the future, MUSE will provide us with about 80 multiple galaxy$-$quasar pairs in two dozen fields.
We present N/O abundance ratios, ionization parameters $q_{\rm ion}$, and oxygen abundance O/H of eleven $z\sim 2$ galaxies determined by the direct temperature $T_{\rm e}$ method with [OIII] $\lambda$4363 and OIII] $\lambda$1665 lines, and investigate galaxy evolution from $z\sim 0$ to $2$ in conjunction with $T_{\rm e}$-method measurements of 208,529 SDSS galaxies and 9 green pea galaxies (GPs). We identify that three out of our eleven $z\sim 2$ galaxies clearly fall beyond the local average of N/O-O/H and N/O-stellar mass ($M_{\star}$) relations, while the rest of the $z\sim 2$ galaxies have N/O ratios comparable with $z\sim 0$ galaxies. The eleven $z\sim 2$ galaxies place the upper limit of N/O ratio $\log ({\rm N/O})\le -1.28$ on average, suggesting that the N/O ratio evolves, if any, by $<0.15$ dex. We find that two of our $z\sim 2$ galaxies with the significant BPT offsets show both N/O ratio and $q_{\rm ion}$ about 0.4-dex higher than the local average, while there exist $z\sim 0$ green-pea galaxies whose N/O ratio and $q_{\rm ion}$ are comparable with these $z\sim 2$ galaxies. Although three (none) of our $z\sim 2$ galaxies have an excess only in $q_{\rm ion}$ (N/O ratio), there are no (many) high-$z$ counterparts so far found in the SDSS and GP galaxies. These missing $z\sim 2$ and local galaxies with the single $q_{\rm ion}$ or N/O-ratio excess may be found after $T_{\rm e}$-method studies proceed for $z\sim 2$ galaxies and rare GP-like objects. Because a large fraction of our $z\sim 2$ galaxies show $q_{\rm ion}$ and/or N/O-ratio excess(es) unlike the average local galaxies, we infer that galaxies with the BPT offsets emerge at $z\sim 2$ and mostly disappear in the present-day universe, except for galaxies such as GPs.
We present results of a photometric study into the cluster population of NGC 1566, a nearby grand design spiral galaxy, sampled out to a Galactocentric radius of $\approx 5.5$ kpc. The shape of the mass-limited age distribution shows negligible variation with radial distance from the centre of the galaxy, and demonstrates three separate sections, with a steep beginning, flat middle and steep end. The luminosity function can be approximated by a power law at lower luminosities with evidence of a truncation at higher luminosity. The power law section of the luminosity function of the galaxy is best fitted by an index $\approx -2$, in agreement with other studies, and is found to agree with a model luminosity function, which uses an underlying Schechter mass function. The recovered power law slope of the mass distribution shows a slight steepening as a function of galactocentric distance, but this is within error estimates. It also displays a possible truncation at the high mass end. Additionally, the cluster formation efficiency ($\Gamma$) and the specific U-band luminosity of clusters ($T_L(U)$) are calculated for NGC 1566 and are consistent with values for similar galaxies. A difference in NGC 1566, however, is that the fairly high star formation rate is in contrast with a low $\Sigma_{SFR}$ and $\Gamma$, indicating that $\Gamma$ can only be said to depend strongly on $\Sigma_{SFR}$, not the star formation rate.
L1551 is chosen because it is relatively isolated in the Taurus molecular cloud shielded from FUV photons, providing an ideal environment for studying the target properties. Our observations cover ~40'x40' with resolution ~30", which are the maps with highest spatial dynamical range to date. We derive the X(13CO)/X(C18O) value on the sub-parsec scales in the range of ~3-27 with a mean value of 8.0+-2.8. Comparing to the visual extinction map derived from the Herschel observations, we found that the abundance ratio reaches its maximum at low Av (i.e., Av ~ 1-4mag), and decreases to the typical solar system value of 5.5 inside L1551 MC. The high X(13CO)/X(C18O) value at the boundary of the cloud is most likely due to the selective FUV photodissociation of C18O. This is in contrast with Orion-A where its internal OB stars keep the abundance ratio at a high level greater than ~10. In addition, we explore the variation of the X-factor, because it is an uncertain but widely used quantity in extragalactic studies. We found that X-factor is proportional to N(H2)^1.0 which is consistent with previous simulations. Excluding the high density region, the average X-factor is similar to the Milky Way average value.
This paper aims at providing aperture corrections for emission lines in a
sample of spiral galaxies from the Calar Alto Legacy Integral Field Area Survey
(CALIFA) database. In particular, we explore the behavior of the
log([OIII]5007/Hbeta)/([NII]6583/Halpha) (O3N2) and log[NII]6583/Halpha (N2)
flux ratios since they are closely connected to different empirical
calibrations of the oxygen abundances in star forming galaxies.
We compute median growth curves of Halpha, Halpha/Hbeta, O3N2 and N2 up to
2.5R_50 and 1.5 disk R_eff. The growth curves simulate the effect of observing
galaxies through apertures of varying radii. The median growth curve of the
Halpha/Hbeta ratio monotonically decreases from the center towards larger
radii, showing for small apertures a maximum value of ~10% larger than the
integrated one. The median growth curve of N2 shows a similar behavior,
decreasing from the center towards larger radii. No strong dependence is seen
with the inclination, morphological type and stellar mass for these growth
curves. Finally, the median growth curve of O3N2 increases monotonically with
radius. However, at small radii it shows systematically higher values for
galaxies of earlier morphological types and for high stellar mass galaxies.
Applying our aperture corrections to a sample of galaxies from the SDSS
survey at 0.02<=z<=0.3 shows that the average difference between fiber-based
and aperture corrected oxygen abundances, for different galaxy stellar mass and
redshift ranges, reaches typically to ~11%, depending on the abundance
calibration used. This average difference is found to be systematically biased,
though still within the typical uncertainties of oxygen abundances derived from
empirical calibrations. Caution must be exercised when using observations of
galaxies for small radii (e.g. below 0.5R_eff) given the high dispersion shown
around the median growth curves.
Radio maps of AGNs often show linear features, called jets, both on VLBI as well as kpc scales. These jets supposedly possess relativistic motion and are oriented close to the line of sight of the observer and accordingly the relativistic Doppler beaming makes them look much brighter than they really are in their respective rest-frames. The flux boosting due to the relativistic beaming is a very sensitive factor of the jet orientation angle, as seen by the observer. Quite often large bends are seen in these jets, with misalignments being $90^\circ$ or more and might imply a change in the orientation angle that could cause a large change in the relativistic beaming factor. But these large bends do not show high contrasts in the brightness of the jets, before and after the misalignments, which seem discordant with the relativistic beaming models. Is something amiss here? It needs to be kept in mind that sometimes a small intrinsic change in the jet angle might appear as a much larger misalignment due to the geometrical projection effects, especially when seen close to the line of sight. Of course what really matters is the final orientation angle of the jet with respect to the observer's line of sight. After taking the geometrical projection effects into account, we show that the assumed relativistic beaming does not seem conversant with the absence of large contrast in the jet brightness before and after the observed misalignments.
We investigate the process of reionisation in a model in which the dark matter is a warm elementary particle such as a sterile neutrino. We focus on models that are consistent with the dark matter decay interpretation of the recently detected line at 3.5 keV in the X-ray spectra of galaxies and clusters. In warm dark matter models the primordial spectrum of density perturbations has a cut-off on the scale of dwarf galaxies. Structure formation therefore begins later than in the standard cold dark matter (CDM) model and very few objects form below the cut-off mass scale. To calculate the number of ionising photons, we use the Durham semi-analytic model of galaxy formation, GALFORM. We find that even the most extreme 7 keV sterile neutrino we consider is able to reionise the Universe early enough to be compatible with the bounds on the epoch of reionisation from Planck. This, perhaps surprising, result arises from the rapid build-up of high redshift galaxies in the sterile neutrino models which is also reflected in a faster evolution of their far-UV luminosity function between $10>z>7$ than in CDM. The dominant sources of ionising photons are systematically more massive in the sterile neutrino models than in CDM. As a consistency check on the models, we calculate the present-day luminosity function of satellites of Milky Way-like galaxies. When the satellites recently discovered in the DES survey are taken into account, strong constraints are placed on viable sterile neutrino models.
We present a detailed X-ray variability study of the low mass Active Galactic Nuclei (AGN) NGC 7314 using the two newly obtained XMM-Newton observations ($140$ and $130$ ks), together with two archival data sets of shorter duration ($45$ and $84$ ks). The relationship between the X-ray variability characteristics and other physical source properties (such as the black hole mass) are still relatively poorly defined, especially for low-mass AGN. We perform a new, fully analytical, power spectral density (PSD) model analysis method, which will be described in detail in a forthcoming paper, that takes into consideration the spectral distortions, caused by red-noise leak. We find that the PSD in the $0.5-10$ keV energy range, can be represented by a bending power-law with a bend around $6.7\times10^{-5}$ Hz, having a slope of $0.51$ and $1.99$ below and above the bend, respectively. Adding our bend time-scale estimate, to an already published ensemble of estimates from several AGN, supports the idea that the bend time-scale depends linearly only on the black hole mass and not on the bolometric luminosity. Moreover, we find that as the energy range increases, the PSD normalization increases and there is a hint that simultaneously the high frequency slope becomes steeper. Finally, the X-ray time-lag spectrum of NGC 7314 shows some very weak signatures of relativistic reflection, and the energy resolved time-lag spectrum, for frequencies around $3\times10^{-4}$ Hz, shows no signatures of X-ray reverberation. We show that the previous claim about ks time-delays in this source, is simply an artefact induced by the minuscule number of points entering during the time-lag estimation in the low frequency part of the time-lag spectrum (i.e. below $10^{-4}$ Hz).
We present a Fourier space method to combine the publicly available Herschel PACS and SPIRE data with the Planck thermal dust emission model. The method effectively combines the Planck large angular scale emission with the small scale \herschel emission, similar to the feathering method used in interferometry and recently implemented by Csengeri et al. to correct the ground-based ATLASGAL data. This method eliminates the pervasive negative fluxes present in the PACS 160 micron archive data while preserving the structure in the background of both the PACS and SPIRE data. We generate column density and temperature maps from data calibrated with this method. We analyze the validity of our new method using star-forming regions in a range of environments: low-mass Bok Globules (B68), nearby star-forming regions (Perseus), and high-mass star forming regions in the Galactic plane (including IRDC G11 and W31). We accurately recover low column density material, comparing well to previous near-infrared extinction methods at intermediate column densities, and providing a more accurate temperature assumption for interpretation of single frequency continuum data such as those obtained by ATLASGAL. The Herschel data, combined with our method, allow studies of the column density and temperature distributions of molecular clouds across a wide range of galactic environments: a key step to study the physics of star and cluster formation.
We present the abundance analysis for a sample of 18 red giant branch stars in the metal-poor globular cluster NGC 4147 based on medium and high resolution spectra. This is the first extensive spectroscopic study of this cluster. We derive abundances of C, N, O, Na, Mg, Al, Si, Ca, Ti, Cr, Fe, Ni, Y, Ba, and Eu. We find a metallicity of [Fe/H]=-1.84+-0.02 and an alpha-enhancement of +0.38+-0.05 (errors on the mean), typical of halo globular clusters in this metallicity regime. A significant spread is observed in the abundances of light elements C, N, O, Na, and Al. In particular we found a Na-O anti-correlation and Na-Al correlation. The cluster contains only 15% of stars that belong to the first generation (Na-poor and O-rich). This implies that it suffered a severe mass loss during its lifetime. Its [Ca/Fe] and [Ti/Fe] mean values agree better with the Galactic Halo trend than with the trend of extragalactic environments at the cluster metallicity. This possibly suggests that NGC 4147 is a genuine Galactic object at odd with what claimed by some author that proposed the cluster to be member of the Sagittarius dwarf galaxy. A anti-relation between the light s-process element Y and Na may also be present.
AX J1845-0258 is a transient X-ray pulsar, with spin period of 6.97s, discovered with the ASCA satellite in 1993. Its soft spectrum and the possible association with a supernova remnant suggest that AX J1845-0258 might be a magnetar, but this has not been confirmed yet. A possible counterpart one order of magnitude fainter, AX J184453-025640, has been found in later X-ray observations, but no pulsations have been detected. In addition, some other X-ray sources are compatible with the pulsar location, which is in a crowded region of the Galactic plane. We have carried out a new investigation of all the X-ray sources in the ASCA error region of AX J1845-0258, using archival data obtained with Chandra in 2007 and 2010, and with XMM-Newton in 2010. We set an upper limit of 6% on the pulsed fraction of AX J184453-025640 and confirmed its rather hard spectrum (power law photon index of 1.2 +\- 0.3). In addition to the other two fainter sources already reported in the literature, we found other X-ray sources positionally consistent with AX J1845-0258. Although many of them are possibly foreground stars likely unrelated to the pulsar, at least another new source, CXOU J184457.5-025823, could be a plausible counterpart of AX J1845-0258. It has a flux of 6x10^{-14} erg cm^{-2} s^{-1} and a spectrum well fit by a power law with photon index ~1.3 and Nh ~ 10^{22} cm^{-2}.
Links to: arXiv, form interface, find, astro-ph, recent, 1605, contact, help (Access key information)
We describe a high-order numerical magnetohydrodynamics (MHD) solver built upon a novel non-linear entropy stable numerical flux function that supports eight travelling wave solutions. By construction the solver conserves mass, momentum, and energy and is entropy stable. The method is designed to treat the divergence-free constraint on the magnetic field in a similar fashion to a hyperbolic divergence cleaning technique. The solver described herein is especially well-suited for flows involving strong discontinuities. Furthermore, we present a new formulation to guarantee positivity of the pressure. We present the underlying theory and implementation of the new solver into the multi-physics, multi-scale adaptive mesh refinement (AMR) simulation code $\texttt{FLASH}$ (this http URL). The accuracy, robustness and computational efficiency is demonstrated with a number of tests, including comparisons to available MHD implementations in $\texttt{FLASH}$.
The hierarchical nature of LCDM suggests that the Magellanic Clouds must have been surrounded by a number of satellites before their infall into the Milky Way. Many of those satellites should still be in close proximity to the Clouds, but some could have dispersed ahead/behind the Clouds along their Galactic orbit. Either way, prior association with the Clouds results in strong restrictions on the present-day positions and velocities of candidate Magellanic satellites: they must lie close to the nearly-polar orbital plane of the Magellanic stream, and their distances and radial velocities must follow the latitude dependence expected for a tidal stream with the Clouds at pericenter. We use a cosmological numerical simulation of the disruption of a massive subhalo in a Milky Way-sized LCDM halo to test whether any of the 20 dwarfs recently-discovered in the DES, SMASH, Pan-STARRS, and ATLAS surveys are truly associated with the Clouds. Of the 6 systems with kinematic data, only Hydra II and Hor 1 have distances and radial velocities consistent with a Magellanic origin. Of the remaining dwarfs, six (Hor 2, Eri 3, Ret 3, Tuc 4, Tuc 5, and Phx 2) have positions and distances consistent with a Magellanic origin, but kinematic data are needed to substantiate that possibility. Conclusive evidence for association would require proper motions to constrain the orbital angular momentum direction, which, for true Magellanic satellites, must coincide with that of the Clouds. We use this result to predict radial velocities and proper motions for all new dwarfs. Our results are relatively insensitive to the assumption of first or second pericenter for the Clouds.
We use numerical simulations to examine the effects of radial migration on the vertical structure of galaxy disks. The simulations follow three exponential disks of different mass but similar circular velocity, radial scalelength, and (constant) scale height. The disks develop different non-axisymmetric patterns, ranging from feeble, long-lived multiple arms to strong, rapidly-evolving few-armed spirals. These fluctuations induce radial migration through secular changes in the angular momentum of disk particles, mixing the disk radially and blurring pre-existing gradients. Migration affects primarily stars with small vertical excursions, regardless of spiral pattern. This "provenance bias" largely determines the vertical structure of migrating stars: inward migrators thin down as they move in, whereas outward migrators do not thicken up but rather preserve the disk scale height at destination. Migrators of equal birth radius thus develop a strong scale-height gradient, not by flaring out as commonly assumed, but by thinning down as they spread inward. Similar gradients have been observed for low-[$\alpha$/Fe] mono-abundance populations (MAPs) in the Galaxy but our results argue against interpreting them as a consequence of radial migration. This is because outward migration does not lead to thickening, implying that the maximum scale height of any population should reflect its value at birth. In contrast, Galactic MAPs have scale heights that increase monotonically outwards, reaching values that greatly exceed those at their presumed birth radii. Given the strong vertical bias affecting migration, a proper assessment of the importance of radial migration in the Galaxy should take carefully into account the strong radial dependence of the scale heights of the various stellar populations.
We discuss the potential of the gravitational microlensing method as a unique tool to detect unambiguous signals caused by intermediate-mass black holes in globular clusters. We select clusters near the line of sight to the Galactic Bulge and the Small Magellanic Cloud, estimate the density of background stars for each of them, and carry out simulations in order to estimate the probabilities of detecting the astrometric signatures caused by black hole lensing. We find that for several clusters, the probability of detecting such an event is significant with available archival data from the Hubble Space Telescope. Specifically, we find that M 22 is the cluster with the best chances of yielding an IMBH detection via astrometric microlensing. If M 22 hosts an IMBH of mass $10^5M_\odot$, then the probability that at least one star will yield a detectable signal over an observational baseline of 20 years is $\sim 86\%$, while the probability of a null result is around $14\%$. For an IMBH of mass $10^6M_\odot$, the detection probability rises to $>99\%$. Future observing facilities will also extend the available time baseline, improving the chance of detections for the clusters we consider.
We use $\sim$83,000 star-forming galaxies at $0.04<z<0.3$ from the Sloan Digital Sky Survey to study the so-called fundamental metallicity relation (FMR) and report on the disappearance of its anti-correlation between metallicity and star formation rate (SFR) when using the new metallicity indicator recently proposed by Dopita et al. In this calibration, metallicity is primarily sensitive to the emission line ratio [NII]$\lambda$6584 / [SII]$\lambda\lambda$6717, 6731 that is insensitive to dilution by pristine, infalling gas that may drive the FMR anti-correlation with SFR. Therefore, we conclude that the apparent disappearance of the FMR (using this new metallicity indicator) does not rule out its existence.
We perform adaptive mesh refinement (AMR) and smoothed particle hydrodynamics (SPH) cosmological zoom simulations of a region around a forming galaxy cluster, comparing the ability of the methods to handle successively more complex baryonic physics. In the simplest, non-radiative case, the two methods are in good agreement with each other, but the SPH simulations generate central cores with slightly lower entropies and virial shocks at slightly larger radii, consistent with what has been seen in previous studies. The inclusion of radiative cooling, star formation, and stellar feedback leads to much larger differences between the two methods. Most dramatically, at z=5, rapid cooling in the AMR case moves the accretion shock well within the virial radius, while this shock remains near the virial radius in the SPH case, due to excess heating, coupled with poorer capturing of the shock width. On the other hand, the addition of feedback from active galactic nuclei (AGN) to the simulations results in much better agreement between the methods. In this case both simulations display halo gas entropies of 100 keV cm^2, similar decrements in the star-formation rate, and a drop in the halo baryon content of roughly 30%. This is consistent with AGN growth being self-regulated, regardless of the numerical method. However, the simulations with AGN feedback continue to differ in aspects that are not self-regulated, such that in SPH a larger volume of gas is impacted by feedback, and the cluster still has a lower entropy central core.
Several `giant' Lyman-$\alpha$ (Ly$\alpha$) nebulae with extent $\gtrsim 300\,$kpc and observed Ly$\alpha$ luminosity of $\gtrsim 10^{44}\,{\rm erg}\,{\rm s}^{-1}\,{\rm cm}^{-2}\,{\rm arcsec}^{-2}$ have recently been detected, and it has been speculated that their presence hints at a substantial cold gas reservoir in small cool clumps not resolved in modern hydro-dynamical simulations. We use the Illustris simulation to predict the Ly$\alpha$ emission emerging from large halos ($M > 10^{11.5}M_{\odot}$) at $z\sim 2$ and thus test this model. We consider both AGN and star driven ionization, and compared the simulated surface brightness maps, profiles and Ly$\alpha$ spectra to a model where most gas is clumped below the simulation resolution scale. We find that while the cold clumps boost the Ly$\alpha$ luminosity especially in the outer regions of the halo -- as expected by previous work -- with Illustris no additional clumping is necessary to explain the extents and luminosities of the `giant Ly$\alpha$ nebulae' observed. Furthermore, the maximal extents of the objects show a wide spread for a given luminosity and do not correlate significantly with any halo properties. We also show how the detected size depends strongly on the employed surface brightness cutoff, and predict that further such objects will be found in the near future.
Aims: We aim at investigating the time-behaviour of a sample of gamma-ray blazars. We present the results from a 13 month-long optical photometry monitoring campaign of the blazars PKS 0048-097, PKS 0754+100, HB89 0827+243, PKS 0851+202, PKS 1253-055, PKS1510-089, PKS 1749+096, PKS 2230+114 and PKS 2251+158. Methods: We analyse the variability of each object, focusing on different time-scales (long term, short term, and microvariability), in an attempt to achieve a statistical comparison of the results. Results: After applying a geometric model to explain the variability results, we found that it is possible that a slight change in the direction of the jet generates the variations detected in some objects during this campaign.
Due to its youth, proximity and richness the Orion Nebula Cloud (ONC) is an ideal testbed to obtain a comprehensive view on the Initial Mass Function (IMF) down to the planetary mass regime. Using the HAWK-I camera at the VLT, we have obtained an unprecedented deep and wide near-infrared JHK mosaic of the ONC (90% completeness at K~19.0mag, 22'x28). Applying the most recent isochrones and accounting for the contamination of background stars and galaxies, we find that ONC's IMF is bimodal with distinct peaks at about 0.25 and 0.025 M_sun separated by a pronounced dip at the hydrogen burning limit (0.08 M_sun), with a depth of about a factor 2-3 below the log-normal distribution. Apart from ~920 low-mass stars (M < 1.4 M_sun) the IMF contains ~760 brown dwarf (BD) candidates and ~160 isolated planetary mass object (IPMO) candidates with M > 0.005 M_sun, hence about ten times more substellar candidates than known before. The substellar IMF peak at 0.025 M_sun could be caused by BDs and IPMOs which have been ejected from multiple systems during the early star-formation process or from circumstellar disks.
We present a novel approach to draw the synthetic color-magnitude diagram of galaxies, which can provide - in principle - a deeper insight in the interpretation and understanding of current observations. In particular, we `light up' the stars of chemical evolution models, according to their initial mass, metallicity and age, to eventually understand how the assumed underlying galaxy formation and evolution scenario affects the final configuration of the synthetic CMD. In this way, we obtain a new set of observational constraints for chemical evolution models beyond the usual photospheric chemical abundances. The strength of our method resides in the very fine grid of metallicities and ages of the assumed database of stellar isochrones. In this work, we apply our photo-chemical model to reproduce the observed CMD of the Sculptor dSph and find that we can reproduce the main features of the observed CMD. The main discrepancies are found at fainter magnitudes in the main sequence turn-off and sub-giant branch, where the observed CMD extends towards bluer colors than the synthetic one; we suggest that this is a signature of metal-poor stellar populations in the data, which cannot be captured by our assumed one-zone chemical evolution model.
I show that a recently discovered star cluster near the center of the ultra-faint dwarf galaxy Eridanus II provides strong constraints on massive compact halo objects (MACHOs) of >~5 M_sun as the main component of dark matter. MACHO dark matter will dynamically heat the cluster, driving it to larger sizes and higher velocity dispersions until it dissolves into its host galaxy. The star cluster has a luminosity of just ~2000 L_sun and is relatively puffy, with a half-light radius of 13 pc, making it much more fragile than other known clusters in dwarf galaxies. For a wide range of plausible dark matter halo properties, Eri II's star cluster combines with existing constraints from microlensing, wide binaries, and disk kinematics to rule out dark matter composed entirely of MACHOs from ~10$^{-7}$ M_sun up to arbitrarily high masses. The cluster in Eri II closes the ~20--100 M_sun window of allowed MACHO dark matter and provides much stronger constraints than wide Galactic binaries for MACHOs of up to thousands of Solar masses.
We used the Revised Flat Galaxy Catalog (RFGC) to select 817 ultra-flat (UF) edge-on disk galaxies with blue and red apparent axial ratios of $(a/b)_B > 10.0$ and $(a/b)_R > 8.5$. The sample covering the whole sky, except the Milky Way zone, contains 490 UF galaxies with measured radial velocities. Our inspection of the neighboring galaxies around them revealed only 30 companions with radial velocity difference of $\mid\Delta V\mid<500$ km s$^{-1}$ inside the projected separation of $R_p < 250$ kpc. Wherein, the wider area around the UF galaxy within $R_p < 750$ kpc contains no other neighbors brighter than the UF galaxy itself in the same velocity span. The resulting sample galaxies mostly belong to the morphological types Sc, Scd, Sd. They have a moderate rotation velocity curve amplitude of about $120$ km s$^{-1}$ and a moderate K-band luminosity of about $10^{10}L_{\odot}$. The median difference of radial velocities of their companions is $87$ km s$^{-1}$, yielding the median orbital mass estimate of about $5\times10^{11}M_{\odot}$. Excluding six probable non-isolated pairs, we obtained a typical halo-mass-to-stellar-mass of UF galaxies of about $30$, what is almost the same one as in the principal spiral galaxies, like M 31 and M 81 in the nearest groups. We also note that ultra-flat galaxies look two times less "dusty" than other spirals of the same luminosity.
New calculations for rotational excitation of cyanoacetylene by collisions with hydrogen molecules are performed to include the lowest 38 rotational levels of HC3N and kinetic temperatures to 300 K. Calculations are based on the interaction potential of Wernli et al. A&A, 464, 1147 (2007) whose accuracy is checked against spectroscopic measurements of the HC3N-H2 complex. The quantum coupled-channel approach is employed and complemented by quasi-classical trajectory calculations. Rate coefficients for ortho-H2 are provided for the first time. Hyperfine resolved rate coefficients are also deduced. Collisional propensity rules are discussed and comparisons between quantum and classical rate coefficients are presented. This collisional data should prove useful in interpreting HC3N observations in the cold and warm ISM, as well as in protoplanetary disks.
We study the nearby lenticular galaxy NGC 3998. This galaxy is known to host a low-power radio AGN with a kpc-size one-sided jet and a large, nearly polar HI disc. It is therefore a good system to study to understand the relation between the availability of cold-gas and the triggering of AGNs in galaxies. Our new WSRT data reveal two faint, S-shaped radio lobes extending out to $\sim$10 kpc from the galaxy centre. Remarkably, we find that the inner HI disc warps back towards the stellar mid-plane in a way that mirrors the warping of the radio lobes. We suggest that the polar HI disc was accreted through a minor merger, and that the torques causing it to warp in the inner regions are also responsible for feeding the AGN. The "S" shape of the radio lobes would then be due to the radio jets adapting to the changing angular momentum of the accreted gas. The extended radio jets are likely poorly collimated, which would explain their quick fading and, therefore, their rarity in galaxies similar to NGC 3998. The fuelling of the central super-massive black hole is likely occurring via "discrete events", suggested by the observed variability of the radio core and the extremely high core dominance, which we attribute to the formation and ejection of a new jet resulting from a recent fuelling event.
Images of dust continuum and CO line emission are powerful tools for deducing structural characteristics of galaxies, such as disk sizes, H$_2$ gas velocity fields and enclosed H$_2$ and dynamical masses. We report on a fundamental constraint set by the cosmic microwave background (CMB) on the observed structural and dynamical characteristics of galaxies, as deduced from dust continuum and CO-line imaging at high redshifts. As the CMB temperature rises in the distant Universe, the ensuing thermal equilibrium between the CMB and the cold dust and H$_2$ gas progressively erases all spatial and spectral contrasts between their brightness distributions and the CMB. For high-redshift galaxies, this strongly biases the recoverable H$_2$ gas and dust mass distributions, scale lengths, gas velocity fields and dynamical mass estimates. This limitation is unique to mm/submm wavelengths and unlike its known effect on the global dust continuum and molecular line emission of galaxies, it cannot be addressed simply. We nevertheless identify a unique signature of CMB-affected continuum brightness distributions, namely an increasing rather than diminishing contrast between such brightness distributions and the CMB when the cold dust in distant galaxies is imaged at frequencies beyond the Raleigh-Jeans limit. For the molecular gas tracers, the same effect makes the atomic carbon (CI) lines maintain a larger contrast than the CO lines against the CMB.
We present a comparison of molecular clouds (MCs) from a simulation of supernova-driven interstellar medium (ISM) turbulence with real MCs from the Outer Galaxy Survey. The radiative transfer calculations to compute synthetic CO spectra are carried out assuming the CO relative abundance depends only on gas density, according to four different models. Synthetic MCs are selected above a threshold brightness temperature value, $T_{\rm B,min}=1.4$ K, of the $J=1-0$ $^{12}$CO line, generating 16 synthetic catalogs (four different spatial resolutions and four CO abundance models), each containing up to several thousands MCs. The comparison with the observations focuses on the mass and size distributions and on the velocity-size and mass-size Larson relations. The mass and size distributions are found to be consistent with the observations, with no significant variations with spatial resolution or chemical model, except in the case of the unrealistic model with constant CO abundance. The velocity-size relation is slightly too steep for some of the models, while the mass-size relation is a bit too shallow for all models only at a spatial resolution $dx\approx 1$ pc. The normalizations of the Larson relations show a clear dependence on spatial resolution, for both the synthetic and the real MCs. The comparison of the velocity-size normalization suggests that the SN rate in the Perseus arm is approximately 70\% or less of the rate adopted in the simulation. Overall, the realistic properties of the synthetic clouds confirm that supernova-driven turbulence can explain the origin and dynamics of MCs.
We use 3035 Herschel-SPIRE 500$\mu$m sources from 20.3 sq deg of sky in the HerMES Lockman, ES1 and XMM-LSS areas to estimate the star-formation rate density at z = 1-6. 500 mu sources are associated first with 350 and 250 mu sources, and then with Spitzer 24 mu sources from the SWIRE photometric redshift catalogue. The infrared and submillimetre data are fitted with a set of radiative-transfer templates corresponding to cirrus (quiescent) and starburst galaxies. Lensing candidates are removed via a set of colour-colour and colour-redshift constraints. Star-formation rates are found to extend from < 1 to 20,000 Mo/yr. Such high values were also seen in the all-sky IRAS Faint Source Survey. Star-formation rate functions are derived in a series of redshift bins from 0-6, combined with earlier far-infrared estimates, where available, and fitted with a Saunders et al (1990) functional form. The star-formation-rate density as a function of redshift is derived and compared with other estimates. There is reasonable agreement with both infrared and ultraviolet estimates for z < 3, but we find higher star-formation-rate densities than ultraviolet estimates at z = 3-6. Given the considerable uncertainties in the submillimetre estimates, we can not rule out the possibility that the ultraviolet estimates are correct. But the possibility that the ultraviolet estimates have seriously underestimated the contribution of dust-shrouded star-formation can also not be excluded.
This work presents the publicly available moving-mesh magnetohydrodynamics code DISCO. DISCO is efficient and accurate at evolving orbital fluid motion in two and three dimensions, especially at high Mach number. DISCO employs a moving-mesh approach utilizing a dynamic cylindrical mesh that can shear azimuthally to follow the orbital motion of the gas. The moving mesh removes diffusive advection errors and allows for longer timesteps than a static grid. Magnetohydrodynamics is implemented in DISCO using an HLLD Riemann solver and a novel constrained transport scheme which is compatible with the mesh motion. DISCO is tested against a wide variety of problems, which are designed to test its stability, accuracy and scalability. In addition, several magnetohydrodynamics tests are performed which demonstrate the accuracy and stability of the new constrained transport approach, including two tests of the magneto-rotational instability (MRI); one testing the linear growth rate and the other following the instability into the fully turbulent regime.
We present a high-precision proper motion study of 873 X-ray and spectroscopically selected stars in the massive OB association Cygnus OB2 as part of the DANCe project. These were calculated from images spanning a 15 year baseline and have typical precisions < 1 mas/yr. We calculate the velocity dispersion in the two axes to be $\sigma_\alpha(c) = 13.0^{+0.8}_{-0.7}$ and $\sigma_\delta(c) = 9.1^{+0.5}_{-0.5}$ km/s, using a 2-component, 2-dimensional model that takes into account the uncertainties on the measurements. This gives a 3-dimensional velocity dispersion of $\sigma_{3D} = 17.8 \pm 0.6$ km/s implying a virial mass significantly larger than the observed stellar mass, confirming that the association is gravitationally unbound. The association appears to be dynamically unevolved, as evidenced by considerable kinematic substructure, non-isotropic velocity dispersions and a lack of energy equipartition. The proper motions show no evidence for a global expansion pattern, with approximately the same amount of kinetic energy in expansion as there is in contraction, which argues against the association being an expanded star cluster disrupted by process such as residual gas expulsion or tidal heating. The kinematic substructures, which appear to be close to virial equilibrium and have typical masses of 40-400 M$_\odot$, also do not appear to have been affected by the expulsion of the residual gas. We conclude that Cyg OB2 was most likely born highly substructured and globally unbound, with the individual subgroups born in (or close to) virial equilibrium, and that the OB association has not experienced significant dynamical evolution since then.
Photoionization produces supra-thermal electrons, electrons with much more energy than is found in a thermalized gas at electron temperatures characteristic of nebulae. The presence of these high energy electrons may solve the long-standing t^2/ADF puzzle, the observations that abundances obtained from recombination and collisionally excited lines do not agree, and that different temperature indicators give different results, if they survive long enough to affect diagnostic emission lines. The presence of these non-Maxwellian distribution electrons is usually designated by the term kappa. Here we use well-established methods to show that the distance over which heating rates change are much longer than the distance supra thermal electrons can travel, and that the timescale to thermalize these electrons are much shorter than the heating or cooling timescales. These estimates establish that supra thermal electrons will have disappeared into the Maxwellian velocity distribution long before they affect the collisionally excited forbidden and recombination lines that are used for deriving abundances relative to hydrogen. The electron velocity distribution in nebulae should be closely thermal.
The understanding of organic content of protoplanetary discs is one of the main goals of the planet formation studies. As an attempt to guide the observational searches for weak lines of complex species in discs, we modelled the (sub-)millimetre spectrum of gaseous methanol (CH$_3$OH), one of the simplest organic molecules, in the representative T Tauri system. We used 1+1D disc physical model coupled to the gas-grain ALCHEMIC chemical model with and without 2D-turbulent mixing. The computed CH$_3$OH abundances along with the CH$_3$OH scheme of energy levels of ground and excited torsional states were used to produce model spectra obtained with the non-local thermodynamic equilibrium (non-LTE) 3D line radiative transfer code LIME. We found that the modelled non-LTE intensities of the CH$_3$OH lines can be lower by factor of $>10$--$100$ than those calculated under assumption of LTE. Though population inversion occurs in the model calculations for many (sub-)millimetre transitions, it does not lead to the strong maser amplification and noticeably high line intensities. We identify the strongest CH$_3$OH (sub-)millimetre lines that could be searched for with the Atacama Large Millimeter Array (ALMA) in nearby discs. The two best candidates are the CH$_{3}$OH $5_0-4_0~A^+$ (241.791 GHz) and $5_{-1}-4_{-1}~E$ (241.767 GHz) lines, which could possibly be detected with the $\sim5\sigma$ signal-to-noise ratio after $\sim3$ hours of integration with the full ALMA array.
We present the effective temperatures ($T_{\rm eff}$), metallicities, and colours in SDSS, 2MASS, and WISE filters, of a sample of 3834 late-K and early-M dwarfs selected from the Sloan Digital Sky Survey APOGEE spectroscopic survey ASPCAP catalog. We confirm that ASPCAP $T_{\rm eff}$ values between 3550 K$<T_{\rm eff}<$4200 K are accurate to $\sim$100 K compared to interferometric $T_{\rm eff}$ values. In that same $T_{\rm eff}$ range, ASPCAP metallicities are accurate to 0.18 dex between $-1.0<$[M/H]$<0.2$. For these cool dwarfs, nearly every colour is sensitive to both $T_{\rm eff}$ and metallicity. Notably, we find that $g-r$ is not a good indicator of metallicity for near-solar metallicity early-M dwarfs. We confirm that $J-K_S$ colour is strongly dependent on metallicity, and find that $W1-W2$ colour is a promising metallicity indicator. Comparison of the late-K and early-M dwarf colours, metallicities, and $T_{\rm eff}$ to those from three different model grids shows reasonable agreement in $r-z$ and $J-K_S$ colours, but poor agreement in $u-g$, $g-r$, and $W1-W2$. Comparison of the metallicities of the KM dwarf sample to those from previous colour-metallicity relations reveals a lack of consensus in photometric metallicity indicators for late-K and early-M dwarfs. We also present empirical relations for $T_{\rm eff}$ as a function of $r-z$ colour combined with either [M/H] or $W1-W2$ colour, and for [M/H] as a function of $r-z$ and $W1-W2$ colour. These relations yield $T_{\rm eff}$ to $\sim$100 K and [M/H] to $\sim$0.18 dex precision with colours alone, for $T_{\rm eff}$ in the range of 3550-4200 K and [M/H] in the range of $-$0.5-0.2.
We report the discovery of an optical Einstein Ring in the Sculptor constellation, IAC J010127-334319, in the vicinity of the Sculptor Dwarf Spheroidal Galaxy. It is an almost complete ring ($\sim 300^{\circ}$) with a diameter of $\sim 4.5\, {\rm arcsec}$. The discovery was made serendipitously from inspecting Dark Energy Camera (DECam) archive imaging data. Confirmation of the object nature has been obtained by deriving spectroscopic redshifts for both components, lens and source, from observations at the $10.4$ m Gran Telescopio CANARIAS (GTC) with the spectrograph OSIRIS. The lens, a massive early-type galaxy, has a redshift of ${\rm z}=0.581$ while the source is a starburst galaxy with redshift of ${\rm z}=1.165$. The total enclosed mass that produces the lensing effect has been estimated to be ${\rm M_{tot}=(1.86 \pm 0.23) \,\cdot 10^{12}\, {\rm M_{\odot}}}$.
EGB 6 is a faint, large, ancient planetary nebula (PN). Its central star, a hot DAOZ white dwarf (WD), is a prototype of a rare class of PN nuclei associated with dense, compact emission-line knots. The central star also shows excess fluxes in both the near- (NIR) and mid-infrared (MIR). In a 2013 paper, we used Hubble Space Telescope (HST) images to show that the compact nebula is a point-like source, located 0".16 (~118 AU) from the WD. We attributed the NIR excess to an M dwarf companion star, which appeared to coincide with the dense emission knot. We now present new ground-based NIR spectroscopy, showing that the companion is actually a much cooler source with a continuous spectrum, apparently a dust-enshrouded low-luminosity star. New HST images confirm common proper motion of the emission knot and red source with the WD. The I-band, NIR, and MIR fluxes are variable, possibly on timescales as short as days. We can fit the spectral-energy distribution with four blackbodies (the WD, a ~1850 K NIR component, and MIR dust at 385 and 175 K). Alternatively, we show that the NIR/MIR SED is very similar to that of Class 0/I young stellar objects. We suggest a scenario in which the EGB 6 nucleus is descended from a wide binary similar to the Mira system, in which a portion of the wind from an AGB star was captured into an accretion disk around a companion star; a remnant of this disk has survived to the present time, and is surrounded by gas photoionized by UV radiation from the WD.
Links to: arXiv, form interface, find, astro-ph, recent, 1605, contact, help (Access key information)