A near pristine atomic cooling halo close to a star forming galaxy offers a natural pathway for forming massive direct collapse black hole (DCBH) seeds which could be the progenitors of the $z>6$ redshift quasars. The close proximity of the haloes enables a sufficient Lyman-Werner flux to effectively dissociate H$_2$ in the core of the atomic cooling halo. A mild background may also be required to delay star formation in the atomic cooling halo, often attributed to distant background galaxies. In this letter we investigate the impact of metal enrichment from both the background galaxies and the close star forming galaxy under extremely unfavourable conditions such as instantaneous metal mixing. We find that within the time window of DCBH formation, the level of enrichment never exceeds the critical threshold (Z$_{cr} \sim 1 \times 10^{-5} \ \rm Z_{\odot})$, and attains a maximum metallicity of Z $\sim 2 \times 10^{-6} \ \rm Z_{\odot}$. As the system evolves, the metallicity eventually exceeds the critical threshold, long after the DCBH has formed.
Similarly to the cosmic star formation history, the black hole accretion rate density of the Universe peaked at 1<z<3. This cosmic epoch is hence best suited for investigating the effects of radiative feedback from AGN. Observational efforts are underway to quantify the impact of AGN feedback, if any, on their host galaxies. Here we present a study of the molecular gas content of AGN hosts at z~1.5 using CO[2-1] line emission observed with ALMA for a sample of 10 AGNs. We compare this with a sample of galaxies without an AGN matched in redshift, stellar mass, and star formation rate. We detect CO in 3 AGNs with $\mathrm{L_{CO} \sim 6.3-25.1\times 10^{9} L_{\odot}}$ which translates to a molecular hydrogen gas mass of $\mathrm{2.5-10\times 10^{10} M_{\odot}}$ assuming conventional conversion factor of $\mathrm{\alpha_{CO}}\sim3.6$. Our results indicate a >99% probability of lower depletion time scales and lower molecular gas fractions in AGN hosts with respect to the non-AGN comparison sample. We discuss the implications of these observations on the impact that AGN feedback may have on star formation efficiency of z>1 galaxies.
We propose an innovative method for measuring the neutral hydrogen (HI) content of an optically-selected spectroscopic sample of galaxies through cross-correlation with HI intensity mapping measurements. We show that the HI-galaxy cross-power spectrum contains an additive shot noise term which scales with the average HI brightness temperature of the optically-selected galaxies, allowing constraints to be placed on the average HI mass per galaxy. This approach can estimate the HI content of populations too faint to directly observe through their 21cm emission over a wide range of redshifts. This cross-correlation, as a function of optical luminosity or colour, can be used to derive HI-scaling relations. We demonstrate that this signal will be detectable by cross-correlating upcoming Australian SKA Pathfinder (ASKAP) observations with existing optically-selected samples. We also use semi-analytic simulations to verify that the HI mass can be successfully recovered by our technique in the range M_HI > 10^8 M_solar, in a manner independent of the underlying power spectrum shape. We conclude that this method is a powerful tool to study galaxy evolution, which only requires a single intensity mapping dataset to infer complementary HI gas information from existing optical and infra-red observations.
CO--0.40--0.22 is a peculiar molecular cloud that is compact and has an extraordinary broad velocity width. It is found in the central molecular zone of our Galaxy. In this direction, there is another cloud with an H$_2$O maser spot at a lower velocity. Collision with this low-velocity cloud could be responsible for the broad velocity width of CO--0.40--0.22. We performed phase-referencing VLBI astrometry with VERA and detected the annual parallax of the H$_2$O maser spot in the low-velocity cloud to be $0.33 \pm 0.14$ mas, which corresponds to a distance of $3.07^{+2.22}_{-0.91}$ kpc from the Sun. This implies that the low-velocity cloud is located in the Galactic disk on the near side of the central molecular zone.
We calculated the polarization degree of hydrogen Balmer broad emission lines from a number of active galactic nuclei (AGNs) with determined virial factors. The objects were selected from the sample presented by Decarli et al.(2008). In our calculations, we used the model of the flattened disc-like structure of the broad-line emission region (BLR). In this model, the expression for the virial factor makes it possible to determine the inclination angle for the flattened BLR, which in turn yields the polarization degree of the broad emission lines. As a result, we obtained the direct relation between the polarization degree and the virial factor. We also compared the determined values of the polarization degree with those obtained in polarimetric observations.
We study the properties of dark matter halos that contain star-forming galaxies at $1.43 \le z \le 1.74$ using the FMOS-COSMOS survey. The sample consists of 516 objects with a detection of the H$\alpha$ emission line, that represent the star-forming population at this epoch having a stellar mass range of $10^{9.57}\le M_\ast/M_\odot \lesssim 10^{11.4}$ and a star formation rate range of $15\lesssim \mathrm{SFR}/(M_\odot \mathrm{yr^{-1}}) \lesssim 600$. We measure the projected two-point correlation function while carefully taking into account observational biases, and find a significant clustering amplitude at scales of $0.04$-$10~h^{-1}~\mathrm{cMpc}$, with a correlation length $r_0 = 5.21^{+0.70}_{-0.67}~h^{-1}~\mathrm{cMpc}$ and a bias $b=2.59^{+0.41}_{-0.34}$. We interpret our clustering measurement using a halo occupation distribution model. The sample galaxies appear to reside in halos with mass $M_\mathrm{h} = 4.6^{+1.1}_{-1.6}\times10^{12}~h^{-1}M_\odot$ on average that will likely become present-day halos of mass $M_\mathrm{h} (z=0) \sim2\times10^{13}~h^{-1}M_\odot$, equivalent to the typical halo mass scale of galaxy groups. We then confirm the decline of the stellar-to-halo mass ratio at $M_\mathrm{h}<10^{12}~M_\odot$, finding $M_\ast/M_\mathrm{h} \approx 5\times10^{-3}$ at $M_\mathrm{h}=10^{11.86}~M_\odot$, which is lower by a factor of 2-4 than those measured at higher masses. Finally, we use our results to illustrate the future capabilities of Subaru's Prime-Focus Spectrograph, a next-generation instrument that will provide strong constraints on the galaxy-formation scenario by obtaining precise measurements of galaxy clustering at $z>1$.
A conclusive model for the formation of dwarf spheroidal (dSph) galaxies still remains elusive. Owing to their proximity to the massive spirals Milky Way (MW) and M31, various environmental processes have been invoked to explain their origin. In this context, the tidal stirring model postulates that interactions with MW-sized hosts can transform rotationally supported dwarfs, resembling present-day dwarf irregular (dIrr) galaxies, into systems with the kinematic and structural properties of dSphs. Using N-body+SPH simulations, we investigate the dependence of this transformation mechanism on the gas fraction, fgas, in the disk of the progenitor dwarf. Our numerical experiments incorporate for the first time the combined effects of radiative cooling, ram-pressure stripping, star formation, supernova (SN) winds, and a cosmic UV background. For a given orbit inside the primary galaxy, rotationally supported dwarfs with gas fractions akin to those of observed dIrrs (fgas >= 0.5), demonstrate a substantially enhanced likelihood and efficiency of transformation into dSphs relative to their collisionless (fgas = 0) counterparts. We argue that the combination of ram-pressure stripping and SN winds causes the gas-rich dwarfs to respond more impulsively to tides, augmenting their transformation. When fgas >= 0.5, disky dwarfs on previously unfavorable low-eccentricity or large-pericenter orbits are still able to transform. On the widest orbits, the transformation is incomplete; the dwarfs retain significant rotational support, a relatively flat shape, and some gas, naturally resembling transition-type systems. We conclude that tidal stirring constitutes a prevalent evolutionary mechanism for shaping the structure of dwarf galaxies within the currently favored CDM cosmological paradigm.
The archival data of 3C 345, a type 1 quasar at $z = 0.5928$, obtained with Suzaku and Swift/BAT are analysed. Though previous studies of this source applied only a simple broken power law model, a heavily obscuring material is found to be required by considering Akaike information criteria. The application of the numerical torus model by Murphy & Yaqoob (2009) surprisingly reveals the existence of Compton thick type 2 nucleus with the line-of-sight hydrogen column density of the torus of $N_{\rm H} = 10^{24.5}$ cm$^{-2}$ and the inclination angle of $\theta_{\rm inc} = 90^{\circ}$. However, this model fails to account for the Eddington ratio obtained with the optical observations by Gu et al. (2001) and Shen et al. (2011), or requires the existence of a supermassive black hole binary, which was suggested by Lobanov & Roland (2005), thus this model is likely to be inappropriate for 3C 345. A partial covering ionized absorber model which accounts for absorption in "hard excess" type 1 AGNs is also applied, and finds a Compton thick absorber with the column density of $N_{\rm H} \simeq 10^{25}$ cm$^{-2}$, the ionization parameter of $\log \xi \gtrsim 2$, and the covering fraction of $75\% \lesssim f_{c} \lesssim 85\%$. Since this model obtains a black hole mass of $\log ( M_{\rm BH} / M_{\odot} ) = 9.8$, which is consistent with the optical observation by Gu et al. (2001), this model is likely to be the best-fitting model of this source. The results suggest that 3C 345 is the most distant and most obscured hard excess AGN at this time.
(Abridged) Aims: We aim to use the progressive heating of the gas caused by the feedback of high-mass young stellar objects (YSOs) to prove the statistical validity of the most common schemes used to define an evolutionary sequence for high-mass clumps, and characterise the sensitivity of different tracers to this process. Methods: From the spectroscopic follow-ups of the ATLASGAL TOP100 sample, we selected several multiplets of CH3CN, CH3CCH, and CH3OH emission lines to derive and compare the physical properties of the gas in the clumps along the evolutionary sequence. Our findings are compared with results obtained from CO isotopologues, dust, and NH3 from previous studies on the same sample. Results: The chemical properties of each species have a major role on the measured physical properties. Low temperatures are traced by NH3, CH3OH, and CO (in the early phases), the warm and dense envelope can be probed with CH3CN, CH3CCH, and, in evolved sources via CO isotopologues. CH3OH and CH3CN are also abundant in the hot cores, and their high-excitation transitions may be good tools to study the kinematics in the hot gas surrounding the YSOs that these clumps are hosting. All tracers show, to different degrees, progressive warming with evolution. The relation between gas temperature and L/M is reproduced by a toy model of a spherical, internally heated clump. Conclusions: The evolutionary sequence defined for the clumps is statistically valid and we could identify the processes dominating in different intervals of L/M. For L/M<2Lsun/Msun a large quantity of gas is still being accumulated and compressed at the bottom of the potential well. Between 2Lsun/Msun<L/M<40Lsun/Msun the YSOs gain mass and increase in L; the first hot cores appear around L/M=10Lsun/Msun. Finally, for L/M>40Lsun/Msun HII regions become common, showing that dissipation of the parental clump dominates.
In this paper, we study the 3-body products (two single stars and a binary) of binary-binary (2+2) scattering interactions. This is done using a combination of analytic methods and numerical simulations of 2+2 scattering interactions, both in isolation and in a homogeneous background potential. We derive analytically a simple formula relating the angle between the velocity vectors of the two ejected single stars and the orbital separation of the remaining binary. We compare our analytic formulation to numerical scattering simulations, and illustrate that the agreement is excellent, both in isolation and in a homogeneous background potential. Our results are ideally suited for application to the GAIA database, which is expected to identify many hundred runaway stars. The analytic relation presented here has the potential to identify runaway stars formed dynamically with high confidence. Finally, by applying our method to the runaways AE Aur and $\mu$ Col, we illustrate that it can be used to constrain the history of the background potential, which was denser than the presently observed density in the case of the Trapezium cluster.
We report on SPT-CLJ2011-5228, a giant system of arcs created by a cluster at $z=1.06$. The arc system is notable for the presence of a bright central image. The source is a Lyman Break galaxy at $z_s=2.39$ and the mass enclosed within the 14 arc second radius Einstein ring is $10^{14.2}$ solar masses. We perform a full light profile reconstruction of the lensed images to precisely infer the parameters of the mass distribution. The brightness of the central image demands that the central total density profile of the lens be shallow. By fitting the dark matter as a generalized Navarro-Frenk-White profile---with a free parameter for the inner density slope---we find that the break radius is $270^{+48}_{-76}$ kpc, and that the inner density falls with radius to the power $-0.38\pm0.04$ at 68 percent confidence. Such a shallow profile is in strong tension with our understanding of relaxed cold dark matter halos; dark matter only simulations predict the inner density should fall as $r^{-1}$. The tension can be alleviated if this cluster is in fact a merger; a two halo model can also reconstruct the data, with both clumps (density going as $r^{-0.8}$ and $r^{-1.0}$) much more consistent with predictions from dark matter only simulations. At the resolution of our Dark Energy Survey imaging, we are unable to choose between these two models, but we make predictions for forthcoming Hubble Space Telescope imaging that will decisively distinguish between them.
Links to: arXiv, form interface, find, astro-ph, recent, 1703, contact, help (Access key information)
We present analyses of the spatial distributions of stars in the young (1 - 3
Myr) star-forming regions IC348 and NGC1333 in the Perseus Giant Molecular
Cloud. We quantify the spatial structure using the $\mathcal{Q}$-parameter and
find that both IC348 and NGC1333 are smooth and centrally concentrated with
$\mathcal{Q}$-parameters of 0.98 and 0.89 respectively. Neither region exhibits
mass segregation ($\Lambda_{\rm MSR} = 1.1^{+0.2}_{-0.3}$ for IC348 and
$\Lambda_{\rm MSR} = 1.2^{+0.4}_{-0.3}$ for NGC1333, where $\Lambda_{\rm MSR}
\sim 1$ corresponds to no mass segregation), nor do the most massive stars
reside in areas of enhanced stellar surface density compared to the average
surface density, according to the $\Sigma_{\rm LDR}$ method.
We then constrain the dynamical histories and hence initial conditions of
both regions by comparing the observed values to $N$-body simulations at
appropriate ages. Stars in both regions likely formed with sub-virial
velocities which contributed to merging of substructure and the formation of
smooth clusters. The initial stellar densities were no higher than $\rho \sim
100 - 500$M$_\odot$pc$^{-3}$ for IC348 and $\rho \sim 500 -
2000$M$_\odot$pc$^{-3}$ for NGC1333. These initial densities, in particular
that of NGC1333, are high enough to facilitate dynamical interactions which
would likely affect $\sim$10 per cent of protoplanetary discs and binary stars.
We study the effect of density fluctuations induced by turbulence on the HI/H$_2$ structure in photodissociation regions (PDRs) both analytically and numerically. We perform magnetohydrodynamic numerical simulations for both subsonic and supersonic turbulent gas and chemical HI/H$_2$ balance calculations. We derive atomic-to-molecular density profiles and the HI column density probability density function (PDF). We find that while the HI/H$_2$ density profiles are strongly perturbed in turbulent gas, the mean HI column density is well approximated by the uniform-density analytic formula of Sternberg et al. (2014). The important effect of supersonic density fluctuations is to produce a wide HI column density PDF which depends on (a) the radiation intensity to mean density ratio, (b) the sonic Mach number and (c) the turbulence decorrelation scale, or driving scale. We derive an analytic model for the HI PDF and demonstrate how our model, combined with 21 cm observations, can be used to constrain the Mach number and driving scale of turbulent gas. We test our model using observations of HI in the Perseus GMC. We find that the narrow observed HI PDF may imply small scale decorrelation, pointing to the potential importance of subcloud-scale turbulence driving in Perseus.
Galaxy clusters are known to induce gas loss in infalling galaxies due to the ram pressure exerted by the intracluster medium over their gas content. In this paper, we investigate this process through a set of simulations of Milky Way like galaxies falling inside idealised clusters of 10$^{14}$ M$_\odot$ and 10$^{15}$ M$_\odot$, containing a cool-core or not, using the adaptive mesh refinement code RAMSES. We use these simulations to constrain how much of the initial mass contained in the gaseous disk of the galaxy will be converted into stars and how much of it will be lost, after a single crossing of the entire cluster. We find that, if the galaxy reaches the central region of a cool-core cluster, it is expected to lose all its gas, independently of its entry conditions and of the cluster's mass. On the other hand, it is expected to never lose all its gas after crossing a cluster without a cool-core just once. Before reaching the centre of the cluster, the SFR of the galaxy is always enhanced, by a factor of 1.5 to 3. If the galaxy crosses the cluster without being completely stripped, its final amount of gas is on average two times smaller after crossing the 10$^{15}$ M$_\odot$ cluster, relative to the 10$^{14}$ M$_\odot$ cluster. This is reflected in the final SFR of the galaxy, which is also two times smaller in the former, ranging from 0.5 $-$ 1 M$_\odot$ yr$^{-1}$, compared to 1 $-$ 2 M$_\odot$ yr$^{-1}$ for the latter.
We study the formation of runaway stars due to binary-binary (2+2) interactions in young star-forming clusters and/or associations. This is done using a combination of analytic methods and numerical simulations of 2+2 scattering interactions, both in isolation and in a homogeneous background potential. We focus on interactions that produce two single stars and a binary, and study the outcomes as a function of the depth of the background potential, within a range typical of cluster cores. As reference parameters for the observational properties, we use those observed for the system of runaway stars AE Aur and $\mu$ Col and binary $\iota$ Ori. We find that the outcome fractions have no appreciable dependence on the depth of the potential, and neither do the velocities of the ejected single stars. However, as the potential gets deeper and a larger fraction of binaries remain trapped, two binary populations emerge, with the escaped component having higher speeds and shorter semi-major axes than the trapped one. Additionally, we find that the relative angles between the ejected products are generally large. In particular, the angle between the ejected fastest star and the escaped binary is typically $\gtrsim 120-135^{\circ}$, with a peak at around $160^{\circ}$. However, as the potential gets deeper, the angle distribution becomes broader. Finally, we discuss the implications of our results for the interpretation of the properties of the runaway stars AE Aur and $\mu$ Col.
Could there be intermediate mass black holes in essentially all old dwarf galaxies? I argue that current observations of Active Galactic Nuclei in dwarfs allow such a radical hypothesis which provides early feedback and potentially provides a unifying explanation for many if not all of the apparent dwarf galaxy anomalies, such as the abundance, core-cusp, "too big to fail", ultra-faint and baryon-fraction issues. I describe the supporting arguments, which are largely circumstantial in nature, and discuss a number of tests. There is no strong motivation for modifying the nature of cold dark matter in order to explain any of the dwarf galaxy "problems".
We analyze the environmental properties of 370 present-day early-type galaxies in the MASSIVE and ATLAS3D surveys, two complementary volume-limited integral-field spectroscopic (IFS) galaxy surveys spanning absolute $K$-band magnitude $-21.5 > M_K > -26.6$, or stellar mass $6\ times 10^{9} < M_* < 2 \times 10^{12} M_\odot$. We find these galaxies to reside in a diverse range of environments measured by four methods: group membership (whether a galaxy is a brightest group/cluster galaxy, satellite, or isolated), halo mass, large-scale mass density (measured over a few Mpc), and local mass density (measured within the $N$th neighbor). The spatially resolved IFS stellar kinematics provide robust measurements of the angular momentum parameter $\lambda_e$ and enable us to examine the relationship among $\lambda_e$, stellar mass, and environment of ETGs. We find a strong correlation between $\lambda_e$ and stellar mass, where the average $\lambda_e$ decreases from $\sim 0.5$ to less than 0.1 with increasing galaxy mass, and the fraction of slow rotators increases correspondingly from $\sim 10$% to 90%. While we see weak trends between environment and both $\lambda_e$ and slow rotator fraction, this kinematic morphology-density relation is fully accounted for by the strong correlation between the slow rotator fraction and stellar mass; due to the correlation between stellar mass and environment, it naturally results in a correlation between slow/fast classification and environment as well. A possible exception is that the increased fraction of slow rotators at high local density is slightly more than expected based only on these joint correlations. Our results suggest that the physical processes responsible for building up the present-day stellar masses of massive galaxies are also very efficient at reducing their angular momentum, in any environment.
Recent observations have revealed massive galactic molecular outflows that may have physical conditions (high gas densities) required to form stars. Indeed, several recent models predict that such massive galactic outflows may ignite star formation within the outflow itself. This star-formation mode, in which stars form with high radial velocities, could contribute to the morphological evolution of galaxies, to the evolution in size and velocity dispersion of the spheroidal component of galaxies, and would contribute to the population of high-velocity stars, which could even escape the galaxy. Such star formation could provide in-situ chemical enrichment of the circumgalactic and intergalactic medium (through supernova explosions of young stars on large orbits), and some models also predict that it may contribute substantially to the global star formation rate observed in distant galaxies. Although there exists observational evidence for star formation triggered by outflows or jets into their host galaxy, as a consequence of gas compression, evidence for star formation occurring within galactic outflows is still missing. Here we report new spectroscopic observations that unambiguously reveal star formation occurring in a galactic outflow at a redshift of 0.0448. The inferred star formation rate in the outflow is larger than 15 Msun/yr. Star formation may also be occurring in other galactic outflows, but may have been missed by previous observations owing to the lack of adequate diagnostics.
We explore the effects of the expected higher cosmic ray (CR) ionization rates $\zeta_{\rm CR}$ on the abundances of carbon monoxide (CO), atomic carbon (C), and ionized carbon (C$^+$) in the H$_2$ clouds of star-forming galaxies. The study of Bisbas et al. (2015) is expanded by: a) using realistic inhomogeneous Giant Molecular Cloud (GMC) structures, b) a detailed chemical analysis behind the CR-induced destruction of CO, and c) exploring the thermal state of CR-irradiated molecular gas. CRs permeating the interstellar medium with $\zeta_{\rm CR}$$\gtrsim 10\times$(Galactic) are found to significantly reduce the [CO]/[H$_2$] abundance ratios throughout the mass of a GMC. CO rotational line imaging will then show much clumpier structures than the actual ones. For $\zeta_{\rm CR}$$\gtrsim 100\times$(Galactic) this bias becomes severe, limiting the utility of CO lines for recovering structural and dynamical characteristics of H$_2$-rich galaxies throughout the Universe, including many of the so-called Main Sequence (MS) galaxies where the bulk of cosmic star formation occurs. Both C$^+$ and C abundances increase with rising $\zeta_{\rm CR}$, with C remaining the most abundant of the two throughout H$_2$ clouds, when $\zeta_{\rm CR}\sim (1-100)\times$(Galactic). C$^+$ starts to dominate for $\zeta_{\rm CR}$$\gtrsim 10^3\times$(Galactic). The thermal state of the gas in the inner and denser regions of GMCs is invariant with $T_{\rm gas}\sim 10\,{\rm K}$ for $\zeta_{\rm CR}\sim (1-10)\times$(Galactic). For $\zeta_{\rm CR}$$\sim 10^3\times$(Galactic) this is no longer the case and $T_{\rm gas}\sim 30-50\,{\rm K}$ are reached. Finally we identify OH as the key species whose $T_{\rm gas}-$sensitive abundance could mitigate the destruction of CO at high temperatures.
HOCO$^+$ is a polar molecule that represents a useful proxy for its parent molecule CO$_2$, which is not directly observable in the cold interstellar medium. This cation has been detected towards several lines of sight, including massive star forming regions, protostars, and cold cores. Despite the obvious astrochemical relevance, protonated CO$_2$ and its deuterated variant, DOCO$^+$, still lack an accurate spectroscopic characterisation. The aim of this work is to extend the study of the ground-state pure rotational spectra of HOCO$^+$ and DOCO$^+$ well into the sub-millimetre region. Ground-state transitions have been recorded in the laboratory using a frequency-modulation absorption spectrometer equipped with a free-space glow-discharge cell. The ions were produced in a low-density, magnetically-confined plasma generated in a suitable gas mixture. The ground-state spectra of HOCO$^+$ and DOCO$^+$ have been investigated in the 213-967 GHz frequency range, with the detection of 94 new rotational transitions. Additionally, 46 line positions taken from the literature have been accurately remeasured. The newly-measured lines have significantly enlarged the available data sets for HOCO$^+$ and DOCO$^+$, thus enabling the determination of highly accurate rotational and centrifugal distortion parameters. Our analysis showed that all HOCO$^+$ lines with Ka $\geq$ 3 are perturbed by a ro-vibrational interaction that couples the ground state with the v$_ 5$ = 1 vibrationally-excited state. This resonance has been explicitly treated in the analysis in order to obtain molecular constants with clear physical meaning. The improved sets of spectroscopic parameters provide enhanced lists of very accurate, sub-millimetre rest-frequencies of HOCO$^+$ and DOCO$^+$ for astrophysical applications. These new data challenges a recent tentative identification of DOCO$^+$ toward a pre-stellar core.
Turbulence models attempt to account for unresolved dynamics and diffusion in hydrodynamical simulations. We develop a common framework for two-equation Reynolds-Averaged Navier-Stokes (RANS) turbulence models, and we implement six models in the Athena code. We verify each implementation with the standard subsonic mixing layer, although the level of agreement depends on the definition of the mixing layer width. We then test the validity of each model into the supersonic regime, showing that compressibility corrections can improve agreement with experiment. For models with buoyancy effects, we also verify our implementation via the growth of the Rayleigh-Taylor instability in a stratified medium. The models are then applied to the ubiquitous astrophysical shock-cloud interaction in three dimensions. We focus on the mixing of shock and cloud material, comparing results from turbulence models to high-resolution simulations (up to 200 cells per cloud radius) and ensemble-averaged simulations. We find that the turbulence models lead to increased spreading and mixing of the cloud, although no two models predict the same result. Increased mixing is also observed in inviscid simulations at resolutions greater than 100 cells per radius, which suggests that the turbulent mixing begins to be resolved.
Follow-up observations at high-angular resolution of submillimeter galaxies showed that the single-dish sources are comprised of a blend of several galaxies. Consequently, number counts derived from low and high angular resolution observations are in disagreement. This demonstrates the importance of resolution effects and the need to have realistic simulations to explore them. We built a new 2deg^2 simulation of the extragalactic sky from the far-infrared to the submillimeter. It is based on an updated version of the two star-formation mode galaxy evolution model. Using global galaxy properties, we use the abundance matching technique to populate a dark-matter lightcone and thus simulate the clustering. We produce maps and extract the sources, and show that the limited angular resolution of single-dish instruments have a strong impact on (sub)millimeter continuum observations. Taking into account these resolution effects, we are reproducing a large set of observables, including number counts, redshift distributions, and cosmic infrared background power spectra. Our simulation describes consistently the number counts from single-dish telescopes and interferometers. In particular, at 350 and 500 um, we find that number counts measured by Herschel between 5 and 50 mJy are biased towards high values by a factor 2, and that redshift distributions are biased towards low z. We also show that the clustering has an important impact on the Herschel pixel histogram used to derive number counts from P(D) analysis. Finally, we demonstrate that the large number density of red Herschel sources found in observations but not in models could be an observational artifact caused by the combination of noise, resolution effects, and the steepness of color and flux density distributions. Our simulation (SIDES) is available at this http URL
We search type 1 AGNs among emission-line galaxies, that are typically classified as type 2 AGNs based on emission line flux ratios if a broad component in the H$\alpha$ line profile is not properly investigated. Using ~24,000 type 2 AGNs at z $<$0.1 initially selected from Sloan Digital Sky Survey Data Release 7 by Bae, et al. 2014, we identify a sample of 611 type 1 AGNs based on the spectral fitting results and visual inspection. These hidden type 1 AGNs have relatively low luminosity with a mean broad \Ha\ luminosity, log L$_{\rm H\alpha }$ $=$ 40.73$\pm$0.32 \ergs\, and low Eddington ratio with a mean log L$_{bol}$/L$_{\rm Edd}$ $=$ -2.04$\pm$0.34, while they do follow the black hole mass - stellar velocity dispersion relation defined by the inactive galaxies and the reverberation-mapped type 1 AGNs. We investigate ionized gas outflows based on the [OIII] $\lambda$5007 kinematics, which show relatively high velocity dispersion and velocity shift, indicating that the line-of-sight velocity and velocity dispersion of the ionized gas in type 1 AGNs is on average larger than that of type 2 AGNs.
We aim at characterizing the large-scale distribution of H2O in G327.3-0.6, a massive star-forming region made of individual objects in different evolutionary phases. We investigate variations of H2O abundance as function of evolution. We present Herschel continuum maps at 89 and 179 $\mu$m of the whole region and an APEX map at 350 {\mu}m of the IRDC. New spectral HIFI maps toward the IRDC region covering low-energy H2O lines at 987 and 1113 GHz are also presented and combined with HIFI pointed observations of the G327 hot core. We infer the physical properties of the gas through optical depth analysis and radiative transfer modeling. The continuum emission at 89 and 179 {\mu}m follows the thermal continuum emission at longer wavelengths, with a peak at the position of the hot core, a secondary peak in the Hii region, and an arch-like layer of hot gas west of the Hii region. The same morphology is observed in the 1113 GHz line, in absorption toward all dust condensations. Optical depths of ~80 and 15 are estimated and correspond to column densities of 10^15 and 2 10^14 cm-2, for the hot core and IRDC position. These values indicate an H2O to H2 ratio of 3 10^-8 toward the hot core; the abundance of H2O does not change along the IRDC with values of some 10^-8. Infall (over ~ 20") is detected toward the hot core position with a rate of 1-1.3 10^-2 M_sun /yr, high enough to overcome the radiation pressure due to the stellar luminosity. The source structure of the hot core region is complex, with a cold outer gas envelope in expansion, situated between the outflow and the observer, extending over 0.32 pc. The outflow is seen face-on and centered away from the hot core. The distribution of H2O along the IRDC is roughly constant with an abundance peak in the more evolved object. These water abundances are in agreement with previous studies in other massive objects and chemical models.
We analyze a sample of $z$-dropout galaxies in the CANDELS GOODS South and UDS fields that have been targeted by a dedicated spectroscopic campaign aimed at detecting their Ly$\alpha$ line. Deep IRAC observations at 3.6 and 4.5 $\mu$m are used to determine the strength of optical emission lines affecting these bands at z$\sim$6.5-6.9 in order to i) investigate possible physical differences between Ly$\alpha$ emitting and non-emitting sources; ii) constrain the escape fraction of ionizing photons; iii) provide an estimate of the specific star-formation rate at high redshifts. We find evidence of strong [OIII]+H$\beta$ emission in the average (stacked) SEDs of galaxies both with and without Ly$\alpha$ emission. The blue IRAC [3.6]-[4.5] color of the stack with detected Ly$\alpha$ line can be converted into a rest-frame equivalent width EW([OIII]+H$\beta$)=1500$^{+530}_{-440}\AA$ assuming a flat intrinsic stellar continuum. This strong optical line emission enables a first estimate of f$_{esc}\lesssim$20% on the escape fraction of ionizing photons from Ly$\alpha$ detected objects. The objects with no Ly$\alpha$ line show less extreme EW([OIII]+H$\beta$)=520$^{+170}_{-150}\AA$ suggesting different physical conditions of the HII regions with respect to Ly$\alpha$-emitting ones, or a larger f$_{esc}$. The latter case is consistent with a combined evolution of f$_{esc}$ and neutral hydrogen fraction as explanation of the lack of bright Ly$\alpha$ emission at z$>$6. A lower limit on the specific star formation rate, SSFR$>$9.1$Gyr^{-1}$ for $M_{star}=2 \times 10^9 M_{\odot}$ galaxies at these redshifts can be derived from the spectroscopically confirmed sample.
The Mid-Infrared Instrument (MIRI) for the {\em James Webb Space Telescope} (JWST) will revolutionize our understanding of infrared stellar populations in the Local Volume. Using the rich {\em Spitzer}-IRS spectroscopic data-set and spectral classifications from the Surveying the Agents of Galaxy Evolution (SAGE)-Spectroscopic survey of over a thousand objects in the Magellanic Clouds, the Grid of Red supergiant and Asymptotic giant branch star ModelS ({\sc grams}), and the grid of YSO models by Robitaille et al. (2006), we calculate the expected flux-densities and colors in the MIRI broadband filters for prominent infrared stellar populations. We use these fluxes to explore the {\em JWST}/MIRI colours and magnitudes for composite stellar population studies of Local Volume galaxies. MIRI colour classification schemes are presented; these diagrams provide a powerful means of identifying young stellar objects, evolved stars and extragalactic background galaxies in Local Volume galaxies with a high degree of confidence. Finally, we examine which filter combinations are best for selecting populations of sources based on their JWST colours.
We use the Sloan Digital Sky Survey Data Release 12, which is the largest available white dwarf catalog to date, to study the evolution of the kinematical properties of the population of white dwarfs in the Galactic disc. We derive masses, ages, photometric distances and radial velocities for all white dwarfs with hydrogen-rich atmospheres. For those stars for which proper motions from the USNO-B1 catalog are available the true three-dimensional components of the stellar space velocity are obtained. This subset of the original sample comprises 20,247 objects, making it the largest sample of white dwarfs with measured three-dimensional velocities. Furthermore, the volume probed by our sample is large, allowing us to obtain relevant kinematical information. In particular, our sample extends from a Galactocentric radial distance $R_{\rm G}=7.8$~kpc to 9.3~kpc, and vertical distances from the Galactic plane ranging from $Z=-0.5$~kpc to 0.5~kpc. We examine the mean components of the stellar three-dimensional velocities, as well as their dispersions with respect to the Galactocentric and vertical distances. We confirm the existence of a mean Galactocentric radial velocity gradient, $\partial\langle V_{\rm R}\rangle/\partial R_{\rm G}=-3\pm5$~km~s$^{-1}$~kpc$^{-1}$. We also confirm North-South differences in $\langle V_{\rm z}\rangle$. Specifically, we find that white dwarfs with $Z>0$ (in the North Galactic hemisphere) have $\langle V_{\rm z}\rangle<0$, while the reverse is true for white dwarfs with $Z<0$. The age-velocity dispersion relation derived from the present sample indicates that the Galactic population of white dwarfs may have experienced an additional source of heating, which adds to the secular evolution of the Galactic disc.
We introduce and present results from the COCOA (Cluster simulatiOn Comparison with ObservAtions) code that has been developed to create idealized mock photometric observations using results from numerical simulations of star cluster evolution. COCOA is able to present the output of realistic numerical simulations of star clusters carried out using Monte Carlo or N-body codes in a way that is useful for direct comparison with photometric observations. In this paper, we describe the COCOA code and demonstrate its different applications by utilizing globular cluster (GC) models simulated with the MOCCA code. COCOA is used to synthetically observe these different GC models with optical telescopes, perform PSF photometry and subsequently produce observed colour magnitude diagrams. We also use COCOA to compare the results from synthetic observations of a cluster model that has the same age and metallicity as the Galactic GC NGC 2808 with a real observation of the same cluster carried out with a 2.2 meter optical telescope. We find that COCOA can effectively simulate realistic observations and recover photometric data. COCOA has numerous scientific applications that maybe be helpful for both theoreticians and observers that work on star clusters. Plans for further improving and developing the code are also discussed in this paper.
To shed light on the time evolution of local star formation episodes in M33, we study the association between 566 Giant Molecular Clouds (GMCs), identified through the CO (J=2-1) IRAM-all-disk survey, and 630 Young Stellar Cluster Candidates (YSCCs), selected via Spitzer-24~$\mu$m emission. The spatial correlation between YSCCs and GMCs is extremely strong, with a typical separation of 17~pc, less than half the CO(2--1) beamsize, illustrating the remarkable physical link between the two populations. GMCs and YSCCs follow the HI filaments, except in the outermost regions where the survey finds fewer GMCs than YSCCs, likely due to undetected, low CO-luminosity clouds. The GMCs have masses between 2$\times 10^4$ and 2$\times 10^6$ M$_\odot$ and are classified according to different cloud evolutionary stages: inactive clouds are 32$\%$ of the total, classified clouds with embedded and exposed star formation are 16$\%$ and 52$\%$ of the total respectively. Across the regular southern spiral arm, inactive clouds are preferentially located in the inner part of the arm, possibly suggesting a triggering of star formation as the cloud crosses the arm. Some YSCCs are embedded star-forming sites while the majority have GALEX-UV and H$\alpha$ counterparts with estimated cluster masses and ages. The distribution of the non-embedded YSCC ages peaks around 5~Myrs with only a few being as old as 8--10~Myrs. These age estimates together with the number of GMCs in the various evolutionary stages lead us to conclude that 14~Myrs is a typical lifetime of a GMC in M33, prior to cloud dispersal. The inactive and embedded phases are short, lasting about 4 and 2~Myrs respectively. This underlines that embedded YSCCs rapidly break out from the clouds and become partially visible in H$\alpha$ or UV long before cloud dispersal.
The $\sigma$ Orionis cluster is important for studying protoplanetary disk evolution, as its intermediate age ($\sim$3-5 Myr) is comparable to the median disk lifetime. We use ALMA to conduct a high-sensitivity survey of dust and gas in 92 protoplanetary disks around $\sigma$ Orionis members with $M_{\ast}\gtrsim0.1 M_{\odot}$. Our observations cover the 1.33 mm continuum and several CO $J=2-1$ lines: out of 92 sources, we detect 37 in the mm continuum and six in $^{12}$CO, three in $^{13}$CO, and none in C$^{18}$O. Using the continuum emission to estimate dust mass, we find only 11 disks with $M_{\rm dust}\gtrsim10 M_{\oplus}$, indicating that after only a few Myr of evolution most disks lack sufficient dust to form giant planet cores. Stacking the individually undetected continuum sources limits their average dust mass to 5$\times$ lower than that of the faintest detected disk, supporting theoretical models that indicate rapid dissipation once disk clearing begins. Comparing the protoplanetary disk population in $\sigma$ Orionis to those of other star-forming regions supports the steady decline in average dust mass and the steepening of the $M_{\rm dust}$-$M_{\ast}$ relation with age; studying these evolutionary trends can inform the relative importance of different disk processes during key eras of planet formation. External photoevaporation from the central O9 star is influencing disk evolution throughout the region: dust masses clearly decline with decreasing separation from the photoionizing source, and the handful of CO detections exist at projected separations $>1.5$ pc. Collectively, our findings indicate that giant planet formation is inherently rare and/or well underway by a few Myr of age.
The latest measurements of CMB electron scattering optical depth reported by Planck significantly reduces the allowed space of HI reionization models, pointing towards a later ending and/or less extended phase transition than previously believed. Reionization impulsively heats the intergalactic medium (IGM) to $\sim10^4$ K, and owing to long cooling and dynamical times in the diffuse gas, comparable to the Hubble time, memory of reionization heating is retained. Therefore, a late ending reionization has significant implications for the structure of the $z\sim5-6$ Lyman-$\alpha$ (ly$\alpha$) forest. Using state-of-the-art hydrodynamical simulations that allow us to vary the timing of reionization and its associated heat injection, we argue that extant thermal signatures from reionization can be detected via the ly$\alpha$ forest power spectrum at $5< z<6$. This arises because the small-scale cutoff in the power depends not only the the IGMs temperature at these epochs, but is also particularly sensitive to the pressure smoothing scale set by the IGMs full thermal history. Comparing our different reionization models with existing measurements of the ly$\alpha$ forest flux power spectrum at $z=5.0-5.4$, we find that models satisfying Planck's $\tau_e$ constraint, favor a moderate amount of heat injection consistent with galaxies driving reionization, but disfavoring quasar driven scenarios. We explore the impact of different reionization histories and heating models on the shape of the power spectrum, and find that they can produce similar effects, but argue that this degeneracy can be broken with high enough quality data. We study the feasibility of measuring the flux power spectrum at $z\simeq 6$ using mock quasar spectra and conclude that a sample of $\sim10$ high-resolution spectra with attainable S/N ratio will allow to discriminate between different reionization scenarios.
We study the intra-cluster magnetic field in the poor galaxy cluster Abell 194 by complementing radio data, at different frequencies, with data in the optical and X-ray bands. We analyze new total intensity and polarization observations of Abell 194 obtained with the Sardinia Radio Telescope (SRT). We use the SRT data in combination with archival Very Large Array observations to derive both the spectral aging and Rotation Measure (RM) images of the radio galaxies 3C40A and 3C40B embedded in Abell 194. The optical analysis indicates that Abell 194 does not show a major and recent cluster merger, but rather agrees with a scenario of accretion of small groups. Under the minimum energy assumption, the lifetimes of synchrotron electrons in 3C40B measured from the spectral break are found to be 157 Myrs. The break frequency image and the electron density profile inferred from the X-ray emission are used in combination with the RM data to constrain the intra-cluster magnetic field power spectrum. By assuming a Kolmogorov power law power spectrum, we find that the RM data in Abell 194 are well described by a magnetic field with a maximum scale of fluctuations of Lambda_max=64 kpc and a central magnetic field strength of <B0>=1.5 microG. Further out, the field decreases with the radius following the gas density to the power of eta=1.1. Comparing Abell 194 with a small sample of galaxy clusters, there is a hint of a trend between central electron densities and magnetic field strengths.
We analyze theoretically the Schrodinger-Poisson equation in two transverse dimensions in the presence of a Kerr term. The model describes the nonlinear propagation of optical beams in thermooptical media and can be regarded as an analogue system for a self-gravitating self-interacting wave. We compute numerically the family of radially symmetric ground state bright stationary solutions for focusing and defocusing local nonlinearity, keeping in both cases a focusing nonlocal nonlinearity. We also analyze excited states and oscillations induced by fixing the temperature at the borders of the material. We provide simulations of soliton interactions, drawing analogies with the dynamics of galactic cores in the scalar field dark matter scenario.
We investigate the formation of circumstellar disks and outflows subsequent to the collapse of molecular cloud cores with the magnetic field and turbulence. Numerical simulations are performed by using an adaptive mesh refinement to follow the evolution up to $\sim 1000$~yr after the formation of a protostar. In the simulations, circumstellar disks are formed around the protostars; those in magnetized models are considerably smaller than those in nonmagnetized models, but their size increases with time. The models with stronger magnetic field tends to produce smaller disks. During evolution in the magnetized models, the mass ratios of a disk to a protostar is approximately constant at $\sim 1-10$\%. The circumstellar disks are aligned according to their angular momentum, and the outflows accelerate along the magnetic field on the $10-100$~au scale; this produces a disk that is misaligned with the outflow. The outflows are classified into two types: a magneto-centrifugal wind and a spiral flow. In the latter, because of the geometry, the axis of rotation is misaligned with the magnetic field. The magnetic field has an internal structure in the cloud cores, which also causes misalignment between the outflows and the magnetic field on the scale of the cloud core. The distribution of the angular momentum vectors in a core also has a non-monotonic internal structure. This should create a time-dependent accretion of angular momenta onto the circumstellar disk. Therefore, the circumstellar disks are expected to change their orientation as well as their sizes in the long-term evolutions.
Links to: arXiv, form interface, find, astro-ph, recent, 1703, contact, help (Access key information)
Recent studies suggest that the quenching properties of galaxies are correlated over several mega-parsecs. The large-scale "galactic conformity" phenomenon around central galaxies has been regarded as a potential signature of "galaxy assembly bias" or "pre-heating", both of which interpret conformity as a result of direct environmental effects acting on galaxy formation. Building on the iHOD halo quenching framework developed in Zu & Mandelbaum (2015, 2016), we discover that our fiducial halo mass quenching model, without any galaxy assembly bias, can successfully explain the overall environmental dependence and the conformity of galaxy colours in SDSS, as measured by the mark correlation functions of galaxy colours and the red galaxy fractions around isolated primaries, respectively. Our fiducial iHOD halo quenching mock also correctly predicts the differences in the spatial clustering and galaxy-galaxy lensing signals between the more vs. less red galaxy subsamples, split by the red-sequence ridge-line at fixed stellar mass. Meanwhile, models that tie galaxy colours fully or partially to halo assembly bias have difficulties in matching all these observables simultaneously. Therefore, we demonstrate that the observed environmental dependence of galaxy colours can be naturally explained by the combination of 1) halo quenching and 2) the variation of halo mass function with environment --- an indirect environmental effect mediated by two separate physical processes.
We present the detection of supermassive black holes (BHs) in two Virgo ultracompact dwarf galaxies (UCDs), VUCD3 and M59cO. We use adaptive optics assisted data from the Gemini/NIFS instrument to derive radial velocity dispersion profiles for both objects. Mass models for the two UCDs are created using multi-band Hubble Space Telescope (HST) imaging, including the modeling of mild color gradients seen in both objects. We then find a best-fit stellar mass-to-light ratio ($M/L$) and BH mass by combining the kinematic data and the deprojected stellar mass profile using Jeans Anisotropic Models (JAM). Assuming axisymmetric isotropic Jeans models, we detect BHs in both objects with masses of $4.4^{+2.5}_{-3.0} \times 10^6$ $M_{\odot}$ in VUCD3 and $5.8^{+2.5}_{-2.8} \times 10^6$ $M_{\odot}$ in M59cO (3$\sigma$ uncertainties). The BH mass is degenerate with the anisotropy parameter, $\beta_z$; for the data to be consistent with no BH requires $\beta_z = 0.4$ and $\beta_z = 0.6$ for VUCD3 and M59cO, respectively. Comparing these values with nuclear star clusters shows that while it is possible that these UCDs are highly radially anisotropic, it seems unlikely. These detections constitute the second and third UCDs known to host supermassive BHs. They both have a high fraction of their total mass in their BH; $\sim$13% for VUCD3 and $\sim$18% for M59cO. They also have low best-fit stellar $M/L$s, supporting the proposed scenario that most massive UCDs host high mass fraction BHs. The properties of the BHs and UCDs are consistent with both objects being the tidally stripped remnants of $\sim$10$^9$ M$_\odot$ galaxies.
Jets from supermassive black holes in the centres of galaxy clusters are a potential candidate for moderating gas cooling and subsequent star formation through depositing energy in the intra-cluster gas. In this work, we simulate the jet-intra-cluster medium interaction using the moving-mesh magnetohydrodynamics code Arepo. Our model injects supersonic, low density, collimated and magnetised outflows in cluster centres, which are then stopped by the surrounding gas, thermalise and inflate low-density cavities filled with cosmic-rays. We perform high-resolution, non-radiative simulations of the lobe creation, expansion and disruption, and find that its dynamical evolution is in qualitative agreement with simulations of idealised low-density cavities that are dominated by a large-scale Rayleigh-Taylor instability. The buoyant rising of the lobe does not create energetically significant small-scale chaotic motion in a volume-filling fashion, but rather a systematic upward motion in the wake of the lobe and a corresponding back-flow perpendicular to it. We find that, overall, 50 per cent of the injected energy ends up in material which is not part of the lobe, and about 25 per cent remains in the inner 100 kpc. We conclude that jet-inflated, buoyantly rising cavities drive systematic gas motions which play an important role in heating the central regions, while mixing of lobe material is sub-dominant. Encouragingly, the main mechanisms responsible for this energy deposition can be modelled already at resolutions within reach in future, high-resolution cosmological simulations of galaxy clusters.
Stars form in cold molecular clouds. However, molecular gas is difficult to observe because the most abundant molecule (H2) lacks a permanent dipole moment. Rotational transitions of CO are often used as a tracer of H2, but CO is much less abundant and the conversion from CO intensity to H2 mass is often highly uncertain. Here we present a new method for estimating the column density of cold molecular gas (Sigma_gas) using optical spectroscopy. We utilise the spatially resolved H-alpha maps of flux and velocity dispersion from the Sydney-AAO Multi-object Integral-field spectrograph (SAMI) Galaxy Survey. We derive maps of Sigma_gas by inverting the multi-freefall star formation relation, which connects the star formation rate surface density (Sigma_SFR) with Sigma_gas and the turbulent Mach number (Mach). Based on the measured range of Sigma_SFR = 0.005-1.5 M_sol/yr/kpc^2 and Mach = 18-130, we predict Sigma_gas = 7-200 M_sol/pc^2 in the star-forming regions of our sample of 260 SAMI galaxies. These values are close to previously measured Sigma_gas obtained directly with unresolved CO observations of similar galaxies at low redshift. We classify each galaxy in our sample as 'Star-forming' (219) or 'Composite/AGN/Shock' (41), and find that in Composite/AGN/Shock galaxies the average Sigma_SFR, Mach, and Sigma_gas are enhanced by factors of 2.0, 1.6, and 1.3, respectively, compared to Star-forming galaxies. We compare our predictions of Sigma_gas with those obtained by inverting the Kennicutt-Schmidt relation and find that our new method is a factor of two more accurate in predicting Sigma_gas, with an average deviation of 32% from the actual Sigma_gas.
We combine Gaia data release 1 astrometry with Sloan Digital Sky Survey (SDSS) images taken some ~10-15 years earlier, to measure proper motions of stars in the halo of our Galaxy. The SDSS-Gaia proper motions have typical statistical errors of 2 mas/yr down to r~20 mag, and are robust to variations with magnitude and colour. Armed with this exquisite set of halo proper motions, we identify RR Lyrae, blue horizontal branch (BHB), and K giant stars in the halo, and measure their net rotation with respect to the Galactic disc. We find evidence for a gently rotating prograde signal (< V_\phi > ~5-25 km/s) in the halo stars, which shows little variation with Galactocentric radius out to 50 kpc. There is also tentative evidence for a kinematic bias with metallicity, whereby the metal richer BHB and K giant stars have slightly stronger prograde rotation than the metal poorer stars. Using the Auriga simulation suite we find that the old (T >10 Gyr) stars in the simulated halos exhibit mild prograde rotation, with little dependence on radius or metallicity, in general agreement with the observations. The weak halo rotation suggests that the Milky Way has a minor in situ halo component, and has undergone a relatively quiet accretion history.
We use $>$9400 $\log(m/M_{\odot})>10$ quiescent and star-forming galaxies at $z\lesssim2$ in COSMOS/UltraVISTA to study the average size evolution of these systems, with focus on the rare, ultra-massive population at $\log(m/M_{\odot})>11.4$. The large 2-square degree survey area delivers a sample of $\sim400$ such ultra-massive systems. Accurate sizes are derived using a calibration based on high-resolution images from the Hubble Space Telescope. We find that, at these very high masses, the size evolution of star-forming and quiescent galaxies is almost indistinguishable in terms of normalization and power-law slope. We use this result to investigate possible pathways of quenching massive $m>M^*$ galaxies at $z<2$. We consistently model the size evolution of quiescent galaxies from the star-forming population by assuming different simple models for the suppression of star-formation. These models include an instantaneous and delayed quenching without altering the structure of galaxies and a central starburst followed by compaction. We find that instantaneous quenching reproduces well the observed mass-size relation of massive galaxies at $z>1$. Our starburst$+$compaction model followed by individual growth of the galaxies by minor mergers is preferred over other models without structural change for $\log(m/M_{\odot})>11.0$ galaxies at $z>0.5$. None of our models is able to meet the observations at $m>M^*$ and $z<1$ with out significant contribution of post-quenching growth of individual galaxies via mergers. We conclude that quenching is a fast process in galaxies with $ m \ge 10^{11} M_\odot$, and that major mergers likely play a major role in the final steps of their evolution.
Many features near the Galactic Center have been called 3-kiloparsec arms. We
reached a point of having too many divergent data, making it difficult to be
constrained by a single physical model. Their differing characteristics suggest
different physical and dynamical objects.
Radial velocity data on the so-called 3-kpc arms do not coincide with radial
velocities of major spiral arms near 3kpc, nor near 2 kpc, nor near 4 kpc from
the Galactic Center (Fig. 1 and 2).
Different 3-kpc arm features may require different models: turbulence around
a shock in a Galactic density wave between 2 and 4 kpc from the Galactic Center
(Table 1), or nuclear rotation between 0 and 2 kpc from the Galactic Center
region (Table 2), or a putative radial expansion between 0 and 4 kpc from the
Galactic Center.
Despite their naming as Near 3-kiloparsec arms or Far 3-kiloparsec arms,
these features are not major arms. Those 3-kpc arm features nearer the Galactic
Center (within 13o of Galactic longitude) may be different than those farther
out (Table 2).
Here we show that the plethora of observed '3-kpc arm' features can be
separated in two: those with Galactic longitude of 13 degrees or more away from
the Galactic Center (Table 1 - some of which are possibly associated with the
observed major spiral arms), and those within 13 degrees from the Galactic
Center (Table 2 - some of which are possibly associated with the observed
central bars; Fig.1 and Fig.2).
We use new deep 21 cm HI observations of the moderately inclined galaxy NGC 4559 in the HALOGAS survey to investigate the properties of extra-planar gas. We use TiRiFiC to construct simulated data cubes to match the HI observations. We find that a thick disk component of scale height $\sim\,2\,\mathrm{kpc}$, characterized by a negative vertical gradient in its rotation velocity (lag) of $\sim13 \pm 5$ km s$^{-1}$ kpc$^{-1}$ is an adequate fit to extra-planar gas features. The tilted ring models also present evidence for a decrease in the magnitude of the lag outside of $R_{25}$, and a radial inflow of $\sim 10$ km s$^{-1}$. We extracted lagging extra-planar gas through Gaussian velocity profile fitting. From both the 3D models and and extraction analyses we conclude that $\sim10-20\%$ of the total {\HI} mass is extra-planar. Most of the extra-planar gas is spatially coincident with regions of star formation in spiral arms, as traced by H$\alpha$ and GALEX FUV images, so it is likely due to star formation processes driving a galactic fountain. We also find the signature of a filament of a kinematically "forbidden" HI, containing $\sim 1.4\times 10^{6}$ M$_{\odot}$ of HI, and discuss its potential relationship to a nearby HI hole. We discover a previously undetected dwarf galaxy in HI located $\sim 0.4^{\circ}$ ($\sim 58$ kpc) from the center of NGC 4559, containing $\sim 4\times10^{5}$ M$_{\odot}$. This dwarf has counterpart sources in SDSS with spectra typical of HII regions, and we conclude it is two merging blue compact dwarf galaxies.
We study how the gas in a sample of galaxies (M* > 10e9 Msun) in clusters, obtained in a cosmological simulation, is affected by the interaction with the intra-cluster medium (ICM). The dynamical state of each elemental parcel of gas is studied using the total energy. At z ~ 2, the galaxies in the simulation are evenly distributed within clusters, moving later on towards more central locations. In this process, gas from the ICM is accreted and mixed with the gas in the galactic halo. Simultaneously, the interaction with the environment removes part of the gas. A characteristic stellar mass around M* ~ 10e10 Msun appears as a threshold marking two differentiated behaviours. Below this mass, galaxies are located at the external part of clusters and have eccentric orbits. The effect of the interaction with the environment is marginal. Above, galaxies are mainly located at the inner part of clusters with mostly radial orbits with low velocities. In these massive systems, part of the gas, strongly correlated with the stellar mass of the galaxy, is removed. The amount of removed gas is sub-dominant compared with the quantity of retained gas which is continuously influenced by the hot gas coming from the ICM. The analysis of individual galaxies reveals the existence of a complex pattern of flows, turbulence and a constant fuelling of gas to the hot corona from the ICM that could make the global effect of the interaction of galaxies with their environment to be substantially less dramatic than previously expected.
The SgrB2 molecular cloud contains several sites forming high-mass stars. SgrB2(N) is one of its main centers of activity. It hosts several compact and UCHII regions, as well as two known hot molecular cores (SgrB2(N1) and SgrB2(N2)), where complex organic molecules are detected. Our goal is to use the high sensitivity of ALMA to characterize the hot core population in SgrB2(N) and shed a new light on the star formation process. We use a complete 3 mm spectral line survey conducted with ALMA to search for faint hot cores in SgrB2(N). We report the discovery of three new hot cores that we call SgrB2(N3), SgrB2(N4), and SgrB2(N5). The three sources are associated with class II methanol masers, well known tracers of high-mass star formation, and SgrB2(N5) also with a UCHII region. The chemical composition of the sources and the column densities are derived by modelling the whole spectra under the assumption of LTE. The H2 column densities are computed from ALMA and SMA continuum emission maps. The H2 column densities of these new hot cores are found to be 16 up to 36 times lower than the one of the main hot core Sgr B2(N1). Their spectra have spectral line densities of 11 up to 31 emission lines per GHz, assigned to 22-25 molecules. We derive rotational temperatures around 140-180 K for the three new hot cores and mean source sizes of 0.4 for SgrB2(N3) and 1.0 for SgrB2(N4) and SgrB2(N5). SgrB2(N3) and SgrB2(N5) show high velocity wing emission in typical outflow tracers, with a bipolar morphology in their integrated intensity maps suggesting the presence of an outflow, like in SgrB2(N1). The associations of the hot cores with class II methanol masers, outflows, and/or UCHII regions tentatively suggest the following age sequence: SgrB2(N4), SgrB2(N3), SgrB2(N5), SgrB2(N1). The status of SgrB2(N2) is unclear. It may contain two distinct sources, a UCHII region and a very young hot core.
The galactocentric distance of quasar absorption outflows are conventionally determined using absorption troughs from excited states, a method hindered by severely saturated or self-blended absorption troughs. We propose a novel method to estimate the size of a broad absorption line (BAL) region which partly obscures an emission line region by assuming virialized gas in the emission region surrounding a supermassive black hole with known mass. When a spiky Lya1216 line emission is present at the flat bottom of the deep Nv1240 absorption trough, the size of BAL region can be estimated. We have found 3 BAL quasars in the SDSS database showing such Lya lines. The scale of their BAL outflows are found to be 3-26 pc, moderately larger than the theoretical scale (0.01-0.1 pc) of trough forming regions for winds originating from accretion discs, but significantly smaller than most outflow sizes derived using the absorption troughs of the excited states of ions. For these three outflows, the lower limits of ratio of kinetic luminosity to Eddington luminosity are 0.02%-0.07%. These lower limits are substantially smaller than that is required to have significant feedback effect on their host galaxies.
We use new X-ray data obtained with the Nuclear Spectroscopic Telescope Array (NuSTAR), near-infrared (NIR) fluxes, and mid-infrared (MIR) spectra of a sample of 24 unobscured type 1 active galactic nuclei (AGN) to study the correlation between various hard X-ray bands between 3 and 80 keV and the infrared (IR) emission. The IR to X-ray correlation spectrum (IRXCS) shows a maximum at ~15-20 micron, coincident with the peak of the AGN contribution to the MIR spectra of the majority of the sample. There is also a NIR correlation peak at ~2 micron, which we associate with the NIR bump observed in some type 1 AGN at ~1-5 micron and is likely produced by nuclear hot dust emission. The IRXCS shows practically the same behaviour in all the X-ray bands considered, indicating a common origin for all of them. We finally evaluated correlations between the X-ray luminosities and various MIR emission lines. All the lines show a good correlation with the hard X-rays (rho>0.7), but we do not find the expected correlation between their ionization potentials and the strength of the IRXCS.
We present new Atacama Large Millimeter/sub-millimeter Array (ALMA) 1.3 mm continuum observations of the SR 24S transition disk with an angular resolution $\lesssim0.18"$ (12 au radius). We perform a multi-wavelength investigation by combining new data with previous ALMA data at 0.45 mm. The visibilities and images of the continuum emission at the two wavelengths are well characterized by a ring-like emission. Visibility modeling finds that the ring-like emission is narrower at longer wavelengths, in good agreement with models of dust trapping in pressure bumps, although there are complex residuals that suggest potentially asymmetric structures. The 0.45 mm emission has a shallower profile inside the central cavity than the 1.3 mm emission. In addition, we find that the $^{13}$CO and C$^{18}$O (J=2-1) emission peaks at the center of the continuum cavity. We do not detect either continuum or gas emission from the northern companion to this system (SR 24N), which is itself a binary system. The upper limit for the dust disk mass of SR 24N is $\lesssim 0.12\,M_{\bigoplus}$, which gives a disk mass ratio in dust between the two components of $M_{\mathrm{dust, SR\,24S}}/M_{\mathrm{dust, SR\,24N}}\gtrsim840$. The current ALMA observations may imply that either planets have already formed in the SR 24N disk or that dust growth to mm-sizes is inhibited there and that only warm gas, as seen by ro-vibrational CO emission inside the truncation radii of the binary, is present.
We present IFS-RedEx, a spectrum and redshift extraction pipeline for integral-field spectrographs. A key feature of the tool is a wavelet-based spectrum cleaner. It identifies reliable spectral features, reconstructs their shapes, and suppresses the spectrum noise. This gives the technique an advantage over conventional methods like Gaussian filtering, which only smears out the signal. As a result, the wavelet-based cleaning allows the quick identification of true spectral features. We test the cleaning technique with degraded MUSE spectra and find that it can detect spectrum peaks down to S/N = 8 while reporting no fake detections. We apply IFS-RedEx to MUSE data of the strong lensing cluster MACSJ1931.8-2635 and extract 54 spectroscopic redshifts. We identify 29 cluster members and 22 background galaxies with z >= 0.4. IFS-RedEx is open source and publicly available.
We identify and investigate known ultracool stars and brown dwarfs that are being observed or indirectly constrained by the Gaia mission. These objects will be the core of the Gaia ultracool dwarf sample composed of all dwarfs later than M7 that Gaia will provide direct or indirect information on. We match known L and T dwarfs to the Gaia first data release, the Two Micron All Sky Survey and the Wide-field Infrared Survey Explorer AllWISE survey and examine the Gaia and infrared colours, along with proper motions, to improve spectral typing, identify outliers and find mismatches. There are 321 L and T dwarfs observed directly in the Gaia first data release, of which 10 are later than L7. This represents 45 % of all the known LT dwarfs with estimated Gaia G magnitudes brighter than 20.3 mag. We determine proper motions for the 321 objects from Gaia and the Two Micron All Sky Survey positions. Combining the Gaia and infrared magnitudes provides useful diagnostic diagrams for the determination of L and T dwarf physical parameters. We then search the Tycho-Gaia astrometric solution Gaia first data release subset to find any objects with common proper motions to known L and T dwarfs and a high probability of being related. We find 15 new candidate common proper motion systems.
The Carina nebula hosts a large number of globulettes. The majority are of planetary mass, but there are also those with masses of several tens up to a few hundred Jupiter masses. We carried out radio observations of molecular line emission in 12CO and 13CO (2-1) and (3-2) of 12 larger objects in addition of positions in adjacent shell structures using APEX. All selected objects were detected with radial velocities shifted relative to the emission from related shell structures and background molecular clouds. Globulettes along the western part of an extended dust shell show a small spread in velocity with small velocity shifts relative to the shell. This system of globulettes and shell structures in the foreground of the bright nebulosity surrounding the cluster Trumpler 14 is expanding with a few km/s relative to the cluster. The Carina globulettes are compact and denser than objects of similar mass studied previously in the Rosette nebula. Some globulettes in Tr 14 are located far from any shell structures. These objects move at a similar speed as the globulettes along the shell. The distribution and velocities of the globulettes studied suggest that they have originated from eroding shells and elephant trunks.
A Type IIn supernova (SN) is dominated by the interaction of SN ejecta with the circumstellar medium (CSM). Some SNe IIn (e.g., SN 2006jd) have episodes of re-brightening ("bumps") in their light curves. We present iPTF13z, a SN IIn discovered by the intermediate Palomar Transient Factory (iPTF) and characterised by several bumps in its light curve. We analyse this peculiar behaviour trying to infer the properties of the CSM and of the SN explosion, as well as the nature of its progenitor star. We obtained multi-band optical photometry for over 1000 days after discovery with the P48 and P60 telescopes at Palomar Observatory. We obtained low-resolution optical spectra in the same period. We did an archival search for progenitor outbursts. We analyse our photometry and spectra, and compare iPTF13z to other SNe IIn. A simple analytical model is used to estimate properties of the CSM. iPTF13z was a SN IIn showing a light curve with five bumps during its decline phase. The bumps had amplitudes between 0.4 and 0.9 mag and durations between 20 and 120 days. The most prominent bumps appeared in all our different optical bands. The spectra showed typical SN IIn characteristics, with emission lines of H$\alpha$ (with broad component FWHM ~$10^{3}-10^{4} ~{\rm ~km ~s^{-1}}$ and narrow component FWHM ~$10^2 \rm ~km ~s^{-1}$) and He I, but also with Fe II, Ca II, Na I D and H$\beta$ P-Cygni profiles (with velocities of ~$10^{3}$ ${\rm ~km ~s^{-1}}$). A pre-explosion outburst was identified lasting $\gtrsim 50$ days, with $M_r \approx -15$ mag around 210 days before discovery. Large, variable progenitor mass-loss rates (~> 0.01 $M_{\odot} \rm ~yr^{-1}$) and CSM densities (~> 10$^{-16}$ g cm$^{-3}$) are derived. We suggest that the light curve bumps of iPTF13z arose from SN ejecta interacting with denser regions in the CSM, possibly produced by the eruptions of a luminous blue variable star.
Links to: arXiv, form interface, find, astro-ph, recent, 1703, contact, help (Access key information)
We determine the physical properties and turbulence driving mode of molecular clouds formed in numerical simulations of a Milky Way-type disc galaxy with parsec-scale resolution. The clouds form through gravitational fragmentation of the gas, leading to average values for mass, radii and velocity dispersion in good agreement with observations of Milky Way clouds. The driving parameter (b) for the turbulence within each cloud is characterised by the ratio of the density contrast (sigma_rho) to the average Mach number (Mach) within the cloud, b = sigma_rho/Mach. As shown in previous works, b ~ 1/3 indicates solenoidal (divergence-free) driving and b ~ 1 indicates compressive (curl-free) driving. We find that the average b value of all the clouds formed in the simulations has a lower limit of b > 0.2. Importantly, we find that b has a broad distribution, covering values from purely solenoidal to purely compressive driving. Tracking the evolution of individual clouds reveals that the b value for each cloud does not vary significantly over their lifetime. Finally, we perform a resolution study with minimum cell sizes of 8, 4, 2 and 1 pc and find that the average b value increases with increasing resolution. Therefore, we conclude that our measured b values are strictly lower limits and that a resolution better than 1 pc is required for convergence. However, regardless of the resolution, we find that b varies by factors of a few in all cases, which means that the effective driving mode alters significantly from cloud to cloud.
We present the VLA-COSMOS 3 GHz Large Project based on 384 hours of observations with the Karl G. Jansky Very Large Array (VLA) at 3 GHz (10 cm) toward the two square degree Cosmic Evolution Survey (COSMOS) field. The final mosaic reaches a median rms of 2.3 uJy/beam over the two square degrees at an angular resolution of 0.75". To fully account for the spectral shape and resolution variations across the broad (2 GHz) band, we image all data with a multiscale, multifrequency synthesis algorithm. We present a catalog of 10,830 radio sources down to 5 sigma, out of which 67 are combined from multiple components. Comparing the positions of our 3 GHz sources with those from the Very Long Baseline Array (VLBA)-COSMOS survey, we estimate that the astrometry is accurate to 0.01" at the bright end (signal-to-noise ratio, S/N_3GHz > 20). Survival analysis on our data combined with the VLA-COSMOS 1.4~GHz Joint Project catalog yields an expected median radio spectral index of alpha=-0.7. We compute completeness corrections via Monte Carlo simulations to derive the corrected 3 GHz source counts. Our counts are in agreement with previously derived 3 GHz counts based on single-pointing (0.087 square degrees) VLA data. In summary, the VLA-COSMOS 3 GHz Large Project simultaneously provides the largest and deepest radio continuum survey at high (0.75") angular resolution to date, bridging the gap between last-generation and next-generation surveys.
(abridged) We study the composition of the faint radio population selected from the VLA-COSMOS 3GHz Large Project. The survey covers a 2.6sq.deg. area with a mean rms of ~2.3uJy/b, cataloging 10830 sources (>5sigma). Combining these radio data with optical, near-infrared (UltraVISTA), mid-infrared (Spitzer/IRAC) data, and X-ray data (Chandra), we find counterparts to radio sources for ~93% of the radio sample (in the areas of the COSMOS field not affected by saturated or bright sources in the optical to NIR bands), reaching out to z<6. We further classify the sources as star forming galaxies or AGN based on various criteria, such as X-ray luminosity, observed MIR color, UV-FIR spectral-energy distribution, rest-frame NUV-optical color corrected for dust extinction, and radio-excess relative to that expected from the hosts' star-formation rate. We separate the AGN into sub-samples dominated by low-to-moderate and moderate-to-high radiative luminosity AGN, candidates for high-redshift analogues to local low- and high-excitation emission line AGN, respectively. We study the fractional contributions of these sub-populations down to radio flux levels of ~11uJy at 3GHz (or ~20uJy at 1.4GHz assuming a spectral index of -0.7). We find that the dominant fraction at 1.4GHz flux densities above ~200uJy is constituted of low-to-moderate radiative luminosity AGN. Below densities of ~100uJy the fraction of star-forming galaxies increases to ~60%, followed by the moderate-to-high radiative luminosity AGN (~20%), and low-to-moderate radiative luminosity AGN (~20%). Based on this observational evidence, we extrapolate the fractions down to sensitivities of the SKA. Our estimates suggest that at the faint flux limits to be reached by the SKA1 surveys, a selection based only on radio flux limits can provide a simple tool to efficiently identify samples highly (>75%) dominated by star-forming galaxies.
We explore the multiwavelength properties of AGN host galaxies for different classes of radio-selected AGN out to z$\lesssim$6 via a multiwavelength analysis of about 7700 radio sources in the COSMOS field. The sources were selected with the Very Large Array (VLA) at 3 GHz (10 cm) within the VLA-COSMOS 3 GHz Large Project, and cross-matched with multiwavelength ancillary data. This is the largest sample of high-redshift (z$\lesssim$6) radio sources with exquisite photometric coverage and redshift measurements available. We constructed a sample of moderate-to-high radiative luminosity AGN (HLAGN) via spectral energy distribution (SED) decomposition combined with standard X-ray and mid-infrared diagnostics. Within the remainder of the sample we further identified low-to-moderate radiative luminosity AGN (MLAGN) via excess in radio emission relative to the star formation rates in their host galaxies. We show that AGN power in HLAGN occurs predominantly in radiative form, while MLAGN display a substantial mechanical AGN luminosity component. We found significant differences in the host properties of the two AGN classes, as a function of redshift. At z$<$1.5, MLAGN appear to reside in significantly more massive and less star-forming galaxies compared to HLAGN. At z$>$1.5, we observed a reversal in the behaviour of the stellar mass distributions with the HLAGN populating the higher stellar mass tail. We interpret this finding as a possible hint of the downsizing of galaxies hosting HLAGN, with the most massive galaxies triggering AGN activity earlier than less massive galaxies, and then fading to MLAGN at lower redshifts. Our conclusion is that HLAGN and MLAGN samples trace two distinct galaxy and AGN populations in a wide range of redshifts, possibly resembling the radio AGN types often referred to as radiative- and jet-mode (or high- and low-excitation), respectively.
We examine the behaviour of the infrared-radio correlation (IRRC) over the range $0<z<6$ using new, highly sensitive 3GHz observations with the Karl G. Jansky Very Large Array (VLA) and infrared data from the Herschel Space Observatory in the 2deg$^{2}$ COSMOS field. We distinguish between objects where emission is believed to arise solely from star-formation, and those where an active galactic nucleus (AGN) is thought to be present. We account for non-detections in the radio or in the infrared using a doubly-censored survival analysis. We find that the IRRC of star-forming galaxies, quantified by the infrared-to-1.4GHz radio luminosity ratio ($q_{\rm TIR}$), decreases with increasing redshift: $q_{\rm TIR}(z)=(2.88\pm0.03)(1+z)^{-0.19\pm0.01}$. Moderate-to-high radiative luminosity AGN do not follow the same $q_{\rm TIR}$$(z)$ trend, having a lower normalisation and steeper decrease with redshift. We cannot rule out the possibility that unidentified AGN contributions only to the radio regime may be steepening the observed $q_{\rm TIR}(z)$ trend of the star-forming population. An increasing fractional contribution to the observed 3GHz flux by free-free emission of star-forming galaxies may also affect the derived evolution. However, we find that the standard (M82-based) assumption of the typical radio spectral energy distribution (SED) for star-forming galaxies is inconsistent with our results. This suggests a more complex shape of the typical radio SED for star-forming galaxies, and that imperfect $K$ corrections in the radio may govern the derived redshift trend of $q_{\rm TIR}$. Lastly, we present a redshift-dependent relation between rest-frame 1.4GHz radio luminosity and star formation rate taking the derived redshift trend into account.
We make use of the deep Karl G. Jansky Very Large Array (VLA) COSMOS radio observations at 3 GHz to infer radio luminosity functions of star-forming galaxies up to redshifts of z~5 based on approximately 6000 detections with reliable optical counterparts. This is currently the largest radio-selected sample available out to z~5 across an area of 2 square degrees with a sensitivity of rms=2.3 ujy/beam. By fixing the faint and bright end shape of the radio luminosity function to the local values, we find a strong redshift trend that can be fitted with a pure luminosity evolution L~(1+z)^{(3.16 +- 0.2)-(0.32 +- 0.07) z}. We estimate star formation rates (SFRs) from our radio luminosities using an infrared (IR)-radio correlation that is redshift dependent. By integrating the parametric fits of the evolved luminosity function we calculate the cosmic SFR density (SFRD) history since z~5. Our data suggest that the SFRD history peaks between 2<z<3 and that the ultraluminous infrared galaxies (ULIRGs; 100 Msol/yr<SFR<1000 Msol/yr) contribute up to ~25% to the total SFRD in the same redshift range. Hyperluminous infrared galaxies (HyLIRGs; SFR>1000 Msol/yr) contribute an additional <2% in the entire observed redshift range. We find evidence of a potential underestimation of SFRD based on ultraviolet (UV) rest-frame observations of Lyman break galaxies (LBGs) at high redshifts (z>4) on the order of 15-20%, owing to appreciable star formation in highly dust-obscured galaxies, which might remain undetected in such UV observations.
We report the detection of a largely ionized very-high velocity cloud (VHVC; $v_{\rm LSR}\sim-350$ km/s) toward M33 with the Hubble Space Telescope/Cosmic Origin Spectrograph. The VHVC is detected in OI, CII, SiII, and SiIII absorption along five sightlines separated by ~0.06-0.4 degree. On sub-degree scales, the velocities and ionic column densities of the VHVC remain relatively smooth with standard deviations of +/-14 km/s and +/-0.15 dex between the sightlines, respectively. The VHVC has a metallicity of [OI/HI]=-0.56+/-0.17 dex (Z=0.28+/-0.11 Z$_{\odot}$). Despite the position-velocity proximity of the VHVC to the ionized Magellanic Stream, the VHVC's higher metallicity makes it unlikely to be associated with the Stream, highlighting the complex velocity structure of this region of sky. We investigate the VHVC's possible origin by revisiting its surrounding HI environment. We find that the VHVC may be: (1) a MW CGM cloud, (2) related to a nearby HI VHVC -- Wright's Cloud, or (3) connected to M33's northern warp. Furthermore, the VHVC could be a bridge connecting Wright's Cloud and M33's northern warp, which would make it a Magellanic-like structure in the halo of M33.
We present an updated version of the mass--metallicity relation (MZR) using integral field spectroscopy data obtained from 734 galaxies observed by the CALIFA survey. These unparalleled spatially resolved spectroscopic data allow us to determine the metallicity at the same physical scale ($\mathrm{R_{e}}$) for different calibrators. We obtain MZ relations with similar shapes for all calibrators, once the scale factors among them are taken into account. We do not find any significant secondary relation of the MZR with either the star formation rate (SFR) or the specific SFR for any of the calibrators used in this study, based on the analysis of the residuals of the best fitted relation. However we do see a hint for a (s)SFR-dependent deviation of the MZ-relation at low masses (M$<$10$^{9.5}$M$_\odot$), where our sample is not complete. We are thus unable to confirm the results by Mannucci et al. (2010), although we cannot exclude that this result is due to the differences in the analysed datasets. In contrast, our results are inconsistent with the results by Lara-Lopez et al. (2010), and we can exclude the presence of a SFR-Mass-Oxygen abundance Fundamental Plane. These results agree with previous findings suggesting that either (1) the secondary relation with the SFR could be induced by an aperture effect in single fiber/aperture spectroscopic surveys, (2) it could be related to a local effect confined to the central regions of galaxies, or (3) it is just restricted to the low-mass regime, or a combination of the three effects.
The abundance of H2 in molecular clouds, relative to the commonly used tracer CO, has only been measured toward a few embedded stars, which may be surrounded by atypical gas. We present observations of near-infrared absorption by H2, CO, and dust toward stars behind molecular clouds, providing a representative sample of these molecules in cold molecular gas, primarily in the Taurus Molecular Cloud. We find N_H2/A_V ~ 1.0x10^21 cm^-2, N_CO/A_V ~ 1.5x10^17 cm^-2 (1.8x10^17 including solid CO), and N_H2/N_CO ~ 6000. The measured N_H2/N_CO ratio is consistent with that toward embedded stars in various molecular clouds, but both are less than that derived from mm-wave observations of CO and star counts. The difference apparently results from the higher directly measured N_CO/A_V ratio.
We aim to characterize the relationship between dust properties. We also aim
to provide equations to estimate accurate dust properties from limited
observational datasets.
We assemble a sample of 1,630 nearby (z<0.1) galaxies-over a large range of
Mstar, SFR - with multi-wavelength observations available from wise, iras,
planck and/or SCUBA. The characterization of dust emission comes from SED
fitting using Draine & Li dust models, which we parametrize using two
components (warm and cold ). The subsample of these galaxies with global
measurements of CO and/or HI are used to explore the molecular and/or atomic
gas content of the galaxies.
The total Lir, Mdust and dust temperature of the cold component (Tc) form a
plane that we refer to as the dust plane. A galaxy's sSFR drives its position
on the dust plane: starburst galaxies show higher Lir, Mdust and Tc compared to
Main Sequence and passive galaxies. Starburst galaxies also show higher
specific Mdust (Mdust/Mstar) and specific Mgas (Mgas/Mstar). The Mdust is more
closely correlated with the total Mgas (atomic plus molecular) than with the
individual components. Our multi wavelength data allows us to define several
equations to estimate Lir, Mdust and Tc from one or two monochromatic
luminosities in the infrared and/or sub-millimeter.
We estimate the dust mass and infrared luminosity from a single monochromatic
luminosity within the R-J tail of the dust emission, with errors of 0.12 and
0.20dex, respectively. These errors are reduced to 0.05 and 0.10 dex,
respectively, if the Tc is used. The Mdust is correlated with the total Mism
(Mism \propto Mdust^0.7). For galaxies with Mstar 8.5<log(Mstar/Msun) < 11.9,
the conversion factor \alpha_850mum shows a large scatter (rms=0.29dex). The SF
mode of a galaxy shows a correlation with both the Mgass and Mdust: high
Mdust/Mstar galaxies are gas-rich and show the highest SFRs.
To constrain models of high-mass star formation, the Herschel/HOBYS KP aims
at discovering massive dense cores (MDCs) able to host the high-mass analogs of
low-mass prestellar cores, which have been searched for over the past decade.
We here focus on NGC6334, one of the best-studied HOBYS molecular cloud
complexes.
We used Herschel PACS and SPIRE 70-500mu images of the NGC6334 complex
complemented with (sub)millimeter and mid-infrared data. We built a complete
procedure to extract ~0.1 pc dense cores with the getsources software, which
simultaneously measures their far-infrared to millimeter fluxes. We carefully
estimated the temperatures and masses of these dense cores from their SEDs.
A cross-correlation with high-mass star formation signposts suggests a mass
threshold of 75Msun for MDCs in NGC6334. MDCs have temperatures of 9.5-40K,
masses of 75-1000Msun, and densities of 10^5-10^8cm-3. Their mid-IR emission is
used to separate 6 IR-bright and 10 IR-quiet protostellar MDCs while their 70mu
emission strength, with respect to fitted SEDs, helps identify 16 starless MDC
candidates. The ability of the latter to host high-mass prestellar cores is
investigated here and remains questionable. An increase in mass and density
from the starless to the IR-quiet and IR-bright phases suggests that the
protostars and MDCs simultaneously grow in mass. The statistical lifetimes of
the high-mass prestellar and protostellar core phases, estimated to be
1-7x10^4yr and at most 3x10^5yr respectively, suggest a dynamical scenario of
high-mass star formation.
The present study provides good mass estimates for a statistically
significant sample, covering the earliest phases of high-mass star formation.
High-mass prestellar cores may not exist in NGC6334, favoring a scenario
presented here, which simultaneously forms clouds and high-mass protostars.
We provide an updated calibration of CIV $\lambda1549$ broad emission line-based single-epoch (SE) black hole (BH) mass estimators for active galactic nuclei (AGNs) using new data for six reverberation-mapped AGNs at redshift $z=0.005-0.028$ with BH masses (bolometric luminosities) in the range $10^{6.5}-10^{7.5}$ $M_{\odot}$ ($10^{41.7}-10^{43.8}$ erg s$^{\rm -1}$). New rest-frame UV-to-optical spectra covering 1150-5700 \AA\ for the six AGNs were obtained with the Hubble Space Telescope (HST). Multi-component spectral decompositions of the HST spectra were used to measure SE emission-line widths for the CIV, MgII, and H$\beta$ lines as well as continuum luminosities in the spectral region around each line. We combine the new data with similar measurements for a previous archival sample of 25 AGNs to derive the most consistent and accurate calibrations of the CIV-based SE BH mass estimators against the H$\beta$ reverberation-based masses, using three different measures of broad-line width: full-width at half maximum (FWHM), line dispersion ($\sigma_{\rm line}$) and mean absolute deviation (MAD). The newly expanded sample at redshift $z=0.005-0.234$ covers a dynamic range in BH mass (bolometric luminosity) of $\log\ M_{\rm BH}/M_{\odot} = 6.5-9.1$ ($\log\ L_{\rm bol}/$erg s$^{\rm -1}=41.7-46.9$), and we derive the new CIV-based mass estimators using a Bayesian linear regression analysis over this range. We generally recommend the use of $\sigma_{\rm line}$ or MAD rather than FWHM to obtain a less biased velocity measurement of the CIV emission line, because its narrow-line component contribution is difficult to decompose from the broad-line profile.
We present a study of the trade-off between depth and resolution using a large number of U-band imaging observations in the GOODS-North field (Giavalisco et al. 2004) from the Large Binocular Camera (LBC) on the Large Binocular Telescope (LBT). Having acquired over 30 hours of data (315 images with 5-6 mins exposures), we generated multiple image mosaics, starting with the best atmospheric seeing images (FWHM $\lesssim$0.8"), which constitute $\sim$10% of the total data set. For subsequent mosaics, we added in data with larger seeing values until the final, deepest mosaic included all images with FWHM $\lesssim$1.8" ($\sim$94% of the total data set). From the mosaics, we made object catalogs to compare the optimal-resolution, yet shallower image to the lower-resolution but deeper image. We show that the number counts for both images are $\sim$90% complete to $U_{AB}$ $\lesssim26$. Fainter than $U_{AB}$$\sim$ 27, the object counts from the optimal-resolution image start to drop-off dramatically (90% between $U_{AB}$ = 27 and 28 mag), while the deepest image with better surface-brightness sensitivity ($\mu^{AB}_{U}$$\lesssim$ 32 mag arcsec$^{-2}$) show a more gradual drop (10% between $U_{AB}$ $\simeq$ 27 and 28 mag). For the brightest galaxies within the GOODS-N field, structure and clumpy features within the galaxies are more prominent in the optimal-resolution image compared to the deeper mosaics. Finally, we find - for 220 brighter galaxies with $U_{AB}$$\lesssim$ 24 mag - only marginal differences in total flux between the optimal-resolution and lower-resolution light-profiles to $\mu^{AB}_{U}$$\lesssim$ 32 mag arcsec$^{-2}$. In only 10% of the cases are the total-flux differences larger than 0.5 mag. This helps constrain how much flux can be missed from galaxy outskirts, which is important for studies of the Extragalactic Background Light.
Decades ago S. Lundquist, S. Chandrasekhar and R. J. Tayler first posed questions about the stability of Taylor-Couette flows of conducting material under the influence of large-scale background magnetic fields. These and many new questions can now be answered with numerical methods where the nonlinear simulations even provide the instability-induced turbulent values of several diffusivities. The cylindrical containers are here often assumed as axially unbounded and the background fields possess axial and/or azimuthal components. The influence of the magnetic Prandtl number Pm on the onset of the instabilities is demonstrated to be substantial. E.g., the potential flow subject to axial fields becomes unstable for a certain value of the averaged Reynolds number $\overline{Rm}=\sqrt{Re\cdot Rm}$, with Reynolds number Re and magnetic Reynolds number Rm. Rotation profiles as flat as the quasi-Keplerian rotation law scale similarly but large Pm. Among the instabilities of azimuthal fields, those where the background field has the same radial profile as the background flow are particularly interesting. They are unstable against nonaxisymmetric perturbations if at least one of the diffusivities is non-zero. For Pm$\gg 1$ the onset of the instability again scales with $\overline{Rm}$, but for Pm$\ll 1$ it scales with Re. Also with superrotation magnetic instabilities exist; they are of a double-diffusive character with Pm$\neq 1$. Several new experiments could be designed with liquid metals as the conducting fluid, starting with the experiment PROMISE which probed the interaction of twisted magnetic fields with differential rotation. Several of these experiments are described together with relevant diversifications of the magnetic instability theory including numerical studies of the kinetic and magnetic energies, the influence of the Hall effect and of axial density stratifications.
In this work we apply the stacking technique to estimate the average fractional polarisation from 30 to 353 GHz of a primary sample of 1560 compact sources - essentially all radio sources - detected in the 30 GHzPlanck all-sky map and listed in the second version of the Planck Catalogue of Compact Sources (PCCS2). We divide our primary sample in two subsamples according to whether the sources lay (679 sources) or not (881 sources) inside the sky region defined by the Planck Galactic mask (f_sky ~ 60%) and the area around the Magellanic Clouds. We find that the average fractional polarisation of compact sources is approximately constant (with frequency) in both samples (with a weighted mean over all the channels of 3.08% outside and 3.54% inside the Planck mask). In the sky region outside the adopted mask, we also estimate the {\mu} and {\sigma} parameters for the log-normal distribution of the mean fractional polarisation, finding a weighted mean value over all the Planck frequency range of 1.0 for {\sigma} and 0.7 for {\mu} (that would imply a weighted mean value for the median fractional polarisation of 1.9%).
We have constructed the database of stars in the local group using the extended version of the SAGA (Stellar Abundances for Galactic Archaeology) database that contains stars in 24 dwarf spheroidal galaxies and ultra faint dwarfs. The new version of the database includes more than 4500 stars in the Milky Way, by removing the previous metallicity criterion of [Fe/H] <= -2.5, and more than 6000 stars in the local group galaxies. We examined a validity of using a combined data set for elemental abundances. We also checked a consistency between the derived distances to individual stars and those to galaxies in the literature values. Using the updated database, the characteristics of stars in dwarf galaxies are discussed. Our statistical analyses of alpha-element abundances show that the change of the slope of the [alpha/Fe] relative to [Fe/H] (so-called "knee") occurs at [Fe/H] = -1.0 for the Milky Way. The knee positions for selected galaxies are derived by applying the same method. Star formation history of individual galaxies are explored using the slope of the cumulative metallicity distribution function. Radial gradients along the four directions are inspected in five galaxies where we find no direction dependence of metallicity gradients along the major and minor axis. The compilation of all the available data shows a lack of CEMP-s population in dwarf galaxies, while there may be some CEMP-no stars at [Fe/H] <~ -3 even in the very small sample. The inspection of the relationship between Eu and Ba abundances confirms an anomalously Ba-rich population in Fornax, which indicates a pre-enrichment of interstellar gas with r-process elements. We do not find any evidence of anti-correlations in O-Na and Mg-Al abundances, which characterises the abundance trends in the Galactic globular clusters.
We present the results of our search for the faint galaxies near the end of the Reionisation Epoch. This has been done using very deep OSIRIS images obtained at the Gran Telescopio Canarias (GTC). Our observations focus around two close, massive Lyman Alpha Emitters (LAEs) at redshift 6.5, discovered in the SXDS field within a large-scale overdense region (Ouchi et al. 2010). The total GTC observing time in three medium band filters (F883w35, F913w25 and F941w33) is over 34 hours covering $7.0\times8.5$ arcmin$^2$ (or $\sim30,000$ Mpc$^3$ at $z=6.5$). In addition to the two spectroscopically confirmed LAEs in the field, we have identified 45 other LAE candidates. The preliminary luminosity function derived from our observations, assuming a spectroscopic confirmation success rate of $\frac{2}{3}$ as in previous surveys, suggests this area is about 2 times denser than the general field galaxy population at $z=6.5$. If confirmed spectroscopically, our results will imply the discovery of one of the earliest protoclusters in the universe, which will evolve to resemble the most massive galaxy clusters today.
We investigate the origin of observed local star formation relations using radiative magnetohydrodynamic simulations with self-consistent star formation and ionising radiation. We compare these clouds to the density distributions of local star-forming clouds and find that the most diffuse simulated clouds match the observed clouds relatively well. We then compute both observationally-motivated and theoretically-motivated star formation efficiencies (SFEs) for these simulated clouds. By including ionising radiation, we can reproduce the observed SFEs in the clouds most similar to nearby Milky Way clouds. For denser clouds, the SFE can approach unity. These observed SFEs are typically 3 to 10 times larger than the "total" SFEs, i.e. the fraction of the initial cloud mass converted to stars. Converting observed to total SFEs is non-trivial. We suggest some techniques for doing so, though estimate up to a factor of ten error in the conversion.
Observations by HAWC and Milagro have detected bright and spatially extended TeV gamma-ray sources surrounding the Geminga and Monogem pulsars. We argue that these observations, along with a substantial population of other extended TeV sources coincident with pulsar wind nebulae, constitute a new morphological class of spatially extended TeV halos. We show that HAWCs wide field-of-view unlocks an expansive parameter space of TeV halos not observable by atmospheric Cherenkov telescopes. Under the assumption that Geminga and Monogem are typical middle-aged pulsars, we show that ten-year HAWC observations should eventually observe 37$^{+17}_{-13}$ middle-aged TeV halos that correspond to pulsars whose radio emission is not beamed towards Earth. Depending on the extrapolation of the TeV halo efficiency to young pulsars, HAWC could detect more than 100 TeV halos from mis-aligned pulsars. These pulsars have historically been difficult to detect with existing multiwavelength observations. TeV halos will constitute a significant fraction of all HAWC sources, allowing follow-up observations to efficiently find pulsar wind nebulae and thermal pulsar emission. The observation and subsequent multi-wavelength follow-up of TeV halos will have significant implications for our understanding of pulsar beam geometries, the evolution of PWN, the diffusion of cosmic-rays near energetic pulsars, and the contribution of pulsars to the cosmic-ray positron excess.
Nuclear starburst discs (NSDs) are star-forming discs that may be residing in the nuclear regions of active galaxies at intermediate redshifts. One dimensional (1D) analytical models developed by Thompson et al. (2005) show that these discs can possess an inflationary atmosphere when dust is sublimated on parsec scales. This make NSDs a viable source for AGN obscuration. We model the two dimensional (2D) structure of NSDs using an iterative method in order to compute the explicit vertical solutions for a given annulus. These solutions satisfy energy and hydrostatic balance, as well as the radiative transfer equation. In comparison to the 1D model, the 2D calculation predicts a less extensive expansion of the atmosphere by orders of magnitude at the parsec/sub-parsec scale, but the new scale-height $h$ may still exceed the radial distance $R$ for various physical conditions. A total of 192 NSD models are computed across the input parameter space in order to predict distributions of a line of sight column density $N_H$. Assuming a random distribution of input parameters, the statistics yield 56% of Type 1, 23% of Compton-thin Type 2s (CN), and 21% of Compton-thick (CK) AGNs. Depending on a viewing angle ($\theta$) of a particular NSD (fixed physical conditions), any central AGN can appear to be Type 1, CN, or CK which is consistent with the basic unification theory of AGNs. Our results show that $\log[N_H(\text{cm}^{-2})]\in$ [23,25.5] can be oriented at any $\theta$ from 0$^\circ$ to $\approx$80$^\circ$ due to the degeneracy in the input parameters.
Ark 564 is an archetypal Narrow line Seyfert 1 that has been well observed in soft X-rays from 0.3-10 keV, revealing a steep spectrum, strong soft excess, iron K emission line and dramatic variability on the order of hours. Because of its very steep spectrum, observations of the source above 10 keV have been sparse. We report here on the first NuSTAR observation of Ark 564. The source was observed for 200 ks with NuSTAR, 50 ks of which were concurrent with Suzaku observations. NuSTAR and Suzaku observed a dramatic flare, in which the hard emission is clearly delayed with respect to the soft emission, consistent with previous detections of a low-frequency hard lag found in XMM-Newton data. The NuSTAR spectrum is well described by a low-temperature Comptonisation continuum (with an electron temperature of 15 +/- 2 keV), which irradiates a highly ionised disc. No further relativistic broadening or ionized absorption is required. These spectral results show that Ark 564 has one of the lowest temperature coronae observed by NuSTAR to date. We discuss possible reasons for low-temperature coronae in high-Eddington sources.
The spin of Wolf-Rayet (WR) stars at low metallicity (Z) is most relevant for our understanding of gravitational wave sources such as GW 150914, as well as the incidence of long-duration gamma-ray bursts (GRBs). Two scenarios have been suggested for both phenomena: one of them involves rapid rotation and quasi-chemical homogeneous evolution (CHE), the other invokes classical evolution through mass loss in single and binary systems. WR spin rates might enable us to test these two scenarios. In order to obtain empirical constraints on black hole progenitor spin, we infer wind asymmetries in all 12 known WR stars in the Small Magellanic Cloud (SMC) at Z = 1/5 Zsun, as well as within a significantly enlarged sample of single and binary WR stars in the Large Magellanic Cloud (LMC at Z = 1/2 Zsun), tripling the sample of Vink (2007). This brings the total LMC sample to 39, making it appropriate for comparison to the Galactic sample. We measure WR wind asymmetries with VLT-FORS linear spectropolarimetry. We report the detection of new line effects in the LMC WN star BAT99-43 and the WC star BAT99-70, as well as the famous WR/LBV HD 5980 in the SMC, which might be evolving chemically homogeneously. With the previous reported line effects in the late-type WNL (Ofpe/WN9) objects BAT99-22 and BAT99-33, this brings the total LMC WR sample to 4, i.e. a frequency of ~10%. Perhaps surprisingly, the incidence of line effects amongst low-Z WR stars is not found to be any higher than amongst the Galactic WR sample, challenging the rotationally-induced CHE model. As WR mass loss is likely Z-dependent, our Magellanic Cloud line-effect WR stars may maintain their surface rotation and fulfill the basic conditions for producing long GRBs, both via the classical post-red supergiant (RSG) or luminous blue variable (LBV) channel, as well as resulting from CHE due to physics specific to very massive stars (VMS).
We present a deep radio search in the Reticulum II dwarf spheroidal (dSph) galaxy performed with the Australia Telescope Compact Array. Observations were conducted at 16 cm wavelength, with an rms sensitivity of 0.01 mJy/beam, and with the goal of searching for synchrotron emission induced by annihilation or decay of weakly interacting massive particles (WIMPs). Data were complemented with observations on large angular scales taken with the KAT-7 telescope. We find no evidence for a diffuse emission from the dSph and we derive competitive bounds on the WIMP properties. In addition, we detect more than 200 new background radio sources. Among them, we show there are two compelling candidates for being the radio counterpart of the possible gamma-ray emission reported by other groups using Fermi-LAT data.
Star-forming galaxies at high redshifts are the ideal targets to probe the hypothetical variation of the fine-structure constant alpha over cosmological time scales. We propose a modification of the alkali doublets method which allows us to search for variation in alpha combining far infrared and submillimeter spectroscopic observations. This variation manifests as velocity offsets between the observed positions of the fine-structure and gross-structure transitions when compared to laboratory wavelengths. Here we describe our method whose sensitivity limit to the fractional changes in alpha is about 5*10^-7. We also demonstrate that current spectral observations of hydrogen and [C II] 158 micron lines provide an upper limit on |Delta alpha/alpha| < 6*10^-5 at redshifts z = 3.1 and z = 4.7.
Links to: arXiv, form interface, find, astro-ph, recent, 1703, contact, help (Access key information)
We report the discovery of the molecular cloud whose kinematic distance is the largest in the Galaxy at the present moment, named G213.042$+$0.003, at $l =$ 213.042$^\circ$ and $b =$ 0.003$^\circ$ in the $^{12}$CO($J =$ 1-0) line using the Nobeyama 45-m telescope and a multi-beam receiver BEARS. This molecular cloud is located at the heliocentric distance of 21$_{-7}^{+12}$ kpc and Galactocentric distance of 29$_{-7}^{+12}$ kpc, which are estimated as the kinematic distances with the Galactic parameters obtained by Reid et al. (2014, ApJ, 783, 130). Its major and minor diameters and line width were measured to be 4.0$_{-1.3}^{+2.3}$ pc, 3.0$_{-1.0}^{+1.7}$ pc, and 1.5 km s$^{-1}$, respectively. The cloud mass was estimated to be 2.5$_{-1.4}^{+3.7}$ $\times$ 10$^2$ $M_\odot$ using the CO-to-H$_2$ conversion factor of 5.6 $\times$ 10$^{20}$ cm$^{-2}$ (K km s$^{-1}$)$^{-1}$ obtained in far outer Galaxy.
We present near-infrared and optical spectroscopic observations of a sample of 450$\mu$m and 850$\mu$m-selected dusty star-forming galaxies (DSFGs) identified in a 400 arcmin$^2$ area in the COSMOS field. Thirty-one sources of the 102 targets were spectroscopically confirmed at $0.2<z<4$, identified primarily in the near-infrared with Keck MOSFIRE and some in the optical with Keck LRIS and DEIMOS. The low rate of confirmation is attributable both to high rest-frame optical obscuration in our targets and limited sensitivity to certain redshift ranges. The high-quality photometric redshifts available in the COSMOS field allow us to test the robustness of photometric redshifts for DSFGs. We find a subset (11/31$\approx35$%) of DSFGs with inaccurate ($\Delta z/(1+z)>0.2$) or non-existent photometric redshifts; these have very distinct spectral energy distributions from the remaining DSFGs, suggesting a decoupling of highly obscured and unobscured components. We present a composite rest-frame 4300--7300\AA\ spectrum for DSFGs, and find evidence of 200$\pm$30 km s$^{-1}$ gas outflows. Nebular line emission for a sub-sample of our detections indicate that hard ionizing radiation fields are ubiquitous in high-z DSFGs, even more so than typical mass or UV-selected high-z galaxies. We also confirm the extreme level of dust obscuration in DSFGs, measuring very high Balmer decrements, and very high ratios of IR to UV and IR to H$\alpha$ luminosities. This work demonstrates the need to broaden the use of wide bandwidth technology in the millimeter to the spectroscopic confirmations of large samples of high-z DSFGs, as the difficulty in confirming such sources at optical/near-infrared wavelengths is exceedingly challenging given their obscuration.
While traditionally associated with active galactic nuclei (AGN), the properties of the CII], CIII] and CIV emission lines are still uncertain as large, unbiased samples of sources are scarce. We present the first blind, statistical study of CII], CIII] and CIV emitters at $z\sim0.68,1.05,1.53$, respectively, uniformly selected down to a flux limit of $\sim4\times10^{-17}$ erg s$^{-1}$ cm$^{-1}$ through a narrow band survey covering an area of $\sim1.4$ deg$^2$ over COSMOS and UDS. We detect 16 CII], 35 CIII] and 17 CIV emitters, whose nature we investigate using optical colours as well as HST, X-ray, radio and far infra-red data. We find that $z\sim0.7$ CII] emitters are consistent with a mixture of blue (UV slope $\beta=-2.0\pm0.4$) star forming galaxies with disky HST structure and AGN with Seyfert-like morphologies. Bright CII] emitters have individual X-ray detections as well as high average black hole accretion rates (BHAR) of $\sim0.1$ $M_{\odot}$ yr$^{-1}$. CIII] emitters at $z\sim1.05$ trace a general population of SF galaxies, with $\beta=-0.8\pm1.1$, a variety of optical morphologies, including isolated and interacting galaxies and low BHAR ($<0.02$ $M_{\odot}$ yr$^{-1}$). Our CIV emitters at $z\sim1.5$ are consistent with young, blue quasars ($\beta\sim-1.9$) with point-like optical morphologies, bright X-ray counterparts and large BHAR ($0.8$ $M_{\odot}$ yr$^{-1}$). We also find some surprising CII], CIII] and CIV emitters with rest-frame equivalent widths which could be as large as $50-100$ {\AA}. AGN or spatial offsets between the UV continuum stellar disk and the line emitting regions may explain the large EW. These bright CII], CIII] and CIV emitters are ideal candidates for spectroscopic follow up to fully unveil their nature.
Utilizing spectroscopic observations taken for the VIMOS Ultra-Deep Survey (VUDS) we report here on the discovery of Cl J1001+0220, a massive proto-cluster located at $z_{spec}\sim4.57$ in the COSMOS field. The proto-cluster was initially detected as a $\sim12\sigma$ overdensity of typical star-forming galaxies in the blind spectroscopic survey of the early universe ($2<z<6$) performed by VUDS. It was further mapped using a new technique developed which statistically combines spectroscopic and photometric redshifts, the latter derived from a recent compilation of deep multi-band imaging. Through various methods, the descendant halo mass of Cl J1001+0220 is estimated to be $\log(M_{h}/M_{\odot})_{z=0}\sim14.5-15$ with a large amount of mass apparently already in place at $z\sim4.57$. Tentative evidence is found for a fractional excess of older and more massive galaxies within the proto-cluster, an observation which suggests the pervasive early onset of vigorous star formation. No evidence is found for the differences in the star formation rates of member and a matched sample of coeval field galaxies either through rest-frame ultraviolet methods or through stacking extremely deep Very Large Array 3 GHz imaging. Additionally, no evidence for pervasive strong active galactic nuclei (AGN) activity is observed. Analysis of Hubble Space Telescope images provides weak evidence for an elevated incidence of galaxy-galaxy interaction within the proto-cluster. The spectral properties of the two samples are compared, with a definite suppression of Ly$\alpha$ seen in the average member galaxy relative to the coeval field ($f_{esc,Ly\alpha}=1.8^{+0.3}_{-1.7}$% and $4.0^{+1.0}_{-0.8}$%, respectively). This observation along with other lines of evidence leads us to infer the possible presence of a large, cool diffuse medium within the proto-cluster environment evocative of a nascent intracluster medium.
We study the correlation between the specific star formation rate of central galaxies and neighbour galaxies, also known as 'galactic conformity', out to several Mpc scales using three different semi-analytic models (SAMs). It has been suggested that SAMs may not show the strong signal of conformity measured in observations. In all the models, when the selection of primary galaxies is based on an isolation criterion in real space, we measure a strong galactic conformity where the mean quenched fraction around quenched primary galaxies is higher than that around star-forming primary galaxies of the same stellar mass. The cumulative signal-to-noise ratio, which is used to evaluate the significance and extension in distance of conformity, depends strongly on the SAM. The two-halo conformity is significant as far as ~5 Mpc/h. The overall signal of galactic conformity decreases when we remove satellites in the selection of primary galaxies. In this case, there is no significant two-halo conformity in two SAMs. For the other SAM, the two-halo conformity is only detected between isolated, central galaxies in relatively low-mass haloes ($M_{\rm halo}$ $< 10^{12.4}$ $h^{-1}~\rm{ M_{\bigodot}}$) and their neighbours. Finally, we explore the case when the secondary galaxies correspond to central galaxies in the SAMs as an attempt to measure the conformity between distinct haloes. Our results show that there is no conformity beyond ~3 Mpc/h. The conformity measured out to scales of 3$-$4 Mpc/h can be related to the environmental influence of large galaxy clusters around near central galaxies.
We present spatially resolved ALMA observations of the CO J=3-2 emission line in two massive galaxies at z=2.5 on the star-forming main sequence. Both galaxies have compact dusty star-forming cores with effective radii of Re=1.3 kpc and Re=1.2 kpc in the 870 um continuum emission. The spatial extent of star-forming molecular gas is also compact with Re=1.9 kpc and Re=2.3 kpc, but more extended than the dust emission. Interpreting the observed position-velocity diagrams with dynamical models, we find the starburst cores to be rotation-dominated with the ratio of the maximum rotation velocity to the local velocity dispersion of v/sigma=7.0 (v=386 km/s) and v/sigma_0=4.1 (v=391 km/s). The comparisons among dynamical, stellar, gas, and dust mass suggest that the starburst CO-to-H2 conversion factor of alpha_CO=0.8 Msun/(K km/s/pc2) is appropriate in the spatially resolved cores. The dense cores are likely to be formed in extreme environments similar to the central regions of local ultraluminous infrared galaxies. Given their short gas depletion timescale of 46 Myr and 79 Myr, they are likely to be the immediate progenitors of bulge-dominated quiescent galaxies. Our work also demonstrates that a combination of medium-resolution CO and high-resolution dust continuum observations is a powerful tool for characterizing the dynamical state of molecular gas in distant galaxies.
Hubble Space Telescope photometry from the ACS/WFC and WFPC2 cameras is used to detect and measure globular clusters (GCs) in the central region of the rich Perseus cluster of galaxies. A detectable population of Intragalactic GCs is found extending out to at least 500 kpc from the cluster center. These objects display luminosity and color (metallicity) distributions that are entirely normal for GC populations. Extrapolating from the limited spatial coverage of the HST fields, we estimate very roughly that the entire Perseus cluster should contain ~50000 or more IGCs, but a targetted wide-field survey will be needed for a more definitive answer. Separate brief results are presented for the rich GC systems in NGC 1272 and NGC 1275, the two largest Perseus ellipticals. For NGC 1272 we find a specific frequency S_N = 8, while for the central giant NGC 1275, S_N ~ 12. In both these giant galaxies, the GC colors are well matched by bimodal distributions, with the majority in the blue (metal-poor) component. This preliminary study suggests that Perseus is a prime target for a more comprehensive deep imaging survey of Intragalactic GCs.
We model the cluster A1689 in two modified MOND frameworks (EMOND and what we call generalised MOND or GMOND) with the aim of determining whether it is possible to explain the inferred acceleration profile, from gravitational lensing, without the aid of dark matter. We also compare our result to predictions from MOG/STVG and Emergent Gravity. By using a baryonic mass model, we determine the total gravitational acceleration predicted by the modified gravitational equations and compare the result to NFW profiles of dark matter studies on A1689 from the literature. Theory parameters are inferred empirically, with the aid of previous work. We are able to reproduce the desired acceleration profile of A1689 for GMOND, EMOND and MOG, but not Emergent Gravity. There is much more work which needs to be conducted with regards to understanding how the GMOND parameters behave in different environments. Furthermore, we show that the exact baryonic profile becomes very important when undertaking modified gravity modelling rather than {\Lambda}CDM.
Using observations from the first two years of the MOSFIRE Deep Evolution Field (MOSDEF) survey, we study 13 AGN-driven outflows detected from a sample of 67 X-ray, IR and/or optically selected AGN at $z \sim 2$. The AGN in the sample have bolometric luminosities spanning $10^{44}-10^{46} \mathrm{erg s^{-1}}$, including both quasars and more common, moderate-luminosity AGN. We detect signatures of blueshifted, ionized gas outflows in the \hbeta , [OIII], \halpha and/or [NII] emission lines of 19\% of the AGN, while only 1.8\% of the MOSDEF galaxies have similar outflows detected. The outflow velocities range from 300 to 1000 km s$^{-1}$. Eight of the 13 detected outflows are spatially extended on similar scales as the host galaxies, with spatial extents of 2.5 to 11.0 kpc along the MOSFIRE slits. Outflows are detected uniformly across the galaxy star-forming main sequence, showing little trend with the stellar mass or SFR of the host galaxy. Line ratio diagnostics in the BPT diagram indicate that the outflowing gas is photoionized by the AGN. We do not find evidence for positive AGN feedback, either in our relatively small MOSDEF sample or in a much larger sample from SDSS. Analysis of the mass and energy outflow rates in the MOSDEF AGN indicates that stellar feedback is likely insufficient to drive these outflows at $z \sim 2$ and that they are very likely to be AGN-driven. The outflows have mass-loading factors of $\sim 1-5$, suggesting that they help regulate star formation in their host galaxies, though they may not be sufficient to fully quench it.
We present observations and analysis of the low-mass star-forming region, Taurus Molecular Cloud-1 (TMC-1). CS ($J$=2-1)/N$_2$H$^+$ ($J$=1-0) and C$^{17}$O ($J$=2-1)/C$^{18}$O ($J$=2-1) were observed with FCRAO (Five College Radio Astronomy Observatory) and SRAO (Seoul Radio Astronomy Observatory), respectively. In addition, Spitzer infrared data and 1.2 mm continuum data observed with MAMBO (Max-Planck Millimetre Bolometer) are used. We also perform chemical modeling to investigate the relative molecular distributions of the TMC-1 filament. Based on Spitzer observations, there is no young stellar object along the TMC-1 filament, while five Class II and one Class I young stellar objects are identified outside the filament. The comparison between column densities calculated from dust continuum and C$^{17}$O 2-1 line emission shows that CO is depleted much more significantly in the ammonia peak than in the cyanopolyyne peak, while the column densities calculated from the dust continuum are similar at the two peaks. N$_2$H$^+$ is not depleted much in either peak. According to our chemical calculation, the differential chemical distribution in the two peaks can be explained by different timescales required to reach the same density, i.e., by different dynamical processes.
A fraction of early-type dwarf galaxies in the Virgo cluster have a disk component and even possess disk features such as bar, lens, and spiral arms. In this study, we construct fifteen galaxy models that resemble VCC856, considered to be an infalling progenitor of disk dwarf galaxies, within observational error ranges, and use $N$-body simulations to study their long-term dynamical evolution in isolation as well as the formation of bar in them. We find that dwarf disk galaxies readily form bars unless they have an excessively concentrated halo or a hot disk. This suggests that infalling dwarf disk galaxies are intrinsically unstable to bar formation, even without any external perturbation, accounting for a population of barred dwarf galaxies in the outskirts of the Virgo cluster. The bars form earlier and stronger in galaxies with a lower fraction of counter-streaming motions, lower halo concentration, lower velocity anisotropy, and thinner disk. Similarly to normal disk galaxies, dwarf disk galaxies also undergo recurrent buckling instabilities. The first buckling instability tends to shorten the bar and to thicken the disk, and drives a dynamical transition in the bar pattern speed as well as mass inflow rate. In nine models, the bars regrow after the mild first buckling instability due to the efficient transfer of disk angular momentum to the halo, and are subject to recurrent buckling instabilities to turn into X-shaped bulges.
The 2.1 GHz radio source catalogue of the 25 deg$^2$ ultimate XMM extragalactic survey south (XXL-S) field, observed with the Australia Telescope Compact Array (ATCA), is presented. The final radio mosaic achieved a resolution of $\sim$$4.8"$ and a median rms noise of $\sigma \approx41$ $\mu$Jy/beam. To date, this is the largest area radio survey to reach this flux density level. A total of 6350 radio components above 5$\sigma$ are included in the component catalogue, 26.4% of which are resolved. Of these components, 111 were merged together to create 48 multiple-component radio sources, resulting in a total of 6287 radio sources in the source catalogue, 25.9% of which were resolved. A survival analysis revealed that the median spectral index of the Sydney University Molonglo Sky Survey (SUMSS) 843 MHz sources in the field is $\alpha$ = $-$0.75, consistent with the values of $-0.7$ to $-0.8$ commonly used to characterise radio spectral energy distributions of active galactic nuclei (AGN). The 2.1 GHz and 1.4 GHz differential radio source counts are presented and compared to other 1.4 GHz radio surveys. The XXL-S source counts show good agreement with the other surveys.
We surveyed a central 0.6 deg^2 region of Abell 1367 cluster for extended ionized gas clouds (EIGs) using the Subaru prime-focus camera (Suprime-Cam) with a narrow-band filter that covers Halpha. We discovered six new EIGs in addition to five known EIGs. We also found that the Halpha tail from the blue infalling group (BIG) is extended to about 330 kpc in projected distance, which is about twice longer than previously reported. Candidates of star-forming blobs in the tail are detected. The properties of the EIG parent galaxies in Abell 1367 basically resemble those in the Coma cluster. A noticeable difference is that the number of detached EIGs is significantly fewer in Abell 1367, while the fraction of blue member galaxies is higher. The results suggest a difference in the evolutionary stage of the clusters; Abell 1367 is at an earlier stage than the Coma cluster.
We measure the mass function for a sample of 840 young star clusters with ages between 10-300 Myr observed by the Panchromatic Hubble Andromeda Treasury (PHAT) survey in M31. The data show clear evidence of a high-mass truncation: only 15 clusters more massive than $10^4$ $M_{\odot}$ are observed, compared to $\sim$100 expected for a canonical $M^{-2}$ pure power-law mass function with the same total number of clusters above the catalog completeness limit. Adopting a Schechter function parameterization, we fit a characteristic truncation mass of $M_c = 8.5^{+2.8}_{-1.8} \times 10^3$ $M_{\odot}$. While previous studies have measured cluster mass function truncations, the characteristic truncation mass we measure is the lowest ever reported. Combining this M31 measurement with previous results, we find that the cluster mass function truncation correlates strongly with the characteristic star formation rate surface density of the host galaxy, where $M_c \propto$ $\langle \Sigma_{\mathrm{SFR}} \rangle^{\sim1.1}$. We also find evidence that suggests the observed $M_c$-$\Sigma_{\mathrm{SFR}}$ relation also applies to globular clusters, linking the two populations via a common formation pathway. If so, globular cluster mass functions could be useful tools for constraining the star formation properties of their progenitor host galaxies in the early Universe.
Observations at 37 GHz, performed at Mets\"ahovi Radio Observatory, are presented for a sample of 78 radio-loud and radio-quiet narrow-line Seyfert 1 (NLS1) galaxies, together with additional lower and higher frequency radio data from RATAN-600, Owens Valley Radio Observatory, and the Planck satellite. Most of the data have been gathered between February 2012 and April 2015 but for some sources even longer lightcurves exist. The detection rate at 37 GHz is around 19%, comparable to other populations of active galactic nuclei presumed to be faint at radio frequencies, such as BL Lac objects. Variability and spectral indices are determined for sources with enough detections. Based on the radio data, many NLS1 galaxies show a blazar-like radio spectra exhibiting significant variability. The spectra at a given time are often inverted or convex. The source of the high-frequency radio emission in NLS1 galaxies, detected at 37 GHz, is most probably a relativistic jet rather than star formation. Jets in NLS1 galaxies are therefore expected to be a much more common phenomenon than earlier assumed.
Galaxy evolution and star formation are two multi-scale problems tightly linked to each other. To understand the interstellar cycle, which triggers galaxy evolution, it is necessary to describe simultaneously the large-scale evolution widely induced by the feedback processes and the details of the gas dynamics that controls the star formation process through gravitational collapse. We perform a set of three-dimensional high-resolution numerical simulations of a turbulent, self-gravitating and magnetized interstellar medium within a $1\ \mathrm{kpc}$ stratified box with supernova feedback correlated with star-forming regions. In particular, we focus on the role played by the magnetic field and the feedback on the galactic vertical structure, the star formation rate (SFR) and the flow dynamics. For this purpose we vary their respective intensities. We extract properties of the dense clouds arising from the turbulent motions and compute power spectra of various quantities. Using a distribution of supernovae sufficiently correlated with the dense gas, we find that supernova explosions can reproduce the observed SFR, particularly if the magnetic field is on the order of a few $\mu G$. The vertical structure, which results from a dynamical and an energy equilibrium is well reproduced by a simple analytical model, which allows us to estimate the coupling between the gas and the supernovae. We found the coupling to be rather low and on the order of 1.5$\%$. Strong magnetic fields may help to increase this coupling by a factor of about 2-3. To characterize the flow we compute the power spectra of various quantities in 3D but also in 2D in order to account for the stratification of the galactic disc.
The signatures for the existence of dark matter are revealed only through its gravitational interaction. Theoretical arguments support that the Weakly Interacting Massive Particle (WIMP) can be a class of dark matter and it can annihilate and/or decay to Standard Model particles, among which neutrino is a favorable candidate. We show that the proposed 50 kt Magnetized Iron CALorimeter (MagICAL) detector under the India-based Neutrino Observatory (INO) project can play an important role in the indirect searches of Galactic diffuse dark matter in the neutrino and antineutrino mode separately. We present the sensitivity of 500 kt$\cdot$yr MagICAL detector to set limits on the velocity-averaged self-annihilation cross-section ($\langle\sigma v\rangle$) and decay lifetime ($\tau$) of dark matter having mass in the range of 2 GeV $\leq m_\chi \leq $ 90 GeV and 4 GeV $\leq m_\chi \leq $ 180 GeV respectively, assuming no excess over the conventional atmospheric neutrino and antineutrino fluxes at the INO site. Our limits for low mass dark matter constrain the parameter space which has not been explored before. We show that MagICAL will be able to set competitive constraints, $\langle\sigma v\rangle\leq 1.87\times 10^{-24}$ cm$^3$ s$^{-1}$ for $\chi\chi\rightarrow\nu\bar\nu$ and $\tau\geq 4.8\times 10^{24}$ s for $\chi\rightarrow\nu\bar\nu$ at 90$\%$ C.L. (1 d.o.f.) for $m_\chi$ = 10 GeV.
Context: Understanding how accretion proceeds is a key question of star
formation, with important implications for both the physical and chemical
evolution of young stellar objects. In particular, very little is known about
the accretion variability in the earliest stages of star formation.
Aims: To characterise protostellar accretion histories towards individual
sources by utilising sublimation and freeze-out chemistry of CO.
Methods: A sample of 24 embedded protostars are observed with the
Submillimeter Array (SMA) in context of the large program "Mass Assembly of
Stellar Systems and their Evolution with the SMA" (MASSES). The size of the
C$^{18}$O emitting region, where CO has sublimated into the gas-phase, is
measured towards each source and compared to the expected size of the region
given the current luminosity. The SMA observations also include 1.3 mm
continuum data, which are used to investigate whether a link can be established
between accretion bursts and massive circumstellar disks.
Results: Depending on the adopted sublimation temperature of the CO ice,
between 20% and 50% of the sources in the sample show extended C$^{18}$O
emission indicating that the gas was warm enough in the past that CO sublimated
and is currently in the process of refreezing; something which we attribute to
a recent accretion burst. Given the fraction of sources with extended C$^{18}$O
emission, we estimate an average interval between bursts of 20000-50000 yr,
which is consistent with previous estimates. No clear link can be established
between the presence of circumstellar disks and accretion bursts, however the
three closest known binaries in the sample (projected separations <20 AU) all
show evidence of a past accretion burst, indicating that close binary
interactions may also play a role in inducing accretion variability.
Supermassive binary black holes (BBHs) on sub-parsec scales are prime targets for gravitational wave experiments. They also provide insights on close binary evolution and hierarchical structure formation. Sub-parsec BBHs cannot be spatially resolved but indirect methods can identify candidates. In 2015 Liu et al. reported an optical-continuum periodicity in the quasar PSO J334.2028+01.4075, with the estimated mass and rest-frame period suggesting an orbital separation of about 0.006 pc (0.7 micro-arcsec). The persistence of the quasar's optical periodicity has recently been disfavored over an extended baseline. However, if a radio jet is launched from a sub-parsec BBH, the binary's properties can influence the radio structure on larger scales. Here, we use the Very Long Baseline Array (VLBA) and Karl G. Jansky Very Large Array (VLA) to study the parsec- and kiloparsec-scale emission energized by the quasar's putative BBH. We find two VLBA components separated by 3.6 mas (30 pc), tentatively identifying one as the VLBA "core" from which the other was ejected. The VLBA components contribute to a point-like, time-variable VLA source that is straddled by lobes spanning 8 arcsec (66 kpc). We classify PSO J334.2028+01.4075 as a lobe-dominated quasar, albeit with an atypically large twist of 39 deg between its elongation position angles on parsec- and kiloparsec-scales. By analogy with 3C 207, a well-studied lobe-dominated quasar with a similarly-rare twist, we speculate that PSO J334.2028+01.4075 could be ejecting jet components over an inner cone that traces a precessing jet in a BBH system.
The fraction of binary stars (fb) is one of most valuable tool to probe the star formation and evolution of multiple systems in the Galaxy. We focus on the relationship between fb and stellar metallicity ([Fe/H]) by employing the differential radial velocity (DRV) method and the large sample observed by the Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST). Main-sequence stars from A- to K-types in the third data release (DR3) of LAMOST are selected to estimate fb. Contributions to a profile of DRV from radial velocity (RV) error of single stars (\sigma) and orbital motion of binary stars are evaluated from the profile of DRV. Finally, we employ 365,911 stars with randomly repeating spectral observations to present a detailed analysis of fb and \sigma in the two-dimensional (2D) space of Teff and [Fe/H]. The A-type stars are more likely to be companions in binary star systems than other stars. Furthermore, the reverse correlation between fb and [Fe/H] can be shown statistically, which suggests that fb is a joint function of Teff and [Fe/H]. At the same time, \sigma of the sample for different Teff and [Fe/H] are fitted. Metal-rich cold stars in our sample have the best RV measurement.
Links to: arXiv, form interface, find, astro-ph, recent, 1703, contact, help (Access key information)