We present an investigation of the globular cluster (GC) systems of NGC 3607
and NGC 3608 as part of the ongoing SLUGGS survey. We use wide-field imaging
data from the Subaru telescope in the g, r, and i filters to analyse the radial
density, colour and azimuthal distributions of both GC systems. With the
complementary kinematic data obtained from the Keck II telescope, we measure
the radial velocities of a total of 81 GCs.
Our results show that the GC systems of NGC 3607 and NGC 3608 have a
detectable spatial extent of ~ 15, and 13 galaxy effective radii, respectively.
Both GC systems show a clear bimodal colour distribution. We detect a
significant radial colour gradient for the GC subpopulations in both galaxies.
NGC 3607 exhibits an overabundance of red GCs on the galaxy minor axis and NGC
3608 shows a misalignment in the GC subpopulation position angles with respect
to the galaxy stellar component.
With the aid of literature data, we discuss several relationships between the
properties of GC systems and their host galaxies. A one-to-one relation between
the ellipticities of red GCs and the galaxy stellar light emphasises the
evolutionary similarities between them. In our sample of four slowly rotating
galaxies with kinematically decoupled cores, we observe a higher ellipticity
for the blue GC subpopulation than their red counterparts. Also, we notice the
flattening of negative colour gradients for the blue GC subpopulations with
increasing galaxy stellar mass. Finally, we discuss the formation scenarios
associated with the blue GC subpopulation.
We measure galaxy sizes on a sample of $\sim1200$ galaxies with confirmed spectroscopic redshifts $2 \leq z_{spec} \leq 4.5$ in the VIMOS Ultra Deep Survey (VUDS), representative of star-forming galaxies with $i_\mathrm{AB} \leq 25$. We first derive galaxy sizes applying a classical parametric profile fitting method using GALFIT. We then measure the total pixel area covered by a galaxy above a given surface brightness threshold, which overcomes the difficulty of measuring sizes of galaxies with irregular shapes. We then compare the results obtained for the equivalent circularized radius enclosing 100\% of the measured galaxy light $r_T^{100}$ to those obtained with the effective radius $r_{e,\mathrm{circ}}$ measured with GALFIT. We find that the sizes of galaxies computed with our non-parametric approach span a large range but remain roughly constant on average with a median value $r_T^{100}\sim2.2$ kpc for galaxies with $2<z<4.5$. This is in stark contrast with the strong downward evolution of $r_e$ with increasing redshift, down to sizes of $<1$ kpc at $z\sim4.5$. We analyze the difference and find that parametric fitting of complex, asymmetric, multi-component galaxies is severely underestimating their sizes. By comparing $r_T^{100}$ with physical parameters obtained through SED fitting we find that the star-forming galaxies that are the largest at any redshift are, on average, more massive and more star-forming. We discover that galaxies present more concentrated light profiles as we move towards higher redshifts. We interpret these results as the signature of several, possibly different, evolutionary paths of galaxies in their early stages of assembly, including major and minor merging or star-formation in multiple bright regions. (abridged)
We aim at improving constraints on the epoch of galaxy formation by measuring the ages of 3597 galaxies with spectroscopic redshifts 2<z<6.5 in the VIMOS Ultra Deep Survey (VUDS). We derive ages and other physical parameters from the simultaneous fitting with the GOSSIP+ software of observed UV rest-frame spectra and photometric data from the u-band up to 4.5 microns using composite stellar population models. We conclude from extensive simulations that at z>2 the joint analysis of spectroscopy and photometry combined with restricted age possibilities when taking into account the age of the Universe substantially reduces systematic uncertainties and degeneracies in the age derivation. We find galaxy ages ranging from very young with a few tens of million years to substantially evolved with ages up to ~1.5-2 Gyr. The formation redshifts z_f derived from the measured ages indicate that galaxies may have started forming stars as early as z_f~15. We produce the formation redshift function (FzF), the number of galaxies per unit volume formed at a redshift z_f, and compare the FzF in increasing redshift bins finding a remarkably constant 'universal' FzF. The FzF is parametrized with (1+z)^\zeta, with \zeta~0.58+/-0.06, indicating a smooth 2 dex increase from z~15 to z~2. Remarkably this observed increase is of the same order as the observed rise in the star formation rate density (SFRD). The ratio of the SFRD with the FzF gives an average SFR per galaxy of ~7-17Msun/yr at z~4-6, in agreement with the measured SFR for galaxies at these redshifts. From the smooth rise in the FzF we infer that the period of galaxy formation extends from the highest possible redshifts that we can probe at z~15 down to redshifts z~2. This indicates that galaxy formation is a continuous process over cosmic time, with a higher number of galaxies forming at the peak in SFRD at z~2 than at earlier epochs. (Abridged)
This paper describes the first data release (DR1) of the VIMOS Ultra Deep Survey (VUDS). The DR1 includes all low-resolution spectroscopic data obtained in 276.9 arcmin2 of the CANDELS-COSMOS and CANDELS-ECFDS survey areas, including accurate spectroscopic redshifts z_spec and individual spectra obtained with VIMOS on the ESO-VLT. A total of 698 objects have a measured redshift, with 677 galaxies, two type-I AGN and a small number of 19 contaminating stars. The targets of the spectroscopic survey are selected primarily on the basis of their photometric redshifts to ensure a broad population coverage. About 500 galaxies have z_spec>2, 48 with z_spec>4, and the highest reliable redshifts reach beyond z_spec=6. This dataset approximately doubles the number of galaxies with spectroscopic redshifts at z>3 in these fields. We discuss the general properties of the sample in terms of the spectroscopic redshift distribution, the distribution of Lyman-alpha equivalent widths, and physical properties including stellar masses M_star and star formation rates (SFR) derived from spectral energy distribution fitting with the knowledge of z_spec. We highlight the properties of the most massive star-forming galaxies, noting the large range in spectral properties, with Lyman-alpha in emission or in absorption, and in imaging properties with compact, multi-component or pair morphologies. We present the catalogue database and data products. All data are publicly available and can be retrieved from a dedicated query-based database available at this http URL
We explore the impact of electron thermal conduction on the evolution of radiatively-cooled cold clouds embedded in flows of hot and fast material, as occur in outflowing galaxies. Performing a parameter study of three-dimensional adaptive mesh refinement hydrodynamical simulations, we show that electron thermal conduction causes cold clouds to evaporate, but it can also extend their lifetimes by compressing them into dense filaments. We distinguish between low column-density clouds, which are disrupted on very short times, and high-column density clouds with much-longer disruption times that are set by a balance between impinging thermal energy and evaporation. We provide fits to the cloud lifetimes and velocities that can be used in galaxy-scale simulations of outflows, in which the evolution of individual clouds cannot be modeled with the required resolution. Moreover, we show that the clouds are only accelerated to a small fraction of the ambient velocity because compression by evaporation causes the clouds to present a small cross-section to the ambient flow. This means that either magnetic fields must suppress thermal conduction, or that the cold clouds observed in galaxy outflows are not formed of cold material carried out from the galaxy.
We present the results of clustering analysis on the $z \sim 2$ star-forming galaxies. By combining our data with data from publicly available archives, we collect $g$-, $\zb / z$-, and $K$-band imaging data over 5.2 deg$^{2}$, which represents the largest area BzK/gzK survey. We apply colour corrections to translate our filter-set to those used in the original BzK selection for the gzK selection. Because of the wide survey area, we obtain a sample of 41,112 star-forming gzK galaxies at $z \sim 2$ (sgzKs) down to $\KAB < 23.0$, and determine high-quality two-point angular correlation functions (ACFs). Our ACFs show an apparent excess from power-law behaviour at small angular scale $(\theta \la 0.01^{\circ})$, which corresponds the virial radius of a dark halo at $z \sim 2$ with a mass of $\sim 10^{13} \Msun$. We find that the correlation lengths are consistent with the previous estimates over all magnitude range; however, our results are evaluated with a smaller margin of error than that in previous studies. The large amount of data enables us to determine ACFs differentially depending on the luminosity of the subset of the data. The mean halo mass of faint sgzKs $(22.0 < K \leq 23.0)$ was found to be $\Mh = (1.32^{+0.09}_{-0.12}) \times 10^{12} h^{-1} \Msun$, whereas bright sgzKs ($18.0 \leq K \leq 21.0)$ were found to reside in dark haloes with a mass of $\Mh = (3.26^{+1.23}_{-1.02}) \times 10^{13} h^{-1} \Msun$.
We present the results of an halo occupation distribution (HOD) analysis of star-forming galaxies at $z \sim 2$. We obtained high-quality angular correlation functions based on a large sgzK sample, which enabled us to carry out the HOD analysis. The mean halo mass and the HOD mass parameters are found to increase monotonically with increasing $K$-band magnitude, suggesting that more luminous galaxies reside in more massive dark haloes. The luminosity dependence of the HOD mass parameters was found to be the same as in the local Universe; however, the masses were larger than in the local Universe over all ranges of magnitude. This implies that galaxies at $z \sim 2$ tend to form in more massive dark haloes than in the local Universe, a process known as downsizing. By analysing the dark halo mass evolution using the extended Press--Schechter formalism and the number evolution of satellite galaxies in a dark halo, we find that faint Lyman break galaxies at $z \sim 4$ could evolve into the faintest sgzKs $(22.0 < K \leq 23.0)$ at $z \sim 2$ and into the Milky-Way-like galaxies or elliptical galaxies in the local Universe, whereas the most luminous sgzKs $(18.0 \leq K \leq 21.0)$ could evolve into the most massive systems in the local Universe. The stellar-to-halo mass ratio (SHMR) of the sgzKs was found to be consistent with the prediction of the model, except that the SHMR of the faintest sgzKs was smaller than the prediction at $z \sim 2$. This discrepancy may be attributed that our samples are confined to star-forming galaxies.
Using the exquisite depth of the Hubble Ultra Deep Field (HUDF12 programme) dataset, we explore the ongoing assembly of the outermost regions of the most massive galaxies (M_{stellar} > 5x10^{10} M_{Sun}) at z < 1. The outskirts of massive objects, particularly Early-Types Galaxies (ETGs), are expected to suffer a dramatic transformation across cosmic time due to continuous accretion of small galaxies. HUDF imaging allows us to study this process at intermediate redshifts in 6 massive galaxies, exploring the individual surface brightness profiles out to 25 effective radii. We find that 10-30% of the total stellar mass for the galaxies in our sample is contained within 10 < R < 50 kpc. These values are in close agreement with numerical simulations, and at least 2-3 times higher than those reported for late-type galaxies. The fraction of stellar mass stored in the outer envelopes/haloes of Massive Early-Type Galaxies increases with decreasing redshift, being 28.7% at < z > = 0.1, 22.6% at < z > = 0.65 and 3.5% at < z > = 2. The fraction of mass in diffuse features linked with ongoing minor merger events is > 1-3%, very similar to predictions based on observed close pair counts. Therefore, our results suggest that the size and mass growth of the most massive galaxies have been solely driven by minor and major merging from z = 1 to today.
We compare APOGEE radial velocities (RVs) of young stars in the Orion A cloud with CO line gas emission and find a correlation between the two at large-scales, in agreement with previous studies. However, at smaller scales we find evidence for the presence of substructure in the stellar velocity field. Using a Friends-of-Friends approach we identify 37 stellar groups with almost identical RVs. These groups are not randomly distributed but form elongated chains or strings of stars with five or more members with low velocity dispersion, across lengths of 1-1.5~pc. The similarity between the kinematic properties of the APOGEE strings and the internal velocity field of the chains of dense cores and fibers recently identified in the dense ISM is striking and suggests that for most of the Orion A cloud, young stars keep memory of the parental gas substructure where they originated.
The northern tile of the wide-area and shallow XMM-XXL X-ray survey field is
used to estimate the average dark matter halo mass of relatively luminous X-ray
selected AGN [$\rm
log\, L_X (\rm 2-10\,keV)= 43.6^{+0.4}_{-0.4}\,erg/s$] in the redshift
interval $z=0.5-1.2$. Spectroscopic follow-up observations of X-ray sources in
the XMM-XXL field by the Sloan telescope are combined with the VIPERS
spectroscopic galaxy survey to determine the cross-correlation signal between
X-ray selected AGN (total of 318) and galaxies (about 20,\,000). We model the
large scales (2-25\,Mpc) of the correlation function to infer a mean dark
matter halo mass of $\log M / (M_{\odot} \, h^{-1}) = 12.50 ^{+0.22} _{-0.30}$
for the X-ray selected AGN sample. This measurement is about 0.5\,dex lower
compared to estimates in the literature of the mean dark matter halo masses of
moderate luminosity X-ray AGN [$L_X (\rm 2-10\,keV)\approx 10^{42} -
10^{43}\,erg/s$] at similar redshifts. Our analysis also links the mean
clustering properties of moderate luminosity AGN with those of powerful
UV/optically selected QSOs, which are typically found in halos with masses few
times $10^{12}\,M_{\odot}$. There is therefore evidence for a negative
luminosity dependence of the AGN clustering. This is consistent with
suggestions that AGN have a broad dark matter halo mass distribution with a
high mass tail that becomes sub-dominant at high accretion luminosities. We
further show that our results are in qualitative agreement with semi-analytic
models of galaxy and AGN evolution, which attribute the wide range of dark
matter halo masses among the AGN population to different triggering mechanisms
and/or black hole fueling modes.
We expand on the dispersion analysis of polarimetry maps toward applications to interferometry data. We show how the filtering of low-spatial frequencies can be accounted for within the idealized Gaussian turbulence model, initially introduced for single-dish data analysis, to recover reliable estimates for correlation lengths of magnetized turbulence, as well as magnetic field strengths (plane-of-the-sky component) using the Davis-Chandrasekhar-Fermi method. We apply our updated technique to TADPOL/CARMA data obtained on W3(OH), W3 Main, and DR21(OH). For W3(OH) our analysis yields a turbulence correlation length $\delta\simeq19$ mpc, a ratio of turbulent-to-total magnetic energy $\left\langle B_{\mathrm{t}}^{2}\right\rangle /\left\langle B^{2}\right\rangle \simeq0.58$, and a magnetic field strength $B_{0}\sim1.1\:\mathrm{mG}$; for W3 Main $\delta\simeq22$ mpc, $\left\langle B_{\mathrm{t}}^{2}\right\rangle /\left\langle B^{2}\right\rangle \simeq0.74$, and $B_{0}\sim0.7\:\mathrm{mG}$; while for DR21(OH) $\delta\simeq12$ mpc, $\left\langle B_{\mathrm{t}}^{2}\right\rangle /\left\langle B^{2}\right\rangle \simeq0.70$, and $B_{0}\sim1.2\:\mathrm{mG}$.
We introduce the Solitary Local Dwarfs Survey (Solo), a wide field photometric study targeting every isolated dwarf galaxy within 3 Mpc of the Milky Way. Solo is based on (u)gi multi-band imaging from CFHT/MegaCam for northern targets, and Magellan/Megacam for southern targets. All galaxies fainter than Mv = -18 situated beyond the nominal virial radius of the Milky Way and M31 (>300 kpc) are included in this volume-limited sample, for a total of 42 targets. In addition to reviewing the survey goals and strategy, we present results for the Sagittarius Dwarf Irregular Galaxy (Sag DIG), one of the most isolated, low mass galaxies, located at the edge of the Local Group. We analyze its resolved stellar populations and their spatial distributions. We provide updated estimates of its central surface brightness and integrated luminosity, and trace its surface brightness profile to a level fainter than 30 mag./sq.arcsec. Sag DIG is well described by a highly elliptical (disk-like) system following a single component Sersic model. However, a low-level distortion is present at the outer edges of the galaxy that, were Sag DIG not so isolated, would likely be attributed to some kind of previous tidal interaction. Further, we find evidence of an extremely low level, extended distribution of stars beyond 5 arcmins (>1.5 kpc) that suggests Sag DIG may be embedded in a very low density stellar halo. We compare the stellar and HI structures of Sag DIG, and discuss results for this galaxy in relation to other isolated, dwarf irregular galaxies in the Local Group.
We present a detailed analysis of the chemistry and kinematics of red giants in the dwarf irregular galaxy NGC 6822. Spectroscopy at 8500 Angstroms was acquired for 72 red giant stars across two fields using FORS2 at the VLT. Line of sight extinction was individually estimated for each target star to accommodate the variable reddening across NGC 6822. The mean radial velocity was found to be v_helio = (52.8 +/- 2.2) km/s with dispersion rms = 24.1 km/s, in agreement with other studies. Ca II triplet equivalent widths were converted into [Fe/H] metallicities using a V magnitude proxy for surface gravity. The average metallicity was [Fe/H] = (-0.84 +/- 0.04) with dispersion rms = 0.31 dex and interquartile range 0.48. Our assignment of individual reddening values makes our analysis more sensitive to spatial variations in metallicity than previous studies. We divide our sample into metal-rich and metal-poor stars; the former are found to cluster towards small radii with the metal-poor stars more evenly distributed across the galaxy. The velocity dispersion of the metal-poor stars is higher than that of the metal-rich stars; combined with the age-metallicity relation this indicates that older populations have either been dynamically heated or were born in a less disclike distribution. The low ratio (v_rot/v_rms) suggests that within the inner 10', NGC 6822's stars are dynamically decoupled from the HI gas, possibly in a thick disc or spheroid.
Using a high-performance computing cluster to mosaic 4,787 pointings, we have imaged the 100 sq. deg. South Pole Telescope (SPT) deep-field at 2.1 GHz using the Australian Telescope Compact Array to an rms of 80 $\mu$Jy and a resolution of 8". Our goal is to generate an independent sample of radio-selected galaxy clusters to study how the radio properties compare with cluster properties at other wavelengths, over a wide range of redshifts in order to construct a timeline of their evolution out to $z \sim 1.3$. A preliminary analysis of the source catalogue suggests there is no spatial correlation between the clusters identified in the SPT-SZ catalogue and our wide-angle tail galaxies.
In the sixth of the series of papers reporting on a large reverberation mapping (RM) campaign of active galactic nuclei (AGNs) with high accretion rates, we present velocity-resolved time lags of H$\beta$ emission lines for nine objects observed in the campaign during 2012$-$2013. In order to correct the line-broadening caused by seeing and instruments before the analysis of velocity-resolved RM, we adopt Richardson-Lucy deconvolution to reconstruct their H$\beta$ profiles. The validity and effectiveness of the deconvolution are checked out by Monte Carlo simulation. Five among the nine objects show clear dependence of time delay on velocity. Mrk 335 and Mrk 486 show signatures of gas inflow whereas the clouds in the broad-line regions (BLRs) of Mrk 142 and MCG +06-26-012 tend to be radial outflowing. Mrk 1044 is consistent with the case of virialized motions. The lags of the rest four are not velocity-resolvable. The velocity-resolved RM of super-Eddington accreting massive black holes (SEAMBHs) shows that they have diversity of the kinematics in their BLRs. Comparing with the AGNs with sub-Eddington accretion rates, we do not find significant differences in the BLR kinematics of SEAMBHs.
We analyze the demographics of black holes (BHs) in the large-volume cosmological hydrodynamical simulation Horizon-AGN. This simulation statistically models how much gas is accreted onto BHs, traces the energy deposited into their environment and, consequently, the back-reaction of the ambient medium on BH growth. The synthetic BHs reproduce a variety of observational constraints such as the redshift evolution of the BH mass density and the mass function. Yet there seem to be too many BHs with mass~ 1e7 Msun at high redshift, and too few BHs with similar mass at z=0 in intermediate-mass galaxies. Strong self-regulation via AGN feedback, weak supernova feedback, and unresolved internal process are likely to be responsible for this, and for a tight BH-galaxy mass correlation. Starting at z~2, tidal stripping creates a small population of BHs over-massive with respect to the halo. The fraction of galaxies hosting a central BH or an AGN increases with stellar mass. The AGN fraction agrees better with multi-wavelength studies, than single-wavelength ones, unless obscuration is taken into account. The most massive halos present BH multiplicity, with additional BHs gained by ongoing or past mergers. In some cases, both a central and an off-center AGN shine concurrently, producing a dual AGN. This dual AGN population dwindles with decreasing redshift, as found in observations. Specific accretion rate and Eddington ratio distributions are in good agreement with observational estimates. The BH population is dominated in turn by fast, slow, and very slow accretors, with transitions occurring at z=3 and z=2 respectively.
The MgII emission line is visible in the optical band for intermediate redshift quasars (0.4<z<1.6) and it is thus an extremely important tool to measure the black hole mass and to understand the structure of the Broad Line Region. We aim to determine the substructure and the variability of the MgII line with the aim to identify which part of the line comes from a medium in Keplerian motion. Using SALT telescope we performed ten spectroscopic observation of a quasar HE 0435-4312 (z=1.2231) over the period of 3 years (Dec 23/24, 2012 to Dec 7/8, 2015). We find that the line is well modeled by two Lorentzian components, and the relative strength of these components vary with time. The line maximum is shifted in a time-dependent way from the position of the Fe II pseudo-continuum, although the effect is not very strong, and the line asymmetry varies in time. We also note very different local conditions in the formation region of Mg II and FeII. The timescale for the line shape variability is of the order of the light travel time to the emitting region, therefore the changes are most likely due to the varying irradiation patterns and the presence of the two components does not imply two distinct emission regions.
In this paper, we test the discrepancy between the stellar mass density and instantaneous star formation rate in redshift range $0<z<8$ using a large observational data sample. We first compile the measurements of the stellar mass densities up to $z\sim 8$. Comparing the observed stellar mass densities with the time-integral of instantaneous star formation history, we find that the observed stellar mass densities are lower than that implied from star formation history at $z<4$. We also use Markov chain monte carlo method to derive the best-fitting star formation history from the observed stellar mass density data. At $0.5<z<6$, the observed star formation rate densities are larger than the best-fitting one, especially at $z\sim2$ where by a factor of about two. However, at lower ($z<0.5$) and higher redshifts ($z>6$), the derived star formation history is consistent with the observations. This is the first time to test the discrepancy between the observed stellar mass density and instantaneous star formation rate up to very high redshift $z\approx8$ using the Markov chain monte carlo method and a varying recycling factor. Several possible reasons for this discrepancy are discussed, such as underestimation of stellar mass density, initial mass function and cosmic metallicity evolution.
Observations of Galactic H I gas for seven intermediate Galactic latitude fields are presented at 1' angular resolution using data from the DRAO Synthesis Telescope (ST) and the Green Bank Telescope (GBT). The DHIGLS data are the most extensive arcminute resolution measurements of the diffuse atomic interstellar medium beyond those in the Galactic plane. The acquisition, reduction, calibration, and mosaicking of the DRAO ST data and the cross calibration and incorporation of the short-spacing information from the GBT are described. The high quality of the DHIGLS data enables a variety of new studies in directions of low Galactic column density. We find evidence for dramatic changes in the structures in channel maps over even small changes in velocity. This narrow line emission has counterparts in absorption spectra against bright background radio sources, quantifying that the gas is cold and dense and can be identified as the cold neutral medium phase. We analyze the angular power spectra of maps of the integrated H I emission (column density) from the mosaics for several distinct velocity ranges. Fitting power spectrum models based on a power law, but including the effects of the synthesized beam and noise at high spatial frequencies, we find exponents ranging from -2.5 to -3.0. Power spectra of maps of the centroid velocity for these components give similar results. These exponents are interpreted as being representative of the 3D density and velocity fields of the atomic gas, respectively. Fully reduced DHIGLS H I data cubes and other data products are available at www.cita.utoronto.ca/DHIGLS.
[Abrigded] We have started a new deep, multi-imaging survey of the Fornax cluster, dubbed Fornax Deep Survey (FDS), at the VLT Survey Telescope. In this paper we present the deep photometry inside two square degrees around the bright galaxy NGC1399 in the core of the cluster. We found a very extended and diffuse envelope surrounding the luminous galaxy NGC1399: we map the surface brightness out to 33 arcmin (~ 192 kpc) from the galaxy center and down to about 31 mag/arcsec^2 in the g band. The deep photometry allows us to detect a faint stellar bridge in the intracluster region between NGC1399 and NGC1387. By analyzing the integrated colors of this feature, we argue that it could be due to the ongoing interaction between the two galaxies, where the outer envelope of NGC1387 on its east side is stripped away. By fitting the light profile, we found that it exists a physical break radius in the total light distribution at R=10 arcmin (~58 kpc) that sets the transition region between the bright central galaxy and the outer exponential stellar halo. We discuss the main implications of this work on the build-up of the stellar halo at the center of the Fornax cluster. By comparing with the numerical simulations of the stellar halo formation for the most massive BCGs, we find that the observed stellar halo mass fraction is consistent with a halo formed through the multiple accretion of progenitors with a stellar mass in the range 10^8 - 10^11 M_sun. This might suggest that the halo of NGC1399 has also gone through a major merging event. The absence of a significant number of luminous stellar streams and tidal tails out to 192 kpc suggests that the epoch of this strong interaction goes back to an early formation epoch. Therefore, differently from the Virgo cluster, the extended stellar halo around NGC1399 is characterised by a more diffuse and well-mixed component, including the ICL.
We present Pipe3D, an analysis pipeline based on the FIT3D fitting tool, devel- oped to explore the properties of the stellar populations and ionized gas of Integral Field Spectroscopy data. Pipe3D was created to provide with coherent, simple to distribute, and comparable dataproducts, independently of the origin of the data, focused on the data of the most recent IFU surveys (e.g., CALIFA, MaNGA, and SAMI), and the last generation IFS instruments (e.g., MUSE). Along this article we describe the different steps involved in the analysis of the data, illustrating them by showing the dataproducts derived for NGC 2916, observed by CALIFA and P-MaNGA. As a practical use of the pipeline we present the complete set of dataproducts derived for the 200 datacubes that comprises the V500 setup of the CALIFA Data Release 2 (DR2), making them freely available through the network (this ftp URL). Finally, we explore the hypothesis that the properties of the stellar populations and ionized gas of galaxies at the effective radius are representative of the overall average ones, finding that this is indeed the case.
The formation of massive black hole binaries (MBHBs) is an unavoidable outcome of galaxy evolution via successive mergers. However, the mechanism that drives their orbital evolution from parsec separations down to the gravitational wave (GW) dominated regime is poorly understood and their final fate is still unclear. If such binaries are embedded in gas-rich and turbulent environments, as observed in remnants of galaxy mergers, the interaction with gas clumps (such as molecular clouds) may efficiently drive their orbital evolution. Using numerical simulations, we test this hypothesis by studying the dynamical evolution of an equal-mass, circular MBHB accreting infalling molecular clouds. We investigate different orbital configurations, modelling a total of 13 systems to explore different possible pericentre distances and relative inclinations of the cloud-binary encounter. We show that the evolution of the binary orbit is dominated by the exchange of angular momentum through gas accretion during the first stages of the interaction for all orbital configurations. Building on these results, we construct a simple model for evolving a MBHB interacting with a sequence of clouds, which are randomly drawn from reasonable populations with different levels of anisotropy in their angular momenta distributions. We show that the binary efficiently evolves down to the GW emission regime within a few hundred million years, overcoming the 'final parsec' problem regardless of the stellar distribution.
A model, which is widely used for inertial rang statistics of supersonic turbulence in the context of molecular clouds and star formation, expresses (measurable) relative scaling exponents Z_p of two-point velocity statistics as a function of two parameters, beta and Delta. The model relates them to the dimension D of the most dissipative structures, D=3-Delta/(1-beta). While this description has proved most successful for incompressible turbulence (beta=Delta=2/3, and D=1), its applicability in the highly compressible regime remains debated. For this regime, theoretical arguments suggest D=2 and Delta=2/3, or Delta=1. Best estimates based on 3D periodic box simulations of supersonic isothermal turbulence yield Delta=0.71 and D=1.9, with uncertainty ranges of Delta in [0.67, 0.78] and D in [2.04,1.60]. With these 5-10\% uncertainty ranges just marginally including the theoretical values of Delta=2/3 and D=2, doubts remain whether the model indeed applies and, if it applies, for what values of beta and Delta. We use a Monte Carlo approach to mimic actual simulation data and examine what factors are most relevant for the fit quality. We estimate that 0.1% (0.05%) accurate Z_p, with p=1...5, should allow for 2% (1%) accurate estimates of beta and Delta in the highly compressible regime, but not in the mildly compressible regime. We argue that simulation-based Z_p with such accuracy are within reach of today's computer resources. If this kind of data does not allow for the expected high quality fit of beta and Delta, then this may indicate the inapplicability of the model for the simulation data. In fact, other models than the one we examine here have been suggested.
Links to: arXiv, form interface, find, astro-ph, recent, 1602, contact, help (Access key information)
The scaling of disk galaxy rotation velocity with baryonic mass (the "Baryonic Tully-Fisher" relation, BTF) has long confounded galaxy formation models. It is steeper than the M ~ V^3 scaling relating halo virial masses and circular velocities and its zero point implies that galaxies comprise a very small fraction of available baryons. Such low galaxy formation efficiencies may in principle be explained by winds driven by evolving stars, but the tightness of the BTF relation argues against the substantial scatter expected from such vigorous feedback mechanism. We use the APOSTLE/EAGLE simulations to show that the BTF relation is well reproduced in LCDM simulations that match the size and number of galaxies as a function of stellar mass. In such models, galaxy rotation velocities are proportional to halo virial velocity and the steep velocity-mass dependence results from the decline in galaxy formation efficiency with decreasing halo mass needed to reconcile the CDM halo mass function with the galaxy luminosity function. Despite the strong feedback, the scatter in the simulated BTF is smaller than observed, even when considering all simulated galaxies and not just rotationally-supported ones. The simulations predict that the BTF should become increasingly steep at the faint end, although the velocity scatter at fixed mass should remain small. Observed galaxies with rotation speeds below ~40 km/s seem to deviate from this prediction. We discuss observational biases and modeling uncertainties that may help to explain this disagreement in the context of LCDM models of dwarf galaxy formation.
We report the results from some of the deepest Keck/MOSFIRE data yet obtained for candidate $z \gtrsim 7$ galaxies. Our data show one significant line detection with 6.5$\sigma$ significance in our combined 10 hours of integration which is independently detected on more than one night, ruling out the possibility that the detection is spurious. The asymmetric line profile and non-detection in the optical bands strongly imply that the detected line is Ly$\alpha$ emission from a galaxy at $z$(Ly$\alpha)=7.6637 \pm 0.0011$, making it the fourth spectroscopically confirmed galaxy at $z>7.5$. This galaxy is bright in the rest-frame ultraviolet (UV; $M_{\rm UV} \sim -21.2$) with a moderately blue UV slope ($\beta=-2.2^{+0.3}_{-0.2}$), and exhibits a rest-frame Ly$\alpha$ equivalent width of EW(Ly$\alpha$) $\sim 15.6^{+5.6}_{-3.6}$ \AA. The non-detection of the 11 other $z \sim$ 7--8 galaxies in our long 10 hr integration, reaching a median 5$\sigma$ sensitivity of 28 \AA\ in the rest-frame EW(Ly$\alpha$), implies a 1.3$\sigma$ deviation from the null hypothesis of a non-evolving distribution in the rest-frame EW(Ly$\alpha$) between $3<z<6$ and $z=$ 7--8. Our results are consistent with previous studies finding a decline in Ly$\alpha$ emission at $z>6.5$, which may signal the evolving neutral fraction in the intergalactic medium at the end of the reionization epoch, although our weak evidence suggests the need for a larger statistical sample to allow for a more robust conclusion.
The Star formation rate (SFR) is a crucial parameter to investigate galaxy evolution. At low redshift the cosmic SFR density declines smoothly, and massive active galaxies become passive, reducing their star formation activity. This implies that the bulk of the SFR density at low redshift is mainly driven by low mass objects. We investigate the properties of a sample of low luminosity Far-Infrared (FIR) sources selected at 250 microns from Pappalardo et al. (2015). We have collected data from Ultraviolet to FIR to perform a multi-wavelengths analysis. The main goal is to investigate the correlation between SFR, stellar mass, and dust mass for a galaxy population with a wide range in dust content and stellar mass, including the low mass regime that most probably dominates the SFR density at low z. We define a main sample of ~800 sources with full Spectral Energy Distribution (SED) coverage between 0.15 < lambda < 500 microns and an extended sample with ~5000 sources in which we remove the constraints on the Ultraviolet and Near-Infrared bands. We analyze both samples with two different SED fitting methods: MAGPHYS and CIGALE. In the SFR versus stellar mass plane our samples occupy a region included between local spirals and higher redshift star forming galaxies. The galaxies subsample with the higher masses (M* > 3e10 Msol) does not lie on the main sequence, but shows a small offset, as a consequence of the decreased star formation. Low mass galaxies (M* < 1e10 Msol) settle in the main sequence with SFR and stellar mass consistent with local spirals. Deep Herschel data allow the identification of a mixed galaxy population, with galaxies still in an assembly phase, or galaxies at the beginning of their passive evolution. We find that the dust luminosity is the parameter that discriminates these two galaxy populations.
Here we fitted a 4-arm spiral structure to the more accurate data on global
arm pitch angle and arm longitude tangents, to get the start of each spiral arm
near the Galactic nucleus. We find that the tangent to the 'start of the
Sagittarius' spiral arm (arm middle) is at l= -17 degrees +/- 0.5 degree, while
the tangent to the 'start of the Norma' spiral arm (arm middle) is at l= +20
degrees +/- 0.5 degree.
Earlier, we published a compilation of observations and analysis of the
tangent to each spiral arm tracer, from longitudes +23 degrees to +340 degrees;
here we cover the arm tracers in the remaining longitudes +340 degrees (=- 20
degrees) to +23 degrees. Our model arm tangents are confirmed through the
recent observed masers data (at the arm's inner edge). Observed arm tracers in
the inner Galaxy show an offset from the mid-arm; this was also found elsewhere
in the Milky Way disk (Vallee 2014c).
In addition, we collated the observed tangents to the so-called '3-kpc-arm'
features; here they are found statistically to be near l= -18 degrees +/- 2
degrees and near l= +21 degrees +/- 2 degrees, after excluding misidentified
spiral arms. We find that the model-computed arm tangents in the inner Galaxcy
are spatially coincident with the mean longitude of the observed tangents to
the '3-kpc-arm' features (same galactic longitudes, within the errors). These
spatial similarities may be suggestive of a contiguous space.
Multiple color selection techniques have been successful in identifying quasars from wide-field broad-band imaging survey data. Among the quasars that have been discovered so far, however, there is a redshift gap at $5 \lesssim {\rm z} \lesssim 5.7$ due to the limitations of filter sets in previous studies. In this work, we present a new selection technique of high redshift quasars using a sequence of medium-band filters: nine filters with central wavelengths from 625 to 1025 nm and bandwidths of 50 nm. Photometry with these medium-bands traces the spectral energy distribution (SED) of a source, similar to spectroscopy with resolution R $\sim$ 15. By conducting medium-band observations of high redshift quasars at 4.7 $\leq$ z $\leq$ 6.0 and brown dwarfs (the main contaminants in high redshift quasar selection) using the SED camera for QUasars in EArly uNiverse (SQUEAN) on the 2.1-m telescope at the McDonald Observatory, we show that these medium-band filters are superior to multi-color broad-band color section in separating high redshift quasars from brown dwarfs. In addition, we show that redshifts of high redshift quasars can be determined to an accuracy of $\Delta{\rm z}/(1+{\rm z}) = 0.002$ -- $0.026$. The selection technique can be extended to z $\sim$ 7, suggesting that the medium-band observation can be powerful in identifying quasars even at the re-ionization epoch.
We model the dust and free-free continuum emission in the high-mass
star-forming region Sagittarius B2 in order to reconstruct the
three-dimensional density and dust temperature distribution, as a crucial input
to follow-up studies of the gas velocity field and molecular abundances. We
employ the three-dimensional radiative transfer program RADMC-3D to calculate
the dust temperature self-consistently, provided a given initial density
distribution. This density distribution of the entire cloud complex is then
recursively reconstructed based on available continuum maps, including both
single-dish and high-resolution interferometric maps covering a wide frequency
range (40 GHz - 4 THz). The model covers spatial scales from 45 pc down to 100
au, i.e. a spatial dynamic range of 10^5.
We find that the density distribution of Sagittarius B2 can be reasonably
well fitted by applying a superposition of spherical cores with Plummer-like
density profiles. In order to reproduce the spectral energy distribution, we
position Sgr B2(N) along the line of sight behind the plane containing Sgr
B2(M). We find that the entire cloud complex comprises a total gas mass of 8.0
x 10^6 Msun within a diameter of 45 pc, corresponding to an averaged gas
density of 170 Msun/pc^3. We estimate stellar masses of 2400 Msun and 20700
Msun and luminosities of 1.8 x 10^6 Lsun and 1.2 x 10^7 Lsun for Sgr B2(N) and
Sgr B2(M), respectively. We report H_2 column densities of 2.9 x 10^24 cm^-2
for Sgr B2(N) and 2.5 x 10^24 cm^-2 for Sgr B2(M) in a 40" beam. For Sgr B2(S),
we derive a stellar mass of 1100 Msun, a luminosity of 6.6 x 10^5 Lsun and a
H_2 column density of 2.2 x 10^24 cm^-2 in a 40" beam. We calculate a star
formation efficiency of 5% for Sgr B2(N) and 50% for Sgr B2(M), indicating that
most of the gas content in Sgr B2(M) has already been converted to stars or
dispersed.
We analyze an optical data cube of the nuclear region of NGC 3621, taken with the integral field unit of the Gemini Multi-object Spectrograph. We found that the previously detected central line emission in this galaxy actually comes from a blob, located at a projected distance of 2.14" +/- 0.08" (70.1 +/- 2.6 pc) from the stellar nucleus. Only diffuse emission was detected in the rest of the field of view, with a deficit of emission at the position of the stellar nucleus. Diagnostic diagram analysis reveals that the off-centered emitting blob has a Seyfert 2 spectrum. We propose that the line-emitting blob may be a "fossil" emission-line region or a light "echo" from an active galactic nucleus (AGN), which was significantly brighter in the past. Our estimates indicate that the bolometric luminosity of the AGN must have decreased by a factor of ~13 - 500 during the last ~230 years. A second scenario to explain the morphology of the line-emitting areas in the nuclear region of NGC 3621 involves no decrease of the AGN bolometric luminosity and establishes that the AGN is highly obscured toward the observer but not toward the line-emitting blob. The third scenario proposed here assumes that the off-centered line-emitting blob is a recoiling supermassive black hole, after the coalescence of two black holes. Finally, an additional hypothesis is that the central X-ray source is not an AGN, but an X-ray binary. This idea is consistent with all the scenarios we proposed.
We present observations of the 1.3 mm continuum emission toward hub-N and hub-S of the infrared dark cloud G14.225-0.506 carried out with the Submillimeter Array, together with observations of the dust emission at 870 and 350 microns obtained with APEX and CSO telescopes. The large scale dust emission of both hubs consists of a single peaked clump elongated in the direction of the associated filament. At small scales, the SMA images reveal that both hubs fragment into several dust condensations. The fragmentation level was assessed under the same conditions and we found that hub-N presents 4 fragments while hub-S is more fragmented, with 13 fragments identified. We studied the density structure by means of a simultaneous fit of the radial intensity profile at 870 and 350 microns and the spectral energy distribution adopting a Plummer-like function to describe the density structure. The parameters inferred from the model are remarkably similar in both hubs, suggesting that density structure could not be responsible in determining the fragmentation level. We estimated several physical parameters such as the level of turbulence and the magnetic field strength, and we found no significant differences between these hubs. The Jeans analysis indicates that the observed fragmentation is more consistent with thermal Jeans fragmentation compared with a scenario that turbulent support is included. The lower fragmentation level observed in hub-N could be explained in terms of stronger UV radiation effects from a nearby HII region, evolutionary effects, and/or stronger magnetic fields at small scales, a scenario that should be further investigated.
The Dragonfish Nebula has been recently claimed to be powered by a superluminous but elusive OB association. Instead, systematic searches in near-infrared photometric surveys have found many other cluster candidates on this sky region. Among these, the first confirmed young massive cluster was Mercer 30, where Wolf-Rayet stars were found. We perform a new characterization of Mercer 30 with unprecedented accuracy, combining NICMOS/HST and VVV photometric data with multi-epoch ISAAC/VLT H- and K-band spectra. Stellar parameters for most of spectroscopically observed cluster members are found through precise non-LTE atmosphere modeling with the CMFGEN code. Our spectrophotometric study for this cluster yields a new, revised distance of d = (12.4 +- 1.7) kpc and a total of Q = 6.70 x 10^50 Lyman ionizing photons. A cluster age of (4.0 +- 0.8) Myr is found through isochrone fitting, and a total mass of (1.6 +- 0.6) x 10^4 Msol is estimated thanks to our extensive knowledge of the post-main-sequence population. As a consequence, membership of Mercer 30 to the Dragonfish star-forming complex is confirmed, allowing us to use this cluster as a probe for the whole complex, which turns out to be extremely large (400 pc across) and located at the outer edge of the Sagittarius-Carina spiral arm (11 kpc from the Galactic Center). The Dragonfish complex hosts 19 young clusters or cluster candidates (including Mercer 30 and a new candidate presented in this work) and an estimated minimum of 9 field Wolf-Rayet stars. The sum of all these contributions accounts for, at least, 73% of the Dragonfish Nebula ionization and leaves little or no room for the alleged superluminous OB association; alternative explanations are discussed.
Our knowledge of the shape of radial metallicity gradients in disc galaxies has recently improved. Conversely, the understanding of their time evolution is more complex, since it requires analysis of stellar populations with different ages, or systematic studies of galaxies at different redshifts. In the Local Universe, Hii regions and planetary nebulae (PNe) are important tools to investigate it. We present an in-depth study of all nearby spiral galaxies (M33, M31, NGC300, and M81) with direct-method nebular abundances of both populations. For the first time, we also evaluate the radial migration of PN populations. We analyse Hii region and PN properties to: determine whether oxygen in PNe is a reliable tracer for past interstellar medium (ISM) composition; homogenise the published datasets; estimate the migration of the oldest stellar populations; determine the overall chemical enrichment and slope evolution. We confirm that oxygen in PNe is a reliable tracer for the past ISM metallicity. We find that PN gradients are flatter than or equal to those of Hii regions. When radial motions are negligible, this result provides a direct measurement of the time evolution of the gradient. For galaxies with dominant radial motions, we provide upper limits on the gradient evolution. Finally, the total metal content increases with time in all target galaxies, with early morphological type having a larger increment Delta(O/H) than late-type galaxies. Our findings provide important constraints to discriminate among different galactic evolutionary scenarios, favouring cosmological models with enhanced feedback from supernovae. The advent of extremely large telescopes will allow us to include galaxies in a wider range of morphologies and environments, thus putting firmer constraints to galaxy formation and evolution scenarios.
We present the results of observations of 22GHz H2O maser emission in NGC 6240 and M51 made with the Karl G. Jansky Very Large Array. Two major H2O maser features and several minor features are detected toward the southern nucleus of NGC 6240. These features are redshifted by about 300 km/s from the galaxy's systemic velocity and remain unresolved at the synthesized beam size. A combination of our two-epoch observations and published data reveals an apparent correlation between the strength of the maser and the 22GHz radio continuum emission, implying that the maser excitation relates to the activity of an active galactic nucleus in the southern nucleus rather than star-forming activity. The star-forming galaxy M51 hosts H2O maser emission in the center of the galaxy; however, the origin of the maser has been an open question. We report the first detection of 22GHz nuclear radio continuum emission in M51. The continuum emission is co-located with the maser position, which indicates that the maser arises from active galactic nucleus-activity and not from star-forming activity in the galaxy.
The Circum-Galactic Medium (CGM) can be probed through the analysis of absorbing systems in the line-of-sight to bright background quasars. We present measurements of the metallicity of a new sample of 15 sub-damped Lyman-$\alpha$ absorbers (sub-DLAs, defined as absorbers with 19.0 < log N(H I) < 20.3) with redshift 0.584 < $\rm z_{abs}$ < 3.104 from the ESO Ultra-Violet Echelle Spectrograph (UVES) Advanced Data Products Quasar Sample (EUADP). We combine these results with other measurements from the literature to produce a compilation of metallicity measurements for 92 sub-DLAs as well as a sample of 362 DLAs. We apply a multi-element analysis to quantify the amount of dust in these two classes of systems. We find that either the element depletion patterns in these systems differ from the Galactic depletion patterns or they have a different nucleosynthetic history than our own Galaxy. We propose a new method to derive the velocity width of absorption profiles, using the modeled Voigt profile features. The correlation between the velocity width delta_V90 of the absorption profile and the metallicity is found to be tighter for DLAs than for sub-DLAs. We report hints of a bimodal distribution in the [Fe/H] metallicity of low redshift (z < 1.25) sub-DLAs, which is unseen at higher redshifts. This feature can be interpreted as a signature from the metal-poor, accreting gas and the metal-rich, outflowing gas, both being traced by sub-DLAs at low redshifts.
We present the results of the most complete ever scan of the parameter space for cosmic ray (CR) injection and propagation. We perform a Bayesian search of the main GALPROP parameters, using the MultiNest nested sampling algorithm, augmented by the BAMBI neural network machine learning package. This is the first such study to separate out low-mass isotopes ($p$, $\bar p$ and He) from the usual light elements (Be, B, C, N, O). We find that the propagation parameters that best fit $p$, $\bar p$, He data are significantly different from those that fit light elements, including the B/C and $^{10}$Be/$^9$Be secondary-to-primary ratios normally used to calibrate propagation parameters. This suggests each set of species is probing a very different interstellar medium, and that the standard approach of calibrating propagation parameters using B/C can lead to incorrect results. We present posterior distributions and best fit parameters for propagation of both sets of nuclei, as well as for the injection abundances of elements from H to Si. The input GALDEF files with these new parameters will be included in an upcoming public GALPROP update.
We present ALMA Cycle 1 observations of the HH 46/47 molecular outflow using combined 12m array and 7m array observations. We use 13CO and C18O emission to correct for the 12CO optical depth, to accurately estimate the outflow mass, momentum and kinetic energy. Applying the optical depth correction increases the mass estimate by a factor of 14, the momentum by a factor of 6, and the kinetic energy by a factor of about 2. The new 13CO(1-0) and C18O(1-0) data also allow us to trace denser and slower outflow material than that traced by 12CO. These species are only detected within about 1~2 km/s from the cloud velocity. The cavity wall of the red lobe appears at very low velocities (~0.2 km/s). Combing the material traced only by 13CO and C18O, the measured total mass of the CO outflow is 1.4 Msun, the total momentum is 1.7 Msun km/s and the total energy is 4.7e43 erg, assuming Tex=15 K. The improved angular resolution and sensitivity in 12CO reveal more details of the outflow structure. Specifically, we find that the outflow cavity wall is composed of multiple shells entrained in a series of jet bow-shock events. The outflow kinetic energy distribution shows that even though the red lobe is mainly entrained by jet bow-shocks, more outflow energy is being deposited into the cloud at the base of the outflow cavity rather than around the heads of the bow shocks. The estimated outflow mass, momentum, and energy indicate that the outflow is capable to disperse the parent core within the typical lifetime of the embedded phase of a low-mass protostar, and regulating a core-to-star efficiency of 1/4~1/3. The 13CO and C18O emission also trace a circumstellar envelope with rotation and infall motions. In CS, we found possible evidence for a slowly-moving rotating outflow, which we believe is entrained not only poloidally but also toroidally by a wind launched from relatively large radii on the disk.
We discuss the science drivers for ALMA Band 2 which spans the frequency range from 67 to 90 GHz. The key science in this frequency range are the study of the deuterated molecules in cold, dense, quiescent gas and the study of redshifted emission from galaxies in CO and other species. However, Band 2 has a range of other applications which are also presented. The science enabled by a single receiver system which would combine ALMA Bands 2 and 3 covering the frequency range 67 to 116 GHz, as well as the possible doubling of the IF bandwidth of ALMA to 16 GHz, are also considered.
Expanding upon our previous work (Rodriguez et al., 2015), we study merging binary black holes formed in globular clusters using our Monte Carlo approach to stellar dynamics. We have created a new set of 52 cluster models with different masses, metallicities, and radii to fully characterize the binary black hole merger rate. These models include all the relevant dynamical processes (such as two-body relaxation, strong encounters, and three-body binary formation) and agree well with detailed direct N-body simulations. In addition, we have enhanced our stellar evolution algorithms with updated metallicity-dependent stellar wind and supernova prescriptions, allowing us to compare our results directly to the most recent population synthesis predictions for merger rates from isolated binary evolution. We explore the relationship between a cluster's global properties and the population of binary black holes that it produces. In particular, we derive a numerically calibrated relationship between the merger times of ejected black hole binaries and a cluster's mass and radius. We explore the masses and mass ratios of these binaries as a function of redshift, and find a merger rate of ~5 Gpc$^{-3}$ yr$^{-1}$ in the local universe, with 90% of sources having total masses from $32M_{\odot}$ to $64M_{\odot}$. Under standard assumptions, approximately 1 out of every 7 binary black hole mergers in the local universe will have originated in a globular cluster, but we also explore the sensitivity of this result to different assumptions for binary stellar evolution. If black holes were born with significant natal kicks, comparable to those of neutron stars, then the merger rate of binary black holes from globular clusters would be comparable to that from the field, with approximately 1/2 of mergers originating in clusters [Abridged].
The white dwarf luminosity function is an important tool to understand the properties of the Solar neighborhood, like its star formation history, and its age. Here we present a population synthesis study of the white dwarf population within 40~pc from the Sun, and compare the results of this study with the properties of the observed sample. We use a state-of-the-art population synthesis code based on Monte Carlo techniques, that incorporates the most recent and reliable white dwarf cooling sequences, an accurate description of the Galactic neighborhood, and a realistic treatment of all the known observational biases and selection procedures. We find a good agreement between our theoretical models and the observed data. In particular, our simulations reproduce a previously unexplained feature of the bright branch of the white dwarf luminosity function, which we argue is due to a recent episode of star formation. We also derive the age of the Solar neighborhood employing the position of the observed cut-off of the white dwarf luminosity function, obtaining ~8.9+-0.2 Gyr. We conclude that a detailed description of the ensemble properties of the population of white dwarfs within 40pc of the Sun allows us to obtain interesting constraints on the history of the Solar neighborhood.
We interpret the stellar population of $\omega$ Centauri by means of a population synthesis analysis, following the most recent observational guidelines for input metallicities, helium and [(C+N+O)/Fe] contents. We deal at the same time with the main sequences, sub-giant and horizontal branch data. The reproduction of the observed colour magnitude features is very satisfying and bears interesting hints concerning the evolutionary history of this peculiar stellar ensemble. Our main results are: 1) no significant spread in age is required to fit the colour-magnitude diagram. Indeed we can use coeval isochrones for the synthetic populations, and we estimate that the ages fall within a $\sim 0.5$ Gyr time interval; in particular the most metal rich population can be coeval (in the above meaning) with the others, if its stars are very helium--rich (Y$\sim$0.37) and with the observed CNO enhancement ([(C+N+O)/Fe] = + 0.7); 2) a satisfactory fit of the whole HB is obtained, consistent with the choice of the populations providing a good reproduction of the main sequence and sub giant data. 3) the split in magnitude observed in the red HB is well reproduced assuming the presence of two stellar populations in the two different sequences observed: a metal poor population made of stars evolving from the blue side (luminous branch) and a metal richer one whose stars are in a stage closer to the zero age HB (dimmer branch). This modelization also fits satisfactorily the period and the [Fe/H] distribution of the RR Lyrae stars.
Using the largest single-survey sample of Type Ia supernovae (SNe Ia) to date, we study the relationship between properties of SNe Ia and those of their host galaxies, focusing primarily on correlations with Hubble residuals (HR). Our sample consists of 345 photometrically-classified or spectroscopically-confirmed SNeIa discovered as part of the SDSS-II Supernova Survey (SDSS-SNS). This analysis utilizes host-galaxy spectroscopy obtained during the SDSS-I/II spectroscopic survey and from an ancillary program on the SDSS-III Baryon Oscillation Spectroscopic Survey (BOSS) that obtained spectra for nearly all host galaxies of SDSS-II SN candidates. In addition, we use photometric host-galaxy properties from the SDSS-SNS data release (Sako et al. 2014) such as host stellar mass and star-formation rate. We confirm the well-known relation between HR and host-galaxy mass and find a 3.6{\sigma} significance of a non-zero linear slope. We also recover correlations between HR and host-galaxy gas-phase metallicity and specific star-formation rate as they are reported in the literature. With our large dataset, we examine correlations between HR and multiple host-galaxy properties simultaneously and find no evidence of a significant correlation. We also independently analyze our spectroscopically-confirmed and photometrically-classified SNe Ia and comment on the significance of similar combined datasets for future surveys.
We have obtained K-band Multi-Object Spectrograph (KMOS) near-IR spectroscopy
for 14 red supergiant stars (RSGs) in the young massive star cluster NGC 2100
in the Large Magellanic Cloud (LMC). Stellar parameters including metallicity
are estimated using the J-band analysis technique, which has been rigorously
tested in the Local Universe. We find an average metallicity for NGC 2100 of
[Z]=$-$0.38$\pm$0.20 dex, in good agreement with estimates from the literature
for the LMC. Comparing our results in NGC 2100 with those for a Galactic
cluster (at Solar-like metallicity) with a similar mass and age we find no
significant difference in the location of RSGs in the Hertzsprung--Russell
diagram. We combine the observed KMOS spectra to form a simulated
integrated-light cluster spectrum and show that, by analysing this spectrum as
a single RSG, the results are consistent with the average properties of the
cluster.
Radial velocities are estimated for the targets and the dynamical properties
are estimated for the first time within this cluster. The data are consistent
with a flat velocity dispersion profile, and with an upper limit of 3.9 \kms,
at the 95\% confidence level, for the velocity dispersion of the cluster.
However, the intrinsic velocity dispersion is unresolved and could, therefore,
be significantly smaller than the upper limit reported here. An upper limit on
the dynamical mass of the cluster is derived as $M_{dyn}$ $\le$
$15.2\times10^{4}M_{\odot}$ assuming virial equilibrium.
Links to: arXiv, form interface, find, astro-ph, recent, 1602, contact, help (Access key information)
The Time-Domain Spectroscopic Survey (TDSS) is an SDSS-IV eBOSS subproject primarily aimed at obtaining identification spectra of ~220,000 optically-variable objects systematically selected from SDSS/Pan-STARRS1 multi-epoch imaging. We present a preview of the science enabled by TDSS, based on TDSS spectra taken over ~320 deg^2 of sky as part of the SEQUELS survey in SDSS-III, which is in part a pilot survey for eBOSS in SDSS-IV. Using the 15,746 TDSS-selected single-epoch spectra of photometrically variable objects in SEQUELS, we determine the demographics of our variability-selected sample, and investigate the unique spectral characteristics inherent in samples selected by variability. We show that variability-based selection of quasars complements color-based selection by selecting additional redder quasars, and mitigates redshift biases to produce a smooth quasar redshift distribution over a wide range of redshifts. The resulting quasar sample contains systematically higher fractions of blazars and broad absorption line quasars than from color-selected samples. Similarly, we show that M-dwarfs in the TDSS-selected stellar sample have systematically higher chromospheric active fractions than the underlying M-dwarf population, based on their H-alpha emission. TDSS also contains a large number of RR Lyrae and eclipsing binary stars with main-sequence colors, including a few composite-spectrum binaries. Finally, our visual inspection of TDSS spectra uncovers a significant number of peculiar spectra, and we highlight a few cases of these interesting objects. With a factor of ~15 more spectra, the main TDSS survey in SDSS-IV will leverage the lessons learned from these early results for a variety of time-domain science applications.
We present a strong lensing system, composed of 4 multiple images of a source at z = 2.387, created by two lens galaxies, G1 and G2, belonging to the galaxy cluster MACS J1115.9+0129 at z = 0.353. We use observations taken as part of the Cluster Lensing and Supernova survey with Hubble, CLASH, and its spectroscopic follow-up programme at the Very Large Telescope, CLASH-VLT, to estimate the total mass distributions of the two galaxies and the cluster through strong gravitational lensing models. We find that the total projected mass values within the half-light radii, R_{e}, of the two lens galaxies are M_{T,G1}(< R_{e,G1}) = (3.6 +/- 0.4) x 10^{10}M_{Sun} and M_{T,G2}(< R_{e,G2}) = (4.2 +/- 1.6) x 10^{10}M_{Sun}. The effective velocity dispersion values of G1 and G2 are (122 +/- 7) km/s and (137 +/- 27) km/s, respectively. We remark that these values are relatively low when compared to those of ~200-300 km/s, typical of lens galaxies found in the field by previous surveys. By fitting the spectral energy distributions of G1 and G2, we measure projected luminous over total mass fractions within R_{e} of 0.11 +/- 0.03, for G1, and 0.73 +/- 0.32, for G2. The fact that the less massive galaxy, G1, is dark-matter dominated in its inner regions raises the question of whether the dark matter fraction in the core of early-type galaxies depends on their mass. Further investigating strong lensing systems will help us understand the influence that dark matter has on the structure and evolution of the inner regions of galaxies.
Core collapse supernovae (CCSNe) are important sources of interstellar dust, potentially capable of producing one solar mass of dust in their explosively expelled ejecta. However, unlike other dust sources, the dust has to survive the passage of the reverse shock, generated by the interaction of the supernova blast wave with its surrounding medium. Knowledge of the net amount of dust produced by CCSNe is crucial for understanding the origin and evolution of dust in the local and high-redshift universe. Our aim is to identify the dust destruction mechanisms in the ejecta, and derive the net amount of dust that survives the passage of the reverse shock. We use analytical models for the evolution of a supernova blast wave and of the reverse shock, with special application to the clumpy ejecta of the remnant of Cassiopeia A. We assume that the dust resides in cool oxygen-rich clumps that are uniformly distributed within the remnant and surrounded by a hot X-ray emitting plasma, and that the dust consists of silicates (MgSiO3) and amorphous carbon grains. The passage of the reverse shock through the clumps gives rise to a relative gas-grain motion and also destroys the clumps. Inside the ejecta clouds, dust is processed via kinetic sputtering, which is terminated either when the grains escape the clumps, or when the clumps are destroyed by the reverse shock. In either case, grain destruction proceeds thereafter by thermal sputtering in the hot ambient gas. We find that 11.8% and 15.9% of, respectively, the silicate and carbon dust survives the passage of the reverse shock by the time the shock has reached the center of the remnant. These fractions depend on the morphology of the ejecta and the medium into which the remnant is expanding, as well as the composition and size distribution of the grains that formed in the ejecta. Results will therefore differ for different types of supernovae.
We investigate the relation between star formation rates ($\dot{M}_{s}$) and AGN properties in optically selected type 1 quasars at $2<z<3$ using data from Herschel and the SDSS. We find that $\dot{\rm{M}}_s$ remains approximately constant with redshift, at $300\pm100~\rm{M}_{\odot}$yr$^{-1}$. Conversely, $\dot{\rm{M}}_s$ increases with AGN luminosity, up to a maximum of $\sim600~\rm{M}_{\odot}$yr$^{-1}$, and with CIV FWHM. In context with previous results, this is consistent with a relation between $\dot{\rm{M}}_s$ and black hole accretion rate ($\dot{\rm{M}}_{bh}$) existing in only parts of the $z-\dot{\rm{M}}_{s}-\dot{\rm{M}}_{bh}$ plane, dependent on the free gas fraction, the trigger for activity, and the processes that may quench star formation. The relations between $\dot{\rm{M}}_s$ and both AGN luminosity and CIV FWHM are consistent with star formation rates in quasars scaling with black hole mass, though we cannot rule out a separate relation with black hole accretion rate. Star formation rates are observed to decline with increasing CIV equivalent width. This decline can be partially explained via the Baldwin effect, but may have an additional contribution from one or more of three factors; $M_i$ is not a linear tracer of L$_{2500}$, the Baldwin effect changes form at high AGN luminosities, and high CIV EW values signpost a change in the relation between $\dot{\rm{M}}_s$ and $\dot{\rm{M}}_{bh}$. Finally, there is no strong relation between $\dot{\rm{M}}_s$ and Eddington ratio, or the asymmetry of the CIV line. The former suggests that star formation rates do not scale with how efficiently the black hole is accreting, while the latter is consistent with CIV asymmetries arising from orientation effects.
We present the first results from our CAlibrating LYMan-$\alpha$ with H$\alpha$ (CALYMHA) pilot survey at the Isaac Newton Telescope. We measure Ly$\alpha$ emission for 488 H$\alpha$ selected galaxies at $z=2.23$ from HiZELS in the COSMOS and UDS fields with a specially designed narrow-band filter ($\lambda_c$ = 3918 {\AA}, $\Delta\lambda$= 52 {\AA}). We find 17 dual H$\alpha$-Ly$\alpha$ emitters ($f_{\rm Ly\alpha} >5\times10^{-17}$ erg s$^{-1}$ cm$^{-2}$, of which 5 are X-ray AGN). For star-forming galaxies, we find a range of Ly$\alpha$ escape fractions (f$_{\rm esc}$, measured with 3$"$ apertures) from $2$\%$-30$\%. These galaxies have masses from $3\times10^8$ M$_{\odot}$ to 10$^{11}$ M$_{\odot}$ and dust attenuations E$(B-V)=0-0.5$. Using stacking, we measure a median escape fraction of $1.6\pm0.5$\% ($4.0\pm1.0$\% without correcting H$\alpha$ for dust), but show that this depends on galaxy properties. The stacked f$_{\rm esc}$ tends to decrease with increasing SFR and dust attenuation. However, at the highest masses and dust attenuations, we detect individual galaxies with f$_{\rm esc}$ much higher than the typical values from stacking, indicating significant scatter in the values of f$_{\rm esc}$. Relations between f$_{\rm esc}$ and UV slope are bimodal, with high f$_{\rm esc}$ for either the bluest or reddest galaxies. We speculate that this bimodality and large scatter in the values of f$_{\rm esc}$ is due to additional physical mechanisms such as outflows facilitating f$_{\rm esc}$ for dusty/massive systems. Ly$\alpha$ is significantly more extended than H$\alpha$ and the UV. f$_{\rm esc}$ continues to increase up to at least 20 kpc (3$\sigma$, 40 kpc [2$\sigma$]) for typical SFGs and thus the aperture is the most important predictor of f$_{\rm esc}$.
The baryonic Tully-Fisher relation (BTFR) is both a valuable observational tool and a critical test of galaxy formation theory. We explore the systematic uncertainty in the slope and the scatter of the observed BTFR utilizing a homogeneously measured dataset of 930 isolated galaxies. We measure a fiducial relation of log_10 M_baryon = 3.24 log_10 V_rot + 3.21 with a scatter of 0.25 dex over the baryonic mass range of 10^7.4 to 10^11.3 M_sun. We then conservatively vary the definitions of M_baryon and V_rot, the sample definition and the linear fitting algorithm used to fit the BTFR. We obtain slopes ranging from 2.64 to 3.46 and scatter measurements ranging from 0.16 to 0.41 dex. We next compare our fiducial slope to literature measurements, where reported slopes range from 3.0 to 4.3 and scatter is either unmeasured, unmeasurable or as large as 0.4 dex. Measurements derived from unresolved HI line-widths tend to produce slopes of 3.2, while measurements derived strictly from resolved asymptotic rotation velocities produce slopes of 4.0. The largest factor affecting the BTFR slope is the definition of rotation velocity. Sample definition, mass range and linear fitting algorithm also significantly affect the measured BTFR. Galaxies with V_rot < 100 km/s are consistent with the BTFR of more massive galaxies, but these galaxies drive most of the scatter in the BTFR. This is most likely due to the diversity in rotation curve shapes of low-mass galaxies and underestimated systematic uncertainties. It is critical when comparing predictions to an observed BTFR that the rotation velocity definition, the sample selection and the fitting algorithm are similarly defined. Fitting a power-law model to the BTFR is an oversimplification and we recommend direct statistical comparisons between datasets with commensurable properties.
Now that ALMA is reaching its full capabilities, observations of sub-mm emission line deep fields become feasible. Deep fields are ideal to study the luminosity function of sub-mm emission lines, ultimately tracing the atomic and molecular gas properties of galaxies. We couple a semi-analytic model of galaxy formation with a radiative transfer code to make predictions for the luminosity function of CO J=1-0 up to CO J=6-5 and [CII] at redshifts z=0-6. We find that: 1) our model correctly reproduces the CO and [CII] emission of low- and high-redshift galaxies and reproduces the available constraints on the CO luminosity function at z<2.75; 2) we find that the CO and [CII] luminosity functions of galaxies increase from z = 6 to z = 4, remain relatively constant till z = 1 and rapidly decrease towards z = 0. The galaxies that are brightest in CO and [CII] are found at z~2; 3) the CO J=3-2 emission line is most favourable to study the CO luminosity and global H2 mass content of galaxies, because of its brightness and observability with currently available sub-mm and radio instruments; 4) the luminosity functions of high-J CO lines show stronger evolution than the luminosity functions of low-J CO lines; 5) our model barely reproduces the available constraints on the CO and [CII] luminosity function of galaxies at z>1.5 and the CO luminosity of individual galaxies at intermediate redshifts. We argue that this is driven by a lack of cold gas in galaxies at intermediate redshifts as predicted by cosmological simulations of galaxy formation. This may lay at the root of other problems theoretical models face at the same redshifts.
Young massive clusters (YMCs) have central stellar mass surface densities exceeding $10^{4} M_{\odot} pc^{-2}$. It is currently unknown whether the stars formed at such high (proto)stellar densities. We compile a sample of gas clouds in the Galaxy which have sufficient gas mass within a radius of a few parsecs to form a YMC, and compare their radial gas mass distributions to the stellar mass distribution of Galactic YMCs. We find that the gas in the progenitor clouds is distributed differently than the stars in YMCs. The mass surface density profiles of the gas clouds are generally shallower than the stellar mass surface density profiles of the YMCs, which are characterised by prominent dense core regions with radii ~ 0.1 pc, followed by a power-law tail. On the scale of YMC core radii, we find that there are no known clouds with significantly more mass in their central regions when compared to Galactic YMCs. Additionally, we find that models in which stars form from very dense initial conditions require surface densities that are generally higher than those seen in the known candidate YMC progenitor clouds. Our results show that the quiescent, less evolved clouds contain less mass in their central regions than in the highly star-forming clouds. This suggests an evolutionary trend in which clouds continue to accumulate mass towards their centres after the onset of star formation. We conclude that a conveyor-belt scenario for YMC formation is consistent with the current sample of Galactic YMCs and their progenitor clouds.
Various unification schemes interpret the complex phenomenology of quasars and luminous active galactic nuclei (AGN) in terms of a simple picture involving a central black hole, an accretion disc and an associated outflow. Here, we continue our tests of this paradigm by comparing quasar spectra to synthetic spectra of biconical disc wind models, produced with our state-of-the-art Monte Carlo radiative transfer code. Previously, we have shown that we could produce synthetic spectra resembling those of observed broad absorption line (BAL) quasars, but only if the X-ray luminosity was limited to $10^{43}$ erg s$^{-1}$. Here, we introduce a simple treatment of clumping, and find that a filling factor of $\sim0.01$ moderates the ionization state sufficiently for BAL features to form in the rest-frame UV at more realistic X-ray luminosities. Our fiducial model shows good agreement with AGN X-ray properties and the wind produces strong line emission in, e.g., Ly \alpha\ and CIV 1550\AA\ at low inclinations. At high inclinations, the spectra possess prominent LoBAL features. Despite these successes, we cannot reproduce all emission lines seen in quasar spectra with the correct equivalent-width ratios, and we find an angular dependence of emission-line equivalent width despite the similarities in the observed emission line properties of BAL and non-BAL quasars. Overall, our work suggests that biconical winds can reproduce much of the qualitative behaviour expected from a unified model, but we cannot yet provide quantitative matches with quasar properties at all viewing angles. Whether disc winds can successfully unify quasars is therefore still an open question.
The "main sequence of galaxies" $-$ defined in terms of the total star formation rate $\psi$ vs. the total stellar mass $M_*$ $-$ is a well-studied tight relation that has been observed at several wavelengths and at different redshifts. All earlier studies have derived this relation from integrated properties of galaxies. We recover the same relation from an analysis of spatially-resolved properties, with integral field spectroscopic (IFS) observations of 306 galaxies from the CALIFA survey. We consider the SFR surface density in units of log(M$_{\odot}$ yr$^{-1}$ Kpc$^{-2}$) and the stellar mass surface density in units of log(M$_{\odot}$ Kpc$^{-2}$) in individual spaxels which probe spatial scales of 0.5-1.5 Kpc. This local relation exhibits a high degree of correlation with small scatter ($\sigma = 0.23$ dex), irrespective of the dominant ionisation source of the host galaxy or its integrated stellar mass. We highlight: $(i)$ the integrated star formation main sequence formed by galaxies whose dominant ionisation process is related to star formation, for which we find a slope of 0.81 $\pm 0.02$; (ii) the spatially-resolved relation obtained with the spaxel analysis, we find a slope of 0.72 $\pm 0.04$; (iii) for the integrated main sequence we identified also a sequence formed by galaxies that are dominated by an old stellar population, which we have called the retired galaxies sequence.
Recent studies found higher galaxy metallicities in richer environments. It is not yet clear, however, whether metallicity-environment dependencies are merely an indirect consequence of environmentally-dependent formation histories, or of environment-related processes affecting metallicity directly. Here we present a first detailed study of metallicity-environment correlations in a cosmological hydrodynamical simulation, in particular the Illustris simulation. Illustris galaxies display similar relations to those observed. Utilizing knowledge of simulated formation histories, and leveraging the large simulation volume, we construct galaxy samples of satellites and centrals that are matched in formation histories. This allows us to find that ~1/3 of the metallicity-environment correlation is due to different formation histories in different environments. This is a combined effect of satellites (in particular in denser environments) having average lower z=0 star-formation rates (SFRs), and of their older stellar ages even at a given z=0 SFR. Most of the difference, ~2/3, is however caused by the higher concentration of the star-forming disks of satellite galaxies, since it biases their SFR-weighted metallicities towards their inner, more metal-rich, parts. With a newly defined quantity, the 'radially-averaged' metallicity, which captures the metallicity profile but is independent of the SFR profile, metallicities of satellites and centrals become environmentally independent, once matched in formation history. We find that circumgalactic metallicity (defined as rapidly inflowing gas around the virial radius), while sensitive to environment, has no measurable effect on the metallicity of the star-forming gas inside the galaxies.
In this paper we present the results of our search for and study of $z \gtrsim 6$ galaxy candidates behind the third Frontier Fields (FF) cluster, MACSJ0717.5+3745, and its parallel field, combining data from Hubble and Spitzer. We select 39 candidates using the Lyman Break technique, for which the clear non-detection in optical make the extreme mid-$z$ interlopers hypothesis unlikely. We also take benefit from $z \gtrsim 6$ samples selected using previous Frontier Fields datasets of Abell 2744 and MACS0416 to improve the constraints on the properties of very high-redshift objects. We compute the redshift and the physical properties, such emission lines properties, star formation rate, reddening, and stellar mass for all Frontier Fields objects from their spectral energy distribution using templates including nebular emission lines. We study the relationship between several physical properties and confirm the trend already observed in previous surveys for evolution of star formation rate with galaxy mass, and between the size and the UV luminosity of our candidates. The analysis of the evolution of the UV Luminosity Function with redshift seems more compatible with an evolution of density. Moreover, no robust $z\ge$8.5 object is selected behind the cluster field, and few $z$$\sim$9 candidates have been selected in the two previous datasets from this legacy survey, suggesting a strong evolution in the number density of galaxies between $z$$\sim$8 and 9. Thanks to the use of the lensing cluster, we study the evolution of the star formation rate density produced by galaxies with L$>$0.03L$^{\star}$, and confirm the strong decrease observed between $z$$\sim$8 and 9.
We study the relationship between stellar mass, star formation rate (SFR),ionization state, and gas-phase metallicity for a sample of 41 normal star-forming galaxies at $3 \lesssim z \lesssim 3.7$. The gas-phase oxygen abundance, ionization parameter, and electron density of ionized gas are derived from rest-frame optical strong emission lines measured on near-infrared spectra obtained with Keck/MOSFIRE. We remove the effect of these strong emission lines in the broad-band fluxes to compute stellar masses via spectral energy distribution fitting, while the SFR is derived from the dust-corrected ultraviolet luminosity. The ionization parameter is weakly correlated with the specific SFR, but otherwise the ionization parameter and electron density do not correlate with other global galaxy properties such as stellar mass, SFR, and metallicity. The mass-metallicity relation (MZR) at $z\simeq3.3$ shows lower metallicity by $\simeq 0.7$ dex than that at $z=0$ at the same stellar mass. Our sample shows an offset by $\simeq 0.3$ dex from the locally defined mass-metallicity-SFR relation, indicating that simply extrapolating such relation to higher redshift may predict an incorrect evolution of MZR. Furthermore, within the uncertainties we find no SFR-metallicity correlation, suggesting a less important role of SFR in controlling the metallicity at high redshift. We finally investigate the redshift evolution of the MZR by using the model by Lilly et al. (2013), finding that the observed evolution from $z=0$ to $z\simeq3.3$ can be accounted for by the model assuming a weak redshift evolution of the star formation efficiency.
Several studies indicate that radio-loud (RL) Active Galactic Nuclei (AGN) are produced only by the most massive black holes (BH), $M_{\rm BH} \sim 10^8$-$10^{10} M_\odot$. This idea has been challenged by the discovery of RL Narrow Line Seyfert 1 (RL NLSy1), having estimated masses of $M_{\rm BH}$$\sim$$10^6$-$10^7$ M$_\odot$. However, these low $M_{\rm BH}$ estimates might be due to projection effects. Spectropolarimetry allows us to test this possibility by looking at RL NLSy1s under a different perspective, i.e., from the viewing angle of the scattering material. We here report the results of a pilot study of VLT spectropolarimetric observations of the RL NLSy1 PKS 2004-447. Its polarization properties are remarkably well reproduced by models in which the scattering occurs in an equatorial structure surrounding its broad line region, seen close to face-on. In particular, we detect a polarized H$\alpha$ line with a width of $\sim$ 9,000 km s$^{-1}$, $\sim 6$ times broader than the width seen in direct light. This corresponds to a revised estimate of $M_{\rm BH}$$\sim$$6\times10^8$ M$_\odot$, well within the typical range of RL AGN. The double-peaked polarized broad H$\alpha$ profile of the target suggests that the rare combination of the orientation effects and a broad line region dominated by the rotation might account for this class of objects, casting doubts on the virial estimates of BH mass for type-I AGN.
We present a detailed study of the shape of the innermost part of the rotation curves of a sample of 29 nearby spiral galaxies, based on high angular and spectral resolution kinematic Halpha Fabry-Perot observations. In particular, we quantify the steepness of the rotation curve by measuring its slope dRvc(0). We explore the relationship between the inner slope and several galaxy parameters, such as stellar mass, maximum rotational velocity, central surface brightness ({\mu}0), bar strength and bulge-to-total ratio. Even with our limited dynamical range, we find a trend for low-mass galaxies to exhibit shallower rotation curve inner slopes than high-mass galaxies, whereas steep inner slopes are found exclusively in high-mass galaxies. This trend may arise from the relationship between the total stellar mass and the mass of the bulge, which are correlated among them. We find a correlation between the inner slope of the rotation curve and the morphological T-type, complementary to the scaling relation between dRvc(0) and {\mu}0 previously reported in the literature. Although we find that the inner slope increases with the Fourier amplitude A2 and decreases with the bar torque Qb, this may arise from the presence of the bulge implicit in both A2 and Qb. As previously noted in the literature, the more compact the mass in the central parts of a galaxy (more concretely, the presence of a bulge), the steeper the inner slopes. We conclude that the baryonic matter dominates the dynamics in the central parts of our sample galaxies.
The all-Galaxy CO survey of Dame, Hartmann, & Thaddeus (2001) is by far the most uniform, large-scale Galactic CO survey. Using a dendrogram-based decomposition of this survey, we present a catalog of 1064 massive molecular clouds throughout the Galactic plane. This catalog contains $2.5 \times 10^8$ solar masses, or $25^{+10.7}_{-5.8} \%$ of the Milky Way's estimated H$_2$ mass. We track clouds in some spiral arms through multiple quadrants. The power index of Larson's first law, the size-linewidth relation, is consistent with 0.5 in all regions - possibly due to an observational bias - but clouds in the inner Galaxy systematically have significantly (~ 30%) higher linewidths at a given size, indicating that their linewidths are set in part by Galactic environment. The mass functions of clouds in the inner Galaxy versus the outer Galaxy are both qualitatively and quantitatively distinct. The inner Galaxy mass spectrum is best described by a truncated power-law with a power index of $\gamma=-1.6\pm0.1$ and an upper truncation mass $M_0 = (1.0 \pm 0.2) \times 10^7 M_\odot$, while the outer Galaxy mass spectrum is better described by a non-truncating power law with $\gamma=-2.2\pm0.1$ and an upper mass $M_0 = (1.5 \pm 0.5) \times 10^6 M_\odot$, indicating that the inner Galaxy is able to form and host substantially more massive GMCs than the outer Galaxy. Additionally, we have simulated how the Milky Way would appear in CO from extragalactic perspectives, for comparison with CO maps of other galaxies.
Recent studies have found that radio-AGN selected by radio-loudness show little difference in terms of their host galaxy properties when compared to non-AGN galaxies of similar stellar mass and redshift. Using new 1.4~GHz VLBI observations of the COSMOS field we find that approximately 49$\pm8$\% of high-mass (M $>$ 10$^{10.5}$ M$_{\odot}$), high luminosity (L$_{1.4}$ $>$ 10$^{24}$ W~Hz$^{-1}$) radio-AGN possess a VLBI detected counterpart. These objects show no discernible bias towards specific stellar masses, redshifts or host properties other than what is shown by the radio-AGN population in general. Radio-AGN that are detected in VLBI observations are not special, but form a representative sample of the radio-loud AGN population.
A blind HI survey of the extragalactic sky behind the southern Milky Way has been conducted with the multibeam receiver on the 64-m Parkes radio telescope. The survey covers the Galactic longitude range 212 < l < 36 and Galactic latitudes |b| < 5, and yields 883 galaxies to a recessional velocity of 12,000 km/s. The survey covers the sky within the HIPASS area to greater sensitivity, finding lower HI-mass galaxies at all distances, and probing more completely the large-scale structures at and beyond the distance of the Great Attractor. Fifty-one percent of the HI detections have an optical/NIR counterpart in the literature. A further 27% have new counterparts found in existing, or newly obtained, optical/NIR images. The counterpart rate drops in regions of high foreground stellar crowding and extinction, and for low-HI mass objects. Only 8% of all counterparts have a previous optical redshift measurement. A notable new galaxy is HIZOA J1353-58, a possible companion to the Circinus galaxy. Merging this catalog with the similarly-conducted northern extension (Donley et al. 2005), large-scale structures are delineated, including those within the Puppis and Great Attractor regions, and the Local Void. Several newly-identified structures are revealed here for the first time. Three new galaxy concentrations (NW1, NW2 and NW3) are key in confirming the diagonal crossing of the Great Attractor Wall between the Norma cluster and the CIZA J1324.7-5736 cluster. Further contributors to the general mass overdensity in that area are two new clusters (CW1 and CW2) in the nearer Centaurus Wall, one of which forms part of the striking 180 deg (100/h Mpc) long filament that dominates the southern sky at velocities of ~3000 km/s, and the suggestion of a further Wall at the Great Attractor distance at slightly higher longitudes.
We report the discovery of significant mass/light offsets in the strong gravitational lensing system SDSS\,J1011$+$0143. We use the high-resolution \textsl{Hubble Space Telescope} (\textsl{HST}) F555W- and F814W-band imaging and SDSS spectroscopy of this system, which consists of a close galaxy pair with a projected separation of $\approx 4.2$ kpc at $z_{\rm lens} \sim 0.331$ lensing a Lyman-alpha emitter (LAE) at $z_{\rm source} = 2.701$. Comparisons between the mass peaks inferred from lens models and light peaks from \textsl{HST} imaging data reveal significant spatial mass/light offsets as large as $(1.72 \pm 0.24 \pm 0.34)$ kpc in both filter bands. Such large mass/light offsets, not seen in isolated field lens galaxies and relaxed galaxy groups, may be related to the interactions between the two lens galaxies. The detected mass/light offsets can potentially serve as an important test for the self-interacting dark matter model. However, other mechanisms such as dynamical friction on spatially differently distributed dark matter and stars could produce similar offsets. Detailed hydrodynamical simulations of galaxy-galaxy interactions with self-interacting dark matter could accurately quantify the effects of different mechanisms. The background LAE is found to contain three distinct star-forming knots with characteristic sizes from 116 pc to 438 pc. It highlights the power of strong gravitational lensing in probing the otherwise too faint and unresolved structures of distance objects below sub-kiloparsec or even 100-parsec scales through its magnification effect.
A precise measure of the mid-infrared interstellar extinction law is crucial to the investigation of the properties of interstellar dust, especially of the grains in the large size end. Based on the stellar parameters derived from the SDSS-III/APOGEE spectroscopic survey, we select a large sample of G- and K-type giants as the tracers of the Galactic mid-infrared extinction. We calculate the intrinsic stellar color excesses from the stellar effective temperatures and use them to determine the mid-infrared extinction for a given line of sight. For the entire sky of the Milky Way surveyed by APOGEE, we derive the extinction (relative to the K$_{\rm S}$ band at wavelength $\lambda=2.16\mu$m) for the four \emph{WISE} bands at 3.4, 4.6, 12 and 22$\mu$m, the four \emph{Spitzer}/IRAC bands at 3.6, 4.5, 5.8 and 8$\mu$m, the \emph{Spitzer}/MIPS24 band at 23.7$\mu$m and for the first time, the \emph{AKARI}/S9W band at 8.23$\mu$m. Our results agree with previous works in that the extinction curve is flat in the ~3--8$\mu$m wavelength range and is generally consistent with the $R_V=5.5$ model curve except our determination exceeds the model prediction in the \emph{WISE}/W4 band. Although some previous works found that the mid-IR extinction law appears to vary with the extinction depth $A_{\rm{K_S}}$, no noticeable variation has been found in this work. The uncertainties are analyzed in terms of the bootstrap resampling method and Monte-Carlo simulation and are found to be rather small.
We present three dimensional relativistic hydrodynamical simulations of a precessing jet interacting with the intracluster medium and compare the simulated jet structure with the observed structure of the Hydra A northern jet. For the simulations, we use jet parameters obtained in the parameter space study of the first paper in this series and probe different values for the precession period and precession angle. We find that for a precession period P = 1 Myr and a precession angle = 20 degree the model reproduces i) the curvature of the jet, ii) the correct number of bright knots within 20 kpc at approximately correct locations, and iii) the turbulent transition of the jet to a plume. The Mach number of the advancing bow shock = 1.85 is indicative of gentle cluster atmosphere heating during the early stages of the AGN's activity.
We explore the role that gravitational instability plays in the active galaxy NGC 1068, and its link with star formation. For this purpose, we use the Romeo-Falstad disc instability diagnostics and data from BIMA SONG, SDSS and SAURON. Our analysis illustrates that NGC 1068 is a gravitationally unstable `monster'. The inner R ~ 2 kpc are subject to strong disc instabilities, which power the vigorous starburst activity observed in this galaxy. AGN/stellar feedback tries to quench such instabilities from inside out by `decapitating' the surface density of molecular gas, but does not succeed. Gravitational instability `wins' because it is driven by the stars via their much higher surface density. This triggers local gravitational collapse/fragmentation in the molecular gas, and thus star formation.
We present a catalogue of dense cores in a $\sim 4^\circ\times2^\circ$ field of the Taurus star-forming region, inclusive of the L1495 cloud, derived from Herschel SPIRE and PACS observations in the 70 $\mu$m, 160 $\mu$m, 250 $\mu$m, 350 $\mu$m, and 500 $\mu$m continuum bands. Estimates of mean dust temperature and total mass are derived using modified blackbody fits to the spectral energy distributions. We detect 528 starless cores of which $\sim10$-20% are gravitationally bound and therefore presumably prestellar. Our census of unbound objects is $\sim85$% complete for $M>0.015\,M_\odot$ in low density regions ($A_V\stackrel{<}{_\sim}5$ mag), while the bound (prestellar) subset is $\sim85$% complete for $M>0.1\,M_\odot$ overall. The prestellar core mass function (CMF) is consistent with lognormal form, resembling the stellar system initial mass function, as has been reported previously. All of the inferred prestellar cores lie on filamentary structures whose column densities exceed the expected threshold for filamentary collapse, in agreement with previous reports. Unlike the prestellar CMF, the unbound starless CMF is not lognormal, but instead is consistent with a power-law form below $0.3\,M_\odot$ and shows no evidence for a low-mass turnover. It resembles previously reported mass distributions for CO clumps at low masses ($M\stackrel{<}{_\sim}0.3\,M_\odot$). The volume density PDF, however, is accurately lognormal except at high densities. It is consistent with the effects of self-gravity on magnetized supersonic turbulence. The only significant deviation from lognormality is a high-density tail which can be attributed unambiguously to prestellar cores.
Accurate estimation of the merger timescale of galaxy clusters is important to understand the cluster merger process and further the formation and evolution of the large-scale structure of the universe. In this paper, we explore a baryonic effect on the merger timescale of galaxy clusters by using hydrodynamical simulations. We find that the baryons play an important role in accelerating the merger process. The merger timescale decreases with increasing the gas fraction of galaxy clusters. For example, the merger timescale is shortened by a factor of up to 3 for merging clusters with gas fractions 0.15, compared with the timescale obtained with zero gas fractions. The baryonic effect is significant for a wide range of merger parameters and especially more significant for nearly head-on mergers and high merging velocities. The baryonic effect on the merger timescale of galaxy clusters is expected to have impacts on the structure formation in the universe, such as the cluster mass function and massive substructures in galaxy clusters, and a bias of "no-gas" may exist in the results obtained from the dark matter-only cosmological simulations.
The dynamical formation of stellar-mass black hole-black hole binaries has long been a promising source of gravitational waves for the Laser Interferometer Gravitational-Wave Observatory (LIGO). Mass segregation, gravitational focusing, and multibody dynamical interactions naturally increase the interaction rate between the most massive black holes in dense stellar systems, eventually leading them to merge. We find that dynamical interactions enhance the merger rate of black hole binaries with total mass M_tot roughly as ~M_tot^beta, with beta >~ 4. We find that this relation holds mostly independently of the initial mass function, but the exact value depends on the degree of mass segregation. The detection rate of such massive black hole binaries is only further enhanced by LIGO's greater sensitivity to massive black hole binaries with M_tot <~ 80 solar masses. We find that for power-law BH mass functions dN/dM ~ M^-alpha with alpha <~ 2, LIGO is most likely to detect black hole binaries with a mass twice that of the maximum initial black hole mass and a mass ratio near one. Repeated mergers of black holes inside the cluster result in about ~5% of mergers being observed between two and three times the maximum initial black hole mass. Using these relations, one may be able to invert the observed distribution to the initial mass function with multiple detections of merging black hole binaries.
We present our investigation into the long-term variability of the X-ray obscuration and optical-UV-X-ray continuum in the Seyfert 1 galaxy NGC 5548. In 2013 and 2014, the Swift observatory monitored NGC 5548 on average every day or two, with archival observations reaching back to 2005, totalling about 670 ks of observing time. Both broadband spectral modelling and temporal rms variability analysis are applied to the Swift data. We disentangle the variability caused by absorption, due to an obscuring weakly-ionised outflow near the disk, from variability of the intrinsic continuum components (the soft X-ray excess and the power-law) originating from the disk and its associated coronae. The spectral model that we apply to this extensive Swift data is the global model that we derived for NGC 5548 from analysis of the stacked spectra from our multi-satellite campaign of summer 2013 (including XMM-Newton, NuSTAR and HST). The results of our Swift study show that changes in the covering fraction of the obscurer is the primary and dominant cause of variability in the soft X-ray band on timescales of 10 days to ~ 5 months. The obscuring covering fraction of the X-ray source is found to range between 0.7 and nearly 1.0. The contribution of the soft excess component to the X-ray variability is often much less than that of the obscurer, but it becomes comparable when the optical-UV continuum flares up. We find that the soft excess is consistent with being the high-energy tail of the optical-UV continuum and can be explained by warm Comptonisation: up-scattering of the disk seed photons in a warm, optically thick corona as part of the inner disk. To this date, the Swift monitoring of NGC 5548 shows that the obscurer has been continuously present in our line of sight for at least 4 years (since at least February 2012).
We reconsider the hypothesis of a vast cometary reservoir surrounding the Solar System - the Oort cloud of comets - within the framework of Milgromian Dynamics (MD or MOND). For this purpose we build a numerical model of the cloud, assuming the theory of modified gravity QUMOND. Adopting popular pair $\nu(x)=[1-\exp(-x^{1/2})]^{-1}$ for the MD interpolating function and $a_{0}=1.2\times10^{-10}$ m s$^{-2}$ for the MD acceleration scale as a standard, we have found that the observationally inferred Milgromian cloud of comets is much more radially compact than its Newtonian counterpart. The comets of the Milgromian cloud stay away from the zone where the Galactic tide can torque their orbits significantly. However, this does not need to be an obstacle for injection of the comets into the inner Solar System as the External Field Effect (EFE) can induce significant change in perihelion distance during one revolution of a comet around the Sun. The efficiency of such injection is further increased by shift of the classical planetary barrier towards the Sun. Adopting constraints on different interpolating function families and revised value of $a_{0}$ as found recently by Hees et al. (2016), the aforementioned qualitative results no longer hold and it can be summarized that the Milgromian cloud is very similar to the Newtonian in its overall size, binding energies of comets and hence operation of the Jupiter-Saturn barrier. However, torquing of perihelia due to EFE still play a significant role in the inner parts of the cloud. Consequently Sedna-like orbits and orbits of large semi-major axis Centaurs are easily comprehensible in MD. In MD they both belong to the same population, just in different modes of their evolution.
Links to: arXiv, form interface, find, astro-ph, recent, 1602, contact, help (Access key information)
We present the line-of-sight (LOS) velocities for 13 distant main sequence Milky Way halo stars with published proper motions. The proper motions were measured using long baseline (5-7 years) multi-epoch HST/ACS photometry, and the LOS velocities were extracted from deep (5-6 hour integrations) Keck II/DEIMOS spectra. We estimate the parameters of the velocity ellipsoid of the stellar halo using a Markov chain Monte Carlo ensembler sampler method. The velocity second moments in the directions of the Galactic $(l,b,$ LOS) coordinate system are $\langle v^2_l \rangle^{1/2} = 138^{+43}_{-26}$ km/s, $\langle v^2_b \rangle^{1/2} = 88^{+28}_{-17}$ km/s, and $\langle v^2_{\rm{LOS}} \rangle^{1/2} = 91^{+27}_{-14}$ km/s. We use these ellipsoid parameters to constrain the velocity anisotropy of the stellar halo. Ours is the first measurement of the anisotropy parameter $\beta$ using 3D kinematics outside of the solar neighborhood. We find $\beta=-0.3^{+0.4}_{-0.9}$, consistent with isotropy and lower than solar neighborhood $\beta$ measurements by 2$\sigma$ ($\beta_{SN} \sim 0.5-0.7$). We identify two stars in our sample that are likely members of the known TriAnd substructure, and excluding these objects from our sample increases our estimate of the anisotropy to $\beta=0.1^{+0.4}_{-1.0}$, which is still lower than solar neighborhood measurements by $1\sigma$. The potential decrease in $\beta$ with Galactocentric radius is inconsistent with theoretical predictions, though consistent with recent observational studies, and may indicate the presence of large, shell-type structure (or structures) at $r \sim 25$ kpc. The methods described in this paper will be applied to a much larger sample of stars with 3D kinematics observed through the ongoing HALO7D program.
Outflows are invoked in co-evolutionary models to link the growth of SMBH and galaxies through feedback phenomena, and from the analysis of both galaxies and Active Galactic Nuclei (AGN) samples at z$\sim1-3$, it is becoming clear that powerful winds are quite common in AGN hosts. High-resolution and high S/N observations are needed in order to uncover the physical properties of the wind through kinematics analysis. We exploited VIMOS, SINFONI and Subaru/IRCS Adaptive Optics data to study the kinematics properties on the scale the host galaxy of XID5395, a luminous, X-ray obscured Starburst/Quasar merging system at z$\sim1.5$ detected in the XMM-COSMOS field, and associated with an extreme [O II] emitter (EW$\sim200$ \AA). We mapped, for the first time, at high resolution the kinematics of the [O III] and H$\alpha$ line complexes and linked them with the [O II] emission. The high spatial resolution achieved allowed us to resolve all the components of the SB-QSO system. Our analysis with a resolution of few kpc reveals complexities and asymmetries in and around the nucleus of XID5395. The velocity field measured via non parametric analysis reveals different kinematic components, with maximum blueshifted and redshifted velocities up to $\simeq1300$ km s$^{-1}$, not spatially coincident with the nuclear core. These extreme values of the observed velocities and the spatial location can be explained by the presence of fast moving material. We also spectroscopically confirm the presence of a merging system at the same redshift of the AGN host. We propose that EW as large as $>150$ \AA\ in X-ray selected AGN may be an efficient criterion to isolate objects associated to the short, transition phase of "feedback" in the AGN-galaxy co-evolutionary path, which will subsequently evolve in an unobscured QSO, as suggested from the different observational evidences we accumulated for XID5395.
We present a new method for assessing the intrinsic 3D shape of prestellar cores from molecular column densities. We have employed hydrodynamic simulations of contracting, isothermal cores considering three intrinsic geometries: spherical, cylindrical/filamentary and disk-like. We have coupled our hydrodynamic simulations with non-equilibrium chemistry. We find that a) when cores are observed very elongated (i.e. for aspect ratios $\le$ 0.15) the intrinsic 3D geometry can be probed by their 2D molecular emission maps, since these exhibit significant qualitative morphological differences between cylindrical and disk-like cores. Specifically, if a disk-like core is observed as a filamentary object in dust emission, then it will be observed as two parallel filaments in $\rm{N_2H^{+}}$; b) for cores with higher aspect ratios (i.e. 0.15 $\sim$ 0.9) we define a metric $\Delta$ that quantifies whether a molecular column density profile is centrally peaked, depressed or flat. We have identified one molecule ($\rm{CN}$) for which $\Delta$ as a function of the aspect ratio probes the 3D geometry of the core; and c) for cores with almost circular projections (i.e. for aspect ratios $\sim$ 1), we have identified three molecules ($\rm{OH}$, $\rm{CO}$ and $\rm{H_2CO}$) that can be used to probe the intrinsic 3D shape by close inspection of their molecular column density radial profiles. We alter the temperature and the cosmic-ray ionization rate and demonstrate that our method is robust against the choice of parameters.
Using three fiducial Nbody+SPH simulations, we follow the merging of two disk galaxies with a hot gaseous halo component each, and examine whether the merger remnant can be a spiral galaxy. The stellar progenitor disks are destroyed by violent relaxation during the merging and most of their stars form a classical bulge, while the remaining form a thick disk and its bar. A new stellar disk forms subsequently and gradually in the remnant from the gas accreted mainly from the halo. It is vertically thin and well extended in its equatorial plane. A bar starts forming before the disk is fully in place, contrary to what is assumed in idealised simulations of isolated bar-forming galaxies. It has morphological features such as ansae and boxy/peanut bulges. Stars of different ages populate different parts of the box/peanut. A disky pseudobulge forms also, so that by the end of the simulation, all three types of bulges coexist. The oldest stars are found in the classical bulge, followed by those of the thick disk, then by those in the thin disk. The youngest stars are in the spiral arms and the disky pseudobulge. The disk surface density profiles are of type II (exponential with downbending), and the circular velocity curves are flat and show that the disks are submaximum in these examples: two clearly so and one near-borderline between maximum and submaximum. On average, only roughly between 10 and 20% of the stellar mass is in the classical bulge of the final models, i.e. much less than in previous simulations.
We report a discovery of a 1.3 Mpc diffuse radio source with extremely steep spectrum fading radio structures in the vicinity of the Abell 449 cluster of galaxies. Its extended diffuse lobes are bright only at low radio frequencies and their synchrotron age is about 160 Myr. The parent galaxy of the extended relic structure, which is the dominant galaxy within the cluster, is starting a new jet activity. There are three weak X-rays sources in the vicinity of the cluster as found in the ROSAT survey, however it is not known if they are connected with this cluster of galaxies. Just a few radio galaxy relics are currently known in the literature, as finding them requires sensitive and high angular resolution low-frequency radio observations. Objects of this kind, which also are starting a new jet activity, are important for understanding the life cycle and evolution of active galactic nuclei. A new 613 MHz map as well as the archival radio data pertaining to this object are presented and analyzed.
Galactic plane surveys of pristine molecular clouds are key for establishing a Galactic-scale view of the earliest stages of star formation. For this reason Peretto & Fuller (2009) built an unbiased sample of IRDCs in the 10 deg < |l| < 65 deg, |b|<1 deg region of the Galactic plane using Spitzer 8micron extinction. However, in absorption studies, intrinsic fluctuations in the mid-infrared background can be mis-interpreted as foreground clouds. The main goal of the study presented here is to disentangle real clouds in the Spitzer Dark Cloud (SDC) catalogue from artefacts due to fluctuations in the mid-infrared background. We constructed H_2 column density maps at ~18" resolution using the 160micron and 250micron data from the Herschel Galactic plane survey Hi-GAL. We also developed an automated detection scheme that confirms the existence of a SDC through its association with a peak on these Herschel column density maps. Detection simulations, along with visual inspection of a small sub-sample of SDCs, have been performed to get better insight into the limitations of our automated identification scheme. Our analysis shows that 76(+/-19)% of the catalogued SDCs are real. This fraction drops to 55(+/-12)% for clouds with angular diameters larger than ~1 arcminute. The contamination of the PF09 catalogue by large spurious sources reflect the large uncertainties associated to the construction of the 8micron background emission, a key stage towards the identification of SDCs. A comparison of the Herschel confirmed SDC sample with the BGPS and ATLASGAL samples shows that SDCs probe a unique range of cloud properties, reaching down to more compact and lower column density clouds than any of these two (sub-)millimetre Galactic plane surveys.
The Draco Nebula is a high Galactic latitude interstellar cloud likely to have been formed by the collision of a Galactic halo cloud entering the disk of the Milky Way. Such conditions are ideal to study the formation of cold and dense gas in colliding flows of warm gas. We present Herschel-SPIRE observations that reveal the fragmented structure of the interface between the infalling cloud and the Galactic layer. This front is characterized by a Rayleigh-Taylor instability structure. From the determination of the typical length of the periodic structure (1.8 pc) we estimated the gas kinematic viscosity and the turbulence dissipation scale (0.1 pc) that is compatible with that expected if ambipolar diffusion is the main mechanism of energy dissipation. The small scale structures of the Nebula are typical of that seen in some molecular clouds. The gas density has a log-normal distribution with an average value of $10^3$ cm$^{-3}$. The size of the structures is 0.1-0.3 pc but this estimate is limited by the resolution of the observations. The mass ranges from 0.2 to 20 M$_{\odot}$ and the distribution of the more massive clumps follows a power law $dN/d\log(M) \sim M^{-1.6}$. We identify a mass-size relation with the same exponent as that found in GMCs ($M\sim L^{2.3}$) but only 15% of the mass of the cloud is in gravitationally bound structures. We conclude that the increase of pressure in the collision is strong enough to trigger the WNM-CNM transition caused by the interplay between turbulence and the thermal instability as self-gravity is not dominating the dynamics.
The radiative cooling of shocked gas with primordial chemical composition is an important process relevant to the formation of the first stars and structures, as well as taking place also in high velocity cloud collisions and supernovae explosions. Among the different processes that need to be considered, the formation kinetics and cooling of molecular hydrogen are of prime interest, since they provide the only way to lower the gas temperature to values well below $\sim$10$^4$~K. In previous works, the internal energy level structure of H$_2$ and its cation has been treated in the approximation of rovibrational ground state at low densities, or trying to describe the dynamics using some arbitrary $v>0$ H$_2$ level that is considered representative of the excited vibrational manifold. In this study, we compute the vibrationally resolved kinetics for the time-dependent chemical and thermal evolution of the post-shock gas in a medium of primordial composition. The calculated non-equilibrium distributions are used to evaluate effects on the cooling function of the gas and on the cooling time. Finally, we discuss the dependence of the results to different initial values of the shock velocity and redshift.
Context. The chemistry of the diffuse interstellar medium rests upon three pillars: exothermic ion-neutral reactions (" cold chemistry "), endothermic neutral-neutral reactions with significant activation barriers (" warm chemistry "), and reactions on the surfaces of dust grains. While warm chemistry becomes important in the shocks associated with turbulent dissipation regions, the main path for the formation of interstellar OH and H2O is that of cold chemistry. Aims. The aim of this study is to observationally confirm the association of atomic oxygen with both atomic and molecular gas phases, and to understand the measured abundances of OH and OH + as a function of the available reservoir of H2. Methods. We obtained absorption spectra of the ground states of OH, OH+ and OI with high-velocity resolution, with GREAT on-board SOFIA, and with the THz receiver at the APEX. We analyzed them along with ancillary spectra of HF and CH from HIFI. To deconvolve them from the hyperfine structure and to separate the blend that is due to various velocity components on the sightline, we fit model spectra consisting of an appropriate number of Gaussian profiles using a method combining simulated annealing with downhill simplex minimization. Together with HF and/or CH as a surrogate for H2, and HI $\lambda$21 cm data, the molecular hydrogen fraction f^N\_H2 = N(H 2)/(N(H) + 2N(H 2)) can be determined. We then investigated abundance ratios as a function of f^N\_H2. Results. The column density of OI is correlated at a high significance with the amount of available molecular and atomic hydrogen, with an atomic oxygen abundance of $3 \times 10 ^{-4}$ relative to H nuclei. While the velocities of the absorption features of OH and OH+ are loosely correlated and reflect the spiral arm crossings on the sightline, upon closer inspection they display an anticorrespondence. The arm-to-interarm density contrast is found to be higher in OH than in OH+. While both species can coexist, with a higher abundance in OH than in OH+, the latter is found less frequently in absence of OH than the other way around, which is a direct consequence of the rapid destruction of OH+ by dissociative recombination when not enough H2 is available. This conjecture has been substantiated by a comparison between the OH/OH+ ratio with f^N\_H2, showing a clear correlation. The hydrogen abstraction reaction chain OH+ (H2,H) H2O+ (H2,H)H3O+ is confirmed as the pathway for the production of OH and H 2 O. Our estimate of the branching ratio of the dissociative recombination of H3O+ to OH and H2O is confined within the interval of 84 to 91%, which matches laboratory measurements (74 to 83%). -- A correlation between the linewidths and column densities of OH+ features is found to be significant with a false-alarm probability below 5%. Such a correlation is predicted by models of interstellar MHD turbulence. For OH the same correlation is found to be insignificant because there are more narrow absorption features. Conclusions. While it is difficult to assess the contributions of warm neutral-neutral chemistry to the observed abundances, it seems fair to conclude that the predictions of cold ion-neutral chemistry match the abundance patterns we observed.
Transitional disks are protoplanetary disks with opacity gaps/cavities in their dust distribution, a feature that may be linked to planet formation. We perform Bayesian modeling of the three transitional disks SZ Cha, CS Cha and T25 including photometry from the Herschel Space Observatory to quantify the improvements added by these new data. We find disk dust masses between 2x10^-5 and 4x10^-4 Msun, and gap radii in the range of 7-18 AU, with uncertainties of ~ one order of magnitude and ~ 4 AU, respectively. Our results show that adding Herschel data can significantly improve these estimates with respect to mid-infrared data alone, which have roughly twice as large uncertainties on both disk mass and gap radius. We also find weak evidence for different density profiles with respect to full disks. These results open exciting new possibilities to study the distribution of disk masses for large samples of disks.
Aims. Study the connection between the masing disk and obscuring torus in
Seyfert 2 galaxies.
Methods. We present a uniform X-ray spectral analysis of the high energy
properties of 14 nearby megamaser Active Galactic Nuclei observed by NuSTAR. We
use a simple analytical model to localize the maser disk and understand its
connection with the torus by combining NuSTAR spectral parameters with
available physical quantities from VLBI mapping.
Results. Most of the sources analyzed are heavily obscured, showing a column
density in excess of $\sim 10^{23}$ cm$^{-2}$. In particular, $79\%$ are
Compton-thick ($N_{\rm H} > 1.5 \times 10^{24}$ cm$^{-2}$). Using column
densities measured by NuSTAR, with the assumption that the torus is the
extension of the maser disk, and further assuming a reasonable density profile,
the torus dimensions can be predicted. They are found to be consistent with
mid-IR interferometry parsec-scale observations of Circinus and NGC 1068. In
this picture, the maser disk is intimately connected to the inner part of the
torus. It is probably made of a large number of molecular clouds connecting the
torus and the outer part of the accretion disk, giving rise to a thin disk
rotating in most cases in Keplerian or sub-Keplerian motion. This toy model
explains the established close connection between water megamaser emission and
nuclear obscuration as a geometric effect.
We identify [Rb IV] 1.5973 and [Cd IV] 1.7204 micron emission lines in high-resolution (R=40,000) near-infrared spectra of the planetary nebulae (PNe) NGC 7027 and IC 5117, obtained with the IGRINS spectrometer on the 2.7-m telescope at McDonald Observatory. We also identify [Ge VI] 2.1930 $\mu$m in NGC 7027. Alternate identifications for these features are ruled out based on the absence of other multiplet members and/or transitions with the same upper levels. Ge, Rb, and Cd can be enriched in PNe by s-process nucleosynthesis during the asymptotic giant branch (AGB) stage of evolution. To determine ionic abundances, we calculate [Rb IV] collision strengths and use approximations for those of [Cd IV] and [Ge VI]. Our identification of [Rb IV] 1.5973 $\mu$m is supported by the agreement between Rb3+/H+ abundances found from this line and the 5759.55 A feature in NGC 7027. Elemental Rb, Cd, and Ge abundances are derived with ionization corrections based on similarities in ionization potential ranges between the detected ions and O and Ne ionization states. Our analysis indicates abundances 2-4 times solar for Rb and Cd in both nebulae. Ge is subsolar in NGC 7027, but its abundance is uncertain due to the large and uncertain ionization correction. The general consistency of the measured relative s-process enrichments with predictions from models appropriate for these PNe (2.0-2.5 M_sun, [Fe/H]= -0.37) demonstrates the potential of using PN compositions to test s-process nucleosynthesis models.
Stellar evolution models predict that internal mixing should cause some sodium overabundance at the surface of red giants more massive than ~ 1.5--2.0 Msun. The surface aluminium abundance should not be affected. Nevertheless, observational results disagree about the presence and/or the degree of the Na and Al overabundances. In addition, Galactic chemical evolution models adopting different stellar yields lead to quite different predictions for the behavior of [Na/Fe] and [Al/Fe] versus [Fe/H]. Overall, the observed trends of these abundances with metallicity are not well reproduced. We readdress both issues, using new Na and Al abundances determined within the Gaia-ESO Survey, using two samples: i) more than 600 dwarfs of the solar neighborhood and of open clusters and ii) low- and intermediate-mass clump giants in six open clusters. Abundances of Na in giants with mass below ~2.0 Msun, and of Al in giants below ~3.0 Msun, seem to be unaffected by internal mixing processes. For more massive giants, the Na overabundance increases with stellar mass. This trend agrees well with predictions of stellar evolutionary models. Chemical evolution models that are able to fit well the observed [Na/Fe] vs. [Fe/H] trend in solar neighborhood dwarfs can not simultaneously explain the run of [Al/Fe] with [Fe/H], and viceversa. The comparison with stellar ages is hampered by severe uncertainties. Indeed, reliable age estimates are available for only a half of the stars of the sample. We conclude that Al is underproduced by the models, except for stellar ages younger than about 7 Gyr. In addition, some significant source of late Na production seems to be missing in the models. Either current Na and Al yields are affected by large uncertainties, and/or some important Galactic source(s) of these elements has not been taken into account up to now. [abridged]
Infrared Herschel imaging and spectroscopic observations of the nebula M1-67 around the Wolf-Rayet star WR 124 have been obtained along with optical imaging observations. The infrared images reveal a clumpy dusty nebula that extends up to 1 pc. The comparison with the optical images shows that the ionized gas nebula coincides with the dust nebula, the dust and the gas being mixed together. A photodissociation region is revealed from the infrared spectroscopic analysis. The analysis of the infrared spectrum of the nebula, where forbidden emission lines of ionized elements were detected, showed that the nebula consists of mildly processed material with the calculated abundance number ratios being N/O = 1.0 +/- 0.5 and C/O = 0.46 +/- 0.27. Based on a radiative transfer model, the dust mass of the nebula was estimated to be 0.22 Msun with a population of large grains being necessary to reproduce the observations. The comparison of the mass-loss rate and the abundance ratios to theoretical models of stellar evolution led to the conclusion that the nebular ejection took place during a RSG/YSG evolutionary phase of a central star with an initial mass of 32 Msun.
Links to: arXiv, form interface, find, astro-ph, recent, 1602, contact, help (Access key information)
Simulating dwarf galaxy halos in a reionizing Universe puts severe constraints on the sub-grid model employed in the simulations. Using the same sub-grid model that works for simulations without a UV-background (UVB) results in gas poor galaxies that stop forming stars very early on, except for halos with high masses. This is in strong disagreement with observed galaxies, which are gas rich and star forming down to a much lower mass range. To resolve this discrepancy, we ran a large suite of isolated dwarf galaxy simulations to explore a wide variety of sub-grid models and parameters, including timing and strength of the UVB, strength of the stellar feedback, and metallicity dependent Pop III feedback. We compared these simulations to observed dwarf galaxies by means of the baryonic Tully-Fisher relation (BTFR), which links the baryonic content of a galaxy to the observationally determined strength of its gravitational potential. We found that the results are robust to changes in the UVB. The strength of the stellar feedback shifts the results on the BTFR, but does not help to form gas rich galaxies at late redshifts. Only by including Pop III feedback are we able to produce galaxies that lie on the observational BTFR and that have neutral gas and ongoing star formation at redshift zero.
We present initial results from the "Ponos" zoom-in numerical simulations of dark matter substructures in massive ellipticals. Two very highly resolved dark matter halos with $M_{\rm vir}=1.2\times 10^{13}$ $M_{\odot}$ and $M_{\rm vir}=6.5\times 10^{12}$ $M_{\odot}$ and different ("violent" vs. "quiescent") assembly histories have been simulated down to $z=0$ in a $\Lambda$CDM cosmology with a total of 921,651,914 and 408,377,544 particles, respectively. Within the virial radius, the total mass fraction in self-bound $M_{\rm sub}>10^6$ $M_{\odot}$ subhalos at the present epoch is 15% for the violent host and 16.5% for the quiescent one. At $z=0.7$, these fractions increase to 19 and 33%, respectively, as more recently accreted satellites are less prone to tidal destruction. In projection, the average fraction of surface mass density in substructure at a distance of $R/R_{\rm vir}=0.02$ ($\sim 5-10$ kpc) from the two halo centers ranges from 0.6% to $\gtrsim 2$%, significantly higher than measured in simulations of Milky Way-sized halos. The contribution of subhalos with $M_{\rm sub} < 10^9$ $M_{\odot}$ to the projected mass fraction is between one fifth and one third of the total, with the smallest share found in the quiescent host. We assess the impact of baryonic effects via twin, lower-resolution hydrodynamical simulations that include metallicity-dependent gas cooling, star formation, and a delayed-radiative-cooling scheme for supernova feedback. Baryonic contraction produces a super-isothermal total density profile and increases the number of massive subhalos in the inner regions of the main host. The host density profiles and projected subhalo mass fractions appear to be broadly consistent with observations of gravitational lenses.
Observations have revealed a relative paucity of red giant (RG) stars within the central 0.5pc in the Galactic Center (GC). Motivated by this finding we investigate the hypothesis that collisions of stars with a fragmenting accretion disk are responsible for the observed dearth of evolved stars. We use 3-dimensional hydrodynamic simulations to model a star with radius $10 R_{\odot}$ and mass $1 M_{\odot}$, representative of the missing population of RGs, colliding with high density clumps. We find that multiple collisions with clumps of column density $\gtrsim10^{8}\, {\rm g\,cm^{-2}}$ can strip a substantial fraction of the star's envelope and in principle render it invisible to observations. Simulations confirm that repeated impacts are particularly efficient in driving mass loss as partially stripped RGs expand and have increased cross sections for subsequent collisions. Because the envelope is unbound on account of the kinetic energy of the star, any significant amount of stripping of the RG population in the GC should be mirrored by a systematic decay of their orbits and possibly by their enhanced rotational velocity. To be viable, this scenario requires that the total mass of the fragmenting disk has been several orders of magnitude higher than that of the early type stars which now form the stellar disk in the GC.
In this Letter we demonstrate that the two claims of $z\simeq 0.03$ OVII K$\alpha$ absorption lines from Warm Hot Intergalactic Medium (WHIM) along the lines of sight to the blazars H~2356-309 (Buote et al., 2009; Fang et al., 2010) and Mkn~501 (Ren, Fang \& Buote, 2014) are likely misidentifications of the $z=0$ OII K$\beta$ line produced by a diffuse Low-Ionization Metal Medium in the Galaxy's Interstellar and Circum-Galactic mediums. We perform detailed modeling of all the available high signal-to-noise Chandra LETG and XMM-Newton RGS spectra of H 2356-309 and Mkn 501 and demonstrate that the $z\simeq 0.03$ WHIM absorption along these two sightlines is statistically not required. Our results, however, do not rule out a small contribution from the $z\simeq 0.03$ OVII K$\alpha$ absorber along the line of sight to H~2356-309. In our model the temperature of the putative $z = 0.031$ WHIM filament is T$= 3\times 10^5$ K and the OVII column density is N$_{OV II} \le 4\times 10^{15}$ cm$^{-2}$, twenty times smaller than the OVII column density previously reported, and now more consistent with the expectations from cosmological hydrodynamical simulations.
We investigate the virilization of the emission lines Hbeta and Mg II in the sample of 287 Type 1 Active Galactic Nuclei taken from the Sloan Digital Sky Survey database. We explore the connections between the intrinsic line shifts and full widths at different levels of maximal intensity. We found that: (i) Hbeta seems to be a good virial estimator of black hole masses, and an intrinsic redshift of Hbeta is dominantly caused by the gravitational effect, (ii) there is an anti-correlation between the redshift and width of the wings of the Mg II line, (iii) the broad Mg II line can be used as virial estimator only at 50% of the maximal intensity, while the widths and intrinsic shifts of the line wings can not be used for this purpose.
The detection of gravitationally redshifted optical emission lines has been reported just for a few active galaxies. In this paper we give a short overview of studies that analyzed or exploited the detection of the gravitational redshift in optical AGN spectra. In addition, we tested the consistency of gravitational redshift as the physical origin of the redward shifts observed in their spectra using a sample of $\approx$ 50 Hamburg-ESO intermediate to high redshift quasars that are among the most luminous quasars known ($10^{47} \lesssim L \lesssim 10^{48}$ erg/s), and are expected to host very massive black holes. To this aim we modeled the line profile with accretion disk models around a black hole.
We perform a series of smoothed particle hydrodynamics simulations of isolated dwarf galaxies to compare different metal mixing models. In particular, we examine the role of diffusion in the production of enriched outflows, and in determining the metallicity distributions of gas and stars. We investigate different diffusion strengths, by changing the pre-factor of the diffusion coefficient, by varying how the diffusion coefficient is calculated from the local velocity distribution, and by varying whether the speed of sound is included as a velocity term. Stronger diffusion produces a tighter [O/Fe]-[Fe/H] distribution in the gas, and cuts off the gas metallicity distribution function at lower metallicities. Diffusion suppresses the formation of low-metallicity stars, even with weak diffusion, and also strips metals from enriched outflows. This produces a remarkably tight correlation between "metal mass-loading" (mean metal outflow rate divided by mean metal production rate) and the strength of diffusion, even when the diffusion coefficient is calculated in different ways. The effectiveness of outflows at removing metals from dwarf galaxies and the metal distribution of the gas is thus dependent on the strength of diffusion. By contrast, we show that the metallicities of stars are not strongly dependent on the strength of diffusion, provided that some diffusion is present.
In this paper, we study the accretion process for fluids flowing near a black hole in the context of $f(T)$ teleparallel gravity. Specifically, by performing a dynamical analysis by a Hamiltonian system, we are able to find the sonic points. After that, we consider different isothermal test fluids in order to study the accretion process when they are falling onto the black hole. We found that these flows can be classified according to the equation of state and the black hole features. Results are compared in $f(T)$ and $f(R)$ gravity.
Some reddened quasars appear to be transitional objects in the merger-induced black hole growth/galaxy evolution paradigm, where a heavily obscured nucleus starts to be unveiled by powerful quasar winds evacuating the surrounding cocoon of dust and gas. Hard X-ray observations are able to peer through this gas and dust, revealing the properties of circumnuclear obscuration. Here, we present NuSTAR and XMM-Newton/Chandra observations of FIRST-2MASS selected red quasars F2M 0830+3759 and F2M 1227+3214. We find that though F2M 0830+3759 is moderately obscured ($N_{\rm H,Z} = 2.1\pm0.2 \times10^{22}$ cm$^{-2}$) and F2M 1227+3214 is mildly absorbed ($N_{\rm H,Z} = 3.4^{+0.8}_{-0.7}\times10^{21}$ cm$^{-2}$) along the line-of-sight, heavier global obscuration may be present in both sources, with $N_{\rm H,S} = 3.7^{+4.1}_{-2.6} \times 10^{23}$ cm$^{-2}$ and $< 5.5\times10^{23}$ cm$^{-2}$, for F2M 0830+3759 and F2M 1227+3214, respectively. F2M 0830+3759 also has an excess of soft X-ray emission below 1 keV which is well accommodated by a model where 7% of the intrinsic AGN X-ray emission is scattered into the line-of-sight. While F2M 1227+3214 has a dust-to-gas ratio ($E(B-V)$/$N_{\rm H}$) consistent with the Galactic value, the $E(B-V)$/$N_{\rm H}$ value for F2M 0830+3759 is lower than the Galactic standard, consistent with the paradigm that the dust resides on galactic scales while the X-ray reprocessing gas originates within the dust-sublimation zone of the broad-line-region. The X-ray and 6.1$\mu$m luminosities of these red quasars are consistent with the empirical relations derived for high-luminosity, unobscured quasars, extending the parameter space of obscured AGN previously observed by NuSTAR to higher luminosities.
In this paper we present new AAOmega spectroscopy of 254 galaxies within a 30' radius around Abell 3888. We combine these data with the existing redshifts measured in a one degree radius around the cluster and performed a substructure analysis. We confirm 71 member galaxies within the core of A3888 and determine a new average redshift and velocity dispersion for the cluster of 0.1535 +\- 0.0009 and 1181 +\- 197 km/s, respectively. The cluster is elongated along an East-West axis and we find the core is bimodal along this axis with two sub-groups of 26 and 41 members detected. Our results suggest that A3888 is a merging system putting to rest the previous conjecture about the morphological status of the cluster derived from X-ray observations. In addition to the results on A3888 we also present six newly detected galaxy over-densities in the field, three of which we classify as new galaxy clusters.
Using the Binary Population and Spectral Synthesis code BPASS, we have calculated the rates, timescales and mass distributions for binary black hole mergers as a function of metallicity. We consider these in the context of the recently reported 1st LIGO event detection. We find that the event has a low probability of arising from a stellar population with initial metallicity mass fraction above $Z=0.010$. Binary black hole merger events with the reported masses are most likely in populations between Z=0.0001 and 0.002 (Z < 0.1Z_sun). The masses inferred for the black holes in the binary progenitor of GW 150914 are close to the predicted peak in the mass distribution for such events. We discuss the implications of our analysis for the electromagnetic follow-up of future LIGO event detections.
Laser Interferometer Gravitational-Wave Observatory, LIGO, found direct evidence of double black hole binaries emitting gravitational waves. Galactic nuclei are expected to harbor the densest population of stellar-mass black holes, accounting for as much as ~2% of the mass of the nuclear stellar cluster. A significant fraction (~30%) of these black holes can reside in binaries. We examine the fate of the black hole binaries in active galactic nuclei, which get trapped in the inner region of the accretion disk around the central supermassive black hole. We show that binary black holes can migrate into and then rapidly merge within the disk well within a Salpeter time. The binaries may also accrete a significant amount of gas from the disk, well above the Eddington rate. This could lead to detectable X-ray or gamma-ray emission, but would require hyper-Eddington accretion with a few % radiative efficiency, comparable to thin disks. We discuss implications for gravitational wave observations and black hole population studies.
Links to: arXiv, form interface, find, astro-ph, recent, 1602, contact, help (Access key information)