Links to: arXiv, form interface, find, astro-ph, recent, 1710, contact, help (Access key information)
Massive black hole (MBH) binaries, formed as a result of galaxy mergers, are expected to harden by dynamical friction and three body stellar scatterings, until emission of gravitational waves (GWs) leads to their final coalescence. According to recent simulations, MBH binaries can efficiently harden via stellar encounters only when the host geometry is triaxial, even if only modestly, as angular momentum diffusion allows an efficient repopulation of the binary loss cone. In this paper, we carry out a suite of N-body simulations of equal-mass galaxy collisions, varying the initial orbits and density profiles for the merging galaxies and running simulations both with and without central MBHs. We find that the presence of a MBH binary in the remnant makes the system maximally triaxial within the binary sphere of influence and nearly oblate, aligned with the galaxy merger plane, within a radius enclosing 100 MBH masses. We never find binary hosts to be prolate on any scale. The decaying MBHs slightly enhance the tangential anisotropy in the centre of the remnant due to angular momentum injection and the slingshot ejection of stars on nearly radial orbits. This latter effect results in about 1% of the remnant stars being expelled from the galactic nucleus. Finally, we do not find any strong connection between the remnant morphology and the binary hardening rate, which depends only on the inner density slope of the remnant galaxy. Our results suggest that MBH binaries are able to coalesce within a few Gyr, even if their host systems retain an axisymmetric shape beyond the binary sphere of influence.
We study the population of supermassive black holes (SMBHs) and their effects on massive central galaxies in the IllustrisTNG cosmological hydrodynamical simulations of galaxy formation. The employed model for SMBH growth and feedback assumes a two-mode scenario in which the feedback from active galactic nuclei occurs through a kinetic, comparatively efficient mode at low accretion rates relative to the Eddington limit, and in the form of a thermal, less efficient mode at high accretion rates. We show that the quenching of massive central galaxies happens coincidently with kinetic-mode feedback, consistent with the notion that active supermassive black cause the low specific star formation rates observed in massive galaxies. However, major galaxy mergers are not responsible for initiating most of the quenching events in our model. Up to black hole masses of about $10^{8.5}\,{\rm M}_\odot$, the dominant growth channel for SMBHs is in the thermal mode. Higher mass black holes stay mainly in the kinetic mode and gas accretion is self-regulated via their feedback, which causes their Eddington ratios to drop, with SMBH mergers becoming the main channel for residual mass growth. As a consequence, the quasar luminosity function is dominated by rapidly accreting, moderately massive black holes in the thermal mode. We show that the associated growth history of SMBHs produces a low-redshift quasar luminosity function and a redshift zero black hole mass-stellar bulge mass relation in good agreement with observations, whereas the simulation tends to over-predict the high-redshift quasar luminosity function.
We present an analysis of 55 central galaxies in clusters and groups with molecular gas masses and star formation rates lying between $10^{8}-10^{11}\ M_{\odot}$ and $0.5-270$ $M_{\odot}\ yr^{-1}$, respectively. We have used Chandra observations to derive profiles of total mass and various thermodynamic variables. Molecular gas is detected only when the central cooling time or entropy index of the hot atmosphere falls below $\sim$1 Gyr or $\sim$35 keV cm$^2$, respectively, at a (resolved) radius of 10 kpc. This indicates that the molecular gas condensed from hot atmospheres surrounding the central galaxies. The depletion timescale of molecular gas due to star formation approaches 1 Gyr in most systems. Yet ALMA images of roughly a half dozen systems drawn from this sample suggest the molecular gas formed recently. We explore the origins of thermally unstable cooling by evaluating whether molecular gas becomes prevalent when the minimum of the cooling to free-fall time ratio ($t_{\rm cool}/t_{\rm ff}$) falls below $\sim10$. We find: 1) molecular gas-rich systems instead lie between $10 < min(t_{\rm cool}/t_{\rm ff}) < 25$, where $t_{\rm cool}/t_{\rm ff}=25$ corresponds approximately to cooling time and entropy thresholds $t_{\rm cool} \lesssim 1$ Gyr and 35 keV~cm$^2$, respectively, 2) $min(t_{\rm cool}/t_{\rm ff}$) is uncorrelated with molecular gas mass and jet power, and 3) the narrow range $10 < min(t_{\rm cool}/t_{\rm ff}) < 25$ can be explained by an observational selection effect. These results and the absence of isentropic cores in cluster atmospheres are in tension with "precipitation" models, particularly those that assume thermal instability ensues from linear density perturbations in hot atmospheres. Some and possibly all of the molecular gas may instead have condensed from atmospheric gas lifted outward either by buoyantly-rising X-ray bubbles or merger-induced gas motions.
The origin and survival of the Smith high-velocity HI cloud has so far defied explanation. This object has several remarkable properties: (i) its prograde orbit is ~100 km/s faster than the underlying Galactic rotation; (ii) its total gas mass ($\gtrsim 4 \times 10^6 ~{\rm M}_{\odot}$) exceeds the mass of all other high-velocity clouds (HVC) outside of the Magellanic Stream; (iii) its head-tail morphology extends to the Galactic HI disc, indicating some sort of interaction. The Smith Cloud's kinetic energy rules out models based on ejection from the disc. We construct a dynamically self-consistent, multi-phase model of the Galaxy with a view to exploring whether the Smith Cloud can be understood in terms of an infalling, compact HVC that has transited the Galactic disc. We show that while a dark-matter (DM) free HVC of sufficient mass and density can reach the disc, it does not survive the transit. The most important ingredient to survival during a transit is a confining DM subhalo around the cloud; radiative gas cooling and high spatial resolution ($\lesssim$ 10 pc) are also essential. In our model, the cloud develops a head-tail morphology within ~10 Myr before and after its first disc crossing; after the event, the tail is left behind and accretes onto the disc within ~400 Myr. In our interpretation, the Smith Cloud corresponds to a gas 'streamer' that detaches, falls back and fades after the DM subhalo, distorted by the disc passage, has moved on. We conclude that subhalos with ${\rm M}_{\rm DM} \lesssim 10^9 ~{\rm M}_{\odot}$ have accreted $\sim10^9 ~{\rm M}_{\odot}$ of gas into the Galaxy over cosmic time - a small fraction of the total baryon budget.
The fraction of galaxies bound in groups in the nearby Universe is high (50% at z~0). Systematic studies of galaxy properties in groups are important in order to improve our understanding of the evolution of galaxies and of the physical phenomena occurring within this environment. We have built a complete spectrophotometric sample of galaxies within X-ray detected, optically spectroscopically confirmed groups and clusters (G&C), covering a wide range of halo masses at z<= 0.6. In the context of the XXL survey, we analyse a sample of 164 G&C in the XXL-North region (XXL-N), at z <= 0.6, with a wide range of virial masses (1.24 x 10^13 <=M_500 M_sun <= 6.63 x 10^14) and X-ray luminosities ( 2.27 x 10^41 <= L^XXL_500 (erg/s)<= 2.15 x10^44). The G&C are X-ray selected and spectroscopically confirmed. We describe the membership assignment and the spectroscopic completeness analysis, and compute stellar masses. As a first scientific exploitation of the sample, we study the dependence of the galaxy stellar mass function (GSMF) on global environment. We present a spectrophotometric characterisation of the G&C and their galaxies. The final sample contains 132 G&C, 22111 field galaxies and 2225 G&C galaxies with r-band magnitude <20. Of the G&C, 95% have at least three spectroscopic members, and 70% at least ten. The shape of the GSMF seems not to depend on environment (field versus G&C) or X-ray luminosity ( used as a proxy for the virial mass of the system). These results are confirmed by the study of the correlation between mean stellar mass of G&C members and L^XXL_500.We release the spectrophotometric catalogue of galaxies with all the quantities computed in this work. As a first homogeneous census of galaxies within X-ray spectroscopically confirmed G&C at these redshifts, this sample will allow environmental studies of the evolution of galaxy properties.
We present high-quality fluid dynamical simulations of isothermal gas flows in a rotating barred potential. We show that a large quantity of gas is driven right into the nucleus of a model galaxy when the potential lacks a central mass concentration, but the inflow stalls at a nuclear ring in comparison simulations that include a central massive object. The radius of the nuclear gas ring increases linearly with the mass of the central object. We argue that bars drive gas right into the nucleus in the early stages of disk galaxy formation, where a nuclear star cluster and perhaps a massive black hole could be created. The process is self-limiting, however, because inflow stalls at a nuclear ring once the mass of gas and stars in the nucleus exceeds ~1% of the disk mass, which shuts off rapid growth of the black hole. We briefly discuss the relevance of these results to the seeding of massive black holes in galaxies, the merger model for quasar evolution, and the existence of massive black holes in disk galaxies that lack a significant classical bulge.
In the next few years the classification of radio sources observed by the large surveys will be a challenging problem, and spectral index is a powerful tool for addressing it. Here we present an algorithm to estimate the spectral index of sources from multiwavelength radio images. We have applied our algorithm to SCORPIO (Umana et al. 2015), a Galactic Plane survey centred around 2.1 GHz carried out with ATCA, and found we can measure reliable spectral indices only for sources stronger than 40 times the rms noise. Above a threshold of 1 mJy, the source density in SCORPIO is 20 percent greater than in a typical extra-galactic field, like ATLAS (Norris et al. 2006), because of the presence of Galactic sources. Among this excess population, 16 sources per square degree have a spectral index of about zero, suggesting optically thin thermal emission such as Hii regions and planetary nebulae, while 12 per square degree present a rising spectrum, suggesting optically thick thermal emission such as stars and UCHii regions.
We simulate the evolution of supernova remnant (SNR) W51C. The simulation shows the existence of a new northeast edge. We present magnetic field structure of the W51 complex (SNR W51C and two HII regions W51A/B) by employing the 11 cm survey data of Effelsberg. This new edge is identified and overlaps with W51A along the line of sight, which gives a new angular diameter of about 37' for the quasi-circular remnant. In addition, we assemble the OH spectral lines (1612/1665/1720 MHz) towards the complex by employing the newly released THOR (The HI OH Recombination line survey of Milky Way) data. We find that the known 1720 MHz OH maser in the W51B/C overlap area is located away from the detected 1612/1665MHz absorption region. The latter is sitting at the peak of the HII region G49.2-0.35 within W51B.
We use Green Bank Ammonia Survey observations of NH$_3$ (1,1) and (2,2) emission with 32'' FWHM resolution from a ~ 10 pc$^{2}$ portion of the Cepheus-L1251 molecular cloud to identify hierarchical dense gas structures. Our dendrogram analysis of the NH$_3$ data results in 22 top-level structures, which reside within 13 lower-level, parent structures. The structures are compact (0.01 pc $\lesssim R_{eff} \lesssim$ 0.1 pc) and are spatially correlated with the highest H$_2$ column density portions of the cloud. We also compare the ammonia data to a catalog of dense cores identified by higher-resolution (18.2'' FWHM) Herschel Space Observatory observations of dust continuum emission from Cepheus-L1251. Maps of kinetic gas temperature, velocity dispersion, and NH$_3$ column density, derived from detailed modeling of the NH$_3$ data, are used to investigate the stability and chemistry of the ammonia-identified and Herschel-identified structures. We show that the dust and dense gas in the structures have similar temperatures, with median $T_{dust}$ and $T_K$ measurements of 11.7 $\pm$ 1.1 K and 10.3 $\pm$ 2.0 K, respectively. Based on a virial analysis, we find that the ammonia-identified structures are gravitationally dominated, yet may be in or near a state of virial equilibrium. Meanwhile, the majority of the Herschel-identified dense cores appear to be not bound by their own gravity and instead confined by external pressure. CCS $(2_0-1_0)$ and HC$_5$N $(9-8)$ emission from the region reveal broader line widths and centroid velocity offsets when compared to the NH$_3$ (1,1) emission in some cases, likely due to these carbon-based molecules tracing the turbulent outer layers of the dense cores.
Interstellar extinction towards the Galactic Center is large and significantly differential. Its reddening and dimming effects in red clump stars in the Galactic Bulge can be exploited to better constrain the extinction law towards the innermost Galaxy. By virtue of a deep and complete catalog of more than 30 million objets at $|l|\le2.7\deg$ and $|b|\le1.55\deg$ obtained from VVV survey observations, we apply the red clump method to infer the selective-to-total extinction ratios in the $Z$, $Y$, $J$, $H$ and $K_s$ broadband near-infrared filters. The measured values are smaller than previously reported, and are not constant, with mean values, e.g., $A_{K_S}/E(J-K_s)=0.428\pm0.005\pm0.04$ and $A_{K_S}/E(H-K_s)=1.104\pm0.022\pm0.2$. We also obtain a ratio $A_Z$:$A_Y$:$A_J$:$A_H$:$A_{K_S}$ of 7.74:5.38:3.30:1.88:1.0, implying extinction towards the Galactic Center to follow a distribution as a function of wavelength steeper than previously reported, consistent with a power law $A_{\lambda}\propto{\lambda}^{-2.47}$ in the near-infrared.
The dynamic activity in massive star forming regions prior to the formation of bright protostars is still not fully investigated. In this work we present observations of HCO+ J=1-0 and N2H+ J=1-0 made with the IRAM 30m telescope towards a sample of 16 Herschel-identified massive 70 micron quiet clumps associated with infrared dark clouds. The clumps span a mass range from 300 M_sun to 2000 M_sun. The N2H+ data show that the regions have significant non-thermal motions with velocity dispersion between 0.28 km s^-1 and 1.5 km s^-1, corresponding to Mach numbers between 2.6 and 11.5. The majority of the 70 micron quiet clumps have asymmetric HCO+ line profiles, indicative of significant dynamical activity. We show that there is a correlation between the degree of line asymmetry and the surface density Sigma of the clumps, with clumps of Sigma>=0.1 g cm^-2 having more asymmetric line profiles, and so are more dynamically active, than clumps with lower Sigma. We explore the relationship between velocity dispersion, radius and Sigma and show how it can be interpreted as a relationship between an acceleration generated by the gravitational field a_G, and the measured kinetic acceleration, a_k, consistent with the majority of the non-thermal motions originating from self-gravity. Finally, we consider the role of external pressure and magnetic fields in the interplay of forces.
Spinning small silicate grains were recently invoked to account for the Galactic foreground anomalous microwave emission. These grains, if present, will absorb starlight in the far ultraviolet (UV). There is also renewed interest in attributing the enigmatic 2175 Angstrom interstellar extinction bump to small silicates. To probe the role of silicon in the UV extinction, we explore the relations between the amount of silicon required to be locked up in silicates [Si/H]_{dust} and the 2175 Angstrom bump or the far-UV extinction rise, based on an analysis of the extinction curves along 46 Galactic sightlines for which the gas-phase silicon abundance [Si/H]_{gas} is known. We derive [Si/H]_{dust} either from {[Si/H]_{ISM} - [Si/H]_{gas}} or from the Kramers-Kronig relation which relates the wavelength-integrated extinction to the total dust volume, where [Si/H]_{ISM} is the interstellar silicon reference abundance and taken to be that of proto-Sun or B stars. We also derive [Si/H]_{dust} from fitting the observed extinction curves with a mixture of amorphous silicates and graphitic grains. We find that in all three cases [Si/H]_{dust} shows no correlation with the 2175 Angstrom bump, while the carbon depletion [C/H]_{dust} tends to correlate with the 2175 Angstrom bump. This supports carbon grains instead of silicates as the possible carrier of the 2175 Angstrom bump. We also find that neither [Si/H]_{dust} nor [C/H]_{dust} alone correlates with the far-UV extinction, suggesting that the far-UV extinction is a combined effect of small carbon grains and silicates.
We revisited barium abundance in a sample of giant stars in the dwarf spheroidal Fornax galaxy. Non-local thermodynamic equilibrium (NLTE) synthesis was used to derive the barium abundance from two Ba II lines. Our new NLTE result for barium shows that in the range of metallicities from -1 to -0.5 the abundance of this element is almost the same as in the stars of the Milky Way. We conclude that the evolution of barium abundance in the dSph Fornax galaxy resembles that of the Milky Way at the epoch of the corresponding metallicity level.
Pulsar timing arrays (PTAs) are used to search for long-wavelength gravitational waves (GWs) by monitoring a set of spin-stable millisecond pulsars. Most theoretical analyses assume that the relevant GW sources are much more distant from Earth than the pulsars comprising the array. Unlike ground- or solar system-based GW detectors, PTAs might well contain embedded GW sources. We derive the PTA response from sources at any distance, with a specific focus on GW bursts with memory (BWMs). We consider supernovae and compact binary mergers as potential Galactic BWM sources and evaluate the signature for an array with pulsars in globular clusters or in the Galactic center. Understanding the response of PTAs to nearby sources of BWM is a step towards investigating other more complex Galactic sources.
We present a Census of the Local Universe (CLU) -- a combination of the largest area H$\alpha$ survey for nearby emission-line galaxies to date as well as a compilation of all publicly available galaxy catalogs out to a distance of 200 Mpc. With the Palomar 48-inch telescope, our H$\alpha$ survey images $\approx$3$\pi$ of the sky (3626 pointings) in 4 narrowband filters to an H$\alpha$ flux limit of $10^{-14}$ erg/s/cm$^2$ at 90% completeness. To characterize the completeness and contamination of the H$\alpha$ survey, we undertook a study of 14 preliminary fields (0.3% of total pointings) where we have obtained spectroscopic follow-up of all galaxy candidates with no previous distance information and a narrow-band color excess greater than 2.5$\sigma$. Comparison of the resulting CLU-H$\alpha$ galaxy catalog to a statistically complete sample of star-forming galaxies in a 11 Mpc volume (the LVL sample), we find that the H$\alpha$ survey would be complete to 85% in star-formation (using H$\alpha$ flux corrected for extinction as a proxy) and 70% in stellar mass at 200 Mpc. The contamination from high-redshift galaxies is lower than 8% (40%) for a 5$\sigma$ (2.5$\sigma$) color excess. Even with just 290 emission-line galaxies in the 14 preliminary fields, we find several interesting objects: 7 newly discovered blue compact dwarfs (aka, blueberries), 1 new green pea, 1 new QSO, and a known planetary nebula. The extreme galaxies (green pea and blueberries) have high star formation rates, low stellar masses, and low metallicities and could be local analogs to higher redshift primordial galaxies. The existence of these objects in our preliminary H$\alpha$ sample exemplifies that the full CLU-H$\alpha$ survey can be used as a discovery machine for a wide variety of objects in our own Galaxy and extreme galaxies out to intermediate redshifts.
Deep exposure imaging of early-type galaxies (ETGs) are revealing the second-order complexity of these objects, which have been long considered uniform, dispersion-supported spheroidals. "Fine structure" features (e.g. ripples, plumes, tidal tails, rings) as well as depleted stellar cores (i.e. central light deficits) characterize a number of massive ETG galaxies, and can be interpreted as the result of galaxy-galaxy interactions. We discuss how the timescale for the evolution of cores and fine structures are comparable, and hence it is expected that they develop in parallel after the major interaction event which shaped the ETG. Using archival data, we compare the "depleted stellar mass" (i.e. the mass missing from the depleted stellar core) against the prominence of the fine structure features, and observe that they correlate inversely. This result confirms our expectation that, while the Super Massive Black Hole (SMBH) binary (constituted by the SMBHs of the merger progenitors) excavates the core via three-body interactions, the gravitational potential of the newborn galaxy relaxes, and the fine structures fade below detection levels. We expect the inverse correlation to hold at least within the first Gyr from the merger which created the SMBH binary; after then, the fine structure evolves independently.
We estimate the population of black holes with luminous stellar companions (BH-LCs) in the Milky Way (MW) observable by Gaia. We evolve a realistic distribution of BH-LC progenitors from zero-age to the current epoch taking into account relevant physics, including binary stellar evolution, BH-formation physics, and star formation rate, to estimate the BH-LC population in the MW today. We predict that Gaia will discover between 3,800 and 12,000 BH-LCs by the end of its 5 yr mission, depending on BH natal kick strength and observability constraints. We find that the overall yield, and distributions of eccentricities and masses of observed BH-LCs can provide important constraints on the strength of BH natal kicks. Gaia-detected BH-LCs are expected to have very different orbital properties compared to those detectable via radio, X-ray, or gravitational wave observations.
The mechanism for producing polarized emission from protostellar disks at (sub)millimeter wavelengths is currently uncertain. Classically, polarization is expected from non-spherical grains aligned with the magnetic field. Recently, two alternatives have been suggested. One polarization mechanism is caused by self-scattering from dust grains of sizes comparable to the wavelength while the other mechanism is due to grains aligned with their short axes along the direction of radiation anisotropy. The latter has recently been shown as a likely mechanism for causing the dust polarization detected in HL Tau at 3 mm. In this Letter, we present ALMA polarization observations of HL Tau for two more wavelengths: 870 $\mu$m and 1.3 mm. The morphology at 870 $\mu$m matches the expectation for self-scattering, while that at 1.3 mm shows a mix between self-scattering and grains aligned with the radiation anisotropy. The observations cast doubt on the ability of (sub)millimeter continuum polarization to probe disk magnetic fields for at least HL Tau. By showing two distinct polarization morphologies at 870 $\mu$m and 3 mm and a transition between the two at 1.3 mm, this Letter provides definitive evidence for two disk polarization mechanisms operating at (sub)millimeter wavelengths. In addition, if the polarization at 870 $\mu$m is due to scattering, the lack of polarization asymmetry along the minor axis of the inclined disk implies that the large grains responsible for the scattering have already settled into a geometrically thin layer, and the presence of asymmetry along the major axis indicates that the HL Tau disk is not completely axisymmetric.
From our recent catalog based on the first Gaia data release (TGAS), we select wide binaries in which both stars have been observed by the Radial Velocity Experiment (RAVE) or the Large Sky Area Multi-Object Fibre Spectroscopic Telescope (LAMOST). Using RAVE and LAMOST metallicities and RAVE Mg, Al, Si, Ti, and Fe abundances, we find that the differences in the metallicities and elemental abundances of components of wide binaries are consistent with being due to observational uncertainties, in agreement with previous results for smaller and more restricted samples. The metallicity and elemental abundance consistency between wide binary components presented in this work confirms their common origin and bolsters the status of wide binaries as "mini-open clusters." Furthermore, this is evidence that wide binaries are effectively co-eval and co-chemical, supporting their use for e.g., constraining age-activity-rotation relations, the initial-final mass relation for white dwarfs, and M-dwarf metallicity indicators. Additionally, we demonstrate that the common proper motion, common parallax pairs in TGAS with the most extreme separations (s $\gtrsim$ 0.1 pc) typically have inconsistent metallicities, radial velocities or both and are therefore likely to be predominantly comprised of random alignments of unassociated stars with similar astrometry, in agreement with our previous results. Finally, we propose that wide binaries form an ideal data set with which to test chemical tagging as a method to identify stars of common origin, particularly because the stars in wide binaries span a wide range of metallicities, much wider than that spanned by nearby open clusters.
It has recently been proposed that massive primordial black holes (PBH) could constitute all of the dark matter, providing a novel scenario of structure formation, with early reionization and a rapid growth of the massive black holes at the center of galaxies and dark matter halos. The scenario arises from broad peaks in the primordial power spectrum that give both a spatially clustered and an extended mass distribution of PBH. The constraints from the observed microlensing events on the extended mass function have already been addressed. Here we study the impact of spatial clustering on the microlensing constraints. We find that the bounds can be relaxed significantly for relatively broad mass distributions if the number of primordial black holes within each cluster is typically above one hundred. On the other hand, even if they arise from individual black holes within the cluster, the bounds from CMB anisotropies are less stringent due to the enhanced black hole velocity in such dense clusters. This way, the window between a few and ten solar masses has opened up for PBH to comprise the totality of the dark matter.
Recently there has been much interest in the spatial distribution of light scalar dark matter, especially axions, throughout the universe. When the local gravitational interactions between the scalar modes are sufficiently rapid, it can cause the field to re-organize into a BEC of gravitationally bound clumps. While these clumps are stable when only gravitation is included, the picture is complicated by the presence of the axion's attractive self-interactions, which can potentially cause the clumps to collapse. Here we perform a detailed stability analysis to determine under what conditions the clumps are stable. In this paper we focus on spherical configurations, leaving aspherical configurations for future work. We identify branches of clump solutions of the axion-gravity-self-interacting system and study their stability properties. We find that clumps that are (spatially) large are stable, while clumps that are (spatially) small are unstable and may collapse. Furthermore, there is a maximum number of particles that can be in a clump. We map out the full space of solutions, which includes quasi-stable axitons, and clarify how a recent claim in the literature of a new ultra-dense branch of stable solutions rests on an invalid use of the non-relativistic approximation. We also consider repulsive self-interactions that may arise from a generic scalar dark matter candidate, finding a single stable branch that extends to arbitrary particle number.
Though theoretically expected, the charge exchange emission from galaxy clusters has not yet been confidently detected. Accumulating hints were reported recently, including a rather marginal detection with the Hitomi data of the Perseus cluster. As suggested in Gu et al. (2015), a detection of charge exchange line emission from galaxy clusters would not only impact the interpretation of the newly-discovered 3.5 keV line, but also open up a new research topic on the interaction between hot and cold matter in clusters. We aim to perform the most systematic search for the O VIII charge exchange line in cluster spectra using the RGS on board XMM. We introduce a sample of 21 clusters observed with the RGS. The dominating thermal plasma emission is modeled and subtracted with a two-temperature CIE component, and the residuals are stacked for the line search. The systematic uncertainties in the fits are quantified by refitting the spectra with a varying continuum and line broadening. By the residual stacking, we do find a hint of a line-like feature at 14.82 A, the characteristic wavelength expected for oxygen charge exchange. This feature has a marginal significance of 2.8 sigma, and the average equivalent width is 2.5E-4 keV. We further demonstrate that the putative feature can be hardly affected by the systematic errors from continuum modelling and instrumental effects, or the atomic uncertainties of the neighbouring thermal lines. Assuming a realistic temperature and abundance pattern, the physical model implied by the possible oxygen line agrees well with the theoretical model proposed previously to explain the reported 3.5 keV line. If the charge exchange source indeed exists, we would expect that the oxygen abundance is potentially overestimated by 8-22% in previous X-ray measurements which assumed pure thermal lines.
The matching of sources between photometric catalogues can lead to cases where objects of differing brightness are incorrectly assumed to be detections of the same source. The rejection of unphysical matches can be achieved through the inclusion of information about the sources' magnitudes. The method described here uses the additional photometric information from both catalogues in the process of accepting or rejecting counterparts, providing approximately a factor 10 improvement in Bayes' factor with its inclusion. When folding in the photometric information we avoid using prior astrophysical knowledge. Additionally, the method allows for the possibility of no counterparts to sources as well as the possibility that sources overlap multiple potential counterparts. We formally describe the probability of two sources being the same astrometric object, allowing systematic effects of astrometric perturbation (by, e.g., contaminant objects) to be accounted for. We apply the method to two cases. First, we test IPHAS-Gaia matches to compare the resulting matches in two catalogues of similar wavelength coverage but differing dynamical ranges. Second, we apply the method to matches between IPHAS and 2MASS and show that the method holds when considering two catalogues with approximately equal astrometric precision. We discuss the importance of including the magnitude information in each case. Additionally, we discuss extending the method to multiple catalogue matches through an iterative matching process. The method allows for the selection of high-quality matches by providing an overall probability for each pairing, giving the flexibility to choose stars known to be good matches.
We present a blind search for doublet intergalactic metal absorption with a method we dub `agnostic stacking'. Using forward-modelling we combine it with direct detections in the literature to measure the overall metal population. Here we apply this novel approach to the search for NeVIII in 26 high-quality COS spectra of QSOs at z>0.7. We probe an unprecedented low limit of log N>12.3 at 0.47<z<1.34 with a total pathlength ${\Delta}$z = 7.36. The method selects absorption without requiring knowledge of its source, be it observing noise, artifacts, or any line transition. Stacking this mixed population with NeVIII absorption dilutes doublet features in composite spectra in a deterministic manner. We stack potential NeVIII absorption in two regimes: absorption too weak to be statistically significant in direct line studies (12.3 < log N< 13.7), and strong absorbers (log N> 13.7). We do not detect NeVIII in either regime, and place upper limits on the population using agnostic stacking alone. Combining our measurements with direct line detections, the NeVIII population is reproduced with a single power law column density distribution of slope \b{eta} = -1.86 and normalisation log f_{13.7} = -13.99, leading to an incidence rate of strong NeVIII absorbers of dn/dz =1.38. Comparing our results with a group of 3 systems in PG1148+549, these have a 0.024% probability of arising by chance. We infer a cosmic mass density for NeVIII in the column density range 12.3 < log N < 15.0 of ${\Omega}$(NeVIII) = 2.2x10^{-8}. We translate this inferred density into an estimate of the baryon density of the NeVIII-bearing gas, and arrive at ${\Omega}b~1.8x10^{-3}$, which constitutes only 4% of the total baryonic mass. The measured NeVIII column density distribution function and cosmic density here are inconsistent with predictions of the EAGLES simulations at ${\sigma}>2.0$ significance. (abridged)
Links to: arXiv, form interface, find, astro-ph, recent, 1710, contact, help (Access key information)
We study radial profiles in H$\alpha$ equivalent width and specific star formation rate (sSFR) derived from spatially-resolved SDSS-IV MaNGA spectroscopy to gain insight on the physical mechanisms that suppress star formation and determine a galaxy's location in the SFR-$\rm M_\star$ diagram. Even within the star-forming blue cloud, the typical sSFR profile depends on stellar mass. Flat radial profiles are observed for $\rm log(M_\star/M_\odot ) < 10.5$, while star-forming galaxies of higher mass show a significant decrease in sSFR in the central regions, a likely consequence of both larger bulges and an inside-out growth history. Our primary focus is the green valley where, at all masses, we find sSFR profiles that are suppressed with respect to normal star-forming galaxies at all galactocentric distances out to 2 effective radii. The responsible quenching mechanism therefore appears to affect the entire galaxy, not simply an expanding central region. However, those galaxies in which central star formation has shut down (classified spectroscopically as central low-ionisation emission-line regions, or cLIERs) show an even stronger suppression in sSFR across the whole disc. These systems are also more suppressed than green valley galaxies with residual, central star-formation. In fact, compared to structural parameters like $\Sigma_1$ (the mass surface density within 1 kpc), central quiescence is a stronger predictor of ongoing sSFR suppression at fixed $\rm M_\star$. The green valley hosts both quiescent bulges and strongly suppressed star forming discs, supporting a scenario in which a slow quenching process affects the entire galaxy, and not just the central regions.
We study the spatially-resolved star formation rates of 1494 galaxies in the SDSSIV-MaNGA Survey. Star formation rates are calculated using a two-step process, exploiting H$_\alpha$ in star-forming regions and $D_n4000$ in zones identified as AGN/LI(N)ER or lineless according to the BPT diagram. This enables a complete view of how radial profiles of specific star formation rate depend on stellar mass, galaxy structure and environment, providing valuable clues about the nature of both secular and environmental quenching processes. We report on the existence of `Centrally Suppressed' galaxies, which have lower SSFR in their cores compared to their disks. The profiles of centrally suppressed and unsuppressed galaxies are distributed in a bimodal way, with very few galaxies having intermediate profile shapes. Galaxies with large stellar masses and central velocity dispersions are much more likely to be centrally-suppressed than low-mass galaxies, a result linked to their morphology and the presence of AGN/LI(N)ER like emission. Centrally suppressed galaxies also display lower star formation at all radii compared to unsuppressed galaxies. Independent of galaxy type, galaxies identified as satellites in larger dark matter halos have noticeably lower SSFR at outer radii compared to central galaxies. This uniform suppression could be a signal of the stripping of hot gas reservoirs through a process like strangulation. We find that satellites are not more likely to be suppressed in their cores than centrals, indicating that core suppression is an entirely secular process. We find no secondary correlation between the local environment density and the profiles of star formation rate surface density.
We tested the performance of photometric redshifts for galaxies in the Hubble Ultra Deep field down to 30th magnitude. We compared photometric redshift estimates from three spectral fitting codes from the literature (EAZY, BPZ and BEAGLE) to high quality redshifts for 1227 galaxies from the MUSE integral field spectrograph. All these codes can return photometric redshifts with bias |Dzn|=|z-z_phot|/(1+z)<0.05 down to F775W=30 and spectroscopic incompleteness is unlikely to strongly modify this statement. We have, however, identified clear systematic biases in the determination of photometric redshifts: in the 0.4<z<1.5 range, photometric redshifts are systematically biased low by as much as Dzn=-0.04 in the median, and at z>3 they are systematically biased high by up to Dzn = 0.05, an offset that can in part be explained by adjusting the amount of intergalactic absorption applied. In agreement with previous studies we find little difference in the performance of the different codes, but in contrast to those we find that adding extensive ground-based and IRAC photometry actually can worsen photo-z performance for faint galaxies. We find an outlier fraction, defined through |Dzn|>0.15, of 8% for BPZ and 10% for EAZY and BEAGLE, and show explicitly that this is a strong function of magnitude. While this outlier fraction is high relative to numbers presented in the literature for brighter galaxies, they are very comparable to literature results when the depth of the data is taken into account. Finally, we demonstrate that while a redshift might be of high confidence, the association of a spectrum to the photometric object can be very uncertain and lead to a contamination of a few percent in spectroscopic training samples that do not show up as catastrophic outliers, a problem that must be tackled in order to have sufficiently accurate photometric redshifts for future cosmological surveys.
The possible detection of C_{24}, a planar graphene, recently reported in several planetary nebulae by Garciaa-Hernandez et al. (2011, 2012) inspires us to explore whether and how much graphene could exist in the interstellar medium (ISM) and how it would reveal its presence through its ultraviolet (UV) extinction and infrared (IR) emission. In principle, interstellar graphene could arise from the photochemical processing of polycyclic aromatic hydrocarbon (PAH) molecules which are abundant in the ISM through a complete loss of their hydrogen atoms and/or from graphite which is thought to be a major dust species in the ISM through fragmentation caused by grain-grain collisional shattering. Both quantum-chemical computations and laboratory experiments have shown that the exciton-dominated electronic transitions in graphene cause a strong absorption band near 2755 Angstrom. We calculate the UV absorption of graphene and place an upper limit of ~5 ppm of C/H (i.e., ~1.9\% of the total interstellar C) on the interstellar graphene abundance. We also model the stochastic heating of graphene C_{24} excited by single starlight photons of the interstellar radiation field in the ISM and calculate its IR emission spectra. We also derive the abundance of graphene in the ISM to be <5 ppm of C/H by comparing the model emission spectra with that observed in the ISM.
We construct average spectra of host galaxies of slower, faster, bluer, and redder Type Ia Supernovae (SNe Ia) from the SDSS-II supernova survey. The average spectrum of slower declining (broader light-curve width or higher stretch) SN Ia hosts shows stronger emission lines compared to the average spectrum of faster declining (narrower light-curve width or lower stretch) SN Ia hosts. Using pPXF, we find that hosts of slower declining SNe Ia have metallicities that are, on average, 0.24 dex higher than average metallicities of faster declining SN Ia hosts. Similarly, redder SN Ia hosts have slightly higher metallicities than bluer SN Ia hosts. Lick index analysis of metallic lines and Balmer lines show that faster declining SN Ia hosts have relatively higher metal content and have relatively older stellar populations compared with slower declining SN Ia hosts. We calculate average $\rm H_{\alpha}$ Star Formation Rate (SFR), stellar mass, and the specific-SFR (sSFR) of host galaxies in these subgroups of SNe Ia. We find that slower declining SN Ia hosts have significantly higher ($>5\sigma$) sSFR than faster declining SN Ia hosts. A Kolmogorov-Smirnov test shows that these two types of hosts originate from different parent distributions. Our results, when compared with the models of \cite{childress14}, indicate that slower declining SNe Ia, being hosted in actively star-forming galaxies, are young (prompt) SNe Ia, originating from similar progenitor age groups.
The intrinsic shape of galactic bulges in nearby galaxies provides crucial information to separate bulge types. We intended to derive accurate constraints to the intrinsic shape of bulges to provide new clues on their formation mechanisms and set new limitations for future simulations. We retrieved the intrinsic shape of a sample of CALIFA bulges using a statistical approach. Taking advantage of GalMer numerical simulations of binary mergers we estimated the reliability of the procedure. Analyzing the $i$-band mock images of resulting lenticular remnants, we studied the intrinsic shape of their bulges at different galaxy inclinations. Finally, we introduced a new ($B/A$, $C/A$) diagram to analyze possible correlations between the intrinsic shape and the properties of bulges. We tested the method on simulated lenticular remnants, finding that for galaxies with inclinations $25^{\circ} < \theta < 65^{\circ}$ we can safely derive the intrinsic shape of their bulges. We found that our CALIFA bulges tend to be nearly oblate systems (66%), with a smaller fraction of prolate spheroids (19%) and triaxial ellipsoids (15%). The majority of triaxial bulges are in barred galaxies (75%). Moreover, we found that bulges with low S\'ersic indices or in galaxies with low bulge-to-total luminosity ratios form a heterogeneous class of objects; additionally, also bulges in late-type galaxies or in less massive galaxies have no preference in being oblate, prolate, or triaxial. On the contrary, bulges with high S\'ersic index, in early-type galaxies, or in more massive galaxies are mostly oblate systems. We concluded that various evolutionary pathways may coexist in galaxies, with merging events and dissipative collapse being the main mechanisms driving the formation of the most massive oblate bulges and bar evolution reshaping the less massive triaxial bulges.
We study a sample of 17 $z>1.5$ C~{\sc i} absorbers observed with the ESO-VLT spectrograph X-SHOOTER. We derive metallicities, depletion onto dust, extinction by dust and analyse the absorption from Mg~{\sc ii}, Mg~{\sc i}, Ca~{\sc ii} and Na~{\sc i} that are redshifted into the near infrared wavelength range. We show that most of them have high metallicity and dust content. We detect 9 Ca~{\sc ii} absorptions with $W$(Ca~{\sc ii}$\lambda$3934) ~$>$~0.23~\AA~ 9 out of 14 systems where we have appropriate wavelength coverage. The observed equivalent widths are similar to what has been measured in other lower redshift surveys of Ca~{\sc ii} systems. We detect 10 Na~{\sc i} absorptions in the 11 systems where we could observe this absorption. The median equivalent width ($W$(Na~{\sc i}$\lambda$5891)~=~0.68~\AA) is larger than what is observed in local clouds with similar H~{\sc i} column densities but also in $z<0.7$ Ca~{\sc ii} systems detected in the SDSS. Most of the systems (12 out of 17) have $W$(Mg~{\sc ii}$\lambda$2796)~$>$~2.5~\AA~ when six of them have log~$N$(H~{\sc i})~$<$~20.3 with the extreme case of J1341+1852 which has log~$N$(H~{\sc i})~=~18.18. The Mg~{\sc ii} absorptions are spread over more than $\Delta v$ $\sim$ 400~km~s$^{-1}$ for half of the systems; three absorbers have $\Delta v$~$>$~500~km~s$^{-1}$. The kinematics is strongly perturbed for most of these systems which cannot arise in quiet disks and must be close to regions with intense star formation activity, contrary to what is usually observed for strong $W > 2$~\AA~ Mg~{\sc ii} systems and/or are part of objects in interaction. All this suggests that a large fraction of the cold gas at high redshift arises in disturbed environments.
We report on interferometric observations of a face-on accretion system around the High-Mass young stellar object, G353.273+0.641. The innermost accretion system of 100 au radius was resolved in a 45 GHz continuum image taken with the Jansky-Very Large Array. Our spectral energy distribution analysis indicated that the continuum could be explained by optically thick dust emission. The total mass of the dusty system is $\sim$ 0.2 $M_{\sun}$ at minimum and up to a few $M_{\sun}$ depending on the dust parameters. 6.7 GHz CH$_{3}$OH masers associated with the same system were also observed with the Australia Telescope Compact Array. The masers showed a spiral-like, non-axisymmetric distribution with a systematic velocity gradient. The line-of-sight velocity field is explained by an infall motion along a parabolic streamline that falls onto the equatorial plane of the face-on system. The streamline is quasi-radial and reaches the equatorial plane at a radius of 16 au. This is clearly smaller than that of typical accretion disks in High-Mass star formation, indicating that the initial angular momentum was very small, or the CH$_{3}$OH masers selectively trace accreting material that has small angular momentum. In the former case, the initial specific angular momentum is estimated to be 8 $\times$ 10$^{20}$ ($M_{*}$$/$10 $M_{\sun}$)$^{0.5}$ cm$^{2}$ s$^{-1}$, or a significant fraction of the initial angular momentum was removed outside of 100 au. The physical origin of such a streamline is still an open question and will be constrained by the higher-resolution ($\sim$ 10 mas) thermal continuum and line observations with ALMA long baselines.
Simulations generally show that non-self-gravitating clouds have a lognormal column density ($\Sigma$) probability distribution function (PDF), while self-gravitating clouds with active star formation develop a distinct power-law tail at high column density. Although the growth of the power law can be attributed to gravitational contraction leading to the formation of condensed cores, it is often debated if an observed lognormal shape is a direct consequence of supersonic turbulence alone, or even if it is really observed in molecular clouds. In this paper we run three-dimensional magnetohydrodynamic simulations including ambipolar diffusion with different initial conditions to see the effect of strong magnetic fields and nonlinear initial velocity perturbations on the evolution of the column density PDFs. Our simulations show that column density PDFs of clouds with supercritical mass-to-flux ratio, with either linear perturbations or nonlinear turbulence, quickly develop a power-law tail such that $dN/d \log \Sigma \propto \Sigma^{-\alpha}$ with index $\alpha \simeq 2$. Interestingly, clouds with subcritical mass-to-flux ratio also proceed directly to a power-law PDF, but with a much steeper index $\alpha \simeq 4$. This is a result of gravitationally-driven ambipolar diffusion. However, for nonlinear perturbations with a turbulent spectrum ($v_{k}^{2} \propto k^{-4}$), the column density PDFs of subcritical clouds do retain a lognormal shape for a major part of the cloud evolution, and only develop a distinct power-law tail with index $\alpha \simeq 2$ at greater column density when supercritical pockets are formed.
We present an analysis of the host-galaxy environment of Swope Supernova Survey 2017a (SSS17a), the discovery of an electromagnetic counterpart to a gravitational wave source, GW170817. SSS17a occurred 1.9 kpc (in projection; 10.2") from the nucleus of NGC 4993, an S0 galaxy at a distance of 40 Mpc. We present a Hubble Space Telescope (HST) pre-trigger image of NGC 4993, Magellan optical spectroscopy of the nucleus of NGC 4993 and the location of SSS17a, and broad-band UV through IR photometry of NGC 4993. The spectrum and broad-band spectral-energy distribution indicate that NGC 4993 has a stellar mass of log (M/M_solar) = 10.49^{+0.08}_{-0.20} and star formation rate of 0.003 M_solar/yr, and the progenitor system of SSS17a likely had an age of >2.8 Gyr. There is no counterpart at the position of SSS17a in the HST pre-trigger image, indicating that the progenitor system had an absolute magnitude M_V > -5.8 mag. We detect dust lanes extending out to almost the position of SSS17a and >100 likely globular clusters associated with NGC 4993. The offset of SSS17a is similar to many short gamma-ray burst offsets, and its progenitor system was likely bound to NGC 4993. The environment of SSS17a is consistent with an old progenitor system such as a binary neutron star system.
This is the first in a series of papers examining the demographics of star-forming galaxies at $0.2<z<2.5$ in CANDELS. We study 9,100 galaxies from GOODS-S and UDS having published values of redshifts, masses, star-formation rates (SFRs), and dust attenuation ($A_V$) derived from UV-optical SED fitting. In agreement with previous works, we find that the $UVJ$ colors of a galaxy are closely correlated with its specific star-formation rate (SSFR) and $A_V$. We define rotated $UVJ$ coordinate axes, termed $S_\mathrm{SED}$ and $C_\mathrm{SED}$, that are parallel and perpendicular to the star-forming sequence and derive a quantitative calibration that predicts SSFR from $C_\mathrm{SED}$ with an accuracy of ~0.2 dex. SFRs from UV-optical fitting and from UV+IR values based on Spitzer/MIPS 24 $\mu\mathrm{m}$ agree well overall, but systematic differences of order 0.2 dex exist at high and low redshifts. A novel plotting scheme conveys the evolution of multiple galaxy properties simultaneously, and dust growth, as well as star-formation decline and quenching, exhibit "mass-accelerated evolution" ("downsizing"). A population of transition galaxies below the star-forming main sequence is identified. These objects are located between star-forming and quiescent galaxies in $UVJ$ space and have lower $A_V$ and smaller radii than galaxies on the main sequence. Their properties are consistent with their being in transit between the two regions. The relative numbers of quenched, transition, and star-forming galaxies are given as a function of mass and redshift.
Star-forming Compact Dwarf Galaxies (CDGs) resemble the expected pristine conditions of the first galaxies in the Universe and are the best systems to test models on primordial galaxy formation and evolution. Here we report on one of such CDGs, Tololo 1214-277, which presents a broad, single peaked, highly symmetric Ly$\alpha$ emission line that had evaded theoretical interpretation so far. In this paper we reproduce for the first time these line features with two different physically motivated kinematic models: an interstellar medium composed by outflowing clumps with random motions and an homogeneous gaseous sphere undergoing solid body rotation. The multiphase model requires a clump velocity dispersion of $54.3\pm 0.6$ km s$^{-1}$ with outflows of $54.3\pm 5.1$ km s$^{-1}$, while the bulk rotation velocity is constrained to be $348^{+75}_{-48}$ km s$^{-1}$. We argue that the results from the multiphase model provide a correct interpretation of the data. In that case the clump velocity dispersion implies a dynamical mass of $2\times 10^{9}$ M$_{\odot}$, ten times its baryonic mass. If future kinematic maps of Tololo 1214-277 confirm the velocities suggested by the multiphase model, it would provide additional support to expect such kinematic state in primordial galaxies, opening the opportunity to use the models and methods presented in this paper to constrain the physics of star formation and feedback in the early generation of Ly-$\alpha$ emitting galaxies.
In Debattista et al. (2015), we proposed that a kiloparsec-scale nuclear disc is responsible for the high-velocity secondary peak in the stellar line-of-sight velocity distributions (LOSVDs) seen at positive longitudes in the bulge by the Apache Point Observatory Galactic Evolution Experiment (APOGEE). Here, we make further qualitative but distinctive predictions of the kinematic properties of a nuclear disc, including for the LOSVDs at negative longitudes (which APOGEE-2 will observe) and examine the proper motions throughout the disc. Since a nuclear ring is also able to produce similar high-velocity LOSVD peaks, we present predictions for the proper motion signatures which distinguish between a nuclear disc and a nuclear ring. We also demonstrate that the stars in a nuclear disc, which would be on x2 orbits perpendicular to the bar, can remain on these orbits for a long time and can therefore be old. We show that such (old) nuclear discs of comparable size exist in external galaxies.
We present the clustering properties and halo occupation distribution (HOD) modeling of very low redshift, hard X-ray-detected active galactic nuclei (AGN) using cross-correlation function measurements with 2MASS galaxies. Spanning a redshift range of $0.007 < z < 0.037$, with a median $z=0.024$, we present a precise AGN clustering study of the most local AGN in the Universe. The AGN sample is drawn from the SWIFT/BAT 70-Month and INTEGRAL/IBIS 8-Year all-sky X-ray surveys and contains both type I and type II AGN. We find a large-scale bias for the full AGN sample of $b=1.04^{+0.10}_{-0.11}$, which corresponds to a typical host dark matter halo mass of $M_{\rm h}^{\rm typ}=12.84^{+0.22}_{-0.30}\,h^{-1} M_{\odot}$. When split into low and high X-ray luminosity and type I and type II AGN subsamples, we detect no statistically significant differences in the large-scale bias parameters. However, there are differences in the small-scale clustering which are reflected in the full HOD model results. We find that low and high X-ray luminosity AGN, as well as type I and type II AGN, occupy dark matter halos differently, with 3.4$\sigma$ and 4.0$\sigma$ differences in their mean halo masses, respectively, when split by luminosity and type. The latter finding contradicts a simple orientation-based AGN unification model. As a by-product of our cross-correlation approach, we also present the first HOD model of 2MASS galaxies.
To investigate the $M_\bullet -\sigma$ relation, we consider realistic elliptical galaxy profiles that are taken to follow a single power law density profile given by $\rho(r) = \rho_{0}(r/ r_{0})^{-\gamma}$ or the Nuker intensity profile. We calculate the density using Abel's formula in the latter case by employing the derived stellar potential in both cases, we derive the distribution function $f(E)$ of the stars in presence of the supermassive black hole (SMBH) at the center and hence compute the line of sight (LOS) velocity dispersion as a function of radius. For the typical range of values for masses of SMBH, we obtain $M_{\bullet} \propto \sigma^{p}$ for different profiles. An analytical relation $p = (2\gamma + 6)/(2 + \gamma)$ is found which is in reasonable agreement with observations (for $\gamma$ = 0.75 - 1.4, $p$ = 3.6 - 5.3). Assuming that a proportionality relation holds between the black hole mass and bulge mass, $ M_{\bullet} =f M_b$, and applying this to several galaxies we find the individual best fit values of $p$ as a function of $f$; also by minimizing $\chi^{2}$, we find the best fit global $p$ and $f$. For Nuker profiles we find that $p$ = $3.81 \pm 0.004$ and $f$ = $(1.23 \pm 0.09)\times 10^{-3}$ which are consistent with the observed ranges.
Observations show that massive star formation quenches first at centers of galaxies. To understand quenching mechanisms, we investigate the thermal and nonthermal energy balance in the central kpc of NGC1097- a prototypical galaxy undergoing quenching- and present a systematic study of the nuclear star formation efficiency and its dependencies. This region is dominated by the nonthermal pressure from the magnetic field, cosmic rays, and turbulence. A comparison of the mass-to-magnetic flux ratio of the molecular clouds shows that most of them are magnetically critical or supported against gravitational collapse needed to form cores of massive stars. Moreover, the star formation efficiency of the clouds drops with the magnetic field strength. Such an anti-correlation holds with neither the turbulent nor the thermal pressure. Hence, a progressive built up of the magnetic field results in high-mass stars forming inefficiently, and it may be the cause of the low-mass stellar population in the bulges of galaxies.
New observations performed at IRAM Plateau de Bure reveal the absence of molecular gas next to Andromeda's black hole. We derived a 3 sigma upper limit on the molecular gas mass of 4300 Msol for the linewidth of 1000 km/s. This is compatible with infra-red observations which reveal a hole in dust emission next to the black hole. Some gas from stellar feedback is expected from the old eccentric stellar disc population, but it is not accreted close to the black hole. This absence of gas explains the absence of stellar formation observed in this region contrary to what is observed next to Sgr A* in the Milky Way. Either the gas has been swallowed by the black hole, or a feedback mechanism has pushed the gas outside the central 1 pc. Nevertheless, we detect a small clump of gas with a very small velocity dispersion at 2.4" from the black hole. It is probable that this clumpy gas is seen in projection, as it does not follow the rotation of the disk surrounding the black hole, its velocity dispersion is ten times smaller than the expected velocity gradient and the tidal shear from the black hole requires a gas density for this clump that is not compatible with our observations.
(Abridged) Context: The assumption of a gas-to-dust mass ratio (\gamma) is a
common approach to estimate the basic properties of molecular clouds, such as
total mass and column density of molecular hydrogen, from (sub)mm continuum
observations of the dust. In the Milky Way a single value is used at all
galactocentric radii, independently of the observed metallicity gradients. Both
models and extragalactic observations suggest that this quantity increases for
decreasing metallicity Z, typical of the outer regions in disks, where fewer
heavy elements are available to form dust grains.
Aims: We aim to investigate the variation of the gas-to-dust ratio as a
function of galactocentric radius and metallicity, to allow a more accurate
characterisation of the quantity of molecular gas across the galactic disk, as
derived from observations of the dust.
Methods: Observations of the optically thin C\$^{18}\$O (2-1) transition were
obtained with the APEX telescope for a sample of 23 massive and dense
star-forming regions in the far outer Galaxy (galactocentric distance greater
than 14 kpc). From the modelling of this line and of the spectral energy
distribution of the selected clumps we computed the gas-to-dust ratio and
compared it to that of well-studied sources from the ATLASGAL TOP100 sample in
the inner galactic disk.
Results: The gradient in gas-to-dust ratio is found to be 0.087 dex/kpc (or
equivalently \gamma\ proportional to Z\$^{-1.4}\$). The dust-to-metal ratio,
decreases with galactocentric radius, which is the most common situation also
for external late-type galaxies. This suggests that grain growth dominates over
destruction. The predicted gas-to-dust ratio is in excellent agreement with the
estimates in Magellanic clouds, for the appropriate value of Z.
We present the analysis of near infrared (NIR), adaptive optics (AO) Subaru and archived HST imaging data of a region near the northern middle lobe (NML) of the Centaurus A (Cen A) jet, at a distance of $\sim15$ kpc north-east (NE) from the center of NGC5128. Low-pass filtering of the NIR images reveals strong -- $>3\sigma$ above the background mean -- signal at the expected position of the brightest star in the equivalent HST field. Statistical analysis of the NIR background noise suggests that the probability to observe $>3\sigma$ signal at the same position, in three independent measurements due to stochastic background fluctuations alone is negligible ($\leq10^{-7}\%$) and, therefore, that this signal should reflect the detection of the NIR counterparts of the brightest HST star. An extensive photometric analysis of this star yields $V-I$, visual-NIR, and NIR colors expected from a yellow supergiant (YSG) with an estimated age $\sim10^{+4}_{-3}$ Myr. Furthermore, the second and third brighter HST stars are, likely, also supergiants in Cen A, with estimated ages $\sim16^{+6}_{-3}$ Myr and $\sim25^{+15}_{-9}$ Myr, respectively. The ages of these three supergiants are in good agreement with the ages of the young massive stars that were previously found in the vicinity and are thought to have formed during the later phases of the jet-HI cloud interaction that appears to drive the star formation (SF) in the region for the past $\sim100$ Myr.
The historic detection of gravitational waves from a binary neutron star merger (GW170817) and its electromagnetic counterpart led to the first accurate (sub-arcsecond) localization of a gravitational-wave event. The transient was found to be $\sim$10" from the nucleus of the S0 galaxy NGC 4993. We report here the luminosity distance to this galaxy using two independent methods. (1) Based on our MUSE/VLT measurement of the heliocentric redshift ($z_{\rm helio}=0.009783\pm0.000023$) we infer the systemic recession velocity of the NGC 4993 group of galaxies in the cosmic microwave background (CMB) frame to be $v_{\rm CMB}=3231 \pm 53$ km s$^{-1}$. Using constrained cosmological simulations we estimate the line-of-sight peculiar velocity to be $v_{\rm pec}=307 \pm 230$ km s$^{-1}$, resulting in a cosmic velocity of $v_{\rm cosmic}=2924 \pm 236$ km s$^{-1}$ ($z_{\rm cosmic}=0.00980\pm 0.00079$) and a distance of $D_z=40.4\pm 3.4$ Mpc assuming a local Hubble constant of $H_0=73.24\pm 1.74$ km s$^{-1}$ Mpc$^{-1}$. (2) Using Hubble Space Telescope measurements of the effective radius (15.5" $\pm$ 1.5") and contained intensity and MUSE/VLT measurements of the velocity dispersion, we place NGC 4993 on the Fundamental Plane (FP) of E and S0 galaxies. Comparing to a frame of 10 clusters containing 226 galaxies, this yields a distance estimate of $D_{\rm FP}=44.0\pm 7.5$ Mpc. The combined redshift and FP distance is $D_{\rm NGC 4993}= 41.0\pm 3.1$ Mpc. This 'electromagnetic' distance estimate is consistent with the independent measurement of the distance to GW170817 as obtained from the gravitational-wave signal ($D_{\rm GW}= 43.8^{+2.9}_{-6.9}$ Mpc) and confirms that GW170817 occurred in NGC 4993.
Some of the heavy elements, such as gold and europium (Eu), are almost exclusively formed by the rapid neutron capture process (r-process). However, it is still unclear which astrophysical site between core-collapse supernovae and neutron star - neutron star (NS-NS) mergers produced most of the r-process elements in the universe. Galactic chemical evolution (GCE) models can test these scenarios by quantifying the frequency and yields required to reproduce the amount of Eu observed in galaxies. Although NS-NS mergers have become popular candidates, their required frequency (or rate) needs to be consistent with that obtained from gravitational wave measurements. Here we address the first NS-NS merger detected by LIGO/Virgo (GW170817) and its associated Gamma-ray burst and analyze their implication on the origin of r-process elements. Among other elements, we find that this event has produced between 15 and 70 Earth masses of gold. The range of NS-NS merger rate densities of 320$-$4740 Gpc$^{-3}$ yr$^{-1}$ provided by LIGO/Virgo is remarkably consistent with the range required by GCE to explain the Eu abundances in the Milky Way with NS-NS mergers, assuming a typical r-process abundance pattern for the ejecta. When using theoretical calculations to derive Eu yields, constraining the role of NS-NS mergers becomes more challenging because of nuclear astrophysics uncertainties. This is the first study that directly combines nuclear physics uncertainties with GCE calculations. If GW170817 is a representative event, NS-NS mergers can produce Eu in sufficient amount and are likely to be the main r-process site.
Dwarf galaxy anomalies, such as their abundance and cusp-core problems, remain a prime challenge in our understanding of galaxy formation. The inclusion of baryonic physics could potentially solve these issues, but the efficiency of stellar feedback is still controversial. We analytically explore the possibility of feedback from Active Galactic Nuclei (AGN) in dwarf galaxies and compare AGN and supernova (SN) feedback. We assume the presence of an intermediate mass black hole within low mass galaxies and standard scaling relations between the relevant physical quantities. We model the propagation and properties of the outflow and explore the critical condition for global gas ejection. Performing the same calculation for SNe, we compare the ability of AGN and SNe to drive gas out of galaxies. We find that a critical halo mass exists below which AGN feedback can remove gas from the host halo and that the critical halo mass for AGN is greater than the equivalent for SNe in a significant part of the parameter space, suggesting that AGN could provide an alternative and more successful source of negative feedback than SNe, even in the most massive dwarf galaxies.
We investigate how Einstein rings and magnified arcs are affected by small-mass dark-matter haloes placed along the line-of-sight to gravitational lens systems. By comparing the gravitational signature of line-of-sight haloes with that of substructures within the lensing galaxy, we derive a mass-redshift relation that allows us to rescale the detection threshold (i.e. lowest detectable mass) for substructures to a detection threshold for line-of-sight haloes at any redshift. We then quantify the line-of-sight contribution to the total number density of low-mass objects that can be detected through strong gravitational lensing. Finally, we assess the degeneracy between substructures and line-of-sight haloes of different mass and redshift to provide a statistical interpretation of current and future detections, with the aim of distinguishing between CDM and WDM. We find that line-of-sight haloes statistically dominate with respect to substructures, by an amount that strongly depends on the source and lens redshifts, and on the chosen dark matter model. Substructures represent about 30 percent of the total number of perturbers for low lens and source redshifts (as for the SLACS lenses), but less than 10 per cent for high redshift systems. We also find that for data with high enough signal-to-noise ratio and angular resolution, the non-linear effects arising from a double-lens-plane configuration are such that one is able to observationally recover the line-of-sight halo redshift with an absolute error precision of 0.15 at the 68 per cent confidence level.
A longstanding problem is identifying the elusive progenitors of Type Ia
supernovae (SNe Ia), which can roughly be split into Chandraksekhar and
sub-Chandrasekhar mass events. An important difference between these two cases
is the nucleosynthetic yield, which is altered by the increased neutron excess
in Chandrasekhar progenitors due to their pre-explosion simmering and high
central density.
From comparison with theoretical nucleosynthesis yields, we show that the
chemical composition of the most metal-rich star in the Ursa Minor dwarf
galaxy, COS 171, is dominated by nucleosynthesis from a low-metallicity,
low-mass, sub-Chandrasekhar mass SN Ia. Key diagnostic abundance ratios include
C/Fe, {\alpha}/Fe, Mn/Fe and Ni/Fe ratios, which could not have been produced
by Chandrasekhar-mass SNe Ia, Core-Collapse Type II supernovae or
Pair-Instsability supernovae. Strong deficiencies of Ni/Fe, Cu/Fe and Zn/Fe
also suggest the absence of alpha-rich freeze-out nucleosynthesis, favoring
low-mass WD progenitor SNe Ia. Based on comparisons of the measured Mn/Fe and
Si/Fe ratios with detonation models, we estimate a WD mass near 0.95Msun. We
also compare Mn/Fe and Ni/Fe ratios to the recent theoretical yields predicted
by Shen et al., finding consistent results. To explain the COS 171 [Fe/H], at
-1.35 dex, requires dilution of the ejecta from a single SNIa event with ~10^4
Msun of material; this is expected for a SN Ia remnant expanding into a warm
interstellar medium with n~1 /cm^3. In the future, finding more stars with the
unique chemical signatures we highlight here will be important for constraining
the rate and environments of sub-Chandrasekhar SNe Ia.
The ground-breaking detections of gravitational waves from black hole mergers by LIGO have rekindled interest in primordial black holes (PBHs) and the possibility of dark matter being composed of PBHs. It has been suggested that PBHs of tens of solar masses could serve as dark matter candidates. Recent analytical studies demonstrated that compact ultra-faint dwarf galaxies can serve as a sensitive test for the PBH dark matter hypothesis, since stars in such a halo-dominated system would be heated by the more massive PBHs, their present-day distribution can provide strong constraints on PBH mass. In this study, we further explore this scenario with more detailed calculations, using a combination of dynamical simulations and Bayesian inference methods. The joint evolution of stars and PBH dark matter is followed with a Fokker-Planck code PhaseFlow. We run a large suite of such simulations for different dark matter parameters, then use a Markov Chain Monte Carlo approach to constrain the PBH properties with observations of ultra-faint galaxies. We find that two-body relaxation between the stars and PBH drives up the stellar core size, and increases the central stellar velocity dispersion. Using the observed half-light radius and velocity dispersion of stars in the compact ultra-faint dwarf galaxies as joint constraints, we infer that these dwarfs may have a cored dark matter halo with the central density in the range of 1-2 $\rm{M_{\odot}/pc^3}$, and that the PBHs may have a mass range of 2-14 $\rm{M_{\odot}}$ if they constitute all or a substantial fraction of the dark matter.
Radio astrometric campaigns using VLBI have provided distances and proper motions for masers associated with young massive stars (BeSSeL survey). The ongoing BAaDE project plans to obtain astrometric information of SiO maser stars located in the inner Galaxy. These stars are associated with evolved, mass-losing stars. By overlapping optical (Gaia), infrared (2MASS, MSX and WISE) and radio (BAaDE) sources, we expect to obtain important clues on the intrinsic properties and population distribution of late-type stars. Moreover, a comparison of the Galactic parameters obtained with Gaia and VLBI can be done using radio observations on different targets: young massive stars (BeSSeL) and evolved stars (BAaDE).
On 2017 August 17, gravitational waves were detected from a binary neutron star merger, GW170817, along with a coincident short gamma-ray burst, GRB170817A. An optical transient source, Swope Supernova Survey 17a (SSS17a), was subsequently identified as the counterpart of this event. We present ultraviolet, optical and infrared light curves of SSS17a extending from 10.9 hours to 18 days post-merger. We constrain the radioactively-powered transient resulting from the ejection of neutron-rich material. The fast rise of the light curves, subsequent decay, and rapid color evolution are consistent with multiple ejecta components of differing lanthanide abundance. The late-time light curve indicates that SSS17a produced at least ~0.05 solar masses of heavy elements, demonstrating that neutron star mergers play a role in r-process nucleosynthesis in the Universe.
We present a near-infrared spectral sequence of the electromagnetic counterpart to the binary neutron star merger GW170817 detected by Advanced LIGO/Virgo. Our dataset comprises seven epochs of J+H spectra taken with FLAMINGOS-2 on Gemini-South between 1.5 and 10.5 days after the merger. In the initial epoch, the spectrum is dominated by a smooth blue continuum due to a high-velocity, lanthanide-poor blue kilonova component. Starting the following night, all subsequent spectra instead show features that are similar to those predicted in model spectra of material with a high concentration of lanthanides, including spectral peaks near 1.07 and 1.55 microns. Our fiducial model with 0.04 M_sun of ejecta, an ejection velocity of v=0.1c, and a lanthanide concentration of X_lan=1e-2 provides a good match to the spectra taken in the first five days, although it over-predicts the late-time fluxes. We also explore models with multiple fitting components, in each case finding that a significant abundance of lanthanide elements is necessary to match the broad spectral peaks that we observe starting at 2.5 d after the merger. These data provide direct evidence that binary neutron star mergers are significant production sites of even the heaviest r-process elements.
The circumstellar envelope of the hypergiant star IRC+10420 has been traced as far out in SiO J=2-1 as in CO J = 1-0 and CO J = 2-1, in dramatic contrast with the centrally condensed (thermal) SiO- but extended CO-emitting envelopes of giant and supergiant stars. Here, we present an observation of the circumstellar envelope in SiO J=1-0 that, when combined with the previous observation in {\sioii}, provide more stringent constraints on the density of the SiO-emitting gas than hitherto possible. The emission in SiO peaks at a radius of $\sim$2\arcsec\ whereas that in SiO J=2-1 emission peaks at a smaller radius of $\sim$1\arcsec, giving rise to their ring-like appearances. The ratio in brightness temperature between SiO J=1-0 and SiO J=2-1 decreases from a value well above unity at the innermost measurable radius to about unity at radius of $\sim$2\arcsec, beyond which this ratio remains approximately constant. Dividing the envelope into three zones as in models for the CO J = 1-0 and CO J = 2-1 emission, we show that the density of the SiO-emitting gas is comparable with that of the CO-emitting gas in the inner zone, but at least an order of magnitude higher by comparison in both the middle and outer zones. The SiO-emitting gas therefore originates from dense clumps, likely associated with the dust clumps seen in scattered optical light, surrounded by more diffuse CO-emitting interclump gas. We suggest that SiO molecules are released from dust grains due to shock interactions between the dense SiO-emitting clumps and the diffuse CO-emitting interclump gas.
We have analysed the [OI]6300 A line in a sample of 131 young stars with discs in the Lupus, Chamaeleon and signa Orionis star forming regions, observed with the X-shooter spectrograph at VLT. The stars have mass accretion rates spanning from 10^{-12} to 10^{-7} Mo/yr. The line profile was deconvolved into a low velocity component (LVC, < 40 km/s) and a high velocity component (HVC, > 40 km/s ), originating from slow winds and high velocity jets, respectively. The LVC is by far the most frequent component, with a detection rate of 77%, while only 30% of sources have a HVC. The [OI]6300 luminosity of both the LVC and HVC, when detected, correlates with stellar and accretion parameters of the central sources (i.e. Lstar , Mstar , Lacc , Macc), with similar slopes for the two components. The line luminosity correlates better with the accretion luminosity than with the stellar luminosity or stellar mass. We suggest that accretion is the main drivers for the line excitation and that MHD disc-winds are at the origin of both components. In the sub-sample of Lupus sources observed with ALMA a relationship is found between the HVC peak velocity and the outer disc inclination angle, as expected if the HVC traces jets ejected perpendicularly to the disc plane. Mass loss rates measured from the HVC span from ~ 10^{-13} to ~10^{-7} Mo/yr. The corresponding Mloss/Macc ratio ranges from ~0.01 to ~0.5, with an average value of 0.07. However, considering the upper limits on the HVC, we infer a ratio < 0.03 in more than 40% of sources. We argue that most of these sources might lack the physical conditions needed for an efficient magneto-centrifugal acceleration in the star-disc interaction region. Systematic observations of populations of younger stars, that is, class 0/I, are needed to explore how the frequency and role of jets evolve during the pre-main sequence phase.
Gaia Data Release 1 allows to recalibrate standard candles such as the Red Clump stars. To use those, they first need to be accurately characterised. In particular, colours are needed to derive the interstellar extinction. As no filter is available for the first Gaia data release and to avoid the atmosphere model mismatch, an empirical calibration is unavoidable. The purpose of this work is to provide the first complete and robust photometric empirical calibration of the Gaia Red Clump stars of the solar neighbourhood, through colour-colour, effective temperature-colour and absolute magnitude-colour relations, from the Gaia, Johnson, 2MASS, Hipparcos, Tycho-2, APASS-SLOAN and WISE photometric systems, and the APOGEE DR13 spectroscopic temperatures. We used a 3D extinction map to select low reddening red giants. To calibrate the colour-colour and the effective temperature-colour relations, we developed a MCMC method which accounts for all variable uncertainties and selects the best model for each photometric relation. We estimate the Red Clump absolute magnitude through the mode of a kernel-based distribution function. We provide 20 colour vs G-Ks relations and the first Teff vs G-Ks calibration. We obtained the Red Clump absolute magnitudes for 15 photometric bands with, in particular, M_Ks = -1.606 +/- 0.009 and M_G = (0.495 +/- 0.009) + (1.121 +/- 0.128) (G-Ks-2.1). We present an unreddened Gaia-TGAS HR diagram and use the calibrations to compare its Red Clump and its Red Giant Branch Bump with the Padova isochrones.
Recently, the optical counterpart of a gravitational wave source GW170817 has been identified in NGC 4993 galaxy. Together with evidence from observations in electromagnetic waves, the event has been suggested as a result of a merger of two neutron stars. We analyze the multi-wavelength data to characterize the host galaxy property and its distance to examine if the properties of NGC 4993 are consistent with this picture. Our analysis shows that NGC 4993 is a bulge-dominated galaxy with reff ~ 2-3 kpc and the Sersic index of n = 3-4 for the bulge component. The spectral energy distribution from 0.15 to 24 micron indicates that this galaxy has no significant ongoing star formation, the mean stellar mass of (0.3 - 1.2) times 10^11 Msun,the mean stellar age greater than ~3 Gyr, and the metallicity of about 20% to 100% of solar abundance. Optical images reveal dust lanes and extended features that suggest a past merging activity. Overall, NGC 4993 has characteristics of normal, but slightly disturbed elliptical galaxies. Furthermore, we derive the distance to NGC 4993 with the fundamental plane relation using 17 parameter sets of 7 different filters and the central stellar velocity dispersion from literature, finding an angular diameter distance of 37.7 +- 8.7 Mpc. NGC 4993 is similar to some host galaxies of short gamma-ray bursts but much different from those of long gamma-ray bursts, supporting the picture of GW170817 as a result of a merger of two NSs.
Using $\sim 100$ X-ray selected clusters in the Dark Energy Survey Science Verification data, we constrain the luminosity function of cluster red sequence galaxies in the redshift range of $0.1<z<1.05$. We develop a hierarchical Bayesian method that simultaneously models redshift evolution and cluster mass dependence. The results from this method are tested by red sequence luminosity function parameters derived in cluster redshift or mass bins. We find a hint that the faint end slope of a Schechter function fit may evolve with redshift at a significance level of $\sim 1.9 \sigma$. Faint cluster red sequence galaxies possibly become more abundant at lower redshift, indicating a different formation time from the bright red sequence galaxies. Optical cluster cosmology analyses may wish to consider this effect when deriving mass proxies. We also constrain the amplitude of the luminosity function with the hierarchical Bayesian method, which strongly correlates with cluster mass and provides an improved estimate of cluster masses. This technique can be applied to a larger sample of X-ray or optically selected clusters from the Dark Energy Survey, significantly improving the sensitivity of the analysis.
Links to: arXiv, form interface, find, astro-ph, recent, 1710, contact, help (Access key information)
We examine the fraction of massive ($M_{*}>10^{10} M_{\odot}$), compact star-forming galaxies (cSFGs) that host an active galactic nucleus (AGN) at $z\sim2$. These cSFGs are likely the direct progenitors of the compact quiescent galaxies observed at this epoch, which are the first population of passive galaxies to appear in large numbers in the early Universe. We identify cSFGs that host an AGN using a combination of Hubble WFC3 imaging and Chandra X-ray observations in four fields: the Chandra Deep Fields, the Extended Groth Strip, and the UKIDSS Ultra Deep Survey field. We find that $39.2^{+3.9}_{-3.6}$\% (65/166) of cSFGs at $1.4<z<3.0$ host an X-ray detected AGN. This fraction is 3.2 times higher than the incidence of AGN in extended star-forming galaxies with similar masses at these redshifts. This difference is significant at the $6.2\sigma$ level. Our results are consistent with models in which cSFGs are formed through a dissipative contraction that triggers a compact starburst and concurrent growth of the central black hole. We also discuss our findings in the context of cosmological galaxy evolution simulations that require feedback energy to rapidly quench cSFGs. We show that the AGN fraction peaks precisely where energy injection is needed to reproduce the decline in the number density of cSFGs with redshift. Our results suggest that the first abundant population of massive, quenched galaxies emerged directly following a phase of elevated supermassive black hole growth and further hints at a possible connection between AGN and the rapid quenching of star formation in these galaxies.
Spectroscopic studies of ices in nearby star-forming regions indicate that ice mantles form on dust grains in two distinct steps, starting with polar ice formation (H2O rich) and switching to apolar ice (CO rich). We test how well the picture applies to more diffuse and quiescent clouds where the formation of the first layers of ice mantles can be witnessed. Medium-resolution near-infrared spectra are obtained toward background field stars behind the Pipe Nebula. The water ice absorption is positively detected at 3.0 micron in seven lines of sight out of 21 sources for which observed spectra are successfully reduced. The peak optical depth of the water ice is significantly lower than those in Taurus with the same visual extinction. The source with the highest water-ice optical depth shows CO ice absorption at 4.7 micron as well. The fractional abundance of CO ice with respect to water ice is 16+7-6 %, and about half as much as the values typically seen in low-mass star-forming regions. A small fractional abundance of CO ice is consistent with some of the existing simulations. Observations of CO2 ice in the early diffuse phase of a cloud play a decisive role in understanding the switching mechanism between polar and apolar ice formation.
We present a comparative analysis of the properties of passive, star-forming and transition (green valley) galaxies, in four discrete environments: field, groups, the outskirts and the core of X-ray clusters. We construct samples of galaxies from the SDSS in these environments so that they are bound to have similar redshift distributions. The classification of galaxies into the three sequences is based on the UV-optical colour $NUV-r$. We study a number of galaxy properties: stellar mass, morphology, specific star formation rate and the history of star formation. The analysis of green valley galaxies reveals that the physical mechanisms responsible for external quenching become more efficient moving from the field to denser environments. We confirm previous findings that green valley galaxies have intermediate morphologies, moreover, we find that this appears to be independent of the environment. Regarding the stellar mass of green valley galaxies, we find that they tend to be more massive in the field than in denser environments. On average, green valley galaxies account for $\sim 20\%$ of all galaxies in groups and X-ray clusters. We find evidence that the field environment is inefficient in transforming low mass galaxies. Green valley galaxies have average star formation histories intermediate between passive and star forming galaxies, and have a clear and consistent dependence on the environment: both, the quenching time, and the amplitude of the star formation rate, decrease towards higher density environments.
Using high-resolution magnetohydrodynamic simulations of idealized, non-cosmological galaxies, we investigate how cooling, star formation, and stellar feedback affect galactic magnetic fields. We find that the amplification histories, saturation values, and morphologies of the magnetic fields vary considerably depending on the baryonic physics employed, primarily because of differences in the gas density distribution. In particular, adiabatic runs and runs with a sub-grid (effective equation of state) stellar feedback model yield lower saturation values and morphologies that exhibit greater large-scale order compared with runs that adopt explicit stellar feedback and runs with cooling and star formation but no feedback. The discrepancies mostly lie in gas denser than the galactic average, which requires cooling and explicit fragmentation to capture. Independent of the baryonic physics included, the magnetic field strength scales with gas density as $B\propto n^{2/3}$, suggesting isotropic flux freezing or equipartition between the magnetic and gravitational energies during the field amplification. We conclude that accurate treatments of cooling, star formation, and stellar feedback are crucial for obtaining the correct magnetic field strength and morphology in dense gas, which, in turn, is essential for properly modeling other physical processes that depend on the magnetic field, such as cosmic ray feedback.
This paper provides a brief overview of the journey of molecules through the Cosmos, from local diffuse interstellar clouds and PDRs to distant galaxies, and from cold dark clouds to hot star-forming cores, protoplanetary disks, planetesimals and exoplanets. Recent developments in each area are sketched and the importance of connecting astronomy with chemistry and other disciplines is emphasized. Fourteen challenges for the field of Astrochemistry in the coming decades are formulated.
Lyman-$\alpha$ emitters (LAEs) are a promising probe of the large-scale structure at high redshift, $z\gtrsim 2$. In particular, the Hobby-Eberly Telescope Dark Energy Experiment aims at observing LAEs at 1.9 $<z<$ 3.5 to measure the Baryon Acoustic Oscillation (BAO) scale and the Redshift-Space Distortion (RSD). However, Zheng et al. (2011) pointed out that the complicated radiative transfer (RT) of the resonant Lyman-$\alpha$ emission line generates an anisotropic selection bias in the LAE clustering on large scales, $s\gtrsim 10$ Mpc. This effect could potentially induce a systematic error in the BAO and RSD measurements. Also, Croft et al. (2016) claims an observational evidence of the effect in the Lyman-$\alpha$ intensity map, albeit statistically insignificant. We aim at quantifying the impact of the Lyman-$\alpha$ RT on the large-scale galaxy clustering in detail. For this purpose, we study the correlations between the large-scale environment and the ratio of an apparent Lyman-$\alpha$ luminosity to an intrinsic one, which we call the `observed fraction', at $2<z<6$. We apply our Lyman-$\alpha$ RT code by post-processing the full Illustris simulations. We simply assume that the intrinsic luminosity of the Lyman-$\alpha$ emission is proportional to the star formation rate of galaxies in Illustris, yielding a sufficiently large sample of LAEs to measure the anisotropic selection bias. We find little correlations between large-scale environment and the observed fraction induced by the RT, and hence a smaller anisotropic selection bias than what was claimed by Zheng et al. (2011). We argue that the anisotropy was overestimated in the previous work due to the insufficient spatial resolution: it is important to keep the resolution such that it resolves the high density region down to the scale of the interstellar medium, $\sim1$ physical kpc. (abridged)
We report the discovery of one of the largest and most distant Giant Radio Galaxy (GRG) in the Lynx field which was discovered using deep Giant Metre-wave Radio Telescope (GMRT) 150 MHz observations. The core is detected at 150 MHz and also in the VLA FIRST survey. Spectroscopic observations carried out using the IUCAA Giravali Observatory(IGO) provided a redshift value of 0.57. This redshift was later confirmed with data from the Sloan Digital Sky Survey (Data Release 12). The angular size of the GRG is 5.5 arcmin and at the redshift of 0.57, its linear size is 2.2 Mpc. At this high redshift, only a few radio sources are known to have such large linear size. In order to estimate the spectral index of the bridge emission as well as the spectral age of the source, we observed this source at L-band, 610 MHz and 325 MHz bands with the GMRT. We present the spectral ageing analysis of the source which puts an upper limit of 20 Myr on the spectral age. The better resolution maps presented here as opposed to the original 150 MHz map shows evidence for a second episode of emission. We also find that the core is detected at all four frequencies with a spectral index of 0.85, which is steeper than normal, hence we speculate that the core may be a compact steep spectrum source (CSS), which makes this giant radio galaxy a candidate triple-double radio galaxy.
We present observations in CO(3-2) that, combined with previous observations in CO(2-1), constrain the physical properties of the filamentary molecular gas in the central $\sim$6.5 kpc of NGC 1275, the central giant elliptical galaxy of the Perseus cluster. We find this molecular gas to have a temperature $\gtrsim 20$ K and a density $\sim$$10^2$-$10^4 {\rm \ cm^{-3}}$, typically warmer and denser than the bulk of Giant Molecular Clouds (GMCs) in the Galaxy. Bathed in the harsh radiation and particle field of the surrounding intracluster X-ray gas, the molecular gas likely has a much higher ionization fraction than that of GMCs. For an ionization fraction of $\sim$$10^{-4}$, similar to that of Galactic diffuse ($\lesssim 250 {\rm \ cm^{-3}}$) partially-molecular clouds that emit in HCN(1-0) and HCO$^+$(1-0), we show that the same gas traced in CO can produce the previously reported emissions in HCN(3-2), HCO$^+$(3-2), and CN(2-1) from NGC 1275; the dominant source of excitation for all the latter molecules is collisions with electrons. To prevent collapse, as evidenced by the lack of star formation in the molecular filaments, they must consist of thin strands that have cross-sectional radii $\lesssim$0.2-2 pc if supported solely by thermal gas pressure; larger radii are permissible if turbulence or poloidal magnetic fields provide additional pressure support. We point out that the conditions required to relate CO luminosities to molecular gas masses in our Galaxy are unlikely to apply in cluster central elliptical galaxies. Rather than being virialized structures analogous to GMCs, we propose that the molecular gas in NGC 1275 comprises pressure-confined structures created by turbulent flows.
The direct dynamical paths leading to Associative Detachment (AD) in the gas-phase, and specifically in the low-temperature regions of the Dark Molecular Clouds (DMC) in the ISM, or in cold trap laboratory experiments, are investigated with quantum chemical methods by using a high-level multi-reference Configuration Interaction (CI) approach that employs single and double excitations plus Davidson perturbative correction [MRSDCI(Q)] and the d-aug-cc-pV5Z basis set. The potential energy curves for H + CN- are constructed for different directions of the H partner approaching the CN- anion within the framework of the Born-Oppenheimer approximation. The present calculations found that the AD energetics at low temperature becomes favorable only along a selected range of approaching directions, thus showing that there is a preferred path of forming HCN at low temperatures, while that of forming its HNC isomer is found to be energetically forbidden. Given the existence in the ISM of different HCN/HNC ratios in different environments, we discuss the implications of our findings for selective formation of either isomer in the low-temperature conditions of the Molecular Cloud Cores.
We review star formation in molecular clouds and describe why magnetic fields may be important and how they can influence filamentary structure and the column density probability distribution function (PDF). We also comment on the origin of the stellar and substellar initial mass function (IMF), which may require explanations beyond a simple Jeans length argument in turbulent molecular clouds. A mathematical model of the modified lognormal power-law (MLP) distribution function provides a framework within which to connect accretion processes with the IMF.
The relationship between the integrated H$\beta$ line luminosity and the velocity dispersion of the ionized gas of HII galaxies and giant HII regions represents an exciting standard candle that presently can be used up to redshifts z ~ 4. Locally it is used to obtain precise measurements of the Hubble constant by combining the slope of the relation obtained from nearby ($z \leq $ 0.2) HII galaxies with the zero point determined from giant HII regions belonging to an `anchor sample' of galaxies for which accurate redshift-independent distance moduli are available. We present new data for 36 giant HII regions in 13 galaxies of the anchor sample that includes the megamaser galaxy NGC 4258. Our data is the result of the first four years of observation of our primary sample of 130 giant HII regions in 73 galaxies with Cepheid determined distances. Our best estimate of the Hubble parameter is $71.0\pm2.8(random)\pm2.1(systematic)$ km /s Mpc This result is the product of an independent approach and, although at present less precise than the latest SNIa results, it is amenable to substantial improvement.
We study the nature of phase transitions between dilute and dense axion stars interpreted as self-gravitating Bose-Einstein condensates. We develop a Newtonian model based on the Gross-Pitaevskii-Poisson equations for a complex scalar field with a self-interaction potential $V(|\psi|^2)$ involving an attractive $|\psi|^4$ term and a repulsive $|\psi|^6$ term. Using a Gaussian ansatz for the wave function, we analytically obtain the mass-radius relation of dilute and dense axion stars for arbitrary values of the self-interaction parameter $\lambda\le 0$. We show the existence of a critical point $|\lambda|_c\sim (m/M_P)^2$ above which a first order phase transition takes place. We qualitatively estimate general relativistic corrections on the mass-radius relation of axion stars. For weak self-interactions $|\lambda|<|\lambda|_c$, a system of self-gravitating axions forms a stable dilute axion star below a general relativistic maximum mass $M_{\rm max,GR}^{\rm dilute}\sim M_P^2/m$ and collapses into a black hole above that mass. For strong self-interactions $|\lambda|>|\lambda|_c$, a system of self-gravitating axions forms a stable dilute axion star below a Newtonian maximum mass $M_{\rm max,N}^{\rm dilute}=5.073 M_P/\sqrt{|\lambda|}$, collapses into a dense axion star above that mass, and collapses into a black hole above a general relativistic maximum mass $M_{\rm max,GR}^{\rm dense}\sim \sqrt{|\lambda|}M_P^3/m^2$. Dense axion stars explode below a Newtonian minimum mass $M_{\rm min,N}^{\rm dense}\sim m/\sqrt{|\lambda|}$ and form dilute axion stars of large size or disperse away. We determine the phase diagram of self-gravitating axions and show the existence of a triple point $(|\lambda|_*,M_*/(M_P^2/m))$ separating dilute axion stars, dense axion stars, and black holes. We make numerical applications for QCD axions and ultralight axions.
Links to: arXiv, form interface, find, astro-ph, recent, 1710, contact, help (Access key information)
The C III] 1907,1909 emission doublet has been proposed as an alternative to Lyman-alpha in redshift confirmations of galaxies at z > 6 since it is not attenuated by the largely neutral intergalactic medium at these redshifts and is believed to be strong in the young, vigorously star-forming galaxies present at these early cosmic times. We present a statistical sample of 17 C III]-emitting galaxies beyond z~1.5 using 30 hour deep VLT/MUSE integral field spectroscopy covering 2 square arcminutes in the Hubble Deep Field South (HDFS) and Ultra Deep Field (UDF), achieving C III] sensitivities of ~2e-17 erg/s/cm^2 in the HDFS and ~7e-18 erg/s/cm^2 in the UDF. The rest-frame equivalent widths range from 2 to 19 Angstroms. These 17 galaxies represent ~3% of the total sample of galaxies found between 1.5 < z < 4. They also show elevated star formation rates, lower dust attenuation, and younger mass-weighted ages than the general population of galaxies at the same redshifts. Combined with deep slitless grism spectroscopy from the HST/WFC3 in the UDF, we can tie the rest-frame ultraviolet C III] emission to rest-frame optical emission lines, namely [O III] 5007, finding a strong correlation between the two. Down to the flux limits that we observe (~1e-18 erg/s/cm^2 with the grism data in the UDF), all objects with a rest-frame [O III] 4959,5007 equivalent width in excess of 250 Angstroms, the so-called Extreme Emission Line Galaxies, have detections of C III] in our MUSE data. More detailed studies of the C III]-emitting population at these intermediate redshifts will be crucial to understand the physical conditions in galaxies at early cosmic times and to determine the utility of C III] as a redshift tracer.
Here I report the scaling relation between the baryonic masses and the scale lengths of stellar discs from $\sim$1000 morphologically late-type galaxies. The baryonic mass-size relation is a single power-law $R_\ast \propto M_b^{0.38}$ across $\sim$3 orders of magnitude in baryonic mass. The scatter in size at fixed baryonic mass is nearly constant and there is essentially no outlier. The baryonic mass-size relation provides a more fundamental description of the structure of the discs than the stellar mass-size relation. The slope and the scatter of the stellar mass-size relation can be understood in the context of the baryonic mass-size relation. For gas-rich galaxies, the stars is no longer a good tracer for the baryons. High baryonic mass, gas-rich galaxies appear to be much larger at fixed stellar mass because most of the baryonic content is gas. The stellar mass-size relation thus deviates from the power law baryonic relation and the scatter increases at the low stellar mass end. Those extremely gas-rich low-mass galaxies can be classified as Ultra Diffuse Galaxies based on the structure.
Supernova driven winds are often invoked to remove chemically enriched gas from dwarf galaxies to match their low observed metallicities. In such shallow potential wells, outflows may produce massive amounts of enriched halo gas (circum-galactic medium or CGM) and pollute the intergalactic medium (IGM). Here, we present a survey of the CGM and IGM around $18$ star-forming field dwarfs with stellar masses of $\log\,M_*/M_\odot\approx8-9$ at $z\approx0.2$. Eight of these have CGM probed by quasar absorption spectra at projected distances, $d$, less than the host virial radius, $R_{\rm h}$. Ten are probed in the surrounding IGM at $d/R_{\rm h}=1-3$. The absorption measurements include neutral hydrogen, the dominant silicon ions for diffuse cool gas ($T\sim10^4$ K; Si II, Si III, and Si IV), moderately ionized carbon (C IV), and highly ionized oxygen (O VI). Metal absorption from the CGM of the dwarfs is less common and $\approx4\times$ weaker compared to massive star-forming galaxies though O VI absorption is still common. None of the dwarfs probed at $d/R_{\rm h}=1-3$ have definitive metal-line detections. Combining the available silicon ions, we estimate that the cool CGM of the dwarfs accounts for only $2-6\%$ of the expected silicon budget from the yields of supernovae associated with past star-formation. The highly ionized O VI accounts for $\approx8\%$ of the oxygen budget. As O VI traces an ion with expected equilibrium ion fractions of $\lesssim 0.2$, the highly ionized CGM may represent a significant metal reservoir even for dwarfs not expected to maintain gravitationally shock heated hot halos.
We use a cosmological hydrodynamic simulation calculated with Enzo and the semi-analytic galaxy formation model (SAM) GAMMA to address the chemical evolution of dwarf galaxies in the early universe. The overall goal is to better understand the origin of metal-poor stars and the formation of dwarf galaxies and the Milky Way halo by cross-validating these theoretical approaches. We combine GAMMA with the merger tree of the most massive galaxy found in the hydrodynamic simulation and compare the star formation rate, the metallicity distribution function (MDF), and the age-metallicity relationship predicted by the two approaches. We found that the SAM can reproduce the global trends of the hydrodynamic simulation. However, there are degeneracies and more constraints need to be extracted from the simulation to isolate the correct semi-analytic solution. Non-uniform mixing in the galaxy's interstellar medium, coming primarily from self-enrichment by local supernovae, causes a broadening in the MDF that can be emulated in the SAM by convolving its predicted MDF with a Gaussian function having a standard deviation of $\sim$0.2 dex. However, stochastic processes such as bursty star formation histories and star formation triggered by supernova explosions cannot be reproduced by the current version of GAMMA. Moreover, we found that massive stars in building-block galaxies tend to explode and generate outflows while falling inside the main galaxy's halo, leading to a complex multiphase circumgalactic medium with a wide range of temperatures, densities, and metallicities, as opposed to the hot, isothermal, well-mixed halo gas in equilibrium assumed in SAMs.
A large population of extended substructures generates a stochastic gravitational field that is fully specified by the function $p({\bf F})$, which defines the probability that a tracer particle experiences a force $\bf F$ within the interval ${\bf F},{\bf F}+ d\bf F$. This paper presents a statistical technique for deriving the spectrum of random fluctuations directly from the number density of substructures with known mass and size functions. Application to the subhalo population found in cold dark matter simulations of Milky Way-sized haloes shows that, while the combined force distribution is governed by the most massive satellites, the fluctuations of the {\it tidal} field are completely dominated by the smallest and most abundant subhaloes. In light of this result we discuss observational experiments that may be sufficiently sensitive to Galactic tidal fluctuations to probe the "dark" low-end of the subhalo mass function and constrain the particle mass of warm and ultra-light axion dark matter models.
We present deep, wide-field Subaru Hyper Suprime-Cam photometry of two recently discovered satellites of the Milky Way (MW): Columba I and Triangulum II. The color magnitude diagrams of both objects point to exclusively old and metal-poor stellar populations. We re-derive structural parameters and luminosities of these satellites, and find $M_{\rm V, Col~I} = -4.2\pm0.2$ for Col I and $M_{\rm V, Tri~II} = -1.2\pm0.4$ for Tri II, with corresponding half-light radii of $r_{\rm h, Col~I} = 117\pm17$ pc and $r_{\rm h, Tri~II} = 21\pm4$ pc. The properties of both systems are consistent with observed scaling relations for MW dwarf galaxies. Based on archival data, we derive upper limits on the neutral gas content of these dwarfs, and find that they lack HI, as do the majority of observed satellites within the MW virial radius. Neither satellite shows evidence of tidal stripping in the form of extensions or distortions in matched-filter stellar density maps or surface density profiles. However, the smaller Tri II system is relatively metal-rich for its luminosity (compared to other MW satellites), possibly because it has been tidally stripped. Through a suite of orbit simulations, we show that Tri II is approaching pericenter of its eccentric orbit, a stage at which tidal debris is unlikely to be seen. In addition, we find that Tri II may be on its first infall into the MW, which helps explain its unique properties among MW dwarfs. Further evidence that Tri II is likely an ultra-faint dwarf comes from its stellar mass function, which is similar to those of other MW dwarfs.
Ice mantles detected on the surface of dust grains toward the coldest regions of the interstellar medium can be photoprocessed by the secondary ultraviolet (UV) field present in dense cloud interiors. In this work, we present UV-irradiation experiments under astrophysically relevant conditions of pure NH3 ice samples in an ultra-high vacuum chamber where solid samples were deposited onto a substrate at 8 K. The ice analogs were subsequently photoprocessed with a microwave-discharged hydrogen- flow lamp. The induced radiation and photochemistry led to the production of H2, N2 and N2H4. In addition, photodesorption to the gas phase of the original ice component, NH3, and two of the three detected photoproducts, H2 and N2, was observed thanks to a quadrupole mass spectrometer (QMS). Calibration of the QMS allowed quantifica- tion of the photodesorption yields, leading to Ypd(NH3)=2.1+2.1-1.0 x 10^-3 molecules per incident photon, which remained constant during the whole experiments, while photodesorption of H2 and N2 increased with fluence, pointing toward an indirect photodesorption mecha- nism involving energy transfer for these species. Photodesorption yield of N2 molecules after a fluence equivalent to that experienced by ice mantles in space was similar to that of the NH3 molecules (Ypd (N2) = 1.7+1.7-0.9 x 10^-3 molecules incident photon ).
Little is known about the portion of the Milky Way lying beyond the Galactic center at distances of more than 9 kilo-parsec from the Sun. These regions are opaque at optical wavelengths due to absorption by interstellar dust, and distances are very large and hard to measure. We report a direct trigonometric parallax distance of 20.4_{-2.2}^{+2.8} kilo-parsec obtained with the Very Long Baseline Array to a water maser source in a region of active star formation. These measurements allow us to shed light on Galactic spiral structure by locating the Scutum-Centaurus spiral arm as it passes through the far side of the Milky Way, and to validate a kinematic method for determining distances in this region based on transverse motions.
We present high-frequency Karl G. Jansky Very Large Array (VLA) continuum and spectral line (NH3, H64${\alpha}$, and H63${\alpha}$) observations of the Galactic Center Radio Arc region, covering the Sickle H II region, the Quintuplet cluster, and molecular clouds M0.20-0.033 and M0.10-0.08. These observations show that the two velocity components of M0.20-0.033 (~25 & 80 km/s), previously thought to be separate clouds along the same line-of-sight, are physically connected in position-velocity space via a third southern component around 50 km/s. Further position-velocity analysis of the surrounding region, using lower-resolution survey observations taken with the Mopra and ATCA telescopes, indicates that both molecular components in M0.20-0.033 are physically connected to the M0.10-0.08 molecular cloud, which is suggested to be located on stream 1 in the Kruijssen et al. (2015) orbital model. The morphology and kinematics of the molecular gas in M0.20-0.033 indicate that the two velocity components in M0.20-0.033 constitute an expanding shell. Our observations suggest that the M0.20-0.033 expanding shell has an expansion velocity of 40 km/s, with a systemic velocity of 53 km/s, comparable to velocities detected in M0.10-0.08. The origin of the expanding shell is located near the Quintuplet cluster, suggesting that the energy and momentum output from this massive stellar cluster may have contributed to the expansion.
We investigate stellar metallicity distribution functions (MDFs), including Fe and ${\alpha}$-element abundances, in dwarf galaxies from the Feedback in Realistic Environments (FIRE) project. We examine both isolated dwarf galaxies and those that are satellites of a Milky Way-mass galaxy. In particular, we study the effects of including a sub-grid turbulent model for the diffusion of metals in gas. Simulations that include diffusion have narrower MDFs and abundance ratio distributions, because diffusion drives individual gas and star particles toward the average metallicity. This effect provides significantly better agreement with observed abundance distributions of dwarf galaxies in the Local Group, including the small intrinsic scatter in [${\alpha}$/Fe] vs. [Fe/H] (less than 0.1 dex). This small intrinsic scatter arises in our simulations because the interstellar medium (ISM) in dwarf galaxies is well-mixed at nearly all cosmic times, such that stars that form at a given time have similar abundances to within 0.1 dex. Thus, most of the scatter in abundances at z = 0 arises from redshift evolution and not from instantaneous scatter in the ISM. We find similar MDF widths and intrinsic scatter for satellite and isolated dwarf galaxies, which suggests that environmental effects play a minor role compared with internal chemical evolution in our simulations. Overall, with the inclusion of metal diffusion, our simulations reproduce abundance distribution widths of observed low-mass galaxies, enabling detailed studies of chemical evolution in galaxy formation.
We study the behavior of eight diffuse interstellar bands (DIBs) in different interstellar environments, as characterized by the fraction of hydrogen in molecular form [$f$(H$_2$)], with comparisons to the corresponding behavior of various known atomic and molecular species. The equivalent widths of the five "normal" DIBs ($\lambda\lambda$5780.5, 5797.1, 6196.0, 6283.8, and 6613.6), normalized to $E(B-V)$, show a "Lambda-shaped" behavior: they increase at low $f$(H$_2$), peak at $f$(H$_2$) ~ 0.3, and then decrease. The similarly normalized column densities of Ca, Ca$^+$, Ti$^+$, and CH$^+$ also decline for $f$(H$_2$) > 0.3. In contrast, the normalized column densities of Na, K, CH, CN, and CO increase monotonically with $f$(H$_2$), and the trends exhibited by the three C$_2$ DIBs ($\lambda\lambda$4726.8, 4963.9, and 4984.8) lie between those two general behaviors. These trends with $f$(H$_2$) are accompanied by cosmic scatter, the dispersion at any given $f$(H$_2$) being significantly larger than the individual errors of measurement. The Lambda-shaped trends suggest the balance between creation and destruction of the DIB carriers differs dramatically between diffuse atomic and diffuse molecular clouds; additional processes besides ionization and shielding are needed to explain those observed trends. Except for several special cases, the highest $W$(5780)/$W$(5797) ratios, characterizing the so-called "sigma-zeta effect", occur only at $f$(H$_2$) < 0.2. We propose a sequence of DIBs based on trends in their pair-wise strength ratios with increasing $f$(H$_2$). In order of increasing environmental density, we find the $\lambda$6283.8 and $\lambda$5780.5 DIBs, the $\lambda$6196.0 DIB, the $\lambda$6613.6 DIB, the $\lambda$5797.1 DIB, and the C$_2$ DIBs.
We have found the atomic gas (HI) reservoirs of the blue ultra diffuse galaxy (UDG) candidates identified by R\'oman and Trujillo in images near Hickson Compact Groups (HCGs). We confirm that all of the objects are indeed UDGs with effective radii R_e > 1.5 kpc. Three of them are likely to be gravitationally bound to the HCG near which they project, one is plausibly gravitationally bound to the nearest HCG, and one is in the background. We measure HI masses and velocity widths for each object directly from the spectra, and use the widths together with the UDG effective radii to estimate dynamical masses and halo spin parameters. The HI and dynamical masses of the blue UDGs span nearly two orders of magnitude; like red UDGs, the blue UDGs may form a heterogeneous class. The location of the blue UDGs in the HI mass - stellar mass plane and their distribution of spin parameter estimators are consistent with those of the broader gas-rich galaxy population. The blue UDGs studied here do not seem to have the high spin parameters predicted by some red UDG progenitor models, but the low-mass object with the largest effective radius is the most gas-rich as predicted by another. The data presented here highlight the potential of single-dish radio observations for measuring the physical properties of blue diffuse objects detected in the optical.
Although now routinely incorporated into hydrodynamic simulations of galaxy evolution, the true importance of the feedback effect of the outflows driven by active galactic nuclei (AGN) remains uncertain from an observational perspective. This is due to a lack of accurate information on the densities, radial scales and level of dust extinction of the outflow regions. Here we use the unique capabilities of VLT/Xshooter to investigate the warm outflows in a representative sample of 9 local (0.06 $<$ z $<$ 0.15) ULIRGs with AGN nuclei and, for the first time, accurately quantify the key outflow properties. We find that the outflows are compact (0.05 $<$ R$_{[OIII]}$ $<$ 1.2 kpc), significantly reddened (median E(B-V)$\sim$0.5 magnitudes), and have relatively high electron densities (3.4 $<$ log$_{10}$ n$_e$ (cm$^{-3}$) $<$ 4.8). It is notable that the latter densities -- obtained using trans-auroral [SII] and [OII] emission-line ratios -- exceed those typically assumed for the warm, emission-line outflows in active galaxies, but are similar to those estimated for broad and narrow absorption line outflow systems detected in some type 1 AGN. Even if we make the most optimistic assumptions about the true (deprojected) outflow velocities, we find relatively modest mass outflow rates ($0.07 < \dot{M} < 11$ M$_{sol}$ yr$^{-1}$) and kinetic powers measured as a fraction of the AGN bolometric luminosities ($4\times10^{-4} < \dot{E}/L_{BOL} <1$%). Therefore, although warm, AGN-driven outflows have the potential to strongly affect the star formation histories in the inner bulge regions ($r \sim$ 1kpc) of nearby ULIRGs, we lack evidence that they have a significant impact on the evolution of these rapidly evolving systems on larger scales.
We use the energy-balance code MAGPHYS to determine stellar and dust masses, and dust corrected star-formation rates for over 200,000 GAMA galaxies, 170,000 G10-COSMOS galaxies and 200,000 3D-HST galaxies. Our values agree well with previously reported measurements and constitute a representative and homogeneous dataset spanning a broad range in stellar mass (10^8---10^12 Msol), dust mass (10^6---10^9 Msol), and star-formation rates (0.01---100 Msol per yr), and over a broad redshift range (0.0 < z < 5.0). We combine these data to measure the cosmic star-formation history (CSFH), the stellar-mass density (SMD), and the dust-mass density (DMD) over a 12 Gyr timeline. The data mostly agree with previous estimates, where they exist, and provide a quasi-homogeneous dataset using consistent mass and star-formation estimators with consistent underlying assumptions over the full time range. As a consequence our formal errors are significantly reduced when compared to the historic literature. Integrating our cosmic star-formation history we precisely reproduce the stellar-mass density with an ISM replenishment factor of 0.50 +/- 0.07, consistent with our choice of Chabrier IMF plus some modest amount of stripped stellar mass. Exploring the cosmic dust density evolution, we find a gradual increase in dust density with lookback time. We build a simple phenomenological model from the CSFH to account for the dust mass evolution, and infer two key conclusions: (1) For every unit of stellar mass which is formed 0.0065---0.004 units of dust mass is also formed; (2) Over the history of the Universe approximately 90 to 95 per cent of all dust formed has been destroyed and/or ejected.
Dwarf spheroidal galaxies are regarded as the basic building blocks in the formation of larger galaxies and are believed to be the most dark matter dominated systems known in the Universe. There are several models that attempt to explain their formation and evolution, but they have problems to model the formation of isolated dwarf spheroidal galaxies. Here, we will explain a possible formation scenario in which star clusters form inside the dark matter halo of a dwarf spheroidal galaxy. Those star clusters suffer from low star formation efficiency and dissolve while orbiting inside the dark matter halo. Thereby, they build the faint luminous components that we observe in dwarf spheroidal galaxies. In this paper we study this model by adding different star formation histories to the simulations to compare the results with our previous work and observational data to show that we can explain the formation of dwarf spheroidal galaxies.
Merging neutron stars offer an exquisite laboratory for simultaneously studying strong-field gravity and matter in extreme environments. We establish the physical association of an electromagnetic counterpart EM170817 to gravitational waves (GW170817) detected from merging neutron stars. By synthesizing a panchromatic dataset, we demonstrate that merging neutron stars are a long-sought production site forging heavy elements by r-process nucleosynthesis. The weak gamma-rays seen in EM170817 are dissimilar to classical short gamma-ray bursts with ultra-relativistic jets. Instead, we suggest that breakout of a wide-angle, mildly-relativistic cocoon engulfing the jet elegantly explains the low-luminosity gamma-rays, the high-luminosity ultraviolet-optical-infrared and the delayed radio/X-ray emission. We posit that all merging neutron stars may lead to a wide-angle cocoon breakout; sometimes accompanied by a successful jet and sometimes a choked jet.
Accurate measurements of the physical structure of protoplanetary discs are critical inputs for planet formation models. These constraints are traditionally established via complex modelling of continuum and line observations. Instead, we present an empirical framework to locate the CO isotopologue emitting surfaces from high spectral and spatial resolution ALMA observations. We apply this framework to the disc surrounding IM Lupi, where we report the first direct, i.e. model independent, measurements of the radial and vertical gradients of temperature and velocity in a protoplanetary disc. The measured disc structure is consistent with an irradiated self-similar disc structure, where the temperature increases and the velocity decreases towards the disc surface. We also directly map the vertical CO snow line, which is located at about one gas scale height at radii between 150 and 300 au, with a CO freeze-out temperature of $21\pm2$ K. In the outer disc ($> 300$ au), where the gas surface density transitions from a power law to an exponential taper, the velocity rotation field becomes significantly sub-Keplerian, in agreement with the expected steeper pressure gradient. The sub-Keplerian velocities should result in a very efficient inward migration of large dust grains, explaining the lack of millimetre continuum emission outside of 300 au. The sub-Keplerian motions may also be the signature of the base of an externally irradiated photo-evaporative wind. In the same outer region, the measured CO temperature above the snow line decreases to $\approx$ 15 K because of the reduced gas density, which can result in a lower CO freeze-out temperature, photo-desorption, or deviations from local thermodynamic equilibrium.
We construct the error distribution of galactic rotation curve ($\Theta$) measurements using 134 data points from the 162 measurements compiled in De Grijs et al. (arXiv:1709.02501), following the same procedures used in previous works by Ratra and collaborators. We determine the weighted mean of these measurements to be $\Theta_{Mean} = 226.73 \pm 0.70$ km/sec and the median estimate is calculated to be $\Theta_{Med} = 234.66\pm 0.58$ km/sec. We also checked if the error distribution (constructed using both the weighted mean and median as the estimate) shows a Gaussian distribution. We find using both the estimates that it has much wider tails than a Gaussian distribution. We also tried to fit the data to four distributions: Gaussian, Cauchy, double-exponential, and Students-t. The best fit is obtained using the Students-$t$ distribution for $n=2$ using the median value as the central estimate, corresponding to $p$-value of 0.19. All other distributions provide poorer fits to the data.
The extragalactic background radiation produced by distant galaxies emitting in the far infrared limits the sensitivity of telescopes operating in this range due to confusion. We have constructed a model of the infrared background based on numerical simulations of the large-scale structure of the Universe and the evolution of dark matter halos. The predictions of this model agree well with the existing data on source counts. We have constructed maps of a sky field with an area of 1 deg$^2$ directly from our simulated observations and measured the confusion limit. At wavelengths $100-300$ $\mu$m the confusion limit for a 10-m telescope has been shown to be at least an order of magnitude lower than that for a 3.5-m one. A spectral analysis of the simulated infrared background maps clearly reveals the large-scale structure of the Universe. The two-dimensional power spectrum of these maps has turned out to be close to that measured by space observatories in the infrared. However, the fluctuations in the number of intensity peaks observed in the simulated field show no clear correlation with superclusters of galaxies; the large-scale structure has virtually no effect on the confusion limit.
We present a study of NGC 4993, the host galaxy of the GW170817 gravitational wave event, the GRB170817A short gamma-ray burst (sGRB) and the AT2017gfo kilonova. We use Dark Energy Camera imaging, AAT spectra and publicly available data, relating our findings to binary neutron star (BNS) formation scenarios and merger delay timescales. NGC4993 is a nearby (40 Mpc) early-type galaxy, with $i$-band S\'ersic index $n=4.0$ and low asymmetry ($A=0.04\pm 0.01$). These properties are unusual for sGRB hosts. However, NGC4993 presents shell-like structures and dust lanes indicative of a recent galaxy merger, with the optical transient located close to a shell. We constrain the star formation history (SFH) of the galaxy assuming that the galaxy merger produced a star formation burst, but find little to no on-going star formation in either spatially-resolved broadband SED or spectral fitting. We use the best-fit SFH to estimate the BNS merger rate in this type of galaxy, as $R_{NSM}^{gal}= 5.7^{+0.57}_{-3.3} \times 10^{-6} {\rm yr}^{-1}$. If star formation is the only considered BNS formation scenario, the expected number of BNS mergers from early-type galaxies detectable with LIGO during its first two observing seasons is $0.038^{+0.004}_{-0.022}$, as opposed to $\sim 0.5$ from all galaxy types. Hypothesizing that the binary system formed due to dynamical interactions during the galaxy merger, the subsequent time elapsed can constrain the delay time of the BNS coalescence. By using velocity dispersion estimates and the position of the shells, we find that the galaxy merger occurred $t_{\rm mer}\lesssim 200~{\rm Myr}$ prior to the BNS coalescence.
Links to: arXiv, form interface, find, astro-ph, recent, 1710, contact, help (Access key information)
A matched filter technique is applied to the Planck all-sky Compton y-parameter map to measure the thermal Sunyaev-Zel'dovich (tSZ) effect produced by galaxy groups of different halo masses selected from large redshift surveys in the low-z Universe. Reliable halo mass estimates are available for all the groups, which allows us to bin groups of similar halo masses to investigate how the tSZ effect depends on halo mass over a large mass range. Filters are simultaneously matched for all groups to minimize projection effects. We find that the integrated y-parameter and the hot gas content it implies are consistent with the predictions of the universal pressure profile model only for massive groups above $10^{14}\,{\rm M}_\odot$, but much lower than the model prediction for low-mass groups. The halo mass dependence found is in good agreement with the predictions of a set of simulations that include strong AGN feedback, but simulations including only supernova feedback significantly over predict the hot gas contents in galaxy groups. Our results suggest that hot gas in galaxy groups is either effectively ejected or in phases much below the virial temperatures of the host halos.
Galactic outflows regulate the amount of gas galaxies convert into stars. However, it is difficult to measure the mass outflows remove because they span a large range of temperatures and phases. Here, we study the rest-frame ultraviolet spectrum of a lensed galaxy at z~2.9 with prominent interstellar absorption lines from O I, tracing neutral gas, up to O VI, tracing transitional phase gas. The O VI profile mimics weak low-ionization profiles at low velocities, and strong saturated profiles at high velocities. These trends indicate that O VI gas is co-spatial with the low-ionization gas. Further, at velocities blueward of -200 km/s the column density of the low-ionization outflow rapidly drops while the O VI column density rises, suggesting that O VI is created as the low-ionization gas is destroyed. Photoionization models do not reproduce the observed O VI, but adequately match the low-ionization gas, indicating that the phases have different formation mechanisms. Photoionized outflows are more massive than O VI outflows for most of the observed velocities, although the O VI mass outflow rate exceeds the photoionized outflow at velocities above the galaxy's escape velocity. Therefore, most gas capable of escaping the galaxy is in a hot outflow phase. We suggest that the O VI absorption is a temporary by-product of conduction transferring mass from the photoionized phase to an unobserved hot wind, and discuss how this mass-loading impacts the observed circum-galactic medium.
Using dark matter haloes identified in a large $N$-body simulation, we study halo assembly bias, with halo formation time, peak maximum circular velocity, concentration, and spin as the assembly variables. Instead of grouping haloes at fixed mass into different percentiles of each assembly variable, we present the joint dependence of halo bias on the {\it values} of halo mass and each assembly variable. In the plane of halo mass and one assembly variable, the joint dependence can be largely described as halo bias increasing outward from a global minimum. We find it unlikely to have a combination of halo variables to absorb all assembly bias effects. We then present the joint dependence of halo bias on two assembly variables at fixed halo mass. The gradient of halo bias does not necessarily follow the correlation direction of the two assembly variables and it varies with halo mass. Therefore in general for two correlated assembly variables one cannot be used as a proxy for the other in predicting halo assembly bias trend. Finally, halo assembly is found to affect the kinematics of haloes. Low-mass haloes formed earlier can have much higher pairwise velocity dispersion than those of massive haloes. In general, halo assembly leads to a correlation between halo bias and halo pairwise velocity distribution, with more strongly clustered haloes having higher pairwise velocity and velocity dispersion. However, the correlation is not tight, and the kinematics of haloes at fixed halo bias still depends on halo mass and assembly variables.
We examine the metallicity trends in the Milky Way (MW) bulge - using APOGEE DR13 data - and explore their origin by comparing two N-body models of isolated galaxies which develop a bar and a boxy/peanut (b/p) bulge. Both models have been proposed as scenarios for reconciling a disc origin of the MW bulge with a negative vertical metallicity gradient. The first model is a superposition of co-spatial, i.e. overlapping, disc populations with different scaleheights, kinematics and metallicities. In this model the thick, metal-poor, and centrally concentrated disc populations contribute significantly to the stellar mass budget in the inner galaxy. The second model is a single disc with an initial steep radial metallicity gradient, which is mapped by the bar into the b/p bulge in such a way that the vertical metallicity gradient of the MW bulge is reproduced -- as shown already in previous works in the literature. However, as we show here, the latter model does not reproduce the positive longitudinal metallicity gradient of the inner disc, nor the metal-poor innermost regions of the Bulge seen in the data. On the other hand, the model with co-spatial thin and thick disc populations reproduces all the aforementioned trends. We therefore see that it is possible to reconcile a (primarily) disc origin for the MW bulge with the observed trends in metallicity by mapping the inner thin and thick discs of the MW into a b/p. For this scenario to reproduce the observations, the $\alpha$-enhanced, metal-poor, thick disc populations must have a significant mass contribution in the inner regions -- as has been suggested for the Milky Way.
The total star formation budget of galaxies consists of the sum of the unobscured star formation, as observed in the rest-frame ultraviolet (UV), together with the obscured component that is absorbed and re-radiated by dust grains in the infrared. We explore how the fraction of obscured star formation depends on stellar mass for mass-complete samples of galaxies at $0<z<2.5$. We combine GALEX and WISE photometry for SDSS-selected galaxies with the 3D-HST treasury program and Spitzer/MIPS 24$\mu$m photometry in the well-studied 5 extragalactic CANDELS fields. We find a strong dependence of the fraction of obscured star formation ($f_{\mathrm{obscured}}$=SFR$_{\mathrm{IR}}$/SFR$_{\mathrm{UV+IR}}$) on stellar mass, with remarkably little evolution in this fraction with redshift out to $z$=2.5. 50\% of star formation is obscured for galaxies with log(M/M$_{\odot}$)=9.4; although unobscured star formation dominates the budget at lower masses, there exists a tail of low mass extremely obscured star-forming galaxies at $z>1$. For log(M/M$_{\odot}$)$>$10.5, $>$90\% of star formation is obscured at all redshifts. We also show that at fixed total SFR, $f_{\mathrm{obscured}}$ is lower at higher redshift. At fixed mass, high-redshift galaxies are observed to have more compact sizes and much higher star formation rates, gas fractions and hence surface densities (implying higher dust obscuration), yet we observe no redshift evolution in $f_{\mathrm{obscured}}$ with stellar mass. This poses a challenge to theoretical models to reproduce, where the observed compact sizes at high redshift seem in tension with lower dust obscuration.
We conducted high-sensitivity, high-velocity resolution CO(1-0) observations in a region containing a portion of the diffuse molecular cloud MBM 53 to determine whether weak CO emission was present. The results of our observations increase the amount of CO-detectable molecular gas in the region by a factor of two. The increased molecular mass for the cloud, if applicable to the molecular clouds in the entire Pegasus-Pisces region, decreases the dark molecular gas content from 58% of the total H2 mass to ~ 30%. If the results for MBM53 are applicable to other diffuse clouds, then the fraction of dark gas directly detectable via sensitive CO(1-0) observations in diffuse molecular clouds is similar to that predicted by models for Giant Molecular Clouds.
The Survey of Water and Ammonia in the Galactic Center (SWAG) covers the Central Molecular Zone (CMZ) of the Milky Way at frequencies between 21.2 and 25.4 GHz obtained at the Australia Telescope Compact Array at $\sim 0.9$ pc spatial and $\sim 2.0$ km s$^{-1}$ spectral resolution. In this paper, we present data on the inner $\sim 250$ pc ($1.4^\circ$) between Sgr C and Sgr B2. We focus on the hyperfine structure of the metastable ammonia inversion lines (J,K) = (1,1) - (6,6) to derive column density, kinematics, opacity and kinetic gas temperature. In the CMZ molecular clouds, we find typical line widths of $8-16$ km s$^{-1}$ and extended regions of optically thick ($\tau > 1$) emission. Two components in kinetic temperature are detected at $25-50$ K and $60-100$ K, both being significantly hotter than dust temperatures throughout the CMZ. We discuss the physical state of the CMZ gas as traced by ammonia in the context of the orbital model by Kruijssen et al. (2015) that interprets the observed distribution as a stream of molecular clouds following an open eccentric orbit. This allows us to statistically investigate the time dependencies of gas temperature, column density and line width. We find heating rates between $\sim 50$ and $\sim 100$ K Myr$^{-1}$ along the stream orbit. No strong signs of time dependence are found for column density or line width. These quantities are likely dominated by cloud-to-cloud variations. Our results qualitatively match the predictions of the current model of tidal triggering of cloud collapse, orbital kinematics and the observation of an evolutionary sequence of increasing star formation activity with orbital phase.
We present and analyse quasi-simultaneous multi-epoch spectral energy distributions (SEDs) of seven reverberation-mapped AGNs for which accurate black hole mass measurements and suitable archival data are available from the 'AGN Watch' monitoring programs. We explore the potential of optical-UV and X-ray data, obtained within 2 days, to provide more accurate SED-based measurements of individual AGN and quantify the impact of source variability on key measurements typically used to characterise the black hole accretion process plus on bolometric correction factors at 5100 {\AA}, 1350 {\AA} and for the 2-10 keV X-ray band, respectively. The largest SED changes occur on long timescales (>1 year). For our small sample, the 1 micron to 10 keV integrated accretion luminosity typically changes by 10% on short time-scales (over 20 days), by ~30% over a year, but can change by 100% or more for individual AGN. The EUV gap is the most uncertain part of the intrinsic SED, introducing a ~25% uncertainty in the accretion-induced luminosity, relative to the model independent interpolation method that we adopt. That aside, our analysis shows that the uncertainty in the accretion-induced luminosity, the Eddington luminosity ratio and the bolometric correction factors can be reduced (by a factor of two or more) by use of SEDs built from data obtained within 20 days. However, mass accretion rate and mass accretion efficiency are mostly limited by the unknown EUV emission and the unknown details of the central engine and our aspect angle.
We study the kinematics of stellar disc populations within the solar neighbourhood in order to find the imprints of the Galactic bar. We carry out the analysis by developing a numerical resolution of the 2D2V Collisionless Boltzmann Equation (CBE) and modelling the stellar motions within the plane of the Galaxy within the solar neighbourhood. We recover similar results to these obtained by different authors using N-body simulations, but we can also numerically identify faint structures thanks to the canceling of the Poisson noise. We find that the ratio of the bar pattern speed to the local circular frequency is in the range $\Omega_B/\Omega$ = 1.77 to 1.91. If the Galactic bar angle orientation is within the range from 24 to 45 degrees, the bar pattern speed is between 46 and 49 km/s/kpc.
The mid-infrared (IR) range contains many spectral features associated with large molecules and dust grains such as polycyclic aromatic hydrocarbons (PAHs) and silicates. These are usually very strong compared to fine-structure gas lines, and thus valuable in studying the spectral properties of faint distant galaxies. In this paper, we evaluate the capability of low-resolution mid-IR spectroscopic surveys of galaxies that could be performed by SPICA. The surveys are designed to address the question how star formation and black hole accretion activities evolved over cosmic time through spectral diagnostics of the physical conditions of the interstellar/circumnuclear media in galaxies. On the basis of results obtained with Herschel far-IR photometric surveys of distant galaxies and Spitzer and AKARI near- to mid-IR spectroscopic observations of nearby galaxies, we estimate the numbers of the galaxies at redshift z > 0.5, which are expected to be detected in the PAH features or dust continuum by a wide (10 deg^2) or deep (1 deg^2) blind survey, both for a given observation time of 600 hours. As by-products of the wide blind survey, we also expect to detect debris disks, through the mid-IR excess above the photospheric emission of nearby main-sequence stars, and we estimate their number. We demonstrate that the SPICA mid-IR surveys will efficiently provide us with unprecedentedly large spectral samples, which can be studied further in the far-IR with SPICA.
This is the first of a series of papers presenting the results from our survey of 25 Galactic globular clusters with the MUSE integral-field spectrograph. In combination with our dedicated algorithm for source deblending, MUSE provides unique multiplex capabilities in crowded stellar fields and allows us to acquire samples of up to 20 000 stars within the half-light radius of each cluster. The present paper focuses on the analysis of the internal dynamics of 22 out of the 25 clusters, using about 500 000 spectra of 200 000 individual stars. Thanks to the large stellar samples per cluster, we are able to perform a detailed analysis of the central rotation and dispersion fields using both radial profiles and two-dimensional maps. The velocity dispersion profiles we derive show a good general agreement with existing radial velocity studies but typically reach closer to the cluster centres. By comparison with proper motion data we derive or update the dynamical distance estimates to 14 clusters. Compared to previous dynamical distance estimates for 47 Tuc, our value is in much better agreement with other methods. We further find significant (>3sigma) rotation in the majority (13/22) of our clusters. Our analysis seems to confirm earlier findings of a link between rotation and the ellipticities of globular clusters. In addition, we find a correlation between the strengths of internal rotation and the relaxation times of the clusters, suggesting that the central rotation fields are relics of the cluster formation that are gradually dissipated via two-body relaxation.
Voronoi grids have been successfully used to represent density structures of gas in astronomical hydrodynamics simulations. While some codes are explicitly built around using a Voronoi grid, others, such as Smoothed Particle Hydrodynamics (SPH), use particle-based representations and can benefit from constructing a Voronoi grid for post-processing their output. So far, calculating the density of each Voronoi cell from SPH data has been done numerically, which is both slow and potentially inaccurate. This paper proposes an alternative analytic method, which is fast and accurate. We derive an expression for the integral of a cubic spline kernel over the volume of a Voronoi cell and link it to the density of the cell. Mass conservation is ensured rigorously by the procedure. The method can be applied more broadly to integrate a spherically symmetric polynomial function over the volume of a random polyhedron.
Links to: arXiv, form interface, find, astro-ph, recent, 1710, contact, help (Access key information)