We analyze the properties of the circumgalactic gas (CGM) around 120 galaxies with stellar and dark matter halo masses similar to that of the Milky Way. We focus on the morphology and kinematics of the neutral hydrogen and how this depends on f_g, the ratio of gas-to-stellar mass within the optical radius. In gas-rich galaxies with f_g > 0.1, gas temperatures rise slowly from center of the halo out to the virial radius and average neutral gas column densities remain above 10^19 atoms cm^-2 out to radii of 50-70 kpc. In gas-poor galaxies with f_g < 0.1, gas temperatures rise quickly outside the edge of the disk to 10^6 K, and then remain fixed out to radii of 100 kpc. The column density of neutral gas quickly drops below 10^19 atoms cm^-2 at radii of 10 kpc. Neutral gas distributions are also more asymmetric in gas-poor galaxies. Most of the differences between gas-poor and gas-rich galaxies in the Illustris simulation can be attributed to the effects of "radio-mode" AGN feedback. In the Illustris simulation, the circumgalactic gas is found to rotate coherently about the center of the galaxy with a maximum rotational velocity of around 200 km/s. In gas-rich galaxies, the average coherence length of the rotating gas is 40 kpc, compared to 10 kpc in gas-poor galaxies. In the very most gas-rich systems, the CGM can rotate coherently over scales of 70-100 kpc. We discuss our results in the context of recent observations of the CGM in low mass galaxies via UV absorption-line spectroscopy and deep 21cm observations of edge-on spiral galaxies.
Using atomic hydrogen (HI) data from the VLA Galactic Plane Survey we measure the HI terminal velocity as a function of longitude for the first quadrant of the Milky Way. We use these data, together with our previous work on the fourth Galactic quadrant, to produce a densely sampled, uniformly measured, rotation curve of the Northern and Southern Milky Way between $3~{\rm kpc} < R < 8~{\rm kpc}$. We determine a new joint rotation curve fit for the first and fourth quadrants, which is consistent with the fit we published in McClure-Griffiths \& Dickey (2007) and can be used for estimating kinematic distances interior to the solar circle. Structure in the rotation curves is now exquisitely well defined, showing significant velocity structure on lengths of $\sim 200$ pc, which is much greater than the spatial resolution of the rotation curve. Furthermore, the shape of the rotation curves for the first and fourth quadrants, even after subtraction of a circular rotation fit shows a surprising degree of correlation with a roughly sinusoidal pattern between $4.2 < R < 7$ kpc.
We examined the reliability of estimates of pseudoisothermal, Burkert and NFW
dark halo parameters for the methods based on the mass-modelling of the
rotation curves.To do it we constructed the $\chi^2$ maps for the grid of the
dark matter halo parameters for a sample of 14 disc galaxies with high quality
rotation curves from THINGS. We considered two variants of models in which: a)
the mass-to-light ratios of disc and bulge were taken as free parameters, b)
the mass-to-light ratios were fixed in a narrow range according to the models
of stellar populations.
To reproduce the possible observational features of the real galaxies we made
tests showing that the parameters of the three halo types change critically in
the cases of a lack of kinematic data in the central or peripheral areas and
for different spatial resolutions.
We showed that due to the degeneracy between the central densities and the
radial scales of the dark haloes there are considerable uncertainties of their
concentrations estimates. Due to this reason it is also impossible to draw any
firm conclusion about universality of the dark halo column density based on
mass-modelling of even a high quality rotation curve. The problem is not solved
by fixing the density of baryonic matter.
In contrast, the estimates of dark halo mass within optical radius are much
more reliable. We demonstrated that one can evaluate successfully the halo mass
using the pure best-fitting method without any restrictions on the
mass-to-light ratios.
The Andromeda galaxy is observed to have a system of two large dwarf ellipticals and ~13 smaller satellite galaxies that are currently co-rotating in a thin plane, in addition to 2 counter-rotating satellite galaxies. We explored the consistency of those observations with a scenario where the majority of the co-rotating satellite galaxies originated from a subhalo group, where NGC 205 was the host and the satellite galaxies occupied dark matter sub-subhalos. We ran N-body simulations of a close encounter between NGC 205 and M31. In the simulations, NGC 205 was surrounded by massless particles to statistically sample the distribution of the sub-subhalos expected in a subhalo that has a mass similar to NGC 205. We made Monte Carlo samplings and found that, using a set of reference parameters, the probability of producing a thinner distribution of sub-subhalos than the observed NGC 205 + 15 smaller satellites (thus including the 2 counter-rotators, but excluding M32) increased from <1e-8 for the initial distribution to ~0.01 at pericentre. The probability of the simulated sub-subhalos occupying the locations of the observed co-rotating satellites in the line of sight velocity versus projected on-sky distance plane is at most 0.002 for 11 out of 13 satellites. Increasing the mass of M31 and the extent of the initial distribution of sub-subhalos gives a maximum probability of 0.004 for all 13 co-rotating satellites, but the probability of producing the thinness would drop to ~ 0.001.
We present eight monitoring seasons of the four brightest images of the gravitational lens SDSS J1004+4112 observed between December 2003 and October 2010. Using measured time delays for the images A, B and C and the model predicted time delay for image D we have removed the intrinsic quasar variability, finding microlensing events of about 0.5 and 0.7 mag of amplitude in the images C and D. From the statistics of microlensing amplitudes in images A, C, and D, we have inferred the half-light radius (at {\lambda} rest = 2407 {\AA}) for the accretion disk using two different methods, $R_{1/2}=8.7^{+18.5}_{-5.5} \sqrt{M/0.3 M_\odot}$ (histograms product) and $R_{1/2} = 4.2^{+3.2}_{-2.2} \sqrt{M/0.3 M_\odot}$ light-days ($\chi^2$). The results are in agreement within uncertainties with the size predicted from the black hole mass in SDSS J1004+4112 using the thin disk theory.
The circular rotation speed of the Milky Way at the solar radius, Theta_o, has been estimated to be 220 km/s by fitting the maximum velocity of HI emission as a function of Galactic longitude. This result is in tension with a recent estimate of Theta_o=240 km/s, based on VLBI parallaxes and proper motions from the BeSSeL and VERA surveys for large numbers of high-mass star forming regions across the Milky Way. We find that the rotation curve best fitted to the VLBI data is slightly curved, and that this curvature results in a biased estimate of Theta_o from the HI data when a flat rotation curve is assumed. This relieves the tension between the methods and favors Theta_o=240 km/s.
The hot intra-cluster medium (ICM) is rich in metals, which are synthesised by supernovae (SNe) and accumulate over time into the deep gravitational potential well of clusters of galaxies. Since most of the elements visible in X-rays are formed by type Ia (SNIa) and/or core-collapse (SNcc) supernovae, measuring their abundances gives us direct information on the nucleosynthesis products of billions of SNe since the epoch of the star formation peak (z~2-3). In this study, we compare the most accurate average X/Fe abundance ratios (compiled in a previous work from XMM-Newton EPIC and RGS observations of 44 galaxy clusters, groups, and ellipticals), representative of the chemical enrichment in the nearby ICM, to various SNIa and SNcc nucleosynthesis models found in the literature. The use of a SNcc model combined to any favoured standard SNIa model (deflagration or delayed-detonation) fails to reproduce our abundance pattern. In particular, the Ca/Fe and Ni/Fe ratios are significantly underestimated by the models. We show that the Ca/Fe ratio can be reproduced better, either by taking a SNIa delayed-detonation model that matches the observations of the Tycho supernova remnant, or by adding a contribution from the Ca-rich gap transient SNe, whose material should easily mix into the hot ICM. On the other hand, the Ni/Fe ratio can be reproduced better by assuming that both deflagration and delayed-detonation SNIa contribute in similar proportions to the ICM enrichment. In either case, the fraction of SNIa over the total number of SNe (SNIa+SNcc) contributing to the ICM enrichment ranges within 29-45%. This fraction is found to be systematically higher than the corresponding SNIa/SNe fraction contributing to the enrichment of the proto-solar environnement (15-25%). We also discuss and quantify two useful constraints on both SNIa and SNcc that can be inferred from the ICM abundance ratios.
The Event Horizon Telescope is a global very-long baseline interferometer capable of probing potential deviations from the Kerr metric, which is believed to provide the unique description of astrophysical black holes. Here we report an updated constraint on the quadrupolar deviation of Sagittarius A* within the context of a radiatively inefficient accretion flow model in a quasi-Kerr background. We also simulate near-future constraints obtainable by the forthcoming eight-station array and show that in this model already a one-day observation can measure the spin magnitude to within 0.005, the inclination to within 0.09{\deg}, the position angle to within 0.04{\deg}, and the quadrupolar deviation to within 0.005 at 3{\sigma} confidence. Thus, we are entering an era of high-precision strong gravity measurements.
We present the measurements of the luminosity-dependent redshift-space three-point correlation functions (3PCFs) for the Sloan Digital Sky Survey (SDSS) DR7 Main galaxy sample. We compare the 3PCF measurements to the predictions from three different halo and subhalo models. One is the halo occupation distribution (HOD) model and the other two are extensions of the subhalo abundance matching (SHAM) model by allowing the central and satellite galaxies to have different occupation distributions in the host halos and subhalos. Parameters in all the models are chosen to best describe the projected and redshift-space two-point correlation functions (2PCFs) of the same set of galaxies. All three model predictions agree well with the 3PCF measurements for the most luminous galaxy sample, while the HOD model better performs in matching the 3PCFs of fainter samples (with luminosity threshold below $L^*$), which is similar in trend to the case of fitting the 2PCFs. The decomposition of the model 3PCFs into contributions from different types of galaxy triplets shows that on small scales the dependence of the 3PCFs on triangle shape is driven by nonlinear redshift-space distortion (and not by the intrinsic halo shape) while on large scales it reflects the filamentary structure. The decomposition also reveals more detailed differences in the three models, which are related to the radial distribution, the mean occupation function, and the velocity distribution of satellite galaxies inside halos. The results suggest that galaxy 3PCFs can further help constrain the above galaxy-halo relation and test theoretical models.
We trace the peripheral magnetic field structure of Bok globule CB130, by estimating the linear polarization of its field stars in the R band. The magnetic field orientation sampled by these stars, aligned on average among themselves and the polarization produced within the cloud has a different direction from that of Galactic plane with an offset of 53$^\circ$. The offset between minor axis and the mean magnetic field of CB130 is found to be 80$^\circ$. The estimated strength of the magnetic field in the plane-of-the-sky is $\sim$ 116$\pm$19 $\mu$G. We constructed the visual extinction map using Near Infrared Color Excess (NICE) method to see the dust distribution around CB130. Contours of Herschel {Herschel is an ESA space observatory with science instruments provided by European-led Principal Investigator consortia and with important participation from NASA} SPIRE 500$\mu$m dust continuum emission maps of this cloud is over-plotted on the visual extinction map which shows the regions having higher optical extinction corresponds to higher densities of dust. Three distinct high dust density cores (named as C1, C2, and C3) are identified in the extinction map. It is observed that the cores C1 and C3 located close to two previously known cores CB130-1 and CB130-2 respectively. Estimate of visual extinction of some moderately obscured stars of CB130 are made utilizing near-infra red photometry. Its observed that there is a feeble dependence of polarization on extinction; and polarization efficiency (defined as $ \textit{p}$/A$ {V}$) of the dust grains decreases with the increase in extinction.
We present a new framework were we simultaneously fit strong lensing (SL) and dynamical data. The SL analysis is based on LENSTOOL, and the dynamical analysis uses MAMPOSSt code, which we have integrated into LENSTOOL. After describing the implementation of this new tool, we apply it on the galaxy group SL2S\,J02140-0535 ($z_{\rm spec}=0.44$), which we have already studied in the past. We use new VLT/FORS2 spectroscopy of multiple images and group members, as well as shallow X-ray data from \xmm. We confirm that the observed lensing features in SL2S\,J02140-0535 belong to different background sources. One of this sources is located at $z_{\rm spec}$ = 1.017 $\pm$ 0.001, whereas the other source is located at $z_{\rm spec}$ = 1.628 $\pm$ 0.001. With the analysis of our new and our previously reported spectroscopic data, we find 24 secure members for SL2S\,J02140-0535. Both data sets are well reproduced by a single NFW mass profile: the dark matter halo coincides with the peak of the light distribution, with scale radius, concentration, and mass equal to $r_s$ =$82^{+44}_{-17}$ kpc , $c_{200}$ = $10.0^{+1.7}_{-2.5}$, and $M_{200}$ = $1.0^{+0.5}_{-0.2}$ $\times$ 10$^{14}$M$_{\odot}$ respectively. These parameters are better constrained when we fit simultaneously SL and dynamical information. The mass contours of our best model agrees with the direction defined by the luminosity contours and the X-ray emission of SL2S\,J02140-0535. The simultaneous fit lowers the error in the mass estimate by 0.34 dex, when compared to the SL model, and in 0.15 dex when compared to the dynamical method.The combination of SL and dynamics tools yields a more accurate probe of the mass profile of SL2S\,J02140-0535 up to $r_{200}$. However, there is tension between the best elliptical SL model and the best spherical dynamical model.
Broad emission lines of active galactic nuclei stem from a spatially extended region (broad-line region; BLR) that are composed of discrete clouds and photoionized by the central ionizing continuum. The temporal behaviors of these emission lines are blurred echoes of the continuum variations (i.e., reverberation mapping; RM) and directly reflect structures and kinematics information of BLRs through the so-called transfer function (also known as velocity-delay map). Based on the previous works of Rybicki & Press (1992) and Zu et al. (2011), we develop an extended, non-parametric approach to determine the transfer function for RM data, in which the transfer function is expressed as a sum of a family of relatively-displaced Gaussian response functions. As such, arbitrary shapes of transfer functions associated with complicated BLR geometry can be seamlessly included, enabling us to relax the presumption of a specified transfer function frequently adopted in previous studies and to let it be determined by observation data. We formulate our approach in the previously well-established framework that incorporates the statistical modeling of the continuum variations as a damped random walk process and takes into account the long-term secular variations which are irrelevant to RM signals. Application to reverberation mapping data shows fidelity of our approach.
Very little information has been so far gathered on atomic jets from young embedded low mass sources (class I stars), especially in the inner jet region. We exploit multiwave spectroscopic observations to infer physical conditions of the inner region of HH34 IRS, a prototypical class I jet. We use a deep X-shooter spectrum ( lambda = 350-2300 nm, R between 8000 and 18000) detecting lines with upper energies from 8000 to 31000 cm-1 . Statistical equilibrium and ionization models are adopted to derive the jet main physical parameters. We derive the physical conditions for the extended high velocity jet (HVC, V about -100 km/s) and for the low velocity and compact gas (LVC, V about -20-50 \kms). At the jet base (< 200 AU) the HVC is mostly neutral (x(e)< 0.1) and very dense (n(H) > 5x10^5 cm-3). The LVC gas has the same density and A_V as the HVC, but it is a factor of two colder. Iron abundance in the HVC is close to solar, while it is 3 times subsolar in the LVC, suggesting an origin of the LVC gas from a dusty disk. We infer a jet mass flux rate > 5x10^{-7} Mo/yr. We find that the relationships between accretion luminosity and line luminosity derived for T Tauri stars cannot be directly extended to class I sources surrounded by reflection nebulae since the permitted lines are seen through scattered light. An accretion rate of about 7\,10^{-6} Mo/yr is instead derived from an empirical relationship connecting the mass accretion rate with L([OI])(HVC). In conclusion, the HH34 IRS jet shares many properties with jets from active CTTs stars. It is however denser as a consequence of the larger mass flux rate. These findings suggest that the acceleration and excitation mechanisms in jets are not influenced by evolution and are similar in CTTs and in still embedded and highly accreting sources.
We present observations at 7 mm that fully resolve the two circumstellar disks, and a reanalyses of archival observations at 3.5 cm that resolve along their major axes the two ionized jets, of the class I binary protostellar system L1551 NE. We show that the two circumstellar disks are better fit by a shallow inner and steep outer power-law than a truncated power-law. The two disks have very different transition radii between their inner and outer regions of $\sim$18.6 AU and $\sim$8.9 AU respectively. Assuming that they are intrinsically circular and geometrically thin, we find that the two circumstellar disks are parallel with each other and orthogonal in projection to their respective ionized jets. Furthermore, the two disks are closely aligned if not parallel with their circumbinary disk. Over an interval of $\sim$10 yr, source B (possessing the circumsecondary disk) has moved northwards with respect to and likely away from source A, indicating an orbital motion in the same direction as the rotational motion of their circumbinary disk. All the aforementioned elements therefore share the same axis for their angular momentum, indicating that L1551 NE is a product of rotationally-driven fragmentation of its parental core. Assuming a circular orbit, the relative disk sizes are compatible with theoretical predictions for tidal truncation by a binary system having a mass ratio of $\sim$0.2, in agreement with the reported relative separations of the two protostars from the center of their circumbinary disk. The transition radii of both disks, however, are a factor of $\gtrsim$1.5 smaller than their predicted tidally-truncated radii.
Links to: arXiv, form interface, find, astro-ph, recent, 1608, contact, help (Access key information)
We investigate the contribution of major mergers (mass ratios $>1:5$) to stellar mass growth and morphological transformations around the epoch of peak cosmic star formation ($z\sim2$). We visually classify a complete sample of massive (M $>$ 10$^{10}$M$_{\odot}$) galaxies at this epoch, drawn from the CANDELS survey, into late-type galaxies, major mergers, spheroids and disturbed spheroids which show morphological disturbances. Given recent simulation work, which indicates that recent ($<$0.3-0.4 Gyr) major-merger remnants exhibit clear tidal features in such images, we use the fraction of disturbed spheroids to probe the role of major mergers in driving morphological transformations. The percentage of blue spheroids (i.e. with ongoing star formation) that show morphological disturbances is only 21 $\pm$ 4%, indicating that major mergers are not the dominant mechanism for spheroid creation at $z\sim2$ - other processes, such as minor mergers or cold accretion are likely to be the main drivers of this process. We also use the rest-frame U-band luminosity as a proxy for star formation to show that only a small fraction of the star formation budget ($\sim$3%) is triggered by major mergers. Taken together, our results show that major mergers are not significant drivers of galaxy evolution at $z\sim2$.
We present the results of a recent reverberation mapping campaign for UGC 06728, a nearby low-luminosity Seyfert 1 in a late-type galaxy. Nightly monitoring in the spring of 2015 allowed us to determine an H$\beta$ time delay of $\tau = 1.4 \pm 0.8$ days. Combined with the width of the variable H$\beta$ line profile, we determine a black hole mass of $M_{\rm BH} = (7.1 \pm 4.0) \times 10^5$ M$_{\odot}$. We also constrain the bulge stellar velocity dispersion from higher-resolution long slit spectroscopy along the galaxy minor axis and find $\sigma_{\star} = 51.6 \pm 4.9$ km s$^{-1}$. The measurements presented here are in good agreement with both the $R_{\rm BLR} - L$ relationship and the $M_{\rm BH}-\sigma_{\star}$ relationship for AGNs. Combined with a previously published spin measurement, our mass determination for UGC 06728 makes it the lowest-mass black hole that has been fully characterized, and thus an important object to help anchor the low-mass end of black hole evolutionary models.
Understanding our Galactic Center is easier with insights from nearby galactic nuclei. Both the star formation activity in nuclear gas disks, driven by bars and nuclear bars, and the fueling of low-luminosity AGN, followed by feedback of jets, driving molecular outflows, were certainly present in our Galactic Center, which appears now quenched. Comparisons and diagnostics are reviewed, in particular of m=2 and m=1 modes, lopsidedness, different disk orientations, and fossil evidences of activity and feedback.
We report the discovery of two new candidate stellar systems in the constellation of Cetus using the data from the first two years of the Dark Energy Survey (DES). The objects, DES J0111-1341 and DES J0225+0304, are located at a heliocentric distance of ~ 24 kpc and appear to have old and metal-poor populations. Their distances to the Sagittarius orbital plane, ~ 1.47 kpc (DES J0111-1341) and ~ 0.51 kpc (DES J0225+0304), indicate that they are possibly associated with the Sagittarius dwarf stream. The half-light radius (r_h ~ 4.10 pc) and luminosity (M_V ~ +0.5) of DES J0111-1341 are consistent with it being an ultra-faint stellar cluster, while the half-light radius (r_h ~ 18.70 pc) and luminosity (M_V ~ -1.2) of DES J0225+0304 place it in an ambiguous region of size-luminosity space between stellar clusters and dwarf galaxies. Determinations of the characteristic parameters of the Sagittarius stream, metallicity spread (-2.18 < [Fe/H] < -0.95) and distance gradient (23 kpc < D_sun < 29 kpc), within the DES footprint in the southern hemisphere, using the same DES data, also indicate a possible association between these systems. If these objects are confirmed through spectroscopic follow-up to be gravitationally bound systems and to share a Galactic trajectory with the Sagittarius stream, DES J0111-1341 and DES J0225+0304 would be the first ultra-faint stellar systems associated with the Sagittarius stream. Furthermore, DES J0225+0304 would also be the first confirmed case of an ultra-faint satellite of a satellite.
Stars ejected from the Galactic centre can be used to place important constraints on the Milky Way potential. Since existing hypervelocity stars are too distant to accurately determine orbits, we have conducted a search for nearby candidates using full three-dimensional velocities. Since the efficacy of such studies are often hampered by deficiencies in proper motion catalogs, we have chosen to utilize the reliable, high-precision SDSS Stripe 82 proper motion catalog. Although we do not find any candidates which have velocities in excess of the escape speed, we identify 226 stars on orbits that are consistent with Galactic centre ejection. This number is significantly larger than what we would expect for halo stars on radial orbits and cannot be explained by disk or bulge contamination. If we restrict ourselves to metal-rich stars, we find 29 candidates with [Fe/H] > -0.8 dex and 10 with [Fe/H] > -0.6 dex. Their metallicities are more consistent with what we expect for bulge ejecta, and so we believe these candidates are especially deserving of further study. We have supplemented this sample using our own radial velocities, developing an algorithm to use proper motions for optimizing candidate selection. This technique provides considerable improvement on the blind spectroscopic sample of SDSS, being able to identify candidates with an efficiency around 20 times better than a blind search.
We study the structure, age and metallicity gradients, and dynamical evolution using a cosmological zoom-in simulation of a Milky Way-mass galaxy from the Feedback in Realistic Environments project. In the simulation, stars older than 6 Gyr were formed in a chaotic, bursty mode and have the largest vertical scale heights (1.5-2.5 kpc) by z=0, while stars younger than 6 Gyr were formed in a relatively calm, stable disk. The vertical scale height increases with stellar age at all radii, because (1) stars that formed earlier were thicker "at birth", and (2) stars were kinematically heated to an even thicker distribution after formation. Stars of the same age are thicker in the outer disk than in the inner disk (flaring). These lead to positive vertical age gradients and negative radial age gradients. The radial metallicity gradient is neg- ative at the mid-plane, flattens at larger disk height |Z|, and turns positive above |Z|~1.5kpc. The vertical metallicity gradient is negative at all radii, but is steeper at smaller radii. These trends broadly agree with observations in the Milky Way and can be naturally understood from the age gradients. The vertical stellar density profile can be well-described by two components, with scale heights 200-500 pc and 1-1.5 kpc, respectively. The thick component is a mix of stars older than 4 Gyr which formed through a combination of several mechanisms. Our results also demonstrate that it is possible to form a thin disk in cosmological simulations even with strong stellar feedback.
The GAIA satellite will provide unprecedented phase-space information for our Galaxy and enable a new era of Galactic dynamics. We may soon see successful realizations of Galactoseismology, i.e., inferring the characteristics of the Galactic potential and sub-structure from a dynamical analysis of observed perturbations in the gas or stellar disk of the Milky Way. Here, we argue that to maximally take advantage of the GAIA data and other complementary surveys, it is necessary to build comprehensive models for both the stars and the gas. We outline several key morphological puzzles of the Galactic disk and proposed solutions that may soon be tested.
We examine the processes triggering star formation and Active Galactic Nucleus (AGN) activity in a sample of 25 low redshift ($z<0.13$) gas-rich galaxy mergers observed at milli-arcsecond resolution with Very Long Baseline Interferometry as part of the mJy Imaging VLBA Exploration at 20cm (mJIVE-20) survey. The high ($>10^7$ K) brightness temperature required for an mJIVE-20 detection allows us to unambiguously identify the radio AGN in our sample. We find three such objects. Our VLBI AGN identifications are classified as Seyferts or LINERs in narrow line optical diagnostic plots; mid-infrared colours of our targets and the comparison of H$\alpha$ star formation rates with integrated radio luminosity are also consistent with the VLBI identifications. We reconstruct star formation histories in our galaxies using optical and UV photometry, and find that these radio AGN are not triggered promptly in the merger process, consistent with previous findings for non-VLBI samples of radio AGN. This delay can significantly limit the efficiency of feedback by radio AGN triggered in galaxy mergers. We find that radio AGN hosts have lower star formation rates than non-AGN radio-selected galaxies at the same starburst age. Conventional and VLBI radio imaging shows these AGN to be compact on arcsecond scales. Our modeling suggests that the actual sizes of AGN-inflated radio lobes may be much larger than this, but these are too faint to be detected in existing observations. Deep radio imaging is required to map out the true extent of the AGN, and to determine whether the low star formation rates in radio AGN hosts are a result of the special conditions required for radio jet triggering, or the effect of AGN feedback.
ESO 243-49 is a high-mass (circular velocity $v_{\rm c}\approx200\,{\rm km\,s^{-1}}$) edge-on S0 galaxy in the Abell 2877 cluster at a distance of $\sim95\,{\rm Mpc}$. To elucidate the origin of its thick disc, we use MUSE science verification data to study its kinematics and stellar populations. The thick disc emits $\sim80\%$ of the light at heights in excess of $3.5^{\prime\prime}$ ($1.6\,{\rm kpc}$). The rotation velocities of its stars lag by $30-40\,{\rm km\,s^{-1}}$ compared to those in the thin disc, which is compatible with the asymmetric drift. The thick disc is found to be more metal-poor than the thin disc, but both discs have old ages. We suggest an internal origin for the thick disc stars in high-mass galaxies. We propose that the thick disc formed either ${\rm a)}$ first in a turbulent phase with a high star formation rate and that a thin disc formed shortly afterwards, or ${\rm b)}$ because of the dynamical heating of a thin pre-existing component. Either way, the star formation in ESO 243-49 was quenched just a few Gyrs after the galaxy was born and the formation of a thin and a thick disc must have occurred before the galaxy stopped forming stars. The formation of the discs was so fast that it could be described as a monolithic collapse where several generations of stars formed in a rapid succession.
Single point observations are presented in NH3 (1,1) and (2,2) inversion transitions using the Effelsberg 100 m telescope for a sample of 100 6.7 GHz methanol masers and mapping observations in the 12CO and 13CO (1-0) transitions using the PMO Delingha 13.7 m telescope for 82 sample sources with detected ammonia. A further 62 sources were selected for either 12CO or 13CO line outflow identification, producing 45 outflow candidates, 29 using 12CO and 16 using 13CO data. Twenty-two of the outflow candidates were newly identified, and 23 had trigonometric parallax distances. Physical properties were derived from ammonia lines and CO outflow parameters calculated. Histograms and statistical correlations for ammonia, CO outflow parameters, and 6.7 GHz methanol maser luminosities are also presented. No significant correlation was found between ammonia and maser luminosity. However, weak correlations were found between outflow properties and maser luminosities, which may indicate that outflows are physically associated with 6.7 GHz masers.
Feedback from Active Galactic Nuclei (AGN) and subsequent jet cocoons and outflow bubbles can have a significant impact on star formation in the host galaxy. To investigate feedback physics on small scales, we perform hydrodynamic simulations of realistically fast AGN winds striking Bonnor-Ebert (BE) spheres and examine gravitational collapse and ablation. We test AGN wind velocities ranging from 300--3,000 km s$^{-1}$ and wind densities ranging from 0.5--10 $m_\mathrm{p}\,\mathrm{cm}^{-3}$. We include heating and cooling of low- and high-temperature gas, self-gravity, and spatially correlated perturbations in the shock, with a maximum resolution of 0.01 pc. We find that the ram pressure is the most important factor that determines the fate of the cloud. High ram pressure winds increase fragmentation and decrease the star formation rate, but also cause star formation to occur on a much shorter time scale and with increased velocities of the newly formed stars. We find a threshold ram pressure of $\sim 2\times10^{-8}$ dyne cm$^{-2}$ above which stars are not formed because the resulting clumps have internal velocities large enough to prevent collapse. Our results indicate that simultaneous positive and negative feedback will be possible in a single galaxy as AGN wind parameters will vary with location within a galaxy.
Tidal tails are created in major mergers involving disk galaxies. How the tidal tails trace the assembly history of massive galaxies remains to be explored. We identify a sample of 461 merging galaxies with long tidal tails from 35076 galaxies mass-complete at $M_\star\ge 10^{9.5}\,M_{\odot}$ and $0.2\leq z\leq1$ based on HST/ACS F814W imaging data and public catalogs of the COSMOS field. The long tails refer to these with length equal to or longer than the diameter of their host galaxies. The mergers with tidal tails are selected using our novel $A_{\rm O}-D_{\rm O}$ technique for strong asymmetric features together with visual examination. Our results show that the fraction of tidal-tailed mergers evolves mildly with redshift, as $\sim (1+z)^{2.0\pm0.4}$, and becomes relatively higher in less massive galaxies out to $z=1$. With a timescale of 0.5 Gyr for the tidal-tailed mergers, we obtain that the occurrence rate of such mergers follows $0.01\pm 0.007\,(1+z)^{2.3\pm 1.4}$ Gyr$^{-1}$ and corresponds to $\sim0.3$ events since $z=1$ and roughly one-third of the total budget of major mergers from the literature. For disk-involved major mergers, nearly half of them have undergone a phase with long tidal tails
We present a comparison of SCUBA-2 850-$\mu$m and Herschel 70--500-$\mu$m observations of the L1495 filament in the Taurus Molecular Cloud with the goal of characterising the SCUBA-2 Gould Belt Survey (GBS) data set. We identify and characterise starless cores in three data sets: SCUBA-2 850-$\mu$m, Herschel 250-$\mu$m, and Herschel 250-$\mu$m spatially filtered to mimic the SCUBA-2 data. SCUBA-2 detects only the highest-surface-brightness sources, principally detecting protostellar sources and starless cores embedded in filaments, while Herschel is sensitive to most of the cloud structure, including extended low-surface-brightness emission. Herschel detects considerably more sources than SCUBA-2 even after spatial filtering. We investigate which properties of a starless core detected by Herschel determine its detectability by SCUBA-2, and find that they are the core's temperature and column density (for given dust properties). For similar-temperature cores, such as those seen in L1495, the surface brightnesses of the cores are determined by their column densities, with the highest-column-density cores being detected by SCUBA-2. For roughly spherical geometries, column density corresponds to volume density, and so SCUBA-2 selects the densest cores from a population at a given temperature. This selection effect, which we quantify as a function of distance, makes SCUBA-2 ideal for identifying those cores in Herschel catalogues that are closest to forming stars. Our results can now be used by anyone wishing to use the SCUBA-2 GBS data set.
We suggest a possible mechanism of ultra diffuse galaxy formation: the UDGs may occur as a result of a central collision of galaxies. If the galaxies are young and contain a lot of gas, the collision may kick all the gas off the systems and thus strongly suppress any farther star formation. As a result, the galaxies now have a very low surface brightness and other properties typical of the ultra diffuse galaxies. We use the Coma cluster (where numerous UDGs were recently discovered) to test the efficiency of the process. The mechanism works very well and can transform a significant fraction of the cluster population into ultra diffuse galaxies. The UDGs formed by the process concentrate towards the center of the cluster, and their globular cluster systems remain undamaged, in accordance with observational results. The projected density of UDGs on the cluster images may help us to recognize the mechanism of the UDG formation that realizes in reality.
The mass loss rates of red supergiants (RSGs) govern their evolution towards supernova and dictate the appearance of the resulting explosion. To study how mass-loss rates change with evolution we measure the mass-loss rates (\mdot) and extinctions of 19 red supergiants in the young massive cluster NGC2100 in the Large Magellanic Cloud. By targeting stars in a coeval cluster we can study the mass-loss rate evolution whilst keeping the variables of mass and metallicity fixed. Mass-loss rates were determined by fitting DUSTY models to mid-IR photometry from WISE and Spitzer/IRAC. We find that the \mdot\ in red supergiants increases as the star evolves, and is well described by \mdot\ prescription of de Jager, used widely in stellar evolution calculations. We find the extinction caused by the warm dust is negligible, meaning the warm circumstellar material of the inner wind cannot explain the higher levels of extinction found in the RSGs compared to other cluster stars. We discuss the implications of this work in terms of supernova progenitors and stellar evolution theory. We argue there is little justification for substantially increasing the \mdot\ during the RSG phase, as has been suggested recently in order to explain the absence of high mass Type IIP supernova progenitors. We also argue that an increase in reddening towards the end of the RSG phase, as observed for the two most evolved cluster stars, may provide a solution to the red supergiant problem.
The failure to find evidence for elementary particles that could serve as the constituents of dark matter brings to mind suggestions that dark matter might consist of massive compact objects (MACHOs). In particular, it has recently been argued that MACHOs with masses > 15 solar masses may have been prolifically produced at the onset of the big bang. Although a variety of astrophysical signatures for primordial MACHOs with masses in this range have been discussed in the literature, we favor a strategy that uses the potential for gravitational microlensing of stars outside our galaxy to directly detect the presence of MACHOs in the halo of our galaxy. We point out that the effect of the motion of the Earth on the shape of the microlensing brightening curves provides a promising approach to confirming over the course of next several years that dark matter consists of MACHOs.
Most existing star-galaxy classifiers use the reduced summary information from catalogs, requiring careful feature extraction and selection. The latest advances in machine learning that use deep convolutional neural networks allow a machine to automatically learn the features directly from data, minimizing the need for input from human experts. We present a star-galaxy classification framework that uses deep convolutional neural networks (ConvNets) directly on the reduced, calibrated pixel values. Using data from the Sloan Digital Sky Survey (SDSS) and the Canada-France-Hawaii Telescope Lensing Survey (CFHTLenS), we demonstrate that ConvNets are able to produce accurate and well-calibrated probabilistic classifications that are competitive with conventional machine learning techniques. Future advances in deep learning may bring more success with current and forthcoming photometric surveys, such as the Dark Energy Survey (DES) and the Large Synoptic Survey Telescope (LSST), because deep neural networks require very little, manual feature engineering.
Links to: arXiv, form interface, find, astro-ph, recent, 1608, contact, help (Access key information)
We investigate how accurately phase space distribution functions (DFs) in galactic models can be reconstructed by a made-to-measure (M2M) method, which constructs $N$-particle models of stellar systems from photometric and various kinematic data. The advantage of the M2M method is that this method can be applied to various galactic models without assumption of the spatial symmetries of gravitational potentials adopted in galactic models, and furthermore, numerical calculations of the orbits of the stars cannot be severely constrained by the capacities of computer memories. The M2M method has been applied to various galactic models. However, the degree of accuracy for the recovery of DFs derived by the M2M method in galactic models has never been investigated carefully. Therefore, we show the degree of accuracy for the recovery of the DFs for the anisotropic Plummer model and the axisymmetric St\"{a}ckel model, which have analytic solutions of the DFs. Furthermore, this study provides the dependence of the degree of accuracy for the recovery of the DFs on various parameters and a procedure adopted in this paper. As a result, we find that the degree of accuracy for the recovery of the DFs derived by the M2M method for the spherical target model is a few percent, and more than ten percent for the axisymmetric target model.
We use the hydrodynamic, cosmological EAGLE simulations to investigate how hot gas in haloes condenses to form and grow galaxies. We select haloes from the simulations that are actively cooling and study the temperature, distribution, and metallicity of their hot, cold, and transitioning `cooling' gas, placing these in context of semi-analytic models. Our selection criteria lead us to focus on Milky Way-like haloes. We find the hot-gas density profiles of the haloes form a progressively stronger core over time, the nature of which can be captured by a beta profile that has a simple dependence on redshift. In contrast, the hot gas that actually cools is broadly consistent with a singular isothermal sphere. We find that cooling gas carries a few times the specific angular momentum of the halo and is offset in spin direction from the rest of the hot gas. The gas loses ~60% of its specific angular momentum during the cooling process, generally remaining greater than that of the halo, and is better aligned with the cold gas already in the disc than anything else. Angular-momentum losses are slightly larger when cooling onto dispersion-supported galaxies. We show that an exponential surface density profile for gas arriving on a disc remains a reasonable approximation, but a cusp is always present, and disc scale radii are larger than predicted by a vanilla Fall & Efstathiou model. These scale radii are still closely correlated with the halo spin parameter, for which we suggest an updated prescription for galaxy formation models.
We model the cosmic distributions in space and time of the formation sites of the first stars that may be the progenitors of supermassive black holes (SMBHs). Pop III.1 stars are defined to form in dark matter minihalos (i.e., with masses $\sim10^6\:M_\odot$) that are isolated from neighboring astrophysical sources by a given isolation distance, $d_{\rm iso}$. We assume these sources are the seeds for the cosmic population of SMBHs, based on a model of protostellar support by dark matter annihilation heating that allows these objects to accrete most of the baryonic content of their minihalos, i.e., $\gtrsim10^5\:M_\odot$. Exploring a range of $d_{\rm iso}$ from 10 to 100~kpc (proper distances), we predict the evolution with redshift of the number density of these Pop III.1 sources and their SMBH remnants. In the context of this model, the local, $z=0$ density of SMBHs constrains $d_{\rm iso}\gtrsim100$~kpc (i.e., a comoving distance of 3~Mpc at $z\simeq30$). In our simulated ($\sim$40.96 $h^{-1}$~Mpc)$^3$ comoving volume, Pop III.1 stars start forming just after $z=40$. Their formation is largely complete by $z\simeq25$ to 20 for $d_{\rm iso}=100$ to 50~kpc. We follow the evolution of these sources down to $z=10$, by which point the SMBHs are expected to reside in halos with $\gtrsim10^8\:M_\odot$. Over this period, there is relatively limited merging of SMBHs for these values of $d_{\rm iso}$. We also predict the clustering properties of the SMBHs in this model at $z=10$: the effective feedback suppression of neighboring sources in the $d_{\rm iso}=100$ model leads to a relatively flat angular correlation function. The general results of this work will be used for population synthesis models of active galactic nuclei, which can be compared with future high redshift observations.
Radiation feedback from stellar clusters is expected to play a key role in setting the rate and efficiency of star formation in giant molecular clouds (GMCs). To investigate how radiation forces influence realistic turbulent systems, we have conducted a series of numerical simulations employing the {\it Hyperion} radiation hydrodynamics solver, considering the regime that is optically thick to ultraviolet (UV) and optically thin to infrared (IR) radiation. Our model clouds cover initial surface densities between $\Sigma_{\rm cl,0} \sim 10-300~M_{\odot}~{\rm pc^{-2}}$, with varying initial turbulence. We follow them through turbulent, self-gravitating collapse, formation of star clusters, and cloud dispersal by stellar radiation. All our models display a lognormal distribution of gas surface density $\Sigma$; for an initial virial parameter $\alpha_{\rm vir,0} = 2$, the lognormal standard deviation is $\sigma_{\rm ln \Sigma} = 1-1.5$ and the star formation rate coefficient $\varepsilon_{\rm ff,\bar\rho} = 0.3-0.5$, both of which are sensitive to turbulence but not radiation feedback. The net star formation efficiency $\varepsilon_\mathrm{final}$ increases with $\Sigma_{\rm cl,0}$ and decreases with $\alpha_{\rm vir,0}$. We interpret these results via a simple conceptual framework, whereby steady star formation increases the radiation force, such that local gas patches at successively higher $\Sigma$ become unbound. Based on this formalism (with fixed $\sigma_{\rm ln \Sigma}$), we provide an analytic upper bound on $\varepsilon_\mathrm{final}$, which is in good agreement with our numerical results. The final star formation efficiency depends on the distribution of Eddington ratios in the cloud and is strongly increased by turbulent compression of gas.
We combine a number of recent studies of the extended main sequence turn-off (eMSTO) phenomenon in intermediate age stellar ($1-2$ Gyr) clusters in the Large Magellanic Cloud (LMC) in order to investigate its origin. By employing the largest sample of eMSTO LMC clusters so far used, we show that cluster core radii, masses, and dynamical state are not related to the genesis of eMSTOs. Indeed, clusters in our sample have core radii, masses and age-relaxation time ratios in the range $\approx$ 2--6 pc, 3.35- 5.50 (log($M_{cls}$/$M_\odot$) and 0.2-8.0, respectively. These results imply that the eMSTO phenomenon is not caused by actual age spreads within the clusters. Furthermore, we confirm from a larger cluster sample recent results including young eMSTO LMC clusters, that the FWHM at the MSTOs correlates most strongly with cluster age, suggesting that a stellar evolutionary effect is the underlying cause.
We discuss 6 GHz JVLA observations covering a volume-limited sample of 178 low redshift ($0.2 < z < 0.3$) optically selected QSOs. Our 176 radio detections fall into two clear categories: (1) About $20$\% are radio-loud QSOs (RLQs) having spectral luminosities $L_6 \gtrsim 10^{\,23.2} \mathrm{~W~Hz}^{-1}$ primarily generated in the active galactic nucleus (AGN) responsible for the excess optical luminosity that defines a \emph{bona fide} QSO. (2) The radio-quiet QSOs (RQQs) have $10^{\,21} \lesssim L_6 \lesssim 10^{\,23.2} \mathrm{~W~Hz}^{-1}$ and radio sizes $\lesssim 10 \mathrm{~kpc}$, and we suggest that the bulk of their radio emission is powered by star formation in their host galaxies. "Radio silent" QSOs ($L_6 \lesssim 10^{\,21} \mathrm{~W~Hz}^{-1}$) are rare, so most RQQ host galaxies form stars faster than the Milky Way; they are not "red and dead" ellipticals. Earlier radio observations did not have the luminosity sensitivity $L_6 \lesssim 10^{\,21} \mathrm{~W~Hz}^{-1}$ needed to distinguish between such RLQs and RQQs. Strong, generally double-sided, radio emission spanning $\gg 10 \mathrm{~kpc}$ was found associated with 13 of the 18 RLQ cores having peak flux densities $S_\mathrm{p} > 5 \mathrm{~mJy~beam}^{-1}$ ($log(L) \gtrsim 24$). The radio luminosity function of optically selected QSOs and the extended radio emission associated with RLQs are both inconsistent with simple "unified" models that invoke relativistic beaming from randomly oriented QSOs to explain the difference between RLQs and RQQs. Some intrinsic property of the AGNs or their host galaxies must also determine whether or not a QSO appears radio loud.
Gas accretion and radial flows are key ingredients of the chemical evolution of spiral galaxies. They are also tightly linked to each other (accretion drives radial flows, due to angular momentum conservation) and should therefore be modelled simultaneously. We summarise an algorithm that can be used to consistently compute accretion profiles, radial flows and abundance gradients under quite general conditions and we describe illustrative applications to the Milky Way. We find that gas-phase abundance gradients strongly depend on the angular momentum of the accreting material and, in the outer regions, they are significantly affected by the choice of boundary conditions.
We analysed the optical spectra of HII regions extracted from a sample of 350 galaxies of the CALIFA survey. We calculated total O/H abundances and, for the first time, N/O ratios using the semi-empirical routine HII-CHI-mistry, which, according to P\'erez-Montero (2014), is consistent with the direct method and reduces the uncertainty in the O/H derivation using [NII] lines owing to the dispersion in the O/H-N/O relation. Then we performed linear fittings to the abundances as a function of the de-projected galactocentric distances. The analysis of the radial distribution both for O/H and N/O in the non-interacting galaxies reveals that both average slopes are negative, but a non-negligible fraction of objects have a flat or even a positive gradient (at least 10\% for O/H and 4\% for N/O). The slopes normalised to the effective radius appear to have a slight dependence on the total stellar mass and the morphological type, as late low-mass objects tend to have flatter slopes. No clear relation is found, however, to explain the presence of inverted gradients in this sample, and there is no dependence between the average slopes and the presence of a bar. The relation between the resulting O/H and N/O linear fittings at the effective radius is much tighter (correlation coefficient $\rho_s$ = 0.80) than between O/H and N/O slopes ($\rho_s$ = 0.39) or for O/H and N/O in the individual \hii\ regions ($\rho_s$ = 0.37). These O/H and N/O values at the effective radius also correlate very tightly (less than 0.03 dex of dispersion) with total luminosity and stellar mass. The relation with other integrated properties, such as star formation rate, colour, or morphology, can be understood only in light of the found relation with mass.
Placing bright sub-millimetre galaxies (SMGs) within the broader context of galaxy formation and evolution requires accurate measurements of their clustering, which can constrain the masses of their host dark matter halos. Recent work has shown that the clustering measurements of these galaxies may be affected by a `blending bias,' which results in the angular correlation function of the sources extracted from single-dish imaging surveys being boosted relative to that of the underlying galaxies. This is due to confusion introduced by the coarse angular resolution of the single-dish telescope and could lead to the inferred halo masses being significantly overestimated. We investigate the extent to which this bias affects the measurement of the correlation function of SMGs when it is derived via a cross-correlation with a more abundant galaxy population. We find that the blending bias is essentially the same as in the auto-correlation case and conclude that the best way to reduce its effects is to calculate the angular correlation function using SMGs in narrow redshift bins. Blending bias causes the inferred host halo masses of the SMGs to be overestimated by a factor of $\sim6$ when a redshift interval of $\delta z=3$ is used. However, this reduces to a factor of $\sim2$ for $\delta z=0.5$. The broadening of photometric redshift probability distributions with increasing redshift can therefore impart a mild halo `downsizing' effect onto the inferred host halo masses, though this trend is not as strong as seen in recent observational studies.
Studying the Milky Way disk structure using stars in narrow bins of [Fe/H] and [alpha/Fe] has recently been proposed as a powerful method to understand the Galactic thick and thin disk formation. It has been assumed so far that these mono-abundance populations (MAPs) are also coeval, or mono-age, populations. Here we study this relationship for a Milky Way chemo-dynamical model and show that equivalence between MAPs and mono-age populations exists only for the high-[alpha/Fe] tail, where the chemical evolution curves of different Galactic radii are far apart. At lower [alpha/Fe]-values a MAP is composed of stars with a range in ages, even for small observational uncertainties and a small MAP bin size. Due to the disk inside-out formation, for these MAPs younger stars are typically located at larger radii, which results in negative radial age gradients that can be as large as 2 Gyr/kpc. Positive radial age gradients can result for MAPs at the lowest [alpha/Fe] and highest [Fe/H] end. Such variations with age prevent the simple interpretation of observations for which accurate ages are not available. Studying the variation with radius of the stellar surface density and scale-height in our model, we find good agreement to recent analyses of the APOGEE red-clump (RC) sample when 1-4 Gyr old stars dominate (as expected for the RC). Our results suggest that the APOGEE data are consistent with a Milky Way model for which mono-age populations flare for all ages. We propose observational tests for the validity of our predictions and argue that using accurate age measurements, such as from asteroseismology, is crucial for putting constraints on the Galactic formation and evolution.
Galaxy cluster mass determinations achieved using X-ray and Sunyaev-Zeldovich data combined with the assumption of hydrostatic equilibrium are generally biased. The bias exists for two main reasons: non-thermal pressure forces are expected to contribute to the overall pressure balance and deviations from spherical symmetry and hydrostatic equilibrium can be present. In this paper, we use a sample of zoom-in hydrodynamical simulations of galaxy clusters to measure the magnitude of hydrostatic bias and the contribution from turbulence to the total pressure. We propose a new empirical model for turbulent pressure based on our simulations that can be applied to observations. We show that our model can be successfully applied to remove most of the bias related to neglection of turbulent pressure, which is usually not included in hydrostatic cluster mass profile reconstructions. The use of this model may significantly improve the calibration of cluster scaling relations that are a key tool for cluster cosmology.
Links to: arXiv, form interface, find, astro-ph, recent, 1608, contact, help (Access key information)
The central region of the Virgo cluster of galaxies contains thousands of globular clusters (GCs), an order of magnitude more than the numbers found in the Local Group. Relics of early star formation epochs in the universe, these GCs also provide ideal targets to test our understanding of the Spectral Energy Distributions (SEDs) of old stellar populations. Based on photometric data from the Next Generation Virgo cluster Survey (NGVS) and its near-infrared counterpart NGVS-IR, we select a robust sample of 1846 GCs with excellent photometry and spanning the full range of colors present in the Virgo core. The selection exploits the well defined locus of GCs in the uiK diagram and the fact that the globular clusters are marginally resolved in the images. We show that the GCs define a narrow sequence in 5-dimensional color space, with limited but real dispersion around the mean sequence. The comparison of these SEDs with the predictions of eleven widely used population synthesis models highlights differences between models, and also shows that no single model adequately matches the data in all colors. We discuss possible causes for some of these discrepancies. Forthcoming papers of this series will examine how best to estimate photometric metallicities in this context, and compare the Virgo globular cluster colors with those in other environments.
We present N-body and hydrodynamical simulations of the response of the Milky Way's baryonic disc to the presence of the Large Magellanic Cloud during a first infall scenario. For a fiducial galactic model reproducing the gross properties of the Galaxy, we explore a set of six initial conditions for the LMC of varying mass which all evolve to fit the measured constraints on its current position and velocity with respect to the Galactic Center. We find that the LMC can produce strong disturbances - warping of the stellar and gaseous discs - in the Galaxy, without violating constraints from the phase-space distribution of stars in the Solar Neighbourhood. All models correctly reproduce the phases of the warp and its anti-symmetrical shape about the disc's mid-plane. If the warp is due to the LMC alone, then the largest mass model is favoured ($2.5\times10^{11}\,\rm{M_{\odot}}$). Still, some quantitative discrepancies remain, including deficits in height of $\Delta Z=0.7 \, \rm{kpc}$ at $R=22 \,\rm{kpc}$ and $\Delta Z=0.7\, \rm{kpc}$ at $R=16 \,\rm{kpc}$. This suggests that even higher infall masses for the LMC's halo are allowed by the data. A comparison with the vertical perturbations induced by a heavy Sagittarius dSph model ($10^{11}\,\rm{M_{\odot}}$) suggest that positive interference with the LMC warp is expected at $R=16 \, \rm{kpc}$. We conclude that the vertical structure of the Galactic disc beyond the Solar Neighbourhood may jointly be shaped by its most massive satellites. As such, the current structure of the MW suggests we are seeing the process of disc heating by satellite interactions in action.
We present mass and mass profile estimates for the Milky Way Galaxy using the Bayesian analysis developed by Eadie et al (2015b) and using globular clusters (GCs) as tracers of the Galactic potential. The dark matter and GCs are assumed to follow different spatial distributions; we assume power-law model profiles and use the model distribution functions described in Evans et al. (1997); Deason et al (2011, 2012a). We explore the relationships between assumptions about model parameters and how these assumptions affect mass profile estimates. We also explore how using subsamples of the GC population beyond certain radii affect mass estimates. After exploring the posterior distributions of different parameter assumption scenarios, we conclude that a conservative estimate of the Galaxy's mass within 125kpc is $5.22\times10^{11} M_{\odot}$, with a $50\%$ probability region of $(4.79, 5.63) \times10^{11} M_{\odot}$. Extrapolating out to the virial radius, we obtain a virial mass for the Milky Way of $6.82\times10^{11} M_{\odot}$ with $50\%$ credible region of $(6.06, 7.53) \times 10^{11} M_{\odot}$ ($r_{vir}=185^{+7}_{-7}$kpc). If we consider only the GCs beyond 10kpc, then the virial mass is $9.02~(5.69, 10.86) \times 10^{11} M_{\odot}$ ($r_{vir}=198^{+19}_{-24}$kpc). We also arrive at an estimate of the velocity anisotropy parameter $\beta$ of the GC population, which is $\beta=0.28$ with a $50\%$ credible region (0.21, 0.35). Interestingly, the mass estimates are sensitive to both the dark matter halo potential and visible matter tracer parameters, but are not very sensitive to the anisotropy parameter.
Reionisation in the early Universe is likely driven by dwarf galaxies. Using cosmological, zoom-in, radiation-hydrodynamic simulations, we study the escape of Lyman continuum (LyC) photons from mini-haloes with $M_{\rm halo} \le 10^8\,M_\odot$. Our simulations include a new thermo-turbulent star formation model, non-equilibrium chemistry, and relevant stellar feedback processes (photoionisation by young massive stars, radiation pressure, and mechanical supernova explosions). We find that the photon number-weighted mean escape fraction in mini-haloes is higher ($\sim20$-$40\%$) than that in atomic-cooling haloes, although the instantaneous fraction in individual haloes varies significantly. The escape fraction from Pop III stars is found to be significant ($\ge10\%$) only when the mass is greater than $\sim$100\,\msun. Because star formation is stochastic and dominated by a few gas clumps, the escape fraction is generally determined by radiation feedback (heating due to photo-ionisation), rather than supernova explosions. We find that the resulting stellar mass of the proto-galaxies in mini-haloes follows the slope and normalisation reported in Kimm \& Cen, which is similar to the empirical stellar mass-to-halo mass relation derived in the local Universe. Based on simple analytic calculations, we show that LyC photons from mini-haloes are, despite their high escape fractions, of minor importance for reionisation, as feedback reduces star formation very efficiently in mini-haloes. We confirm previous claims that stars in atomic-cooling haloes with masses $10^8\,M_\odot\le M_{\rm halo} \le 10^{11}\,M_\odot$ are likely to be the most important source of reionisation.
We report idealized simulations that mimic the growth of galaxy disks embedded in responsive halos and bulges. The disks manifested an almost overwhelming tendency to form strong bars that we found very difficult to prevent. We found that fresh bars formed in growing disks after we had destroyed the original, indicating that bar formation also afflicts continued galaxy evolution, and not just the early stages of disk formation. This behavior raises still more insistently the previously unsolved question of how some galaxies avoid bars. Since our simulations included only collisionless star and halo particles, our findings may apply to gas-poor galaxies only; however the conundrum persists for the substantial unbarred fraction of those galaxies. Our original objective was to study how internal dynamics rearranged the distribution of mass in the disk as a generalization of our earlier study with rigid spherical components. With difficulty, we were able to construct some models that were not strongly influenced by bars, and found that halo compression, and angular momentum exchange with the disk did not alter our earlier conclusion that spiral activity is largely responsible for creating smooth density profiles and rotation curves.
We develop a model of dust evolution in a multiphase, inhomogeneous ISM including dust growth and destruction processes. The physical conditions for grain evolution are taken from hydrodynamical simulations of giant molecular clouds in a Milky Way-like spiral galaxy. We improve the treatment of dust growth by accretion in the ISM to investigate the role of the temperature-dependent sticking coefficient and ion-grain interactions. From detailed observational data on the gas-phase Si abundances [Si/H]_{gas} measured in the local Galaxy, we derive a relation between the average [Si/H]_{gas} and the local gas density n(H) which we use as a critical constraint for the models. This relation requires a sticking coefficient that decreases with the gas temperature. The synthetic relation constructed from the spatial dust distribution reproduces the slope of -0.5 of the observed relation in cold clouds. This slope is steeper than that for the warm medium and is explained by the dust growth. We find that it occurs for all adopted values of the minimum grain size a_{min} from 1 to 5nm. For the classical cut-off of a_{min}=5 nm, the ion-grain interactions result in longer growth timescales and higher [Si/H]_{gas} than the observed values. For a_{min} below 3 nm, the ion-grain interactions enhance the growth rates, steepen the slope of [Si/H]_{gas}-n(H) relation and provide a better match to observations. The rates of dust re-formation in the ISM by far exceed the rates of dust production by stellar sources as expected from simple evolution models. After the cycle of matter in and out of dust reaches a steady state, the dust growth balances the destruction operating on similar timescales of 350 Myr.
We have analysed the chemical and kinematic properties of the 20 and 50 km s$^{-1}$ molecular clouds in the Central Molecular Zone of the Milky Way Galaxy, as well as those of the molecular ridge bridging these two clouds. Our work has utilized 37 molecular transitions in the 0.65, 3 and 7-mm wavebands, from the Mopra and NANTEN2 telescopes. The 0.65-mm NANTEN2 data highlights a dense condensation of emission within the western part of the 20 km s$^{-1}$ cloud, visible in only four other transitions, which are 3-mm H$^{13}$CN (1--0), H$^{13}$CO$^{+}$ (1--0), HNC (1--0) and N$_{2}$H$^{+}$ (1--0), suggesting that the condensation is moderately optically thick and cold. We find that while the relative chemical abundances between both clouds are alike in many transitions, suggesting little variation in the chemistry between both clouds; the 20 km s$^{-1}$, cold cloud is brighter than the 50 km s$^{-1}$ cloud in shock and high density tracers. The spatial distribution of enhanced emission is widespread in the 20 km s$^{-1}$ cloud, as shown via line ratio maps. The position velocity diagrams across both clouds indicate that the gas is well mixed. We show that the molecular ridge is most likely part of the 20 km s$^{-1}$ cloud and that both of them may possibly extend to include the 50 km s$^{-1}$ cloud, as part of one larger cloud. Furthermore, we expect that the 20 km s$^{-1}$ cloud is being tidally sheared as a result of the gravitational potential from Sgr A*.
Pegasus III (Peg III) is one of the few known ultra-faint stellar systems in the outer halo of the Milky Way. We present the results from a follow-up campaign with Magellan/IMACS and Keck/DEIMOS. Deep stellar photometry down to $r_0\approx 26$ mag has allowed accurate measurements of its photometric and structural properties. The color-magnitude diagram of Peg III confirms that the stellar system is well described by an old ($\sim$13.5 Gyr), metal-poor ($\langle\lbrack$Fe/H$\rbrack\rangle\sim -2.5$ dex) and $\alpha$-enhanced ([$\alpha$/Fe]$\sim+0.4$ dex.) stellar population at a heliocentric distance of $215\pm12$ kpc. The revised half-light radius $r_{h}=53\pm14$ pc, ellipticity $\epsilon=0.38^{+0.22}_{-0.38}$, and total luminosity $M_{V}=-3.4\pm0.4$ are in good agreement with the values quoted in our previous paper. We further report on the spectroscopic identification of seven, possibly eight member stars of Peg III. Peg III has a systemic velocity of $-222.9 \pm 2.6$ km s$^{-1}$ and a velocity dispersion of $5.4^{+3.0}_{-2.5}$ km s$^{-1}$. The inferred dynamical mass within the half-light radius is $1.4^{+3.0}_{-1.1} \times 10^6\rm{M_{\odot}}$ and the mass-to-light ratio $\rm{M/L}_{V} = 1470^{+5660}_{-1240}$ $\rm{M_{\odot}/L_{\odot}}$, providing further evidence that Peg III is a dwarf galaxy satellite. We find that Peg III and another distant dwarf satellite Pisces II lie relatively close to each other ($\Delta d_{spatial}\approx47$ kpc) and share similar radial velocities in the Galactic standard-of-rest frame ($\Delta v_{GSR}\approx10$ km$s^{-1}$). This suggests that they may be a physically bound pair.
Using combined asteroseismic and spectroscopic observations of 418 red-giant stars close to the Galactic disc plane (6 kpc $<R_{\rm Gal}\lesssim13$ kpc, $|Z_{\rm Gal}|<0.3$ kpc), we measure the age dependence of the radial metallicity distribution in the Milky Way's thin disc over cosmic time. The radial metallicity gradient is constant, at a value of $\sim-0.07$ dex/kpc, for stellar populations younger than $\sim4$ Gyr, and flattens to reach a value of $\sim-0.02$ dex/kpc for stars with ages between 6 and 10 Gyr. After a long-lasting observational debate, our data are able to rule out certain chemical-evolution scenarios, especially models in which the radial abundance gradient flattens with time. Our results are in excellent agreement with a state-of-the-art chemo-dynamical Milky-Way model in which the evolution of the abundance gradient and its scatter can be entirely explained by a non-varying negative metallicity gradient in the interstellar medium, together with stellar radial mixing and migration. We also offer an explanation for why old open clusters in the Solar Neighbourhood can be more metal rich than young clusters: Already within 2 Gyr, radial migration can bring metal-rich clusters from the innermost regions of the disc to Galactocentric radii of 5 to 8 kpc. In the near future, asteroseismic data from the K2 mission will allow for improved statistics and a better coverage of the inner-disc regions, thereby providing tighter constraints on the evolution of the central parts of the Milky Way.
We present the first detection of a jet in the far-IR [O I] lines from an intermediate mass protostar. We have carried out a Herschel/PACS spectral mapping study in the [O I] lines of OMC-2 FIR 3 and FIR 4, two of the most luminous protostars in Orion outside of the Orion Nebula. The spatial morphology of the fine structure line emission reveals the presence of an extended photodissociation region (PDR) and a narrow, but intense jet connecting the two protostars. The jet seen in [O I] emission is spatially aligned with the Spitzer/IRAC 4.5 micron jet and the CO (6-5) molecular outflow centered on FIR 3. The mass loss rate derived from the total [O I] 63 micron line luminosity of the jet is 7.7 x 10^-6 M_sun/yr, more than an order of magnitude higher than that measured for typical low mass class 0 protostars. The implied accretion luminosity is significantly higher than the observed bolometric luminosity of FIR 4, indicating that the [O I] jet is unlikely to be associated with FIR 4. We argue that the peak line emission seen toward FIR 4 originates in the terminal shock produced by the jet driven by FIR 3. The higher mass-loss rate that we find for FIR 3 is consistent with the idea that intermediate mass protostars drive more powerful jets than their low-mass counterparts. Our results also call into question the nature of FIR 4.
We investigate the evolution of far-IR CO emission from protostars observed with Herschel/PACS for 50 sources from the combined sample of HOPS and DIGIT Herschel key programs. From the uniformly sampled spectral energy distributions, we computed $L_{\rm{bol}}$, $T_{\rm{bol}}$ and $L_{\rm {bol}}/L_{\rm {smm}}$ for these sources to search for correlations between far-IR CO emission and protostellar properties. We find a strong and tight correlation between far-IR CO luminosity ($L^{\rm fir}_{\rm CO}$) and the bolometric luminosity ($L_{\rm{bol}}$) of the protostars with $L^{\rm fir}_{\rm CO}$ $\propto$ $L_{\rm{bol}}^{0.7}$. We, however, do not find a strong correlation between $L^{\rm fir}_{\rm CO}$ and protostellar evolutionary indicators, $T_{\rm{bol}}$ and $L_{\rm {bol}}/L_{\rm {smm}}$. FIR CO emission from protostars traces the currently shocked gas by jets/outflows, and $L^{\rm fir}_{\rm CO}$ is proportional to the instantaneous mass loss rate, $\dot{M}_{\rm{out}}$. The correlation between $L^{\rm fir}_{\rm CO}$ and $L_{\rm{bol}}$ is indicative of instantaneous $\dot{M}_{\rm{out}}$ tracking instantaneous $\dot{M}_{\rm{acc}}$. The lack of correlation between $L^{\rm fir}_{\rm CO}$ and evolutionary indicators $T_{\rm{bol}}$ and $L_{\rm {bol}}/L_{\rm {smm}}$ suggests that $\dot{M}_{\rm{out}}$ and, therefore, $\dot{M}_{\rm{acc}}$ do not show any clear evolutionary trend. These results are consistent with mass accretion/ejection in protostars being episodic. Taken together with the previous finding that the time-averaged mass ejection/accretion rate declines during the protostellar phase (e.g., Bontemps et al. 1996), our results suggest that the instantaneous accretion/ejection rate of protostars is highly time variable and episodic, but the amplitude and/or frequency of this variability decreases with time such that the time averaged accretion/ejection rate declines with system age.
We present and discuss initial selection criteria and first results in M33 from a systematic search for extragalactic symbiotic stars. We show that the presence of diffuse interstellar gas emission can significantly contaminate the spectra of symbiotic star candidates. This important effect forces upon us a more stringent working definition of an extragalactic symbiotic star. We report the first detections and spectroscopic characterisation of 12 symbiotic binaries in M33. We found that four of our systems contain carbon-rich giants. In another two of them the giant seems to be a Zr-enhanced MS star, while the remaining six objects host M-type giants. The high number ratio of C to M giants in these binaries is consistent with the low metallicity of M33. The spatial and radial velocity distributions of these new symbiotic binaries are consistent with a wide range of progenitor star ages.
Fast jets are thought to be a crucial ingredient of star formation because
they might extract angular momentum from the disk and thus allow mass accretion
onto the star. However, it is unclear whether jets are ubiquitous, and
likewise, their contribution to mass and angular momentum extraction during
protostar formation remains an open question.
Our aim is to investigate the ejection process in the low-mass Class 0
protostar L1157. This source is associated with a spectacular bipolar outflow,
and the recent detection of high-velocity SiO suggests the occurrence of a jet.
Observations of CO 2-1 and SiO 5-4 at 0.8" resolution were obtained with the
IRAM Plateau de Bure Interferometer as part of the CALYPSO large program. The
jet and outflow structure were fit with a precession model. We derived the
column density of CO and SiO, as well as the jet mass-loss rate and mechanical
luminosity.
High-velocity CO and SiO emission resolve for the first time the first 200 au
of the outflow-driving molecular jet. The jet is strongly asymmetric, with the
blue lobe 0.65 times slower than the red lobe. This suggests that the
large-scale asymmetry of the outflow is directly linked to the jet velocity and
that the asymmetry in the launching mechanism has been at work for the past
1800 yr. Velocity asymmetries are common in T Tauri stars, which suggests that
the jet formation mechanism from Class 0 to Class II stages might be similar.
Our model simultaneously fits the inner jet and the clumpy 0.2 pc scale outflow
by assuming that the jet precesses counter-clockwise on a cone inclined by 73
degree to the line of sight with an opening angle of 8 degree on a period of
1640 yr. The estimated jet mass flux and mechanical luminosity are 7.7e-7
Msun/yr, and 0.9 Lsun, indicating that the jet could extract at least 25% of
the gravitational energy released by the forming star.
The reported observations of an unidentified X-ray line feature at $\sim$3.5 keV have driven a lively discussion about its possible dark matter origin. Motivated by this, we have measured the \emph{K}-shell X-ray spectra of highly ionized bare sulfur ions following charge exchange with gaseous molecules in an electron beam ion trap, as a source of or a contributor to this X-ray line. We produce $\mathrm{S}^{16+}$ and $\mathrm{S}^{15+}$ ions and let them capture electrons in collision with those molecules with the electron beam turned off while recording X-ray spectra. We observed a charge-exchanged-induced X-ray feature at the Lyman series limit (3.47 $\pm$ 0.06 keV). The inferred X-ray energy is in full agreement with the reported astrophysical observations and supports the novel scenario proposed by Gu and Kaastra (A \& A \textbf{584}, {L11} (2015)).
Deep observations of galaxy outskirts reveal faint extended stellar components (ESCs) of streams, shells, and halos, which are ghostly remnants of the tidal disruption of satellite galaxies. We use cosmological galaxy formation simulations in Cold Dark Matter (CDM) and Warm Dark Matter (WDM) models to explore how the dark matter model influences the spatial, kinematic, and orbital properties of ESCs. These reveal that the spherically averaged stellar mass density at large galacto-centric radius can be depressed by up to a factor of 10 in WDM models relative to the CDM model, reflecting the anticipated suppressed abundance of satellite galaxies in WDM models. However, these differences are much smaller in WDM models that are compatible with observational limits, and are comparable in size to the system-to-system variation we find within the CDM model. This suggests that it will be challenging to place limits on dark matter using only the unresolved ESC.
The putative black holes which may constitute all the dark matter are described by a Kerr metric with only two parameters, mass M and angular momentum J. There has been little discussion of J since it plays no role in the upcoming attempt at detection by microlensing. Nevertheless J does play a central role in understanding the previous lack of detection, especially of CMB distortion. We explain why bounds previously derived from lack of CMB distortion are too strong for primordial black holes with J non-vanishing. Almost none of the dark matter black holes can be from stellar collapse, and nearly all are primordial, to avoid excessive CMB distortion.
Links to: arXiv, form interface, find, astro-ph, recent, 1608, contact, help (Access key information)
We present the new single dish CO (3-2) emission data obtained toward 19 early stage and 7 late stage nearby merging galaxies using the Atacama Submillimeter Telescope Experiment (ASTE). Combining with the single dish and interferometric data of galaxies observed in previous studies, we investigate the relation between the CO (3-2) luminosity (L'CO(3-2)) and the far Infrared luminosity (LFIR) in a sample of 29 early stage and 31 late stage merging galaxies, and 28 nearby isolated spiral galaxies. We find that normal isolated spiral galaxies and merging galaxies have different slopes (alpha) in the log L'CO(3-2) - log LFIR plane (alpha ~ 0.79 for spirals and ~ 1.12 for mergers). The large slope (alpha > 1) for merging galaxies can be interpreted as an evidence for increasing Star Formation Efficiency (SFE=LFIR/L'CO(3-2)) as a function of LFIR. Comparing our results with sub-kpc scale local star formation and global star-burst activity in the high-z Universe, we find deviations from the linear relationship in the log L'CO(3-2) - log LFIR plane for the late stage mergers and high-z star forming galaxies. Finally, we find that the average SFE gradually increases from isolated galaxies, merging galaxies, and to high-z submillimeter galaxies / quasi-stellar objects (SMGs/QSOs). By comparing our findings with the results from numerical simulations, we suggest; (1) inefficient star-bursts triggered by disk-wide dense clumps occur in the early stage of interaction and (2) efficient star-bursts triggered by central concentration of gas occur in the final stage. A systematic high spatial resolution survey of diffuse and dense gas tracers is a key to confirm this scenario.
We studied AGN activity in twelve LSSs in the ORELSE survey, at 0.65<z<1.28, using a combination of Chandra observations, optical and NIR imaging and spectroscopy. We located a total of 61 AGNs across our sample that were successfully matched to optical counterparts in the LSSs. Seeking to study AGN triggering mechanisms, we examined the spatial distribution of the AGNs and their average spectral properties. We found that AGN populations across our sample had less time since the last starburst than the overall galaxy populations. We did not find any relation between AGN activity and location within the LSSs, suggesting triggering mechanisms which depend on global environment are at most sub-dominant. To focus on differences between our AGNs, we grouped them into four sub-samples based on the spectral properties of their parents LSSs. We found one of the sub-samples, SG0023 & SC1604, stood out from the others: AGNs in this sample were disproportionately luminous, their average time since the last starburst event was the smallest, despite the fact that this was not true of the overall galaxy population in those LSSs, and both the AGNs and the overall galaxy population had the largest fraction of close kinematic pairs, which indicates a higher rate of galaxy mergers and interactions. These results suggest that major mergers are driving AGN activity in SG0023 & SC1604, while other processes are likely triggering less luminous AGNs in the rest of our sample. Additionally, minor mergers are unlikely to play a significant role, since the same conditions that lead to more major mergers should should also lead to more minor mergers, which is not observed in SG0023 & SC1604.
We have observed the dust continuum of ten z=3.1 Lyman Break Galaxies with the Atacama Large Millimeter/Submillimeter Array at ~450 mas resolution in Band 7. We detect and resolve the 870um emission in one of the targets with an integrated flux density of S(870)=(192+/-57) uJy, and measure a stacked 3-sigma signal of S(870)=(67+/-23) uJy for the remaining nine. The total infrared luminosities estimated from full spectral energy distribution fits are L(8-1000um)=(8.4+/-2.3)x10^10 Lsun for the detection and L(8-1000um)=(2.9+/-0.9)x10^10 Lsun for the stack. With HST ACS I-band imaging we map the rest-frame UV emission on the same scale as the dust, effectively resolving the 'infrared excess' (IRX=L_FIR/L_UV) in a normal galaxy at z=3. Integrated over the galaxy we measure IRX=0.56+/-0.15, and the galaxy-averaged UV slope is beta=-1.25+/-0.03. This puts the galaxy a factor of ~10 below the IRX-beta relation for local starburst nuclei of Meurer et al. (1999). However, IRX varies by more than a factor of 3 across the galaxy, and we conclude that the complex relative morphology of the dust relative to UV emission is largely responsible for the scatter in the IRX-beta relation at high-z. A naive application of a Meurer-like dust correction based on the UV slope would dramatically over-estimate the total star formation rate, and our results support growing evidence that when integrated over the galaxy, the typical conditions in high-z star-forming galaxies are not analogous to those in the local starburst nuclei used to establish the Meurer relation.
We report the discovery of a very diverse set of five low-surface brightness (LSB) dwarf galaxy candidates in Hickson Compact Group 90 (HCG 90) detected in deep U- and I-band images obtained with VLT/VIMOS. These are the first LSB dwarf galaxy candidates found in a compact group of galaxies. We measure spheroid half-light radii in the range $0.7\!\lesssim\! r_{\rm eff}/{\rm kpc}\! \lesssim\! 1.5$ with luminosities of $-11.65\!\lesssim\! M_U\! \lesssim\! -9.42$ and $-12.79\!\lesssim\! M_I\! \lesssim\! -10.58$ mag, corresponding to a color range of $(U\!-\!I)_0\!\simeq\!1.1\!-\!2.2$ mag and surface brightness levels of $\mu_U\!\simeq\!28.1\,{\rm mag/arcsec^2}$ and $\mu_I\!\simeq\!27.4\,{\rm mag/arcsec^2}$. Their colours and luminosities are consistent with a diverse set of stellar population properties. Assuming solar and 0.02 Z$_\odot$ metallicities we obtain stellar masses in the range $M_*|_{Z_\odot} \simeq 10^{5.7-6.3} M_{\odot}$ and $M_*|_{0.02\,Z_\odot}\!\simeq\!10^{6.3-8}\,M_{\odot}$. Three dwarfs are older than 1 Gyr, while the other two significantly bluer dwarfs are younger than $\sim 2$ Gyr at any mass/metallicity combination. Altogether, the new LSB dwarf galaxy candidates share properties with dwarf galaxies found throughout the Local Volume and in nearby galaxy clusters such as Fornax. We find a pair of candidates with $\sim\!2$ kpc projected separation, which may represent one of the closest dwarf galaxy pairs found. We also find a nucleated dwarf candidate, with a nucleus size of $r_{\rm eff}\!\simeq\!46\!-\!63$ pc and magnitude M$_{U,0}=-7.42$ mag and $(U\!-\!I)_0\!=\!1.51$ mag, which is consistent with a nuclear stellar disc with a stellar mass in the range $10^{4.9-6.5}\,M_\odot$.
We test for galactic conformity at $0.2<z<1.0$ to a projected distance of 5 Mpc using spectroscopic redshifts from the PRism MUlti-object Survey (PRIMUS). Our sample consists of $\sim60,000$ galaxies in five separate fields covering a total of $\sim5.5$ square degrees, which allows us to account for cosmic variance. We identify star-forming and quiescent "isolated primary" (i.e., central) galaxies using isolation criteria and cuts in specific star formation rate. We match the redshift and stellar mass distributions of these samples, to control for correlations between quiescent fraction and redshift and stellar mass. We detect a significant $(>3\sigma)$ one-halo conformity signal, or an excess of star-forming neighbors around star-forming central galaxies, of $\sim5$% on scales of 0-1 Mpc and a $2.5\sigma$ two-halo signal of $\sim1$% on scales of 1-3 Mpc. These signals are weaker than those detected in SDSS and are consistent with galactic conformity being the result of large-scale tidal fields and reflecting assembly bias. We also measure the star-forming fraction of central galaxies at fixed stellar mass as a function of large-scale environment and find that central galaxies are more likely to be quenched in overdense environments, independent of stellar mass. However, we find that environment does not affect the star formation efficiency of central galaxies, as long as they are forming stars. We test for redshift and stellar mass dependence of the conformity signal within our sample and show that large volumes and multiple fields are required at intermediate redshift to adequately account for cosmic variance.
The ultra-faint satellite galaxy Hercules has a strongly elongated and irregular morphology with detections of tidal features up to 1.3 deg (3 kpc) from its center. This suggests that Hercules may be dissolving under the Milky Way's gravitational influence, and hence could be a tidal stream in formation rather than a bound, dark-matter dominated satellite. Using Bayesian inference in combination with N-body simulations, we show that Hercules has to be on a very eccentric orbit (epsilon~0.95) within the Milky Way in this scenario. On such an orbit, Hercules "explodes" as a consequence of the last tidal shock at pericenter 0.5 Gyr ago. It is currently decelerating towards apocenter of its orbit with a velocity of V=157 km/s -- of which 99% is directed radially outwards. Due to differential orbital precession caused by the non-spherical nature of the Galactic potential, its debris fans out nearly perpendicular to its orbit. This explains why Hercules has an elongated shape without showing a distance gradient along its main body: it is in fact a stream that is significantly broader than it is long. In other words, it is moving perpendicular to its apparent major axis. In this scenario, there is a spike in the radial velocity profile created by the dominant debris component that formed through the last pericenter passage. This is similar to kinematic substructure that is observed in the real Hercules. Modeling a satellite on such a highly eccentric orbit is strongly dependent on the form of the Galactic potential. We therefore propose that detailed kinematic investigation of Hercules and other exploding satellite candidates can yield strong constraints on the potential of the Milky Way.
We compare predictions of a number of empirical models and numerical simulations of galaxy formation to the conditional stellar mass functions (CSMF) of galaxies in groups of different masses obtained recently by Lan et al. to test how well different models accommodate the data. Among all the models considered, only the model of Lu et al. can match the observational data; all other models fail to reproduce the faint-end upturn seen in the observation. The CSMFs are used to update the halo-based empirical model of Lu et al., and the model parameters obtained are very similar to those inferred by Lu et al. from a completely different set of observational constraints. The observational data clearly prefer a model in which star formation in low-mass halos changes behavior at a characteristic redshift $z_c \sim 2$. There is also tentative evidence that this characteristic redshift depends on environments, becoming $z_c \sim 4$ in regions that eventually evolve into rich clusters of galaxies. The constrained model is used to understand how galaxies form and evolve in dark matter halos, and to make predictions for other statistical properties of the galaxy population, such as the stellar mass functions of galaxies at high z, the star formation and stellar mass assembly histories in dark matter halos. A comparison of our model predictions with those of other empirical models shows that different models can make vastly different predictions, even though all of them are tuned to match the observed stellar mass functions of galaxies.
We present the catalogue of the Mg II absorption systems detected at a high significance level using an automated search algorithm in the spectra of quasars from the twelfth data release of the Sloan Digital Sky Survey. A total of 266,433 background quasars were searched for the presence of absorption systems in their spectra. The continuum modelling for the quasar spectra was performed using a mean filter. A pseudo-continuum derived using a median filter was used to trace the emission lines. The absorption system catalogue contains 39,694 Mg II systems detected at a 6.0, 3.0$\sigma$ level respectively for the two lines of the doublet. The catalogue was constrained to an absorption line redshift of 0.35 $\le$ z$_{2796}$ $\le$ 2.3. The rest-frame equivalent width of the $\lambda$2796 line ranges between 0.2 $\le$ W$_r$ $\le$ 6.2 \AA. Using Gaussian-noise only simulations we estimate a false positive rate of 7.7 per cent in the catalogue. We measured the number density $\partial N^{2796}/\partial z$ of Mg II absorbers and find evidence for steeper evolution of the systems with W$_r \ge$ 1.2 \AA\ at low redshifts (z$_{2796}$ $\le$ 1.0), consistent with other earlier studies. A suite of null tests over the redshift range 0.5 $\le$ z$_{2796}$ $\le$ 1.5 was used to study the presence of systematics and selection effects like the dependence of the number density evolution of the absorption systems on the properties of the background quasar spectra. The null tests do not indicate the presence of any selection effects in the absorption catalogue if the quasars with spectral signal-to-noise level less than 5.0 are removed. The resultant catalogue contains 36,981 absorption systems. The Mg II absorption catalogue is publicly available.
We present the X-ray spectral analysis of the 1855 extragalactic sources in the Chandra COSMOS-Legacy survey catalog having more than 30 net counts in the 0.5-7 keV band. 38% of the sources are optically classified Type 1 active galactic nuclei (AGN), 60% are Type 2 AGN and 2% are passive, low-redshift galaxies. We study the distribution of AGN photon index and of the intrinsic absorption N(H,z) based on the sources optical classification: Type 1 have a slightly steeper mean photon index than Type 2 AGN, which on the other hand have average intrinsic absorption ~3 times higher than Type 1 AGN. We find that ~15% of Type 1 AGN have N(H,z)>1E22 cm^(-2), i.e., are obscured according to the X-ray spectral fitting; the vast majority of these sources have L(2-10keV)>$1E44 erg/s. The existence of these objects suggests that optical and X-ray obscuration can be caused by different phenomena, the X-ray obscuration being for example caused by dust-free material surrounding the inner part of the nuclei. ~18% of Type 2 AGN have N(H,z)<1E22 cm^(-2), and most of these sources have low X-ray luminosities (L(2-10keV)<$1E43 erg/s). We expect a part of these sources to be low-accretion, unobscured AGN lacking of broad emission lines. Finally, we also find a direct proportional trend between N(H,z) and host galaxy mass and star formation rate, although part of this trend is due to a redshift selection effect.
Gas-grain and gas-phase reactions dominate the formation of molecules in the interstellar medium (ISM). Gas-grain reactions require a substrate (e.g. a dust or ice grain) on which the reaction is able to occur. The formation of molecular hydrogen (H$_2$) in the ISM is the prototypical example of a gas-grain reaction. In these reactions, an atom of hydrogen will strike a surface, stick to it, and diffuse across it. When it encounters another adsorbed hydrogen atom, the two can react to form molecular hydrogen and then be ejected from the surface by the energy released in the reaction. We perform in-depth classical molecular dynamics (MD) simulations of hydrogen atoms interacting with an amorphous water-ice surface. This study focuses on the first step in the formation process; the sticking of the hydrogen atom to the substrate. We find that careful attention must be paid in dealing with the ambiguities in defining a sticking event. The technical definition of a sticking event will affect the computed sticking probabilities and coefficients. Here, using our new definition of a sticking event, we report sticking probabilities and sticking coefficients for nine different incident kinetic energies of hydrogen atoms [5 K - 400 K] across seven different temperatures of dust grains [10 K - 70 K]. We find that probabilities and coefficients vary both as a function of grain temperature and incident kinetic energy over the range of 0.99 - 0.22.
By following the Kazantsev theory and taking into account both microscopic and turbulent diffusion of magnetic fields, we develop a unified treatment of the kinematic and nonlinear stages of turbulent dynamo, and study the dynamo process for a full range of magnetic Prandtl number Pm and ionization fractions. We find a striking similarity between the dependence of dynamo behavior on Pm in a conducting fluid and R (a function of ionization fraction) in partially ionized gas. In a weakly ionized medium, the kinematic stage is largely extended, including not only exponential growth but a new regime of dynamo characterized by linear-in-time growth of magnetic field strength, and the resulting magnetic energy is much higher than the kinetic energy carried by viscous-scale eddies. Unlike the kinematic stage, the subsequent nonlinear stage is unaffected by microscopic diffusion processes and has a universal linear-in-time growth of magnetic energy with the growth rate as a constant fraction $3/38$ of the turbulent energy transfer rate, showing good agreement with earlier numerical results. Applying the analysis to the first stars and galaxies, we find that the kinematic stage is able to generate a field strength only an order of magnitude smaller than the final saturation value. But the generation of large-scale magnetic fields can only be accounted for by the relatively inefficient nonlinear stage and requires longer time than the free-fall time. It suggests that magnetic fields may not have played a dynamically important role during the formation of the first stars.
Recent studies based on the integrated light of distant galaxies suggest that the initial mass function (IMF) might not be universal. Variations of the IMF with galaxy type and/or formation time may have important consequences for our understanding of galaxy evolution. We have developed a new stellar population synthesis (SPS) code specifically designed to reconstruct the IMF. We implement a novel approach combining regularization with hierarchical Bayesian inference. Within this approach we use a parametrized IMF prior to regulate a direct inference of the IMF. This direct inference gives more freedom to the IMF and allows the model to deviate from parametrized models when demanded by the data. We use Markov Chain Monte Carlo sampling techniques to reconstruct the best parameters for the IMF prior, the age, and the metallicity of a single stellar population. We present our code and apply our model to a number of mock single stellar populations with different ages, metallicities, and IMFs. When systematic uncertainties are not significant, we are able to reconstruct the input parameters that were used to create the mock populations. Our results show that if systematic uncertainties do play a role, this may introduce a bias on the results. Therefore, it is important to objectively compare different ingredients of SPS models. Through its Bayesian framework, our model is well-suited for this.
The majority of recent hydrodynamical simulations indicate the creation of central "cores" in the mass profiles of low-mass halos, a process that is attributed to star formation-related baryonic feedback. Core creation is regarded as one of the most promising solutions to potential issues faced by LambdaCDM cosmology on small scales. For example, the reduced dynamical mass enclosed by cores can explain the low rotational velocities measured for nearby dwarf galaxies, thus possibly lifting the seeming contradiction with the LambdaCDM expectations (the so-called "too big to fail" problem). Here we test core creation as a solution of cosmological issues by using a sample of dwarfs with measurements of their atomic hydrogen (HI) kinematics extending to large radii. Using the NIHAO hydrodynamical simulation as an example, we show that core creation can successfully reproduce the kinematics of dwarfs with small kinematic radii, R <~ 1.5 kpc. However, the agreement with observations becomes poor once galaxies with kinematic measurements extending beyond the core region, R ~ 1.5 - 4 kpc, are considered. This result illustrates the importance of testing the predictions of hydrodynamical simulations that are relevant for cosmology against a broad range of observational samples. We would like to stress that our result is valid only under the following set of assumptions: i) that our sample of dwarfs with HI kinematics is representative of the overall population of field dwarfs, ii) that there are no severe measurement biases in the observational parameters of our HI dwarfs (e.g., related to inclination estimates), and iii) that the HI velocity fields of dwarfs are regular enough to allow the recovery of the true enclosed dynamical mass.
Stars form in dense, dusty structures, which are embedded in larger clumps of molecular clouds often showing a clear filamentary structure on large scales (> 1pc). One of the best-studied regions in the Hi-GAL survey can be observed toward the l=224deg field. Here, a filamentary region has been studied and it has been found that protostellar clumps are mostly located along the main filament, whereas starless clumps are detected off this filament and are instead found on secondary, less prominent filaments. We want to investigate this segregation effect and how it may affect the clumps properties. We mapped the 12CO(1-0) line and its main three isotopologues toward the two most prominent filaments observed toward the l=224deg field using the Mopra radio telescope, in order to set observational constraints on the dynamics of these structures and the associated starless and protostellar clumps. Compared to the starless clumps, the protostellar clumps are more luminous, more turbulent and lie in regions where the filamentary ambient gas shows larger linewidths. We see evidence of gas flowing along the main filament, but we do not find any signs of accretion flow from the filament onto the Hi-GAL clumps. We analyze the radial column density profile of the filaments and their gravitational stability. The more massive and highly fragmented main filament appears to be thermally supercritical and gravitationally bound, assuming that all of the non-thermal motion is contributing thermal-like support, suggesting a later stage of evolution compared to the secondary filament. The status and evolutionary phase of the Hi-GAL clumps would then appear to correlate with that of the host filament.
The Breakthrough Starshot initiative aims to launch a gram-scale spacecraft to a speed of $v\sim 0.2$c, capable of reaching the nearest star system, $\alpha$ Centauri, in about 20 years. However, a critical challenge for the initiative is the damage to the spacecraft by interstellar gas and dust during the journey. In this paper, we quantify the interaction of a relativistic spacecraft with gas and dust in the interstellar medium. For gas bombardment, we find that damage by track formation due to heavy elements is an important effect. We find that gas bombardment can potentially damage the surface of the spacecraft to a depth of $\sim 0.1$ mm for quartz material after traversing a gas column of $N_{\rm H}\sim 2\times 10^{18}\rm cm^{-2}$ along the path to $\alpha$ Centauri, whereas the effect is much weaker for graphite material. The effect of dust bombardment erodes the spacecraft surface and produces numerous craters due to explosive evaporation of surface atoms. For a spacecraft speed $v=0.2c$, we find that dust bombardment can erode a surface layer of $\sim 0.5$ mm thickness after the spacecraft has swept a column density of $N_{\rm H}\sim 3\times 10^{17}\rm cm^{-2}$, assuming the standard gas-to-dust ratio of the interstellar medium. Dust bombardment also damages the spacecraft surface by modifying the material structure through melting. We calculate the equilibrium surface temperature due to collisional heating by gas atoms as well as the temperature profile as a function of depth into the spacecraft. Our quantitative results suggest methods for damage control, and we highlight possibilities for shielding strategies and protection of the spacecraft.
We present large-field (3x2 deg^2) and high-sensitivity CO(1-0) molecular line observations toward the Tycho's supernova remnant, using the 13.7-meter radio telescope of the Purple Mountain Observatory. Based on the CO observations, we discover a large cavity around the remnant, with radii of about 0.3x0.6 deg (or ~13x27 pc at a distance of 2.5 kpc), which is further supported by the complementary infrared images from the space telescopes. The observed CO line broadenings and asymmetries in the surrounding clouds, the infrared pillar-like structures found around the remnant, in concert with enhanced 12CO(2-1)/(1-0) intensity ratio detected in previous studies, indicate strong interaction of the large cavity with a wind in the region. After excluding the scenario of a large bubble produced by bright massive stars, we consider that the large cavity could be most likely explained by the accretion wind from the progenitor system of the Tycho's supernova. The CO gas kinematics indicates that the large cavity is expanding at a velocity of about 4 km/s. The estimated velocity (~1000 km/s, with a mass-loss rate of ~10^(-6)*M_sun*yr^(-1)) and timescale (~4x10^6 yr) of the wind needed for creating such a cavity are consistent with the predictions from the wind-regulated accretion model. We conclude that Tycho's supernova, the prototypical Type-Ia supernova in the Milky Way, arose from accretion onto a white dwarf.
Using synthetic absorption lines generated from 3D hydro-dynamical simulations we explore how the velocity of a starburst-driven galactic wind correlates with the star formation rate (SFR) and SFR density. We find strong correlations until the scaling relations flatten abruptly at a point set by the mass loading of the starburst. Below this point the scaling relation depends on the temperature regime being probed by the absorption line, not on the mass loading. The exact scaling relation depends on whether the maximum or mean velocity of the absorption line is used. We find that the outflow velocity of neutral gas is four to five times lower than the average velocity of the hottest gas, with the difference in velocity between the neutral and ionized gas increasing with gas ionization. Thus, absorption lines of neutral or low ionized gas will underestimate the outflow velocity of hot gas, severely underestimating outflow energetics.
The local Galactic HI gas was found to contain cold neutral medium (CNM) filaments that are aligned with polarized dust emission. These filaments appear to be dominated by the magnetic field and in this case turbulence is expected to show distinct anisotropies. We use the Galactic Effelsberg--Bonn HI Survey (EBHIS) to derive 2D turbulence spectra for the HI distribution in direction to 3C196 and two more comparison fields. Prior to Fourier transform we apply a rotational symmetric 50% Tukey window to apodize the data. We derive average as well as position angle dependent power spectra. Anisotropies in the power distribution are defined as the ratio of the spectral power in orthogonal directions. We find strong anisotropies. For a narrow range in position angle, in direction perpendicular to the filaments and the magnetic field, the spectral power is on average more than an order of magnitude larger than parallel. In the most extreme case the anisotropy reaches locally a factor of 130. Anisotropies increase on average with spatial frequency as predicted by Goldreich and Sridhar, at the same time the Kolmogorov spectral index remains almost unchanged. The strongest anisotropies are observable for a narrow range in velocity and decay with a power law index close to --8/3, almost identical to the average isotropic spectral index of $-2.9 < \gamma < -2.6$. HI filaments, associated with linear polarization structures in LOFAR observations in direction to 3C196, show turbulence spectra with marked anisotropies. Decaying anisotropies appear to indicate that we witness an ongoing shock passing the HI and affecting the observed Faraday depth.
Links to: arXiv, form interface, find, astro-ph, recent, 1608, contact, help (Access key information)