The recent detection of two faint and extended star clusters in the central regions of two Local Group dwarf galaxies, Eridanus II and Andromeda XXV, raises the question of whether clusters with such low densities can survive the tidal field of cold dark matter haloes with central density cusps. Using both analytic arguments and a suite of collisionless N-body simulations, I show that these clusters are extremely fragile and quickly disrupted in the presence of central cusps $\rho\sim r^{-\alpha}$ with $\alpha\gtrsim 0.2$. Furthermore, the scenario in which the clusters where originally more massive and sank to the center of the halo requires extreme fine tuning and does not naturally reproduce the observed systems. In turn, these clusters are long lived in cored haloes, whose central regions are safe shelters for $\alpha\lesssim 0.2$. The only viable scenario for hosts that have preserved their primoridal cusp to the present time is that the clusters formed at rest at the bottom of the potential, which is easily tested by measurement of the clusters proper velocity within the host. This offers means to readily probe the central density profile of two dwarf galaxies as faint as $L_V\sim5\times 10^5 L_\odot$ and $L_V\sim6\times10^4 L_\odot$, in which stellar feedback is unlikely to be effective.
We present dynamical measurements from the KMOS (K-band Multi-Object Spectrograph) Deep Survey (KDS), which is comprised of 77 typical star-forming galaxies at z $\simeq$ 3.5 in the mass range 9.0 < log(M$_{\star}$/M$_{\odot}$) < 10.5. These measurements constrain the internal dynamics, the intrinsic velocity dispersions ({\sigma}$_{int}$) and rotation velocities (V$_{C}$) of galaxies in the high redshift Universe. The mean velocity dispersion of the galaxies in our sample is {\sigma}$_{int}$ = $70.8^{+3.3}_{-3.1}$ km s$^{-1}$, revealing that the increasing average {\sigma}$_{int}$ with increasing redshift, reported for z $\lesssim2$, continues out to z $\simeq$ 3.5. Only 34 $\pm$ 8% of our galaxies are rotation-dominated (V$_{C}$/{\sigma}$_{int}$ > 1), with the sample average V$_{C}$/{\sigma}$_{int}$ value much smaller than at lower redshift. After carefully selecting comparable star-forming samples at multiple epochs, we find that the rotation-dominated fraction evolves with redshift with a z$^{-0.2}$ dependence. The rotation-dominated KDS galaxies show no clear offset from the local rotation velocity-stellar mass (i.e. V$_{C}$-M$_{\star}$) relation, although a smaller fraction of the galaxies are on the relation due to the increase in the dispersion-dominated fraction. These observations are consistent with a simple equilibrium model picture, in which random motions are boosted in high redshift galaxies by a combination of the increasing gas fractions, accretion efficiency, specific star-formation rate and stellar feedback and which may provide significant pressure support against gravity on the galactic disk scale.
The velocity anisotropy parameter, beta, is a measure of the kinematic state of orbits in the stellar halo which holds promise for constraining the merger history of the Milky Way (MW). We determine global trends for beta as a function of radius from three suites of simulations, including accretion only and cosmological hydrodynamic simulations. We find that both types of simulations are consistent and predict strong radial anisotropy (<beta>~0.7) for Galactocentric radii greater than 10 kpc. Previous observations of beta for the MW's stellar halo claim a detection of an isotropic or tangential "dip" at r~20 kpc. Using the N-body+SPH simulations, we investigate the temporal persistence, population origin, and severity of "dips" in beta. We find dips in the in situ stellar halo are long-lived, while dips in the accreted stellar halo are short-lived and tied to the recent accretion of satellite material. We also find that a major merger as early as z~1 can result in a present day low (isotropic to tangential) value of beta over a wide range of radii. While all of these mechanisms are plausible drivers for the beta dip observed in the MW, in the simulations, each mechanism has a unique metallicity signature associated with it, implying that future spectroscopic surveys could distinguish between them. Since an accurate knowledge of beta(r) is required for measuring the mass of the MW halo, we note significant transient dips in beta could cause an overestimate of the halo's mass when using spherical Jeans equation modeling.
Low radio frequency surveys are important for testing unified models of radio-loud quasars and radio galaxies. Intrinsically similar sources that are randomly oriented on the sky will have different projected linear sizes. Measuring the projected linear sizes of these sources provides an indication of their orientation. Steep-spectrum isotropic radio emission allows for orientation-free sample selection at low radio frequencies. We use a new radio survey of the Bo\"otes field at 150 MHz made with the Low Frequency Array (LOFAR) to select a sample of radio sources. We identify 44 radio galaxies and 16 quasars with powers $P>10^{25.5}$ W Hz$^{-1}$ at 150 MHz using cross-matched multi-wavelength information from the AGN and Galaxy Evolution Survey (AGES), which provides spectroscopic redshifts. We find that LOFAR-detected radio sources with steep spectra have projected linear sizes that are on average 4.4$\pm$1.4 larger than those with flat spectra. The projected linear sizes of radio galaxies are on average 3.1$\pm$1.0 larger than those of quasars (2.0$\pm$0.3 after correcting for redshift evolution). Combining these results with three previous surveys, we find that the projected linear sizes of radio galaxies and quasars depend on redshift but not on power. The projected linear size ratio does not correlate with either parameter. The LOFAR data is consistent within the uncertainties with theoretical predictions of the correlation between the quasar fraction and linear size ratio, based on an orientation-based unification scheme.
We present a study of the low-frequency radio properties of star forming (SF) galaxies and active galactic nuclei (AGN) up to redshift $z=2.5$. The new spectral window probed by the Low Frequency Array (LOFAR) allows us to reconstruct the radio continuum emission from 150 MHz to 1.4 GHz to an unprecedented depth for a radio-selected sample of $1542$ galaxies in $\sim 7~ \rm{deg}^2$ of the LOFAR Bo\"otes field. Using the extensive multi-wavelength dataset available in Bo\"otes and detailed modelling of the FIR to UV spectral energy distribution (SED), we are able to separate the star-formation (N=758) and the AGN (N=784) dominated populations. We study the shape of the radio SEDs and their evolution across cosmic time and find significant differences in the spectral curvature between the SF galaxy and AGN populations. While the radio spectra of SF galaxies exhibit a weak but statistically significant flattening, AGN SEDs show a clear trend to become steeper towards lower frequencies. No evolution of the spectral curvature as a function of redshift is found for SF galaxies or AGN. We investigate the redshift evolution of the infrared-radio correlation (IRC) for SF galaxies and find that the ratio of total infrared to 1.4 GHz radio luminosities decreases with increasing redshift: $ q_{\rm 1.4GHz} = (2.45 \pm 0.04) \times (1+z)^{-0.15 \pm 0.03} $. Similarly, $q_{\rm 150MHz}$ shows a redshift evolution following $ q_{\rm 150GHz} = (1.72 \pm 0.04) \times (1+z)^{-0.22 \pm 0.05}$. Calibration of the 150 MHz radio luminosity as a star formation rate tracer suggests that a single power-law extrapolation from $q_{\rm 1.4GHz}$ is not an accurate approximation at all redshifts.
Does the environment of a galaxy directly influence the quenching history of a galaxy? Here we investigate the detailed morphological structures and star formation histories of a sample of SDSS group galaxies with both classifications from Galaxy Zoo 2 and NUV detections in GALEX. We use the optical and NUV colours to infer the quenching time and rate describing a simple exponentially declining SFH for each galaxy, along with a control sample of field galaxies. We find that the time since quenching and the rate of quenching do not correlate with the relative velocity of a satellite but are correlated with the group potential. This quenching occurs within an average quenching timescale of $\sim2.5~\rm{Gyr}$ from star forming to complete quiescence, during an average infall time (from $\sim 10R_{200}$ to $0.01R_{200}$) of $\sim 2.6~\rm{Gyr}$. Our results suggest that the environment does play a direct role in galaxy quenching through quenching mechanisms which are correlated with the group potential, such as harassment, interactions or starvation. Environmental quenching mechanisms which are correlated with satellite velocity, such as ram pressure stripping, are not the main cause of quenching in the group environment. We find that no single mechanism dominates over another, except in the most extreme environments or masses. Instead an interplay of mergers, mass & morphological quenching and environment driven quenching mechanisms dependent on the group potential drive galaxy evolution in groups.
We study the orientations of satellite galaxies in redMaPPer clusters constructed from the Sloan Digital Sky Survey at $0.1<z<0.35$ to determine whether there is any preferential tendency for satellites to point radially toward cluster centers. We analyze the satellite alignment (SA) signal based on three shape measurement methods (re-Gaussianization, de Vaucouleurs, and isophotal shapes), which trace galaxy light profiles at different radii. While no net SA signal is detected using re-Gaussianization shapes across the entire sample, the observed SA signal reaches a statistically significant level when using a subsample of satellites with higher luminosity. We detect the strongest SA signals using isophotal shapes, followed by de Vaucouleurs shapes, and investigate the impact of noise, systematics, and real physical effects such as isophotal twisting in the comparison between the results based on different shape measurement methods. After studying the correlation of the SA signal with a total of 17 galaxy and cluster properties, we find that the measured SA signal is strongest for satellites with the following characteristics: higher luminosity, smaller distance to the cluster center, rounder in shape, higher bulge fraction in the light profile, distributed preferentially along the major axis directions of their centrals, and residing in clusters with less luminous centrals. Finally, we provide physical explanations for the identified dependences, and discuss the connection to theories of SA.
We attempt to constrain the kinematics of the thin and thick discs using the Besancon population synthesis model together with RAVE DR4 and Gaia first data release (TGAS). The RAVE fields were simulated applying a detailed target selection function and the kinematics was computed using velocity ellipsoids depending on age in order to study the secular evolution. We accounted for the asymmetric drift computed from fitting a St\"ackel potential to orbits. Model parameters such as velocity dispersions, mean motions, velocity gradients, were adjusted using an ABC-MCMC method. We made use of the metallicity to enhanced the separation between thin and thick discs. We show that this model is able to reproduce the kinematics of the local discs in great detail. The disc follows the expected secular evolution, in very good agreement with previous studies of the thin disc. The new asymmetric drift formula, fitted to the St\"ackel potential described in Bienaym\'e et al (2015), fairly well reproduces the velocity distribution in a wide solar neighbourhood. The U and W components of the Solar motion determined with this method agree well with previous studies. However we find a smaller V component than previously thought, essentially due to the inclusion of the variation of the asymmetric drift with distance to the plane. The thick disc is represented by a long period of formation (at least 2 Gyr) during which it is shown that the mean velocity increases with time, while the scale height and scale length decrease, very consistently with a collapse phase with conservation of angular momentum. This new Galactic dynamical model is able to reproduce the observed velocities in a wide solar neighborhood at the level of quality of the TGAS-RAVE sample, allowing to constrain the thin and thick disc dynamical evolution, as well as determining the Solar motion.
Motivated by observational and theoretical work which both suggest very small scale ($\lesssim 1\,$pc) structure in the circum-galactic medium of galaxies and in other environments, we study Lyman-$\alpha$ (Ly$\alpha$) radiative transfer in an extremely clumpy medium with many "clouds" of neutral gas along the line of sight. While previous studies have typically considered radiative transfer through sightlines intercepting $\lesssim 10$ clumps, we explore the limit of a very large number of clumps per sightline (up to $f_{\mathrm{c}} \sim 1000$). Our main finding is that, for covering factors greater than some critical threshold, a multiphase medium behaves similar to a homogeneous medium in terms of the emergent Ly$\alpha$ spectrum. The value of this threshold depends on both the clump column density and on the movement of the clumps. We estimate this threshold analytically and compare our findings to radiative transfer simulations with a range of covering factors, clump column densities, radii, and motions. Our results suggest that (i) the success in fitting observed Ly$\alpha$ spectra using homogeneous "shell models" (and the corresponding failure of multiphase models) hints towards the presence of very small-scale structure in neutral gas, in agreement within a number of other observations; and (ii) the recurrent problems of reproducing realistic line profiles from hydrodynamical simulations may be due to their inability to resolve small-scale structure, which causes simulations to underestimate the effective covering factor of neutral gas clouds.
We investigate the relation between kinematic morphology, intrinsic colour and stellar mass of galaxies in the EAGLE cosmological hydrodynamical simulation. We calculate the intrinsic u-r colours and measure the fraction of kinetic energy invested in ordered corotation of 3562 galaxies at z=0 with stellar masses larger than $10^{10}M_{\odot}$. We perform a visual inspection of gri-composite images and find that our kinematic morphology correlates strongly with visual morphology. EAGLE produces a galaxy population for which morphology is tightly correlated with the location in the colour- mass diagram, with the red sequence mostly populated by elliptical galaxies and the blue cloud by disc galaxies. Satellite galaxies are more likely to be on the red sequence than centrals, and for satellites the red sequence is morphologically more diverse. These results show that the connection between mass, intrinsic colour and morphology arises from galaxy formation models that reproduce the observed galaxy mass function and sizes.
Determining the velocity distribution of halo stars is essential for estimating the mass of the Milky Way and for inferring its formation history. Since the stellar halo is a dynamically hot system, the velocity distribution of halo stars is well described by the 3-dimensional velocity dispersions $(\sigma_r, \sigma_\theta, \sigma_\phi)$, or by the velocity anisotropy parameter $\beta=1-(\sigma_\theta^2+\sigma_\phi^2)/(2\sigma_r^2)$. Direct measurements of $(\sigma_r, \sigma_\theta, \sigma_\phi)$ consistently suggest $\beta =0.5$-$0.7$ for nearby halo stars. In contrast, the value of $\beta$ at large Galactocentric radius $r$ is still controversial, since reliable proper motion data are available for only a handful of stars. In the last decade, several authors have tried to estimate $\beta$ for distant halo stars by fitting the observed line-of-sight velocities at each radius with simple velocity distribution models (local fitting methods). Some results of local fitting methods imply $\beta<0$ at $r \gtrsim 20 \;\rm{kpc}$, which is inconsistent with recent predictions from cosmological simulations. Here we perform mock-catalogue analyses to show that the estimates of $\beta$ based on local fitting methods are reliable only at $r \leq 15 \;\rm{kpc}$ with the current sample size ($\sim10^3$ stars at a given radius). As $r$ increases, the line-of-sight velocity (corrected for the Solar reflex motion) becomes increasingly closer to the Galactocentric radial velocity, so that it becomes increasingly more difficult to estimate tangential velocity dispersion $(\sigma_\theta, \sigma_\phi)$ from line-of-sight velocity distribution. Our results suggest that the forthcoming Gaia data will be crucial for understanding the velocity distribution of halo stars at $r \geq 20\;\rm{kpc}$.
We present an overview of the first data release (DR1) and first-look science from the Green Bank Ammonia Survey (GAS). GAS is a Large Program at the Green Bank Telescope to map all Gould Belt star-forming regions with $A_V \gtrsim 7$ mag visible from the northern hemisphere in emission from NH$_3$ and other key molecular tracers. This first release includes the data for four regions in Gould Belt clouds: B18 in Taurus, NGC 1333 in Perseus, L1688 in Ophiuchus, and Orion A North in Orion. We compare the NH$_3$ emission to dust continuum emission from Herschel, and find that the two tracers correspond closely. NH$_3$ is present in over 60\% of lines-of-sight with $A_V \gtrsim 7$ mag in three of the four DR1 regions, in agreement with expectations from previous observations. The sole exception is B18, where NH$_3$ is detected toward ~ 40\% of lines-of-sight with $A_V \gtrsim 7$ mag. Moreover, we find that the NH$_3$ emission is generally extended beyond the typical 0.1 pc length scales of dense cores. We produce maps of the gas kinematics, temperature, and NH$_3$ column densities through forward modeling of the hyperfine structure of the NH$_3$ (1,1) and (2,2) lines. We show that the NH$_3$ velocity dispersion, ${\sigma}_v$, and gas kinetic temperature, $T_K$, vary systematically between the regions included in this release, with an increase in both the mean value and spread of ${\sigma}_v$ and $T_K$ with increasing star formation activity. The data presented in this paper are publicly available.
We present a study of the hierarchical clustering of the young stellar clusters in six local (3--15 Mpc) star-forming galaxies using Hubble Space Telescope broad band WFC3/UVIS UV and optical images from the Treasury Program LEGUS (Legacy ExtraGalactic UV Survey). We have identified 3685 likely clusters and associations, each visually classified by their morphology, and we use the angular two-point correlation function to study the clustering of these stellar systems. We find that the spatial distribution of the young clusters and associations are clustered with respect to each other, forming large, unbound hierarchical star-forming complexes that are in general very young. The strength of the clustering decreases with increasing age of the star clusters and stellar associations, becoming more homogeneously distributed after ~40--60 Myr and on scales larger than a few hundred parsecs. In all galaxies, the associations exhibit a global behavior that is distinct and more strongly correlated from compact clusters. Thus, populations of clusters are more evolved than associations in terms of their spatial distribution, traveling significantly from their birth site within a few tens of Myr whereas associations show evidence of disruption occurring very quickly after their formation. The clustering of the stellar systems resembles that of a turbulent interstellar medium that drives the star formation process, correlating the components in unbound star-forming complexes in a hierarchical manner, dispersing shortly after formation, suggestive of a single, continuous mode of star formation across all galaxies.
We present optical integral field spectroscopy $-$ obtained with the Gemini Multi-Object Spectrograph $-$ of the inner $4.0 \times 5.8$ kpc$^2$ of the narrow line radio galaxy 3C 33 at a spatial resolution of 0.58 kpc. The gas emission shows three brightest structures: a strong knot of nuclear emission and two other knots at $\approx 1.4$ kpc south-west and north-east of the nucleus along the ionization axis. We detect two kinematic components in the emission lines profiles, with a "broader component" (with velocity dispersion $\sigma > 150$ km s$^{-1}$) being dominant within a $\sim$ 1 kpc wide strip ("the nuclear strip") running from the south-east to the north-west, perpendicular to the radio jet, and a narrower component ($\sigma < 100$ km s$^{-1}$) dominating elsewhere. Centroid velocity maps reveal a rotation pattern with velocity amplitudes reaching $\sim \pm 350$ km s$^{-1}$ in the region dominated by the narrow component, while residual blueshifts and redshifts relative to rotation are observed in the nuclear strip, where we also observe the highest values of the [N II]/H{\alpha}, [S II]/H{\alpha} and [O I]/H{\alpha} line ratios, and an increase of the gas temperature ($\sim 18000$ K), velocity dispersion and electron density ($\sim 500$ cm$^{-3}$). We interpret these residuals and increased line ratios as due to a lateral expansion of the ambient gas in the nuclear strip due to shocks produced by the passage of the radio jet. The effect of this expansion in the surrounding medium is very small, as its estimated kinetic power represents only $2.6 - 3.0 \times 10^{-5}$ of the AGN bolometric luminosity. A possible signature of inflow is revealed by an increase in the [O I]/H{\alpha} ratio values and velocity dispersions in the shape of two spiral arms extending to 2.3 kpc north-east and south-west from the nucleus.
During the Space Telescope and Optical Reverberation Mapping Project (STORM) observations of NGC 5548, the continuum and emission-line variability became de-correlated during the second half of the 6-month long observing campaign. Here we present Swift and Chandra X-ray spectra of NGC 5548 obtained as a part of the campaign. The Swift spectra show that excess flux (relative to a power-law continuum) in the soft X-ray band appears before the start of the anomalous emission-line behavior, peaks during the period of the anomaly, and then declines. This is a model-independent result suggesting that the soft excess is related to the anomaly. We divide the Swift data into on- and off-anomaly spectra to characterize the soft excess via spectral fitting. The cause of the spectral differences is likely due to a change in the intrinsic spectrum rather than being due to variable obscuration or partial covering. The Chandra spectra have lower signal-to-noise ratios, but are consistent with Swift data. Our preferred model of the soft excess is emission from an optically thick, warm Comptonizing corona, the effective optical depth of which increases during the anomaly. This model simultaneously explains all the three observations: the UV emission line flux decrease, the soft-excess increase, and the emission line anomaly.
We have modified the iterative procedure introduced by Lin et al. (2016), to systematically combine the submm images taken from ground based (e.g., CSO, JCMT, APEX) and space (e.g., Herschel, Planck) telescopes. We applied the updated procedure to observations of three well studied Infrared Dark Clouds (IRDCs): G11.11-0.12, G14.225-0.506 and G28.34+0.06, and then performed single-component, modified black-body fits to derive $\sim$10$"$ resolution dust temperature and column density maps. The derived column density maps show that these three IRDCs exhibit complex filamentary structures embedding with rich clumps/cores. We compared the column density probability distribution functions (N-PDFs) and two-point correlation (2PT) functions of the column density field between these IRDCs with several OB cluster-forming regions. Based on the observed correlation and measurements, and complementary hydrodynamical simulations for a 10$^{4}$ $\rm M_{\odot}$ molecular cloud, we hypothesize that cloud evolution can be better characterized by the evolution of the (column) density distribution function and the relative power of dense structures as a function of spatial scales, rather than merely based on the presence of star-forming activity. Based on the small analyzed sample, we propose four evolutionary stages, namely: {\it cloud integration, stellar assembly, cloud pre-dispersal and dispersed-cloud.} The initial {\it cloud integration} stage and the final {\it dispersed cloud} stage may be distinguished from the two intermediate stages by a steeper than $-$4 power-law index of the N-PDF. The {\it cloud integration} stage and the subsequent {\it stellar assembly} stage are further distinguished from each other by the larger luminosity-to-mass ratio ($>$40 $\rm L_{\odot}/M_{\odot}$) of the latter.
We present 3D simulations of the formation process of two molecular clouds (MCs) within their larger-scale galactic environment. Using adaptive mesh refinement, we model the two MCs within the SILCC project with an unprecedented resolution of 0.06 pc combined with a chemical network for the formation of H$_2$ and CO including (self-) shielding and important thermal processes. The MCs form within a few Myr with mass growth rates of up to 10$^{-2}$ M$_\rm{sun}$ yr$^{-1}$ and final masses of $\sim$ 50000 M$_\rm{sun}$. We show that the usage of different definitions for MCs by thresholds in density, H$_2$ or CO mass fraction significantly change the inferred cloud properties. While CO traces well the evolution of dense gas with $n \geq$ 300 cm$^{-3}$, H$_2$ is also found in gas with lower number density ($n \lesssim$ 30 cm$^{-3}$) due to turbulent mixing. The CO-to-H$_2$ ratio increases within the first 2 Myr reaching a value of $\sim$ 1.8 $\times$ 10$^{-4}$ at later stages. The $X_\rm{CO}$ factor, however, is rather time-independent with values of 1 - 4 $\times$ 10$^{20}$ cm$^{-2}$ (K km s$^{-1}$)$^{-1}$. We show that a spatial resolution of $\sim$ 0.1 pc is required to accurately model the chemical, dynamical, and structural evolution of MCs. At a coarser resolution the mass, velocity dispersion, and chemical abundances of the clouds are underestimated. Furthermore, we show that the progressive increase of resolution has to occur over a time of 1 - 1.5 Myr. This ensures that the maximum refinement level is reached within the free-fall time of the densest structures and avoids the spurious formation large-scale, rotating objects by unresolved turbulent flows. In addition, the accelerated formation of chemical species in dense, turbulent environments is captured properly. Finally, we demonstrate that $\gtrsim$ 200 time steps should be spent on each refinement level to avoid grid artefacts.
We present clustering properties from 579,492 Lyman break galaxies (LBGs) at $z\sim4-6$ over the 100 deg$^2$ sky (corresponding to a 1.4 Gpc$^3$ volume) identified in early data of the Hyper Suprime-Cam (HSC) Subaru strategic program survey. We derive angular correlation functions (ACFs) of the HSC LBGs with unprecedentedly high statistical accuracies at $z\sim4-6$, and compare them with the halo occupation distribution (HOD) models. We clearly identify significant ACF excesses in $10"<\theta<90"$, the transition scale between 1- and 2-halo terms, suggestive of the existence of the non-linear halo bias effect. Combining the HOD models and previous clustering measurements of faint LBGs at $z\sim4-7$, we investigate dark-matter halo mass ($M_\mathrm{h}$) of the $z\sim4-7$ LBGs and its correlation with various physical properties including the star-formation rate (SFR), the stellar-to-halo mass ratio (SHMR), and the dark-matter mass accretion rate ($\dot{M}_\mathrm{h}$) over a wide-mass range of $M_\mathrm{h}/M_\odot=4\times10^{10}-4\times10^{12}$. We find that the SHMR increases from $z\sim4$ to $7$ by a factor of $\sim4$ at $M_\mathrm{h}\simeq1\times10^{11}\ M_\odot$, while the SHMR shows no strong evolution in the similar redshift range at $M_\mathrm{h}\simeq1\times10^{12}\ M_\odot$. Interestingly, we identify a tight relation of $SFR/\dot{M}_\mathrm{h}-M_\mathrm{h}$ showing no significant evolution beyond 0.15 dex in this wide-mass range over $z\sim4-7$. This weak evolution suggests that the $SFR/\dot{M}_\mathrm{h}-M_\mathrm{h}$ relation is a fundamental relation in high-redshift galaxy formation whose star-formation activities are regulated by the dark-matter mass assembly.
The radiative transport of photons through arbitrary three-dimensional (3D) structures of dust is a challenging problem due to the anisotropic scattering of dust grains and strong coupling between different spatial regions. The radiative transfer problem in 3D is solved using Monte Carlo or Ray Tracing techniques as no full analytic solution exists for the true 3D structures. We provide the first 3D dust radiative transfer benchmark composed of a slab of dust with uniform density externally illuminated by a star. This simple 3D benchmark is explicitly formulated to provide tests of the different components of the radiative transfer problem including dust absorption, scattering, and emission. This benchmark includes models with a range of dust optical depths fully probing cases that are optically thin at all wavelengths to optically thick at most wavelengths. This benchmark includes solutions for the full dust emission including single photon (stochastic) heating as well as two simplifying approximations: One where all grains are considered in equilibrium with the radiation field and one where the emission is from a single effective grain with size-distribution-averaged properties. A total of six Monte Carlo codes and one Ray Tracing code provide solutions to this benchmark. Comparison of the results revealed that the global SEDs are consistent on average to a few percent for all but the scattered stellar flux at very high optical depths. The image results are consistent within 10%, again except for the stellar scattered flux at very high optical depths. The lack of agreement between different codes of the scattered flux at high optical depths is quantified for the first time. We provide the first 3D dust radiative transfer benchmark and validate the accuracy of this benchmark through comparisons between multiple independent codes and detailed convergence tests.
We reconsider the lower bound on the mass of a fermionic dark matter (DM) candidate resulting from the existence of known small Dwarf Spheroidal galaxies, in the hypothesis that their DM halo is constituted by degenerate fermions, with phase-space density limited by the Pauli exclusion principle. By relaxing the common assumption that the DM halo scale radius is tied to that of the luminous stellar component and by marginalizing on the unknown stellar velocity dispersion anisotropy, we prove that observations lead to rather weak constraints on the DM mass, that could be as low as tens of eV. In this scenario, however, the DM halos would be quite large and massive, so that a bound stems from the requirement that the time of orbital decay due to dynamical friction in the hosting Milky Way DM halo is longer than their lifetime. The smallest and nearest satellites Segue I and Willman I lead to a final lower bound of $m\gtrsim100$ eV, still weaker than previous estimates but robust and independent on the model of DM formation and decoupling. We thus show that phase space constraints do not rule out the possibility of sub-keV fermionic DM.
We present three-dimensional simulations of a protoplanetary disk subject to the effect of a nearby (0.3pc distant) supernova, using a time-dependent flow from a one dimensional numerical model of the supernova remnant (SNR), in addition to constant peak ram pressure simulations. Simulations are performed for a variety of disk masses and inclination angles. We find disk mass-loss rates that are typically 1e-7 to 1e-6 Msol/yr (but peak near 1e-5 Msol/yr during the "instantaneous" stripping phase) and are sustained for around 200 yr. Inclination angle has little effect on the mass loss unless the disk is close to edge-on. Inclined disks also strip asymmetrically with the trailing edge ablating more easily. Since the interaction lasts less than one outer rotation period, there is not enough time for the disk to restore its symmetry, leaving the disk asymmetrical after the flow has passed. Of the low-mass disks considered, only the edge-on disk is able to survive interaction with the SNR (with 50% of its initial mass remaining). At the end of the simulations, disks that survive contain fractional masses of SN material up to 5e-6. This is too low to explain the abundance of short-lived radionuclides in the early solar system, but a larger disk and the inclusion of radiative cooling might allow the disk to capture a higher fraction of SN material.
Different combinations of input parameters to filament identification algorithms, such as Disperse and FilFinder, produce numerous different output skeletons. The skeletons are a one pixel wide representation of the filamentary structure in the original input image. However, these output skeletons may not necessarily be a good representation of that structure. Furthermore, a given skeleton may not be as good a representation as another. Previously there has been no mathematical `goodness-of-fit' measure to compare output skeletons to the input image. Thus far this has been assessed visually, introducing visual bias. We propose the application of the mean structural similarity index (MSSIM) as a mathematical goodness-of-fit measure. We describe the use of the MSSIM to find the output skeletons most mathematically similar to the original input image (the optimum, or `best', skeletons) for a given algorithm, and independently of the algorithm. This measure makes possible systematic parameter studies, aimed at finding the subset of input parameter values returning optimum skeletons. It can also be applied to the output of non-skeleton based filament identification algorithms, such as the Hessian matrix method. The MSSIM removes the need to visually examine thousands of output skeletons, and eliminates the visual bias, subjectivity, and limited reproducibility inherent in that process, representing a major improvement on existing techniques. Importantly, it also allows further automation in the post-processing of output skeletons, which is crucial in this era of `big data'.
The profiles of the broad emission lines of active galactic nuclei (AGNs) and the time delays in their response to changes in the ionizing continuum ("lags") give information about the structure and kinematics of the inner regions of AGNs. Line profiles are also our main way of estimating the masses of the supermassive black holes (SMBHs). However, the profiles often show ill-understood, asymmetric structure and velocity-dependent lags vary with time. Here we show that partial obscuration of the BLR by outflowing, compact, dusty clumps produces asymmetries and velocity-dependent lags similar to those observed. Our model explains previously inexplicable changes in the ratios of the hydrogen lines with time and velocity, the lack of correlation of changes in line profiles with variability of the central engine, the velocity dependence of lags, and the change of lags with time. We propose that changes on timescales longer than the light-crossing time do not come from dynamical changes in the BLR, but are a natural result of the effect of outflowing dusty clumps driven by radiation pressure acting on the dust. The effects of this dust complicate the study of the structure and kinematics of the BLR and the search for sub-parsec SMBH binaries. Partial obscuration of the accretion disc can also explain the puzzling flux gradients implied by microlensing.
An infinite family of exact solutions of the electrovacuum Einstein-Maxwell equations is presented. The family is static, axially symmetric and describe thin disks composed by electrically polarized material in a conformastatic spacetime. The form of the conformastatic metric allows us to write down the metric functions and the electromagnetic potentials in terms of a solution of the Laplace equation. We find a general expression for the surface energy density of the disk, the pressure, the polarization vector, the electromagnetic fields and the velocity rotation for circular orbits. As an example, we present the first model of the family and show the behavior of the different physical variables.
Links to: arXiv, form interface, find, astro-ph, recent, 1704, contact, help (Access key information)
Radio relics in galaxy clusters are believed to trace merger shock fronts. If cosmological structure formation determines the luminosity, size and shape distributions of radio relics then merger shocks need to be lighted up in a homogeneous way. We investigate if a mock relic sample, obtained from zoomed galaxy cluster simulations, is able to match the properties of relics measured in the NRAO VLA Sky Survey (NVSS). We compile a list of all radio relics known to date and homogeneously measure their parameters in all NVSS images and apply the same procedure to relics in our simulations. Number counts in the mock relic sample increase more steeply towards lower relic flux densities, suggesting an incompleteness of NVSS in this regime. Overall, we find that NVSS and mock samples show similar properties. However, large simulated relics tend to be somewhat smaller and closer to the cluster centre than observed ones. Besides this, the mock sample reproduces very well-known correlations for radio relics, in particular those relating the radio luminosity with the largest linear size and the X-ray luminosity. We show that these correlations are largely governed by the sensitivity of the NVSS observations. Mock relics show a similar orientation with respect to the direction to the cluster centre as the NVSS sample. Moreover, we find that their maximum radio luminosity roughly correlates with cluster mass, although displaying a large scatter. The overall good agreement between NVSS and the mock sample suggests that properties of radio relics are indeed governed by merger shock fronts, emitting in a homogeneous fashion. Our study demonstrates that the combination of mock observations and data from upcoming radio surveys will allow to shed light on both the origin of radio relics and the nature of the intracluster medium.
We present ALMA $^{12}$CO (J=1-0, 3-2 and 6-5), $^{13}$CO (J=1-0) and C$^{18}$O (J=1-0) observations of the local Ultra Luminous Infrared Galaxy, IRAS 13120-5453 (dubbed "The Yo-yo"). The morphologies of the three isotopic species differ, where $^{13}$CO shows a hole in emission towards the center. We measure integrated brightness temperature line ratios of $^{12}$CO/$^{13}$CO $\geq$ 60 (exceeding 200) and $^{13}$CO/C$^{18}$O $\leq$ 1 in the central region. Assuming optical thin emission, C$^{18}$O is more abundant than $^{13}$CO in several regions. The abundances within the central 500 pc are consistent with enrichment of the ISM via a young starburst ($<$7Myr), a top-heavy initial mass function or a combination of both.
We report the discovery of a diffuse stellar cloud with an angular extent $\gtrsim30^{\prime\prime}$, which we term "Sumo Puff", in data from the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP). While we do not have a redshift for this object, it is in close angular proximity to a post-merger galaxy at redshift $z=0.0431$ and is projected within a few virial radii (assuming similar redshifts) of two other ${\sim}L_\star$ galaxies, which we use to bracket a potential redshift range of $0.0055 < z < 0.0431$. The object's light distribution is flat, as characterized by a low Sersic index ($n\sim0.3$). It has a low central $g$-band surface brightness of ${\sim}26.4$ mag arcsec$^{-2}$, large effective radius of ${\sim}13^{\prime\prime}$ (${\sim}11$ kpc at $z=0.0431$ and ${\sim}1.5$ kpc at $z=0.0055$), and an elongated morphology ($b/a\sim0.4$). Its red color ($g-i\sim1$) is consistent with a passively evolving stellar population and similar to the nearby post-merger galaxy, and we may see tidal material connecting Sumo Puff with this galaxy. We offer two possible interpretations for the nature of this object: (1) it is an extreme, galaxy-size tidal feature associated with a recent merger event, or (2) it is a foreground dwarf galaxy with properties consistent with a quenched, disturbed ultra-diffuse galaxy. We present a qualitative comparison with simulations that demonstrates the feasibility of forming a structure similar to this object in a merger event. Follow-up spectroscopy and/or deeper imaging to confirm the presence of the bridge of tidal material will be necessary to reveal the true nature of this object.
We use state-of-the-art chemical models to track the cosmic evolution of the CNO isotopes in the interstellar medium (ISM) of galaxies, yielding powerful constraints on their stellar initial mass function (IMF). We re-assess the relative roles of massive stars, asymptotic giant branch (AGB) stars and novae in the production of rare isotopes such as 13C, 15N, 17O and 18O, along with 12C, 14N and 16O. The CNO isotope yields of super-AGB stars, novae and fast-rotating massive stars are included. Having reproduced the available isotope enrichment data in the solar neighbourhood, and across the Galaxy, and having assessed the sensitivity of our models to the remaining uncertainties, e.g. nova yields and star-formation history, we show that we can meaningfully constrain the stellar IMF in galaxies using C, O and N isotope abundance ratios. In starburst galaxies, where data for multiple isotopologue lines are available, we find compelling new evidence for a top-heavy stellar IMF, with profound implications for their star-formation rates and efficiencies, perhaps also their stellar masses. Neither chemical fractionation nor selective photodissociation can significantly perturb globally-averaged isotopologue abundance ratios away from the corresponding isotope ones, as both these processes will typically affect only small mass fractions of molecular clouds in galaxies. Thus the Atacama Large Millimetre Array now stands ready to probe the stellar IMF, and even the ages of specific starburst events in star-forming galaxies across cosmic time unaffected by the dust obscuration effects that plague optical/near-infrared studies.
Several authors have claimed that the less luminous active galactic nuclei (AGN) are not capable of sustaining the dusty torus structure. Thus, a gradual re-sizing of the torus is expected when the AGN luminosity decreases. Our aim is to confront mid-infrared observations of local AGN of different luminosities with this scenario. We decomposed about ~100 IRS/Spitzer spectra of LLAGN and powerful Seyferts in order to decontaminate the torus component from other contributors. We have used the affinity propagation (AP) method to cluster the data into five groups within the sample according to torus contribution to the 5-15 um range (Ctorus) and bolometric luminosity. The AP groups show a progressively higher torus contribution and an increase of the bolometric luminosity, from Group 1 (Ctorus~ 0% and logLbol ~ 41) and up to Group 5 (Ctorus ~80% and log(Lbol) ~44). We have fitted the average spectra of each of the AP groups to clumpy models. The torus is no longer present in Group 1, supporting the disappearance at low-luminosities. We were able to fit the average spectra for the torus component in Groups 3 (Ctorus~ 40% and log(Lbol)~ 42.6), 4 (Ctorus~ 60% and log(Lbol)~ 43.7), and 5 to Clumpy torus models. We did not find a good fitting to Clumpy torus models for Group 2 (Ctorus~ 18% and log(Lbol)~ 42). This might suggest a different configuration and/or composition of the clouds for Group 2, which is consistent with a different gas content seen in Groups 1, 2, and 3, according to the detections of H2 molecular lines. Groups 3, 4, and 5 show a trend to decrease of the width of the torus (which yields to a likely decrease of the geometrical covering factor), although we cannot confirm it with the present data. Finally, Groups 3, 4, and 5 show an increase on the outer radius of the torus for higher luminosities, consistent with a re-sizing of the torus according to the AGN luminosity.
The spectra of 413 star-forming (or HII) regions in M33 (NGC 598) were observed by using the multifiber spectrograph of Hectospec at the 6.5-m Multiple Mirror Telescope (MMT). By using this homogeneous spectra sample, we measured the intensities of emission lines and some physical parameters, such as electron temperatures, electron densities, and metallicities. Oxygen abundances were derived via the direct method (when available) and two empirical strong-line methods, namely, O3N2 and N2. In the high-metallicity end, oxygen abundances derived from O3N2 calibration were higher than those derived from N2 index, indicating an inconsistency between O3N2 and N2 calibrations. We presented a detailed analysis of the spatial distribution of gas-phase oxygen abundances in M33 and confirmed the existence of the axisymmetric global metallicity distribution widely assumed in literature. Local variations were also observed and subsequently associated with spiral structures to provide evidence of radial migration driven by arms. Our O/H gradient fitted out to 1.1 $R_{25}$ resulted in slopes of $-0.17\pm0.03$, $-0.19\pm0.01$, and $-0.16\pm0.17$ dex $R_{25}^{-1}$ utilizing abundances from O3N2, N2 diagnostics, and direct method, respectively.
We study the global consequences in the halos of spiral galaxies of the steady, axially symmetric, mean field dynamo. We use the classical theory but add the possibility of using the velocity field components as parameters in addition to the helicity and diffusivity. The analysis is based on the simplest version of the theory and uses scale-invariant solutions. The velocity field (subject to restrictions) is a scale invariant field in a `pattern' frame, in place of a full dynamical theory. The `pattern frame' of reference may either be the systemic frame or some rigidly rotating spiral pattern frame. One type of solution for the magnetic field yields off-axis, spirally wound, magnetic field lines. These predict sign changes in the Faraday screen rotation measure in every quadrant of the halo of an edge-on galaxy. Such rotation measure oscillations have been observed in the CHANG-ES survey.
We seek the conditions for a {\it steady} mean field galactic dynamo. The parameter set is reduced to those appearing in the $\alpha^2$ and $\alpha/\omega$ dynamo, namely velocity amplitudes, and the ratio of sub-scale helicity to diffusivity. The parameters can be allowed to vary on conical spirals. We analyze the mean field dynamo equations in terms of scale invariant logarithmic spiral modes and special exact solutions. Compatible scale invariant gravitational spiral arms are introduced and illustrated in an appendix, but the detailed dynamical interaction with the magnetic field is left for another work. As a result of planar magnetic spirals `lifting' into the halo, multiple sign changes in average rotation measures forming a regular pattern on each side of the galactic minor axis, are predicted. Such changes have recently been detected in the CHANG-ES survey.
Galactic star formation scaling relations show increased scatter from kpc to
sub-kpc scales. Investigating this scatter may hold important clues to how the
star formation process evolves in time and space. Here, we combine different
molecular gas tracers, different star formation indicators probing distinct
populations of massive stars, and knowledge on the evolutionary state of each
star forming region to derive star formation properties of $\sim$150 star
forming complexes over the face of the Large Magellanic Cloud. We find that the
rate of massive star formation ramps up when stellar clusters emerge and boost
the formation of subsequent generations of massive stars. In addition, we
reveal that the star formation efficiency of individual GMCs declines with
increasing cloud gas mass ($M_\mathrm{cloud}$). This trend persists in Galactic
star forming regions, and implies higher molecular gas depletion times for
larger GMCs.
We compare the star formation efficiency per freefall time
($\epsilon_\mathrm{ff}$) with predictions from various widely-used analytical
star formation models. We show that while these models can produce large
dispersions in $\epsilon_\mathrm{ff}$ similar to observations, the origin of
the model-predicted scatter is inconsistent with observations. Moreover, all
models fail to reproduce the observed decline of $\epsilon_\mathrm{ff}$ with
increasing $M_\mathrm{cloud}$ in the LMC and the Milky Way. We conclude that
analytical star formation models idealizing global turbulence levels, cloud
densities, and assuming a stationary SFR are inconsistent with observations
from modern datasets tracing massive star formation on individual cloud scales.
Instead, we reiterate the importance of local stellar feedback in shaping the
properties of GMCs and setting their massive star formation rate.
The range of currently proposed active galactic nucleus (AGN) far-infrared templates results in uncertainties in retrieving host galaxy information from infrared observations and also undermines constraints on the outer part of the AGN torus. We discuss how to test and reconcile these templates. Physically, the fraction of the intrinsic AGN IR-processed luminosity compared with that from the central engine should be consistent with the dust-covering factor. In addition, besides reproducing the composite spectral energy distributions (SEDs) of quasars, a correct AGN IR template combined with an accurate library of star-forming galaxy templates should be able to reproduce the IR properties of the host galaxies, such as the luminosity-dependent SED shapes and aromatic feature strengths. We develop tests based on these expected behaviors and find that the shape of the AGN intrinsic far-IR emission drops off rapidly starting at $\sim20~\mu$m and can be matched by an Elvis et al. (1994)-like template with minor modification. Despite the variations in the near- to mid-IR bands, AGNs in quasars and Seyfert galaxies have remarkably similar intrinsic far-IR SEDs at $\lambda \sim 20$-$100~\mu$ m, suggesting similar emission character of the outermost region of the circumnuclear torus. The variations of the intrinsic AGN IR SEDs among the type-1 quasar population can be explained by the changing relative strengths of four major dust components with similar characteristic temperatures, and there is evidence for compact AGN-heated dusty structures at sub-kpc scales in the far-IR.
The Extragalactic Background Light (EBL) captures the total integrated emission from stars and galaxies throughout the cosmic history. The amplitude of the near-infrared EBL from space absolute photometry observations has been controversial and depends strongly on the modeling and subtraction of the Zodiacal light foreground. We report the first measurement of the diffuse background spectrum at 0.8-1.7 um from the CIBER experiment. The observations were obtained with an absolute spectrometer over two flights in multiple sky fields to enable the subtraction of Zodiacal light, stars, terrestrial emission, and diffuse Galactic light. After subtracting foregrounds and accounting for systematic errors, we find the nominal EBL brightness, assuming the Kelsall Zodiacal light model, is 42.7+11.9/-10.6 nW/m2/sr at 1.4 um. We also analyzed the data using the Wright Zodiacal light model, which results in a worse statistical fit to the data and an unphysical EBL, falling below the known background light from galaxies at <1.3 um. Using a model-independent analysis based on the minimum EBL brightness, we find an EBL brightness of 28.7+5.1/-3.3 nW/m2/sr at 1.4 um. While the derived EBL amplitude strongly depends on the Zodiacal light model, we find that we cannot fit the spectral data to Zodiacal light, Galactic emission, and EBL from solely integrated galactic light from galaxy counts. The results require a new diffuse component, such as an additional foreground or an excess EBL with a redder spectrum than that of Zodiacal light.
The ongoing surveys of galaxies and those for the next generation of telescopes will demand the execution of high-CPU consuming machine codes for recovering detailed star formation histories (SFHs) and hence age-metallicity relationships (AMRs). We present here an expeditive method which provides quick-look AMRs on the basis of representative ages and metallicities obtained from colour-magnitude diagram (CMD) analyses. We have tested its perfomance by generating synthetic CMDs for a wide variety of galaxy SFHs. The representative AMRs turn out to be reliable down to a magnitude limit with a photometric completeness factor higher than $\sim$ 85 per cent, and trace the chemical evolution history for any stellar population (represented by a mean age and an intrinsic age spread) with a total mass within ~ 40 per cent of the more massive stellar population in the galaxy.
The supermassive black holes in active galactic nuclei (AGN) and stellar-mass black holes in X-ray binaries (XRBs) are believed to work in a similar way. While XRBs evolve rapidly and several sources have undergone a few complete cycles from quiescence to an outburst and back, most AGN remain in the same state over periods of decades, due to their longer characteristic timescale proportional to their size. However, the study of the AGN spectral states is still possible with a large sample of sources. Multi-wavelength observations are needed for this purpose since the AGN thermal disc emission dominates in the ultraviolet energy range, while the up-scattered hot-corona emission is detected in X-rays. We compared simultaneous UV and X-ray measurements of AGN obtained by the XMM-Newton satellite. The non-thermal flux was constrained from the 2-12 keV X-ray luminosity, while the thermal disc component was estimated from the UV flux at 2900A. The hardness (ratio between the X-ray and UV plus X-ray luminosity) and the total luminosity were used to construct the AGN state diagrams. For sources with reliable mass measurements, the Eddington ratio was used instead of the total luminosity. The state diagrams show that the radio-loud sources have on average higher hardness, due to the lack of the thermal disc emission in the UV band, and have flatter intrinsic X-ray spectra. In contrast, the sources with high luminosity and low hardness are radio-quiet AGN. The hardness-Eddington ratio diagram reveals that the average radio-loudness is stronger for low-accreting sources, while it decreases when the accretion rate is close to the Eddington limit. Our results indicate that the general properties of AGN accretion states are similar to those of X-ray binaries. This suggests that the AGN radio dichotomy of radio-loud and radio-quiet sources can be explained by the evolution of the accretion states.
We utilize elemental-abundance information for Galactic red giant stars in five open clusters (NGC 7789, NGC 6819, M67, NGC 188, and NGC 6791) from the Apache Point Observatory Galactic Evolution Experiment (APOGEE) DR13 dataset to age-date the chemical evolution of the high- and low-$\alpha$ element sequences of the Milky Way. Key to this time-stamping is the cluster NGC 6791, whose stellar members have mean abundances that place it in the high-$\alpha$, high-[Fe/H] region of the [$\alpha$/Fe]-[Fe/H] plane. Based on the cluster's age ($\sim 8$ Gyr), Galactocentric radius, and height above the Galactic plane, as well as comparable chemistry reported for APOGEE stars in Baade's Window, we suggest that the two most likely origins for NGC 6791 are as an original part of the thick-disk, or as a former member of the Galactic bulge. Moreover, because NGC 6791 lies at the \textit{high metallicity end} ([Fe/H] $\sim 0.4$) of the high-$\alpha$ sequence, the age of NGC 6791 places a limit on the \textit{youngest age} of stars in the high-metallicity, high-$\alpha$ sequence for the cluster's parent population (i.e., either the bulge or the disk). In a similar way, we can also use the age and chemistry of NGC 188 to set a limit of $\sim 7$ Gyr on the \textit{oldest age} of the low-$\alpha$ sequence of the Milky Way. Therefore, NGC 6791 and NGC 188 are potentially a pair of star clusters that bracket both the timing and the duration of an important transition point in the chemical history of the Milky Way.
The formaldehyde MegaMaser emission has been mapped for the three host galaxies IC\,860. IRAS\,15107$+$0724, and Arp\,220. Elongated emission components are found at the nuclear centres of all galaxies with an extent ranging between 30 to 100 pc. These components are superposed on the peaks of the nuclear continuum. Additional isolated emission components are found superposed in the outskirts of the radio continuum structure. The brightness temperatures of the detected features ranges from 0.6 to 13.4 $\times 10^{4}$ K, which confirms their masering nature. The masering scenario is interpreted as amplification of the radio continuum by foreground molecular gas that is pumped by far-infrared radiation fields in these starburst environments of the host galaxies.
We investigate the dependence of black-hole accretion rate (BHAR) on host-galaxy star formation rate (SFR) and stellar mass ($M_*$) in the CANDELS/GOODS-South field in the redshift range of $0.5\leq z < 2.0$. Our sample consists of $\approx 18000$ galaxies, allowing us to probe galaxies with $0.1 \lesssim \mathrm{SFR} \lesssim 100\ M_\odot\ \mathrm{yr}^{-1}$ and/or $10^8 \lesssim M_* \lesssim 10^{11}\ M_{\odot}$. We use sample-mean BHAR to approximate long-term average BHAR. Our sample-mean BHARs are derived from the $Chandra$ Deep Field-South 7 Ms observations, while the SFRs and $M_*$ have been estimated by the CANDELS team through SED fitting. The average BHAR is correlated positively with both SFR and $M_*$, and the BHAR-SFR and BHAR-$M_*$ relations can both be described acceptably by linear models with a slope of unity. However, BHAR appears to be correlated more strongly with $M_*$ than SFR. This result indicates that $M_*$ is the primary host-galaxy property related to black-hole growth, and the apparent BHAR-SFR relation is largely a secondary effect due to the star-forming main sequence. Among our sources, massive galaxies ($M_* \gtrsim 10^{10} M_{\odot}$) have significantly higher BHAR/SFR ratios than less-massive galaxies, indicating the former have higher black-hole fueling efficiency and/or higher SMBH occupation fraction than the latter. Our results can naturally explain the observed proportionality between $M_{\rm BH}$ and $M_*$ for local giant ellipticals, and suggest their $M_{\rm BH}/M_*$ is higher than that of local star-forming galaxies. Among local star-forming galaxies, massive systems might have higher $M_{\rm BH}/M_*$ compared to dwarfs.
A key goal of the Stage IV dark energy experiments Euclid, LSST and WFIRST is to measure the growth of structure with cosmic time from weak lensing analysis over large regions of the sky. Weak lensing cosmology will be challenging: in addition to highly accurate galaxy shape measurements, statistically robust and accurate photometric redshift (photo-z) estimates for billions of faint galaxies will be needed in order to reconstruct the three-dimensional matter distribution. Here we present an overview of and initial results from the Complete Calibration of the Color-Redshift Relation (C3R2) survey, designed specifically to calibrate the empirical galaxy color-redshift relation to the Euclid depth. These redshifts will also be important for the calibrations of LSST and WFIRST. The C3R2 survey is obtaining multiplexed observations with Keck (DEIMOS, LRIS, and MOSFIRE), the Gran Telescopio Canarias (GTC; OSIRIS), and the Very Large Telescope (VLT; FORS2 and KMOS) of a targeted sample of galaxies most important for the redshift calibration. We focus spectroscopic efforts on under-sampled regions of galaxy color space identified in previous work in order to minimize the number of spectroscopic redshifts needed to map the color-redshift relation to the required accuracy. Here we present the C3R2 survey strategy and initial results, including the 1283 high confidence redshifts obtained in the 2016A semester and released as Data Release 1.
We briefly present the spectroscopic evolution of the recent outburst of the classical nova V1369 Cen, and the presence of a narrow absorption line identified as due to the resonance of neutral lithium at 6708 \AA. We also discuss the consequences for the chemical evolution of lithium in the Galaxy.
We report calculations of energy levels, radiative rates, collision strengths, and effective collision strengths for transitions among the lowest 25 levels of the n <= 5 configurations of He~II. The general-purpose relativistic atomic structure package (GRASP) and Dirac atomic R-matrix code (DARC) are adopted for the calculations. Radiative rates, oscillator strengths, and line strengths are reported for all electric dipole (E1), magnetic dipole (M1), electric quadrupole (E2), and magnetic quadrupole (M2) transitions among the 25 levels. Furthermore, collision strengths and effective collision strengths are listed for all 300 transitions among the above 25 levels over a wide energy (temperature) range up to 9 Ryd (10**5.4 K). Comparisons are made with earlier available results and the accuracy of the data is assessed.
Links to: arXiv, form interface, find, astro-ph, recent, 1704, contact, help (Access key information)
The discovery of multiple stellar populations in the Milky Way globular clusters (GCs) has stimulated various follow-up studies on helium-enhanced stellar populations. Here, we present the evolutionary population synthesis models for spectro-photometric evolution of simple stellar populations (SSPs) with varying initial helium abundance ($Y_{ini}$). We show that $Y_{ini}$ brings about dramatic changes in spectro-photometric properties of SSPs. Like the normal-helium SSPs, the integrated spectro-photometric evolution of helium-enhanced SSPs is also dependent on metallicity and age for a given $Y_{ini}$. We discuss the implications and the prospects for the helium-enhanced populations in relation to the second generation populations found in the Milky Way GCs. All models are available at \url{this http URL}.
PBC J2333.9-2343 is a giant radio galaxy which shows different characteristics at different wavebands that are difficult to explain within the actual generic schemes of unification of active galactic nuclei (AGN), thus being a good candidate to host different phases of nuclear activity. We aim at disentangling the nature of this AGN by using simultaneous multiwavelength data. We obtained data in 2015 from the Very Long Baseline Array (VLBA), the San Pedro M\'{a}rtir telescope, and the XMM-Newton}observatories. This allows the study of the nuclear parts of the galaxy through its morphology and spectra, as well as the analysis of the spectral energy distribution (SED). We also reanalysed optical data from the San Pedro M\'{a}rtir telescope from 2009 previously presented in Parisi et al. (2012) for a homogeneous comparison. At X-ray frequencies, the source is unabsorbed. The optical spectra are of a type 1.9 AGN, both in 2009 and 2015, although showing a broader component in 2015. The VLBA radio images show an inverted spectrum with self-absorbed, optically thick compact core and steep spectrum, optically thin jet. The SED resembles that of typical blazars and is best represented by an external Compton (EC) model with a viewing angle of $\sim$ 3-6 degrees. The apparent size of the large scale structure of PBC J2333.9-2343 must correspond to an intrinsic deprojected value of $\sim 7$ Mpc for $\theta_v<10^\circ$, and to $> 13$ Mpc for $\theta_v<5^\circ$, a value much larger than the bigger giant radio galaxy known (4.5 Mpc). The above arguments suggest that PBC J2333.9-2343 has undergone through a new episode of nuclear activity and that the direction of the new jet has changed, in the plane of the sky and is now pointing towards us, making this source from being a radio galaxy to become a blazar, a very exceptional case of restarting activity.
We characterize the near-infrared (NIR) dust attenuation for a sample of ~5500 local (z<0.1) star-forming galaxies and obtain an estimate of their average total-to-selective attenuation $k(\lambda)$. We utilize data from the United Kingdom Infrared Telescope (UKIRT) and the Two Micron All-Sky Survey (2MASS), which is combined with previously measured UV-optical data for these galaxies. The average attenuation curve is slightly lower in the far-UV than local starburst galaxies, by roughly 15%, but appears similar at longer wavelengths with a total-to-selective normalization at V-band of $R_V=3.67\substack{+0.44 \\ -0.35}$. Under the assumption of energy balance, the total attenuated energy inferred from this curve is found to be broadly consistent with the observed infrared dust emission ($L_{\rm{TIR}}$) in a small sample of local galaxies for which far-IR measurements are available. However, the significant scatter in this quantity among the sample may reflect large variations in the attenuation properties of individual galaxies. We also derive the attenuation curve for sub-populations of the main sample, separated according to mean stellar population age (via $D_n4000$), specific star formation rate, stellar mass, and metallicity, and find that they show only tentative trends with low significance, at least over the range which is probed by our sample. These results indicate that a single curve is reasonable for applications seeking to broadly characterize large samples of galaxies in the local Universe, while applications to individual galaxies would yield large uncertainties and is not recommended.
Broad absorption lines (BALs) found in a significant fraction of quasar spectra identify high-velocity outflows that might be present in all quasars and could be a major factor in feedback to galaxy evolution. Understanding the nature of these flows requires further constraints on their physical properties, including their column densities, for which well-studied BALs, such as CIV 1548,1551, typically provide only a lower limit because of saturation effects. Low-abundance lines, such as PV 1118,1128, indicate large column densities, implying outflows more powerful than measurements of CIV alone would indicate. We search through a sample of 2694 BAL quasars from the SDSS-III/BOSS DR9 quasar catalog for such absorption, and we identify 81 `definite' and 86 `probable' detections of PV broad absorption, yielding a firm lower limit of 3.0-6.2% for the incidence of such absorption among BAL quasars. The PV-detected quasars tend to have stronger CIV and SiIV absorption, as well as a higher incidence of LoBAL absorption, than the overall BAL quasar population. Many of the PV-detected quasars have CIV troughs that do not reach zero intensity (at velocities where PV is detected), confirming that the outflow gas only partially covers the UV continuum source. PV appears significantly in a composite spectrum of non-PV-detected BAL quasars, indicating that PV absorption (and large column densities) are much more common than indicated by our search results. Our sample of PV detections significantly increases the number of known PV detections, providing opportunities for follow-up studies to better understand BAL outflow energetics.
We present the SILVERRUSH program strategy and clustering properties investigated with $2,354$ Ly$\alpha$ emitters at $z=5.7$ and $6.6$ found in the early data of the Hyper Suprime-Cam (HSC) Subaru Strategic Program survey exploiting the carefully designed narrowband filters. We derive angular correlation functions with the unprecedentedly large samples of LAEs at $z=6-7$ over the large total area of $14-21$ deg$^2$ corresponding to $0.3-0.5$ comoving Gpc$^2$. We obtain the average large-scale bias values of $b_{\rm avg}=4.1\pm 0.2$ ($4.5\pm 0.6$) at $z=5.7$ ($z=6.6$) for $\gtrsim L^*$ LAEs, indicating the weak evolution of LAE clustering from $z=5.7$ to $6.6$. We compare the LAE clustering results with two independent theoretical models that suggest an increase of an LAE clustering signal by the patchy ionized bubbles at the epoch of reionization (EoR), and estimate the neutral hydrogen fraction to be $x_{\rm HI}=0.3^{+0.1}_{-0.2}$ at $z=6.6$. Based on the halo occupation distribution models, we find that the $\gtrsim L^*$ LAEs are hosted by the dark-matter halos with the average mass of $\log (\left < M_{\rm h} \right >/M_\odot) =11.1^{+0.2}_{-0.4}$ ($10.8^{+0.3}_{-0.5}$) at $z=5.7$ ($6.6$) with a Ly$\alpha$ duty cycle of 1 % or less, where the results of $z=6.6$ LAEs may be slightly biased, due to the increase of the clustering signal at the EoR. Our clustering analysis reveals the low-mass nature of $\gtrsim L^*$ LAEs at $z=6-7$, and that these LAEs probably evolve into massive super-$L^*$ galaxies in the present-day universe.
We present the results of our investigation of the star-forming potential in the Perseus star-forming complex. We build on previous starless core, protostellar core, and young stellar object (YSO) catalogs from Spitzer, Herschel, and SCUBA observations in the literature. We place the cores and YSOs within seven star-forming clumps based on column densities greater than 5x10^21 cm^-2. We calculate the mean density and free-fall time for 69 starless cores as 5.55x10^-19 gcm^-3 and 0.1 Myr,respectively, and we estimate the star formation rate for the near future as 150 Msun Myr^-1. According to Bonnor Ebert stability analysis, we find that majority of starless cores in Perseus are unstable. Broadly, these cores can collapse to form the next generation of stars. We found a relation between starless cores and YSOs, where the numbers of young protostars (Class 0 + Class I) are similar to the numbers of starless cores. This similarity, which shows a one-to-one relation, suggests that these starless cores may form the next generation of stars with approximately the same formation rate as the current generation, as identified by the Class 0 and Class I protostars. It follows that if such a relation between starless cores and any YSO stage exists, the SFR values of these two populations must be nearly constant. In brief, we propose that this one-to-one relation is an important factor in better understanding the star formation process within a cloud.
The distribution of N/O abundance ratios calculated by the detailed modelling of different galaxy spectra at z<4 is investigated. Supernova (SN) and long gamma-ray-burst (LGRB) host galaxies cover different redshift domains. N/O in SN hosts increases due to secondary N production towards low z (0.01) accompanying the growing trend of active galaxies (AGN, LINER). N/O in LGRB hosts decreases rapidly between z>1 and z ~0.1 following the N/H trend and reach the characteristic N/O ratios calculated for the HII regions in local and nearby galaxies. The few short period GRB (SGRB) hosts included in the galaxy sample show N/H <0.04 solar and O/H solar. They seem to continue the low bound N/H trend of SN hosts at z<0.3. The distribution of N/O as function of metallicity for SN and LGRB hosts is compared with star chemical evolution models. The results show that several LGRB hosts can be explained by star multi-bursting models when 12+log(O/H) <8.5, while some objects follow the trend of continuous star formation models. N/O in SN hosts at log(O/H)+12 <8.5 are not well explained by stellar chemical evolution models calculated for starburst galaxies. At 12+log(O/H) >8.5 many different objects are nested close to O/H solar with N/O ranging between the maximum corresponding to starburst galaxies and AGN and the minimum corresponding to HII regions and SGRB.
We report HST/ACS observations of two nearby gas-rich dwarf galaxies: DDO 161 and UGCA 319. Their distances determined via the Tip of the Red Giant Branch are 6.03 (-0.21/+0.29) Mpc and 5.75+-0.18 Mpc, respectively. The galaxies form an isolated pair dynamically well separated from the nearest neighbors: KK 176 (7.28+-0.29 Mpc) and NGC 5068 (5.16+-0.21 Mpc). All four galaxies have a bulk spatial peculiar velocity towards the Virgo cluster of ~158+-17 km/s in the Local Group rest frame and ~330 km/s with respect to the cluster center.
We report the discovery of tidal tails around the two outer halo globular clusters, Eridanus and Palomar 15, based on $gi$-band images obtained with DECam at the CTIO 4-m Blanco Telescope. The tidal tails are among the most remote stellar streams presently known in the Milky Way halo. Cluster members have been determined from the color-magnitude diagrams and used to establish the radial density profiles, which show, in both cases, a strong departure in the outer regions from the best-fit King profile. Spatial density maps reveal tidal tails stretching out on opposite sides of both clusters, extending over a length of $\sim$760 pc for Eridanus and $\sim$1160 pc for Palomar 15. The great circle projected from the Palomar 15 tidal tails encompasses the Galactic Center, while that for Eridanus passes close to four dwarf satellite galaxies, one of which (Sculptor) is at a comparable distance to that of Eridanus.
NGC\,2403, NGC\,300 and M33 are three nearby pure-disc galaxies with similar stellar mass in different environments, they are benchmarks for understanding late-type spiral galaxies in different environments. The chemical evolution and growth of their discs are investigated by using the simple chemical evolution model, in which their discs are assumed to originate and grow through the accretion of the primordial gas, and the gas outflow process is also taken into account. Through comparative study of the best-fitting model predicted star formation histories for them, we hope to derive a picture of the local environment on the evolution and star formation histories of galaxies and whether or not the isolated galaxies follow similar evolution history. Our results show that these three galaxies accumulated more than 50 percent of their stellar mass at $z\,<\,1$. It can be also found that the metallicity gradients in isolated spiral galaxies NGC\,2403 and NGC\,300 are similar and obviously steeper than that in M33, when the metallicity gradients are expressed in ${\rm dex}\,R_{\rm 25}^{-1}$. The similar metallicity gradients in NGC\,2403 and NGC\,300 indicate that they may experience similar chemical evolutionary histories. The principal epoch of star formation on the discs of NGC\,2403 and NGC\,300 is earlier than that on the disc of M33, and the mean age of stellar populations along the discs of both NGC\,2403 and NGC\,300 is older than that of M33. Our results indicates that the evolution and star formation history of a galaxy indeed depends on its local environment, at least for galaxies with stellar mass of $10^{9.2}\,\rm M_{\odot}\sim10^{9.7}\,\rm M_{\odot}$.
We present a new catalog of narrow-line Seyfert 1 (NLSy1) galaxies from the Sloan Digital Sky Survey Data Release 12 (SDSS DR12). This was obtained by a systematic analysis through modeling of the continuum and emission lines of the spectra of all the 68,859 SDSS DR12 objects that are classified as "QSO" by the SDSS spectroscopic pipeline with z < 0.8 and a median signal-to-noise ratio (S/N) > 2 per pixel. This catalog contains a total of 11,101 objects, which is about 5 times larger than the previously known NLSy1 galaxies. Their monochromatic continuum luminosity at 5100 A is found to be strongly correlated with H-beta, H-alpha, and [O III] emission line luminosities. The optical Fe II strength in NLSy1 galaxies is about two times larger than the broad- line Seyfert 1 (BLSy1) galaxies. About 5% of the catalog sources are detected in the FIRST survey. The Eddington ratio (XEdd) of NLSy1 galaxies has an average of log XEdd of -0.34, much higher than -1.03 found for BLSy1 galaxies. Their black hole masses (MBH) have an average of log MBH of 6.9 Msun, which is less than BLSy1 galaxies, which have an average of log MBH of 8.0 Msun. The MBH of NLSy1 galaxies is found to be correlated with their host galaxy velocity dispersion. Our analysis suggests that geometrical effects playing an important role in defining NLSy1 galaxies and their MBH deficit is perhaps due to their lower inclination compared to BLSy1 galaxies.
We reinvestigate a claimed sample of 22 X-ray detected active galactic nuclei (AGN) at redshifts z > 4, which has reignited the debate as to whether young galaxies or AGN reionized the Universe. These sources lie within the GOODS-S/CANDELS field, and we examine both the robustness of the claimed X-ray detections (within the Chandra 4Ms imaging) and perform an independent analysis of the photometric redshifts of the optical/infrared counterparts. We confirm the reality of only 15 of the 22 reported X-ray detections, and moreover find that only 12 of the 22 optical/infrared counterpart galaxies actually lie robustly at z > 4. Combining these results we find convincing evidence for only 7 X-ray AGN at z > 4 in the GOODS-S field, of which only one lies at z > 5. We recalculate the evolving far-UV (1500 Angstrom) luminosity density produced by AGN at high redshift, and find that it declines rapidly from z = 4 to z = 6, in agreement with several other recent studies of the evolving AGN luminosity function. The associated rapid decline in inferred hydrogen-ionizing emissivity contributed by AGN falls an order-of-magnitude short of the level required to maintain hydrogen ionization at z ~ 6. We conclude that all available evidence continues to favour a scenario in which young galaxies reionized the Universe, with AGN making, at most, a very minor contribution to cosmic hydrogen reionization.
We present a statistical study on the [C I]($^{3} \rm P_{1} \rightarrow {\rm ^3 P}_{0}$), [C I] ($^{3} \rm P_{2} \rightarrow {\rm ^3 P}_{1}$) lines (hereafter [C I] (1$-$0) and [C I] (2$-$1), respectively) and the CO (1$-$0) line for a sample of (ultra)luminous infrared galaxies [(U)LIRGs]. We explore the correlations between the luminosities of CO (1$-$0) and [C I] lines, and find that $L'_\mathrm{CO(1-0)}$ correlates almost linearly with both $L'_ \mathrm{[CI](1-0)}$ and $L'_\mathrm{[CI](2-1)}$, suggesting that [C I] lines can trace total molecular gas mass at least for (U)LIRGs. We also investigate the dependence of $L'_\mathrm{[CI](1-0)}$/$L'_\mathrm{CO(1-0)}$, $L'_\mathrm{[CI](2-1)}$/$L'_\mathrm{CO(1-0)}$ and $L'_\mathrm{[CI](2-1)}$/$L'_\mathrm{[CI](1-0)}$ on the far-infrared color of 60-to-100 $\mu$m, and find non-correlation, a weak correlation and a modest correlation, respectively. Under the assumption that these two carbon transitions are optically thin, we further calculate the [C I] line excitation temperatures, atomic carbon masses, and the mean [C I] line flux-to-H$_2$ mass conversion factors for our sample. The resulting $\mathrm{H_2}$ masses using these [C I]-based conversion factors roughly agree with those derived from $L'_\mathrm{CO(1-0)}$ and CO-to-H$_2$ conversion factor.
We use two catalogues, a Herschel catalogue selected at 500 mu (HerMES) and
an IRAS catalogue selected at 60 mu (RIFSCz), to contrast the sky at these two
wavelengths.
Both surveys demonstrate the existence of extreme starbursts, with
star-formation rates (SFRs) > 5000 Msun/yr. There appears to be a maximum
(un-lensed) star-formation rate of 30,000 Msun/yr. IRAS sources with estimates
higher than this are either lensed systems, blazars, or erroneous photometric
redshifts.
At redshifts 3 to 5, the time-scale for the Herschel galaxies to make their
current mass of stars at their present rate of formation ~ 10^8 yrs, so these
galaxies are making a significant fraction of their stars in the current
star-formation episode. Using dust mass as a proxy for gas mass, the Herschel
galaxies at redshift 3 to 5 have gas masses comparable to their mass in stars.
We have plotted here the individual spectral energy distributions (SEDs) for
the 58 extreme starbursts in our Herschel survey, for which we have more
complete SED information. Over 50% are QSOs or have an AGN dust torus, i.e. are
Type 2 AGN, but in all cases the infrared luminosity is dominated by a
starburst component. We derive a mean covering factor for AGN dust as a
function of redshift and derive black hole masses and black hole accretion
rates. There is a universal ratio of black-hole mass to stellar mass, ~
10^{-3}, driven by the strong period of star-formation and black-hole growth at
z = 1-5.
We analyze three-dimensional hydrodynamical simulations of the interaction of jets and the bubbles they inflate with the intra-cluster medium (ICM), and show that the heating of the ICM by mixing hot bubble gas with the ICM operates over tens of millions of years, and hence can smooth the sporadic activity of the jets. The inflation process of hot bubbles by propagating jets forms many vortices, and these vortices mix the hot bubble gas with the ICM. The mixing, hence the heating of the ICM, starts immediately after the jets are launched, but continues for tens of millions of years. We suggest that the smoothing of the active galactic nucleus (AGN) sporadic activity by the long-lived vortices accounts for the recent finding of a gentle energy coupling between AGN heating and the ICM.
We present the first self-consistent chemodynamical model fitted to reproduce data for the galactic bulge, bar and inner disk. We extend the Made-to-Measure method to an augmented phase-space including the metallicity of stars, and show its first application to the bar region of the Milky Way. Using data from the ARGOS and APOGEE (DR12) surveys, we adapt the recent dynamical model from Portail et al. to reproduce the observed spatial and kinematic variations as a function of metallicity, thus allowing the detailed study of the 3D density distributions, kinematics and orbital structure of stars in different metallicity bins. We find that metal-rich stars with [Fe/H] > -0.5 are strongly barred and have dynamical properties that are consistent with a common disk origin. Metal-poor stars with [Fe/H] < -0.5 show strong kinematic variations with metallicity, indicating varying contributions from the underlying stellar populations. Outside the central kpc, metal-poor stars are found to have the density and kinematics of a thick disk while in the inner kpc, evidence for an extra concentration of metal-poor stars is found. Finally, the combined orbit distributions of all metallicities in the model naturally reproduce the observed vertex deviations in the bulge. This paper demonstrates the power of Made-to-Measure chemodynamical models, that when extended to other chemical dimensions will be very powerful tools to maximize the information obtained from large spectroscopic surveys such as APOGEE, GALAH and MOONS.
We interpret the recent discovery of a preferable VLBI/Gaia offset direction for active galactic nuclei (AGNs) along the parsec-scale radio jet as a manifestation of their optical structure on scales of 1 to 100 milliarcseconds. The extended jet structure affects the Gaia position stronger than the VLBI position due to the difference in observing techniques. Gaia detects total power while VLBI measures the correlated quantity, visibility. The synergy of VLBI that is sensitive to the position of the most compact source component, usually associated with the opaque radio core, and Gaia that is sensitive to the centroid of optical emission, opens a window of opportunity to study optical jets at milliarcsecond resolution, two orders of magnitude finer than the resolution of most existing optical instruments. We have demonstrated that strong variability of optical jets is able to cause a jitter comparable to the VLBI/Gaia offsets at a quiet state, i.e. several milliarcseconds. We have shown that VLBI/Gaia position jitter correlation against the AGN optical light curve may help to locate the region where the flare occurred, estimate its distance from the super-massive black hole and the ratio of the flux density in the flaring region to the total flux density.
Gravitational waves from binary black hole mergers provide a new probe of massive-star evolution and the formation channels of binary compact objects. By coupling the growing sample of binary black hole systems with realistic population synthesis models, we can begin to constrain the parameters of such models and glean unprecedented knowledge of the inherent physical processes that underpin binary stellar evolution. In this study, we apply a hierarchical Bayesian model on mass measurements from a synthetic gravitational-wave sample to constrain the physical prescriptions in population models and the relative fractions of systems generated from various channels. We employ population models of two canonical formation scenarios in our analysis --- isolated binary evolution involving common-envelope evolution and dynamical formation within globular clusters --- with model variations for different black-hole natal kick magnitude prescriptions. We show that solely with chirp mass measurements, it is possible to constrain natal kick prescriptions and the relative fraction of systems detected from each formation channel with $\mathcal{O}(100)$ of confident detections. This framework can be extended to include additional formation scenarios, model parameters, and measured properties of the compact binary.
We have constructed the most-comprehensive catalog of photometry and proper motions ever assembled for a globular cluster (GC). The core of $\omega$Cen has been imaged over 650 times through WFC3's UVIS and IR channels for the purpose of detector calibration. There exist from 4 to over 60 exposures through each of 26 filters, stretching continuously from F225W in the UV to F160W in the infrared. Furthermore, the 11-year baseline between these data and a 2002 ACS survey has allowed us to more than double the proper-motion accuracy and triple the number of well-measured stars compared to our previous groundbreaking effort. This totally unprecedented complete spectral coverage for over 470,000 stars within the cluster's core, from the tip of the red-giant branch down to the white dwarfs, provides the best astro-photometric observational data base yet to understand the multiple-population phenomenon in any GC. In this first paper of the series we describe in detail the data-reduction processes and deliver the astro-photometric catalog to the astronomical community.
The propagation of charged cosmic rays through the Galactic environment
influences all aspects of the observation at Earth. Energy spectrum,
composition and arrival directions are changed due to deflections in magnetic
fields and interactions with the interstellar medium. Today the transport is
simulated with different simulation methods either based on the solution of a
transport equation (multi-particle picture) or a solution of an equation of
motion (single-particle picture).
We developed a new module for the publicly available propagation software
CRPropa 3.1, where we implemented an algorithm to solve the transport equation
using stochastic differential equations. This technique allows us to use a
diffusion tensor which is anisotropic with respect to an arbitrary magnetic
background field.
In this paper, we present the new low-energy propagation code together with
validation procedures that are developed to proof the accuracy of the new
implementation. Furthermore, we show first examples of the cosmic ray density
evolution, which depends strongly on the ratio of the parallel
$\kappa_\parallel$ and perpendicular $\kappa_\perp$ diffusion coefficients.
This dependency is systematically examined as well the influence of the
particle rigidity on the diffusion process.
In molecular outflows from forming low-mass protostars, most oxygen is expected to be locked up in water. However, Herschel observations have shown that typically an order of magnitude or more of the oxygen is still unaccounted for. To test if the oxygen is instead in atomic form, SOFIA-GREAT observed the R1 position of the bright molecular outflow from NGC1333-IRAS4A. The [OI] 63 um line is detected and spectrally resolved. From an intensity peak at +15 km/s, the intensity decreases until +50 km/s. The profile is similar to that of high-velocity (HV) H2O and CO 16-15, the latter observed simultaneously with [OI]. A radiative transfer analysis suggests that ~15% of the oxygen is in atomic form toward this shock position. The CO abundance is inferred to be ~10^-4 by a similar analysis, suggesting that this is the dominant oxygen carrier in the HV component. These results demonstrate that a large portion of the observed [OI] emission is part of the outflow. Further observations are required to verify whether this is a general trend.
We present a homogeneous set of accurate atmospheric parameters for a complete sample of very and extremely metal-poor stars in the dwarf spheroidal galaxies (dSphs) Sculptor, Ursa Minor, Sextans, Fornax, Bo\"otes I, Ursa Major II, and Leo IV. We also deliver a Milky Way (MW) comparison sample of giant stars covering the -4 < [Fe/H] < -1.7 metallicity range. We show that, in the [Fe/H] > -3.5 regime, the non-local thermodynamic equilibrium (NLTE) calculations with non-spectroscopic effective temperature (Teff) and surface gravity (log~g) based on the photometric methods and known distance provide consistent abundances of the Fe I and Fe II lines. This justifies the Fe I/Fe II ionisation equilibrium method to determine log g for the MW halo giants with unknown distance. The atmospheric parameters of the dSphs and MW stars were checked with independent methods. In the [Fe/H] > -3.5 regime, the Ti I/Ti II ionisation equilibrium is fulfilled in the NLTE calculations. In the log~g - Teff plane, all the stars sit on the giant branch of the evolutionary tracks corresponding to [Fe/H] = -2 to -4, in line with their metallicities. For some of the most metal-poor stars of our sample, we hardly achieve consistent NLTE abundances from the two ionisation stages for both iron and titanium. We suggest that this is a consequence of the uncertainty in the Teff-colour relation at those metallicities. The results of these work provide the base for a detailed abundance analysis presented in a companion paper.
Context. Globular clusters host stars with chemical peculiarities. The associated helium enrichment is expected to affect the evolution of stars, in general, and of low-mass stars, and in particular the progenitors of white dwarfs (WDs). Aims. We investigate the effects of different initial helium contents on the properties of white dwarfs such as their masses, compositions, and the time since their formation. Results. In a globular cluster at a given age, the He enrichment implies lower initial masses for stars at a given phase. Thus it leads to a decrease of the masses of white dwarfs reaching the cooling sequence. In addition the He enrichment increases the total mass and number of white dwarfs and eventually allows the presence of He white dwarf from single progenitors. Conclusions. The low He enrichment determined in most globular clusters with different methods results in negligible effects on the white dwarf properties. However, in the few globular clusters that display a high He enrichment, this may significantly affect the characteristics of the white dwarfs. In NGC 2808 and {\omega} Centauri the high He enrichment even leads to the formation of He white dwarfs from single He-rich progenitors. Therefore investigating the white dwarf mass domain in globular clusters with a high He enrichment would provide an additional indirect way to measure and constrain the He enrichment degree.
Links to: arXiv, form interface, find, astro-ph, recent, 1704, contact, help (Access key information)
We present a family of self-consistent axisymmetric rotating globular cluster models which are fitted to spectroscopic data for NGC 362, NGC 1851, NGC 2808, NGC 4372, NGC 5927 and NGC 6752 to provide constraints on their physical and kinematic properties, including their rotation signals. They are constructed by flattening Modified Plummer profiles, which have the same asymptotic behaviour as classical Plummer models, but can provide better fits to young clusters due to a slower turnover in the density profile. The models are in dynamical equilibrium as they depend solely on the action variables. We employ a fully Bayesian scheme to investigate the uncertainty in our model parameters (including mass-to-light ratios and inclination angles) and evaluate the Bayesian evidence ratio for rotating to non-rotating models. We find convincing levels of rotation only in NGC 2808. In the other clusters, there is just a hint of rotation (in particular, NGC 4372 and NGC 5927), as the data quality does not allow us to draw strong conclusions. Where rotation is present, we find that it is confined to the central regions, within radii of $R \leq 2 r_h$. As part of this work, we have developed a novel q-Gaussian basis expansion of the line-of-sight velocity distributions, from which general models can be constructed via interpolation on the basis coefficients.
The spatial and kinematic distribution of warm gas in and around the Coma Cluster is presented through observations of Lyman-alpha absorbers using background QSOs. Updates to the Lyman-alpha absorber distribution found in Yoon et al. (2012) for the Virgo Cluster are also presented. At 0.2-2.0 R_vir of Coma we identify 14 Lyman-alpha absorbers (N_HI = 10^{12.8-15.9} cm^-2) towards 5 sightlines and no Lyman-alpha absorbers along 3 sightlines within 3\sigmav_coma. For both Coma and Virgo, most Lyman-alpha absorbers are found outside the virial radius or beyond 1\sigmav consistent with them largely representing the infalling intergalactic medium. The few exceptions in the central regions can be associated with galaxies. The Lyman-alpha absorbers avoid the hot ICM, consistent with the infalling gas being shock-heated within the cluster. The massive dark matter halos of clusters do not show the increasing column density with decreasing impact parameter relationship found for the smaller mass galaxy halos. In addition, while the covering fraction within R_vir is lower for clusters than galaxies, beyond R_vir the covering fraction is somewhat higher for clusters. The velocity dispersion of the absorbers compared to the galaxies is higher for Coma, consistent with the absorbers tracing additional turbulent gas motions in the cluster outskirts. The results are overall consistent with cosmological simulations, with the covering fraction being high in the observations standing out as the primary discrepancy.
We investigate the distribution of central velocity dispersions for quiescent galaxies in the SDSS at $0.03 \leq z \leq 0.10$. To construct the field velocity dispersion function (VDF), we construct a velocity dispersion complete sample of quiescent galaxies with Dn4000$ > 1.5$. The sample consists of galaxies with central velocity dispersion larger than the velocity dispersion completeness limit of the SDSS survey. Our VDF measurement is consistent with previous field VDFs for $\sigma > 200$ km s$^{-1}$. In contrast with previous results, the VDF does not decline significantly for $\sigma < 200$ km s$^{-1}$. The field and the similarly constructed cluster VDFs are remarkably flat at low velocity dispersion ($\sigma < 250$ km s$^{-1}$). The cluster VDF exceeds the field for $\sigma < 250$ km s$^{-1}$ providing a measure of the relatively larger number of massive subhalos in clusters. The VDF is a probe of the dark matter halo distribution because the measured central velocity dispersion may be directly proportional to the dark matter velocity dispersion. Thus the VDF provides a potentially powerful test of simulations for models of structure formation.
The recent discovery of thousands of ultra diffuse galaxies (UDGs) in nearby galaxy clusters has opened a new window into the process of galaxy formation and evolution. Several scenarios have been proposed to explain the formation history of UDGs, and their ability to survive in the harsh cluster environments. A key requirement to distinguish between these scenarios is a measurement of their halo masses which, due to their low surface brightnesses, has proven difficult if one relies on stellar tracers of the potential. We exploit weak gravitational lensing, a technique that does not depend on these baryonic tracers, to measure the average subhalo mass of 784 UDGs selected in 18 clusters at $z\leq0.09$. Our sample of UDGs has a median stellar mass $\langle m_\star\rangle=2\times10^8\,\mathrm{M}_\odot$ and a median effective radius $\langle r_\mathrm{eff}\rangle=2.8$ kpc. We constrain the average mass of subhaloes within 30 kpc to $\log m_\mathrm{UDG}(r<30\,\mathrm{kpc})/\mathrm{M}_\odot\leq10.99$ at 95 per cent credibility, implying an effective virial mass $\log m_{200}/\mathrm{M}_\odot\leq11.80$, and a lower limit on the stellar mass fraction within 10 kpc of 1.0 per cent. Such mass is consistent with a simple extrapolation of the subhalo-to-stellar mass relation of typical satellite galaxies in massive clusters. However, our analysis is not sensitive to scatter about this mean mass; the possibility remains that extreme UDGs reside in haloes as massive as the Milky Way.
Accretion of gas and the interaction of matter and radiation are at the heart of many questions pertaining to black hole (BH) growth and the coevolution of massive BHs and their host galaxies. In order to answer them it is critical to quantify how the ionizing radiation that emanates from the innermost regions of the BH accretion flow couples to the surrounding medium and how it regulates the BH fueling. In this work we use high resolution 3-dimensional (3D) radiation-hydrodynamic simulations with the code Enzo, equipped with adaptive ray tracing module Moray, to investigate BH accretion of cold gas regulated by radiative feedback. Our simulations reproduce findings from an earlier generation of 1D and 2D simulations, that the accretion powered UV and X-ray radiation forms a highly ionized bubble, which leads to suppression of BH accretion rate characterized by quasi-periodic outbursts. A new feature revealed by the 3D simulations is the highly turbulent nature of the gas flow in vicinity of the ionization front. Because turbulence is efficient in replenishing the gas in the low-density region inside the ionization front, the 3D simulations show oscillations in the accretion rate of only ~2-3 orders of magnitude, significantly smaller than 1D and 2D models. We calculate the energy budget of the gas flow and find that turbulence is the main contributor to the kinetic energy of the gas but corresponds to less than 10% of its thermal energy. The turbulence therefore does not contribute significantly to the pressure support of the gas but is the key factor that sets the level of BH fueling during quiescent periods between the accretion outbursts.
We compare the existent methods including the minimum spanning tree based method and the local stellar density based method in measuring mass segregation of star clusters. We find that the minimum spanning tree method reflects more the compactness, which represents the global spatial distribution of massive stars, while the local stellar density method reflects more the crowdedness, which provides the local gravitational potential information. It is suggested to measure the local and the global mass segregation simultaneously. We also develop a hybrid method which takes both aspects into account. This hybrid method balances the local and the global mass segregation in the sense that if the predominant one is caused either by dynamical evolution or purely accidental, especially when the knowledge is unknown a priori. In addition, we test our prescriptions with numerical models and show the impact of binaries, in estimating the mass segregation value. As an application, we apply these methods on the Orion Nebula Cluster and the Taurus Cluster observations. We find that the Orion Nebula Cluster is significantly mass segregated down to the 20th most massive stars. In contrast, the massive stars of the Taurus cluster are sparsely distributed in many different sub-clusters and show a low degree of compactness. The massive stars of Taurus are also found to be distributed in the high-density regions of sub-clusters, showing significant mass segregation at sub-cluster scale. Finally, we apply these methods to discuss the possible mechanisms of the dynamical evolution of simulated sub-structured star clusters.
The paucity of hypervelocity stars (HVSs) known to date has severely hampered their potential to investigate the stellar population of the Galactic Centre and the Galactic Potential. The first Gaia data release gives an opportunity to increase the current sample. The challenge is of course the disparity between the expected number of hypervelocity stars and that of bound background stars (around 1 in $10^6$). We have applied a novel data mining algorithm based on machine learning techniques, an artificial neural network, to the Tycho-Gaia astrometric solution (TGAS) catalogue. With no pre-selection of data, we could exclude immediately $\sim 99 \%$ of the stars in the catalogue and find 80 candidates with more than $90\%$ predicted probability to be HVSs, based only on their position, proper motions, and parallax. We have cross-checked our findings with other spectroscopic surveys, determining radial velocities for 30 and spectroscopic distances for 5 candidates. In addition, follow-up observations have been carried out at the Isaac Newton Telescope for 22 stars, for which we obtained radial velocities and distance estimates. We discover 14 stars with a total velocity in the Galactic rest frame > 400 km/s, and 5 of these have a probability $>50\%$ of being unbound from the Milky Way. Tracing back their orbits in different Galactic potential models we find one possible unbound HVS with velocity $\sim$ 520 km/s, 5 bound HVSs, and, notably, 5 runaway stars with median velocity between 400 and 780 km/s. At the moment, uncertainties in the distance estimates and ages are too large to confirm the nature of our candidates by narrowing down their ejection location, and we wait for future Gaia releases to validate the quality of our sample. This test successfully demonstrates the feasibility of our new data mining routine.
Starting from a summary of detection statistics of our recent X-shooter campaign, we review the major surveys, both space and ground based, for emission counterparts of high-redshift damped Ly$\alpha$ absorbers (DLAs) carried out since the first detection 25 years ago. We show that the detection rates of all surveys are precisely reproduced by a simple model in which the metallicity and luminosity of the galaxy associated to the DLA follow a relation of the form, ${\rm M_{UV}} = -5 \times \left(\,[{\rm M/H}] + 0.3\, \right) - 20.8$, and the DLA cross-section follows a relation of the form $\sigma_{DLA} \propto L^{0.8}$. Specifically, our spectroscopic campaign consists of 11 DLAs preselected based on their equivalent width of SiII $\lambda1526$ to have a metallicity higher than [Si/H] > -1. The targets have been observed with the X-shooter spectrograph at the Very Large Telescope to search for emission lines around the quasars. We observe a high detection rate of 64% (7/11), significantly higher than the typical $\sim$10% for random, HI-selected DLA samples. We use the aforementioned model, to simulate the results of our survey together with a range of previous surveys: spectral stacking, direct imaging (using the `double DLA' technique), long-slit spectroscopy, and integral field spectroscopy. Based on our model results, we are able to reconcile all results. Some tension is observed between model and data when looking at predictions of Ly$\alpha$ emission for individual targets. However, the object to object variations are most likely a result of the significant scatter in the underlying scaling relations as well as uncertainties in the amount of dust which affects the emission.
We present an unprecedentedly large catalog consisting of 2,354 >~ L^* Lya emitters (LAEs) at z=5.7 and 6.6 on the 13.8 and 21.2 deg^2 sky, respectively, that are identified by the SILVERRUSH program with the first narrowband imaging data of the Hyper Suprime-Cam (HSC) survey. We confirm that the LAE catalog is reliable on the basis of 97 LAEs whose spectroscopic redshifts are already determined by this program and the previous studies. This catalogue is also available on-line. Based on this catalogue, we derive the rest-frame Lya equivalent-width distributions of LAEs at z=5.7 and 6.6 that are reasonably explained by the exponential profiles with the scale lengths of 72+/-19 and 119+/-4 A, respectively, showing the increase trend towards high-z. We find that ~700 LAEs with a large equivalent width (LEW) of >~ 240 A are candidates of young-metal poor galaxies and AGNs. We also find that the fraction of LEW LAEs to all ones is moderately large, ~30%. Our LAE catalog includes 11 Lya blobs (LABs) that are LAEs with spatially extended Lya emission whose profile is clearly distinguished from those of stellar objects at the >~ 3sigma level. The number density of the LABs at z=6-7 is ~10^{-7}-10^{-6} Mpc-3, being ~10-100 times lower than those claimed for LABs at z~2-3, suggestive of disappearing LABs at z>~6, albeit with the different selection methods and criteria for the low and high-z LABs.
We present a near infrared study of the spectral components of the continuum in the inner 500$\times$500 pc$^2$ of the nearby Seyfert galaxy Mrk573 using adaptive optics near-infrared integral field spectroscopy with the instrument NIFS of the Gemini North Telescope at a spatial resolution of $\sim$50 pc. We performed spectral synthesis using the {\sc starlight} code and constructed maps for the contributions of different age components of the stellar population: young ($age\leq100$ Myr), young-intermediate ($100<age\leq700$ Myr), intermediate-old ($700$ Myr $<age\leq2$ Gyr) and old ($age>2$ Gyr) to the near-IR K-band continuum, as well as their contribution to the total stellar mass. We found that the old stellar population is dominant within the inner 250 pc, while the intermediate age components dominate the continuum at larger distances. A young stellar component contributes up to $\sim$20% within the inner $\sim$70 pc, while hot dust emission and featureless continuum components are also necessary to fit the nuclear spectrum, contributing up to 20% of the K-band flux there. The radial distribution of the different age components in the inner kiloparsec of Mrk573 is similar to those obtained by our group for the Seyfert galaxies Mrk1066, Mrk1157 and NGC1068 in previous works using a similar methodology. Young stellar populations ($\leq$100 Myr) are seen in the inner 200-300 pc for all galaxies contributing with $\ge$20% of the K-band flux, while the near-IR continuum is dominated by the contribution of intermediate-age stars ($t=$100 Myr-2 Gyr) at larger distances. Older stellar populations dominate in the inner 250 pc.
We mine the Tycho-{\it Gaia} astrometric solution (TGAS) catalog for wide stellar binaries by matching positions, proper motions, and astrometric parallaxes. We separate genuine binaries from unassociated stellar pairs through a Bayesian formulation that includes correlated uncertainties in the proper motions and parallaxes. Rather than relying on assumptions about the structure of the Galaxy, we calculate Bayesian priors and likelihoods based on the nature of Keplerian orbits and the TGAS catalog itself. We calibrate our method using radial velocity measurements and obtain 6196 high-confidence candidate wide binaries with projected separations $s\lesssim1$ pc. The normalization of this distribution suggests that at least 0.6\% of TGAS stars have an associated, distant TGAS companion in a wide binary. We demonstrate that {\it Gaia}'s astrometry is precise enough that it can detect projected orbital velocities in wide binaries with orbital periods as large as 10$^6$ yr. For pairs with $s\ \lesssim\ 4\times10^4$~AU, characterization of random alignments indicate our contamination to be $\approx$5\%. For $s \lesssim 5\times10^3$~AU, our distribution is consistent with \"{O}pik's Law. At larger separations, the distribution is steeper and consistent with a power-law $P(s)\propto s^{-1.6}$; there is no evidence in our data of any bimodality in this distribution for $s \lesssim$ 1 pc. Using radial velocities, we demonstrate that at large separations, i.e., of order $s \sim$ 1 pc and beyond, any potential sample of genuine wide binaries in TGAS cannot be easily distinguished from ionized former wide binaries, moving groups, or contamination from randomly aligned stars.
We present the highest spatial resolution ALMA observations to date of the Class I protostar WL 17 in the $\rho$ Ophiuchus L1688 molecular cloud complex, which show that it has a 12 AU hole in the center of its disk. We consider whether WL 17 is actually a Class II disk being extincted by foreground material, but find that such models do not provide a good fit to the broadband SED and also require such high extinction that it would presumably arise from dense material close to the source such as a remnant envelope. Self-consistent models of a disk embedded in a rotating collapsing envelope can nicely reproduce both the ALMA 3 mm observations and the broadband SED of WL 17. This suggests that WL 17 is a disk in the early stages of its formation, and yet even at this young age the inner disk has been depleted. Although there are multiple pathways for such a hole to be created in a disk, if this hole were produced by the formation of planets it could place constraints on the timescale for the growth of planetesimals in protoplanetary disks.
We discuss how contemporary multiwavelength observations of young OB-dominated clusters address long-standing astrophysical questions: Do clusters form rapidly or slowly with an age spread? When do clusters expand and disperse to constitute the field star population? Do rich clusters form by amalgamation of smaller subclusters? What is the pattern and duration of cluster formation in massive star forming regions (MSFRs)? Past observational difficulties in obtaining good stellar censuses of MSFRs have been alleviated in recent studies that combine X-ray and infrared surveys to obtain rich, though still incomplete, censuses of young stars in MSFRs. We describe here one of these efforts, the MYStIX project, that produced a catalog of 31,784 probable members of 20 MSFRs. We find that age spread within clusters are real in the sense that the stars in the core formed after the cluster halo. Cluster expansion is seen in the ensemble of (sub)clusters, and older dispersing populations are found across MSFRs. Direct evidence for subcluster merging is still unconvincing. Long-lived, asynchronous star formation is pervasive across MSFRs.
The formation process of massive stars is still poorly understood. Massive young stellar objects (mYSOs) are deeply embedded in their parental clouds, they are rare and thus typically distant, and their reddened spectra usually preclude the determination of their photospheric parameters. M17 is one of the best studied HII regions in the sky, is relatively nearby, and hosts a young stellar population. With X-shooter on the ESO Very Large Telescope we have obtained optical to near-infrared spectra of candidate mYSOs, identified by Hanson et al. (1997), and a few OB stars in this region. The large wavelength coverage enables a detailed spectroscopic analysis of their photospheres and circumstellar disks. We confirm the pre-main sequence (PMS) nature of six of the stars and characterise the O stars. The PMS stars have radii consistent with being contracting towards the main sequence and are surrounded by a remnant accretion disk. The observed infrared excess and the (double-peaked) emission lines provide the opportunity to measure structured velocity profiles in the disks. We compare the observed properties of this unique sample of young massive stars with evolutionary tracks of massive protostars by Hosokawa & Omukai (2009), and propose that these mYSOs near the western edge of the HII region are on their way to become main-sequence stars ($\sim 6 - 20$ $M_{\odot}$) after having undergone high mass-accretion rates (${\dot{M}_{\rm acc}} \sim 10^{-4} - 10^{-3}$ $M_{\odot}$ $\rm yr^{-1}$). Their spin distribution upon arrival at the zero age main sequence (ZAMS) is consistent with that observed for young B stars, assuming conservation of angular momentum and homologous contraction.
The discovery of radionuclides like 60Fe with half-lives of million years in deep-sea crusts and sediments offers the unique possibility to date and locate nearby supernovae. We want to quantitatively establish that the 60Fe enhancement is the result of several supernovae which are also responsible for the formation of the Local Bubble, our Galactic habitat. We performed three-dimensional hydrodynamic adaptive mesh refinement simulations (with resolutions down to subparsec scale) of the Local Bubble and the neighbouring Loop I superbubble in different homogeneous, self-gravitating environments. For setting up the Local and Loop I superbubble, we took into account the time sequence and locations of the generating core-collapse supernova explosions, which were derived from the mass spectrum of the perished members of certain stellar moving groups. The release of 60Fe and its subsequent turbulent mixing process inside the superbubble cavities was followed via passive scalars, where the yields of the decaying radioisotope were adjusted according to recent stellar evolution calculations. The models are able to reproduce both the timing and the intensity of the 60Fe excess observed with rather high precision, provided that the external density does not exceed 0.3 cm-3 on average. Thus the two best-fit models presented here were obtained with background media mimicking the classical warm ionised and warm neutral medium. We also found that 60Fe (which is condensed onto dust grains) can be delivered to Earth via two physical mechanisms: either through individual fast-paced supernova blast waves, which cross the Earth's orbit sometimes even twice as a result of reflection from the Local Bubble's outer shell, or, alternatively, through the supershell of the Local Bubble itself, injecting the 60Fe content of all previous supernovae at once, but over a longer time range.
Links to: arXiv, form interface, find, astro-ph, recent, 1704, contact, help (Access key information)
We quantify the gas-phase abundance of deuterium in cosmological zoom-in simulations from the Feedback In Realistic Environments project. The cosmic deuterium fraction decreases with time, because mass lost from stars is deuterium-free. At low metallicity, our simulations confirm that the deuterium abundance is very close to the primordial value. The deuterium abundance decreases towards higher metallicity, with very small scatter between the deuterium and oxygen abundance. We compare our simulations to existing high-redshift observations in order to determine a primordial deuterium fraction of (2.549 +/- 0.033) x 10^-5 and stress that future observations at higher metallicity can also be used to constrain this value. At fixed metallicity, the deuterium fraction decreases slightly with decreasing redshift, due to the increased importance of mass loss from intermediate-mass stars. We find that the evolution of the average deuterium fraction in a galaxy correlates with its star formation history. Our simulations are consistent with observations of the Milky Way's interstellar medium: the deuterium fraction at the solar circle is 83-92% of the primordial deuterium fraction. We use our simulations to make predictions for future observations. In particular, the deuterium abundance is lower at smaller galactocentric radii and in higher mass galaxies, showing that stellar mass loss is more important for fuelling star formation in these regimes (and can even dominate). Gas accreting onto galaxies has a deuterium fraction above that of the galaxies' interstellar medium, but below the primordial fraction, because it is a mix of gas accreting from the intergalactic medium and gas previously ejected or stripped from galaxies.
We characterize the contribution from accreted material to the galactic discs of the Auriga Project, a set of high resolution magnetohydrodynamic cosmological simulations of late-type galaxies performed with the moving-mesh code AREPO. Our goal is to explore whether a significant accreted (or ex-situ) stellar component in the Milky Way disc could be hidden within the near-circular orbit population, which is strongly dominated by stars born in-situ. One third of our models shows a significant ex-situ disc but this fraction would be larger if constraints on orbital circularity were relaxed. Most of the ex-situ material ($\gtrsim 50\%$) comes from single massive satellites ($> 6 \times 10^{10}~M_{\odot}$). These satellites are accreted with a wide range of infall times and inclination angles (up to $85^{\circ}$). Ex-situ discs are thicker, older and more metal-poor than their in-situ counterparts. They show a flat median age profile, which differs from the negative gradient observed in the in-situ component. As a result, the likelihood of identifying an ex-situ disc in samples of old stars on near-circular orbits increases towards the outskirts of the disc. We show three examples that, in addition to ex-situ discs, have a strongly rotating dark matter component. Interestingly, two of these ex-situ stellar discs show an orbital circularity distribution that is consistent with that of the in-situ disc. Thus, they would not be detected in typical kinematic studies.
The Infrared Dark Cloud (IRDC) G028.23-00.19 hosts a massive (1,500 Msun), cold (12 K), and 3.6-70 um IR dark clump (MM1) that has the potential to form high-mass stars. We observed this prestellar clump candidate with the SMA (~3.5" resolution) and JVLA (~2.1" resolution) in order to characterize the early stages of high-mass star formation and to constrain theoretical models. Dust emission at 1.3 mm wavelength reveals 5 cores with masses <15 Msun. None of the cores currently have the mass reservoir to form a high-mass star in the prestellar phase. If the MM1 clump will ultimately form high-mass stars, its embedded cores must gather a significant amount of additional mass over time. No molecular outflows are detected in the CO (2-1) and SiO (5-4) transitions, suggesting that the SMA cores are starless. By using the NH3 (1,1) line, the velocity dispersion of the gas is determined to be transonic or mildly supersonic (DeltaV_nt}/DeltaV_th ~1.1-1.8). The cores are not highly supersonic as some theories of high-mass star formation predict. The embedded cores are 4 to 7 times more massive than the clump thermal Jeans mass and the most massive core (SMA1) is 9 times less massive than the clump turbulent Jeans mass. These values indicate that neither thermal pressure nor turbulent pressure dominates the fragmentation of MM1. The low virial parameters of the cores (0.1-0.5) suggest that they are not in virial equilibrium, unless strong magnetic fields of ~1-2 mG are present. We discuss high-mass star formation scenarios in a context based on IRDC G028.23-00.19, a study case believed to represent the initial fragmentation of molecular clouds that will form high-mass stars.
We present $Suzaku$ off-center observations of two poor galaxy groups, NGC 3402 and NGC 5129, with temperatures below 1 keV. Through spectral decomposition, we measure their surface brightnesses and temperatures out to 330 and 680 times the critical density of the universe for NGC 3402 and NGC 5129, respectively. These quantities are consistent with extrapolations from existing inner measurements of the two groups. With the refined X-ray luminosities, both groups prefer $L_X-T$ relations without a break in the group regime. Furthermore, we measure the electron number densities and hydrostatic masses at these radii. We find that the electron number density profiles require three $\beta$ model components, with nearly flat slopes in the 3$^{rd}$ $\beta$ component for both groups. However, we find the effective slope in the outskirts to be $\beta_{out}$ = 0.59 and 0.49 for NGC 3402 and NGC 5129, respectively. Adding the gas mass measured from the X-ray data and stellar mass from group galaxy members, we measure baryon fractions of $f_b$ = 0.113 $\pm$ 0.013 and 0.091 $\pm$ 0.006 for NGC 3402 and NGC 5129, respectively. Combining other poor groups with well measured X-ray emission to the outskirts, we find an average baryon fraction of $f_{b,ave}$ = 0.100 $\pm$ 0.004 for X-ray bright groups with temperatures between 0.8$-$1.3 keV, extending existing constraints to lower mass systems.
We used the Gemini Multi-Object Spectrograph Integral Field Unit to map the gas distribution, excitation and kinematics within the inner kiloparsec of four nearby low-luminosity active galaxies: NGC3982, NGC4501, NGC2787 and NGC4450. The observations cover the spectral range 5600-7000{\AA} at a velocity resolution of 120km/s and spatial resolution ranging from 50 to 70pc at the galaxies. Extended emission in H{\alpha}, [NII]{\lambda}{\lambda}6548,6583, [SII]{\lambda}{\lambda}6716,6730 over most of the field-of-view is observed for all galaxies, while only NGC3982 shows [OI]{\lambda}6300 extended emission. The H{\alpha} equivalent widths combined with the [NII]/H{\alpha} line ratios reveal that NGC3982 and NGC4450 harbor Seyfert nuclei surrounded by regions with LINER excitation, while NGC2787 and NGC4501 harbor LINER nuclei. NGC3982 shows a partial ring of recent star-formation at 500pc from the nucleus, while in NGC4501 a region at 500pc west of the nucleus shows LINER excitation but has been interpreted as an aging HII region with the gas excitation dominated by shocks from supernovae. The line-of-sight velocity field of the gas shows a rotation pattern for all galaxies, with deviations from pure disk rotation observed in NGC3982, NGC4501 and NGC4450. For NGC4501 and NGC4450, many of these deviations are spatially coincident with dust structures seen in optical continuum images, leading to the interpretation that the deviations are due to shocks in the gas traced by the dust. A speculation is that these shocks lead to loss of angular momentum, allowing the gas to be transferred inwards to feed the AGN. In the case of NGC2787, instead of deviations in the rotation field, we see a misalignment of 40{^\circ} between the orientation of the line of nodes of the gas rotation and the photometric major axis of the galaxy. Evidence of compact nuclear outflows are seen in NGC4501 and NGC4450.
We study winds in 12 X-ray AGN host galaxies at z ~ 1. We find, using the low-ionization Fe II 2586 absorption in the stacked spectra, that the probability distribution function (PDF) of the centroid velocity shift in AGN has a median, 16th and 84th percentiles of (-87, -251, +86) km/s respectively. The PDF of the velocity dispersion in AGN has a median, 84th and 16th percentile of (139, 253, 52) km/s respectively. The centroid velocity and the velocity dispersions are obtained from a two component (ISM+wind) absorption line model. The equivalent width PDF of the outflow in AGN has median, 84th and 16th percentiles of (0.4, 0.8, 0.1) Angstrom. There is a strong ISM component in Fe II 2586 absorption with (1.2, 1.5, 0.8) Angstrom, implying presence of substantial amount cold gas in the host galaxies. For comparison, star-forming and X-ray undetected galaxies at a similar redshift, matched roughly in stellar mass and galaxy inclination, have a centroid velocity PDF with percentiles of (-74, -258, +90) km/s, and a velocity dispersion PDF percentiles of (150, 259, 57) km/s. Thus, winds in the AGN are similar to star-formation-driven winds, and are too weak to escape and expel substantial cool gas from galaxies. Our sample doubles the previous sample of AGN studied at z ~ 0.5 and extends the analysis to z ~ 1. A joint reanalysis of the z ~ 0.5 AGN sample and our sample yields consistent results to the measurements above.
We report new spectroscopic observations obtained with the Michigan/Magellan Fiber System of 308 red giants (RGs) located in two fields near the photometric center of the bar of the Large Magellanic Cloud. This sample consists of 131 stars observed in previous studies (in one field) and 177 newly-observed stars (in the second field) selected specifically to more reliably establish the metallicity and age distributions of the bar. For each star, we measure its heliocentric line-of-sight velocity, surface gravity and metallicity from its high-resolution spectrum (effective temperatures come from photometric colors). The spectroscopic Hertzsprung-Russell diagrams---modulo small offsets in surface gravities---reveal good agreement with model isochrones. The mean metallicity of the 177-RG sample is $\rm [Fe/H]=-0.76\pm0.02$ with a metallicity dispersion $\sigma=0.28\pm0.03$. The corresponding metallicity distribution---corrected for selection effects---is well fitted by two Gaussian components: one metal-rich with a mean $-0.66\pm0.02$ and a standard deviation $0.17\pm0.01$, and the other metal-poor with $-1.20\pm0.24$ and $0.41\pm0.06$. The metal-rich and metal-poor populations contain approximately 85% and 15% of stars, respectively. We also confirm the velocity dispersion in the bar center decreases significantly from $31.2\pm4.3$ to $18.7\pm1.9$ km s$^{-1}$ with increasing metallicity over the range $-2.09$ to $-0.38$. Individual stellar masses are estimated using the spectroscopic surface gravities and the known luminosities. We find that lower mass hence older RGs have larger metallicity dispersion and lower mean metallicity than the higher-mass, younger RGs. The estimated masses, however, extend to implausibly low values ($\rm \sim 0.1~M_{\odot}$) making it impossible to obtain an absolute age-metallicity or age distribution of the bar.
We study the physical and dynamical properties of the ionized gas in the prototypical HII galaxy Henize 2-10 using MUSE integral field spectroscopy. The large scale dynamics is dominated by extended outflowing bubbles, probably the results of massive gas ejection from the central star forming regions. We derive a mass outflow rate dMout/dt~0.30 Msun/yr, corresponding to mass loading factor eta~0.4, in range with similar measurements in local LIRGs. Such a massive outflow has a total kinetic energy that is sustainable by the stellar winds and Supernova Remnants expected in the galaxy. We use classical emission line diagnostic to study the dust extinction, electron density and ionization conditions all across the galaxy, confirming the extreme nature of the highly star forming knots in the core of the galaxy, which show high density and high ionization parameter. We measure the gas phase metallicity in the galaxy taking into account the strong variation of the ionization parameter, finding that the external parts of the galaxy have abundances as low as 12 + log(O/H)~8.3, while the central star forming knots are highly enriched with super solar metallicity. We find no sign of AGN ionization in the galaxy, despite the recent claim of the presence of a super massive active Black Hole in the core of He~2-10. We therefore reanalyze the X-ray data that were used to propose the presence of the AGN, but we conclude that the observed X-ray emission can be better explained with sources of a different nature, such as a Supernova Remnant.
We estimated the ortho-{\rm{H$_2$O}} abundances of G267.9--1.1, G268.4--0.9, G333.1--0.4 and G336.5--1.5, four of the brightest ortho-{\rm{H$_2$O}} sources in the southern sky observed by the Submillimeter Wave Astronomy Satellite (ortho-{\rm{H$_2$O}} 1$_{10}$ -- 1$_{01}$ line, 556.936~GHz). The typical molecular clumps in our sample have H$_2$ column densities of $10 ^{22}$ to $10 ^{23}${\,}cm$^{-2}$ and ortho-{\rm{H$_2$O}} abundances of 10$^{-10}$. Compared with previous studies, the ortho-{\rm{H$_2$O}} abundances are at a low level, which can be caused by the low temperatures of these clumps. To estimate the ortho-{\rm{H$_2$O}} abundances, we used the CS $J = 2 \to 1$ line (97.98095~GHz) and CS $J = 5 \to 4$ (244.93556~GHz) line observed by{ the} Swedish-ESO 15\,m Submillimeter Telescope (SEST) to calculate the temperatures of the clumps and the 350~$\upmu$m dust continuum observed by{ the} Caltech Submillimeter Observatory (CSO) telescope to estimate the H$_2$ column densities. The observations of {\rm{N$_2$H$^+$}} ($J = 1 \to 0$) for these clumps were also acquired by SEST and the corresponding abundances were estimated. The {\rm{N$_2$H$^+$}} abundance in each clump shows a common decreasing trend toward the center and {a} typical abundance range from 10$^{-11}$ to 10$^{-9}$.
Recently, large samples of visually classified early-type galaxies (ETGs) containing dust have been identified using space-based infrared observations with the Herschel Space Telescope. The presence of large quantities of dust in massive ETGs is peculiar as X-ray halos of these galaxies are expected to destroy dust in 10 Myr (or less). This has sparked a debate regarding the origin of the dust: is it internally produced by asymptotic giant branch (AGB) stars, or is it accreted externally through mergers? We examine the 2D stellar and ionised gas kinematics of dusty ETGs using IFS observations from the SAMI galaxy survey, and integrated star-formation rates, stellar masses, and dust masses from the GAMA survey. Only 8% (4/49) of visually-classified ETGs are kinematically consistent with being dispersion-supported systems. These "dispersion-dominated galaxies" exhibit discrepancies between stellar and ionised gas kinematics, either offsets in the kinematic position angle or large differences in the rotational velocity, and are outliers in star-formation rate at a fixed dust mass compared to normal star-forming galaxies. These properties are suggestive of recent merger activity. The remaining 90% of dusty ETGs have low velocity dispersions and/or large circular velocities, typical of "rotation-dominated galaxies". These results, along with the general evidence of published works on X-ray emission in ETGs, suggest that they are unlikely to host hot, X-ray gas consistent with their low stellar mass when compared to dispersion-dominated galaxies. This means dust will be long lived and thus these galaxies do not require external scenarios for the origin of their dust content.
The observation of counter rotation in galaxies (i.e. gas that rotates in the opposite direction to the stellar component or two co-spatial stellar populations with opposite rotation) is becoming more commonplace with modern integral field spectroscopic surveys. In this paper we explore the emergence of counter-rotation (both stellar and gaseous) in S0 galaxies from smoothed-particle hydrodynamics simulations of 1/10 mass ratio minor mergers between a 10^10.8 solar mass disk galaxy with a bulge-to-total ration of 0.17 and a gas rich companion (gas-to-stellar mass fraction of 5.0). These simulations include a self-consistent treatment of gas dynamics, star formation, the production/destruction of H2 and dust, and the time evolution of the interstellar radiation field. We explore the effect of retrograde versus prograde obits, gas and bulge mass fractions of the primary galaxy, and orbital parameters of the companion. The key requirement for producing counter rotation in stars or gas in a merger remnant is a retrograde primary, while the relative spin of the companion affects only the radial extent of the accreted gas. We also find that including a significant amount of gas in the primary can prevent the emergence of counter-rotating gas, although accreted stars retain counter-rotation. Bulge mass and orbit have a secondary effect, generally influencing the final distribution of accreted stars and gas within the framework outlined above. In addition to our primary focus of counter-rotating components in galaxies, we also make some predictions regarding the SFRs, H2 distributions, and dust in minor merger remnants.
We examine the clustering of quasars over a wide luminosity range, by utilizing 901 quasars at $\overline{z}_{\rm phot}\sim3.8$ with $-24.73<M_{\rm 1450}<-22.23$ photometrically selected from the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) S16A Wide2 date release and 342 more luminous quasars at $3.4<z_{\rm spec}<4.6$ having $-28.0<M_{\rm 1450}<-23.95$ from the Sloan Digital Sky Survey (SDSS) that fall in the HSC survey fields. We measure the bias factors of two quasar samples by evaluating the cross-correlation functions (CCFs) between the quasar samples and 25790 bright $z\sim4$ Lyman Break Galaxies (LBGs) in $M_{\rm 1450}<-21.25$ photometrically selected from the HSC dataset. Over an angular scale of \timeform{10.0"} to \timeform{1000.0"}, the bias factors are $5.93^{+1.34}_{-1.43}$ and $2.73^{+2.44}_{-2.55}$ for the low and high luminosity quasars, respectively, indicating no luminosity dependence of quasar clustering at $z\sim4$. It is noted that the bias factor of the luminous quasars estimated by the CCF is smaller than that estimated by the auto-correlation function (ACF) over a similar redshift range, especially on scales below \timeform{40.0"}. Moreover, the bias factor of the less-luminous quasars implies the minimal mass of their host dark matter halos (DMHs) is $0.3$-$2\times10^{12}h^{-1}M_{\odot}$, corresponding to a quasar duty cycle of $0.001$-$0.06$.
We present the first results from the B-fields In STar-forming Region Observations (BISTRO) survey, using the Sub-millimetre Common-User Bolometer Array 2 (SCUBA-2) camera, with its associated polarimeter (POL-2), on the James Clerk Maxwell Telescope (JCMT) in Hawaii. We discuss the survey's aims and objectives. We describe the rationale behind the survey, and the questions which the survey will aim to answer. The most important of these is the role of magnetic fields in the star formation process on the scale of individual filaments and cores in dense regions. We describe the data acquisition and reduction processes for POL-2, demonstrating both repeatability and consistency with previous data. We present a first-look analysis of the first results from the BISTRO survey in the OMC 1 region. We see that the magnetic field lies approximately perpendicular to the famous 'integral filament' in the densest regions of that filament. Furthermore, we see an 'hour-glass' magnetic field morphology extending beyond the densest region of the integral filament into the less-dense surrounding material, and discuss possible causes for this. We also discuss the more complex morphology seen along the Orion Bar region. We examine the morphology of the field along the lower-density north-eastern filament. We find consistency with previous theoretical models that predict magnetic fields lying parallel to low-density, non-self-gravitating filaments, and perpendicular to higher-density, self-gravitating filaments.
We report high angular resolution (4.9" x 3.0") images of reactive ions SH+, HOC+, and SO+ toward the Orion Bar photodissociation region (PDR). We used ALMA-ACA to map several rotational lines at 0.8 mm, complemented with multi-line observations obtained with the IRAM 30m telescope. The SH+ and HOC+ emission is restricted to a narrow layer of 2"- to 10"-width (~800 to 4000 AU depending on the assumed PDR geometry) that follows the vibrationally excited H2^* emission. Both ions efficiently form very close to the H/H2 transition zone, at a depth of A_V < 1 mag into the neutral cloud, where abundant C+, S+, and H2^* coexist. SO+ peaks slightly deeper into the cloud. The observed ions have low rotational temperatures (T_rot~10-30 K << T_k) and narrow line-widths (~2-3 km/s), a factor of ~2 narrower that those of the lighter reactive ion CH+. This is consistent with the higher reactivity and faster radiative pumping rates of CH+ compared to the heavier ions, which are driven relatively faster toward smaller velocity dispersion by elastic collisions and toward lower T_rot by inelastic collisions. We estimate column densities and average physical conditions from a non-LTE excitation model (n(H2)~10^5-10^6 cm^-3, n(e^-)~10 cm^-3, and T_k~200 K). Regardless of the excitation details, SH+ and HOC+ clearly trace the most exposed layers of the UV-irradiated molecular cloud surface, whereas SO+ arises from slightly more shielded layers.
We present Rotation Measures (RM) of the diffuse Galactic synchrotron emission from the Canadian Galactic Plane Survey (CGPS) and compare them to RMs of extragalactic sources in order to study the large-scale reversal in the Galactic magnetic field (GMF). Using Stokes Q, U and I measurements of the Galactic disk collected with the Synthesis Telescope at the Dominion Radio Astrophysical Observatory, we calculate RMs over an extended region of the sky, focusing on the low longitude range of the CGPS (l=52deg to l=72deg). We note the similarity in the structures traced by the compact sources and the extended emission and highlight the presence of a gradient in the RM map across an approximately diagonal line, which we identify with the well-known field reversal of the Sagittarius-Carina arm. We suggest that the orientation of this reversal is a geometric effect resulting from our location within a GMF structure arising from current sheets that are not perpendicular to the Galactic plane, as is required for a strictly radial field reversal, but that have at least some component parallel to the disk. Examples of models that fit this description are the three-dimensional dynamo-based model of Gressel et al. (2013) and a Galactic scale Parker spiral (Akasofu & Hakamada 1982), although the latter may be problematic in terms of Galactic dynamics. We emphasize the importance of constructing three-dimensional models of the GMF to account for structures like the diagonal RM gradient observed in this dataset.
The nature of ultraluminous X-ray sources (ULXs) -- off-nuclear extra-galactic sources with luminosity, assumed isotropic, $\gtrsim 10^{39}$ erg s$^{-1}$ -- is still debated. One possibility is that ULXs are stellar black holes accreting beyond the Eddington limit. This view has been recently reinforced by the discovery of ultrafast outflows at $\sim 0.1$-$0.2c$ in the high resolution spectra of a handful of ULXs, as predicted by models of supercritical accretion discs. Under the assumption that ULXs are powered by super-Eddington accretion onto black holes, we use the properties of the observed outflows to self-consistently constrain their masses and accretion rates. We find masses $\lesssim 100$ M$_{\odot}$ and typical accretion rates $\sim 10^{-5}$ M$_{\odot}$ yr$^{-1}$, i.e. $\approx 10$ times larger than the Eddington limit calculated with a radiative efficiency of 0.1. However, the emitted luminosity is only $\approx 10\%$ beyond the Eddington luminosity, because most of the energy released in the inner part of the accretion disc is used to accelerate the wind, which implies radiative efficiency $\sim 0.01$. Our results are consistent with a formation model where ULXs are black hole remnants of massive stars evolved in low-metallicity environments.
The wide-area XMM-XXL X-ray survey is used to explore the fraction of obscured AGN at high accretion luminosities, $L_X (\rm 2-10 \, keV) > 10^{44} \, erg \,s ^{-1}$, and out to redshift $z\approx1.5$. The sample covers an area of about $\rm14\,deg^2$ and provides constraints on the space density of powerful AGN over a wide range of neutral hydrogen column densities extending beyond the Compton-thick limit, $\rm N_H\approx10^{24}\,cm^{-2}$. The fraction of obscured Compton-thin ($\rm N_H=10^{22}-10^{24}\,cm^{-2}$) AGN is estimated to be $\approx0.35$ for luminosities $L_X(\rm 2-10\,keV)>10^{44}\,erg\,s^{-1}$ independent of redshift. For less luminous sources the fraction of obscured Compton-thin AGN increases from $0.45\pm0.10$ at $z=0.25$ to $0.75\pm0.05$ at $z=1.25$. Studies that select AGN in the infrared via template fits to the observed Spectral Energy Distribution of extragalactic sources estimate space densities at high accretion luminosities consistent with the XMM-XXL constraints. There is no evidence for a large population of AGN (e.g. heavily obscured) identified in the infrared and missed at X-ray wavelengths. We further explore the mid-infrared colours of XMM-XXL AGN as a function of accretion luminosity, column density and redshift. The fraction of XMM-XXL sources that lie within the mid-infrared colour wedges defined in the literature to select AGN is primarily a function of redshift. This fraction increases from about 20-30% at z=0.25 to about 50-70% at $z=1.5$.
Fragmentation of filaments into dense cores is thought to be an important step in forming stars. The bar-mode instability of spherically collapsing cores found in previous linear analysis invokes a possibility of re-fragmentation of the cores due to their ellipsoidal (prolate or oblate) deformation. To investigate this possibility, here we perform three-dimensional self-gravitational hydrodynamics simulations that follow all the way from filament fragmentation to subsequent core collapse. We assume the gas is polytropic with index \gamma, which determines the stability of the bar-mode. For the case that the fragmentation of isolated hydrostatic filaments is triggered by the most unstable fragmentation mode, we find the bar mode grows as collapse proceeds if \gamma < 1.1, in agreement with the linear analysis. However, it takes more than ten orders-of-magnitude increase in the central density for the distortion to become non-linear. In addition to this fiducial case, we also study non-fiducial ones such as the fragmentation is triggered by a fragmentation mode with a longer wavelength and it occurs during radial collapse of filaments and find the distortion rapidly grows. In most of astrophysical applications, the effective polytropic index of collapsing gas exceeds 1.1 before ten orders-of-magnitude increase in the central density. Thus, supposing the fiducial case of filament fragmentation, re-fragmentation of dense cores would not be likely and their final mass would be determined when the filaments fragment.
A superbubble which advances in a symmetric Navarro--Frenk--White density profile or in an auto-gravitating density profile generates a thick shell with a radius that can reach 10 kpc. The application of the symmetric and asymmetric image theory to this thick 3D shell produces a ring in the 2D map of intensity and a characteristic `U' shape in the case of 1D cut of the intensity. A comparison of such a ring originating from a superbubble is made with the Einstein's ring. A Taylor approximation of order 10 for the angular diameter distance is derived in order to deal with high values of the redshift.
Links to: arXiv, form interface, find, astro-ph, recent, 1704, contact, help (Access key information)