X-ray emitting gaseous coronae around massive galaxies are a basic prediction of galaxy formation models. Although the coronae around spiral galaxies offer a fundamental test of these models, observational constraints on their characteristics are still scarce. While the presence of extended hot coronae has been established around a handful of massive spiral galaxies, the short X-ray observations only allowed to measure the basic characteristics of the coronae. In this work, we utilize deep XMM-Newton observations of NGC 6753 to explore its extended X-ray corona in unprecedented detail. Specifically, we establish the isotropic morphology of the hot gas, suggesting that it resides in hydrostatic equilibrium. The temperature profile of the gas shows a decrease with increasing radius: it drops from $kT\approx0.7$ keV in the innermost parts to $kT\approx0.4$ keV at 50 kpc radius. The temperature map reveals the complex temperature structure of the gas. We study the metallicity distribution of the gas, which is uniform at $Z\approx0.1$ Solar. This value is about an order of magnitude lower than that obtained for elliptical galaxies with similar dark matter halo mass, hinting that the hot gas in spiral galaxies predominantly originates from external gas inflows rather than from internal sources. By extrapolating the density profile of the hot gas out to the virial radius, we estimate the total gas mass and derive the total baryon mass of NGC 6753. We conclude that the baryon mass fraction is $f_{\rm b} \approx 0.06$, implying that about half of the baryons are missing.
Previous attempts at segmenting molecular line maps of molecular clouds have focused on using position-position-velocity data cubes of a single line to separate the spatial components of the cloud. In contrast, wide field spectral imaging with large spectral bandwidth in the (sub)mm domain now allows to combine multiple molecular tracers to understand the different physical and chemical phases that constitute giant molecular clouds. We aim at using multiple tracers (sensitive to different physical processes) to segment a molecular cloud into physically/chemically similar regions (rather than spatially connected components). We use a machine learning clustering method (the Meanshift algorithm) to cluster pixels with similar molecular emission, ignoring spatial information. Simple radiative transfer models are used to interpret the astrophysical information uncovered by the clustering. A clustering analysis based only on the J=1-0 lines of 12CO, 13CO and C18O reveals distinct density/column density regimes (nH~100, 500, and >1000 cm-3), closely related to the usual definitions of diffuse, translucent and high-column-density regions. Adding two UV-sensitive tracers, the (1-0) lines of HCO+ and CN, allows us to distinguish two clearly distinct chemical regimes, characteristic of UV-illuminated and UV-shielded gas. The UV-illuminated regime shows overbright HCO+ and CN emission, which we relate to photochemical enrichment. We also find a tail of high CN/HCO+ intensity ratio in UV-illuminated regions. Finer distinctions in density classes (nH~7E3, and 4E4 cm-3) for the densest regions are also identified, likely related to the higher critical density of the CN and HCO+ (1-0) lines. The association of simultaneous multi-line, wide-field mapping and powerful machine learning methods such as the Meanshift algorithm reveals how to decode the complex information available in molecular tracers.
We introduce Project MEGaSaURA: The Magellan Evolution of Galaxies Spectroscopic and Ultraviolet Reference Atlas. MEGaSaURA comprises medium-resolution, rest-frame ultraviolet spectroscopy of N=15 bright gravitationally lensed galaxies at redshifts of 1.68$<$z$<$3.6, obtained with the MagE spectrograph on the Magellan telescopes. The spectra cover the observed-frame wavelength range $3200 < \lambda_o < 8280$ \AA ; the verage spectral resolving power is R=3300. The median spectrum has a signal-to-noise ratio of $SNR=21$ per resolution element at 5000 \AA . As such, the MEGaSaURA spectra have superior signal-to-noise-ratio and wavelength coverage compared to what COS/HST provides for starburst galaxies in the local universe. This paper describes the sample, the observations, and the data reduction. We compare the measured redshifts for the stars, the ionized gas as traced by nebular lines, and the neutral gas as traced by absorption lines; we find the expected bulk outflow of the neutral gas, and no systemic offset between the redshifts measured from nebular lines and the redshifts measured from the stellar continuum. We provide the MEGaSaURA spectra to the astronomical community through a data release.
Spatially resolved kinematics of nearby galaxies has shown that the ratio of dynamical- to stellar population-based estimates of the mass of a galaxy (M$_*^{JAM}$ / M$_*$) correlates with $\sigma_e$, the light-weighted velocity dispersion within its half-light radius, if M$_*$ is estimated using the same Initial Mass Function (IMF) for all galaxies. This correlation may indicate that, in fact, the IMF is more bottom-heavy or dwarf-rich for galaxies with large $\sigma$. We use this correlation to estimate a dynamical or IMF-corrected stellar mass, M$_*^{\alpha JAM}$, from M$_*$ and $\sigma_e$ for a sample of 6x10$^5$ SDSS galaxies for which spatially resolved kinematics is not available. We also compute the `virial' mass estimate $k(n,R) R_e \sigma_R^2/G$, where $n$ is the Sersic index, in the SDSS and ATLAS-3D samples. We show that an $n$-dependent correction must be applied to the $k(n,R)$ values provided by Prugniel & Simien (1997). Our analysis also shows that the shape of the velocity dispersion profile in the ATLAS-3D sample varies weakly with $n$: $\sigma_R/\sigma_e=(R/R_e)^{-\gamma(n)}$. The resulting stellar mass functions, based on M$_*^{\alpha JAM}$ and the recalibrated virial mass, are in good agreement. If the M$_*^{\alpha JAM}$ / M$_*$ - $\sigma_e$ correlation is indeed due to the IMF, then our $\Phi$(M$_*^{\alpha JAM}$) is the first estimate of the stellar mass function in which $\sigma_e$-dependent variations in the IMF across the population have been accounted for. Using a Fundamental Plane based observational proxy for $\sigma_e$ produces comparable results. The use of direct measurements for estimating the IMF-dependent stellar mass is prohibitively expensive for a large sample of galaxies. Our analysis should enable a more accurate census of the mass in stars, especially at high redshift, at a fraction of the cost. Our results are provided in tabular form.
We construct $z\sim6-7$, 8, and 9 faint Lyman break galaxy samples (334, 61, and 37 galaxies, respectively) with accurate size measurements with the software $\texttt{glafic}$ from the complete Hubble Frontier Fields (FF) cluster and parallel fields data. These are the largest samples hitherto and reach down to the faint ends of recently obtained deep luminosity functions. At faint magnitudes, however, these samples are highly incomplete for galaxies with large sizes, implying that derivation of the luminosity function (LF) sensitively depends on the intrinsic size--luminosity (RL) relation. We thus conduct simultaneous maximum-likelihood estimation of LF and RL relation parameters from the observed distribution of galaxies on the RL plane with help of a completeness map as a function of size and luminosity. At $z\sim6-7$, we find that the intrinsic RL relation expressed as $r_\textrm{e} \propto L^\beta$ has a notably steeper slope of $\beta=0.46^{+0.08}_{-0.09}$ than those at lower redshifts, which in turn implies that the LF has a relatively shallow faint-end slope of $\alpha=-1.86^{+0.17}_{-0.18}$. This steep $\beta$ can be reproduced by a simple analytical model in which smaller galaxies have lower specific angular momenta. The $\beta$ and $\alpha$ values for the $z\sim8$ and 9 samples are consistent with those for $z\sim6-7$ but with larger errors. For all three samples there is a large, positive covariance between $\beta$ and $\alpha$, implying that the simultaneous determination of these two parameters is important. We also provide new strong lens mass models of Abell S1063 and Abell 370 as well as updated mass models of Abell 2744 and MACS J0416.1$-$2403.
We present the analysis of the XMM-Newton data of the Circum-Galactic Medium of MASsive Spirals (CGM-MASS) sample of six extremely massive spiral galaxies in the local Universe. All the CGM-MASS galaxies have diffuse X-ray emission from hot gas detected above the background extending $\sim(30-100)\rm~kpc$ from the galactic center. This doubles the existing detection of such extended hot CGM around massive spiral galaxies. The radial soft X-ray intensity profile of hot gas can be fitted with a $\beta$-function with the slope typically in the range of $\beta=0.35-0.55$. This range, as well as those $\beta$ values measured for other massive spiral galaxies, including the Milky Way (MW), are in general consistent with X-ray luminous elliptical galaxies of similar hot gas luminosity and temperature, and with those predicted from a hydrostatic isothermal gaseous halo. Hot gas in such massive spiral galaxy tends to have temperature comparable to its virial value, indicating the importance of gravitational heating. This is in contrast to lower mass galaxies where hot gas temperature tends to be systematically higher than the virial one. The ratio of the radiative cooling to free fall timescales of hot gas is much larger than the critical value of $\sim10$ throughout the entire halos of all the CGM-MASS galaxies, indicating the inefficiency of gas cooling and precipitation in the CGM. The hot CGM in these massive spiral galaxies is thus most likely in a hydrostatic state, with the feedback material mixed with the CGM, instead of escaping out of the halo or falling back to the disk. We also homogenize and compare the halo X-ray luminosity measured for the CGM-MASS galaxies and other galaxy samples and discuss the "missing" galactic feedback detected in these massive spiral galaxies.
Linearly polarized emission is described, in general, in terms of the Stokes parameters $Q$ and $U$, from which the polarization intensity and polarization angle can be determined. Although the polarization intensity and polarization angle provide an intuitive description of the polarization, they are affected by the limitations of interferometric data, such as missing single-dish data in the u-v plane, from which radio frequency interferometric data is visualized. To negate the effects of these artefacts, it is desirable for polarization diagnostics to be rotationally and translationally invariant in the $Q$-$U$ plane. One rotationally and translationally invariant quantity, the polarization gradient, has been shown to provide a unique view of spatial variations in the turbulent interstellar medium when applied to diffuse radio frequency synchrotron emission. In this paper we develop a formalism to derive additional rotationally and translationally invariant quantities. We present new diagnostics that can be applied to diffuse or point-like polarized emission in any waveband, including a generalization of the polarization gradient, the polarization directional curvature, polarization wavelength derivative, and polarization wavelength curvature. In Paper II we will apply these diagnostics to observed and simulated images of diffuse radio frequency synchrotron emission.
We study the dynamical state of cores by using a simple analytical model, a
sample of observational massive cores, and numerical simulations of collapsing
massive cores. From the model, we find that, if cores are formed from turbulent
compressions, they evolve from small to large column densities, increasing
their velocity dispersion as they collapse, while they tend to equipartition
between gravity and kinetic energy.
From the observational sample, we find that: (a) cores with substantially
different column densities in the sample do not follow a Larson-like
linewidth-size relation. Instead, cores with higher column densities tend to be
located in the upper-left corner of the Larson velocity dispersion-size
diagram, a result predicted previously (Ballesteros-Paredes et al. 2011a). (b)
The data exhibit cores with overvirial values.
Finally, in the simulations we reproduce the behavior depicted by the model
and observational sample: cores evolve towards larger velocity dispersions and
smaller sizes as they collapse and increase their column density. However,
collapsing cores appear to approach overvirial states within a free-fall time.
The cause of this apparent excess of kinetic energy is an underestimation of
the gravitational energy, due to the assumption that the gravitational energy
is given by the energy of an isolated sphere of constant column density. This
excess disappears when the gravitational energy is correctly calculated from
the actual spatial mass distribution, where inhomogeneities, as well as the
potential due to the mass outside of the core, also contribute to the
gravitational energy. We conclude that the observed energy budget of cores in
surveys is consistent with their non-thermal motions being driven by their
self-gravity and in a hierarchical and chaotic collapse.
As part of our ongoing HII Region Discovery Survey (HRDS), we report the Green Bank Telescope detection of 148 new angularly-large Galactic HII regions in radio recombination line (RRL) emission. Our targets are located at a declination greater than -45deg., which corresponds to 266deg. > l > -20deg. at b = 0deg. All sources were selected from the WISE Catalog of Galactic HII Regions, and have infrared angular diameters >260''. The Galactic distribution of these "large" HII regions is similar to that of the previously-known sample of Galactic HII regions. The large HII region RRL line width and peak line intensity distributions are skewed toward lower values compared with that of previous HRDS surveys. We discover 7 sources with extremely narrow RRLs <10 km/s. If half the line width is due to turbulence, these 7 sources have thermal plasma temperatures <1100 K. These temperatures are lower than any measured for Galactic HII regions, and the narrow line components may arise instead from partially ionized zones in the HII region photo-dissociation regions. We discover G039.515+00.511, one of the most luminous HII regions in the Galaxy. We also detect the RRL emission from three HII regions with diameters >100 pc, making them some of the physically largest known HII regions in the Galaxy. This survey completes the HRDS HII region census in the Northern sky, where we have discovered 887 HII regions and more than doubled the previously-known census of Galactic HII regions.
By performing an orbital analysis within a Galactic model including a bar, we found that it is plausible that the newly discovered stars that show enhanced Al and N accompanied by Mg underabundances may have formed in the outer halo, or were brought in by satellites field possibly accreted a long time ago. However, another subsample of three N- and Al-rich stars with Mg-deficiency are kinematically consistent with the inner stellar halo. A speculative scenario to explain the origin of the atypical chemical composition of these stars in the inner halo is that they migrated to the inner stellar halo as unbound stars due to the mechanism of bar-induced resonant trapping.
[Abridged] The Lupus I cloud is found between the Upper-Scorpius and the Upper-Centaurus-Lupus sub-groups, where the expanding USco HI shell appears to interact with a bubble currently driven by the winds of the remaining B-stars of UCL. We investigate if the Lupus I molecular could have formed in a colliding flow, and how the kinematics of the cloud might have been influenced by the larger scale gas dynamics. We performed APEX 13CO and C18O observations of three parts of Lupus. We compare these results to the atomic hydrogen data from the GASS HI survey and our dust emission results presented in the previous paper. Based on the velocity information, we present a geometric model for the interaction zone between the USco shell and the UCL wind bubble. We present evidence that the molecular gas of Lupus I is tightly linked to the atomic material of the USco shell. The CO emission in Lupus I is found mainly at velocities in the same range as the HI velocities. Thus, the molecular cloud is co-moving with the expanding USco atomic Hi shell. The gas in the cloud shows a complex kinematic structure with several line-of-sight components that overlay each other. The non-thermal velocity dispersion is in the transonic regime in all parts of the cloud and could be injected by external compression. Our observations and the derived geometric model agree with a scenario where Lupus I is located in the interaction zone between the USco shell and the UCL wind bubble. The kinematics observations are consistent with a scenario where the Lupus I cloud formed via shell instabilities. The particular location of Lupus I between USco and UCL suggests that counter-pressure from the UCL wind bubble and pre-existing density enhancements, perhaps left over from the gas stream that formed the stellar subgroups, may have played a role in its formation.
Local Group Analogs (LGA) are galaxy associations dominated by few bright Spirals, reminiscent of the LG. The NGC3447/NGC3447A system, member of the LGG 225 group, a nearby LGA, is considered a physical pair: an intermediate luminosity late type spiral, NGC3447, and an irregular companion, NGC3447A, linked by a faint filament of matter. A ring-like structure in the NGC3447 outskirts is emphasised by UV observations. This work aims to contribute to the understanding of galaxy evolution in low density environments, favourable habitat to highly effective encounters. We performed a multi-wavelength analysis of the surface photometry of this system to derive spectral energy distribution and structural properties using UV and optical images. We also characterised the velocity field of the pair using new kinematic observations. All these data are used to constrain smooth particle hydrodynamic simulations with chemo-photometric implementation to shed light on the evolution of this system. Luminosity profiles are all consistent with the presence of a disc extending and including NGC3447A. The overall velocity field does not emphasise any significant rotation pattern, rather a small velocity gradient between NGC3447 and NGC3447A. Our simulation, detached from a large grid explored to best-fit the global properties of the system, suggests that this arises from an encounter between two halos of equal mass. NGC3447 and NGC3447A belong to the same halo, NGC3447A being a substructure of the same disk as NGC3447. The halo gravitational instability, enhanced by the encounter, fuels a long lived instability in this dark matter dominated disk, driving its morphology. This system may warn about a new class of "false pairs" and the potential danger of a misunderstanding of such objects in pair surveys that could produce a severe underestimate of the total mass of the system. (abridged)
The star-forming ability of a molecular cloud depends on the fraction of gas it can cycle into the dense-phase. Consequently, one of the crucial questions in reconciling star-formation in clouds is to understand the factors that control this process. While it is widely accepted that the variation in ambient conditions can alter significantly the ability of a cloud to spawn stars, the observed variation in the star-formation rate in nearby clouds that experience similar ambient conditions, presents an interesting question. In this work we attempted to reconcile this variation within the paradigm of colliding flows. To this end we develop self-gravitating, hydrodynamic realisations of identical flows, but allowed to collide off-centre. Typical observational diagnostics such as the gas-velocity dispersion, the fraction of dense-gas, the column density distribution ({\small N-PDF}), the distribution of gas mass as a function of $K$-band extinction and the strength of compressional/solenoidal modes in the post-collision cloud were deduced for different choices of the impact parameter of collision. We find that a strongly sheared cloud is terribly inefficient in cycling gas into the dense phase and that such a cloud can possibly reconcile the sluggish nature of star-formation reported for some clouds. Within the paradigm of cloud-formation via colliding flows this is possible in case of flows colliding with a relatively large impact parameter. We conclude that compressional modes - though probably essential - are insufficient to ensure a relatively higher star-formation efficiency in a cloud.
We stack the rest-frame ultraviolet spectra of N=14 highly magnified gravitationally lensed galaxies at redshifts 1.6<z<3.6. The resulting new composite spans $900< \lambda_{rest} < 3000$ \AA, with a peak signal-to-noise ratio of 103 per spectral resolution element ($\sim$100 km/s). It is the highest signal-to-noise ratio, highest spectral resolution composite spectrum of $z\sim2$--3 galaxies yet published. The composite reveals numerous weak nebular emission lines and stellar photospheric absorption lines that can serve as new physical diagnostics, particularly at high redshift with the James Webb Space Telescope (JWST). We report equivalent widths to aid in proposing for and interpreting JWST spectra. We examine the velocity profiles of strong absorption features in the composite, and in a matched composite of $z\sim0$ COS/HST galaxy spectra. We find remarkable similarity in the velocity profiles at $z\sim 0$ and $z\sim2$, suggesting that similar physical processes control the outflows across cosmic time. While the maximum outflow velocity depends strongly on ionization potential, the absorption-weighted mean velocity does not. As such, the bulk of the high-ionization absorption traces the low-ionization gas, with an additional blueshifted absorption tail extending to at least $-2000$ km/s . We interpret this tail as arising from the stellar wind and photospheres of massive stars. Starburst99 models are able to replicate this high-velocity absorption tail. However, these theoretical models poorly reproduce several of the photospheric absorption features, indicating that improvements are needed to match observational constraints on the massive stellar content of star-forming galaxies at $z \sim 2$. We publicly release our composite spectra.
Supernovae represent significant sources of dust in the interstellar medium. In this work, deep far-infrared (FIR) observations of the Crab Nebula are studied to provide a new and reliable constraint on the amount of dust present in this supernova remnant. Deep exposures between 70 and 500 $\mu$m taken by PACS and SPIRE instruments on-board the Herschel Space Telescope, compiling all observations of the nebula including PACS observing mode calibration, are refined using advanced processing techniques, thus providing the most accurate data ever generated by Herschel on the object. We carefully find the intrinsic flux of each image by masking the source and creating a 2D polynomial fit to deduce the background emission. After subtracting the estimated non-thermal synchrotron component, two modified blackbodies were found to best fit the remaining infrared continuum, the cold component with T$_c$ = 8.3 $\pm$ 3.0 K and M$_d$ = 0.27 $\pm$ 0.05 M$_{\odot}$ and the warmer component with T$_w$ = 27.2 $\pm$ 1.3 K and M$_d$ = (1.3 $\pm$ 0.4) $\times$10$^{-3}$ M$_{\odot}$.
We present new evidence for AGN feedback in a subset of 69 quenched low-mass galaxies ($M_{\star} \lesssim 5\times10^{9}$ M$_{\odot}$, $M_{\rm{r}} > -19$) selected from the first two years of the SDSS-IV MaNGA survey. The majority (85 per cent) of these quenched galaxies appear to reside in a group environment. We find 6 galaxies in our sample that appear to have an active AGN that is preventing on-going star-formation; this is the first time such a feedback mechanism has been observed in this mass range. Interestingly, five of these six galaxies have an ionised gas component that is kinematically offset from their stellar component, suggesting the gas is either recently accreted or outflowing. We hypothesise these six galaxies are low-mass equivalents to the "red geysers" observed in more massive galaxies. Of the other 62 galaxies in the sample, we find 8 do appear for have some low-level, residual star formation, or emission from hot, evolved stars. The remaining galaxies in our sample have no detectable ionised gas emission throughout their structures, consistent with them being quenched. This work shows the potential for understanding the detailed physical properties of dwarf galaxies through spatially resolved spectroscopy.
We investigate radial gradients in the recent star formation history (SFH) using 1917 galaxies with $0.01<z<0.14$ and integral-field spectroscopy from the ongoing MaNGA survey. For each galaxy, we obtain two-dimensional maps and radial profiles for three spectroscopically-measured parameters that are sensitive to the recent SFH: D$_n$(4000) (the 4000\AA\ break), EW(H$\delta_A$) (equivalent width of the H$\delta$ absorption line), and EW(H$\alpha$) (equivalent width of the H$\alpha$ emission line). We find the majority of the spaxels in these galaxies are consistent with models with continuously declining star formation rate, indicating that starbursts occur rarely in local galaxies. We classify the galaxies into three classes: fully star-forming (SF), partly quenched (PQ) and totally quenched (TQ), according to the fraction of quenched area within 1.5 times the effective radius, $f_Q(1.5R_e)$. We find that galaxies less massive than $10^{10}$M$_{\odot}$ present at most weak radial gradients in all the diagnostic parameters. In contrast, galaxies with stellar mass above $10^{10}$M$_{\odot}$ present significant gradients in all the three diagnostic parameters if classified as SF or PQ, but show weak gradients in D$_n$(4000) and EW(H$\delta_A$) and no gradients in EW(H$\alpha$) if in the TQ class. This implies the existence of a critical stellar mass, at $\sim10^{10}$M$_{\odot}$, above which the star formation in a galaxy gets shutdown from the inside out. Galaxies tend to evolve synchronously from inner to outer regions before their mass reaches the critical value. We have further divided the sample at fixed mass by both bulge-to-total luminosity ratio and morphological type, finding that our conclusions hold regardless of these factors: it appears that the presence of a central dense object is not a driving parameter, but rather a byproduct of the star formation cessation process.
We analyze methanol excitation in the absence of external radiation and consider LTE methods for probing interstellar gas. We show that rotation diagrams correctly estimate the gas kinetic temperature only if they are built from lines with the upper levels located in the same K-ladders, such as the J_0-J_{-1}E lines at 157~GHz, the J_1-J_0E lines at 165~GHz or the J_2-J_1E lines at 25~GHz. The gas density should be no less than 10^7~cm^{-3}. Rotation diagrams built from lines with different K values of the upper levels (2_K-1_K at 96~GHz, 3_K-2_K at 145~GHz, or 5_K-4_K at 241~GHz) significantly underestimate the temperature but allow a density estimation. In addition, the diagrams based on the 2_K-1_K lines make possible methanol column density estimates within a factor of about 2--5. We suggest that rotation diagrams should be used in the following manner. First, one should build two rotation diagrams, one from the lines at 96, 145, or 241~GHz, and another from the lines at 157, 165, or 25~GHz. The former diagram is used to estimate the gas density. If the density is about 10^7~cm^{-3} or higher, the latter diagram reproduces the temperature fairly well. If the density is around 10^6~cm^{-3}, the temperature obtained from the latter diagram should be multiplied by a factor of 1.5--2. If the density is about 10^5~cm^{-3} or lower, then the latter diagram yields a temperature that is lower than the kinetic temperature by a factor of three or larger and should be used only as a lower limit on the kinetic temperature. Errors of methanol column density determined from the integrated intensity of a single line may be larger than an order of magnitude even when the gas temperature is well-known. However, if the J_0-(J-1)_0E lines, as well as the J_1-(J-1)_1A^{+} or A^{-} lines are used, the relative error of the column density proves to be no larger than several units.
Radiation in the extreme ultraviolet (EUV) and soft X-ray holds clues to the location of the missing baryons, the energetics in stellar feedback processes, and the cosmic enrichment history. Additionally, EUV and soft X-ray photons help determine the ionization state of most intergalactic and circumgalactic metals, shaping the rate at which cosmic gas cools. Unfortunately, this band is extremely difficult to probe observationally due to absorption from the Galaxy. In this paper, we model the contributions of various sources to the cosmic EUV and soft X-ray backgrounds. We bracket the contribution from (1) quasars, (2) X-ray binaries, (3) hot interstellar gas, (4) circumgalactic gas, and (5) virialized gas, developing models that extrapolate into these bands using both empirical and theoretical inputs. While quasars are traditionally assumed to dominate these backgrounds, we discuss the substantial uncertainty in their contribution. Furthermore, we find that hot intrahalo gases likely emit an O(1) fraction of this radiation at low redshifts, and that interstellar and circumgalactic emission potentially contribute tens of percent to these backgrounds at all redshifts. We estimate that uncertainties in the angular-averaged background intensity impact the ionization corrections for common circumgalactic and intergalactic metal absorption lines by ~0.3-1 dex, and we show that local emissions are comparable to the cosmic background only at r_prox = 10-100 kpc from Milky Way-like galaxies.
Links to: arXiv, form interface, find, astro-ph, recent, 1710, contact, help (Access key information)
The reionization of cosmic hydrogen marks a critical juncture in the history of structure formation in the universe. Here we present a new formulation of the standard reionization equation for the evolution of the volume-averaged HII fraction that is more consistent with the accepted conceptual model of inhomogeneous intergalactic absorption. The revised equation retains the basic terminology and simplicity of the classic calculation but explicitly accounts for the presence of the optically thick "Lyman-limit systems" that are known to determine the mean free path of ionizing radiation after overlap. Integration of this equation provides a better characterization of the timing of reionization by smoothly linking the pre-overlap with the post-overlap phases of such process. We confirm the validity of the quasi-instantaneous approximation as predictor of reionization completion/maintenance, and discuss new insights on the sources of cosmic reionization using the improved formalism. A constant emission rate into the intergalactic medium (IGM) of 3 Lyman continuum (LyC) photons per atom per Gyr leads to a reionization history that is consistent with a number of observational constraints on the ionization state of the z=5-9 universe and with the reduced Thomson scattering optical depth recently reported by the Planck Collaboration. While star-forming galaxies can dominate the reionization process if the luminosity-weighted fraction of LyC photons that escape into the IGM, f_esc, exceeds 15% (for a faint magnitude cut-off of the galaxy UV luminosity function of M_lim=-13 and a LyC photon yield per unit 1500 AA luminosity of xi_ion=10^{25.3} Hz/erg, simple models where the product of the two unknowns f_esc xi_ion is not evolving with redshift fail to reproduce the changing neutrality of the IGM observed at these epochs.
We study the connection between the observed star formation rate-stellar mass (SFR-$M_*$) relation and the evolution of the stellar mass function (SMF) by means of a Subhalo Abundance Matching technique coupled to merger trees extracted from a N-body simulation. Our approach consist of forcing the model to match the observed SMF at redshift $z \sim 2.3$, and let it evolve down to $z \sim 0.3$ according to a $\tau$ model, an exponentially declining functional form which describes the star formation rate decay of both satellite and central galaxies. In this study, we use three different sets of SMFs: ZFOURGE data from \cite{tomczak14}; UltraVISTA data from \cite{ilbert13} and COSMOS data from \cite{davidzon17}. We also build a mock survey combining UltraVISTA with ZFOURGE. Our modelling of quenching time scales is consistent with the evolution of the SMF down to $z \sim 0.3$, with different accuracy depending on the particular survey used for calibration. We tested our model against the observed SMFs at low redshift and it predicts residuals (observation-model) within $1\sigma$ observed scatter along most of the stellar mass range investigated, and with mean residuals below 0.1 dex in the range $\sim [10^{8.7}-10^{11.7}] M_{\odot}$. We then compare the SFR-$M_*$ relation predicted by the model with the observed one at different redshifts. The predicted SFR-$M_*$ relation underpredicts the median SFR at fixed stellar mass relative to observations at all redshifts. Nevertheless, the shapes are consitent with the observed relations up to intermediate mass galaxies, followed by a rapid decline for massive galaxies.
The spatial distribution of satellite galaxies around pairs of galaxies in the Sloan Digital Sky Survey (SDSS) have been found to bulge significantly towards the respective partner. Highly anisotropic, planar distributions of satellite galaxies are in conflict with expectations derived from cosmological simulations. Does the lopsided distribution of satellite systems around host galaxy pairs constitute a similar challenge to the standard model of cosmology? We investigate whether such satellite distributions are present around stacked pairs of hosts extracted from the $\Lambda$CDM simulations Millennium-I, Millennium-II, ELVIS, and Illustris-1. By utilizing this set of simulations covering different volumes, resolutions, and physics, we implicitly test whether a lopsided signal exists for different ranges of satellite galaxy masses, and whether the inclusion of hydrodynamical effects produces significantly different results. All simulations display a lopsidedness similar to the observed situation. The signal is highly significant for simulations containing a sufficient number of hosts and resolved satellite galaxies (up to $5\,\sigma$ for Millennium-II). We find a projected signal that is up to twice as strong as that reported for the SDSS systems for certain opening angles ($\sim16\%$ more satellites in the direction between the pair than expected for uniform distributions). Considering that the SDSS signal is a lower limit owing to likely back- and foreground contamination, the $\Lambda$CDM simulations appear to be consistent with this particular empirical property of galaxy pairs.
Active galactic nuclei (AGN) influence the formation and evolution of galaxies and their environments. Since the large scale environment can also affect AGN, this work studies how their formation and properties depend on the environment. We use a reconstructed three-dimensional high-resolution density field obtained from a Bayesian large-scale structure reconstruction method applied to the 2M++ galaxy sample. A web-type classification relying on the shear tensor is used to identify different structures on the cosmic web, defining voids, sheets, filaments, and clusters. We confirm that the environmental density affects the AGN formation and their properties. We found that the AGN clustering is equivalent to the galaxy clustering, indicating that active and non-active galaxies reside in similar dark matter halos. However, the clustering and the occurrence rate are different for each spectral type and accretion rate. These differences are consistent with the AGN evolutionary sequence suggested by previous authors, Seyferts and Transition objects transforming into LINERs (Low-Ionization Nuclear Emission Line Regions), the weaker counterpart of Seyferts. AGN properties show correlations with the environmental density that can be explained by mergers and interactions at high densities. More powerful starbursts and younger stellar populations are found in high densities, where interactions and mergers are more likely. AGN hosts show smaller masses in clusters for Seyferts and Transition objects, which might be due to gas stripping. In voids, the AGN population is dominated by the most massive galaxy hosts.
The radio source J1819+3845 underwent a period of extreme interstellar scintillation between circa 1999 and 2007. The plasma structure responsible for this scintillation was determined to be just $1$-$3\,$pc from the solar system and to posses a density of $n_e\sim 10^2\,$cm$^{-3}$ that is three orders of magnitude higher than the ambient interstellar density (de Bruyn & Macquart 2015). Here we present radio-polarimetric images of the field towards J1819+3845 at wavelengths of 0.2, 0.92 and 2$\,$m. We detect an elliptical plasma globule of approximate size $1^\circ \times \gtrsim 2^\circ$ (major-axis position angle of $\approx -40^\circ$), via its Faraday-rotation imprint ($\approx 15\,$rad$\,$m$^{-2}$) on the diffuse Galactic synchrotron emission. The extreme scintillation of J1819+3845 was most likely caused at the turbulent boundary of the globule (J1819+3845 is currently occulted by the globule). The origin and precise nature of the globule remain unknown. Our observations are the first time plasma structures that likely cause extreme scintillation have been directly imaged.
Using infrared data from the Herschel Space Observatory and Karl G. Jansky Very Large Array (VLA) 3 GHz observations in the COSMOS field, we investigate the redshift evolution of the infrared-radio correlation (IRRC) for star-forming galaxies (SFGs) we classify as either spheroid- or disc-dominated based on their morphology. The sample predominantly consists of disc galaxies with stellar mass ${\gtrsim}10^{10}\,M_{\odot}$, and residing on the star-forming main sequence (MS). After the removal of AGN using standard approaches, we observe a significant difference between the redshift-evolution of the median IR/radio ratio $\overline{q}_{\mathrm{TIR}}$ of (i) a sample of ellipticals, plus discs with a substantial bulge component (`spheroid-dominated' SFGs) and, (ii) virtually pure discs and irregular systems (`disc-dominated' SFGs). The spheroid-dominated population follows a declining $\overline{q}_{\mathrm{TIR}}$ vs. $z$ trend similar to that measured in recent evolutionary studies of the IRRC. However, for disc-dominated galaxies, where radio and IR emission should be linked to star formation in the most straightforward way, we measure very little change in $\overline{q}_{\mathrm{TIR}}$. This suggests that low-redshift calibrations of radio emission as an SFR-tracer may remain valid out to at least $z\,{\simeq}\,1\,{-}\,1.5$ for pure star-forming systems. We find that the different redshift-evolution of $q_{\rm TIR}$ for the spheroid- and disc-dominated sample is mainly due to an increasing radio excess for spheroid-dominated galaxies at $z\,{\gtrsim}\,$0.8, hinting at some residual AGN activity in these systems. This finding demonstrates that in the absence of AGN the IRRC is independent of redshift, and that radio observations can therefore be used to estimate SFRs at all redshifts for genuinely star-forming galaxies.
We present spatially resolved stellar kinematic maps, for the first time, for a sample of 17 intermediate redshift galaxies (0.2 < z < 0.8). We used deep MUSE/VLT integral field spectroscopic observations in the Hubble Deep Field South (HDFS) and Hubble Ultra Deep Field (HUDF), resulting from ~30h integration time per field, each covering 1'x1' field of view, with ~0.65" spatial resolution. We selected all galaxies brighter than 25mag in the I band and for which the stellar continuum is detected over an area that is at least two times larger than the spatial resolution. The resulting sample contains mostly late-type disk, main-sequence star-forming galaxies with 10^8.5 - 10^10.5 Msun. Using a full-spectrum fitting technique, we derive two-dimensional maps of the stellar and gas kinematics, including the radial velocity V and velocity dispersion sigma. We find that most galaxies in the sample are consistent with having rotating stellar disks with roughly constant velocity dispersions and that the Vrms=sqrt{V^2+sigma^2} of the gas and stars, a scaling proxy for the galaxy gravitational potential, compare well to each other. These spatially resolved observations of intermediate redshift galaxies suggest that the regular stellar kinematics of disk galaxies that is observed in the local Universe was already in place 4 - 7 Gyr ago and that their gas kinematics traces the gravitational potential of the galaxy, thus is not dominated by shocks and turbulent motions. Finally, we build dynamical axisymmetric Jeans models constrained by the derived stellar kinematics for two specific galaxies and derive their dynamical masses. These are in good agreement (within 25%) with those derived from simple exponential disk models based on the gas kinematics. The obtained mass-to-light ratios hint towards dark matter dominated systems within a few effective radii.
We mapped the filamentary infrared dark cloud (IRDC) G304.74+01.32 at 350 $\mu$m with the SABOCA bolometer. The new SABOCA data have a factor of 2.2 times higher resolution than our previous LABOCA 870 $\mu$m map of the cloud. We also employed the Herschel far-IR and submillimetre, and WISE IR data available for G304.74. The SABOCA data show that G304.74 is composed of a dense filamentary structure with a mean width of only $0.18\pm0.05$ pc. The percentage of LABOCA clumps that are found to be fragmented into SABOCA cores is $36\%\pm16\%$. The WISE data suggest that $65\%\pm18\%$ of the SABOCA cores host young stellar objects (YSOs). The mean dust temperature of the clumps, derived by comparing the Herschel/SPIRE flux densities, was found to be $15.0 \pm 0.8$ K. The mean mass, beam-averaged H$_2$ column density, and H$_2$ number density of the LABOCA clumps are estimated to be $55\pm10$ M$_{\odot}$, $(2.0\pm0.2)\times10^{22}$ cm$^{-2}$, and $(3.1\pm0.2)\times10^4$ cm$^{-3}$. The corresponding values for the SABOCA cores are $29\pm3$ M$_{\odot}$, $(2.9\pm0.3)\times10^{22}$ cm$^{-2}$, and $(7.9\pm1.2)\times10^4$ cm$^{-3}$. The G304.74 filament is estimated to be thermally supercritical by a factor of $\gtrsim3.5$ on the scale probed by LABOCA, and by a factor of $ \gtrsim1.5$ for the SABOCA filament. Our data strongly suggest that G304.74 has undergone hierarchical fragmentation. The IRDC G304.74 has a seahorse-like morphology in the Herschel images, and the filament appears to be attached by elongated, perpendicular striations. Besides the presence of perpendicularly oriented, dusty striations and potential embedded intermediate-mass YSOs, G304.74 is a relatively nearby ($d\sim2.5$ kpc) IRDC, which makes it a useful target for future star formation studies. Owing to its observed morphology, we propose that G304.74 could be nicknamed the Seahorse Nebula.
We show that the mass of a dark matter halo can be inferred from the dynamical status of its satellite galaxies. Using 9 dark-matter simulations of halos like the Milky Way (MW), we find that the present-day substructures in each halo follow a characteristic distribution in the phase space of orbital binding energy and angular momentum, and that this distribution is similar from halo to halo but has an intrinsic dependence on the halo formation history. We construct this distribution directly from the simulations for a specific halo and extend the result to halos of similar formation history but different masses by scaling. The mass of an observed halo can then be estimated by maximizing the likelihood in comparing the measured kinematic parameters of its satellite galaxies with these distributions. We test the validity and accuracy of this method with mock samples taken from the simulations. Using the positions, radial velocities, and proper motions of 9 tracers and assuming observational uncertainties comparable to those of MW satellite galaxies, we find that the halo mass can be recovered to within $\sim$40%. The accuracy can be improved to within $\sim$25% if 30 tracers are used. However, the dependence of the phase-space distribution on the halo formation history sets a minimum uncertainty of $\sim$20% that cannot be reduced by using more tracers. We believe that this minimum uncertainty also applies to any mass determination for a halo when the phase space information of other kinematic tracers is used.
We present the public release of the MultiDark-Galaxies: three distinct galaxy catalogues derived from one of the Planck cosmology MultiDark simulations (i.e. MDPL2, with a volume of (1 Gpc/$h$)$^{3}$ and mass resolution of $1.5 \times 10^{9} M_{\odot}/h$) by applying the semi-analytic models GALACTICUS, SAG, and SAGE to it. We compare the three models and their conformity with observational data for a selection of fundamental properties of galaxies like stellar mass function, star formation rate, cold gas fractions, and metallicities - noting that they sometimes perform differently reflecting model designs and calibrations. We have further selected galaxy subsamples of the catalogues by number densities in stellar mass, cold gas mass, and star formation rate in order to study the clustering statistics of galaxies. We show that despite different treatment of orphan galaxies, i.e. galaxies that lost their dark-matter host halo due to the finite mass resolution of the N-body simulation or tidal stripping, the clustering signal is comparable, and reproduces the observations in all three models - in particular when selecting samples based upon stellar mass. Our catalogues provide a powerful tool to study galaxy formation within a volume comparable to those probed by on-going and future photometric and redshift surveys. All model data consisting of a range of galaxy properties - including broad-band SDSS magnitudes - are publicly available.
Event Horizon Telescope (EHT) observations at 230 GHz are combined with Very Long Baseline Interferometry (VLBI) observations at 86 GHz and high resolution Hubble Space Telescope optical observations in order to constrain the broadband spectrum of the emission from the base of the jet in M87. The recent VLBI observations of Hada et al provide much stricter limits on the 86 GHz luminosity and component acceleration in the jet base than was available to previous modelers. They reveal an almost hollow jet on sub-mas scales. Thus, tubular models of the jet base emanating from the innermost accretion disk are considered within the region responsible for the EHT correlated flux. There is substantial synchrotron self absorbed opacity at 86 GHz. A parametric analysis indicates that the jet dimensions and power depend strongly on the 86 GHz flux density and the black hole spin, but weakly on other parameters such as jet speed, 230 GHz flux density and optical flux. The entire power budget of the M87 jet, $\lesssim 10^{44}\rm{ergs/sec}$, can be accommodated by the tubular jet. No invisible, powerful spine is required. Even though this analysis never employs the resolution of the EHT, the spectral shape implies a dimension transverse to the jet direction of 12-21 $\mu \rm{as}$ ($\sim$24-27$\, \mu \rm{as}$) for $0.99 > a/M > 0.95$ ($a/M\sim 0.7 $), where $M$ is the mass and $a$ is the angular momentum per unit mass of the central black hole.
We present MUSE observations of the field of the quasar Q0152$-$020 whose spectrum shows a Lyman limit system (LLS) at redshift $z_{\rm abs} = 0.38$, with a metallicity Z $\gtrsim 0.06$ Z$_\odot$. The low ionization metal lines associated with the LLS present two narrow distinct absorption components with a velocity separation of 26 km ${\rm s}^{-1}$. We detect six galaxies within 600 km ${\rm s}^{-1}$ from the absorption redshift; their projected distances from the quasar sightline range from 60 to 200 kpc. The optical spectra of five of these galaxies exhibit prominent nebular emission lines, from which we deduce extinction-corrected star formation rates in the range SFR = 0.06-1.3 M$_\odot$~yr$^{-1}$, and metallicities between 0.2 Z$_\odot$ and Z$_\odot$. The sixth galaxy is only detected in the stellar continuum. By combining our data with archival Keck/HIRES spectroscopy of the quasar and HST/WFPC2 imaging of the field, we can relate absorption line and galaxy kinematics; we conclude that the LLS is most likely associated with the galaxy closest to the quasar sight-line (galaxy "a"). Our morphokinematic analysis of galaxy "a" combined with the absorption line kinematics supports the interpretation that one of the absorption components originates from an extension of the stellar disk of galaxy "a", while the other component may arise in accreting gas in a warped disk with specific angular momentum $\sim 3$ times larger than the specific angular momentum of the galaxy halo. Such warped disks are common features in hydrodynamical simulations of cold-flow accretion onto galaxies; the data presented here provide observational evidence in favour of this scenario.
Recent X-ray observations of merger shocks in galaxy clusters have shown that the post-shock plasma is two-temperature, with the protons hotter than the electrons. By means of two-dimensional particle-in-cell simulations, we study the physics of electron irreversible heating in perpendicular low Mach number shocks, for a representative case with sonic Mach number of 3 and plasma beta of 16. We find that two basic ingredients are needed for electron entropy production: (i) an electron temperature anisotropy, induced by field amplification coupled to adiabatic invariance; and (ii) a mechanism to break the electron adiabatic invariance itself. In shocks, field amplification occurs at two major sites: at the shock ramp, where density compression leads to an increase of the frozen-in field; and farther downstream, where the shock-driven proton temperature anisotropy generates strong proton cyclotron and mirror modes. The electron temperature anisotropy induced by field amplification exceeds the threshold of the electron whistler instability. The growth of whistler waves breaks the electron adiabatic invariance, and allows for efficient entropy production. We find that the electron heating efficiency displays only a weak dependence on mass ratio (less than 30 percent drop, as we increase the mass ratio from 49 up to 1600). We develop an analytical model of electron irreversible heating and show that it is in excellent agreement with our simulation results.
The youngest low-mass protostars are known to be chemically rich, accreting matter most vigorously, and producing the most powerful outflows. Molecules are unique tracers of these phenomena. We use ALMA to study several outflow sources in the Serpens Main region. The most luminous source, Ser-SMM1, shows the richest chemical composition, but some complex molecules are also present in S68N. No emission from complex organics is detected toward Ser-emb 8N, which is the least luminous in the sample. We discuss whether these differences reflect an evolutionary effect or whether they are due to different physical structures. We also analyze the outflow structure from these young protostars by comparing emission of CO and SiO. EHV molecular jets originating from SMM1-a,b and Ser-emb 8N contrast with no such activity from S68N, which on the other hand presents a complex outflow structure.
Links to: arXiv, form interface, find, astro-ph, recent, 1710, contact, help (Access key information)
Using observed stellar mass functions out to $z=5$, we measure the main progenitor stellar mass growth of descendant galaxies with masses of $\log{M_{*}/M_{\odot}}=11.5,11.0,10.5,10.0$ at $z\sim0.1$ using an evolving cumulative number density selection. From these mass growth histories, we are able to measure the time at which half the total stellar mass of the descendant galaxy was assembled, $t_{a}$, which, in order of decreasing mass corresponds to redshifts of $z_{a}=1.28, 0.92, 0.60$ and $0.51$. We compare this to the median light-weighted stellar age $t_{*}$ ($z_{*} = 2.08, 1.49, 0.82$ and $0.37$) of a sample of low redshift SDSS galaxies (from the literature) and find the timescales are consistent with more massive galaxies forming a higher fraction of their stars ex-situ compared to lower mass descendants. We find that both $t_{*}$ and $t_{a}$ strongly correlate with mass which is in contrast to what is found in the EAGLE hydrodynamical simulation which shows a flat relationship between $t_{a}$ and $M_{*}$. However, the semi-analytic model of \citet{henriques2015} is consistent with the observations in both $t_{a}$ and $t_{*}$ with $M_{*}$, showing the most recent semi-analytic models are better able to decouple the evolution of the baryons from the dark matter in lower-mass galaxies.
We present new HST/WFC3 grism observations and re-analyse VLT data to unveil the continuum, variability and rest-frame UV lines of the three UV components of the most luminous Ly-alpha (Lya) emitter at z=6.6, COSMOS Redshift 7 (CR7). Our re-reduced, flux calibrated X-SHOOTER spectra of CR7 reveal a tentative detection of HeII with F(HeII)=$(1.8\pm0.7)\times10^{-17}$erg s$^{-1}$cm$^{-2}$ and we identify the signal (~2.6$\sigma$) as coming only from observations obtained along the major axis of Lya emission. There is a change of +0.2-0.5mag in UltraVISTA J band data for CR7 from DR2 to DR3, which virtually eliminates the strong J-band excess previously interpreted as being caused by HeII. Our WFC3 grism spectra provide a significant detection of the UV continuum of CR7's clump A, yielding an excellent fit to a power law with $\beta=-2.4\pm0.4$ and $M_{UV}=-21.7\pm0.3$, consistent with no variability. HST grism data fail to detect any rest-frame UV line in clump A above 3$\sigma$, yielding F(HeII)<$0.5\times10^{-17}$erg s$^{-1}$cm$^{-2}$ (EW$_0$<10A) at a 95% confidence level. Clump C is tentatively identified as a potential variable and high ionisation source with F(HeII)=$(1.0\pm0.4)\times10^{-17}$erg s$^{-1}$cm$^{-2}$. We perform CLOUDY modelling to constrain the metallicity and the ionising nature of CR7, and also make emission-line predictions for JWST/NIRSpec. CR7 seems to be actively forming stars without any clear AGN activity in clumps A and B, consistent with a metallicity of ~0.05-0.2 Z$_{\odot}$ and with component A experiencing the most massive starburst. Component C may host a high ionisation source/AGN. Our results highlight the need for spatially resolved information to study the complex formation and assembly of early galaxies within the epoch of re-ionisation.
We present a comprehensive study of the interstellar NaI $\lambda$5890, 5895 (NaD) resonant lines in a complete spectroscopic sample of $\sim 600,000$ passive, star-forming and starburst galaxies drawn from SDSS DR7 in order to look for cold-gas outflows in the local Universe. Individual galaxy spectra are stacked in bins of stellar mass and SFR and the dependence of galactic winds, with respect to the galaxies position in the SFR-$M_{\star}$ plane is investigated. While in most cases the interstellar medium (ISM) lines are fixed at the galaxy systemic velocity, at the higher SFR tail (SFR$>12.5 M_{\odot} yr^{-1}$), we find evidence of blue-shifted NaD absorption profiles, which we interpret as evidence of neutral outflowing gas. We explore the properties of the ISM in these galaxies with high SFR, in particular relating the absorption NaD line shape with the galaxy geometry in galaxies with different ionization mechanisms: AGN and star-formation. We find that: a) the ISM NaD absorption lines show a clear transition from a strong disk-like component, perfectly centered to the systemic velocity, in the edge-on system ($i > 50^\circ$ of the disk rotation axis), to an outflow, blue-shifted, component in face-on galaxies ($i < 50^\circ$); b) these trends are observed in galaxies classified as "purely" SF and AGN dominated objects. We compare the kinematics of the neutral gas with the kinematics of the ionized gas as traced by the [OIII]$\lambda$5007 emission lines. We find that, in these high SFR galaxies, the perturbations of the [OIII] emission line are present only in AGN or composite systems. In conclusion, we find that, in the local Universe, galactic winds show two faces which are related to two different ejection mechanisms, namely the neutral outflowing gas phase related to the star formation rate along the galaxy disk and the ionized winds related to the AGN feedback.
We measure the parsec-scale relationship between integrated CO intensity (I_CO) and visual extinction (A_V) in 24 local molecular clouds using maps of CO emission and dust optical depth from Planck. This relationship informs our understanding of CO emission across environments, but clean Milky Way measurements remain scarce. We find uniform I_CO for a given A_V, with the results bracketed by previous studies of the Pipe and Perseus clouds. Our measured I_CO-A_V relation broadly agrees with the standard Galactic CO-to-H2 conversion factor, the relation found for the Magellanic clouds at coarser resolution, and numerical simulations by Glover & Clark (2016). This supports the idea that CO emission primarily depends on shielding, which protects molecules from dissociating radiation. Evidence for CO saturation at high A_V and a threshold for CO emission at low A_V varies remains uncertain due to insufficient resolution and ambiguities in background subtraction. Resolution of order 0.1 pc may be required to measure these features. We use this I_CO-AV relation to predict how the CO-to-H2 conversion factor (X_CO) would change if the Solar Neighborhood clouds had different dust-to-gas ratio (metallicity). The calculations highlight the need for improved observations of the CO emission threshold and HI shielding layer depth. They are also sensitive to the shape of the column density distribution. Because local clouds collectively show a self-similar distribution, we predict a shallow metallicity dependence for X_CO down to a few tenths of solar metallicity. However, our calculations also imply dramatic variations in cloud-to-cloud X_CO at subsolar metallicity.
The life-cycle of cosmic dust grains is far from being understood and the origin and evolution of interstellar medium (ISM) grains is still under debate. In the ISM, the cosmic dust destruction rate is faster than the production rate by stellar sources. However, observations of ISM refractory matter suggest that to maintain a steady amount of cosmic grains, some supplementary production mechanism takes place. In this context, we aimed to study possible re-formation mechanisms of cosmic grains taking place at low temperature directly in the ISM. The low temperature condensation of carbonaceous materials has been investigated in experiments mimicking the ISM conditions. Gas-phase carbonaceous precursors created by laser ablation of graphite were forced to accrete on cold substrates (T about 10 K) representing surviving dust grains. The growing and evolution of the condensing carbonaceous precursors have been monitored by MIR and UV spectroscopy under a number of experimental scenarios. It is demonstrated, for the first time, the possibility to form ISM carbonaceous grains in "situ". The condensation process is governed by carbon chains that first condense into small carbon clusters and finally into more stable carbonaceous materials, which structural characteristics are comparable to the material formed in gas-phase condensation experiments at very high temperature. We also show that the so-formed fullerene-like carbonaceous material is transformed into a more ordered material under VUV processing. The cold condensation mechanisms here discussed can give fundamental clues to fully understand the balance between the timescale for dust injection, destruction and re-formation in the ISM.
We report the first high spectral resolution study of 17 M giants kinematically confirmed to lie within a few parsecs of the Galactic Center, using R=24,000 spectroscopy from Keck/NIRSPEC and a new linelist for the infrared K band. We consider their luminosities and kinematics, which classify these stars as members of the older stellar population and the central cluster. We find a median metallicity of <[Fe/H]>=-0.16 and a large spread from approximately -0.3 to +0.3 (quartiles). We find that the highest metallicities are [Fe/H]<+0.6, with most of the stars being at or below the Solar iron abundance. The abundances and the abundance distribution strongly resembles that of the Galactic bulge rather than disk or halo; in our small sample we find no statistical evidence for a dependence of velocity dispersion on metallicity.
The stellar velocity ellipsoid of the solar neighbour is re-examined using intermediate-old mono-abundances stellar groups with high quality chemistry data together with parallaxes and proper motions from Gaia DR1. We find the average velocity dispersion values for the three space velocity components for the thin and thick disc of (\sigma_{U},\sigma_{V},\sigma_{W})_{thin} = (33 \pm 4, 28 \pm 2, 23 \pm 2) and (\sigma_{U},\sigma_{V},\sigma_{W})_{thick} = (57 \pm 6, 38 \pm 5, 37 \pm 4) km s^{-1}, respectively. The mean values of the ratio between the semi-axes of the velocity ellipsoid for the thin disc are found to be, \sigma_{V}/\sigma_{U} = 0.70 \pm 0.13 and \sigma_{W}/\sigma_{U} is 0.64 \pm 0.08, while for the thick disc \sigma_{V}/\sigma_{U} = 0.67 \pm 0.11 and \sigma_{W}/\sigma_{U} is 0.66 \pm 0.11. Inputting these dispersions into the linear Str\"omberg relation for the thin disc groups, we find the Sun's velocity with respect to the LSR in Galactic rotation to be V_{\sun} = 13.9 \pm 3.4 km s^{-1}. A relation is found between the vertex deviation and the chemical abundances for the thin disc, ranging from -5 to +40^{\circ} as iron-abundance increases. For the thick disc we find a vertex deviation of l_{uv} \sim -15^{\circ}. The tilt angle (l_{uw}) in the U-W plane for the thin disc groups ranges from -10 to +15^\circ, but there is no evident relation between l_{uw} and the mean abundances. However we find a weak relation for l_{uw} as a function of iron abundances and \alpha-elements for most of the groups in the thick disc, where the tilt angle decreases from -5 to -20^\circ when [Fe/H] decreases and [\alpha/Fe] increases. The velocity anisotropy parameter is independent of the chemical group abundances and its value is nearly constant for both discs (\beta \sim 0.5), suggesting that the combined disc is dynamically relaxed.
We carried out an optical polarimetric study in the direction of the RCW95 star forming region in order to probe the sky-projected magnetic field structure by using the distribution of linear polarization segments which seem to be well aligned with the more extended cloud component. A mean polarization angle of $\theta=49.8^o\pm7.7^o$ was derived. Through the spectral dependence analysis of polarization it was possible to obtain the total-to-selective extinction ratio ($R_V$) by fitting the Serkowski function, resulting in a mean value of $R_V=2.93\pm0.47$. The foreground polarization component was estimated and is in agreement with previous studies in this direction of the Galaxy. Further, near-infrared images from Vista Variables in the Via L\'actea (VVV) survey were collected to improve the study of the stellar population associated with the HII region. The Automated Stellar Cluster Analysis (ASteCA) algorithm was employed to derive structural parameters for two clusters in the region, and a set of PAdova and TRieste Stellar Evolution Code (PARSEC) isochrones was superimposed on the decontaminated colour-magnitude diagrams (CMDs) to estimate an age of about 3 Myr for both clusters. Finally, from the near-infrared photometry study combined with spectra obtained with the Ohio State Infrared Imager and Spectrometer (OSIRIS) mounted at the Southern Astrophysics Research Telescope (SOAR) we derived the spectral classification of the main ionizing sources in the clusters associated with IRAS 15408$-$5356 and IRAS 15412$-$5359, both objects classified as O4 V stars.
The 5th RAVE data release is based on 520,781 spectra ($R\approx7500$ in the CaT region at $8410$ - $8795$\AA) of 457,588 unique stars. RAVE DR5 provides radial velocities, stellar parameters and individual abundances for up to seven elements and distances found using isochrones for a considerable subset of these objects. In particular, RAVE DR5 has 255,922 stellar observations that also have parallaxes and proper motions from the Tycho-Gaia astrometric solution (TGAS) in Gaia DR1. The combination of RAVE and TGAS thus provides the currently largest overlap of spectroscopic and space-based astrometric data and thus can serve as a formidable preview of what Gaia is going to deliver in coming data releases. Basic properties of the RAVE+TGAS survey and its derived data products are presented as well as first applications w.r.t wave-like patterns in the disk structure. An outlook to the 6th RAVE data release is given.
We present the discovery of strong Balmer line absorption in H$\alpha$ to H$\gamma$ in two luminous low-ionization broad absorption line quasars (LoBAL QSOs) at z~1.5, with black hole masses around $10^{10}$ $M_\odot$ from near-IR spectroscopy. There are only two previously known quasars at z>1 showing Balmer line absorption. SDSS J1019+0225 shows blueshifted absorption by ~1400 km/s with an H$\alpha$ rest-frame equivalent width of 13 \AA{}. In SDSS J0859+4239 we find redshifted absorption by ~500 km/s with an H$\alpha$ rest-frame equivalent width of 7 \AA{}. The redshifted absorption could indicate an inflow of high density gas onto the black hole, while we cannot rule out alternative interpretations. The Balmer line absorption in both objects appears to be saturated, indicating partial coverage of the background source by the absorber. We estimate the covering fractions and optical depth of the absorber and derive neutral hydrogen column densities, $N_{\rm{HI}}\sim1.3\times 10^{18}$ cm$^{-2}$ for SDSS J1019+0225 and $N_{\rm{HI}}\sim9\times 10^{17}$ cm$^{-2}$ for SDSS J0859+4239, respectively. In addition, the optical spectra reveal also absorption troughs in HeI* $\lambda3889$ and $\lambda3189$ in both objects.
We report an observational study of the galactic giant molecular cloud (GMC) associated with the three HII regions including the infrared ring N35 and two nearby HII regions G024.392+00.072 (HII region A) and G024.510-00.060 (HII region B). As a part of the FOREST Unbiased Galactic Plane Imaging survey with the Nobeyama 45-m telescope (FUGIN) project, we obtained the CO J=1-0 data covering the entirety of the GMC at a spatial resolution of 21arcsec. Our CO data revealed that the GMC having a total molecular mass of 1.3x10^6Mo, has two velocity components over 10km/s. The majority of the molecular gas in the GMC is distributed in the lower velocity component (LVC) at 110-114km/s, while the higher velocity components (HVCs) around 118-126km/s comprises three molecular clouds distributed around the HII regions. LVC and HVCs are separated by the CO emission with intermediate intensities in the velocity space, and comparisons of spatial distribution between LVC and HVCs show complementary distributions of gas. Based on the analytic model of HII region expansion, we discussed that the observed gas distribution cannot be explained with expansion of the HII regions. Alternatively, we proposed a cloud-cloud collision scenario to interpret the observed properties of GMC. We discussed that collisions between LVC and HVCs happened 0.2-1.0Myr ago at relative velocities of ~10km/s and triggered the high-mass star formation which ionizes the three HII regions in the present GMC. The intermediate velocity features between LVC and HVCs can be understood as broad bridge features, which indicate turbulent motion of gas at the interfaces of the collisions, and the spatially complementary distributions between LVC and HVCs are understood as the cavities created on LVC by the collisions.
To investigate the connection between radio activity and AGN outflows, we present a study of ionized gas kinematics based on [O III] $\lambda$5007 emission line along the large-scale radio jet for six radio AGNs. These AGNs are selected based on the radio activity (i.e., $\mathrm{L_{1.4GHz}}$ $\geqslant$ 10$^{39.8}$ erg s$^{-1}$) as well as optical emission line properties as type 2 AGNs. Using the Red Channel Cross Dispersed Echellette Spectrograph at the Multiple Mirror Telescope, we investigate in detail the [O III] and stellar kinematics. We spatially resolve and probe the central AGN-photoionization sizes, which is important in understanding the structures and evolutions of galaxies. We find that the typical central AGN-photoionization radius of our targets are in range of 0.9$-$1.6 kpc, consistent with the size-luminosity relation of [O III] in the previous studies. We investigate the [O III] kinematics along the large-scale radio jets to test whether there is a link between gas outflows in the narrow-line region and extended radio jet emissions. Contrary to our expectation, we find no evidence that the gas outflows are directly connected to the large scale radio jets.
We study the role of cold gas in quenching star formation in the green valley by analysing ALMA $^{12}$CO (1-0) observations of three galaxies with resolved optical spectroscopy from the MaNGA survey. We present resolution-matched maps of the star formation rate and molecular gas mass. These data are used to calculate the star formation efficiency (SFE) and gas fraction ($f_{\rm~gas}$) for these galaxies separately in the central `bulge' regions and outer disks. We find that, for the two galaxies whose global specific star formation rate (sSFR) deviates most from the star formation main sequence, the gas fraction in the bulges is significantly lower than that in their disks, supporting an `inside-out' model of galaxy quenching. For the two galaxies where SFE can be reliably determined in the central regions, the bulges and disks share similar SFEs. This suggests that a decline in $f_{\rm~gas}$ is the main driver of lowered sSFR in bulges compared to disks in green valley galaxies. Within the disks, there exist common correlations between the sSFR and SFE and between sSFR and $f_{\rm~gas}$ on kpc scales -- the local SFE or $f_{\rm~gas}$ in the disks declines with local sSFR. Our results support a picture in which the sSFR in bulges is primarily controlled by $f_{\rm~gas}$, whereas both SFE and $f_{\rm~gas}$ play a role in lowering the sSFR in disks. A larger sample is required to confirm if the trend established in this work is representative of green valley as a whole.
Galaxies are surrounded by extended gaseous halos which store significant fractions of chemical elements. These are syntethized by the stellar populations and later ejected into the circumgalactic medium (CGM) by different mechanism, of which supernova feedback is considered one of the most relevant. We explore the properties of this metal reservoir surrounding star-forming galaxies in a cosmological context aiming to investigate the chemical loop between galaxies and their CGM, and the ability of the subgrid models to reproduce observational results. Using cosmological hydrodynamical simulations, we analyse the gas-phase chemical contents of galaxies with stellar masses in the range $10^{9} - 10^{11}\,{\rm M}_{\odot}$. We estimate the fractions of metals stored in the different CGM phases, and the predicted OVI and SiIII column densities within the virial radius. We find roughly $10^{7}\,{\rm M}_{\odot}$ of oxygen in the CGM of simulated galaxies having $M_{\star}{\sim}10^{10}\,{\rm M}_{\odot}$, in fair agreement with the lower limits imposed by observations. The $M_{\rm oxy}$ is found to correlate with $M_{\star}$, at odds with current observational trends but in agreement with other numerical results. The estimated profiles of OVI column density reveal a substantial shortage of that ion, whereas SiIII, which probes the cool phase, is overpredicted. The analysis of the relative contributions of both ions from the hot, warm and cool phases suggests that the warm gas ($ 10^5~{\rm K} < T < 10^6~{\rm K}$) should be more abundant in order to bridge the mismatch with the observations, or alternatively, that more metals should be stored in this gas-phase. Adittionally, we find that the X-ray coronae around the simulated galaxies have luminosities and temperatures in decent agreement with the available observational estimates. [abridged]
Understanding of the role of X-rays for driving the thermal evolution of the intergalactic medium (IGM) at high redshifts is one of important questions in astrophysics. High-mass X-ray binaries (HMXBs) in early stellar populations are prime X-ray source; however, their formation efficiency is not well understood. Using $N$-body simulations, we estimate the HMXB formation rate via mutual gravitational interactions of nascent, small groups of the Population~III stars. We find that HMXBs form at a rate of one per $\gtrsim 10^{4}M_{\odot}$ in newly born stars, and that they emit with a power of $\sim 10^{41} {\rm erg}~{\rm s}^{-1}$ in the $2-10$ keV band per star formation rate (SFR). This value is a factor $\sim 10^{2}$ larger than what is observed in star forming galaxies at lower redshifts; the X-ray production from early HMXBs would have been even more copious, if they also formed \textit{in situ} or via migration in protostellar disks. Combining our results with earlier studies suggests that early HMXBs were highly effective at heating the IGM and leaving a strong 21 cm signature. We discuss broader implications of our results, such as the rate of long gamma-ray bursts from Population~III stars and the direct collapse channel for massive black hole formation.
Stars and their corresponding protoplanetary disks form in diverse environments. To account for these natural variations, we investigate the formation process around nine solar mass stars with a maximum resolution of 2 AU in a Giant Molecular Cloud of (40 pc)$^3$ in volume by using the adaptive mesh refinement code \ramses. The magnetohydrodynamic simulations reveal that the accretion process is heterogeneous in time, in space, and among protostars of otherwise similar mass. During the first roughly 100 kyr of a protostar evolving to about a solar mass, the accretion rates peak around $10^{-5}$ to $10^{-4}$ M$_{\odot}$ yr$^{-1}$ shortly after its birth, declining with time after that. The different environments also affect the spatial accretion, and infall of material to the star-disk system is mostly through filaments and sheets. Furthermore, the formation and evolution of disks varies significantly from star to star. We interpret the variety in disk formation as a consequence of the differences in the combined effects of magnetic fields and turbulence that may cause differences in the efficiency of magnetic braking, as well as differences in the strength and distribution of specific angular momentum.
The thermodynamic structure of hot gas in galaxy clusters is sensitive to astrophysical processes and typically difficult to model with galaxy formation simulations. We explore the fraction of cool-core (CC) clusters in a large sample of $370$ clusters from IllustrisTNG, examining six common CC definitions. IllustrisTNG produces continuous CC criteria distributions, the extremes of which are classified as CC and non-cool-core (NCC), and the criteria are increasingly correlated for more massive clusters. At $z=0$ the CC fraction is systematically lower than observed for the complete sample, selecting massive systems increases the CC fraction for $3$ criteria and reduces it for others. This result is partly driven by systematic differences between the simulated and observed gas fraction profiles. The simulated CC fraction increases more rapidly with redshift than observed, independent of mass or redshift range, and the CC fraction is overpredicted at $z\geq1$. The conversion of CCs to NCCs begins later and acts more rapidly in the simulations. Examining the fraction of CCs and NCCs defined as relaxed we find no evidence that CCs are more relaxed, suggesting that mergers are not solely responsible for disrupting CCs. A comparison of the median thermodynamic profiles defined by different CC criteria shows that the extent to which they evolve in the cluster core is dependent on the CC criteria. We conclude that the thermodynamic structure of galaxy clusters in IllustrisTNG shares many similarities with observations, but achieving better agreement most likely requires modifications of the underlying galaxy formation model.
Far-infrared and (sub)millimeter fluxes can be used to study dust in protoplanetary disks, the building blocks of planets. Here, we combine observations from the Herschel Space Observatory with ancillary data of 284 protoplanetary disks in the Taurus, Chamaeleon I, and Ophiuchus star-forming regions, covering from the optical to mm/cm wavelengths. We analyze their spectral indices as a function of wavelength and determine their (sub)millimeter slopes when possible. Most disks display observational evidence of grain growth, in agreement with previous studies. No correlation is found between other tracers of disk evolution and the millimeter spectral indices. A simple disk model is used to fit these sources, and we derive posterior distributions for the optical depth at 1.3 mm and 10 au, the disk temperature at this same radius, and the dust opacity spectral index. We find the fluxes at 70 microns to correlate strongly with disk temperatures at 10 au, as derived from these simple models. We find tentative evidence for spectral indices in Chamaeleon I being steeper than those of disks in Taurus/Ophiuchus, although more millimeter observations are needed to confirm this trend and identify its possible origin. Additionally, we determine the median spectral energy distribution of each region and find them to be similar across the entire wavelength range studied, possibly due to the large scatter in disk properties and morphologies.
We have studied the chemical enrichment history of the interstellar medium through an analysis of the n-dimensional stellar abundances space. This work is a non-parametric analysis of the stellar chemical abundance space. The main goal is to study the stars from their organization within this abundance space. Within this space, we seek to find clusters (in a statistical sense), that is, stars likely to share similar chemo-evolutionary history, using two methods: the hierarchical clustering and the principal component analysis. We analysed some selected abundance surveys available in the literature. For each sample, we labelled the group of stars according to its average abundance curve. In all samples, we identify the existence of a main enrichment pattern of the stars, which we call chemical enrichment flow. This flow is set by the structured and well defined mean rate at which the abundances of the interstellar medium increase, resulting from the mixture of the material ejected from the stars and stellar mass loss and interstellar medium gas. One of the main results of our analysis is the identification of subgroups of stars with peculiar chemistry. These stars are situated in regions outside of the enrichment flow in the abundance space. These peculiar stars show a mismatch in the enrichment rate of a few elements, such as Mg, Si, Sc and V, when compared to the mean enrichment rate of the other elements of the same stars. We believe the existence of these groups of stars with peculiar chemistry may be related to the accretion of planetary material onto stellar surfaces or may be due to production of the same chemical element by different nucleosynthetic sites.
We have developed a new X-ray absorption model, called {\tt IONeq}, which computes the optical depth $\tau(E)$ simultaneously for ions of all abundant elements, assuming ionization equilibrium and taking into account turbulent broadening. We use this model to analyze the interstellar medium (ISM) absorption features in the Milky Way for a sample of 18 galactic (LMXBs) and 42 extragalactic sources (mainly Blazars). The absorbing ISM was modeled as a combination of three components/phases - neutral ($T\lesssim1\times10^{4}$ K), warm ($T\sim5\times 10^{4}$ K) and hot ($T\sim2\times10^{6}$ K). We found that the spatial distribution of both, neutral and warm components, are difficult to describe using smooth profiles due to nonuniform distribution of the column densities over the sky. For the hot phase we used a combination of a flattened disk and a halo, finding comparable column densities for both spatial components, in the order of $\sim 6-7\times10^{18}\;{\rm cm^{-2}}$, although this conclusion depends on the adopted parametrization. If the halo component has sub-solar abundance $Z$, then the column density has be scaled up by a factor $\frac{Z_\odot}{Z}$. The vertically integrated column densities of the disk components suggests the following mass fractions for these three ISM phases in the Galactic disk: neutral $\sim~89\%$, warm $\sim 8\%$ and hot $\sim 3\%$ components, respectively. The constraints on the radial distribution of the halo component of the hot component are weak.
Complex organic molecules (COMs) have been observed in comets, hot cores and cold dense regions of the interstellar medium. It is generally accepted that these COMs form on icy dust grain through the recombination reaction of radicals triggered by either energetic UV- photon or non-energetic H-atom addition processing. In this work, we present for the first time laboratory studies that allow for quantitative comparison of hydrogenation and UV-induced reactions as well as their cumulative effect in astronomically relevant CO:CH3OH=4:1 ice analogues. The formation of glycolaldehyde (GA) and ethylene glycol (EG) is confirmed in pure hydrogenation experiments at 14 K, except methyl formate (MF), which is only clearly observed in photolysis. The fractions for MF:GA:EG are 0 : (0.2-0.4) : (0.8-0.6) for pure hydrogenation, and 0.2 : 0.3 : 0.5 for UV involving experiments and can offer a diagnostic tool to derive the chemical origin of these species. The GA/EG ratios in the laboratory (0.3-1.5) compare well with observations toward different objects.
Links to: arXiv, form interface, find, astro-ph, recent, 1710, contact, help (Access key information)
Using observed stellar mass functions out to $z=5$, we measure the main progenitor stellar mass growth of descendant galaxies with masses of $\log{M_{*}/M_{\odot}}=11.5,11.0,10.5,10.0$ at $z\sim0.1$ using an evolving cumulative number density selection. From these mass growth histories, we are able to measure the time at which half the total stellar mass of the descendant galaxy was assembled, $t_{a}$, which, in order of decreasing mass corresponds to redshifts of $z_{a}=1.28, 0.92, 0.60$ and $0.51$. We compare this to the median light-weighted stellar age $t_{*}$ ($z_{*} = 2.08, 1.49, 0.82$ and $0.37$) of a sample of low redshift SDSS galaxies (from the literature) and find the timescales are consistent with more massive galaxies forming a higher fraction of their stars ex-situ compared to lower mass descendants. We find that both $t_{*}$ and $t_{a}$ strongly correlate with mass which is in contrast to what is found in the EAGLE hydrodynamical simulation which shows a flat relationship between $t_{a}$ and $M_{*}$. However, the semi-analytic model of \citet{henriques2015} is consistent with the observations in both $t_{a}$ and $t_{*}$ with $M_{*}$, showing the most recent semi-analytic models are better able to decouple the evolution of the baryons from the dark matter in lower-mass galaxies.
We present new HST/WFC3 grism observations and re-analyse VLT data to unveil the continuum, variability and rest-frame UV lines of the three UV components of the most luminous Ly-alpha (Lya) emitter at z=6.6, COSMOS Redshift 7 (CR7). Our re-reduced, flux calibrated X-SHOOTER spectra of CR7 reveal a tentative detection of HeII with F(HeII)=$(1.8\pm0.7)\times10^{-17}$erg s$^{-1}$cm$^{-2}$ and we identify the signal (~2.6$\sigma$) as coming only from observations obtained along the major axis of Lya emission. There is a change of +0.2-0.5mag in UltraVISTA J band data for CR7 from DR2 to DR3, which virtually eliminates the strong J-band excess previously interpreted as being caused by HeII. Our WFC3 grism spectra provide a significant detection of the UV continuum of CR7's clump A, yielding an excellent fit to a power law with $\beta=-2.4\pm0.4$ and $M_{UV}=-21.7\pm0.3$, consistent with no variability. HST grism data fail to detect any rest-frame UV line in clump A above 3$\sigma$, yielding F(HeII)<$0.5\times10^{-17}$erg s$^{-1}$cm$^{-2}$ (EW$_0$<10A) at a 95% confidence level. Clump C is tentatively identified as a potential variable and high ionisation source with F(HeII)=$(1.0\pm0.4)\times10^{-17}$erg s$^{-1}$cm$^{-2}$. We perform CLOUDY modelling to constrain the metallicity and the ionising nature of CR7, and also make emission-line predictions for JWST/NIRSpec. CR7 seems to be actively forming stars without any clear AGN activity in clumps A and B, consistent with a metallicity of ~0.05-0.2 Z$_{\odot}$ and with component A experiencing the most massive starburst. Component C may host a high ionisation source/AGN. Our results highlight the need for spatially resolved information to study the complex formation and assembly of early galaxies within the epoch of re-ionisation.
We present a comprehensive study of the interstellar NaI $\lambda$5890, 5895 (NaD) resonant lines in a complete spectroscopic sample of $\sim 600,000$ passive, star-forming and starburst galaxies drawn from SDSS DR7 in order to look for cold-gas outflows in the local Universe. Individual galaxy spectra are stacked in bins of stellar mass and SFR and the dependence of galactic winds, with respect to the galaxies position in the SFR-$M_{\star}$ plane is investigated. While in most cases the interstellar medium (ISM) lines are fixed at the galaxy systemic velocity, at the higher SFR tail (SFR$>12.5 M_{\odot} yr^{-1}$), we find evidence of blue-shifted NaD absorption profiles, which we interpret as evidence of neutral outflowing gas. We explore the properties of the ISM in these galaxies with high SFR, in particular relating the absorption NaD line shape with the galaxy geometry in galaxies with different ionization mechanisms: AGN and star-formation. We find that: a) the ISM NaD absorption lines show a clear transition from a strong disk-like component, perfectly centered to the systemic velocity, in the edge-on system ($i > 50^\circ$ of the disk rotation axis), to an outflow, blue-shifted, component in face-on galaxies ($i < 50^\circ$); b) these trends are observed in galaxies classified as "purely" SF and AGN dominated objects. We compare the kinematics of the neutral gas with the kinematics of the ionized gas as traced by the [OIII]$\lambda$5007 emission lines. We find that, in these high SFR galaxies, the perturbations of the [OIII] emission line are present only in AGN or composite systems. In conclusion, we find that, in the local Universe, galactic winds show two faces which are related to two different ejection mechanisms, namely the neutral outflowing gas phase related to the star formation rate along the galaxy disk and the ionized winds related to the AGN feedback.
We measure the parsec-scale relationship between integrated CO intensity (I_CO) and visual extinction (A_V) in 24 local molecular clouds using maps of CO emission and dust optical depth from Planck. This relationship informs our understanding of CO emission across environments, but clean Milky Way measurements remain scarce. We find uniform I_CO for a given A_V, with the results bracketed by previous studies of the Pipe and Perseus clouds. Our measured I_CO-A_V relation broadly agrees with the standard Galactic CO-to-H2 conversion factor, the relation found for the Magellanic clouds at coarser resolution, and numerical simulations by Glover & Clark (2016). This supports the idea that CO emission primarily depends on shielding, which protects molecules from dissociating radiation. Evidence for CO saturation at high A_V and a threshold for CO emission at low A_V varies remains uncertain due to insufficient resolution and ambiguities in background subtraction. Resolution of order 0.1 pc may be required to measure these features. We use this I_CO-AV relation to predict how the CO-to-H2 conversion factor (X_CO) would change if the Solar Neighborhood clouds had different dust-to-gas ratio (metallicity). The calculations highlight the need for improved observations of the CO emission threshold and HI shielding layer depth. They are also sensitive to the shape of the column density distribution. Because local clouds collectively show a self-similar distribution, we predict a shallow metallicity dependence for X_CO down to a few tenths of solar metallicity. However, our calculations also imply dramatic variations in cloud-to-cloud X_CO at subsolar metallicity.
The life-cycle of cosmic dust grains is far from being understood and the origin and evolution of interstellar medium (ISM) grains is still under debate. In the ISM, the cosmic dust destruction rate is faster than the production rate by stellar sources. However, observations of ISM refractory matter suggest that to maintain a steady amount of cosmic grains, some supplementary production mechanism takes place. In this context, we aimed to study possible re-formation mechanisms of cosmic grains taking place at low temperature directly in the ISM. The low temperature condensation of carbonaceous materials has been investigated in experiments mimicking the ISM conditions. Gas-phase carbonaceous precursors created by laser ablation of graphite were forced to accrete on cold substrates (T about 10 K) representing surviving dust grains. The growing and evolution of the condensing carbonaceous precursors have been monitored by MIR and UV spectroscopy under a number of experimental scenarios. It is demonstrated, for the first time, the possibility to form ISM carbonaceous grains in "situ". The condensation process is governed by carbon chains that first condense into small carbon clusters and finally into more stable carbonaceous materials, which structural characteristics are comparable to the material formed in gas-phase condensation experiments at very high temperature. We also show that the so-formed fullerene-like carbonaceous material is transformed into a more ordered material under VUV processing. The cold condensation mechanisms here discussed can give fundamental clues to fully understand the balance between the timescale for dust injection, destruction and re-formation in the ISM.
We report the first high spectral resolution study of 17 M giants kinematically confirmed to lie within a few parsecs of the Galactic Center, using R=24,000 spectroscopy from Keck/NIRSPEC and a new linelist for the infrared K band. We consider their luminosities and kinematics, which classify these stars as members of the older stellar population and the central cluster. We find a median metallicity of <[Fe/H]>=-0.16 and a large spread from approximately -0.3 to +0.3 (quartiles). We find that the highest metallicities are [Fe/H]<+0.6, with most of the stars being at or below the Solar iron abundance. The abundances and the abundance distribution strongly resembles that of the Galactic bulge rather than disk or halo; in our small sample we find no statistical evidence for a dependence of velocity dispersion on metallicity.
The stellar velocity ellipsoid of the solar neighbour is re-examined using intermediate-old mono-abundances stellar groups with high quality chemistry data together with parallaxes and proper motions from Gaia DR1. We find the average velocity dispersion values for the three space velocity components for the thin and thick disc of (\sigma_{U},\sigma_{V},\sigma_{W})_{thin} = (33 \pm 4, 28 \pm 2, 23 \pm 2) and (\sigma_{U},\sigma_{V},\sigma_{W})_{thick} = (57 \pm 6, 38 \pm 5, 37 \pm 4) km s^{-1}, respectively. The mean values of the ratio between the semi-axes of the velocity ellipsoid for the thin disc are found to be, \sigma_{V}/\sigma_{U} = 0.70 \pm 0.13 and \sigma_{W}/\sigma_{U} is 0.64 \pm 0.08, while for the thick disc \sigma_{V}/\sigma_{U} = 0.67 \pm 0.11 and \sigma_{W}/\sigma_{U} is 0.66 \pm 0.11. Inputting these dispersions into the linear Str\"omberg relation for the thin disc groups, we find the Sun's velocity with respect to the LSR in Galactic rotation to be V_{\sun} = 13.9 \pm 3.4 km s^{-1}. A relation is found between the vertex deviation and the chemical abundances for the thin disc, ranging from -5 to +40^{\circ} as iron-abundance increases. For the thick disc we find a vertex deviation of l_{uv} \sim -15^{\circ}. The tilt angle (l_{uw}) in the U-W plane for the thin disc groups ranges from -10 to +15^\circ, but there is no evident relation between l_{uw} and the mean abundances. However we find a weak relation for l_{uw} as a function of iron abundances and \alpha-elements for most of the groups in the thick disc, where the tilt angle decreases from -5 to -20^\circ when [Fe/H] decreases and [\alpha/Fe] increases. The velocity anisotropy parameter is independent of the chemical group abundances and its value is nearly constant for both discs (\beta \sim 0.5), suggesting that the combined disc is dynamically relaxed.
We carried out an optical polarimetric study in the direction of the RCW95 star forming region in order to probe the sky-projected magnetic field structure by using the distribution of linear polarization segments which seem to be well aligned with the more extended cloud component. A mean polarization angle of $\theta=49.8^o\pm7.7^o$ was derived. Through the spectral dependence analysis of polarization it was possible to obtain the total-to-selective extinction ratio ($R_V$) by fitting the Serkowski function, resulting in a mean value of $R_V=2.93\pm0.47$. The foreground polarization component was estimated and is in agreement with previous studies in this direction of the Galaxy. Further, near-infrared images from Vista Variables in the Via L\'actea (VVV) survey were collected to improve the study of the stellar population associated with the HII region. The Automated Stellar Cluster Analysis (ASteCA) algorithm was employed to derive structural parameters for two clusters in the region, and a set of PAdova and TRieste Stellar Evolution Code (PARSEC) isochrones was superimposed on the decontaminated colour-magnitude diagrams (CMDs) to estimate an age of about 3 Myr for both clusters. Finally, from the near-infrared photometry study combined with spectra obtained with the Ohio State Infrared Imager and Spectrometer (OSIRIS) mounted at the Southern Astrophysics Research Telescope (SOAR) we derived the spectral classification of the main ionizing sources in the clusters associated with IRAS 15408$-$5356 and IRAS 15412$-$5359, both objects classified as O4 V stars.
The 5th RAVE data release is based on 520,781 spectra ($R\approx7500$ in the CaT region at $8410$ - $8795$\AA) of 457,588 unique stars. RAVE DR5 provides radial velocities, stellar parameters and individual abundances for up to seven elements and distances found using isochrones for a considerable subset of these objects. In particular, RAVE DR5 has 255,922 stellar observations that also have parallaxes and proper motions from the Tycho-Gaia astrometric solution (TGAS) in Gaia DR1. The combination of RAVE and TGAS thus provides the currently largest overlap of spectroscopic and space-based astrometric data and thus can serve as a formidable preview of what Gaia is going to deliver in coming data releases. Basic properties of the RAVE+TGAS survey and its derived data products are presented as well as first applications w.r.t wave-like patterns in the disk structure. An outlook to the 6th RAVE data release is given.
We present the discovery of strong Balmer line absorption in H$\alpha$ to H$\gamma$ in two luminous low-ionization broad absorption line quasars (LoBAL QSOs) at z~1.5, with black hole masses around $10^{10}$ $M_\odot$ from near-IR spectroscopy. There are only two previously known quasars at z>1 showing Balmer line absorption. SDSS J1019+0225 shows blueshifted absorption by ~1400 km/s with an H$\alpha$ rest-frame equivalent width of 13 \AA{}. In SDSS J0859+4239 we find redshifted absorption by ~500 km/s with an H$\alpha$ rest-frame equivalent width of 7 \AA{}. The redshifted absorption could indicate an inflow of high density gas onto the black hole, while we cannot rule out alternative interpretations. The Balmer line absorption in both objects appears to be saturated, indicating partial coverage of the background source by the absorber. We estimate the covering fractions and optical depth of the absorber and derive neutral hydrogen column densities, $N_{\rm{HI}}\sim1.3\times 10^{18}$ cm$^{-2}$ for SDSS J1019+0225 and $N_{\rm{HI}}\sim9\times 10^{17}$ cm$^{-2}$ for SDSS J0859+4239, respectively. In addition, the optical spectra reveal also absorption troughs in HeI* $\lambda3889$ and $\lambda3189$ in both objects.
We report an observational study of the galactic giant molecular cloud (GMC) associated with the three HII regions including the infrared ring N35 and two nearby HII regions G024.392+00.072 (HII region A) and G024.510-00.060 (HII region B). As a part of the FOREST Unbiased Galactic Plane Imaging survey with the Nobeyama 45-m telescope (FUGIN) project, we obtained the CO J=1-0 data covering the entirety of the GMC at a spatial resolution of 21arcsec. Our CO data revealed that the GMC having a total molecular mass of 1.3x10^6Mo, has two velocity components over 10km/s. The majority of the molecular gas in the GMC is distributed in the lower velocity component (LVC) at 110-114km/s, while the higher velocity components (HVCs) around 118-126km/s comprises three molecular clouds distributed around the HII regions. LVC and HVCs are separated by the CO emission with intermediate intensities in the velocity space, and comparisons of spatial distribution between LVC and HVCs show complementary distributions of gas. Based on the analytic model of HII region expansion, we discussed that the observed gas distribution cannot be explained with expansion of the HII regions. Alternatively, we proposed a cloud-cloud collision scenario to interpret the observed properties of GMC. We discussed that collisions between LVC and HVCs happened 0.2-1.0Myr ago at relative velocities of ~10km/s and triggered the high-mass star formation which ionizes the three HII regions in the present GMC. The intermediate velocity features between LVC and HVCs can be understood as broad bridge features, which indicate turbulent motion of gas at the interfaces of the collisions, and the spatially complementary distributions between LVC and HVCs are understood as the cavities created on LVC by the collisions.
To investigate the connection between radio activity and AGN outflows, we present a study of ionized gas kinematics based on [O III] $\lambda$5007 emission line along the large-scale radio jet for six radio AGNs. These AGNs are selected based on the radio activity (i.e., $\mathrm{L_{1.4GHz}}$ $\geqslant$ 10$^{39.8}$ erg s$^{-1}$) as well as optical emission line properties as type 2 AGNs. Using the Red Channel Cross Dispersed Echellette Spectrograph at the Multiple Mirror Telescope, we investigate in detail the [O III] and stellar kinematics. We spatially resolve and probe the central AGN-photoionization sizes, which is important in understanding the structures and evolutions of galaxies. We find that the typical central AGN-photoionization radius of our targets are in range of 0.9$-$1.6 kpc, consistent with the size-luminosity relation of [O III] in the previous studies. We investigate the [O III] kinematics along the large-scale radio jets to test whether there is a link between gas outflows in the narrow-line region and extended radio jet emissions. Contrary to our expectation, we find no evidence that the gas outflows are directly connected to the large scale radio jets.
We study the role of cold gas in quenching star formation in the green valley by analysing ALMA $^{12}$CO (1-0) observations of three galaxies with resolved optical spectroscopy from the MaNGA survey. We present resolution-matched maps of the star formation rate and molecular gas mass. These data are used to calculate the star formation efficiency (SFE) and gas fraction ($f_{\rm~gas}$) for these galaxies separately in the central `bulge' regions and outer disks. We find that, for the two galaxies whose global specific star formation rate (sSFR) deviates most from the star formation main sequence, the gas fraction in the bulges is significantly lower than that in their disks, supporting an `inside-out' model of galaxy quenching. For the two galaxies where SFE can be reliably determined in the central regions, the bulges and disks share similar SFEs. This suggests that a decline in $f_{\rm~gas}$ is the main driver of lowered sSFR in bulges compared to disks in green valley galaxies. Within the disks, there exist common correlations between the sSFR and SFE and between sSFR and $f_{\rm~gas}$ on kpc scales -- the local SFE or $f_{\rm~gas}$ in the disks declines with local sSFR. Our results support a picture in which the sSFR in bulges is primarily controlled by $f_{\rm~gas}$, whereas both SFE and $f_{\rm~gas}$ play a role in lowering the sSFR in disks. A larger sample is required to confirm if the trend established in this work is representative of green valley as a whole.
Galaxies are surrounded by extended gaseous halos which store significant fractions of chemical elements. These are syntethized by the stellar populations and later ejected into the circumgalactic medium (CGM) by different mechanism, of which supernova feedback is considered one of the most relevant. We explore the properties of this metal reservoir surrounding star-forming galaxies in a cosmological context aiming to investigate the chemical loop between galaxies and their CGM, and the ability of the subgrid models to reproduce observational results. Using cosmological hydrodynamical simulations, we analyse the gas-phase chemical contents of galaxies with stellar masses in the range $10^{9} - 10^{11}\,{\rm M}_{\odot}$. We estimate the fractions of metals stored in the different CGM phases, and the predicted OVI and SiIII column densities within the virial radius. We find roughly $10^{7}\,{\rm M}_{\odot}$ of oxygen in the CGM of simulated galaxies having $M_{\star}{\sim}10^{10}\,{\rm M}_{\odot}$, in fair agreement with the lower limits imposed by observations. The $M_{\rm oxy}$ is found to correlate with $M_{\star}$, at odds with current observational trends but in agreement with other numerical results. The estimated profiles of OVI column density reveal a substantial shortage of that ion, whereas SiIII, which probes the cool phase, is overpredicted. The analysis of the relative contributions of both ions from the hot, warm and cool phases suggests that the warm gas ($ 10^5~{\rm K} < T < 10^6~{\rm K}$) should be more abundant in order to bridge the mismatch with the observations, or alternatively, that more metals should be stored in this gas-phase. Adittionally, we find that the X-ray coronae around the simulated galaxies have luminosities and temperatures in decent agreement with the available observational estimates. [abridged]
Understanding of the role of X-rays for driving the thermal evolution of the intergalactic medium (IGM) at high redshifts is one of important questions in astrophysics. High-mass X-ray binaries (HMXBs) in early stellar populations are prime X-ray source; however, their formation efficiency is not well understood. Using $N$-body simulations, we estimate the HMXB formation rate via mutual gravitational interactions of nascent, small groups of the Population~III stars. We find that HMXBs form at a rate of one per $\gtrsim 10^{4}M_{\odot}$ in newly born stars, and that they emit with a power of $\sim 10^{41} {\rm erg}~{\rm s}^{-1}$ in the $2-10$ keV band per star formation rate (SFR). This value is a factor $\sim 10^{2}$ larger than what is observed in star forming galaxies at lower redshifts; the X-ray production from early HMXBs would have been even more copious, if they also formed \textit{in situ} or via migration in protostellar disks. Combining our results with earlier studies suggests that early HMXBs were highly effective at heating the IGM and leaving a strong 21 cm signature. We discuss broader implications of our results, such as the rate of long gamma-ray bursts from Population~III stars and the direct collapse channel for massive black hole formation.
Stars and their corresponding protoplanetary disks form in diverse environments. To account for these natural variations, we investigate the formation process around nine solar mass stars with a maximum resolution of 2 AU in a Giant Molecular Cloud of (40 pc)$^3$ in volume by using the adaptive mesh refinement code \ramses. The magnetohydrodynamic simulations reveal that the accretion process is heterogeneous in time, in space, and among protostars of otherwise similar mass. During the first roughly 100 kyr of a protostar evolving to about a solar mass, the accretion rates peak around $10^{-5}$ to $10^{-4}$ M$_{\odot}$ yr$^{-1}$ shortly after its birth, declining with time after that. The different environments also affect the spatial accretion, and infall of material to the star-disk system is mostly through filaments and sheets. Furthermore, the formation and evolution of disks varies significantly from star to star. We interpret the variety in disk formation as a consequence of the differences in the combined effects of magnetic fields and turbulence that may cause differences in the efficiency of magnetic braking, as well as differences in the strength and distribution of specific angular momentum.
The thermodynamic structure of hot gas in galaxy clusters is sensitive to astrophysical processes and typically difficult to model with galaxy formation simulations. We explore the fraction of cool-core (CC) clusters in a large sample of $370$ clusters from IllustrisTNG, examining six common CC definitions. IllustrisTNG produces continuous CC criteria distributions, the extremes of which are classified as CC and non-cool-core (NCC), and the criteria are increasingly correlated for more massive clusters. At $z=0$ the CC fraction is systematically lower than observed for the complete sample, selecting massive systems increases the CC fraction for $3$ criteria and reduces it for others. This result is partly driven by systematic differences between the simulated and observed gas fraction profiles. The simulated CC fraction increases more rapidly with redshift than observed, independent of mass or redshift range, and the CC fraction is overpredicted at $z\geq1$. The conversion of CCs to NCCs begins later and acts more rapidly in the simulations. Examining the fraction of CCs and NCCs defined as relaxed we find no evidence that CCs are more relaxed, suggesting that mergers are not solely responsible for disrupting CCs. A comparison of the median thermodynamic profiles defined by different CC criteria shows that the extent to which they evolve in the cluster core is dependent on the CC criteria. We conclude that the thermodynamic structure of galaxy clusters in IllustrisTNG shares many similarities with observations, but achieving better agreement most likely requires modifications of the underlying galaxy formation model.
Far-infrared and (sub)millimeter fluxes can be used to study dust in protoplanetary disks, the building blocks of planets. Here, we combine observations from the Herschel Space Observatory with ancillary data of 284 protoplanetary disks in the Taurus, Chamaeleon I, and Ophiuchus star-forming regions, covering from the optical to mm/cm wavelengths. We analyze their spectral indices as a function of wavelength and determine their (sub)millimeter slopes when possible. Most disks display observational evidence of grain growth, in agreement with previous studies. No correlation is found between other tracers of disk evolution and the millimeter spectral indices. A simple disk model is used to fit these sources, and we derive posterior distributions for the optical depth at 1.3 mm and 10 au, the disk temperature at this same radius, and the dust opacity spectral index. We find the fluxes at 70 microns to correlate strongly with disk temperatures at 10 au, as derived from these simple models. We find tentative evidence for spectral indices in Chamaeleon I being steeper than those of disks in Taurus/Ophiuchus, although more millimeter observations are needed to confirm this trend and identify its possible origin. Additionally, we determine the median spectral energy distribution of each region and find them to be similar across the entire wavelength range studied, possibly due to the large scatter in disk properties and morphologies.
We have studied the chemical enrichment history of the interstellar medium through an analysis of the n-dimensional stellar abundances space. This work is a non-parametric analysis of the stellar chemical abundance space. The main goal is to study the stars from their organization within this abundance space. Within this space, we seek to find clusters (in a statistical sense), that is, stars likely to share similar chemo-evolutionary history, using two methods: the hierarchical clustering and the principal component analysis. We analysed some selected abundance surveys available in the literature. For each sample, we labelled the group of stars according to its average abundance curve. In all samples, we identify the existence of a main enrichment pattern of the stars, which we call chemical enrichment flow. This flow is set by the structured and well defined mean rate at which the abundances of the interstellar medium increase, resulting from the mixture of the material ejected from the stars and stellar mass loss and interstellar medium gas. One of the main results of our analysis is the identification of subgroups of stars with peculiar chemistry. These stars are situated in regions outside of the enrichment flow in the abundance space. These peculiar stars show a mismatch in the enrichment rate of a few elements, such as Mg, Si, Sc and V, when compared to the mean enrichment rate of the other elements of the same stars. We believe the existence of these groups of stars with peculiar chemistry may be related to the accretion of planetary material onto stellar surfaces or may be due to production of the same chemical element by different nucleosynthetic sites.
We have developed a new X-ray absorption model, called {\tt IONeq}, which computes the optical depth $\tau(E)$ simultaneously for ions of all abundant elements, assuming ionization equilibrium and taking into account turbulent broadening. We use this model to analyze the interstellar medium (ISM) absorption features in the Milky Way for a sample of 18 galactic (LMXBs) and 42 extragalactic sources (mainly Blazars). The absorbing ISM was modeled as a combination of three components/phases - neutral ($T\lesssim1\times10^{4}$ K), warm ($T\sim5\times 10^{4}$ K) and hot ($T\sim2\times10^{6}$ K). We found that the spatial distribution of both, neutral and warm components, are difficult to describe using smooth profiles due to nonuniform distribution of the column densities over the sky. For the hot phase we used a combination of a flattened disk and a halo, finding comparable column densities for both spatial components, in the order of $\sim 6-7\times10^{18}\;{\rm cm^{-2}}$, although this conclusion depends on the adopted parametrization. If the halo component has sub-solar abundance $Z$, then the column density has be scaled up by a factor $\frac{Z_\odot}{Z}$. The vertically integrated column densities of the disk components suggests the following mass fractions for these three ISM phases in the Galactic disk: neutral $\sim~89\%$, warm $\sim 8\%$ and hot $\sim 3\%$ components, respectively. The constraints on the radial distribution of the halo component of the hot component are weak.
Complex organic molecules (COMs) have been observed in comets, hot cores and cold dense regions of the interstellar medium. It is generally accepted that these COMs form on icy dust grain through the recombination reaction of radicals triggered by either energetic UV- photon or non-energetic H-atom addition processing. In this work, we present for the first time laboratory studies that allow for quantitative comparison of hydrogenation and UV-induced reactions as well as their cumulative effect in astronomically relevant CO:CH3OH=4:1 ice analogues. The formation of glycolaldehyde (GA) and ethylene glycol (EG) is confirmed in pure hydrogenation experiments at 14 K, except methyl formate (MF), which is only clearly observed in photolysis. The fractions for MF:GA:EG are 0 : (0.2-0.4) : (0.8-0.6) for pure hydrogenation, and 0.2 : 0.3 : 0.5 for UV involving experiments and can offer a diagnostic tool to derive the chemical origin of these species. The GA/EG ratios in the laboratory (0.3-1.5) compare well with observations toward different objects.
Links to: arXiv, form interface, find, astro-ph, recent, 1710, contact, help (Access key information)
Super-Massive Black Holes weighing up to $\sim 10^9 \, \mathrm{M_{\odot}}$ are in place by $z \sim 7$, when the age of the Universe is $\lesssim 1 \, \mathrm{Gyr}$. This implies a time crunch for their growth since such high masses cannot be easily reached in standard accretion scenarios. Here, we explore the physical conditions that would lead to optimal growth wherein stable super-Eddington accretion would be permitted. Our analysis suggests that the preponderance of optimal conditions depends on two key parameters: the black hole mass and the host galaxy central gas density. In the high-efficiency region of this parameter space, a continuous stream of gas can accrete onto the black hole from large to small spatial scales, assuming a global isothermal profile for the host galaxy. By using analytical initial mass functions for black hole seeds, we find an enhanced probability of high-efficiency growth for seeds with initial masses $\gtrsim 10^4 \, \mathrm{M_{\odot}}$. Our picture suggests that a large population of high-$z$ lower-mass black holes that formed in the low-efficiency region, with low duty cycles and accretion rates, might remain undetectable as quasars, since we predict their bolometric luminosities to be $\lesssim 10^{41} \, \mathrm{erg \, s^{-1}}$. The presence of these sources might be revealed only via gravitational wave detections of their mergers.
The formation and evolution of gravitationally bound, star forming substructures in tidal tails of interacting galaxies, called tidal dwarf galaxies (TDG), has been studied, until now, only in idealised simulations of individual pairs of interacting galaxies for pre-determined orbits, mass ratios, and gas fractions. Here, we present the first identification of TDG candidates in fully cosmological simulations, specifically the high-resolution simulations of the EAGLE suite. The finite resolution of the simulation limits their ability to predict the exact formation rate and survival timescale of TDGs, but we show that gravitationally bound baryonic structures in tidal arms already form in current state-of-the-art cosmological simulations. In this case, the orbital parameter, disc orientations as well as stellar and gas masses and the specific angular momentum of the TDG forming galaxies are a direct consequence of cosmic structure formation. We identify TDG candidates in a wide range of environments, such as multiple galaxy mergers, clumpy high-redshift (up to z = 2) galaxies, high-speed encounters, and tidal interactions with gas-poor galaxies. We present selection methods, the properties of the identified TDG candidates and a roadmap for more quantitative analyses using future high-resolution simulations.
After changing optical AGN type from 1.9 to 1 in 1984, the AGN Mrk 1018 recently reverted back to its type 1.9 state. Our ongoing monitoring now reveals that the AGN has halted its dramatic dimming, reaching a minimum around October 2016. The minimum was followed by an outburst rising with $\sim$0.25 U-band mag/month. The rebrightening lasted at least till February 2017 as confirmed by joint Chandra and Hubble observations. Monitoring was resumed in July 2017 after the source emerged from sunblock, at which point the AGN was found only $\sim$0.4 mag brighter than its minimum. The intermittent outburst was accompanied with the appearance of a red wing asymmetry in broad-line shape, indicative of an inhomogeneous broad line region. Mrk 1018's current flickering brightness following its rapid fade suggests either that the source has reignited, remains variable at a low level, or may continue dimming over the next few years. Discriminating between these possibilities demands continual multi-wavelength monitoring.
In this paper, we investigate 2727 galaxies observed by MaNGA as of June 2016 to develop spatially resolved techniques for identifying signatures of active galactic nuclei (AGN). We identify 303 AGN candidates. The additional spatial dimension imposes challenges in identifying AGN due to contamination from diffuse ionized gas, extra-planar gas and photoionization by hot stars. We show that the combination of spatially-resolved line diagnostic diagrams and additional cuts on H$\alpha$ surface brighness and H$\alpha$ equivalent width can distinguish between AGN-like signatures and high-metallicity galaxies with LINER-like spectra. Low mass galaxies with high specific star formation rates are particularly difficult to diagnose and routinely show diagnostic line ratios outside of the standard star-formation locus. We develop a new diagnostic -- the distance from the standard diagnostic line in the line-ratios space -- to evaluate the significance of the deviation from the star-formation locus. We find 173 galaxies that would not have been selected as AGN candidates based on single-fibre spectral measurements but exhibit photoionization signatures suggestive of AGN activity in the MaNGA resolved observations, underscoring the power of large integral field unit (IFU) surveys. A complete census of these new AGN candidates is necessary to understand their nature and probe the complex co-evolution of supermassive black holes and their hosts.
We presents results from a spectroscopic survey of $z\sim5$ quasars in the CFHT Legacy Survey (CFHTLS). Using both optical color selection and a likelihood method we select 97 candidates over an area of 105 deg$^2$ to a limit of $i_{\rm AB} < 23.2$, and 7 candidates in the range $23.2 < i_{\rm AB} < 23.7$ over an area of 18.5 deg$^2$. Spectroscopic observations for 43 candidates were obtained with Gemini, MMT, and LBT, of which 37 are $z>4$ quasars. This sample extends measurements of the quasar luminosity function $\sim$1.5 mag fainter than our previous work in SDSS Stripe 82. The resulting luminosity function is in good agreement with our previous results, and suggests that the faint end slope is not steep. We perform a detailed examination of our survey completeness, particularly the impact of the Ly$\alpha$ emission assumed in our quasar spectral models, and find hints that the observed Ly$\alpha$ emission from faint $z\sim5$ quasars is weaker than for $z\sim3$ quasars at a similar luminosity. Our results strongly disfavor a significant contribution of faint quasars to the hydrogen-ionizing background at $z=5$.
SubHalo Abundance Matching (SHAM) assumes that one (sub)halo property, such as mass Mvir or peak circular velocity Vpeak, determines properties of the galaxy hosted in each (sub)halo such as its luminosity or stellar mass. This assumption implies that the dependence of Galaxy Luminosity Functions (GLFs) and the Galaxy Stellar Mass Function (GSMF) on environmental density is determined by the corresponding halo density dependence. In this paper, we test this by determining from an SDSS sample the observed dependence with environmental density of the ugriz GLFs and GSMF for all galaxies, and for central and satellite galaxies separately. We then show that the SHAM predictions are in remarkable agreement with these observations, even when the galaxy population is divided between central and satellite galaxies. However, we show that SHAM fails to reproduce the correct dependence between environmental density and g-r color for all galaxies and central galaxies, although it better reproduces the color dependence on environmental density of satellite galaxies.
Previous studies suggest that the growth of supermassive black holes (SMBHs) may be fundamentally related to host-galaxy stellar mass ($M_\star$). To investigate this SMBH growth-$M_\star$ relation in detail, we calculate long-term SMBH accretion rate as a function of $M_\star$ and redshift [$\overline{\rm BHAR}(M_\star, z)$] over ranges of $\log(M_\star/M_\odot)=\text{9.5--12}$ and $z=\text{0.4--4}$. Our $\overline{\rm BHAR}(M_\star, z)$ is constrained by high-quality survey data (GOODS-South, GOODS-North, and COSMOS), and by the stellar mass function and the X-ray luminosity function. At a given $M_\star$, $\overline{\rm BHAR}$ is higher at high redshift. This redshift dependence is stronger in more massive systems (for $\log(M_\star/M_\odot)\approx 11.5$, $\overline{\rm BHAR}$ is three decades higher at $z=4$ than at $z=0.5$), possibly due to AGN feedback. Our results indicate that the ratio between $\overline{\rm BHAR}$ and average star formation rate ($\overline{\rm SFR}$) rises toward high $M_\star$ at a given redshift. This $\overline{\rm BHAR}/\overline{\rm SFR}$ dependence on $M_\star$ does not support the scenario that SMBH and galaxy growth are in lockstep. We calculate SMBH mass history [$M_{\rm BH}(z)$] based on our $\overline{\rm BHAR}(M_\star, z)$ and the $M_\star(z)$ from the literature, and find that the $M_{\rm BH}$-$M_\star$ relation has weak redshift evolution since $z\approx 2$. The $M_{\rm BH}/M_\star$ ratio is higher toward massive galaxies: it rises from $\approx 1/5000$ at $\log M_\star\lesssim 10.5$ to $\approx 1/500$ at $\log M_\star \gtrsim 11.2$. Our predicted $M_{\rm BH}/M_\star$ ratio at high $M_\star$ is similar to that observed in local giant ellipticals, suggesting that SMBH growth from mergers is unlikely to dominate over growth from accretion.
We present Lightning, a new spectral energy distribution (SED) fitting procedure, capable of quickly and reliably recovering star formation history (SFH) and extinction parameters. The SFH is modeled as discrete steps in time. In this work, we assumed lookback times of 0-10 Myr, 10-100 Myr, 0.1-1 Gyr, 1-5 Gyr, and 5-13.6 Gyr. Lightning consists of a fully vectorized inversion algorithm to determine SFH step intensities and combines this with a grid-based approach to determine three extinction parameters. We apply our procedure to the extensive FUV-to-FIR photometric data of M51, convolved to a common spatial resolution and pixel scale, and make the resulting maps publicly available. We recover, for M51a, a peak star formation rate (SFR) between 0.1 and 5 Gyr ago, with much lower star formation activity over the last 100 Myr. For M51b, we find a declining SFR toward the present day. In the outskirt regions of M51a, which includes regions between M51a and M51b, we recover a SFR peak between 0.1 and 1 Gyr ago, which corresponds to the effects of the interaction between M51a and M51b. We utilize our results to (1) illustrate how UV+IR hybrid SFR laws vary across M51, and (2) provide first-order estimates for how the IR luminosity per unit stellar mass varies as a function of the stellar age. From the latter result, we find that IR emission from dust heated by stars is not always associated with young stars, and that the IR emission from M51b is primarily powered by stars older than 5 Gyr.
We have obtained the time and space-resolved star formation history (SFH) of M51a (NGC 5194) by fitting GALEX, SDSS, and near infrared pixel-by-pixel photometry to a comprehensive library of stellar population synthesis models drawn from the Synthetic Spectral Atlas of Galaxies (SSAG). We fit for each space-resolved element (pixel) an independent model where the SFH is averaged in 137 age bins, each one 100 Myr wide. We used the Bayesian Successive Priors (BSP) algorithm to mitigate the bias in the present-day spatial mass distribution. We test BSP with different prior probability distribution functions (PDFs); this exercise suggests that the best prior PDF is the one concordant with the spatial distribution of the stellar mass as inferred from the near infrared images. We also demonstrate that varying the implicit prior PDF of the SFH in SSAG does not affects the results. By summing the contributions to the global star formation rate of each pixel, at each age bin, we have assembled the resolved star formation history of the whole galaxy. According to these results, the star formation rate of M51a was exponentially increasing for the first 10 Gyr after the Big Bang, and then turned into an exponentially decreasing function until the present day. Superimposed, we find a main burst of star formation at t 11.9 Gyr after the Big Bang.
We present a new technique for empirically calibrating how the X-ray luminosity function (XLF) of X-ray binary (XRB) populations evolves following a star-formation event. We first utilize detailed stellar population synthesis modeling of far-UV to far-IR photometry of the nearby face-on spiral galaxy M51 to construct maps of the star-formation histories (SFHs) on subgalactic (~400 pc) scales. Next, we use the ~850 ks cumulative Chandra exposure of M51 to identify and isolate 2-7 keV detected point sources within the galaxy, and we use our SFH maps to recover the local properties of the stellar populations in which each X-ray source is located. We then divide the galaxy into various subregions based on their SFH properties (e.g., star-formation rate [SFR] per stellar mass [M*] and mass-weighted stellar age) and group the X-ray point sources according to the characteristics of the regions in which they are found. Finally, we construct and fit a parameterized XLF model that quantifies how the XLF shape and normalization evolves as a function of the XRB population age. Our best-fit model indicates the XRB XLF per unit stellar mass declines in normalization, by ~3-3.5 dex, and steepens in slope from ~10 Myr to ~10 Gyr. We find that our technique recovers results from past studies of how XRB XLFs and XRB luminosity scaling relations vary with age and provides a self-consistent picture for how the XRB XLF evolves with age.
We investigate the contribution of outer HI disks to the observable population of merging black hole binaries. Like dwarf galaxies, the outer HI disks of spirals have low star formation rates and lower metallicities than the inner disks of spirals. Since low-metallicity star formation can produce more detectable compact binaries than typical star formation, the environments in the outskirts of spiral galaxies may be conducive to producing a rich population of massive binary black holes. We consider here both detailed controlled simulations of spirals and cosmological simulations, as well as the current range of observed values for metallicity and star formation in outer disks. We find that outer HI disks contribute at least as much as dwarf galaxies do to the observed LIGO/Virgo detection rates. Identifying the host galaxies of merging massive black holes should provide constraints on cosmological parameters and insights into the formation channels of binary mergers.
We present results of hydrodynamic simulations of massive star forming regions with and without protostellar jets. We show that jets change the normalization of the stellar mass accretion rate, but do not strongly affect the dynamics of star formation. In particular, $M_*(t) \propto f^2 (t-t_*)^2$ where $f = 1 - f_{\rm jet}$ is the fraction of mass accreted onto the protostar, $f_{\rm jet}$ is the fraction ejected by the jet, and $(t-t_*)^2$ is the time elapsed since the formation of the first star. The star formation efficiency is nonlinear in time. We find that jets have only a small effect (of order 25\%) on the accretion rate onto the protostellar disk (the "raw" accretion rate). We show that the small scale structure -- the radial density, velocity, and mass accretion profiles are very similar in the jet and no-jet cases. Finally, we show that the inclusion of jets does drive turbulence but only on small (parsec) scales.
Massive black hole binaries (BHBs) are expected to form as the result of galaxy mergers; they shrink via dynamical friction and stellar scatterings, until gravitational waves (GWs) bring them to the final coalescence. It has been argued that BHBs may stall at a parsec scale and never enter the GW stage if stars are not continuously supplied to the BHB loss cone. Here we perform several N-body experiments to study the effect of an 80,000 solar masses stellar cluster (SC) infalling on a parsec-scale BHB. We explore different orbital elements for the SC and we perform runs both with and without accounting for the influence of a rigid stellar cusp (modelled as a rigid Dehnen potential). We find that the semi-major axis of the BHB shrinks by more than 10 per cent if the SC is on a nearly radial orbit; the shrinking is more efficient when a Dehnen potential is included and the orbital plane of the SC coincides with that of the BHB. In contrast, if the SC orbit has non-zero angular momentum, only a few stars enter the BHB loss cone and the resulting BHB shrinking is negligible. Our results indicate that SC disruption might significantly contribute to the shrinking of a parsec-scale BHB only if the SC approaches the BHB on a nearly radial orbit.
We address in detail the radiation forces on spherical dust envelopes around luminous stars, and numerical solutions for these forces, as a first step toward more general dust geometries. Two physical quantities, a normalized force and a force-averaged radius, suffice to capture the overall effects of radiation forces. In addition to the optically thin and thick regimes, the wavelength dependence of dust opacity allows for an inter- mediate case in which starlight is easily trapped but infrared radiation readily escapes. Scattering adds a non-negligible force in this intermediate regime. We address all three regimes analytically and provide approximate formulae for the force parameters, for arbitrary optical depth and inner dust temperature. Turning to numerical codes, we examine the convergence properties of DUSTY and of the Monte Carlo code Hyperion. We find that Monte Carlo codes tend to underestimate the radiation force when the mean free path of starlight is not well resolved, as this causes the inner dust temperature, and therefore the inner Rosseland opacity, to be too low. We briefly discuss implications for more complicated radiation transfer problems.
Observations of luminous quasars have detected fast (~1000 km s$^{-1}$) outflows of molecular gas. We recently used a series of hydro-chemical simulations to demonstrate that these molecular outflows can be explained by the formation of molecules within the AGN wind material itself. However, such simulations are computationally expensive, which limits the physical parameter space that we can explore. We have therefore developed an analytic model to follow the radiative cooling of the shocked ISM layer of a spherically symmetric AGN wind. We use this model to explore a wide range of ambient densities ($1-10^{4}$ cm$^{-3}$), density profile slopes ($0-1.5$), AGN luminosities ($10^{44}-10^{47}$ erg s$^{-1}$), and metallicities ($0.1-3$ Z$_{\odot}$). The analytic models mostly cool within ~1 Myr, and are therefore likely to produce observable molecular outflows. The momentum boost, calculated from the instantaneous rate of change of the outflow momentum, initially increases as the outflow decelerates, as expected for an energy-conserving outflow. However, it reaches a maximum of $\approx$20, due to the work done against the gravitational potential of the black hole and host galaxy. A time-averaged observational estimate of the momentum boost, calculated from the outflow velocity, radius and molecular mass, reaches a maximum of $\approx 1-2$. This is partly due to our assumed molecular fraction, 0.2, but we also show that the instantaneous and observationally time-averaged definitions are not equivalent. Thus recent observations that estimated momentum boosts of order unity do not necessarily rule out an energy-driven outflow.
We present a comprehensive and precise description of the Sagittarius (Sgr) stellar stream's 3D geometry as traced by its old stellar population. This analysis draws on the sample of ${\sim}44,000$ RR Lyrae (RRab) stars from the Pan-STARRS1 (PS1) 3$\pi$ survey (Hernitschek et al. 2016,Sesar et al. 2017b), which is ${\sim}80\%$ complete and ${\sim}90\%$ pure within 80~kpc, and extends to ${\gtrsim} 120$~kpc with a distance precision of ${\sim} 3\%$. A projection of RR Lyrae stars within $|\tilde{B}|_{\odot}<9^\circ$ of the Sgr stream's orbital plane reveals the morphology of both the leading and the trailing arms at very high contrast, across much of the sky. In particular, the map traces the stream near-contiguously through the distant apocenters. We fit a simple model for the mean distance and line-of-sight depth of the Sgr stream as a function of the orbital plane angle $\tilde{\Lambda}_{\odot}$, along with a power-law background-model for the field stars. This modeling results in estimates of the mean stream distance precise to ${\sim}1\%$ and it resolves the stream's line-of-sight depth. These improved geometric constraints can serve as new constraints for dynamical stream models.
Radio jets play an important role in quasar feedback, but direct observations showing how the jets interact with the multi-phase interstellar medium of galaxy disks are few and far between. In this work, we provide new millimeter interferometric observations of PG 1700+518 in order to investigate the effect of its radio jet on the surrounding molecular gas. PG 1700 is a radio-quiet, low-ionization broad absorption line quasar whose host galaxy has a nearby interacting companion. On sub-kiloparsec scales, the ionized gas is driven to high velocities by a compact radio jet that is identified by radio interferometry. We present observations from the NOrthern Extended Millimeter Array (NOEMA) interferometer with a 3.8 arcsec (15.96 kpc) synthesized beam where we detect the CO (1-0) emission line at $30\sigma$ significance with a total flux of $3.12\pm0.02$ Jy km s$^{-1}$ and a typical velocity dispersion of $125\pm5$ km s$^{-1}$. Despite the outflow in ionized gas, we find no concrete evidence that the CO gas is being affected by the radio jet on size scales of a kiloparsec or more. However, a $\sim\!1$ arcsec drift in the spatial centroid of the CO emission as a function of velocity across the emission line and the compact nature of the jet hint that higher spatial resolution observations may reveal a signal of interaction between the jet and molecular gas.
We present Keck/OSIRIS adaptive optics observations with 150-400 pc spatial sampling of 7 turbulent, clumpy disc galaxies from the DYNAMO sample ($0.07<z<0.2$). DYNAMO galaxies have previously been shown to be well matched in properties to main sequence galaxies at $z\sim1.5$. Integral field spectroscopy observations using adaptive optics are subject to a number of systematics including a variable PSF and spatial sampling, which we account for in our analysis. We present gas velocity dispersion maps corrected for these effects, and confirm that DYNAMO galaxies do have high gas velocity dispersion ($\sigma=40-80$\kms), even at high spatial sampling. We find statistically significant structure in 6 out of 7 galaxies. The most common distance between the peaks in velocity dispersion and emission line peaks is $\sim0.5$~kpc, we note this is very similar to the average size of a clump measured with HST H$\alpha$ maps. This could suggest that the peaks in velocity dispersion in clumpy galaxies likely arise due to some interaction between the clump and the surrounding ISM of the galaxy, though our observations cannot distinguish between outflows, inflows or velocity shear. Observations covering a wider area of the galaxies will be needed to confirm this result.
We present rest-frame ultraviolet and optical spectroscopy of the brightest lensed galaxy yet discovered, at redshift z = 2.4. This source reveals a characteristic, triple-peaked Lyman {\alpha} profile which has been predicted by various theoretical works but to our knowledge has not been unambiguously observed previously. The feature is well fit by a superposition of two components: a double-peak profile emerging from substantial radiative transfer, and a narrow, central component resulting from directly escaping Lyman {\alpha} photons; but is poorly fit by either component alone. We demonstrate that the feature is unlikely to contain contamination from nearby sources, and that the central peak is unaffected by radiative transfer effects apart from very slight absorption. The feature is detected at signal-to-noise ratios exceeding 80 per pixel at line center, and bears strong resemblance to synthetic profiles predicted by numerical models.
We present an analysis of stellar populations in passive galaxies in seven massive X-ray clusters at z=0.19-0.89. Based on absorption line strengths measured from our high signal-to-noise spectra, the data support primarily passive evolution of the galaxies. We use the scaling relations between velocity dispersions and the absorption line strengths to determine representative mean line strengths for the clusters. From the age determinations based on the line strengths (and stellar population models), we find a formation redshift of z_form=1.96(-0.19,+0.24). Based on line strength measurements from high signal-to-noise composite spectra of our data, we establish the relations between velocity dispersion, ages, metallicities [M/H] and abundance ratios [alpha/Fe] as a function of redshift. The [M/H]-velocity dispersion and [alpha/Fe]-velocity dispersion relations are steep and tight. The age-velocity dispersion relation is flat, with zero point changes reflecting passive evolution. The scatter in all three parameters are within 0.08-0.15 dex at fixed velocity dispersions, indicating a large degree of synchronization in the evolution of the galaxies. We find indication of cluster-to-cluster differences in metallicities and abundance ratios. However, variations in stellar populations with the cluster environment can only account for a very small fraction of the intrinsic scatter in the scaling relations. Thus, within these very massive clusters the main driver of the properties of the stellar populations in passive galaxies appears to be the galaxy velocity dispersion.
We here report an identification of SDSS J090152.04+624342.6 as a new "overlapping-trough" iron low-ionization broad absorption line quasar at redshift of $z\sim2.1$. No strong variation of the broad absorption lines can be revealed through the two spectra taken by the Sloan Digital Sky Survey with a time interval of $\sim6$yr. Further optical and infrared spectroscopic study on this object is suggested.
Previously, it has been established that axion dark matter (DM) is clustered to form clumps (axion miniclusters) with masses $M\sim10^{-12}M_\odot$. The passages of such clumps through the Earth are very rare events occurring once in $10^5$ years. It has also been shown that the Earth's passage through DM streams, which are the remnants of clumps destroyed by tidal gravitational forces from Galactic stars, is a much more probable event occurring once in several years. In this paper we have performed details calculations of the destruction of miniclusters by taking into account their distribution in orbits in the Galactic halo. We have investigated two DM halo models, the Navarro-Frenk-White and isothermal density profiles. Apart from the Galactic disk stars, we have also taken into account the halo and bulge stars. We show that about 2-5% of the axion miniclusters are destroyed when passing near stars and transform into axion streams, while the clump destruction efficiency depends on the DM halo model. The expected detection rate of streams with an overdensity exceeding an order of magnitude is 1-2 in 20 years. The possibility of detecting streams by their tidal gravitational effect on gravitational-wave interferometers is also considered.
Aims: To reveal the morphology, chemical composition, kinematics and to
establish the main processes prevalent in the gas at the foot points of the
giant molecular loops (GMLs) in the Galactic center region
Methods: Using the 22-m Mopra telescope, we mapped the M$-3.8+0.9$ molecular
cloud, placed at the foot points of a giant molecular loop, in 3-mm range
molecular lines. To derive the molecular hydrogen column density, we also
observed the $^{13}$CO $(2-1)$ line at 1 mm using the 12-m APEX telescope. From
the 3 mm observations 12 molecular species were detected, namely HCO$^+$, HCN,
H$^{13}$CN, HNC, SiO, CS, CH$_3$OH, N$_2$H$^+$, SO, HNCO, OCS, and HC$_3$N.
Results: Maps revealing the morphology and kinematics of the M$-3.8+0.9$
molecular cloud in different molecules are presented. We identified six main
molecular complexes. We derive fractional abundances in 11 selected positions
of the different molecules assuming local thermodynamical equilibrium.
Conclusions: Most of the fractional abundances derived for the M$-3.8+0.9$
molecular cloud are very similar over the whole cloud. However, the fractional
abundances of some molecules show significant difference with respect to those
measured in the central molecular zone (CMZ). The abundances of the shock
tracer SiO are very similar between the GMLs and the CMZ. The methanol emission
is the most abundant specie in the GMLs. This indicates that the gas is likely
affected by moderate $\sim $ 30 km s$^{-1}$ or even high velocity (50 km
s$^{-1}$) shocks, consistent with the line profile observed toward one of the
studied position. The origin of the shocks is likely related to the flow of the
gas throughout the GMLs towards the foot points.
ALMA is providing us essential information on where certain molecules form. Observing where these molecules emission arises from, the physical conditions of the gas, and how this relates with the presence of other species allows us to understand the formation of many species, and to significantly improve our knowledge of the chemistry that occurs in the space. We studied the molecular distribution of NaCN around IRC +10216, a molecule detected previously, but whose origin is not clear. High angular resolution maps allow us to model the abundance distribution of this molecule and check suggested formation paths. We modeled the emission of NaCN assuming local thermal equilibrium (LTE) conditions. These profiles were fitted to azimuthal averaged intensity profiles to obtain an abundance distribution of NaCN. We found that the presence of NaCN seems compatible with the presence of CN, probably as a result of the photodissociation of HCN, in the inner layers of the ejecta of IRC +10216. However, similar as for CH 3 CN, current photochemical models fail to reproduce this CN reservoir. We also found that the abundance peak of NaCN appears at a radius of 3 x 10 15 cm, approximately where the abundance of NaCl, suggested to be the parent species, starts to decay. However, the abundance ratio shows that the NaCl abundance is lower than that obtained for NaCN. We expect that the LTE assumption might result in NaCN abundances higher than the real ones. Updated photochemical models, collisional rates, and reaction rates are essential to determine the possible paths of the NaCN formation.
We present a deep, low-frequency radio continuum study of the nearby Fanaroff--Riley class I (FR I) radio galaxy 3C 31 using a combination of LOw Frequency ARray (LOFAR; 30--85 and 115--178 MHz), Very Large Array (VLA; 290--420 MHz), Westerbork Synthesis Radio Telescope (WSRT; 609 MHz) and Giant Metre Radio Telescope (GMRT; 615 MHz) observations. Our new LOFAR 145-MHz map shows that 3C 31 has a largest physical size of $1.1$ Mpc in projection, which means 3C 31 now falls in the class of giant radio galaxies. We model the radio continuum intensities with advective cosmic-ray transport, evolving the cosmic-ray electron population and magnetic field strength in the tails as functions of distance to the nucleus. We find that if there is no in-situ particle acceleration in the tails, then decelerating flows are required that depend on radius $r$ as $v\propto r^{\beta}$ ($\beta\approx -1$). This then compensates for the strong adiabatic losses due to the lateral expansion of the tails. We are able to find self-consistent solutions in agreement with the entrainment model of Croston & Hardcastle, where the magnetic field provides $\approx$$1/3$ of the pressure needed for equilibrium with the surrounding intra-cluster medium (ICM). We obtain an advective time-scale of $\approx$$190$ Myr, which, if equated to the source age, would require an average expansion Mach number ${\cal M} \approx 5$ over the source lifetime. Dynamical arguments suggest that instead, either the outer tail material does not represent the oldest jet plasma or else the particle ages are underestimated due to the effects of particle acceleration on large scales.
We obtain a more accurate statistical estimation of the mass of double galaxies moving in circular orbits, including confidence intervals for different confidence levels.
Links to: arXiv, form interface, find, astro-ph, recent, 1710, contact, help (Access key information)