In the interstellar medium of galaxies and the intracluster gas of galaxy clusters, the charged particles making up cosmic rays are moving almost exclusively along (but not across) magnetic field lines. The resulting anisotropic transport of cosmic rays in the form of diffusion or streaming not only affects the gas dynamics but also rearranges the magnetic fields themselves. The coupled dynamics of magnetic fields and cosmic rays can thus impact the formation and evolution of galaxies and the thermal evolution of galaxy clusters in critical ways. Numerically studying these effects requires solvers for anisotropic diffusion that are accurate, efficient, and robust, requirements that have proven difficult to satisfy in practice. Here, we present an anisotropic diffusion solver on an unstructured moving mesh that is conservative, does not violate the entropy condition, allows for semi-implicit time integration with individual timesteps, and only requires solving a single linear system of equations per timestep. We apply our new scheme to a large number of test problems and show that it works as well or better than previous implementations. Finally, we demonstrate for a numerically demanding simulation of the formation of an isolated disk galaxy that our local time-stepping scheme reproduces the results obtained with global time-stepping at a fraction of the computational cost.
We report new properties of the 11 and 12.7 {\mu}m emission complexes of polycyclic aromatic hydrocarbons (PAHs) by applying a Gaussian-based decomposition technique. Using high-resolution \textit{Spitzer} Space Telescope data, we study in detail the spectral and spatial characteristics of the 11 and 12.7 {\mu}m emission bands in maps of reflection nebulae NGC 7023 and NGC 2023 (North and South) and the star-forming region M17. Profile variations are observed in both the 11 and 12.7 {\mu}m emission bands. We identify a neutral contribution to the traditional 11.0 {\mu}m PAH band and a cationic contribution to the traditional 11.2 {\mu}m band, the latter of which affects the PAH class of the 11.2 {\mu}m emission in our sample. The peak variations of the 12.7 {\mu}m complex are explained by the competition between two underlying blended components. The spatial distributions of these components link them to cations and neutrals. We conclude that the 12.7 {\mu}m emission originates in both neutral and cationic PAHs, lending support to the use of the 12.7/11.2 intensity ratio as a charge proxy.
We present deep imaging of the most distant dwarf discovered by the Dark Energy Survey, Eridanus II (Eri II). Our Magellan/Megacam stellar photometry reaches $\sim$$3$ mag deeper than previous work, and allows us to confirm the presence of a stellar cluster whose position is consistent with Eri II's center. This makes Eri II, at $M_V=-7.1$, the least luminous galaxy known to host a (possibly central) cluster. The cluster is partially resolved, and at $M_V=-3.5$ it accounts for $\sim$$4\%$ of Eri II's luminosity. We derive updated structural parameters for Eri II, which has a half-light radius of $\sim$$280$ pc and is elongated ($\epsilon$$\sim$$0.48$), at a measured distance of $D$$\sim$$370$ kpc. The color-magnitude diagram displays a blue, extended horizontal branch, as well as a less populated red horizontal branch. The presence of the latter, together with a central concentration of stars brighter than the old main sequence turnoff, hints at a possible intermediate-age ($\sim$$3$ Gyr) population. Alternatively, these sources could be blue straggler stars. A deep Green Bank Telescope observation of Eri II reveals no associated atomic gas.
The confirmation of a globular cluster (GC) in the recently discovered ultrafaint galaxy Eridanus II (Eri II) motivated us to examine the question posed in the title. After estimating the halo mass of Eri II using a published stellar mass - halo mass relation, the one GC in this galaxy supports extending the relationship between the number of GCs hosted by a galaxy and the galaxy's total mass about two orders of magnitude in stellar mass below the previous limit. For this empirically determined specific frequency of between 0.06 and 0.39 globular clusters per 10$^9$ $M_\odot$ of total mass, the surviving Milky Way (MW) subhalos with masses smaller than $10^{10} M_\odot$ could host as many as 5 to 31 GCs, broadly consistent with the actual population of outer halo MW GCs, although matching the radial distribution in detail remains a challenge. Using a subhalo mass function from published high resolution numerical simulations and a Poissonian model for populating those halos with the aforementioned empirically constrained frequency, we find that about 90$\%$ of these GCs lie in lower-mass subhalos than that of Eri II. From what we know about the stellar mass-halo mass function, the subhalo mass function, and the mass-normalized GC specific frequency, we conclude that some of the MW's outer halo GCs are likely to be hosted by undetected subhalos with extremely modest stellar populations.
Current observational evidence suggests that the star formation rate (SFR) efficiency of neutral atomic hydrogen gas measured in Damped Ly-alpha Systems (DLAs) at z~3 is more than 10 times lower than predicted by the Kennicutt-Schmidt (KS) relation. To understand the origin of this deficit, and to investigate possible evolution with redshift and galaxy properties, we measure the SFR efficiency of atomic gas at z~1, z~2, and z~3 around star-forming galaxies. We use new robust photometric redshifts in the Hubble Ultra Deep Field to create galaxy stacks in these three redshift bins, and measure the SFR efficiency by combining DLA absorber statistics with the observed rest-frame UV emission in the galaxies' outskirts. We find that the SFR efficiency of HI gas at z>1 is ~1-3% of that predicted by the KS relation. Contrary to simulations and models that predict a reduced SFR efficiency with decreasing metallicity and thus with increasing redshift, we find no significant evolution in the SFR efficiency with redshift. Our analysis instead suggests that the reduced SFR efficiency is driven by the low molecular content of this atomic-dominated phase, with metallicity playing a secondary effect in regulating the conversion between atomic and molecular gas. This interpretation is supported by the similarity between the observed SFR efficiency and that observed in local atomic-dominated gas, such as in the outskirts of local spiral galaxies and local dwarf galaxies.
We study how the cosmic environment affects galaxy evolution in the Universe by comparing the metallicities of dwarf galaxies in voids with dwarf galaxies in more dense regions. Ratios of the fluxes of emission lines, particularly those of the forbidden [OIII] and [SII] transitions, provide estimates of a region's electron temperature and number density. From these two quantities and the emission line fluxes [OII] 3727, [OIII] 4363, and [OIII] 4959,5007, we estimate the abundance of oxygen with the Direct Te method. We estimate the metallicity of 37 void dwarf galaxies and 75 dwarf galaxies in more dense regions using spectroscopic observations from the Sloan Digital Sky Survey Data Release 7, as re-processed in the MPA-JHU value-added catalog. We find very little difference between the two sets of galaxies, indicating little influence from the large-scale environment on their chemical evolution. Of particular interest are a number of extremely metal-poor dwarf galaxies that are equally abundant in both voids and denser regions.
When dark matter halos are accreted by massive host clusters, strong gravitational tidal forces begin stripping mass from the accreted subhalos. This stripping eventually removes all mass beyond a subhalo's tidal radius, but the unbound mass remains in the vicinity of the satellite for at least a dynamical time t_dynamical. The N-body subhalo study of Chamberlain et al. verified this picture and pointed out a useful observational consequence: measurements of subhalo correlations beyond the tidal radius are sensitive to the infall time, t_infall, of the subhalo onto its host. We perform this cross-correlation measurement using ~ 160,000 red satellite galaxies in SDSS redMaPPer clusters and find evidence that subhalo correlations do persist well beyond the tidal radius, suggesting that many of the observed satellites fell into their current host less than a dynamical time ago, t_infall < t_dyn. Combined with estimated dynamical times t_dynamical ~ 3-5 Gyr and SED fitting results for the time at which satellites stopped forming stars, t_quench ~ 6 Gyr, we infer that for a significant fraction of the satellites, star formation quenched before those satellites entered their current hosts. The result holds for red satellites over a large range of cluster-centric distances 0.1 - 0.6 Mpc/h. We discuss the implications of this result for models of galaxy formation.
Chemical evolution models are powerful tools for interpreting stellar abundance surveys and understanding galaxy evolution. However, their predictions depend heavily on the treatment of inflow, outflow, star formation efficiency (SFE), the IMF, the SNIa delay time distribution, stellar yields, and mixing of stellar populations. Using flexCE, a new, flexible one-zone chemical evolution code, we investigate the effects of individual parameters and the trade-offs between them. Two of the most important parameters are the SFE and outflow mass-loading parameter, which shift the knee in [O/Fe]-[Fe/H] and the equilibrium abundances, respectively. One-zone models with simple star formation histories follow narrow tracks in [O/Fe]-[Fe/H] that do not match the observed bimodality in this plane. A mix of one-zone models with variations in their inflow timescales and outflow mass-loading parameters, as motivated by the inside-out galaxy formation scenario with radial mixing, reproduces the high- and low-alpha sequences better than a single model with two infall epochs. We present [X/Fe]-[Fe/H] tracks for 20 elements assuming three different SN yield models and find some significant discrepancies with observations, especially for elements with strongly metallicity-dependent yields. Analyzing the high dimensional abundance space probed by surveys like APOGEE, GALAH, and Gaia-ESO will require more advanced statistical techniques. We applied one such technique, principal component abundance analysis, to the simulations and data to reveal the main correlations amongst abundances and quantify their contributions to variation in abundance space. PC1 and PC2 of the stellar population mixing scenario are dominated by alpha-elements and elements with metallicity-dependent yields, respectively, and they collectively explain 99% of the variance. flexCE is available at https://github.com/bretthandrews/flexCE.
We present new observations of the field containing the z=3.786 protocluster, PC217.96+32.3. We confirm that it is one of the largest and most overdense high-redshift structures known. Such structures are rare even in the largest cosmological simulations. We used the Mayall/MOSAIC1.1 imaging camera to image a 1.2x0.6 deg area (~150x75 comoving Mpc) surrounding the protocluster's core and discovered 165 candidate Lyman Alpha emitting galaxies (LAEs) and 788 candidate Lyman Break galaxies (LBGs). There are at least 2 overdense regions traced by the LAEs, the largest of which shows an areal overdensity in its core (i.e., within a radius of 2.5 comoving Mpc) of 14+/-7 relative to the average LAE spatial density in the imaged field. Further, the average LAE spatial density in the imaged field is twice that derived by other field LAE surveys. Spectroscopy with Keck/DEIMOS yielded redshifts for 164 galaxies (79 LAEs and 85 LBGs); 65 lie at a redshift of 3.785+/-0.010. The velocity dispersion of galaxies near the core is 350+/-40 km/s, a value robust to selection effects. The overdensities are likely to collapse into systems with present-day masses of >10^{15} solar masses and >6x10^{14} solar masses. The low velocity dispersion may suggest a dynamically young protocluster. We find a weak trend between narrow-band (Lyman Alpha) luminosity and environmental density: the Lyman Alpha luminosity is enhanced on average by 1.35X within the protocluster core. There is no evidence that the Lyman Alpha equivalent width depends on environment. These suggest that star-formation and/or AGN activity is enhanced in the higher density regions of the structure. PC217.96+32.3 is a Coma cluster analog, witnessed in the process of formation.
The outskirts of galaxies - especially the very extended HI disks of galaxies - are strongly affected by their local environment. I highlight the giant 2X-HI disks of nearby galaxies (M 83, NGC 3621, and NGC 1512), studied as part of the Local Volume HI Survey (LVHIS), their kinematics and relation to XUV disks, signatures of tidal interactions and accretion events, the MHI - DHI relation as well as the formation of tidal dwarf galaxies. - Using multi-wavelength data, I create 3D visualisations of the gas and stars in galaxies, with the shape of their warped disks obtained through kinematic modelling of their HI velocity fields.
Massive black-hole binaries (MBHBs) are thought to be the main source of gravitational waves (GWs) in the low-frequency domain surveyed by ongoing and forthcoming Pulsar Timing Array campaigns and future space-borne missions, such as {\it eLISA}. However, many low-redshift MBHBs in realistic astrophysical environments may not reach separations small enough to allow significant GW emission, but rather stall on (sub)pc-scale orbits. This "last-parsec problem" can be eased by the appearance of a third massive black hole (MBH) -- the "intruder" -- whose action can force, under certain conditions, the inner MBHB on a very eccentric orbit, hence allowing intense GW emission eventually leading to coalescence. A detailed assessment of the process, ultimately driven by the induced Kozai-Lidov oscillations of the MBHB orbit, requires a general relativistic treatment and the inclusion of external factors, such as the Newtonian precession of the intruder orbit in the galactic potential and its hardening by scattering off background stars. In order to tackle this problem, we developed a three-body Post-Newtonian (PN) code framed in a realistic galactic potential, including both non-dissipative 1PN and 2PN terms, and dissipative terms such as 2.5PN effects, orbital hardening of the outer binary, and the effect of the dynamical friction on the early stages of the intruder dynamics. In this first paper of a series devoted at studing the dynamics of MBH triplets from a cosmological perspective, we describe, test and validate our code.
Here I present results from individual galaxy studies and galaxy surveys in the Local Universe with particular emphasis on the spatially resolved properties of neutral hydrogen gas. The 3D nature of the data allows detailed studies of the galaxy morphology and kinematics, their relation to local and global star formation as well as galaxy environments. I use new 3D visualisation tools to present multi-wavelength data, aided by tilted-ring models of the warped galaxy disks. Many of the algorithms and tools currently under development are essential for the exploration of upcoming large survey data, but are also highly beneficial for the analysis of current galaxy surveys.
The HI and CO components of the interstellar medium (ISM) are usually used to derive the dynamical mass M_dyn of nearby galaxies. Both components become too faint to be used as a tracer in observations of high-redshift galaxies. In those cases, the 158 $\mu$m line of atomic carbon [CII] may be the only way to derive M_dyn. As the distribution and kinematics of the ISM tracer affects the determination of M_dyn, it is important to quantify the relative distributions of HI, CO and [CII]. HI and CO are well-characterised observationally, however, for [CII] only very few measurements exist. Here we compare observations of CO, HI, and [CII] emission of a sample of nearby galaxies, drawn from the HERACLES, THINGS and KINGFISH surveys. We find that within R_25, the average [CII] exponential radial profile is slightly shallower than that of the CO, but much steeper than the HI distribution. This is also reflected in the integrated spectrum ("global profile"), where the [CII] spectrum looks more like that of the CO than that of the HI. For one galaxy, a spectrally resolved comparison of integrated spectra was possible; other comparisons were limited by the intrinsic line-widths of the galaxies and the coarse velocity resolution of the [CII] data. Using high-spectral-resolution SOFIA [CII] data of a number of star forming regions in two nearby galaxies, we find that their [CII] linewidths agree better with those of the CO than the HI. As the radial extent of a given ISM tracer is a key input in deriving M_dyn from spatially unresolved data, we conclude that the relevant length-scale to use in determining M_dyn based on [CII] data, is that of the well-characterised CO distribution. This length scale is similar to that of the optical disk.
Signal eigenvectors and components analyser (RF Utility model patent 116242) was used to explore the time-series of the electric field vertical component Ez in the Earth atmosphere boundary layer. There have been detected non-coherent complex-periodic components localized at the frequencies of gravity-wave impact of binary stars and at the frequency of axion-photon interaction. These components cannot be detected by means of quadrature scheme of spectral analysis and have RMS values from 0.05 V/m to 0.5 V/m at binary stars gravity-wave impact frequencies and from 0.7 V/m to 2.7 V/m at axion-photon interaction frequency. It was also demonstrated that the axion-photon interaction frequency modulates the amplitude
Context. The APOGEE survey has obtained high-resolution infrared spectra of more than 100,000 stars. Deriving chemical abundances patterns of these stars is paramount to piecing together the structure of the Milky Way. While the derived chemical abundances have been shown to be precise for most stars, some calibration problems have been reported, in particular for more metal- poor stars. Aims. In this paper, we aim to (1) re-determine the chemical abundances of the APOGEE+Kepler stellar sample (APOKASC) with an independent procedure, line list and line selection, and high quality surface gravity information from astroseismology, and (2) extend the abundance catalogue by including abundances that are not currently reported in the most recent APOGEE release (DR12). Methods. We fixed the Teff and log g to those determined using spectrophotometric and asteroseismic techniques, respectively. We made use of the Brussels Automatic Stellar Parameter (BACCHUS) code to derive the metallicity and broadening parameters for the APOKASC sample. In addition, we derived differential abundances with respect to Arcturus. Results. We have validated the BACCHUS code on APOGEE data using several well-known stars, and stars from open and globular clusters. We also provide the abundances of C, N, O, Mg, Ca, Si, Ti, S, Al, Na, Ni, Mn, Fe, K, P, Cr, Co, Cu, Rb, Yb and V for every star, line, and show the impact of line selection on the final abundances. These include abundances of five new elements and improved abundances for Si, Ti, S, and V. Conclusions. In this paper, we present an independent analysis of the APOKASC sample and provide abundances of up to 21 elements. This catalogue can be used not only to study chemical abundance patterns of the Galaxy but also to train data driven spectral approaches which can improve the abundance precision in a restricted dataset, but also full APOGEE sample.
The study of the evolution of the morphological distribution of galaxies in different environments can provide important information about the effects of the environment and the physical mechanisms responsible for the morphological transformations. As part of a complete analysis of the young cluster RXJ1257+4738 at z$\sim$0.9, we studied in this work the morphological properties of its galaxies. We used non-parametric methods of morphological classification, as implemented in the galSVM code. The classification with the applied method was possible even using ground-based observations: r'-band imaging from OSIRIS/GTC. We defined very conservative probability limits, taking into account the probability errors, in order to obtain a trustworthy classification. In this way we were able to classify about the 30% of all cluster members, and to separate between LT and ET galaxies. Additionally, when analysing the colour-magnitude diagram, we observed a significant population of blue ET galaxies between the classified ones. We discussed possible explanations for the finding of this population. Moreover, we studied different physical properties of LT, ET, and blue ET galaxies. They turn out to be comparable, with the exception of the stellar mass that shows that the red ET population is more massive. We also analysed the morphology-density and morphology-radius relations, observing that, only when considering the morphological separation between ET and LT galaxies, a mild classical behaviour is obtained. RXJ1257+4738 is a young galaxy cluster, showing a clumpy structure and being still in the process of formation, which could explain the lack of some of the standard morphological relations. This makes this cluster a very attractive case for obtaining the higher resolution data and for studying in more details the morphological properties of the entire cluster and relation with the environment.
We present the HI content of galaxies in nearby groups and clusters as measured by the 70% complete Arecibo Legacy Fast-ALFA (ALFALFA) survey, including constraints from ALFALFA detection limits. Our sample includes 22 systems at distances between 70-160 Mpc over the mass range 12.5<log M/M_sun<15.0, for a total of 1986 late-type galaxies. We find that late-type galaxies in the centers of groups lack HI at fixed stellar mass relative to the regions surrounding them. Larger groups show evidence of a stronger dependence of HI properties on environment, despite a similar dependence of color on environment at fixed stellar mass. We compare several environment variables to determine which is the best predictor of galaxy properties; group-centric distance r and r/R_200 are similarly effective predictors, while local density is slightly more effective and group size and halo mass are slightly less effective. While both central and satellite galaxies in the blue cloud exhibit a significant dependence of HI content on local density, only centrals show a strong dependence on stellar mass, and only satellites show a strong dependence on halo mass. Finally, we see evidence that HI is deficient for blue cloud galaxies in denser environments even when both stellar mass and color are fixed. This is consistent with a picture where HI is removed or destroyed, followed by reddening within the blue cloud. Our results support the existence of pre-processing in isolated groups, along with an additional rapid mechanism for gas removal within larger groups and clusters, perhaps ram-pressure stripping.
In this paper, we investigate numerically the stochastic ABC model, a toy model in the theory of astrophysical kinematic dynamos, within the recently proposed supersymmetric theory of stochastics (STS). STS characterises stochastic differential equations (SDEs) by the spectrum of the stochastic evolution operator (SEO) on elements of the exterior algebra or differentials forms over the system's phase space, X. STS can thereby classify SDEs as chaotic or non-chaotic by identifying the phenomenon of stochastic chaos with the spontaneously broken topological supersymmetry that all SDEs possess. We demonstrate the following three properties of the SEO, deduced previously analytically and from physical arguments: the SEO spectra for zeroth and top degree forms never break topological supersymmetry, all SDEs possesses pseudo-time-reversal symmetry, and each de Rahm cohomology class provides one supersymmetric eigenstate. Our results also suggests that the SEO spectra for forms of complementary degrees, i.e., k and dim X -k, may be isospectral.
Links to: arXiv, form interface, find, astro-ph, recent, 1605, contact, help (Access key information)
The "too big to fail" problem is revisited by studying the tidal evolution of populations of dwarf satellites with different density profiles. The high resolution cosmological $\rm \Lambda CDM$ "ErisMod" set of simulations is used. These simulations can model both the stellar and dark matter components of the satellites, and their evolution under the action of the tides of a MW-sized host halo at a force resolution better than 10 pc. The stronger tidal mass loss and re-shaping of the mass distribution induced in satellites with $\gamma=0.6$ dark matter density distributions, as those resulting from the effect of feedback in hydrodynamical simulations of dwarf galaxy formation, is sufficient to bring the circular velocity profiles in agreement with the kinematics of MW's dSphs. In contrast, in simulations in which the satellites retain cusps at $z=0$ there are several "massive failures" with circular velocities in excess of the observational constraints. Various sources of deviations in the conventionally adopted relation between the circular velocity at the half light radius and the one dimensional line-of-sight velocity dispersions are found. Such deviations are caused by the response of circular velocity profiles to tidal effects, which also varies depending on the initially assumed inner density profile, and by the complexity of the stellar kinematics, which include residual rotation and anisotropy. In addition tidal effects naturally induce large deviations in the stellar mass-halo mass relation for halo masses below $\rm 10^9 ~ M_{\odot}$, preventing any reliable application of the abundance matching technique to dwarf galaxy satellites.
We investigate the multiplicity of extragalactic sources detected by the Herschel Space Observatory in the COSMOS field. Using 3.6- and 24-$\mu$m catalogues, in conjunction with 250-$\mu$m data from Herschel, we seek to determine if a significant fraction of Herschel sources are composed of multiple components emitting at 250 $\mu$m. We use the XID+ code, using Bayesian inference methods to produce probability distributions of the possible contributions to the observed 250-$\mu$m flux for each potential component. The fraction of Herschel flux assigned to the brightest component is highest for sources with total 250-$\mu$m fluxes < 45 mJy; however, the flux in the brightest component is still highest in the brightest Herschel sources. The faintest 250-$\mu$m sources (30-45 mJy) have the majority of their flux assigned to a single bright component; the second brightest component is typically significantly weaker, and contains the remainder of the 250-$\mu$m source flux. At the highest 250-$\mu$m fluxes (45-110 mJy), the brightest and second brightest components are assigned roughly equal fluxes, and together are insufficient to reach 100 per cent of the 250-$\mu$m source flux. This indicates that additional components are required, beyond the brightest two components, to reproduce the observed flux. 95 per cent of the sources in our sample have a second component that contains more than 10 per cent of the total source flux. Particularly for the brightest Herschel sources, assigning the total flux to a single source may overestimate the flux contributed by around 150 per cent.
Cumulative number density matching of galaxies is a method to observationally connect descendent galaxies to their typical main progenitors at higher redshifts and thereby to assess the evolution of galaxy properties. The accuracy of this method is limited due to galaxy merging and scatter in the stellar mass growth history of individual galaxies. Behroozi et al. (2013) have introduced a refinement of the method, based on abundance matching of observed galaxies to the Bolshoi dark-matter-only simulation. The EAGLE cosmological hydro-simulation is well suited to test this method, because it reproduces the observed evolution of the galaxy stellar mass function and has a representative sample of passive/active galaxies. We find agreement with the Behroozi et al. (2013) method for the complete sample of main progenitors of z = 0 galaxies, but we also find a strong dependence on the current star formation rate. Passive galaxies with a stellar mass up to 10^10.75 Msun have a completely different median mass history than active galaxies of the same mass. This difference persists if we only select central galaxies. This means that the cumulative number density method should be applied separately to active and passive galaxies. Even then, the typical main progenitor of a z = 0 galaxy already spans two orders of magnitude in stellar mass at z = 2. Any observational analysis using cumulative number density matching therefore needs to be tested on a range of simulations.
We present deep number counts at 450 $\mu$m and 850 $\mu$m using the SCUBA-2 camera on the James Clerk Maxwell Telescope. We combine data for six lensing cluster fields and three blank fields to measure the counts over a wide flux range at each wavelength. Thanks to the lensing magnification, our measurements extends to fluxes fainter than 1 mJy and 0.2 mJy at 450 $\mu$m and 850 $\mu$m, respectively. Our combined data highly constrain the faint end of the number counts. Integrating our counts shows that the majority of the extragalactic background light (EBL) at each wavelength is contributed by faint sources with $L_{\rm IR} < 10^{12} L_{\odot }$, corresponding to luminous infrared galaxies (LIRGs) or normal galaxies. By comparing our result with the 500 $\mu$m stacking of $K$-selected sources from the literature, we conclude that the $K$-selected LIRGs and normal galaxies still cannot fully account for the EBL that originates from sources with $L_{\rm IR} < 10^{12} L_{\odot }$. This suggests that many faint submillimeter galaxies may not be included in the UV star formation history. We also explore the submillimeter flux ratio between the two bands for our 450 $\mu$m and 850 $\mu$m selected sources. At 850 $\mu$m, we find a clear relation between the flux ratio and the observed flux. This relation can be explained by a redshift evolution, where galaxies at higher redshifts have higher luminosities and star formation rates. In contrast, at 450 $\mu$m, we do not see a clear relation between the flux ratio and the observed flux.
We present full-polarization observations of the compact, steep-spectrum radio quasar 3C~286 made with the ALMA at 1.3~mm. These are the first full-polarization ALMA observations, which were obtained in the framework of Science Verification. A bright core and a south-west component are detected in the total intensity image, similar to previous centimeter images. Polarized emission is also detected toward both components. The fractional polarization of the core is about 17\%, this is higher than the fractional polarization at centimeter wavelengths, suggesting that the magnetic field is even more ordered in the millimeter radio core than it is further downstream in the jet. The observed polarization position angle (or EVPA) in the core is $\sim$\,$39^{\circ}$, which confirms the trend that the EVPA slowly increases from centimeter to millimeter wavelengths. With the aid of multi-frequency VLBI observations, we argue that this EVPA change is associated with the frequency-dependent core position. We also report a serendipitous detection of a sub-mJy source in the field of view, which is likely to be a submillimeter galaxy.
For the first time, we obtain chemical evolution models (CEMs) for Tucana and Cetus, two isolated dwarf spheroidal galaxies (dSphs) of the Local Group. The CEMs have been built from the star formation histories (SFHs) and the metallicity histories, both obtained independently by the LCID project from deep color-magnitude dia- grams. Based on our models, we find that the chemical histories were complex and can be divided into different epochs and scenarios. In particular, during 75 percent of the SFH, the galaxies behaved as closed boxes and, during the remaining 25 percent, either received a lot of primordial gas by accretion or they lost metals through metal-rich winds. In order to discriminate between these two scenarios, abundances ratios in old stars are needed. At t approximately 4.5 Gyr, the galaxies lost most of their gas due to a short-strong, well-mixed wind. We obtain very similar CEMs for both galaxies, although Cetus is twice as massive as Tucana. We conclude that the star formation in both galaxies began with only 1.5 percent of the baryonic mass fraction predicted by lambda CDM.
Based on the Sloan Digital Sky Survey (SDSS), we develop a new monte-carlo based method to estimate the photometric metallicity distribution function (MDF) for stars in the Milky Way. Compared with other photometric calibration methods, this method enables a more reliable determination of the MDF, in particular at the metal-poor and metal-rich ends. We present a comparison of our new method with a previous polynomial-based approach, and demonstrate its superiority. As an example, we apply this method to main-sequence stars with $0.2<g-r<0.6$, $6$ kpc$<R<9$ kpc, and in different intervals in height above the plane, $|Z|$. The MDFs for the selected stars within two relatively local intervals ($0.8$ kpc$<|Z|<1.2$ kpc, $1.5$ kpc$<|Z|<2.5$ kpc) can be well-fit by two Gaussians, with peaks at [Fe/H] $\approx-0.6$ and $-1.2$ respectively, one associated with the disk system, the other with the halo. The MDFs for the selected stars within two more distant intervals ($3$ kpc$<|Z|<5$ kpc, $6$ kpc$<|Z|<9$ kpc) can be decomposed into three Gaussians, with peaks at [Fe/H] $\approx-0.6$, $-1.4$ and $-1.9$ respectively, where the two lower peaks may provide evidence for a two-component model of the halo: the inner halo and the outer halo. The number ratio between the disk component and halo component(s) decreases with vertical distance from the Galactic plane, consistent with the previous literature.
We present spatial correlations of galaxies and IGM neutral hydrogen HI in the COSMOS/UltraVISTA 1.62 deg$^2$ field. Our data consist of 13415 photo-$z$ galaxies at $z\sim2-3$ with $K_s<23.4$ and the Ly$\alpha$ forest absorptions in the background quasar spectra selected from SDSS data with no signature of damped Ly$\alpha$ system contamination. We estimate a galaxy overdensity $\delta_{gal}$ in an impact parameter of 2.5 pMpc, and calculate the Ly$\alpha$ forest fluctuations $\delta_{\langle F\rangle}$ whose negative values correspond to the strong Ly$\alpha$ forest absorptions. We identify a weak anti-correlation between $\delta_{gal}$ and $\delta_{\langle F\rangle}$ with a Spearman's rank correlation coefficient of $-0.39$ suggesting that the galaxy overdensities and the Ly$\alpha$ forest absorptions positively correlate in space at the $\sim90\%$ confidence level. This positive correlation indicates that high-$z$ galaxies exist around an excess of HI gas in the Ly$\alpha$ forest. We find four cosmic volumes, dubbed $A_{obs}$-$D_{obs}$, that have extremely large (small) values of $\delta_{gal} \simeq0.8$ ($-1$) and $\delta_{\langle F\rangle}$ $\simeq0.1$ ($-0.4$), three out of which, $B_{obs}$-$D_{obs}$, significantly depart from the correlation, and weaken the correlation signal. We perform cosmological hydrodynamical simulations, and compare with our observational results. Our simulations reproduce the correlation, agreeing with the observational results. Moreover, our simulations have model counterparts of $A_{obs}$-$D_{obs}$, and suggest that the observations pinpoint, by chance, a galaxy overdensity like a proto-cluster, gas filaments lying on the sightline, a large void, and orthogonal low-density filaments. Our simulations indicate that the significant departures of $B_{obs}$-$D_{obs}$ are produced by the filamentary large-scale structures and the observation sightline effects.
We present a detailed study of the neutral and ionised gas phases in the galactic wind for the nearby starburst galaxy NGC 5394 based on new integral field spectroscopy obtained with the INTEGRAL fibre system at the William Herschel Telescope. The neutral gas phase in the wind is detected via the interstellar NaI D doublet absorption. After a careful removal of the stellar contribution to these lines, a significant amount of neutral gas (~10^7 Msun) is detected in a central region of ~1.75 kpc size. This neutral gas is blueshifted by ~165 km/s with respect to the underlying galaxy. The mass outflow of neutral gas is comparable to the star formation rate of the host galaxy. Simultaneously, several emission lines (Ha, [NII], [SII]) are also analysed looking for the ionised warm phase counterpart of the wind. A careful kinematic decomposition of the line profiles reveals the presence of a secondary, broader, kinematic component. This component is found roughly in the same region where the NaI D absorption is detected. It presents higher [NII]/Ha and [SII]/Ha line ratios than the narrow component at the same locations, indicative of contamination by shock ionization. This secondary component also presents blueshifted velocities, although smaller than those measured for the neutral gas, averaging to ~ -30 km/s. The mass and mass outflow rate of the wind is dominated by the neutral gas, of which a small fraction might be able to escape the gravitational potential of the host galaxy. The observations in this system can be readily understood within a bipolar gas flow scenario.
In this paper we selected a sample of massive star forming regions from the Red MSX Source (RMS) survey, to study star formation activities (mainly outflow and inflow signatures).
The supermassive black hole, Sagittarius A* (Sgr A*), at the centre of the Milky Way undergoes regular flaring activity which is thought to arise from the innermost region of the accretion flow. We performed the monitoring observations of the Galactic Centre to study the flux-density variations at 3mm using the Australia Telescope Compact Array (ATCA) between 2010 and 2014. We obtain the light curves of Sgr A* by subtracting the contributions from the extended emission around it, and the elevation and time dependent gains of the telescope. We perform structure function analysis and the Bayesian blocks representation to detect flare events. The observations detect six instances of significant variability in the flux density of Sgr A* in three observations, with variations between 0.5 to 1.0 Jy, which last for 1.5 $-$ 3 hours. We use the adiabatically expanding plasmon model to explain the short time-scale variations in the flux density. We derive the physical quantities of the modelled flare emission, such as the source expansion speed $v_{\mathrm{exp}}$, source sizes, spectral indices, and the turnover frequency. These parameters imply that the expanding source components are either confined to the immediate vicinity of Sgr A* by contributing to the corona or the disc, or have a bulk motion greater than $v_{\mathrm{exp}}$. No exceptional flux density variation on short flare time-scales was observed during the approach and the flyby of the dusty S-cluster object (DSO/G2). This is consistent with its compactness and the absence of a large bow shock.
The Gaia-ESO Survey was designed to target all major Galactic components (i.e., bulge, thin and thick discs, halo and clusters), with the goal of constraining the chemical and dynamical evolution of the Milky Way. This paper presents the methodology and considerations that drive the selection of the targeted, allocated and successfully observed Milky Way field stars. The detailed understanding of the survey construction, specifically the influence of target selection criteria on observed Milky Way field stars is required in order to analyse and interpret the survey data correctly. We present the target selection process for the Milky Way field stars observed with VLT/FLAMES and provide the weights that characterise the survey target selection. The weights can be used to account for the selection effects in the Gaia-ESO Survey data for scientific studies. We provide a couple of simple examples to highlight the necessity of including such information in studies of the stellar populations in the Milky Way.
We have incorporated our experimentally derived thermal rate coefficients for C + H$_3^+$ forming CH$^+$ and CH$_2^+$ into a commonly used astrochemical model. We find that the Arrhenius-Kooij equation typically used in chemical models does not accurately fit our data and use instead a more versatile fitting formula. At a temperature of 10 K and a density of 10$^4$ cm$^{-3}$, we find no significant differences in the predicted chemical abundances, but at higher temperatures of 50, 100, and 300 K we find up to factor of 2 changes. Additionally, we find that the relatively small error on our thermal rate coefficients, $\sim15\%$, significantly reduces the uncertainties on the predicted abundances compared to those obtained using the currently implemented Langevin rate coefficient with its estimated factor of 2 uncertainty.
We have analyzed rotational spectral line emission of OCS, CH3OH, HCOOCH3, and H2CS observed toward the low-mass Class 0 protostellar source IRAS 16293-2422 Source A at a sub-arcsecond resolution (~0".6 x 0".5) with ALMA. Significant chemical differentiation is found at a 50 AU scale. The OCS line is found to well trace the infalling-rotating envelope in this source. On the other hand, the CH3OH and HCOOCH3 distributions are found to be concentrated around the inner part of the infalling-rotating envelope. With a simple ballistic model of the infalling-rotating envelope, the radius of the centrifugal barrier (a half of the centrifugal radius) and the protostellar mass are evaluated from the OCS data to be from 40 to 60 AU and from 0.5 to 1.0 Msun, respectively, assuming the inclination angle of the envelope/disk structure to be 60 degrees (90 degrees for the edge-on configuration). Although the protostellar mass is correlated with the inclination angle, the radius of the centrifugal barrier is not. This is the first indication of the centrifugal barrier of the infalling-rotating envelope in a hot corino source. CH3OH and HCOOCH3 may be liberated from ice mantles due to weak accretion shocks around the centrifugal barrier, and/or due to protostellar heating. The H2CS emission seems to come from the disk component inside the centrifugal barrier in addition to the envelope component. The centrifugal barrier plays a central role not only in the formation of a rotationally-supported disk but also in the chemical evolution from the envelope to the protoplanetary disk.
We aim to constrain the temperature and velocity structures, and H2O abundances in the winds of a sample of M-type AGB stars. We further aim to determine the effect of H2O line cooling on the energy balance in the inner circumstellar envelope. We use two radiative-transfer codes to model molecular emission lines of CO and H2O towards four M-type AGB stars. We focus on spectrally resolved observations of CO and H2O from HIFI. The observations are complemented by ground-based CO observations, and spectrally unresolved CO and H2O observations with PAC. The observed line profiles constrain the velocity structure throughout the circumstellar envelopes (CSEs), while the CO intensities constrain the temperature structure in the CSEs. The H2O observations constrain the o-H2O and p-H2O abundances relative to H2. Finally, the radiative-transfer modelling allows to solve the energy balance in the CSE, in principle including also H2O line cooling. The fits to the line profiles only set moderate constraints on the velocity profile, indicating shallower acceleration profiles in the winds of M-type AGB stars than predicted by dynamical models, while the CO observations effectively constrain the temperature structure. Including H2O line cooling in the energy balance was only possible for the low-mass-loss-rate objects in the sample, and required an ad hoc adjustment of the dust velocity profile in order to counteract extreme cooling in the inner CSE. H2O line cooling was therefore excluded from the models. The constraints set on the temperature profile by the CO lines nevertheless allowed us to derive H2O abundances. The derived H2O abundances confirm previous estimates and are consistent with chemical models. However, the uncertainties in the derived abundances are relatively large, in particular for p-H2O, and consequently the derived o/p-H2O ratios are not well constrained.
Links to: arXiv, form interface, find, astro-ph, recent, 1605, contact, help (Access key information)
The physics of cosmic rays (CR) is a promising candidate for explaining the driving of galactic winds and outflows. Recent galaxy formation simulations have demonstrated the need for active CR transport either in the form of diffusion or streaming to successfully launch winds in galaxies. However, due to computational limitations, most previous simulations have modeled CR transport isotropically. Here, we discuss high resolution simulations of isolated disk galaxies in a $10^{11}\rm{M_\odot}$ halo with the moving mesh code {\sc Arepo} that include injection of CRs from supernovae, advective transport, CR cooling, and CR transport through isotropic or anisotropic diffusion. We show that either mode of diffusion leads to the formation of strong bipolar outflows. However, they develop significantly later in the simulation with anisotropic diffusion compared to the simulation with isotropic diffusion. Moreover, we find that isotropic diffusion allows most of the CRs to quickly diffuse out of the disk, while in the simulation with anisotropic diffusion, most CRs remain in the disk once the magnetic field becomes dominated by its azimuthal component, which occurs after $\sim 300\,{\rm Myrs}$. This has important consequences for the gas dynamics in the disk. In particular, we show that isotropic diffusion strongly suppresses the amplification of the magnetic field in the disk compared to anisotropic or no diffusion models. We therefore conclude that reliable simulations which include CR transport inevitably need to account for anisotropic diffusion.
We present the results of ALMA Cycle 2 observations of the ultraluminous infrared galaxy, IRAS 20551-4250, at HCN/HCO+/HNC J=3-2 lines at both vibrational-ground (v=0) and vibrationally excited (v2=1) levels. This galaxy contains a luminous buried active galactic nucleus (AGN), in addition to starburst activity, and our ALMA Cycle 0 data revealed a tentatively detected vibrationally excited HCN v2=1f J=4-3 emission line. In our ALMA Cycle 2 data, the HCN/HCO+/HNC J=3-2 emission lines at v=0 are clearly detected. The HCN and HNC v2=1f J=3-2 emission lines are also detected, but the HCO+ v2=1f J=3-2 emission line is not. Given the high-energy level of v2=1 and the resulting difficulty of collisional excitation, we compared these results with those of the calculation of infrared radiative pumping, using the available infrared 5-35 micron spectrum. We found that all of the observational results were reproduced, if the HCN abundance was significantly higher than that of HCO+ and HNC. The flux ratio and excitation temperature between v2=1f and v=0, after correction for possible line opacity, suggests that infrared radiative pumping affects rotational (J-level) excitation at v=0 at least for HCN and HNC. The HCN-to-HCO+ v=0 flux ratio is higher than those of starburst-dominated regions, and will increase even more when thederived high HCN opacity is corrected. The enhanced HCN-to-HCO+ flux ratio in this AGN-hosting galaxy can be explained by the high HCN-to-HCO+ abundance ratio and sufficient HCN excitation at up to J=4, rather than the significantly higher efficiency of infrared radiative pumping for HCN than HCO+.
Self-gravity and stellar feedback are capable of driving turbulence and transporting mass and angular momentum in disk galaxies, but the balance between them is not well understood. In the previous paper in this series, we showed that gravity alone can drive turbulence in galactic disks, regulate their Toomre $Q$ parameters to $\sim$ 1, and transport mass inwards at a rate sufficient to fuel star formation in the centers of present-day galaxies. In this paper we extend our models to include the effects of star formation feedback. We show that feedback suppresses galaxies' star formation rates by a factor of $\sim$ 5 and leads to the formation of a multi-phase atomic and molecular ISM. Both the star formation rate and the phase balance produced in our simulations agree well with observations of nearby spirals. After our galaxies reach steady state, we find that the inclusion of feedback actually lowers the gas velocity dispersion slightly compared to the case of pure self-gravity, and also slightly reduces the rate of inward mass transport. Nevertheless, we find that, even with feedback included, our galactic disks self-regulate to $Q$ $\sim$ 1, and transport mass inwards at a rate sufficient to supply a substantial fraction of the inner disk star formation. We argue that gravitational instability is therefore likely to be the dominant source of turbulence and transport in galactic disks, and that it is responsible for fueling star formation in the inner parts of galactic disks over cosmological times.
We present the new semi-analytic model of galaxy evolution, DARK SAGE, a heavily modified version of the publicly available SAGE code. The model is designed for detailed evolution of galactic discs. We evolve discs in a series of annuli with fixed specific angular momentum, which allows us to make predictions for the radial and angular-momentum structure of galaxies. Most physical processes, including all channels of star formation and associated feedback, are performed in these annuli. We present the surface density profiles of our model spiral galaxies, both as a function of radius and specific angular momentum, and find the discs naturally build a pseduobulge-like component. Our main results are focussed on predictions relating to the integrated mass--specific angular momentum relation of stellar discs. The model produces a distinct sequence between these properties in remarkable agreement with recent observational literature. We investigate the impact Toomre disc instabilities have on shaping this sequence and find they are crucial for regulating both the mass and spin of discs. Without instabilities, high-mass discs would be systematically deficient in specific angular momentum by a factor of ~2.5, with increased scatter. Instabilities also appear to drive the direction in which the mass--spin sequence of spiral galaxy discs evolves. With them, we find galaxies of fixed mass have higher specific angular momentum at later epochs.
Diffuse striations in molecular clouds are preferentially aligned with local magnetic fields whereas dense filaments tend to be perpendicular to them. When and why this transition occurs remain uncertain. To explore the physics behind this transition, we compute the histogram of relative orientation (HRO) between the density gradient and the magnetic field in 3D MHD simulations of prestellar core formation in shock-compressed regions within GMCs. We find that, in the magnetically-dominated (sub-Alfv\'enic) post-shock region, the gas structure is preferentially aligned with the local magnetic field. For overdense sub-regions with super-Alfv\'enic gas, their elongation becomes preferentially perpendicular to the local magnetic field instead. The transition occurs when self-gravitating gas gains enough kinetic energy from the gravitational acceleration to overcome the magnetic support against the cross-field contraction, which results in a power-law increase of the field strength with density. Similar results can be drawn from HROs in projected 2D maps with integrated column densities and synthetic polarized dust emission. We quantitatively analyze our simulated polarization properties, and interpret the reduced polarization fraction at high column densities as the result of increased distortion of magnetic field directions in trans- or super-Alfv\'enic gas. Furthermore, we introduce measures of the inclination and tangledness of the magnetic field along the line of sight as the controlling factors of the polarization fraction. Observations of the polarization fraction and angle dispersion can therefore be utilized in studying local magnetic field morphology in star-forming regions.
Photoelectric heating has long been recognized as the primary source of heating for the neutral interstellar medium. Simulations of spiral galaxies found some indication that photoelectric heating could suppress star formation. However, simulations that include photoelectric heating have typically found that it has little effect on the rate of star formation in either spiral galaxies or dwarfs suggesting that supernovae and not photoelectric heating are responsible for setting the star formation law in galaxies. This result is in tension with recent work indicating that a star formation law that depends on galaxy metallicity, as expected for photoelectric heating but not for supernovae, reproduces the present-day galaxy population better than a metallicity-independent one. Here we report a series of simulations of dwarf galaxies, where the effects of both photoelectric heating and supernovae are expected to be strongest. We simultaneously include space- and time-dependent photoelectric heating, and we resolve the Sedov phase of every supernova blast wave, allowing us to make a direct measurement of the relative importance of momentum injection by supernovae and dust heating by far ultraviolet (FUV) photons in suppressing star formation. We find that supernovae are unable to account for the long observed gas depletion times in dwarf galaxies. Instead, ordinary photoelectric heating is the dominant means by which dwarf galaxies regulate their star formation rate at any given time, suppressing the star formation rate by more than an order of magnitude relative to simulations with only supernovae.
Previous studies have shown that WISE-selected hyperluminous, hot dust-obscured galaxies (Hot DOGs) are powered by highly dust-obscured, possibly Compton-thick AGNs. High obscuration provides us a good chance to study the host morphology of the most luminous AGNs directly. We analyze the host morphology of 18 Hot DOGs at $z\sim3$ using Hubble Space Telescope/WFC3 imaging. We find that Hot DOGs have a high merger fraction ($62\pm 14 \%$). By fitting the surface brightness profiles, we find that the distribution of S\'ersic indices in our Hot DOG sample peaks around 2, which suggests that most of Hot DOGs have transforming morphologies. We also derive the AGN bolometric luminosity ($\sim10^{14}L_\odot$) of our Hot DOG sample by using IR SEDs decomposition. The derived merger fraction and AGN bolometric luminosity relation is well consistent with the variability-based model prediction (Hickox et al. 2014). Both the high merger fraction in IR-luminous AGN sample and relatively low merger fraction in UV/optical-selected, unobscured AGN sample can be expected in the merger-driven evolutionary model. Finally, we conclude that Hot DOGs are merger-driven and may represent a transit phase during the evolution of massive galaxies, transforming from the dusty starburst dominated phase to the unobscured QSO phase.
The interstellar medium (ISM) consists of highly ionized and neutral atomic, as well as molecular, components. Knowledge of their distribution is important for tracing the structure and lifecycle of the ISM. Here we determine the properties of the highly ionized and neutral weakly ionized gas in the Galaxy traced by the fine-structure lines of ionized nitrogen, [N II], and ionized carbon, [C II]. To analyze the ionized ISM we utilize [C II] 158 micron and [N II] 205 micron lines taken with the high spectral resolution Heterodyne Instrument in the Far-Infrared (HIFI) on the Herschel Space Observatory along ten lines of sight towards the inner Galaxy. [N II] emission can be used to estimate the contribution of the highly ionized gas to the [C II] emission and separate the highly ionized and weakly ionized neutral gas. We find that [N II] has strong emission in distinct spectral features along all lines of sight associated with strong [C II] emission. The [N II] arises from moderate density extended HII regions or ionized boundary layers of clouds. Comparison of the [N II] and [C II] spectra in 31 separate kinematic features shows that many of the [C II] spectra are affected by absorption from low excitation gas associated with molecular clouds, sometimes strongly so. The apparent fraction of the [C II] associated with the [N II] gas is unrealistically large in many cases, most likely due to the reduction of [C II] by absorption. In a few cases the foreground absorption can be modeled to determine the true source intensity. In these sources we find that the foreground absorbing gas layer has C$^+$ column densities of order 10$^{18}$ cm$^{-2}$. [C II] emission arising from strong sources of [N II] emission is frequently absorbed by low excitation foreground gas complicating the interpretation of the properties of the ionized and neutral gas components that give rise to [C II] emission.
New near- and far-ultraviolet (NUV and FUV) HST spectra of Mrk 231, the nearest quasar known, are combined with ground-based optical spectra to study the remarkable dichotomy between the FUV and NUV-optical spectral regions in this object. The FUV emission-line features are faint, broad, and highly blueshifted (up to ~7000 km/s), with no significant accompanying absorption. In contrast, the profiles of the NUV absorption features resemble those of the optical Na I D, He I, and Ca II H and K lines, exhibiting broad blue-shifted troughs that overlap in velocity space with the FUV emission-line features and indicate a dusty, high-density and patchy broad absorption line (BAL) screen covering ~90% of the observed continuum source at a distance less than ~2 - 20 pc. The FUV continuum emission does not show the presence of any obvious stellar features and is remarkably flat compared with the steeply declining NUV continuum. The NUV (FUV) features and continuum emission have not varied significantly over the past ~22 (3) years and are unresolved on scales ~40 (170) pc. These results favor an AGN origin for the NUV - FUV line and continuum emission. The observed FUV line emission is produced in the outflowing BAL cloud system, while the Balmer lines arise primarily from the standard broad line region seen through the dusty BAL screen. Our data are inconsistent with the recently proposed binary black hole model. We argue instead that Mrk~231 is the nearest example of weak-lined "wind-dominated" quasars with high Eddington ratios and geometrically thick ("slim") accretion disks; these quasars are likely more common in the early universe.
The surface-brightness profiles of galaxies I(R) and the density profiles of dark-matter halos rho(r) are well represented by the same analytic function, named after either S\'ersic, I~exp[-(R/R*)^(1/m)], or Einasto, rho~[exp[-(r/r*)^alpha], where R* and r* are characteristic radii. Systems with high S\'ersic index m (or low Einasto index alpha) have steep central profiles and shallow outer profiles, while systems with low m (or high alpha) have shallow central profiles and steep profiles in the outskirts. We present the results of idealized numerical experiments which suggest that the origin of these profiles can be traced back to the initial density fluctuation field: high-alpha (low-m) systems form in smooth regions via few mergers, while low-alpha (high-m) systems form in clumpy regions via several mergers.
The current IAU Symposium is closely connected to the EU-funded network DAGAL (Detailed Anatomy of Galaxies), with the final annual network meeting of DAGAL being at the core of this international symposium. In this short paper, we give an overview of DAGAL, its training activities, and some of the scientific advances that have been made under its umbrella.
We present maps of the Large and Small Magellanic Clouds from combined South Pole Telescope (SPT) and Planck data. Both instruments are designed to make measurements of the cosmic microwave background but are sensitive to any source of millimeter-wave (mm-wave) emission. The Planck satellite observes in nine mm-wave bands, while the SPT data used in this work were taken with the three-band SPT-SZ camera. The SPT-SZ bands correspond closely to three of the nine Planck bands, namely those centered at 1.4, 2.1, and 3.0 mm. The angular resolution of the Planck data in these bands ranges from 5 to 10 arcmin, while the SPT resolution in these bands ranges from 1.0 to 1.7 arcmin. The combined maps take advantage of the high resolution of the SPT data and the long-timescale stability of the space-based Planck observations to deliver high signal-to-noise and robust brightness measurements on scales from the size of the maps down to ~1 arcmin. In each of the three bands, we first calibrate and color-correct the SPT data to match the Planck data, then we use noise estimates from each instrument and knowledge of each instrument's beam, or point-spread function, to make the inverse-variance-weighted combination of the two instruments' data as a function of angular scale. We create maps assuming a range of underlying emission spectra (for the color correction) and at a range of final resolutions. We perform several consistency tests on the combined maps and estimate the expected noise in measurements of features in the maps. We compare the maps of the Large Magellanic Cloud (LMC) from this work to maps from the Herschel HERITAGE survey, finding general consistency between the datasets. The broad wavelength coverage provides evidence of different emission mechanisms at work in different environments in the LMC. [Abridged]
We present new Chandra X-ray observations of seven low-mass black holes (~1e6 Msun) accreting at low Eddington ratios between -2.0<log L/Ledd<-1.5. We compare the X-ray properties of these seven low-mass active galactic nuclei (AGN) to a total of 73 other low-mass AGN in the literature with published Chandra observations (with Eddington ratios extending from -2.0<log L/Ledd<-0.1). We do not find any statistical differences between low- and high-Eddington ratio low-mass AGN in the distributions of their X-ray to ultraviolet luminosity ratios (aox), or in their X-ray spectral shapes. Furthermore, the aox distribution of low-L/Ledd AGN displays an X-ray weak tail that is also observed within high-L/Ledd objects. Our results indicate that between -2<log L/Ledd<-0.1, there is no systematic change in the structure of the accretion flow for active galaxies hosting 1e6 Msun black holes. We examine the accuracy of current bolometric luminosity estimates for our low-L/Ledd objects with new Chandra observations, and it is plausible that their Eddington ratios could be underestimated by up to an order of magnitude. If so, then in analogy with weak emission line quasars, we suggest that accretion from a geometrically thick, radiatively inefficient `slim disk' could explain their diverse properties in aox. Alternatively, if current Eddington ratios are in fact correct (or overestimated), then the X-ray weak tail would imply that there is diversity in disk/corona couplings among individual low-mass objects. Finally, we conclude by noting that the aox distribution for low-mass black holes may have favorable consequences for the epoch of cosmic reionization being driven by AGN.
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)
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)