We examine the relation between stellar mass and central stellar velocity dispersion-the M-sigma relation-for massive quiescent galaxies at z<0.7. We measure the local relation from the Sloan Digital Sky Survey and the intermediate redshift relation from the Smithsonian Hectospec Lensing Survey. Both samples are highly complete (>85%) and we consistently measure the stellar mass and velocity dispersion for the two samples. The M-sigma relation and its scatter are independent of redshift with sigma ~ M^0.3 for M > 10^10.3 solar masses. The measured slope of the M-sigma relation is the same as the scaling between the total halo mass and the dark matter halo velocity dispersion obtained by N-body simulations. This consistency suggests that massive quiescent galaxies are virialized systems where the central dark matter concentration is either a constant or negligible fraction of the stellar mass. The relation between the total galaxy mass (stellar + dark matter) and the central stellar velocity dispersion is consistent with the observed relation between the total mass of a galaxy cluster and the velocity dispersion of the cluster members. This result suggests that the central stellar velocity dispersion is directly proportional to the velocity dispersion of the dark matter halo. Thus the central stellar velocity dispersion is a fundamental, directly observable property of galaxies that may robustly connect galaxies to dark matter halos in N-body simulations. To interpret the results further in the context of Lambda-CDM, it would be useful to analyze the relationship between the velocity dispersion of stellar particles and the velocity dispersion characterizing their dark matter halos in high-resolution cosmological hydrodynamic simulations.
We study how outflows of gas launched from a central galaxy undergoing repeated starbursts propagate through the circumgalactic medium (CGM), using the simulation code RAMSES. We assume that the outflow from the disk can be modelled as a rapidly moving bubble of hot gas at $\mathrm{\sim1\;kpc}$ above disk, then ask what happens as it moves out further into the halo around the galaxy on $\mathrm{\sim 100\;kpc}$ scales. To do this we run 60 two-dimensional simulations scanning over parameters of the outflow. Each of these is repeated with and without radiative cooling, assuming a primordial gas composition to give a lower bound on the importance of cooling. In a large fraction of radiative-cooling cases we are able to form rapidly outflowing cool gas from in situ cooling of the flow. We show that the amount of cool gas formed depends strongly on the 'burstiness' of energy injection; sharper, stronger bursts typically lead to a larger fraction of cool gas forming in the outflow. The abundance ratio of ions in the CGM may therefore change in response to the detailed historical pattern of star formation. For instance, outflows generated by star formation with short, intense bursts contain up to 60 per cent of their gas mass at temperatures $<5 \times 10^4\,\mathrm{K}$; for near-continuous star formation the figure is $\lesssim$ 5 per cent. Further study of cosmological simulations, and of idealised simulations with e.g., metal-cooling, magnetic fields and/or thermal conduction, will help to understand the precise signature of bursty outflows on observed ion abundances.
We present a new exploration of the cosmic star-formation history and dust obscuration in massive galaxies at redshifts $0.5< z<6$. We utilize the deepest 450 and 850$\mu$m imaging from SCUBA-2 CLS, covering 230arcmin$^2$ in the AEGIS, COSMOS and UDS fields, together with 100-250$\mu$m imaging from Herschel. We demonstrate the capability of the T-PHOT deconfusion code to reach below the confusion limit, using multi-wavelength prior catalogues from CANDELS/3D-HST. By combining IR and UV data, we measure the relationship between total star-formation rate (SFR) and stellar mass up to $z\sim5$, indicating that UV-derived dust corrections underestimate the SFR in massive galaxies. We investigate the relationship between obscuration and the UV slope (the IRX-$\beta$ relation) in our sample, which is similar to that of low-redshift starburst galaxies, although it deviates at high stellar masses. Our data provide new measurements of the total SFR density (SFRD) in $M_\ast>10^{10}M_\odot$ galaxies at $0.5<z<6$. This is dominated by obscured star formation by a factor of $>10$. One third of this is accounted for by 450$\mu$m-detected sources, while one fifth is attributed to UV-luminous sources (brighter than $L^\ast_{UV}$), although even these are largely obscured. By extrapolating our results to include all stellar masses, we estimate a total SFRD that is in good agreement with previous results from IR and UV data at $z\lesssim3$, and from UV-only data at $z\sim5$. The cosmic star-formation history undergoes a transition at $z\sim3-4$, as predominantly unobscured growth in the early Universe is overtaken by obscured star formation, driven by the build-up of the most massive galaxies during the peak of cosmic assembly.
We present a dynamical friction model based on Chandrasekhar's formula that reproduces the fast inspiral and stalling experienced by satellites orbiting galaxies with a large constant density core. We show that the fast inspiral phase does not owe to resonance. Rather, it owes to the background velocity distribution function for the constant density cores being dissimilar from the usually-assumed Maxwellian distribution. Using the correct background velocity distribution function and the semi-analytic model from Petts et al. (2015), we are able to correctly reproduce the infall rate in both cored and cusped potentials. However, in the case of large cores, our model is no longer able to correctly capture core-stalling. We show that this stalling owes to the tidal radius of the satellite approaching the size of the core. By switching off dynamical friction when rt(r) = r (where rt is the tidal radius at the satellite's position) we arrive at a model which reproduces the N-body results remarkably well. Since the tidal radius can be very large for constant density background distributions, our model recovers the result that stalling can occur for Ms/Menc << 1, where Ms and Menc are the mass of the satellite and the enclosed galaxy mass, respectively. Finally, we include the contribution to dynamical friction that comes from stars moving faster than the satellite. This next-to-leading order effect becomes the dominant driver of inspiral near the core region, prior to stalling.
We attempt to interpret existing data on the evolution of the UV luminosity function and UV colours, $\beta$, of galaxies at $5 \leq z \leq 8$, to improve our understanding of their dust content and ISM properties. To this aim, we post-process the results of a cosmological hydrodynamical simulation with a chemical evolution model, which includes dust formation by supernovae and intermediate mass stars, dust destruction in supernova shocks, and grain growth by accretion of gas-phase elements in dense gas. We find that observations require a steep, Small Magellanic Cloud-like extinction curve and a clumpy dust distribution, where stellar populations younger than 15 Myr are still embedded in their dusty natal clouds. Investigating the scatter in the colour distribution and stellar mass, we find that the observed trends can be explained by the presence of two populations: younger, less massive galaxies where dust enrichment is mainly due to stellar sources, and massive, more chemically evolved ones, where efficient grain growth provides the dominant contribution to the total dust mass. Computing the IR-excess - UV color relation we find that all but the dustiest model galaxies follow a relation shallower than the Meurer et al. (1999) one, usually adopted to correct the observed UV luminosities of high-$z$ galaxies for the effects of dust extinction. As a result, their total star formation rates might have been over-estimated. Our study illustrates the importance to incorporate a proper treatment of dust in simulations of high-$z$ galaxies, and that massive, dusty, UV-faint galaxies might have already appeared at $z \lesssim 7$.
Pop III stars are typically massive stars of primordial composition forming at the centers of the first collapsed dark matter structures. Here we estimate the optimal X-ray emission in the early universe for promoting the formation of Pop III stars. This is important in determining the number of dwarf galaxies formed before reionization and their fossils in the local universe, as well as the number of intermediate-mass seed black holes. A mean X-ray emission per source above the optimal level reduces the number of Pop III stars because of the increased Jeans mass of the intergalactic medium (IGM), while a lower emission suppresses the formation rate of H2 preventing or delaying star formation in dark matter minihalos above the Jeans mass. The build up of the H2 dissociating background is slower than the X-ray background due to the shielding effect of resonant hydrogen Lyman lines. Hence, the nearly unavoidable X-ray emission from supernova remnants of Pop III stars is sufficient to boost their number to few tens per comoving Mpc^3 by redshift z~15. We find that there is a critical X-ray to UV energy ratio emitted per source that produces a universe where the number of Pop III stars is largest: 400 per comoving- Mpc^3. This critical ratio is very close to the one provided by 20-40 M_sun Pop III stars exploding as hypernovae. High mass X-ray binaries in dwarf galaxies are far less effective at increasing the number of Pop III stars than normal supernova remnants, we thus conclude that supernovae drove the formation of Pop III stars.
In this paper, the first in a series on galaxy formation before reionization, we focus on understanding what determines the size and morphology of stellar objects in the first low mass galaxies, using parsec- scale cosmological simulations performed with an adaptive mesh hydrodynamics code. Although the dense gas in which stars are formed tends to have a disk structure, stars are found in spheroids with little rotation. Halos with masses between 10^6 M_sun and 5x10^8 M_sun form stars stochastically, with stellar masses in the range 10^4 M_sun to 2x10^6 M_sun. Nearly independent of stellar mass, we observe a large range of half-light radii for the stars, from a few parsecs to a few hundred parsecs and surface brightnesses and mass-to-light ratios ranging from those typical of globular clusters to ultra-faint dwarfs. In our simulations, stars form in dense stellar clusters with high gas-to-star conversion efficiencies and rather uniform metallicities. A fraction of these clusters remain bound after the gas is removed by feedback, but others are destroyed, and their stars, which typically have velocity dispersions of 20 to 40 km/s, expand until they become bound by the dark matter halo. We thus speculate that the stars in ultra-faint dwarf galaxies may show kinematic and chemical signatures consistent with their origin in a few distinct stellar clusters. On the other hand, some globular clusters may form at the center of primordial dwarf galaxies and may contain dark matter, perhaps detectable in the outer parts.
We use the Jeans equations for an ensemble of collisionless particles to describe the distribution of broad-line region (BLR) cloud in three classes: (A) non disc (B) disc-wind (c) pure disc structure. We propose that clumpy structures in the brightest quasars belong to class A, fainter quasars and brighter Seyferts belong to class B, and dimmer Seyfert galaxies and all low-luminosity AGNs (LLAGNs) belong to class C. We derive the virial factor, $f$, for disc-like structures and find a negative correlation between the inclination angle, $\theta_{0}$, and $f$. We find similar behaviour for $f$ as a function of the FWHM and $\sigma_{z}$, the $z$ component of velocity dispersion. For different values of $\theta_{0}$ we find that $ 1.0 \lesssim f \lesssim 9.0 $ in type1 AGNs and $ 0.5 \lesssim f \lesssim 1.0 $ in type2 AGNs. Moreover we have $ 0.5 \lesssim f \lesssim 6.5 $ for different values of FWHM and $ 1.4 \lesssim f \lesssim 1.8 $ for different values of $ \sigma_{z} $. We also find that $ f $ is relatively insensitive to the variations of bolometric luminosity and column density of each cloud and the range of variation of $ f $ is in order of 0.01. Considering wide range of $ f $ we see the use of average virial factor $ \langle f \rangle $ is not very safe. Therefore we propose AGN community to divide a sample into a few subsamples based on the value of $ \theta_{0}$ and FWHM of members and calculate $ \langle f \rangle $ for each group separately to reduce uncertainty in black hole mass estimation.
AGN exhibit rapid, high amplitude stochastic flux variability across the entire electromagnetic spectrum on timescales ranging from hours to years. The cause of this variability is poorly understood. We present a new method for using variability to (1) measure the time-scales on which flux perturbations evolve and (2) characterize the driving flux perturbations. We model the observed light curve of an AGN as a linear differential equation driven by stochastic impulses. Physically, the impulses could be local `hot-spots' in the accretion disk---the linear differential equation then governs how the hot spots evolve and dissipate. The impulse-response function of the accretion disk material is given by the Green's function of the linear differential equation. The timescales on which the hot-spots radiate energy is characterized by the powerspectrum of the driving stochastic impulses. We analyze the light curve of the \Kepler AGN Zw 229-15 and find that the observed variability behavior can be modeled as a damped harmonic oscillator perturbed by a colored noise process. The model powerspectrum turns over on time-scale $385$~d. On shorter time-scales, the log-powerspectrum slope varies between $2$ and $4$, explaining the behavior noted by previous studies. We recover and identify both the $5.6$~d and $67$~d timescales reported by previous work. These timescales represent the time-scale on which flux perturbations grow, and the time-scale on which flux perturbations decay back to the steady-state flux level respectively. We make the software package used to study light curves using our method, \textsc{k\={a}l\={i}}, available to the community.
This paper presents Gemini-$gri'$ high quality photometry for cluster candidates in the field of NGC 1316 (Fornax A) as part of a study that also includes GMOS spectroscopy. A preliminary discussion of the photometric data indicates the presence of four stellar cluster populations with distinctive features in terms of age, chemical abundance and spatial distribution. Two of them seem to be the usually old (metal poor and metal rich) populations typically found in elliptical galaxies. In turn, an intermediate-age (5 Gyr) globular cluster population is the dominant component of the sample (as reported by previous papers). We also find a younger cluster population with a tentative age of $\approx$ 1 Gyr.
We study possibility of obtaining velocity spectra by studying turbulence in an optically thick medium using velocity centroids (VCs).We find that the regime of universal, i.e. independent of underlying turbulence statistics, fluctuations discovered originally within the velocity channel analysis (VCA) carries over to the statistics of VCs. In other words, for large absorptions the VC lose their ability to reflect the spectra of turbulence. Combining our present study with the earlier studies of centroids in Esquivel & Lazarian, we conclude that centroids are applicable for studies subsonic/transsonic turbulence for the range of scales that is limited by the absorption effects. We also consider VCs based on absorption lines and define the range of their applicability. We address the problem of analytical description of spectra and anisotropies of fluctuations that are available through studies using VC. We obtain spectra and anisotropy of VC fluctuations arising from Alfv\'en, slow and fast modes that constitute the compressible MHD cascade to address the issue of anisotropy of VC statistics, and show how the VC anisotropy can be used to find the media magnetization as well as to identify and separate contributions from Alfv\'en, slow and fast modes. Our study demonstrates that VCs are complementary to the tools provided by the VCA. In order to study turbulent volume for which the resolution of single dish telescopes is insufficient, we demonstrate how the studies of anisotropy can be performed using interferometers.
Low resolution (4.5 to 5 Angstroem) spectra of 58 blue supergiant stars distributed over the disk of the Magellanic spiral galaxy NGC 55 in the Sculptor group are analyzed by means of non-LTE techniques to determine stellar temperatures, gravities and metallicities (from iron peak and alpha-elements). A metallicity gradient of -0.22 +/- 0.06$ dex/R_25 is detected. The central metallicity on a logarithmic scale relative to the Sun is [Z] = -0.37 +\- 0.03. A chemical evolution model using the observed distribution of stellar and interstellar medium gas mass column densities reproduces the observed metallicity distribution well and reveals a recent history of strong galactic mass accretion and wind outflows with accretion and mass-loss rates of the order of the star formation rate. There is an indication of spatial inhomogeneity in metallicity. In addition, the relatively high central metallicity of the disk confirms that two extra-planar metal poor HII regions detected in previous work 1.13 to 2.22 kpc above the galactic plane are ionized by massive stars formed in-situ outside the disk. For a sub-sample of supergiants, for which Hubble Space Telescope photometry is available, the flux-weighted gravity--luminosity relationship is used to determine a distance modulus of 26.85 +\- 0.10 mag.
Recent observations have detected molecular outflows in a few nearby starburst nuclei. We discuss the physical processes at work in such an environment in order to outline a scenario that can explain the observed parameters of the phenomenon, such as the molecular mass, speed and size of the outflows. We show that outflows triggered by OB associations, with $N_{OB}\ge 10^5$ (corresponding to a star formation rate (SFR)$\ge 1$ M$_{\odot}$ yr$^{-1}$ in the nuclear region), in a stratified disk with mid-plane density $n_0\sim 200\hbox{--}1000$ cm$^{-3}$ and scale height $z_0\ge 200 (n_0/10^2 \, {\rm cm}^{-3})^{-3/5}$ pc, can form molecules in a cool dense and expanding shell. The associated molecular mass is $\ge 10^7$ M$_\odot$ at a distance of a few hundred pc, with a speed of several tens of km s$^{-1}$. We show that a SFR surface density of $10 \le \Sigma_{SFR} \le 50$ M$_\odot$ yr$^{-1}$ kpc$^{-2}$ favours the production of molecular outflows, consistent with observed values.
The EAGLE cosmological simulations reproduce the observed galaxy stellar mass function and many galaxy properties. In this work, we study the dust-related properties of present-day EAGLE galaxies through mock observations in the far-infrared and submm wavelength ranges obtained with the 3D dust radiative transfer code SKIRT. To prepare an EAGLE galaxy for radiative transfer processing, we derive a diffuse dust distribution from the gas particles and we re-sample the star-forming gas particles and the youngest star particles into star-forming regions that are assigned dedicated emission templates. We select a set of redshift-zero EAGLE galaxies that matches the K-band luminosity distribution of the galaxies in the Herschel Reference Survey (HRS), a volume-limited sample of about 300 normal galaxies in the Local Universe. We find overall agreement of the EAGLE dust scaling relations with those observed in the HRS, such as the dust-to-stellar mass ratio versus stellar mass and versus NUV-r colour relations. A discrepancy in the f_250/f_350 versus f_350/f_500 submm colour-colour relation implies that part of the simulated dust is insufficiently heated, likely because of limitations in our sub-grid model for star-forming regions. We also investigate the effect of adjusting the metal-to-dust ratio and the covering factor of the photodissociation regions surrounding the star-forming cores. We are able to constrain the important dust-related parameters in our method, informing the calculation of dust attenuation for EAGLE galaxies in the UV and optical domain.
This is the second paper in a series on a new luminous z ~ 5 quasar survey using optical and near-infrared colors. Here we present a new determination of the bright end of the quasar luminosity function (QLF) at z ~ 5. Combined our 45 new quasars with previously known quasars that satisfy our selections, we construct the largest uniform luminous z ~ 5 quasar sample to date, with 99 quasars in the range 4.7 <= z < 5.4 and -29 < M1450 <= -26.8, within the Sloan Digital Sky Survey (SDSS) footprint. We use a modified 1/Va method including flux limit correction to derive a binned QLF, and we model the parametric QLF using maximum likelihood estimation. With the faint-end slope of the QLF fixed as alpha = -2.03 from previous deeper samples, the best fit of our QLF gives a flatter bright end slope beta = -3.58+/-0.24 and a fainter break magnitude M*1450 = -26.98+/-0.23 than previous studies at similar redshift. Combined with previous work at lower and higher redshifts, our result is consistent with a luminosity evolution and density evolution (LEDE) model. Using the best fit QLF, the contribution of quasars to the ionizing background at z ~ 5 is found to be 18% - 45% with a clumping factor C of 2 - 5. Our sample suggests an evolution of radio loud fraction with optical luminosity but no obvious evolution with redshift.
Clusters of galaxies are the most massive gravitationally-bound objects in the Universe and are still forming. They are thus important probes of cosmological parameters and a host of astrophysical processes. Knowledge of the dynamics of the pervasive hot gas, which dominates in mass over stars in a cluster, is a crucial missing ingredient. It can enable new insights into mechanical energy injection by the central supermassive black hole and the use of hydrostatic equilibrium for the determination of cluster masses. X-rays from the core of the Perseus cluster are emitted by the 50 million K diffuse hot plasma filling its gravitational potential well. The Active Galactic Nucleus of the central galaxy NGC1275 is pumping jetted energy into the surrounding intracluster medium, creating buoyant bubbles filled with relativistic plasma. These likely induce motions in the intracluster medium and heat the inner gas preventing runaway radiative cooling; a process known as Active Galactic Nucleus Feedback. Here we report on Hitomi X-ray observations of the Perseus cluster core, which reveal a remarkably quiescent atmosphere where the gas has a line-of-sight velocity dispersion of 164+/-10 km/s in a region 30-60 kpc from the central nucleus. A gradient in the line-of-sight velocity of 150+/-70 km/s is found across the 60 kpc image of the cluster core. Turbulent pressure support in the gas is 4% or less of the thermodynamic pressure, with large scale shear at most doubling that estimate. We infer that total cluster masses determined from hydrostatic equilibrium in the central regions need little correction for turbulent pressure.
The "torus" is the central element of the most popular theory unifying various classes of AGNs, but it is usually described as "putative" because it has not been imaged yet. Since it is too small to be resolved with single-dish telescopes, one can only make indirect assumptions about its structure using models. Using infrared interferometry, however, we were able to resolve the circum-nuclear dust distributions for several nearby AGNs and achieved constraints on some further two dozen sources. We discovered circum-nuclear dust on parsec scales in all sources and, in two nearby sources, were able to dissect this dust into two distinct components. The compact component, a very thin disk, appears to be connected to the maser disk and the extended one, which is responsible for most of the mid-IR flux, is oriented perpendicularly to the circum-nuclear gas disks. What may come as a surprise when having in mind the standard unification cartoon actually connects well to observations on larger scales. Optically thin dust in the polar region, perhaps driven by a disk wind, could solve both the scale height problem of the torus and explain the missing anisotropy in the mid-IR - X-ray relation.
We investigate an infrared-excess source called G2 or Dusty S-cluster Object (DSO) moving on a highly eccentric orbit around the Galaxy's central black hole, Sgr A*. We use, for the first time, near-infrared polarimetric imaging data to determine the nature and the properties of the DSO, and obtain an improved K_s-band identification of this source in median polarimetry images of different observing years. The source starts to deviate from the stellar confusion in 2008, and it does not show any flux density variability over the years we analyzed it. We measure the polarization degree and angle of the DSO between 2008 and 2012 and conclude, based on the significance analysis on polarization parameters, that it is an intrinsically polarized source (>20%) with a varying polarization angle as it approaches Sgr A* position. DSO shows a near-infrared excess of K_s-L' > 3 that remains compact close to the pericenter of its orbit. Its observed parameters and the significant polarization obtained in this work show that the DSO might be a dust-enshrouded young star, forming a bow shock as it approaches the super massive black hole. The significantly high measured polarization degree indicates that it has a non-spherical geometry and it can be modelled as a combination of a bow shock with a bipolar wind of the star. We use a 3D radiative transfer model that can reproduce the observed properties of the source such as the total flux density and the polarization degree. We obtain that the change of the polarization angle can be due to an intrinsic change in the source structure. Accretion disc precession of the young star in the gravitational field of the black hole can lead to the change of the bipolar outflow and therefore the polarization angle variation. It might also be the result of the source interaction with the ambient medium.
We use the SDSS and WISE surveys to investigate the real nature of galaxies defined as LINERs in the BPT diagram. After establishing a mid-infrared colour W2-W3 = 2.5 as the optimal separator between galaxies with and without star formation, we investigate the loci of different galaxy classes in the W_{Ha} versus W2-W3 space. We find that: (1) A large fraction of LINER-like galaxies are emission-line retired galaxies, i.e galaxies which have stopped forming stars and are powered by hot low-mass evolved stars (HOLMES). Their W2-W3 colours show no sign of star formation and their Ha equivalent widths, W_{Ha}, are consistent with ionization by their old stellar populations. (2) Another important fraction have W2-W3 indicative of star formation. This includes objects located in the supposedly `pure AGN' zone of the BPT diagram. (3) A smaller fraction of LINER-like galaxies have no trace of star formation from W2-W3 and a high W_{Ha}, pointing to the presence of an AGN. (4) Finally, a few LINERs tagged as retired by their W_{Ha} but with W2-W3 values indicative of star formation are late-type galaxies whose SDSS spectra cover only the old `retired' bulge. This reinforces the view that LINER-like galaxies are a mixed bag of objects involving different physical phenomena and observational effects thrusted into the same locus of the BPT diagram.
Aims. In this work we rigorously show the shortcomings of various
simplifications that are used to calculate the total internal partition
function. These shortcomings can lead to errors of up to 40 percent or more in
the estimated partition function. These errors carry on to calculations of
thermodynamic quantities. Therefore a more complicated approach has to be
taken.
Methods. Seven possible simplifications of various complexity are described,
together with advantages and disadvantages of direct summation of experimental
values. These were compared to what we consider the most accurate and most
complete treatment (case 8). Dunham coefficients were determined from
experimental and theoretical energy levels of a number of electronically
excited states of H$_2$ . Both equilibrium and normal hydrogen was taken into
consideration.
Results. Various shortcomings in existing calculations are demonstrated, and
the reasons for them are explained. New partition functions for equilibrium,
normal, and ortho and para hydrogen are calculated and thermodynamic quantities
are reported for the temperature range 1 - 20000 K. Our results are compared to
previous estimates in the literature. The calculations are not limited to the
ground electronic state, but include all bound and quasi-bound levels of
excited electronic states. Dunham coefficients of these states of H$_2$ are
also reported.
Conclusions. For most of the relevant astrophysical cases it is strongly
advised to avoid using simplifications, such as a harmonic oscillator and rigid
rotor or ad hoc summation limits of the eigenstates to estimate accurate
partition functions and to be particularly careful when using polynomial fits
to the computed values. Reported internal partition functions and thermodynamic
quantities in the present work are shown to be more accurate than previously
available data.
We present the MAssive ClusterS and Intercluster Structures (MACSIS) project, a suite of 390 clusters simulated with baryonic physics that yields realistic massive galaxy clusters capable of matching a wide range of observed properties. MACSIS extends the recent BAHAMAS simulation to higher masses, enabling robust predictions for the redshift evolution of cluster properties and an assessment of the effect of selecting only the hottest systems. We study the observable-total mass and X-ray luminosity-temperature scaling relations across the complete observed cluster mass range, finding the slope of the scaling relations and the evolution of their normalization with redshift to depart significantly from self-similar predictions. This is driven by the impact of AGN feedback, the presence of non-thermal pressure support and biased X-ray temperatures. For a sample of hot clusters with core-excised temperatures $k_{\rm{B}}T\geq5\,\rm{keV}$ the normalization and slope of the observable-mass relations and their evolution are significantly closer to self-similar. The exception is the temperature-mass relation, for which the increased importance of non-thermal pressure support and biased X-ray temperatures leads to a greater departure from self-similarity in the hottest systems. We also demonstrate that this affects the slope and evolution of the normalization in the luminosity-temperature relation. The median hot gas profiles also show good agreement with observational data at $z=0$ and $z=1$, with their evolution again departing significantly from the self-similar prediction. However, selecting a hot sample of clusters yields profiles that evolve significantly closer to the self-similar prediction.
Links to: arXiv, form interface, find, astro-ph, recent, 1607, contact, help (Access key information)
Here is considered the full evolution of a spherical supernova remnant. We start by calculating the early time ejecta dominated stage and continue through the different phases of interaction with the circumstellar medium, and end with the dissipation and merger phase. The physical connection between the phases reveals new results. One is that the blast wave radius during the adiabatic phase is significantly smaller than it would be, if one does notaccount for the blast wave interaction with the ejecta.
NGC 4660, in the Virgo cluster, is a well-studied elliptical galaxy which has a strong disk component (D/T about 0.2-0.3). The central regions including the disk component have stellar populations with ages about 12-13 Gyr from SAURON studies. However we report the discovery of a long narrow tidal filament associated with the galaxy in deep co-added Schmidt plate images and deep CCD frames, implying that the galaxy has undergone a tidal interaction and merger within the last few Gyr. The relative narrowness of the filament implies a wet merger with at least one spiral galaxy involved, but the current state of the system has little evidence for this. However a 2-component photometric fit using GALFIT shows much bluer B-V colours for the disk component than for the elliptical component, which may represent a residual trace of enhanced star formation in the disk caused by the interaction 1-2 Gyr ago. There are brighter concentrations within the filament which resemble Tidal Dwarf Galaxies, although they are at least 40 times fainter. These may represent faint, evolved versions of these galaxies. A previously detected stripped satellite galaxy south of the nucleus is seen in our residual image and may imply that the filament is a tidal stream produced by perigalactic passages of this satellite.
We present measurements of the clustering of galaxies as a function of their stellar mass in the Baryon Oscillation Spectroscopic Survey. We compare the clustering of samples using 12 different methods for estimating stellar mass, isolating the method that has the smallest scatter at fixed halo mass. In this test, the stellar mass estimate with the smallest errors yields the highest amplitude of clustering at fixed number density. We find that the PCA stellar masses of Chen etal (2012) clearly have the tightest correlation with halo mass. The PCA masses use the full galaxy spectrum, differentiating them from other estimates that only use optical photometric information. Using the PCA masses, we measure the large-scale bias as a function of Mgal for galaxies with logMgal>=11.4, correcting for incompleteness at the low-mass end of our measurements. Using the abundance-matching ansatz to connect dark matter halo mass to stellar mass, we construct theoretical models of b(Mgal) that match the same stellar mass function but have different amounts of scatter in stellar mass at fixed halo mass, sigma_logM. Using this approach, we find sigma_logM=0.18^{+0.01}_{-0.02}. This value includes both intrinsic scatter as well as random errors in the stellar masses. To partially remove the latter, we use repeated spectra to estimate statistical errors on the stellar masses, yielding an upper limit to the intrinsic scatter of 0.16 dex.
We present a statistical method for measuring the average HI spin temperature in distant galaxies using the expected detection yields from future wide-field 21-cm absorption surveys. As a demonstrative case study we consider a simulated survey with the Australian Square Kilometre Array Pathfinder for intervening HI absorbers at intermediate cosmological redshifts between $z = 0.4$ and 1. If such a survey yielded $\sim 1000$ absorbers we would infer a harmonic-mean spin temperature of $\overline{T}_\mathrm{spin} \sim 100$K for the population of Damped Lyman-$\alpha$ absorbers at these redshifts, indicating that more than $50$ percent of the neutral gas in these systems is in a cold neutral medium (CNM). Conversely, a lower yield of only 100 detections would imply $\overline{T}_\mathrm{spin} \sim 1000$K and a correspondingly lower CNM fraction. We propose that this method can be used to provide independent verification of the spin temperature evolution reported in recent 21-cm surveys of known DLAs at high redshift and for measuring the spin temperature at intermediate redshifts below $z \approx 1.7$, where the Lyman-$\alpha$ line is inaccessible using ground-based observatories. Increasingly more sensitive and larger surveys with the Square Kilometre Array should provide stronger statistical constraints on the average spin temperature. However, these will ultimately be limited by the accuracy to which we can determine the HI column density frequency distribution, the covering factor and the redshift distribution of the background radio source population.
The Fermi bubbles are two lobes filled with non-thermal particles that emit gamma rays, extend $\approx$10 kpc vertically from the Galactic center, and formed from either nuclear star formation or accretion activity on Sgr A*. Simulations predict a range of shock strengths as the bubbles expand into the surrounding hot gas halo distribution ($T_{halo} \approx 2 \times 10^6$ K), but with significant uncertainties in the energetics, age, and thermal gas structure. The bubbles should contain thermal gas with temperatures between $10^6$ and $10^8$ K, with potential X-ray signatures. In this work, we constrain the bubbles' thermal gas structure by modeling the OVII and OVIII emission line strengths from archival XMM-Newton and Suzaku data. Our emission model includes a hot thermal volume-filled bubble component cospatial with the gamma-ray region, and a shell of compressed material. We find that a bubble/shell model with $n \approx 1 \times 10^{-3}$ cm$^{-3}$ and with log($T$) $\approx$ 6.60-6.70 is consistent with the observed line intensities. In the framework of a continuous Galactic outflow, we infer a bubble expansion rate, age, and energy injection rate of $490_{-77}^{+230}$ km s$^{-1}$, $4.3_{-1.4}^{+0.8}$ Myr, and $2.3_{-0.9}^{+5.1} \times 10^{42}$ erg s$^{-1}$. These estimates are consistent with the bubbles forming from a Sgr A* accretion event rather than from nuclear star formation.
We present a polarization catalog of 533 extragalactic radio sources with 2.3 GHz total intensity above 420 mJy from the S-band Polarization All Sky Survey, S-PASS, with corresponding 1.4 GHz polarization information from the NRAO VLA Sky Survey, NVSS. We studied selection effects and found that fractional polarization, $\pi$, of radio objects at both wavelengths depends on the spectral index, source magnetic field disorder, source size and depolarization. The relationship between depolarization, spectrum and size shows that depolarization occurs primarily in the source vicinity. The median $\pi_{2.3}$ of resolved objects in NVSS is approximately two times larger than that of unresolved sources. Sources with little depolarization are $\sim2$ times more polarized than both highly depolarized and re-polarized sources. This indicates that intrinsic magnetic field disorder is the dominant mechanism responsible for the observed low fractional polarization of radio sources at high frequencies. We predict that number counts from polarization surveys will be similar at 1.4 GHz and at 2.3 GHz, for fixed sensitivity, although $\sim$10% of all sources may be currently missing because of strong depolarization. Objects with $\pi_{1.4}\approx \pi_{2.3} \ge 4\%$ typically have simple Faraday structures, so are most useful for background samples. Almost half of flat spectrum ($\alpha \ge -0.5$) and $\sim$25% of steep spectrum objects are re-polarized. Steep spectrum, depolarized sources show a weak negative correlation of depolarization with redshift in the range 0 $<$ z $<$ 2.3. Previous non-detections of redshift evolution are likely due the inclusion of re-polarized sources as well.
It is believed that, in radio-loud active galactic nuclei (AGN), the core radio flux density can be normalized to the flux density of the extended lobe emission to infer the orientation of a radio source. However very little is known about the reliability and precision of this method, and we are unaware of any robust conversion recipe to infer the inclination from the core dominance. Investigating whether or not the radio core dominance parameter R separates the quasars from the radio-galaxies in the $z \ge$~1 3CRR catalog, we found excellent agreement of R with optical type, infrared flux ratios and optical polarization. This indicates that probably both R and optical classification are very good orientation indicators, and the unified model is strongly predictive for these objects. The relative number densities indicate half-opening angles close to 60$^\circ$, as expected from large surveys. The separations of optical types according to radio core dominance as well as NIR/MIR ratios, which are essentially perfect, means that there can be only a small dispersion of torus half-opening angles. Also, even though torus dust is thought to be clumpy, there is an almost zero probability to see a type-1 source at high inclination. Finally, using only the Copernican Principle, i,e, the assumption that solid angle is filled uniformly with source axis orientations, we estimated a semi-empirical relation between core dominance and AGN inclination. This makes it possible to use R to infer the inclination of a source to an accuracy of $\sim$ 10 degrees or less, at least for this type of object.
We present the analysis of 35.5 square degrees of images in the 1-0S(1) line
of H2 from the UK Widefield Infrared Survey for H2 (UWISH2) towards Cassiopeia
and Auriga. We have identified 98 Molecular Hydrogen emission-line Objects
(MHOs) driven by Young Stellar Objects, 60% of which are bipolar outflows and
all are new discoveries. We estimate that the UWISH2 extended emission object
catalogue contains fewer than % false positives and is complete at the 95%
level for jets and outflows brighter than the UWISH2 detection limit. We
identified reliable driving source candidates for three quarters of the
detected outflows, 40% of which are associated with groups and clusters of
stars. The driving source candidates are 20% protostars, the remainder are
CTTSs. We also identified 15 new star cluster candidates near MHOs in the
survey area.
We find that the typical outflow identified in the sample has the following
characteristics: the position angles are randomly orientated; bipolar outflows
are straight within a few degrees; the two lobes are slightly asymmetrical in
length and brightness; the length and brightness of the lobes are not
correlated; typical time gaps between major ejections of material are 1-3kyr,
hence FU-Ori or EX-Ori eruptions are most likely not the cause of these, but we
suggest MNors as a possible source. Furthermore, we find that outflow lobe
length distributions are statistically different from the widely used total
length distributions. There are a larger than expected number of bright
outflows indicating that the flux distribution does not follow a power law.
We report the results of our radio, optical and infrared studies of a peculiar radio source 4C 35.06, an extended radio-loud AGN found at the center of galaxy cluster Abell 407 (z=0.047). The central region also hosts a remarkably tight ensemble of nine galaxies within an ~1' (53 kpc) region, surrounded by a very faint, diffuse stellar halo. This remarkable system (named here the 'Zwicky's Nonet') provides a unique and compelling evidence for an ongoing formation of a giant cD galaxy at the cluster center. Multifrequency deep radio observations carried out with Giant Meterwave Radio Telescope at 610, 235 and 150 MHz reveal a system of 400 kpc scale helically twisted and kinked radio jets and outer relic lobes associated with 4C 35.06 in 'Zwicky's Nonet'. The outer extremities of jets show extremely steep spectrum (spectral index -1.7 to -2.5) relic or fossil radio plasma with a relatively short spectral age of few times 10^6 - 10^7 years. Such ultra steep spectrum relic radio plasma lobes without definitive hot-spots are very rare and they provide an interesting opportunity towards understanding the life-cycle of relativistic jets and physics of black hole mergers. We discuss our observations of this complex radio galaxy in the context of growth of its central black hole, triggering of its AGN activity and jet oscillations; presumably all caused by mergers in this dense galactic system. A slow conical precession of the jet axis due to gravitational perturbation effects between interacting black holes is invoked to understand the peculiar jet morphology. The observed close resemblance of the morphology of 4C 35.06 with the precessing relativistic jets of the galactic microquasar SS 433 is noted, and scale-invariance of the disk-jet coupling processes in two systems over a large black hole mass range is discussed.
The role of the Hall term on large scale dynamo action is investigated by means of the First Order Smoothing Approximation. It is shown that the standard $\alpha$ coefficient is altered, and is zero when a specific double Beltrami state is attained, in contrast to the Alfv\'enic state for MHD dynamos. The $\beta$ coefficient is no longer positive definite, and thereby enables dynamo action even if $\alpha$-quenching were to operate. The similarities and differences with the (magnetic) shear-current effect are pointed out, and a mechanism that may be potentially responsible for $\beta < 0$ is advanced. The results are compared against previous studies, and their astrophysical relevance is also highlighted.
We have carried out a high-precision astrometric analysis of two very-long-baseline-interferometry (VLBI) epochs of observation of the 13 extragalactic radio sources in the complete S5 polar cap sample. The VLBI epochs span a time baseline of 10 years and enable us to achieve precisions in the proper motions of the source cores up to a few micro-arcseconds per year. The observations were performed at 14.4 GHz and 43.1 GHz, and enable us to estimate the frequency core-shifts in a subset of sources, for which the spectral-index distributions can be computed. We study the source-position stability by analysing the changes in the relative positions of fiducial source points (the jet cores) over a decade. We find motions of 0.1-0.9 mas among close-by sources between the two epochs, which imply drifts in the jet cores of approximately a few tens of micro-as per year. These results have implications for the standard Active Galactic Nucleus (AGN) jet model (where the core locations are supposed to be stable in time). For one of our sources, 0615+820, the morphological and spectral properties in year 2010, as well as the relative astrometry between years 2000 and 2010, suggest the possibility of either a strong parsec-scale interaction of the AGN jet with the ISM, a gravitational lens with ~1 mas diameter, or a resolved massive binary black hole.
I derive a new MOND relation for pure-disc galaxies: The `dynamical' central surface density, $\Sigma^0_D$, deduced from the measured velocities, is a universal function of only the true, `baryonic' central surface density, $\Sigma^0_B$: $\Sigma^0_D=\Sigma_M \mathcal{S}(\Sigma^0_B/\Sigma_M)$, where $\Sigma_M\equiv a_0/2\pi G$ is the MOND surface density constant. This surprising result is shown to hold in both existing, nonrelativistic MOND theories (the nonlinear Poisson formulation, and QUMOND). $\mathcal{S}(y)$ is derived, giving in the two limits: $\Sigma^0_D=\Sigma^0_B$ for very high arguments, and $\Sigma^0_D=(4\Sigma_M\Sigma^0_B)^{1/2}$ for $\Sigma^0_B/\Sigma_M\ll 1$. This study was prompted by the recent finding of a correlation between related attributes in a large sample of disc galaxies by Lelli et al. (2016). The MOND relation is shown to agree very well with these results.
We present an analysis of the complex gas hydrodynamics in the X-ray luminous galaxy cluster RXJ1347.5-1145 caught in the act of merging with a subcluster to its southeast using a combined $186$ ks Chandra exposure, $2.5$ times greater than previous analyses. The primary cluster hosts a sloshing cold front spiral traced by four surface brightness edges $5.^{\prime \prime}85^{+0.04}_{-0.03}$ west, $7.^{\prime \prime}10^{+0.07}_{-0.03}$ southeast, $11.^{\prime \prime}5^{+1.3}_{-1.2}$ east, and $16.^{\prime \prime}7^{+0.3}_{-0.5}$ northeast from the primary central dominant galaxy, suggesting the merger is in the plane of the sky. We measure temperature and density ratios across these edges, confirming they are sloshing cold fronts. We observe the eastern edge of the subcluster infall shock, confirming the observed subcluster is traveling from the southwest to the northeast in a clockwise orbit. We measure a shock density contrast of $1.38^{+0.16}_{-0.15}$ and infer a Mach number $1.25\pm0.08$ and a shock velocity of $2810^{+210}_{-240}$ km s$^{-1}$. Temperature and entropy maps show cool, low entropy gas trailing the subcluster in a southwestern tail, consistent with core shredding. Simulations suggest a perturber in the plane of the sky on a clockwise orbit would produce a sloshing spiral winding counterclockwise, opposite to that observed. The most compelling solution to this discrepancy is that the observed southeastern subcluster is on its first passage, shock heating gas during its clockwise infall, while the main cluster's clockwise cold front spiral formed from earlier encounters with a second perturber orbiting counterclockwise.
The CONT14 campaign with state-of-the-art VLBI data has observed the source 0642+449 with about one thousand observables each day during a continuous observing period of fifteen days, providing tens of thousands of closure delays---the sum of the delays around a closed loop of baselines. The closure delay is independent of the instrumental and propagation delays and provides valuable additional information about the source structure. We demonstrate the use of this new "observable" for the determination of the structure in the radio source 0642+449. This source, as one of the defining sources in the second realization of the International Celestial Reference Frame (ICRF2), is found to have two point-like components with a relative position offset of -426 microarcseconds in right ascension and -66 microarcseconds in declination. The two components are almost equally bright with a flux-density ratio of 0.92. The standard deviation of closure delays for source 0642+449 was reduced from 139 ps to 90 ps by using this two-component model. Closure delays larger than one nanosecond are found to be related to the source structure, demonstrating that structure effects for a source with this simple structure could be up to tens of nanoseconds. The method described in this paper does not rely on a priori source structure information, such as knowledge of source structure determined from direct (Fourier) imaging of the same observations or observations at other epochs. We anticipate our study to be a starting point for more effective determination of the structure effect in VLBI observations.
Correlations between intrinsic galaxy shapes on large-scales arise due to the effect of the tidal field of the large-scale structure. Anisotropic primordial non-Gaussianity induces a distinct scale-dependent imprint in these tidal alignments on large scales. Motivated by the observational finding that the alignment strength of luminous red galaxies depends on how galaxy shapes are measured, we study the use of two different shape estimators as a multi-tracer probe of intrinsic alignments. We show, by means of a Fisher analysis, that this technique promises a significant improvement on anisotropic non-Gaussianity constraints over a single-tracer method. For future weak lensing surveys, the uncertainty in the anisotropic non-Gaussianity parameter, $A_2$, is forecast to be $\sigma(A_2)\approx 50$, $\sim 40\%$ smaller than currently available constraints from the bispectrum of the Cosmic Microwave Background. This corresponds to an improvement of a factor of $4-5$ over the uncertainty from a single-tracer analysis.
Links to: arXiv, form interface, find, astro-ph, recent, 1607, contact, help (Access key information)
Using high-resolution simulations with explicit treatment of stellar feedback physics based on the FIRE (Feedback in Realistic Environments) project, we study how galaxy formation and the interstellar medium (ISM) are affected by magnetic fields, anisotropic Spitzer-Braginskii conduction and viscosity, and sub-grid turbulent metal diffusion. We consider controlled simulations of isolated (non-cosmological) galaxies but also a limited set of cosmological "zoom-in" simulations. Although simulations have shown significant effects from these physics with weak or absent stellar feedback, the effects are much weaker than those of stellar feedback when the latter is modeled explicitly. The additional physics have no systematic effect on galactic star formation rates (SFRs). In contrast, removing stellar feedback leads to SFRs being over-predicted by factors of $\sim 10 -100$. Without feedback, neither galactic winds nor volume filling hot-phase gas exist, and discs tend to runaway collapse to ultra-thin scale-heights with unphysically dense clumps congregating at the galactic center. With stellar feedback, a multi-phase, turbulent medium with galactic fountains and winds is established. At currently achievable resolutions, the additional physics investigated here (MHD, conduction, viscosity, metal diffusion) have only weak ($\sim10\%$-level) effects on these properties and do not significantly alter the balance of phases, outflows, or the energy in ISM turbulence, consistent with simple equipartition arguments. We conclude that galactic star formation and the ISM are primarily governed by a combination of turbulence, gravitational instabilities, and feedback.
The exponential scale length ($L_d$) of the Milky Way's (MW's) disk is a critical parameter for describing the global physical size of our Galaxy, important both for interpreting other Galactic measurements and helping us to understand how our Galaxy fits into extragalactic contexts. Unfortunately, current estimates span a wide range of values and often are statistically incompatible with one another. Here, we aim to determine an improved, aggregate estimate for $L_d$ by utilizing a hierarchical Bayesian (HB) meta-analysis technique that accounts for the possibility that any one measurement has not properly accounted for all statistical or systematic errors. Within this machinery we explore a variety of ways of modeling the nature of problematic measurements, and then use a Bayesian model averaging technique to derive net posterior distributions that incorporate any model-selection uncertainty. Our meta-analysis combines 29 different (15 visible and 14 infrared) photometric measurements of $L_d$ available in the literature; these involve a broad assortment of observational datasets, MW models and assumptions, and methodologies, all tabulated herein. Analyzing the visible and infrared measurements separately yields estimates for $L_d$ of $2.71^{+0.22}_{-0.20}$ kpc and $2.51^{+0.15}_{-0.13}$ kpc, respectively, whereas considering them all combined yields $2.64\pm0.13$ kpc. The ratio between the visible and infrared scale lengths determined here is very similar to that measured in external spiral galaxies. We use these results to updated the model of the Galactic disk from our previous work, constraining its stellar mass to be $4.8^{+1.5}_{-1.1}\times10^{10}$ $\mathrm{M}_\odot$, corresponding to a total stellar mass of $5.7^{+1.5}_{-1.1}\times10^{10}$ $\mathrm{M}_\odot$.
We use spectra from the ALFALFA, GASS and COLD GASS surveys to quantify variations in the mean atomic and molecular gas mass fractions throughout the SFR-M* plane and along the main sequence (MS) of star-forming galaxies. Although galaxies well below the MS tend to be undetected in the Arecibo and IRAM observations, reliable mean atomic and molecular gas fractions can be obtained through a spectral stacking technique. We find that the position of galaxies in the SFR-M* plane can be explained mostly by their global cold gas reservoirs as observed in the HI line, with in addition systematic variations in the molecular-to-atomic ratio and star formation efficiency. When looking at galaxies within +/-0.4 dex of the MS, we find that as stellar mass increases, both atomic and molecular gas mass fractions decrease, stellar bulges become more prominent, and the mean stellar ages increase. Both star formation efficiency and molecular-to-atomic ratios vary little for massive main sequence galaxies, indicating that the flattening of the MS is due to the global decrease of the cold gas reservoirs of galaxies rather than to bottlenecks in the process of converting cold atomic gas to stars.
The average bulge [Fe/H] and [Mg/H] are +0.06 and +0.17 dex, respectively, in Baade's Window, roughly 0.2 dex higher than the thin disk and ~0.7 dex higher than the local thick disk metallicity. This suggests a higher effective yield in the bulge, perhaps due to more efficient retention of supernova ejecta. The bulge vertical [Fe/H] gradient, at ~0.5 dex/kpc, appears to be due to a changing mixture of sub-populations (near +0.3 dex and -0.3 dex and one possibly near -0.7 dex) with latitude. The bulge is enhanced in O, Mg, Si, Ca, Ti, and Al relative to the sun, with [alpha/Fe]=+0.15 dex at [Fe/H]=0.0 dex. Below [Fe/H]~-0.5 dex, the bulge and local thick disk compositions are very similar, but small [Mg/Fe] and possibly [<SiCaTi>/Fe] enhancements, low [La/Eu] ratios and large [Cu/Fe], relative to the thick disk suggest slightly higher SFR in the bulge. However, these composition differences could simply be due to measurement errors and non-LTE effects. Unfortunately, comparison with the thick disk near solar [Fe/H] suffers considerable confusion, due to poor identification of the local thick disk. Unusual zig-zag abundance trends of [Cu/Fe] and [Na/Fe] suggest nucleosynthesis dominated by core-collapse supernovae with metallicity-dependent yields, in the Type~Ia supernova time-delay scenario. Thus, the bulge sub-population compositions resemble the local thin and thick disks, but at higher [Fe/H], suggesting a radial [Fe/H] gradient in both the thin and thick disks, within the solar circle, near -0.04 to -0.05 dex/kpc. If the bulge was built through accretion of inner disk stars by a bar, it appears that the inner thin and thick disk stars, incorporated into the bulge, retained vertical scale heights characteristic of their kinematic origin, and resulting in the vertical [Fe/H] gradient and [alpha/Fe] trends seen today.
We present the spatially resolved star formation history (SFH) of the Carina dwarf spheroidal galaxy, obtained from deep, wide-field g,r imaging and a metallicity distribution from the literature. Our photometry covers $\sim2$ deg$^2$, reaching up to $\sim10$ times the half-light radius of Carina with a completeness higher than $50\%$ at $g\sim24.5$, more than one magnitude fainter than the oldest turnoff. This is the first time a combination of depth and coverage of this quality has been used to derive the SFH of Carina, enabling us to trace its different populations with unprecedented accuracy. We find that Carina's SFH consists of two episodes well separated by a star formation temporal gap. These episodes occurred at old ($>10$ Gyr) and intermediate ($2$-$8$ Gyr) ages. Our measurements show that the old episode comprises the majority of the population, accounting for $54\pm5\%$ of the stellar mass within $1.3$ times the King tidal radius, while the total stellar mass derived for Carina is $1.60\pm0.09\times 10^{6} M_{\rm{\odot}}$, and the stellar mass-to-light ratio $1.8\pm0.2$. The SFH derived is consistent with no recent star formation which hints that the observed blue plume is due to blue stragglers. We conclude that the SFH of Carina evolved independently of the tidal field of the Milky Way, since the frequency and duration of its star formation events do not correlate with its orbital parameters. This result is supported by the age/metallicity relation observed in Carina, and the gradients calculated indicating that outer regions are older and more metal poor.
We study the behaviour of the dynamical and stellar mass inside the effective radius (re) of early-type galaxies (ETGs). We use several samples of ETGs -ranging from 19 000 to 98 000 objects- from the ninth data release of the Sloan Digital Sky Survey. We consider Newtonian dynamics, different light profiles and different Initial Mass Functions (IMF) to calculate the dynamical and stellar mass. We assume that any difference between these two masses is due to dark matter and/or a non Universal IMF. The main results for galaxies in the redshift range 0.0024 < z < 0.3500 and in the dynamical mass range 9.5 < log(M) < 12.5 are: i) A significant part of the intrinsic dispersion of the distribution of dynamical vs. stellar mass is due to redshift. ii) The difference between dynamical and stellar mass increases as a function of dynamical mass and decreases as a function of redshift. iii) The difference between dynamical and stellar mass goes from approximately 0% to 70% of the dynamical mass depending on mass and redshift. iv) These differences could be due to dark matter or a non Universal IMF or a combination of both. v) The amount of dark matter inside ETGs would be equal to or less than the difference between dynamical and stellar mass depending on the impact of the IMF on the stellar mass estimation. vi) The previous results go in the same direction of some results of the Fundamental Plane (FP) found in the literature in the sense that they could be interpreted as an increase of dark matter along the FP and a dependence of the FP on redshift.
In a previous study of the L1157 B1 shocked cavity, a comparison between NH$_3$(1$_0$-$0_0$) and H$_2$O(1$_{\rm 10}$--1$_{\rm 01}$) transitions showed a striking difference in the profiles, with H$_2$O emitting at definitely higher velocities. This behaviour was explained as a result of the high-temperature gas-phase chemistry occurring in the postshock gas in the B1 cavity of this outflow. If the differences in behaviour between ammonia and water are indeed a consequence of the high gas temperatures reached during the passage of a shock, then one should find such differences to be ubiquitous among chemically rich outflows. In order to determine whether the difference in profiles observed between NH$_3$ and H$_2$O is unique to L1157 or a common characteristic of chemically rich outflows, we have performed Herschel-HIFI observations of the NH$_3$(1$_0$-0$_0$) line at 572.5 GHz in a sample of 8 bright low-mass outflow spots already observed in the H$_2$O(1$_{\rm 10}$--1$_{\rm 01}$) line within the WISH KP. We detected the ammonia emission at high-velocities at most of the outflows positions. In all cases, the water emission reaches higher velocities than NH$_3$, proving that this behaviour is not exclusive of the L1157-B1 position. Comparisons with a gas-grain chemical and shock model confirms, for this larger sample, that the behaviour of ammonia is determined principally by the temperature of the gas.
The halo of the Milky-Way circumgalactic gas extends up to the virial radius of the Galaxy, ~250 kpc. The halo properties may be deduced from X-ray spectroscopic observations and from studies of the ram-pressure stripping of satellite dwarf galaxies. The former method is more precise but its results depend crucially on the assumed metallicity of the circumgalactic gas; the latter one does not need these assumptions. Here, the information from both approaches is combined to constrain observationally the gas metallicity and density as functions of the galactocentric distance. It is demonstrated that the metallicity decreases in the outer parts of the extended halo.
(abridged) The fairly recent detection of a variety of anions in the Interstellar Molecular Clouds have underlined the importance of realistically modeling the processes governing their abundance. To this aim, our earlier calculations for the radiative electron attachment (REA) rates for C4H-, C6H-, and C8H- are employed to generate the corresponding column density ratios of anion/neutral (A/N) relative abundances. The latter are then compared with those obtained from observational measurements. The calculations involved the time-dependent solutions of a large network of chemical processes over an extended time interval and included a series of runs in which the values of REA rates were repeatedly scaled. Macroscopic parameters for the clouds' modeling were also varied to cover a broad range of physical environments. It was found that, within the range and quality of the processes included in the present network,and selected from state-of-the-art astrophysical databases, the REA values required to match the observed A/N ratios needed to be reduced by orders of magnitude for C4H- case, while the same rates for C6H- and C8H- only needed to be scaled by much smaller factors. The results suggest that the generally proposed formation of interstellar anions by REA mechanism is overestimated by current models for the C4H- case, for which is likely to be an inefficient path to formation. This path is thus providing a rather marginal contribution to the observed abundances of C4H-, the latter being more likely to originate from other chemical processes in the network, as we discuss in some detail in the present work.Possible physical reasons for the much smaller differences against observations found instead for the values of the (A/N) ratios in two other, longer members of the series are put forward and analyzed within the evolutionary modeling discussed in the present work.
We use the Sloan Digital Sky Survey (SDSS) database to explore the effect of the disk inclination angle on the derived star formation rate (SFR), hence on the slope and width of the Main Sequence (MS) relation for star-forming galaxies. We find that SFRs for nearly edge-on disks are underestimated by factors ranging from $\sim$ 0.2 dex for low mass galaxies up to $\sim$ 0.4 dex for high mass galaxies. This results in a substantially flatter MS relation for high-inclination disks compared to that for less inclined ones, though the global effect over the whole sample of star-forming galaxies is relatively minor, given the small fraction of high-inclination disks. However, we also find that galaxies with high-inclination disks represent a non negligible fraction of galaxies populating the so-called green valley, with derived SFRs intermediate between the MS and those of quenched, passively evolving galaxies.
We develop models of two-component spherical galaxies to establish scaling relations linking the properties of spheroids at $z=0$ (total stellar masses, effective radii $R_e$ and velocity dispersions within $R_e$) to the properties of their dark-matter halos at both $z=0$ and higher redshifts. . Our main motivation is the widely accepted idea that the accretion-driven growth of supermassive black holes (SMBHs) in protogalaxies is limited by quasar-mode feedback and gas blow-out. The SMBH masses, $M_{\rm{BH}}$, should then be connected to the dark-matter potential wells at the redshift $z_{\rm{qso}}$ of the blow-out. We specifically consider the example of a power-law dependence on the maximum circular speed in a protogalactic dark-matter halo: $M_{\rm{BH}}\propto V^4_{\rm{d,pk}}$, as could be expected if quasar-mode feedback were momentum-driven. For halos with a given $V_{\rm{d,pk}}$ at a given $z_{\rm{qso}}\ge 0$, our model scaling relations give a typical stellar velocity dispersion $\sigma_{\rm{ap}}(R_e)$ at $z=0$. Thus, they transform a theoretical $M_{\rm{BH}}$-$V_{\rm{d,pk}}$ relation into a prediction for an observable $M_{\rm{BH}}$-$\sigma_{\rm{ap}}(R_e)$ relation. We find the latter to be distinctly non-linear in log-log space. Its shape depends on the generic redshift-evolution of halos in a {$\Lambda$}CDM cosmology and the systematic variation of stellar-to-dark matter mass fraction at $z=0$, in addition to any assumptions about the physics underlying the $M_{\rm{BH}}$-$V_{\rm{d,pk}}$ relation. Despite some clear limitations of the form we use for $M_{\rm{BH}}$ versus $V_{\rm{d,pk}}$, and even though we do not include any SMBH growth through dry mergers at low redshift, our results for $M_{\rm{BH}}$-$\sigma_{\rm{ap}}(R_e)$ compare well to data for local early types if we take $z_{\rm{qso}} \sim$ 2-4.
Molecular cloud structure is regulated by stellar feedback in various forms. Two of the most important feedback processes are UV photoionisation and supernovae from massive stars. However, the precise response of the cloud to these processes, and the interaction between them, remains an open question. In particular, we wish to know under which conditions the cloud can be dispersed by feedback, which in turn can give us hints as to how feedback regulates the star formation inside the cloud. We perform a suite of radiative magnetohydrodynamic simulations of a 10^5 solar mass cloud with embedded sources of ionising radiation and supernovae, including multiple supernovae and a hypernova model. A UV source corresponding to 10% of the mass of the cloud is required to disperse the cloud, suggesting that the star formation efficiency should be on the order of 10%. A single supernova is unable to significantly affect the evolution of the cloud. However, energetic hypernovae and multiple supernovae are able to add significant quantities of momentum to the cloud, approximately 10^{43} g cm/s of momentum per 10^{51} ergs of supernova energy. This is on the lower range of estimates in other works, since dense gas clumps that remain embedded inside the HII region cause rapid cooling in the supernova blast. We argue that supernovae alone are unable to regulate star formation in molecular clouds, and that strong pre-supernova feedback is required to allow supernova blastwaves to propagate efficiently into the interstellar medium
We present an analysis of the gas-phase oxygen abundances of a sample of 28 galaxies in the local Universe (z < 0.02) hosting Type Ia Supernovae (SNe Ia). The data were obtained with the 4.2m William Herschel Telescope (WHT). We derive local oxygen abundances for the regions where the SNe Ia exploded by calculating oxygen gradients through each galaxy (when possible) or assuming the oxygen abundance of the closest HII region. The sample selection only considered galaxies for which distances not based on the the SN Ia method are available. Then, we use a principal component analysis to study the dependence of the absolute magnitudes on the color of the SN Ia, the oxygen abundances of the region where they exploded, and the stretch of the SN light curve. We demonstrate that our previous result suggesting a metallicity-dependence on the SN Ia luminosity for not-reddened SNe Ia (Moreno-Raya et al. 2016) can be extended to our whole sample. These results reinforce the need of including a metallicity proxy, such as the oxygen abundance of the host galaxy, to minimize the systematic effect induced by the metallicity-dependence of the SN Ia luminosity in future studies of SNe Ia at cosmological distances.
We study the evolution of galactic magnetic fields using 3D smoothed particle magnetohydrodynamics (SPMHD) simulations of galaxies with an imposed spiral potential. We consider the appearance of reversals of the field, and amplification of the field. We find magnetic field reversals occur when the velocity jump across the spiral shock is above $\approx$20km s$^{-1}$, occurring where the velocity change is highest, typically at the inner Lindblad resonance (ILR) in our models. Reversals also occur at corotation, where the direction of the velocity field reverses in the co-rotating frame of a spiral arm. They occur earlier with a stronger amplitude spiral potential, and later or not at all with weaker or no spiral arms. The presence of a reversal at a radii of around 4--6 kpc in our fiducial model is consistent with a reversal identified in the Milky Way, though we caution that alternative Galaxy models could give a similar reversal. We find that relatively high resolution, a few million particles in SPMHD, is required to produce consistent behaviour of the magnetic field. Amplification of the magnetic field occurs in the models, and while some may be genuinely attributable to differential rotation or spiral arms, some may be a numerical artefact. We check our results using Athena, finding reversals but less amplification of the field, suggesting that some of the amplification of the field with SPMHD is numerical.
This paper describes an automatic isophotal fitting procedure that succeeds, without the support of any visual inspection of neither the images nor the ellipticity/P.A. radial profiles, at extracting a fairly pure sample of barred LTGs among thousands of optical images from the SDSS. The procedure relies on the methods described in Consolandi et al. (2016) to robustly extract the photometrical properties of a large sample of local SDSS galaxies and is tailored to extract bars on the basis of their well-known peculiarities in their P.A. and ellipticity profiles. It has been run on a sample of 5853 galaxies in the Coma and Local supercluster. The procedure extracted for each galaxy a color, an ellipticity and a position angle radial profile of the ellipses fitted to the isophotes. Examining automatically the profiles of 922 face-on late-type galaxies (B/A >0.7) the procedure found that ~ 36 % are barred. The local bar fraction strongly increases with stellar mass. The sample of barred galaxies is used to construct a set of template radial color profiles in order to test the impact of the barred galaxy population on the average color profiles shown by Consolandi et al. (2016) and to test the bar-quenching scenario proposed in Gavazzi et al. (2015). The radial color profile of barred galaxy shows that bars are on average redder than their surrounding disk producing an outside-in gradient toward red in correspondence of their corotation radius. The distribution of the extension of the deprojected length of the bar suggests that bars have strong impacts on the gradients of averaged color profiles. The dependence of the profiles on the mass is consistent with the bar-quenching scenario, i.e. more massive barred galaxies have redder colors (hence older stellar population and suppressed star formation) inside their corotation radius with respect to their lower mass counterparts.
We compare observed far infra-red/sub-millimetre (FIR/sub-mm) galaxy spectral energy distributions (SEDs) of massive galaxies ($M_{\star}\gtrsim10^{10}$ $h^{-1}$M$_{\odot}$) derived through a stacking analysis with predictions from a new model of galaxy formation. The FIR SEDs of the model galaxies are calculated using a self-consistent model for the absorption and re-emission of radiation by interstellar dust based on radiative transfer calculations and global energy balance arguments. Galaxies are selected based on their position on the specific star formation rate (sSFR) - stellar mass ($M_{\star}$) plane. We identify a main sequence of star-forming galaxies in the model, i.e. a well defined relationship between sSFR and $M_\star$, up to redshift $z\sim6$. The scatter of this relationship evolves such that it is generally larger at higher stellar masses and higher redshifts. There is remarkable agreement between the predicted and observed average SEDs across a broad range of redshifts ($0.5\lesssim z\lesssim4$) for galaxies on the main sequence. However, the agreement is less good for starburst galaxies at $z\gtrsim2$, selected here to have elevated sSFRs$>10\times$ the main sequence value. We find that the predicted average SEDs are robust to changing the parameters of our dust model within physically plausible values. We also show that the dust temperature evolution of main sequence galaxies in the model is driven by star formation on the main sequence being more burst-dominated at higher redshifts.
We construct a Bayesian framework to perform inference of dim or overlapping point sources. The method involves probabilistic cataloging, where samples are taken from the posterior probability distribution of catalogs consistent with an observed photon count map. By implementing across-model jumps between point source models of different dimensionality, we collect a representative ensemble of catalogs consistent with the image. In order to validate our method we sample random catalogs of the gamma-ray sky in the direction of the North Galactic Pole (NGP) by binning the data in energy and PSF (Point Spread Function) classes. Using three energy bins between $0.3 - 1$, $1 - 3$ and $3 - 10$ GeV, we identify $270\substack{+30 -10}$ point sources inside a $40^\circ \times 40^\circ$ region around the NGP above our point-source inclusion limit of $3 \times 10^{-11}$/cm$^2$/s/sr/GeV at the $1-3$ GeV energy bin. Most of these point sources are time-variable blazars. Modeling the flux distribution as a single power law, we infer the slope to be $-1.92\substack{+0.07 -0.05}$ and estimate the contribution of point sources (resolved and unresolved) to the total emission as $18\substack{+2 -2}$\%. Further analyses that rely on the ensemble of sample catalogs instead of only the most likely catalog, can perform reliable marginalization over uncertainties in the number as well as spatial and spectral properties of the point sources. This marginalization allows a robust test of whether the apparently isotropic emission in an image is due to unresolved point sources or of truly diffuse origin. With the increase in the availability of computational resources in the near future, probabilistic cataloging can potentially be applied to full sky datasets or optical images and replace the standard data reduction pipelines for crowded fields.
We present a systematic study of the intensity mapping technique using updated models for the different emission lines from galaxies and identify which ones are more promising for cosmological studies of the post reionization epoch. We consider the emission of ${\rm Ly\alpha}$, ${\rm H\alpha}$, H$\beta$, optical and infrared oxygen lines, nitrogen lines, CII and the CO rotational lines. We then identify that ${\rm Ly\alpha}$, ${\rm H\alpha}$, OII, CII and the lowest rotational CO lines are the best candidates to be used as IM probes. These lines form a complementary set of probes of the galaxies emission spectra. We then use reasonable experimental setups from current, planned or proposed experiments to access the detectability of the power spectrum of each emission line. Intensity mapping of ${\rm Ly\alpha}$ emission from $z=2$ to 3 will be possible in the near future with HETDEX, while far-infrared lines require new dedicated experiments. We also show that the proposed SPHEREx satellite can use OII and ${\rm H\alpha}$ IM to study the large-scale distribution of matter in intermediate redshifts of 1 to 4. We found that submilimeter experiments with bolometers can have similar performances at intermediate redshifts using CII and CO(3-2).
We present the first scientific results from the luminous red galaxy sample (LRG) of the extended Baryon Oscillation Spectroscopic Survey (eBOSS). We measure the small and intermediate scale clustering from a sample of more than 61,000 galaxies in the redshift range $0.6 < z < 0.9$. We interpret these measurements in the framework of the Halo Occupation Distribution. The bias of eBOSS LRGs is $2.30 \pm 0.03$, with a satellite fraction of $13\pm3$\% and a mean halo mass of $2.5\times10^{13}h^{-1}M_{\odot}$. These results are consistent with expectations, demonstrating that eBOSS galaxies will be reliable tracers of large scale structure at $z\sim 0.7$. The eBOSS galaxy bias implies a scatter of luminosity at fixed halo mass, $\sigma_{\log L}$, of 0.19 dex. Using the clustering of massive galaxies from BOSS-CMASS, BOSS-LOWZ, and SDSS, we find that $\sigma_{\log L}=0.19$ is consistent with observations over the full redshift range that these samples cover. The addition of eBOSS to previous surveys allows investigation of the evolution of massive galaxies over the past $\sim 7$ Gyr.
We present an observational estimate of the fraction and distribution of dark mass in the innermost region of the two Galactic globular clusters NGC 6218 (M12) and NGC 288. Such an assessment has been made by comparing the dynamical and luminous mass profiles derived from an accurate analysis of the most extensive spectroscopic and photometric surveys performed on these stellar systems. We find that non-luminous matter constitutes more than 60% of the total mass in the region probed by our data (R<1.6 arcmin~r_h) in both clusters. We have carefully analyzed the effects of binaries and tidal heating on our estimate and ruled out the possibility that our result is a spurious consequence of these effects. The dark component appears to be more concentrated than the most massive stars suggesting that it is likely composed of dark remnants segregated in the cluster core.
Links to: arXiv, form interface, find, astro-ph, recent, 1607, contact, help (Access key information)
Using high-resolution simulations with explicit treatment of stellar feedback physics based on the FIRE (Feedback in Realistic Environments) project, we study how galaxy formation and the interstellar medium (ISM) are affected by magnetic fields, anisotropic Spitzer-Braginskii conduction and viscosity, and sub-grid turbulent metal diffusion. We consider controlled simulations of isolated (non-cosmological) galaxies but also a limited set of cosmological "zoom-in" simulations. Although simulations have shown significant effects from these physics with weak or absent stellar feedback, the effects are much weaker than those of stellar feedback when the latter is modeled explicitly. The additional physics have no systematic effect on galactic star formation rates (SFRs). In contrast, removing stellar feedback leads to SFRs being over-predicted by factors of $\sim 10 -100$. Without feedback, neither galactic winds nor volume filling hot-phase gas exist, and discs tend to runaway collapse to ultra-thin scale-heights with unphysically dense clumps congregating at the galactic center. With stellar feedback, a multi-phase, turbulent medium with galactic fountains and winds is established. At currently achievable resolutions, the additional physics investigated here (MHD, conduction, viscosity, metal diffusion) have only weak ($\sim10\%$-level) effects on these properties and do not significantly alter the balance of phases, outflows, or the energy in ISM turbulence, consistent with simple equipartition arguments. We conclude that galactic star formation and the ISM are primarily governed by a combination of turbulence, gravitational instabilities, and feedback.
The exponential scale length ($L_d$) of the Milky Way's (MW's) disk is a critical parameter for describing the global physical size of our Galaxy, important both for interpreting other Galactic measurements and helping us to understand how our Galaxy fits into extragalactic contexts. Unfortunately, current estimates span a wide range of values and often are statistically incompatible with one another. Here, we aim to determine an improved, aggregate estimate for $L_d$ by utilizing a hierarchical Bayesian (HB) meta-analysis technique that accounts for the possibility that any one measurement has not properly accounted for all statistical or systematic errors. Within this machinery we explore a variety of ways of modeling the nature of problematic measurements, and then use a Bayesian model averaging technique to derive net posterior distributions that incorporate any model-selection uncertainty. Our meta-analysis combines 29 different (15 visible and 14 infrared) photometric measurements of $L_d$ available in the literature; these involve a broad assortment of observational datasets, MW models and assumptions, and methodologies, all tabulated herein. Analyzing the visible and infrared measurements separately yields estimates for $L_d$ of $2.71^{+0.22}_{-0.20}$ kpc and $2.51^{+0.15}_{-0.13}$ kpc, respectively, whereas considering them all combined yields $2.64\pm0.13$ kpc. The ratio between the visible and infrared scale lengths determined here is very similar to that measured in external spiral galaxies. We use these results to updated the model of the Galactic disk from our previous work, constraining its stellar mass to be $4.8^{+1.5}_{-1.1}\times10^{10}$ $\mathrm{M}_\odot$, corresponding to a total stellar mass of $5.7^{+1.5}_{-1.1}\times10^{10}$ $\mathrm{M}_\odot$.
We use spectra from the ALFALFA, GASS and COLD GASS surveys to quantify variations in the mean atomic and molecular gas mass fractions throughout the SFR-M* plane and along the main sequence (MS) of star-forming galaxies. Although galaxies well below the MS tend to be undetected in the Arecibo and IRAM observations, reliable mean atomic and molecular gas fractions can be obtained through a spectral stacking technique. We find that the position of galaxies in the SFR-M* plane can be explained mostly by their global cold gas reservoirs as observed in the HI line, with in addition systematic variations in the molecular-to-atomic ratio and star formation efficiency. When looking at galaxies within +/-0.4 dex of the MS, we find that as stellar mass increases, both atomic and molecular gas mass fractions decrease, stellar bulges become more prominent, and the mean stellar ages increase. Both star formation efficiency and molecular-to-atomic ratios vary little for massive main sequence galaxies, indicating that the flattening of the MS is due to the global decrease of the cold gas reservoirs of galaxies rather than to bottlenecks in the process of converting cold atomic gas to stars.
The average bulge [Fe/H] and [Mg/H] are +0.06 and +0.17 dex, respectively, in Baade's Window, roughly 0.2 dex higher than the thin disk and ~0.7 dex higher than the local thick disk metallicity. This suggests a higher effective yield in the bulge, perhaps due to more efficient retention of supernova ejecta. The bulge vertical [Fe/H] gradient, at ~0.5 dex/kpc, appears to be due to a changing mixture of sub-populations (near +0.3 dex and -0.3 dex and one possibly near -0.7 dex) with latitude. The bulge is enhanced in O, Mg, Si, Ca, Ti, and Al relative to the sun, with [alpha/Fe]=+0.15 dex at [Fe/H]=0.0 dex. Below [Fe/H]~-0.5 dex, the bulge and local thick disk compositions are very similar, but small [Mg/Fe] and possibly [<SiCaTi>/Fe] enhancements, low [La/Eu] ratios and large [Cu/Fe], relative to the thick disk suggest slightly higher SFR in the bulge. However, these composition differences could simply be due to measurement errors and non-LTE effects. Unfortunately, comparison with the thick disk near solar [Fe/H] suffers considerable confusion, due to poor identification of the local thick disk. Unusual zig-zag abundance trends of [Cu/Fe] and [Na/Fe] suggest nucleosynthesis dominated by core-collapse supernovae with metallicity-dependent yields, in the Type~Ia supernova time-delay scenario. Thus, the bulge sub-population compositions resemble the local thin and thick disks, but at higher [Fe/H], suggesting a radial [Fe/H] gradient in both the thin and thick disks, within the solar circle, near -0.04 to -0.05 dex/kpc. If the bulge was built through accretion of inner disk stars by a bar, it appears that the inner thin and thick disk stars, incorporated into the bulge, retained vertical scale heights characteristic of their kinematic origin, and resulting in the vertical [Fe/H] gradient and [alpha/Fe] trends seen today.
We present the spatially resolved star formation history (SFH) of the Carina dwarf spheroidal galaxy, obtained from deep, wide-field g,r imaging and a metallicity distribution from the literature. Our photometry covers $\sim2$ deg$^2$, reaching up to $\sim10$ times the half-light radius of Carina with a completeness higher than $50\%$ at $g\sim24.5$, more than one magnitude fainter than the oldest turnoff. This is the first time a combination of depth and coverage of this quality has been used to derive the SFH of Carina, enabling us to trace its different populations with unprecedented accuracy. We find that Carina's SFH consists of two episodes well separated by a star formation temporal gap. These episodes occurred at old ($>10$ Gyr) and intermediate ($2$-$8$ Gyr) ages. Our measurements show that the old episode comprises the majority of the population, accounting for $54\pm5\%$ of the stellar mass within $1.3$ times the King tidal radius, while the total stellar mass derived for Carina is $1.60\pm0.09\times 10^{6} M_{\rm{\odot}}$, and the stellar mass-to-light ratio $1.8\pm0.2$. The SFH derived is consistent with no recent star formation which hints that the observed blue plume is due to blue stragglers. We conclude that the SFH of Carina evolved independently of the tidal field of the Milky Way, since the frequency and duration of its star formation events do not correlate with its orbital parameters. This result is supported by the age/metallicity relation observed in Carina, and the gradients calculated indicating that outer regions are older and more metal poor.
We study the behaviour of the dynamical and stellar mass inside the effective radius (re) of early-type galaxies (ETGs). We use several samples of ETGs -ranging from 19 000 to 98 000 objects- from the ninth data release of the Sloan Digital Sky Survey. We consider Newtonian dynamics, different light profiles and different Initial Mass Functions (IMF) to calculate the dynamical and stellar mass. We assume that any difference between these two masses is due to dark matter and/or a non Universal IMF. The main results for galaxies in the redshift range 0.0024 < z < 0.3500 and in the dynamical mass range 9.5 < log(M) < 12.5 are: i) A significant part of the intrinsic dispersion of the distribution of dynamical vs. stellar mass is due to redshift. ii) The difference between dynamical and stellar mass increases as a function of dynamical mass and decreases as a function of redshift. iii) The difference between dynamical and stellar mass goes from approximately 0% to 70% of the dynamical mass depending on mass and redshift. iv) These differences could be due to dark matter or a non Universal IMF or a combination of both. v) The amount of dark matter inside ETGs would be equal to or less than the difference between dynamical and stellar mass depending on the impact of the IMF on the stellar mass estimation. vi) The previous results go in the same direction of some results of the Fundamental Plane (FP) found in the literature in the sense that they could be interpreted as an increase of dark matter along the FP and a dependence of the FP on redshift.
In a previous study of the L1157 B1 shocked cavity, a comparison between NH$_3$(1$_0$-$0_0$) and H$_2$O(1$_{\rm 10}$--1$_{\rm 01}$) transitions showed a striking difference in the profiles, with H$_2$O emitting at definitely higher velocities. This behaviour was explained as a result of the high-temperature gas-phase chemistry occurring in the postshock gas in the B1 cavity of this outflow. If the differences in behaviour between ammonia and water are indeed a consequence of the high gas temperatures reached during the passage of a shock, then one should find such differences to be ubiquitous among chemically rich outflows. In order to determine whether the difference in profiles observed between NH$_3$ and H$_2$O is unique to L1157 or a common characteristic of chemically rich outflows, we have performed Herschel-HIFI observations of the NH$_3$(1$_0$-0$_0$) line at 572.5 GHz in a sample of 8 bright low-mass outflow spots already observed in the H$_2$O(1$_{\rm 10}$--1$_{\rm 01}$) line within the WISH KP. We detected the ammonia emission at high-velocities at most of the outflows positions. In all cases, the water emission reaches higher velocities than NH$_3$, proving that this behaviour is not exclusive of the L1157-B1 position. Comparisons with a gas-grain chemical and shock model confirms, for this larger sample, that the behaviour of ammonia is determined principally by the temperature of the gas.
The halo of the Milky-Way circumgalactic gas extends up to the virial radius of the Galaxy, ~250 kpc. The halo properties may be deduced from X-ray spectroscopic observations and from studies of the ram-pressure stripping of satellite dwarf galaxies. The former method is more precise but its results depend crucially on the assumed metallicity of the circumgalactic gas; the latter one does not need these assumptions. Here, the information from both approaches is combined to constrain observationally the gas metallicity and density as functions of the galactocentric distance. It is demonstrated that the metallicity decreases in the outer parts of the extended halo.
(abridged) The fairly recent detection of a variety of anions in the Interstellar Molecular Clouds have underlined the importance of realistically modeling the processes governing their abundance. To this aim, our earlier calculations for the radiative electron attachment (REA) rates for C4H-, C6H-, and C8H- are employed to generate the corresponding column density ratios of anion/neutral (A/N) relative abundances. The latter are then compared with those obtained from observational measurements. The calculations involved the time-dependent solutions of a large network of chemical processes over an extended time interval and included a series of runs in which the values of REA rates were repeatedly scaled. Macroscopic parameters for the clouds' modeling were also varied to cover a broad range of physical environments. It was found that, within the range and quality of the processes included in the present network,and selected from state-of-the-art astrophysical databases, the REA values required to match the observed A/N ratios needed to be reduced by orders of magnitude for C4H- case, while the same rates for C6H- and C8H- only needed to be scaled by much smaller factors. The results suggest that the generally proposed formation of interstellar anions by REA mechanism is overestimated by current models for the C4H- case, for which is likely to be an inefficient path to formation. This path is thus providing a rather marginal contribution to the observed abundances of C4H-, the latter being more likely to originate from other chemical processes in the network, as we discuss in some detail in the present work.Possible physical reasons for the much smaller differences against observations found instead for the values of the (A/N) ratios in two other, longer members of the series are put forward and analyzed within the evolutionary modeling discussed in the present work.
We use the Sloan Digital Sky Survey (SDSS) database to explore the effect of the disk inclination angle on the derived star formation rate (SFR), hence on the slope and width of the Main Sequence (MS) relation for star-forming galaxies. We find that SFRs for nearly edge-on disks are underestimated by factors ranging from $\sim$ 0.2 dex for low mass galaxies up to $\sim$ 0.4 dex for high mass galaxies. This results in a substantially flatter MS relation for high-inclination disks compared to that for less inclined ones, though the global effect over the whole sample of star-forming galaxies is relatively minor, given the small fraction of high-inclination disks. However, we also find that galaxies with high-inclination disks represent a non negligible fraction of galaxies populating the so-called green valley, with derived SFRs intermediate between the MS and those of quenched, passively evolving galaxies.
We develop models of two-component spherical galaxies to establish scaling relations linking the properties of spheroids at $z=0$ (total stellar masses, effective radii $R_e$ and velocity dispersions within $R_e$) to the properties of their dark-matter halos at both $z=0$ and higher redshifts. . Our main motivation is the widely accepted idea that the accretion-driven growth of supermassive black holes (SMBHs) in protogalaxies is limited by quasar-mode feedback and gas blow-out. The SMBH masses, $M_{\rm{BH}}$, should then be connected to the dark-matter potential wells at the redshift $z_{\rm{qso}}$ of the blow-out. We specifically consider the example of a power-law dependence on the maximum circular speed in a protogalactic dark-matter halo: $M_{\rm{BH}}\propto V^4_{\rm{d,pk}}$, as could be expected if quasar-mode feedback were momentum-driven. For halos with a given $V_{\rm{d,pk}}$ at a given $z_{\rm{qso}}\ge 0$, our model scaling relations give a typical stellar velocity dispersion $\sigma_{\rm{ap}}(R_e)$ at $z=0$. Thus, they transform a theoretical $M_{\rm{BH}}$-$V_{\rm{d,pk}}$ relation into a prediction for an observable $M_{\rm{BH}}$-$\sigma_{\rm{ap}}(R_e)$ relation. We find the latter to be distinctly non-linear in log-log space. Its shape depends on the generic redshift-evolution of halos in a {$\Lambda$}CDM cosmology and the systematic variation of stellar-to-dark matter mass fraction at $z=0$, in addition to any assumptions about the physics underlying the $M_{\rm{BH}}$-$V_{\rm{d,pk}}$ relation. Despite some clear limitations of the form we use for $M_{\rm{BH}}$ versus $V_{\rm{d,pk}}$, and even though we do not include any SMBH growth through dry mergers at low redshift, our results for $M_{\rm{BH}}$-$\sigma_{\rm{ap}}(R_e)$ compare well to data for local early types if we take $z_{\rm{qso}} \sim$ 2-4.
Molecular cloud structure is regulated by stellar feedback in various forms. Two of the most important feedback processes are UV photoionisation and supernovae from massive stars. However, the precise response of the cloud to these processes, and the interaction between them, remains an open question. In particular, we wish to know under which conditions the cloud can be dispersed by feedback, which in turn can give us hints as to how feedback regulates the star formation inside the cloud. We perform a suite of radiative magnetohydrodynamic simulations of a 10^5 solar mass cloud with embedded sources of ionising radiation and supernovae, including multiple supernovae and a hypernova model. A UV source corresponding to 10% of the mass of the cloud is required to disperse the cloud, suggesting that the star formation efficiency should be on the order of 10%. A single supernova is unable to significantly affect the evolution of the cloud. However, energetic hypernovae and multiple supernovae are able to add significant quantities of momentum to the cloud, approximately 10^{43} g cm/s of momentum per 10^{51} ergs of supernova energy. This is on the lower range of estimates in other works, since dense gas clumps that remain embedded inside the HII region cause rapid cooling in the supernova blast. We argue that supernovae alone are unable to regulate star formation in molecular clouds, and that strong pre-supernova feedback is required to allow supernova blastwaves to propagate efficiently into the interstellar medium
We present an analysis of the gas-phase oxygen abundances of a sample of 28 galaxies in the local Universe (z < 0.02) hosting Type Ia Supernovae (SNe Ia). The data were obtained with the 4.2m William Herschel Telescope (WHT). We derive local oxygen abundances for the regions where the SNe Ia exploded by calculating oxygen gradients through each galaxy (when possible) or assuming the oxygen abundance of the closest HII region. The sample selection only considered galaxies for which distances not based on the the SN Ia method are available. Then, we use a principal component analysis to study the dependence of the absolute magnitudes on the color of the SN Ia, the oxygen abundances of the region where they exploded, and the stretch of the SN light curve. We demonstrate that our previous result suggesting a metallicity-dependence on the SN Ia luminosity for not-reddened SNe Ia (Moreno-Raya et al. 2016) can be extended to our whole sample. These results reinforce the need of including a metallicity proxy, such as the oxygen abundance of the host galaxy, to minimize the systematic effect induced by the metallicity-dependence of the SN Ia luminosity in future studies of SNe Ia at cosmological distances.
We study the evolution of galactic magnetic fields using 3D smoothed particle magnetohydrodynamics (SPMHD) simulations of galaxies with an imposed spiral potential. We consider the appearance of reversals of the field, and amplification of the field. We find magnetic field reversals occur when the velocity jump across the spiral shock is above $\approx$20km s$^{-1}$, occurring where the velocity change is highest, typically at the inner Lindblad resonance (ILR) in our models. Reversals also occur at corotation, where the direction of the velocity field reverses in the co-rotating frame of a spiral arm. They occur earlier with a stronger amplitude spiral potential, and later or not at all with weaker or no spiral arms. The presence of a reversal at a radii of around 4--6 kpc in our fiducial model is consistent with a reversal identified in the Milky Way, though we caution that alternative Galaxy models could give a similar reversal. We find that relatively high resolution, a few million particles in SPMHD, is required to produce consistent behaviour of the magnetic field. Amplification of the magnetic field occurs in the models, and while some may be genuinely attributable to differential rotation or spiral arms, some may be a numerical artefact. We check our results using Athena, finding reversals but less amplification of the field, suggesting that some of the amplification of the field with SPMHD is numerical.
This paper describes an automatic isophotal fitting procedure that succeeds, without the support of any visual inspection of neither the images nor the ellipticity/P.A. radial profiles, at extracting a fairly pure sample of barred LTGs among thousands of optical images from the SDSS. The procedure relies on the methods described in Consolandi et al. (2016) to robustly extract the photometrical properties of a large sample of local SDSS galaxies and is tailored to extract bars on the basis of their well-known peculiarities in their P.A. and ellipticity profiles. It has been run on a sample of 5853 galaxies in the Coma and Local supercluster. The procedure extracted for each galaxy a color, an ellipticity and a position angle radial profile of the ellipses fitted to the isophotes. Examining automatically the profiles of 922 face-on late-type galaxies (B/A >0.7) the procedure found that ~ 36 % are barred. The local bar fraction strongly increases with stellar mass. The sample of barred galaxies is used to construct a set of template radial color profiles in order to test the impact of the barred galaxy population on the average color profiles shown by Consolandi et al. (2016) and to test the bar-quenching scenario proposed in Gavazzi et al. (2015). The radial color profile of barred galaxy shows that bars are on average redder than their surrounding disk producing an outside-in gradient toward red in correspondence of their corotation radius. The distribution of the extension of the deprojected length of the bar suggests that bars have strong impacts on the gradients of averaged color profiles. The dependence of the profiles on the mass is consistent with the bar-quenching scenario, i.e. more massive barred galaxies have redder colors (hence older stellar population and suppressed star formation) inside their corotation radius with respect to their lower mass counterparts.
We compare observed far infra-red/sub-millimetre (FIR/sub-mm) galaxy spectral energy distributions (SEDs) of massive galaxies ($M_{\star}\gtrsim10^{10}$ $h^{-1}$M$_{\odot}$) derived through a stacking analysis with predictions from a new model of galaxy formation. The FIR SEDs of the model galaxies are calculated using a self-consistent model for the absorption and re-emission of radiation by interstellar dust based on radiative transfer calculations and global energy balance arguments. Galaxies are selected based on their position on the specific star formation rate (sSFR) - stellar mass ($M_{\star}$) plane. We identify a main sequence of star-forming galaxies in the model, i.e. a well defined relationship between sSFR and $M_\star$, up to redshift $z\sim6$. The scatter of this relationship evolves such that it is generally larger at higher stellar masses and higher redshifts. There is remarkable agreement between the predicted and observed average SEDs across a broad range of redshifts ($0.5\lesssim z\lesssim4$) for galaxies on the main sequence. However, the agreement is less good for starburst galaxies at $z\gtrsim2$, selected here to have elevated sSFRs$>10\times$ the main sequence value. We find that the predicted average SEDs are robust to changing the parameters of our dust model within physically plausible values. We also show that the dust temperature evolution of main sequence galaxies in the model is driven by star formation on the main sequence being more burst-dominated at higher redshifts.
We construct a Bayesian framework to perform inference of dim or overlapping point sources. The method involves probabilistic cataloging, where samples are taken from the posterior probability distribution of catalogs consistent with an observed photon count map. By implementing across-model jumps between point source models of different dimensionality, we collect a representative ensemble of catalogs consistent with the image. In order to validate our method we sample random catalogs of the gamma-ray sky in the direction of the North Galactic Pole (NGP) by binning the data in energy and PSF (Point Spread Function) classes. Using three energy bins between $0.3 - 1$, $1 - 3$ and $3 - 10$ GeV, we identify $270\substack{+30 -10}$ point sources inside a $40^\circ \times 40^\circ$ region around the NGP above our point-source inclusion limit of $3 \times 10^{-11}$/cm$^2$/s/sr/GeV at the $1-3$ GeV energy bin. Most of these point sources are time-variable blazars. Modeling the flux distribution as a single power law, we infer the slope to be $-1.92\substack{+0.07 -0.05}$ and estimate the contribution of point sources (resolved and unresolved) to the total emission as $18\substack{+2 -2}$\%. Further analyses that rely on the ensemble of sample catalogs instead of only the most likely catalog, can perform reliable marginalization over uncertainties in the number as well as spatial and spectral properties of the point sources. This marginalization allows a robust test of whether the apparently isotropic emission in an image is due to unresolved point sources or of truly diffuse origin. With the increase in the availability of computational resources in the near future, probabilistic cataloging can potentially be applied to full sky datasets or optical images and replace the standard data reduction pipelines for crowded fields.
We present a systematic study of the intensity mapping technique using updated models for the different emission lines from galaxies and identify which ones are more promising for cosmological studies of the post reionization epoch. We consider the emission of ${\rm Ly\alpha}$, ${\rm H\alpha}$, H$\beta$, optical and infrared oxygen lines, nitrogen lines, CII and the CO rotational lines. We then identify that ${\rm Ly\alpha}$, ${\rm H\alpha}$, OII, CII and the lowest rotational CO lines are the best candidates to be used as IM probes. These lines form a complementary set of probes of the galaxies emission spectra. We then use reasonable experimental setups from current, planned or proposed experiments to access the detectability of the power spectrum of each emission line. Intensity mapping of ${\rm Ly\alpha}$ emission from $z=2$ to 3 will be possible in the near future with HETDEX, while far-infrared lines require new dedicated experiments. We also show that the proposed SPHEREx satellite can use OII and ${\rm H\alpha}$ IM to study the large-scale distribution of matter in intermediate redshifts of 1 to 4. We found that submilimeter experiments with bolometers can have similar performances at intermediate redshifts using CII and CO(3-2).
We present the first scientific results from the luminous red galaxy sample (LRG) of the extended Baryon Oscillation Spectroscopic Survey (eBOSS). We measure the small and intermediate scale clustering from a sample of more than 61,000 galaxies in the redshift range $0.6 < z < 0.9$. We interpret these measurements in the framework of the Halo Occupation Distribution. The bias of eBOSS LRGs is $2.30 \pm 0.03$, with a satellite fraction of $13\pm3$\% and a mean halo mass of $2.5\times10^{13}h^{-1}M_{\odot}$. These results are consistent with expectations, demonstrating that eBOSS galaxies will be reliable tracers of large scale structure at $z\sim 0.7$. The eBOSS galaxy bias implies a scatter of luminosity at fixed halo mass, $\sigma_{\log L}$, of 0.19 dex. Using the clustering of massive galaxies from BOSS-CMASS, BOSS-LOWZ, and SDSS, we find that $\sigma_{\log L}=0.19$ is consistent with observations over the full redshift range that these samples cover. The addition of eBOSS to previous surveys allows investigation of the evolution of massive galaxies over the past $\sim 7$ Gyr.
We present an observational estimate of the fraction and distribution of dark mass in the innermost region of the two Galactic globular clusters NGC 6218 (M12) and NGC 288. Such an assessment has been made by comparing the dynamical and luminous mass profiles derived from an accurate analysis of the most extensive spectroscopic and photometric surveys performed on these stellar systems. We find that non-luminous matter constitutes more than 60% of the total mass in the region probed by our data (R<1.6 arcmin~r_h) in both clusters. We have carefully analyzed the effects of binaries and tidal heating on our estimate and ruled out the possibility that our result is a spurious consequence of these effects. The dark component appears to be more concentrated than the most massive stars suggesting that it is likely composed of dark remnants segregated in the cluster core.
Links to: arXiv, form interface, find, astro-ph, recent, 1607, contact, help (Access key information)
We present a new model of the nebular emission from star-forming galaxies in a wide range of chemical compositions, appropriate to interpret observations of galaxies at all cosmic epochs. The model relies on the combination of state-of-the-art stellar population synthesis and photoionization codes to describe the ensemble of HII regions and the diffuse gas ionized by young stars in a galaxy. A main feature of this model is the self-consistent yet versatile treatment of element abundances and depletion onto dust grains, which allows one to relate the observed nebular emission from a galaxy to both gas-phase and dust-phase metal enrichment. We show that this model can account for the rest-frame ultraviolet and optical emission-line properties of galaxies at different redshifts and find that ultraviolet emission lines are more sensitive than optical ones to parameters such as C/O abundance ratio, hydrogen gas density, dust-to-metal mass ratio and upper cutoff of the stellar initial mass function. We also find that, for gas-phase metallicities around solar to slightly sub-solar, widely used formulae to constrain oxygen ionic fractions and the C/O ratio from ultraviolet and optical emission-line luminosities are reasonable faithful. However, the recipes break down at non-solar metallicities, making them inappropriate to study chemically young galaxies. In such cases, a fully self-consistent model of the kind presented in this paper is required to interpret the observed nebular emission.
We present a panoramic map of the entire Milky Way halo north of dec~-30 degrees (~30,000 deg^2), constructed by applying the matched-filter technique to the Pan-STARRS1 3Pi Survey dataset. Using single-epoch photometry reaching to g~22, we are sensitive to stellar substructures with heliocentric distances between 3.5 and ~35 kpc. We recover almost all previously-reported streams in this volume and demonstrate that several of these are significantly more extended than earlier datasets have indicated. In addition, we also report five new candidate stellar streams. One of these features appears significantly broader and more luminous than the others and is likely the remnant of a dwarf galaxy. The other four streams are consistent with a globular cluster origin, and three of these are rather short in projection (<10 degrees), suggesting that streams like Ophiuchus may not be that rare. Finally, a significant number of more marginal substructures are also revealed by our analysis; many of these features can also be discerned in matched-filter maps produced by other authors from SDSS data, and hence they are very likely to be genuine. However, the extant 3Pi data is currently too shallow to determine their properties or produce convincing CMDs. The global view of the Milky Way provided by Pan-STARRS1 provides further evidence for the important role of both globular cluster disruption and dwarf galaxy accretion in building the Milky Way's stellar halo.
We report the detection of morphology-dependent stellar age in massive quenched galaxies at redshift z~1.2. The sense of the dependence is that compact quenched galaxies are 0.5-2 Gyr older than normal-sized ones. The evidence comes from three different age indicators, Dn4000, H$\delta_{A}$, and fits to spectral synthesis models, applied to the stacked optical spectra of massive, quenched galaxies observed during the CANDELS project. All age indicators consistently show that the dominant stellar populations of compact passive galaxies are older than those of their normally-sized counterparts, with the Dn4000 and model fitting in excellent quantitative agreement. We detect weak [OII] emission in a fraction of passive galaxies, and the strength of the line, when present, is similar between the two samples; however, compact galaxies have a significantly lower frequency of [OII] emission than normal ones. A fraction of both compact and normal galaxies are also individually detected in the 4 Ms Chandra X-ray images, with luminosity in the range 10$^{40}$-10$^{41}$ erg/sec; while both [OII] emitters and non-emitters are found to be X-ray sources among normal galaxies, no compact galaxy with [OII] emission is an X-ray source, arguing against an AGN powering the line in compact galaxies. We interpret the [OII] properties as further evidence that compact galaxies are older and are further along into the process of quenching star formation and suppressing gas accretion. Finally, we argue that the older stellar age of the compact passive galaxies is further evidence of progenitor bias: compact passive galaxies are older because they reflect the smaller sizes of galaxies at their earlier quenching epoch, with stellar density most likely having nothing directly to do with the cessation of their star formation.
Using the Sloan Digital Sky Survey, we adopt the sSFR-$\Sigma_{1kpc}$ diagram as a diagnostic tool to understand the nature of quenching in different environments. sSFR is the specific star formation rate, and $\Sigma_{1kpc}$ is the stellar surface density in the inner kpc. Although both the host halo mass and group-centric distance affect the satellite population, we find that these two properties can be characterised by a single number, the quenched fraction, such that key features of the sSFR-$\Sigma_{1kpc}$ diagram vary smoothly with this proxy for the "environment". Particularly, the sSFR of star-forming galaxies decreases smoothly with the quenched fraction of a given environment. Furthermore, the location of the transition galaxies (i.e., the "green valley" or GV) in the sSFR-$\Sigma_{1kpc}$ diagram also varies smoothly with the environment, $\Sigma_{1kpc}$ being lower for satellites than the field, and lower for satellites in larger halos and at smaller radial distances within the same-mass halos. We interpret this shift as indicating the relative importance in different environments of today's field quenching track vs. the cluster quenching track. These environmental effects in the sSFR-$\Sigma_{1kpc}$ diagram are most significant in our lowest mass range ($9.75 < \log M_{*}/M_{\odot} < 10$). One feature of the sSFR-$\Sigma_{1kpc}$ diagram that is shared between all environments is that at a given $M_{*}$ quenched galaxies have higher $\Sigma_{1kpc}$ than the star-forming population. These results disfavour scenarios that quench satellites without any (subsequent) change in $\Sigma_{1kpc}$ or $M_{*}$. We discuss possible scenarios to explain the difference in $\Sigma_{1kpc}$ between GV and quenched satellites.
We use the scatter in the stellar-to-halo mass relation to constrain galaxy evolution models. If the efficiency of converting accreted baryons into stars varies with time, halos of the same present-day mass but different formation histories will have different z=0 galaxy stellar mass. This is one of the sources of scatter in stellar mass at fixed halo mass, $\sigma_{\log M\ast}$. For massive halos that undergo rapid quenching of star formation at z~2, different mechanisms that trigger this quenching yield different values of $\sigma_{\log M\ast}$. We use this framework to test various models in which quenching begins after a galaxy crosses a threshold in one of the following physical quantities: redshift, halo mass, stellar mass, and stellar-to-halo mass ratio. Our model is highly idealized, with other sources of scatter likely to arise as more physics is included. Thus, our test is whether a model can produce scatter lower than observational bounds, leaving room for other sources. Recent measurements find $\sigma_{\log M\ast}=0.16$ dex for 10^11 Msol galaxies. Under the assumption that the threshold is constant with time, such a low value of $\sigma_{\log M\ast}$ rules out all of these models with the exception of quenching by a stellar mass treshold. Most physical quantities, such as metallicity, will increase scatter if they are uncorrelated with halo formation history. Thus, to decrease the scatter of a given model, galaxy properties would correlate tightly with formation history, creating testable predictions for their clustering. Understanding why $\sigma_{\log M\ast}$ is so small may be key to understanding the physics of galaxy formation.
We use the halo occupation distribution (HOD) framework to characterise the predictions from two independent galaxy formation models for the galactic content of dark matter haloes and its evolution with redshift. Our galaxy samples correspond to a range of fixed number densities defined by stellar mass and span $0 \le z \le 3$. We find remarkable similarities between the model predictions. Differences arise at low galaxy number densities which are sensitive to the treatment of heating of the hot halo by active galactic nuclei. The evolution of the form of the HOD can be described in a relatively simple way, and we model each HOD parameter using its value at $z=0$ and an additional evolutionary parameter. In particular, we find that the ratio between the characteristic halo masses for hosting central and satellite galaxies can serve as a sensitive diagnostic for galaxy evolution models. Our results can be used to test and develop empirical studies of galaxy evolution and can facilitate the construction of mock galaxy catalogues for future surveys.
(abridged) We report rest-frame submillimeter H2O emission line observations of 11 HyLIRGs/ULIRGs at z~2-4 selected among the brightest lensed galaxies discovered in the Herschel-ATLAS. Using the IRAM NOEMA, we have detected 14 new H2O emission lines. The apparent luminosities of the H2O emission lines are $\mu L_{\rm{H_2O}} \sim 6-21 \times 10^8 L_\odot$, with velocity-integrated line fluxes ranging from 4-15 Jy km s$^{-1}$. We have also observed CO emission lines using EMIR on the IRAM 30m telescope in seven sources. The velocity widths for CO and H2O lines are found to be similar. With almost comparable integrated flux densities to those of the high-J CO line, H2O is found to be among the strongest molecular emitters in high-z Hy/ULIRGs. We also confirm our previously found correlation between luminosity of H2O ($L_{\rm{H_2O}}$) and infrared ($L_{\rm{IR}}$) that $L_{\rm{H_2O}} \sim L_{\rm{IR}}^{1.1-1.2}$, with our new detections. This correlation could be explained by a dominant role of far-infrared (FIR) pumping in the H2O excitation. Modelling reveals the FIR radiation fields have warm dust temperature $T_\rm{warm}$~45-75 K, H2O column density per unit velocity interval $N_{\rm{H_2O}}/\Delta V \gtrsim 0.3 \times 10^{15}$ cm$^{-2}$ km$^{-1}$ s and 100 $\mu$m continuum opacity $\tau_{100} > 1$ (optically thick), indicating that H2O is likely to trace highly obscured warm dense gas. However, further observations of $J\geq4$ H2O lines are needed to better constrain the continuum optical depth and other physical conditions of the molecular gas and dust. We have also detected H2O+ emission in three sources. A tight correlation between $L_{\rm{H_2O}}$ and $L_{\rm{H_2O^+}}$ has been found in galaxies from low to high redshift. The velocity-integrated flux density ratio between H2O+ and H2O suggests that cosmic rays generated by strong star formation are possibly driving the H2O+ formation.
We summarize the results obtained from our suite of chemical evolution models for spiral disks, computed for different total masses and star formation efficiencies. Once the gas, stars and star formation radial distributions are reproduced, we analyze the Oxygen abundances radial profiles for gas and stars, in addition to stellar averaged ages and global metallicity. We examine scenarios for the potential origin of the apparent flattening of abundance gradients in the outskirts of disk galaxies, in particular the role of molecular gas formation prescriptions.
We calculate the absorption efficiencies of composite silicate grains with inclusions of graphite and silicon carbide in the spectral range 5--25$\rm \mu m$. We study the variation in absorption profiles with volume fractions of inclusions. In particular we study the variation in the wavelength of peak absorption at 10 and 18$\rm \mu m$. We also study the variation of the absorption of porous silicate grains. We use the absorption efficiencies to calculate the infrared flux at various dust temperatures and compare with the observed infrared emission flux from the circumstellar dust around some M-Type \& AGB stars obtained from IRAS and a few stars from Spitzer satellite. We interpret the observed data in terms of the circumstellar dust grain sizes; shape; composition and dust temperature.
The CO+ reactive ion is thought to be a tracer of the boundary between a HII region and the hot molecular gas. In this study, we present the spatial distribution of the CO+ rotational emission toward the Mon R2 star-forming region. The CO+ emission presents a clumpy ring-like morphology, arising from a narrow dense layer around the HII region. We compare the CO+ distribution with other species present in photon-dominated regions (PDR), such as [CII] 158 mm, H2 S(3) rotational line at 9.3 mm, polycyclic aromatic hydrocarbons (PAHs) and HCO+. We find that the CO+ emission is spatially coincident with the PAHs and [CII] emission. This confirms that the CO+ emission arises from a narrow dense layer of the HI/H2 interface. We have determined the CO+ fractional abundance, relative to C+ toward three positions. The abundances range from 0.1 to 1.9x10^(-10) and are in good agreement with previous chemical model, which predicts that the production of CO+ in PDRs only occurs in dense regions with high UV fields. The CO+ linewidth is larger than those found in molecular gas tracers, and their central velocity are blue-shifted with respect to the molecular gas velocity. We interpret this as a hint that the CO+ is probing photo-evaporating clump surfaces.
There is strong observational evidence indicating a time lag of order of some 100 Myr between the onset of starburst and AGN activity in galaxies. Dynamical time lags have been invoked to explain this. We extend this approach by introducing a viscous time lag the gas additionally needs to flow through the AGN's accretion disc before it reaches the central black hole. Our calculations reproduce the observed time lags and are in accordance with the observed correlation between black hole mass and stellar velocity dispersion.
Ultra-compact dwarf galaxies (UCDs) share many properties with globular clusters (GCs) and are found in similar environments. A large sample of UCDs and GCs around NGC 1399, the central giant elliptical of the Fornax galaxy cluster, is used to infer their formation history and also that of NGC 1399. We assumed that all GCs and UCDs in our sample are star clusters (SCs) and used them as tracers of past star formation activities. After correcting our GC/UCD sample for mass loss, we interpreted their overall mass function to be a superposition of SC populations that formed coevally during different times. The SC masses of each population were distributed according to the embedded cluster mass function (ECMF), a pure power law with the slope $-\beta$ and a stellar upper mass limit, $M_{\mathrm{max}}$, which depended on the star formation rate (SFR). We decomposed the observed GC/UCD mass function into individual SC populations and converted $M_{\mathrm{max}}$ of each SC population to an SFR. The overall distribution of SFRs reveals how the GC/UCD sample formed. Considering the age of the GCs/UCDs and the present stellar mass of NGC 1399, we found that the formation of the GCs/UCDs can be well explained for $\beta<2.3$. This agrees very well with the observation in young SCs where $\beta\approx2.0$ is usually found. Even if taking into account that some of the most massive objects might not be genuine SCs and applying different corrections for the mass loss, the outcome is not influenced much. We found peak SFRs between approximately 300 and 3000 $M_{\odot}\mathrm{yr}^{-1}$, which matches the SFRs observed in massive high-$z$ sub-mm galaxies and an SFR estimate inferred from NGC 1399 based on "downsizing", i.e. more massive galaxies must have formed over shorter times. Our results suggest that NGC 1399 and its GC/UCD system formed in an early, short, and intense starburst.
Large-scale gaseous filaments with length up to the order of 100 pc are on the upper end of the filamentary hierarchy of the Galactic interstellar medium. Their association with respect to the Galactic structure and their role in Galactic star formation are of great interest from both observational and theoretical point of view. Previous "by-eye" searches, combined together, have started to uncover the Galactic distribution of large filaments, yet inherent bias and small sample size limit conclusive statistical results to be drawn. Here, we present (1) a new, automated method to identify large-scale velocity-coherent dense filaments, and (2) the first statistics and the Galactic distribution of these filaments. We use a customized minimum spanning tree algorithm to identify filaments by connecting voxels in the position-position-velocity space, using the Bolocam Galactic Plane Survey spectroscopic catalog. In the range of $7.^{\circ}5 \le l \le 194^{\circ}$, we have identified 54 large-scale filaments and derived mass ($\sim 10^3 - 10^5 \, M_\odot$), length (10-276 pc), linear mass density (54-8625 $M_\odot \, \rm{pc}^{-1}$), aspect ratio, linearity, velocity gradient, temperature, fragmentation, Galactic location and orientation angle. The filaments concentrate along major spiral arms. They are widely distributed across the Galactic disk, with 50% located within $\pm$20 pc from the Galactic mid-plane and 27% run in the center of spiral arms (aka "bones"). An order of 1% of the molecular ISM is confined in large filaments. Massive star formation is more favorable in large filaments compared to elsewhere. This is the first comprehensive catalog of large filaments useful for a quantitative comparison with spiral structures and numerical simulations.
We used the newly commissioned 50 cm Binocular Network (50BiN) telescope at Qinghai Station of Purple Mountain Observatory (Chinese Academy of Sciences) to observe the old open cluster NGC 188 in V and R as part of a search for variable objects. Our time-series data span a total of 36 days. Radial velocity and proper-motion selection resulted in a sample of 532 genuine cluster members. Isochrone fitting was applied to the cleaned cluster sequence, yielding a distance modulus of (m - M)0V = 11.35 \pm 0.10 mag and a total foreground reddening of E(V - R) = 0.062 \pm 0.002 mag. Light-curve solutions were obtained for eight W Ursae Majoris eclipsing-binary systems (W UMas) and their orbital parameters were estimated. Using the latter parameters, we estimate a distance to the W UMas which is independent of the host cluster's physical properties. Based on combined fits to six of the W UMas (EP Cep, EQ Cep, ES Cep, V369 Cep, and--for the first time--V370 Cep and V782 Cep), we obtain an average distance modulus of (m - M)0V = 11.31 \pm 0.08 mag, which is comparable with that resulting from our isochrone fits. These six W UMas exhibit an obvious period-luminosity relation. We derive more accurate physical parameters for the W UMa systems and discuss their initial masses and ages. The former show that these W UMa systems have likely undergone angular-momentum evolution within a convective envelope (W-type evolution). The ages of the W UMa systems agree well with the cluster's age.
We investigate anew the distribution of absolute carbon abundance, $A$(C) $= \log\,\epsilon $(C), for carbon-enhanced metal-poor (CEMP) stars in the halo of the Milky Way, based on high-resolution spectroscopic data for a total sample of 301 CEMP stars. The sample includes 147 CEMP-$s$ (and CEMP-r/s) stars, 124 CEMP-no stars, and 30 CEMP stars that are unclassified, based on the currently employed [Ba/Fe] criterion. We confirm previous claims that the distribution of $A$(C) for CEMP stars is (at least) bimodal, with newly determined peaks centered on $A$(C)$=7.93$ (the high-C region) and $A$(C)$ =6.28$ (the low-C region). A very high fraction of CEMP-$s$ (and CEMP-r/s) stars belong to the high-C region, while the great majority of CEMP-no stars reside in the low-C region. However, there exists complexity in the morphology of the $A$(C)-[Fe/H] space for the CEMP-no stars, a first indication that more than one class of first-generation stellar progenitors may be required to account for their observed abundances. The two groups of CEMP-no stars we identify exhibit clearly different locations in the $A$(Na)-$A$(C) and $A$(Mg)-$A$(C) spaces, also suggesting multiple progenitors. The clear distinction in $A$(C) between the CEMP-$s$ (and CEMP-$r/s$) stars and the CEMP-no stars appears to be $as\ successful$, and $likely\ more\ astrophysically\ fundamental$, for the separation of these sub-classes as the previously recommended criterion based on [Ba/Fe] (and [Ba/Eu]) abundance ratios. This result opens the window for its application to present and future large-scale low- and medium-resolution spectroscopic surveys.
Links to: arXiv, form interface, find, astro-ph, recent, 1607, contact, help (Access key information)