Links to: arXiv, form interface, find, astro-ph, recent, 1709, contact, help (Access key information)
We study low-density axisymmetric accretion flows onto black holes (BHs) with two-dimensional hydrodynamical simulations, adopting the $\alpha$-viscosity prescription. When the gas angular momentum is low enough to form a rotationally supported disk within the Bondi radius ($R_{\rm B}$), we find a global steady accretion solution. The solution consists of a rotational equilibrium distribution at $r\sim R_{\rm B}$, where the density follows $\rho \propto (1+R_{\rm B}/r)^{3/2}$, surrounding a geometrically thick and optically thin accretion disk at the centrifugal radius, where thermal energy generated by viscosity is transported via strong convection. Physical properties of the inner solution agree with those expected in convection-dominated accretion flows (CDAF; $\rho \propto r^{-1/2}$). In the inner CDAF solution, the gas inflow rate decreases towards the center due to convection ($\dot{M}\propto r$), and the net accretion rate (including both inflows and outflows) is strongly suppressed by several orders of magnitude from the Bondi accretion rate $\dot{M}_{\rm B}$ The net accretion rate depends on the viscous strength, following $\dot{M}/\dot{M}_{\rm B}\propto (\alpha/0.01)^{0.6}$. This solution holds for low accretion rates of $\dot{M}_{\rm B}/\dot{M}_{\rm Edd}< 10^{-3}$ having minimal radiation cooling, where $\dot{M}_{\rm Edd}$ is the Eddington rate. In a hot plasma at the bottom ($r<10^{-3}~R_{\rm B}$), thermal conduction would dominate the convective energy flux. Since suppression of the accretion by convection ceases, the final BH feeding rate is found to be $\dot{M}/\dot{M}_{\rm B} \sim 10^{-3}-10^{-2}$. This rate is as low as $\dot{M}/\dot{M}_{\rm Edd} \sim 10^{-7}-10^{-6}$ inferred for SgrA$^*$ and the nuclear BHs in M31 and M87, and can explain the low luminosities in these sources, without invoking any feedback mechanism.
Most giant spiral galaxies have pseudo or disk-like bulges that are considered to be the result of purely secular processes. This may challenge the hierarchical scenario predicting about one major merger per massive galaxy ($>$$3\times 10^{10} M_{\odot}$) since the last $\sim$ 9 billion years. Here we verify whether or not the association between pseudo-bulges and secular processes is irrevocable. Using GADGET2 N-body/SPH simulations, we have conducted a systematic study of remnants of major mergers which progenitors have been selected (1) to follow the gas richness-look back time relationship, and (2) with a representative distribution of orbits and spins in a cosmological frame. Analyzing the surface-mass density profile of both nearby galaxies and merger remnants with two components, we find that most of them show pseudo-bulges or bar dominated centers. Even if some orbits lead to classical bulges just after the fusion, the contamination by the additional gas that gradually accumulates to the center and forming stars later on, leads to remnants apparently dominated by pseudo-bulges. We also found that simple SPH simulations should be sufficient to form realistic spiral galaxies as remnants of ancient gas-rich mergers without need for specifically tuned feedback conditions. We then conclude that pseudo-bulges and bars in spiral galaxies are natural consequences of major mergers when they are realized in a cosmological context, i.e., with gas-rich progenitors as expected when selected in the distant Universe.
We examine the z = 0 group-integrated stellar and cold baryonic (stars + cold atomic gas) mass functions (group SMF and CBMF) and the baryonic collapse efficiency (group cold baryonic to dark matter halo mass ratio) using the RESOLVE and ECO survey galaxy group catalogs and a galform semi-analytic model (SAM) mock catalog. The group SMF and CBMF fall off more steeply at high masses and rise with a shallower low-mass slope than the theoretical halo mass function (HMF). The transition occurs at group-integrated cold baryonic mass M_coldbary ~ 10^11 Msun. The SAM, however, has significantly fewer groups at the transition mass ~ 10^11 Msun and a steeper low-mass slope than the data, suggesting that feedback is too weak in low-mass halos and conversely too strong near the transition mass. Using literature prescriptions to include hot halo gas and potential unobservable galaxy gas produces a group BMF with slope similar to the HMF even below the transition mass. Its normalization is lower by a factor of ~2, in agreement with estimates of warm-hot gas making up the remaining difference. We compute baryonic collapse efficiency with the halo mass calculated two ways, via halo abundance matching (HAM) and via dynamics (extended all the way to three-galaxy groups using stacking). Using HAM, we find that baryonic collapse efficiencies reach a flat maximum for groups across the halo mass range of M_halo ~ 10^11.4-12 Msun, which we label "nascent groups." Using dynamics, however, we find greater scatter in baryonic collapse efficiencies, likely indicating variation in group hot-to-cold baryon ratios. Similarly, we see higher scatter in baryonic collapse efficiencies in the SAM when using its true groups and their group halo masses as opposed to friends-of-friends groups and HAM masses.
We present VLA images and HST/STIS spectra of sources within the center of the brightest cluster galaxy (BCG) in Abell 2261. These observations were obtained to test the hypothesis that its extremely large, flat core reflects the ejection of its supermassive black hole. Spectra of three of the four most luminous "knots" embedded in the core were taken to test whether one may represent stars bound to a displaced massive black hole. The three knots have radial velocity offsets dV < ~150 km/s from the BCG. Knots 2 and 3 show kinematics, colors, and stellar masses consistent with infalling low-mass galaxies or larger stripped cluster members. Large errors in the stellar velocity dispersion of Knot 1, however, mean that we cannot rule out the hypothesis that it hosts a high-mass black hole. A2261-BCG has a compact, relic radio-source offset by 6.5 kpc (projected) from the optical core's center, but no active radio core that would pinpoint the galaxy's central black hole to a tight 10 GHz flux limit <3.6 uJy. Its spectrum and morphology are suggestive of an active galactic nucleus that switched off >48 Myr ago, with an equipartition condition magnetic field of 15 uG. These observations are still consistent with the hypothesis that the nuclear black hole has been ejected from its core, but the critical task of locating the supermassive black hole or demonstrating that A2261-BCG lacks one remains to be done.
The physics of star formation at its extreme, in the nuclei of the largest and densest star clusters in the universe - potential massive black hole nurseries - has for decades eluded scrutiny. Spectroscopy of these systems has been scarce, whereas theoretical claims have been made that radiation pressure on dust grains somehow inhibits star formation. Here, we harness an accelerated Monte Carlo radiation transport scheme to report a radiation hydrodynamical simulation of super star cluster formation in turbulent clouds. We find that radiation pressure reduces the global star formation efficiency by 30-35%, and the star formation rate by 15-50%, both relative to a radiation-free control run. Overall, radiation pressure is ineffective in limiting the gas supply for star formation and the final stellar mass of the most massive cluster is $\sim 1.3\times10^6\,M_\odot$. The limited impact as compared to that implied by idealized theoretical models is attributed to a radiation-matter anti-correlation in the supersonically turbulent, gravitationally collapsing medium. In isolated regions outside massive clusters, where the gas distribution is less disturbed, radiation pressure is more effective in limiting star formation. The resulting stellar density at the cluster core is $\ge 10^{8}\,M_{\odot}\,\textrm{pc}^{-3}$, with stellar velocity dispersion $\gtrsim 70\,\text{km}\,\text{s}^{-1}$. We conclude that the super star cluster nucleus is propitious to the formation of very massive stars via dynamical core collapse and stellar merging. We speculate that the very massive star may avoid the claimed catastrophic mass loss by continuing to accrete dense gas condensing from a gravitationally-confined ionized phase.
We explore the circumgalactic metal content traced by commonly observed low ion absorbers, including C II, Si II, Si III, Si IV, and Mg II. We use a set of cosmological hydrodynamical zoom simulations run with the EAGLE model and including a non-equilibrium ionization and cooling module that follows 136 ions. The simulations of z~0.2 L* (M_200=10^11.7-10^12.3 Msol) haloes hosting star-forming galaxies and group-sized (M_200=10^12.7-10^13.3 Msol) haloes hosting mainly passive galaxies reproduce key trends observed by the COS-Halos survey-- low ion column densities show 1) little dependence on galaxy specific star formation rate, 2) a patchy covering fraction indicative of 10^4 K clouds with a small volume filling factor, and 3) a declining covering fraction as impact parameter increases from 20-160 kpc. Simulated Si II, Si III, Si IV, C II, and C III column densities show good agreement with observations, while Mg II is under-predicted. Low ions trace a significant metal reservoir, ~10^8 Msol, residing primarily at 10-100 kpc from star-forming and passive central galaxies. These clouds tend to flow inwards and most will accrete onto the central galaxy within the next several Gyr, while a small fraction are entrained in strong outflows. A two-phase structure describes the inner CGM (<0.5 R_200) with low-ion metal clouds surrounded by a hot, ambient medium. This cool phase is separate from the O VI observed by COS-Halos, which arises from the outer CGM (>0.5 R_200) tracing virial temperature gas around L* galaxies. Physical parameters derived from standard photo-ionization modelling of observed column densities (e.g. aligned Si II/Si III absorbers) are validated against our simulations. Our simulations therefore support previous ionization models indicating that cloud covering factors decline while densities and pressures show little variation with increasing impact parameter.
The cosmological origin of carbon, the fourth most abundant element in the Universe, is not well known and matter of heavy debate. We investigate the behavior of C/O to O/H in order to constrain the production mechanism of carbon. We measured emission-line intensities in a spectral range from 1600 to 10000 \AA\ on Space Telescope Imaging Spectrograph (STIS) long-slit spectra of 18 starburst galaxies in the local Universe. We determined chemical abundances through traditional nebular analysis and we used a Markov Chain Monte Carlo (MCMC) method to determine where our carbon and oxygen abundances lie in the parameter space. We conclude that our C and O abundance measurements are sensible. We analyzed the behavior of our sample in the [C/O] vs. [O/H] diagram with respect to other objects such as DLAs, neutral ISM measurements, and disk and halo stars, finding that each type of object seems to be located in a specific region of the diagram. Our sample shows a steeper C/O vs. O/H slope with respect to other samples, suggesting that massive stars contribute more to the production of C than N at higher metallicities, only for objects where massive stars are numerous; otherwise intermediate-mass stars dominate the C and N production.
We use observations made with the Giant Metrewave Radio Telescope (GMRT) to probe the neutral hydrogen (HI) gas content of field galaxies in the VIMOS VLT Deep Survey (VVDS) 14h field at $z \approx 0.32$. Because the HI emission from individual galaxies is too faint to detect at this redshift, we use an HI spectral stacking technique using the known optical positions and redshifts of the 165 galaxies in our sample to co-add their HI spectra and thus obtain the average HI mass of the galaxies. Stacked HI measurements of 165 galaxies show that 95 per cent of the neutral gas is found in blue, star-forming galaxies. Among these galaxies, those having lower stellar mass are more gas-rich than more massive ones. We apply a volume correction to our HI measurement to evaluate the HI gas density at $z \approx 0.32$ as $\Omega_{HI}=(0.50\pm0.18) \times 10^{-3}$ in units of the cosmic critical density. This value is in good agreement with previous results at z < 0.4, suggesting no evolution in the neutral hydrogen gas density over the last $\sim 4$ Gyr. However the $z \approx 0.32$ gas density is lower than that at $z \sim 5$ by at least a factor of two.
Abridged. Understanding the detailed structure of the interstellar gas is essential for our knowledge of the star formation process. The small-scale structure of the interstellar medium (ISM) is a direct consequence of the galactic scales and making the link between the two is essential. We perform adaptive mesh simulations that aim to bridge the gap between the intermediate galactic scales and the self-gravitating prestellar cores. For this purpose we use stratified supernova regulated ISM magneto-hydrodynamical (MHD) simulations at the kpc scale to set up the initial conditions. We then zoom, performing a series of concentric uniform refinement and then refining on the Jeans length for the last levels. This allows us to reach a spatial resolution of a few $10^{-3}$ pc. The cores are identified using a clump finder and various criteria based on virial analysis. Their most relevant properties are computed and, due to the large number of objects formed in the simulations, reliable statistics are obtained. The cores properties show encouraging agreements with observations. The mass spectrum presents a clear powerlaw at high masses with an exponent close to $\simeq -1.3$ and a peak at about 1-2 $M_\odot$. The velocity dispersion and the angular momentum distributions are respectively a few times the local sound speed and a few $10^{-2}$ pc km s$^{-1}$. We also find that the distribution of thermally supercritical cores present a range of magnetic mass-to-flux over critical mass-to-flux ratio which typically ranges between $\simeq$0.3 and 3.
We mapped the kinetic temperature structure of the Orion molecular cloud 1 with para-H2CO(303-202, 322-221, and 321-220) using the APEX 12m telescope. This is compared with the temperatures derived from the ratio of the NH3(2,2)/(1,1) inversion lines and the dust emission. Using the RADEX non-LTE model, we derive the gas kinetic temperature modeling the measured averaged line ratios of para-H2CO 322-221/303-202 and 321-220/303-202. The gas kinetic temperatures derived from the para-H2CO line ratios are warm, ranging from 30 to >200 K with an average of 62 K at a spatial density of 10$^5$ cm$^{-3}$. These temperatures are higher than those obtained from NH3(2,2)/(1,1) and CH3CCH(6-5) in the OMC-1 region. The gas kinetic temperatures derived from para-H2CO agree with those obtained from warm dust components measured in the mid infrared (MIR), which indicates that the para-H2CO(3-2) ratios trace dense and warm gas. The cold dust components measured in the far infrared (FIR) are consistent with those measured with NH3(2,2)/(1,1) and the CH3CCH(6-5) line series. With dust at MIR wavelengths and para-H2CO(3-2) on one side and dust at FIR wavelengths, NH3(2,2)/(1,1), and CH3CCH(6-5) on the other, dust and gas temperatures appear to be equivalent in the dense gas of the OMC-1 region, but provide a bimodal distribution, one more directly related to star formation than the other. The non-thermal velocity dispersions of para-H2CO are positively correlated with the gas kinetic temperatures in regions of strong non-thermal motion (Mach number >2.5) of the OMC-1, implying that the higher temperature traced by para-H2CO is related to turbulence on a 0.06 pc scale. Combining the temperature measurements with para-H2CO and NH3(2,2)/(1,1) line ratios, we find direct evidence for the dense gas along the northern part of the OMC-1 10 km s$^{-1}$ filament heated by radiation from the central Orion nebula.
In an earlier paper we modeled the far-infrared emission from a star-forming galaxy using the photoionisation code CLOUDY and presented metallicity sensitive diagnostics based on far-infrared fine structure line ratios. Here, we focus on the applicability of the [OIII]88/[NII]122 microns line ratio as a gas phase metallicity indicator in high redshift submillimetre luminous galaxies. The [OIII]88/[NII]122 microns ratio is strongly dependent on the ionization parameter (which is related to the total number of ionizing photons) as well as the gas electron density. We demonstrate how the ratio of 88/$122 continuum flux measurements can provide a reasonable estimate of the ionization parameter while the availability of the [NII]205 microns line can constrain the electron density. Using the [OIII]88/[NII]122 microns line ratios from a sample of nearby normal and star-forming galaxies we measure their gas phase metallicities and find that their mass metallicity relation is consistent with the one derived using optical emission lines. Using new, previously unpublished, Herschel spectroscopic observations of key far-infrared fine structure lines of the z~3 galaxy HLSW-01 and additional published measurements of far-infrared fine structure lines of high-z submillimetre luminous galaxies we derive gas phase metallicities using their [OIII]88/[NII]122 microns line ratio. We find that the metallicities of these z~3 submm luminous galaxies are consistent with solar metallicities and that they appear to follow the mass-metallicity relation expected for z~3 systems.
The origin of the observed morphological and kinematic substructure of young star forming regions is a matter of debate. We offer a new analysis of data from simulations of globally gravitationally collapsing clouds of progenitor gas to answer questions about sub-structured star formation in the context of cold collapse. As a specific example, we compare our models to recent radial velocity survey data from the IN-SYNC survey of Orion and new observations of dense gas kinematics, and offer possible interpretations of kinematic and morphological signatures in the region. In the context of our model, we find the frequently-observed hub-filament morphology of the gas naturally arises during gravitational evolution, as well as the dynamically-distinct kinematic substructure of stars. We emphasize that the global and not just the local gravitational potential plays an important role in determining the dynamics of both clusters and filaments.
Cosmic Rays (CR) process the matter of the Interstellar Medium. Such
energetic processing not only modifies the interstellar matter but also injects
chemical species in the gas phase. In this work, we study the effect of the CR
on the astrophysical polycyclic aromatic hydrocarbons (PAH). For events in
which many electrons are stripped out from the PAH by interaction with a heavy
cosmic ray particle, coulomb explosion takes place and carbon chains are
produced. The fragments production rates of carbon chains are of particular
interest for astrophysical models. We computed PAH multi-ionization cross
sections with an Independent Atom and Electron collisional model. We introduced
and used a model to predict the fragmentation pattern for the coulomb
explosion. Experimental measurements on small hydrocarbons, C$_{60}$ and PAHs
were used to set confidence intervals on the calculations results.
The carbon chains production rates were calculated using different CR fluxes
and elemental compositions, to account for the variations expected in various
astrophysical environments. A range of PAH sizes and compactness were also
explored. The PAH lifetime with respect to a standard interstellar CR flux
(corresponding to an H$_2$ ionization rate of $\zeta \approx$
6.10$^{-17}$s$^{-1}$) is found to be in the order of a few billion years. The
production rates of interstellar carbon chains containing around 5-15 carbon
atoms are in the order of few to many tens of percent of the H$_2$ ionization
rate $\zeta$. The exact rate value relies on the nature of the PAH and on the
CR composition. In diffuse medium, with ten percent of the available cosmic
carbon locked in PAHs, this process leads to carbon chain fractional abundances
at steady state, in the range of $10^{-15}$-$10^{-14}$, with a confidence
interval of about one order of magnitude. It reaches $10^{-13}$ in quiescent
dense clouds.
FU Orionis-type objects are young, low-mass stars with large outbursts in visible light that last for several years or decades. They are thought to represent an evolutionary phase during the life of every young star when accretion from the circumstellar disk is enhanced during recurring time periods. These outbursts are able to rapidly build up the star while affecting the circumstellar disk and thus the ongoing or future planet formation. In many models infall from a circumstellar envelope seems to be necessary to trigger the outbursts. We observed the J=1$-$0 rotational transition of $^{13}$CO and C$^{18}$O towards eight northern FU Orionis-type stars (V1057 Cyg, V1515 Cyg, V2492 Cyg, V2493 Cyg, V1735 Cyg, V733 Cep, RNO 1B and RNO 1C) and derive temperatures and envelope masses and discuss the morphology and kinematics of the circumstellar material. We detected extended CO emission associated with all our targets. Smaller scale CO clumps were found to be associated with five objects with radii of 2000$-$5000 AU and masses of 0.02$-$0.5 $M_{\odot}$; these are clearly heated by the central stars. Three of these envelopes are also strongly detected in the 2.7 mm continuum. No central CO clumps were detected around V733 Cep and V710 Cas but there are many other clumps in their environments. Traces of outflow activity were observed towards V1735 Cyg, V733 Cep and V710 Cas. The diversity of the observed envelopes enables us to set up an evolutionary sequence between the objects. We find their evolutionary state to range from early, embedded Class I stage to late, Class II-type objects with very low-mass circumstellar material. The results reinforce the idea of FU Orionis-type stars as representatives of a transitory stage between embedded Class I young stellar objects and classical T-Tauri stars.
The formation epoch of protostellar disks is debated because of the competing roles of rotation, turbulence, and magnetic fields in the early stages of low-mass star formation. Magnetohydrodynamics simulations of collapsing cores predict that rotationally supported disks may form in strongly magnetized cores through ambipolar diffusion or misalignment between the rotation axis and the magnetic field orientation. Detailed studies of individual sources are needed to cross check the theoretical predictions. We present 0.06-0.1" resolution images at 350 GHz toward B1b-N and B1b-S, which are young class 0 protostars, possibly first hydrostatic cores. The images have been obtained with ALMA, and we compare these data with magnetohydrodynamics simulations of a collapsing turbulent and magnetized core. The submillimeter continuum emission is spatially resolved by ALMA. Compact structures with optically thick 350 GHz emission are detected toward both B1b-N and B1b-S, with 0.2 and 0.35" radii (46 and 80 au at the Perseus distance of 230 pc), within a more extended envelope. The flux ratio between the compact structure and the envelope is lower in B1b-N than in B1b-S, in agreement with its earlier evolutionary status. The size and orientation of the compact structure are consistent with 0.2" resolution 32 GHz observations obtained with the Very Large Array as a part of the VANDAM survey, suggesting that grains have grown through coagulation. The morphology, temperature, and densities of the compact structures are consistent with those of disks formed in numerical simulations of collapsing cores. Moreover, the properties of B1b-N are consistent with those of a very young protostar, possibly a first hydrostatic core. These observations provide support for the early formation of disks around low-mass protostars.
We study the asymmetry in the two-point cross-correlation function of two populations of galaxies focusing in particular on the relativistic effects that include the gravitational redshift. We derive the cross-correlation function on small and large scales using two different approaches: General Relativistic and Newtonian perturbation theory. Following recent work by Bonvin et al., Gaztanaga et al. and Croft, we calculate the dipole and the shell estimator with the two procedures and we compare our results. We find that while General Relativistic Perturbation Theory (GRPT) is able to make predictions of relativistic effects on very large, obviously linear scales (r > 50 Mpc/h), the presence of non-linearities physically occurring on much smaller scales (down to those describing galactic potential wells) can strongly affect the asymmetry estimators. These can lead to cancellations of the relativistic terms, and sign changes in the estimators on scales up to r ~ 50 Mpc/h. On the other hand, with an appropriate non-linear gravitational potential, the results obtained using Newtonian theory can successfully describe the asymmetry on smaller, non-linear scales (r < 20 Mpc/h) where gravitational redshift is the dominant term. On larger scales the asymmetry is much smaller in magnitude, and measurement is not within reach of current observations. This is in agreement with the observational results obtained by Gaztnaga et al. and the first detection of relativistic effects (on (r < 20 Mpc/h) scales) by Alam et al.
General relativistic effects have long been predicted to subtly influence the observed large-scale structure of the universe. The current generation of galaxy redshift surveys have reached a size where detection of such effects is becoming feasible. In this paper, we report the first detection of the redshift asymmetry from the cross-correlation function of two galaxy populations which is consistent with relativistic effects. The dataset is taken from the Sloan Digital Sky Survey DR12 CMASS galaxy sample, and we detect the asymmetry at the $2.7\sigma$ level by applying a shell-averaged estimator to the cross-correlation function. Our measurement dominates at scales around $10$ h$^{-1}$Mpc, larger than those over which the gravitational redshift profile has been recently measured in galaxy clusters, but smaller than scales for which linear perturbation theory is likely to be accurate. The detection significance varies by 0.5$\sigma$ with the details of our measurement and tests for systematic effects. We have also devised two null tests to check for various survey systematics and show that both results are consistent with the null hypothesis. We measure the dipole moment of the cross-correlation function, and from this the asymmetry is also detected, at the $2.8 \sigma$ level. The amplitude and scale-dependence of the clustering asymmetries are approximately consistent with the expectations of General Relativity and a biased galaxy population, within large uncertainties. We explore theoretical predictions using numerical simulations in a companion paper.
In a galaxy redshift survey the objects to be targeted for spectra are selected from a photometrically observed sample. The observed magnitudes and colours of galaxies in this parent sample will be affected by their peculiar velocities, through relativistic Doppler and relativistic beaming effects. In this paper we compute the resulting expected changes in galaxy photometry. The magnitudes of the relativistic effects are a function of redshift, stellar mass, galaxy velocity and velocity direction. We focus on the CMASS sample from the Sloan Digital Sky Survey (SDSS), Baryon Oscillation Spectroscopic Survey (BOSS), which is selected on the basis of colour and magnitude. We find that 0.10\% of the sample ($\sim 585$ galaxies) has been scattered into the targeted region of colour-magnitude space by relativistic effects, and conversely 0.09\% of the sample ($\sim 532$ galaxies) has been scattered out. Observational consequences of these effects include an asymmetry in clustering statistics, which we explore in a companion paper. Here we compute a set of weights which can be used to remove the effect of modulations introduced into the density field inferred from a galaxy sample. We conclude by investigating the possible effects of these relativistic modulation on large scale clustering of the galaxy sample.
Large redshift surveys of galaxies and clusters are providing the first opportunities to search for distortions in the observed pattern of large-scale structure due to such effects as gravitational redshift. We focus on non-linear scales and apply a quasi-Newtonian approach using N-body simulations to predict the small asymmetries in the cross-correlation function of two galaxy different populations. Following recent work by Bonvin et al., Zhao and Peacock and Kaiser on galaxy clusters, we include effects which enter at the same order as gravitational redshift: the transverse Doppler effect, light-cone effects, relativistic beaming, luminosity distance perturbation and wide-angle effects. We find that all these effects cause asymmetries in the cross-correlation functions. Quantifying these asymmetries, we find that the total effect is dominated by the gravitational redshift and luminosity distance perturbation at small and large scales, respectively. By adding additional subresolution modelling of galaxy structure to the large-scale structure information, we find that the signal is significantly increased, indicating that structure on the smallest scales is important and should be included. We report on comparison of our simulation results with measurements from the SDSS/BOSS galaxy redshift survey in a companion paper.
Links to: arXiv, form interface, find, astro-ph, recent, 1709, contact, help (Access key information)
RGG 118 (SDSS 1523+1145) is a nearby ($z=0.0243$), dwarf disk galaxy ($M_{\ast}\approx2\times10^{9} M_{\odot}$) found to host an active $\sim50,000$ solar mass black hole at its core (Baldassare et al. 2015). RGG 118 is one of a growing collective sample of dwarf galaxies known to contain active galactic nuclei -- a group which, until recently, contained only a handful of objects. Here, we report on new \textit{Hubble Space Telescope} Wide Field Camera 3 UVIS and IR imaging of RGG 118, with the main goal of analyzing its structure. Using 2-D parametric modeling, we find that the morphology of RGG 118 is best described by an outer spiral disk, inner component consistent with a pseudobulge, and central PSF. The luminosity of the PSF is consistent with the central point source being dominated by the AGN. We measure the luminosity and mass of the "pseudobulge" and confirm that the central black hole in RGG 118 is under-massive with respect to the $M_{BH}-M_{\rm bulge}$ and $M_{BH}-L_{\rm bulge}$ relations. This result is consistent with a picture in which black holes in disk-dominated galaxies grow primarily through secular processes.
In this paper, we stack neutral atomic hydrogen (HI) spectra for 9,720 star forming galaxies along the mass-metallicity relation. The sample is selected according to stellar mass (10$^9 \leq$ M$_{\star}$/M$_{\odot}\leq$10$^{11}$) and redshift ($0.02 \leq z \leq 0.05$) from the overlap of the Sloan Digital Sky Survey and Arecibo Legacy Fast ALFA survey. We confirm and quantify the strong anti-correlation between HI mass and gas-phase metallicity at fixed stellar mass. Furthermore, we show for the first time that the relationship between gas content and metallicity is consistent between different metallicity estimators, contrary to the weaker trends found with star formation which are known to depend on the observational techniques used to derive oxygen abundances and star formation rates. When interpreted in the context of theoretical work, this result supports a scenario where galaxies exist in an evolving equilibrium between gas, metallicity and star formation. The fact that deviations from this equilibrium are most strongly correlated with gas mass suggests that the scatter in the mass-metallicity relation is primarily driven by fluctuations in gas accretion.
We investigate the population of high-redshift ($3\leq z < 6$) AGN selected in the two deepest X-ray surveys, the 7 Ms \textit{Chandra} Deep Field-South and 2 Ms \textit{Chandra} Deep Field-North. Their outstanding sensitivity and spectral characterization of faint sources allow us to focus on the sub-$L_*$ regime (log$L_{\mathrm{X}}\lesssim44$), poorly sampled by previous works using shallower data, and the obscured population. Taking fully into account the individual photometric-redshift probability distribution functions, the final sample consists of $\approx102$ X-ray selected AGN at $3\leq z < 6$. The fraction of AGN obscured by column densities log$N_{\mathrm{H}}>23$ is $\sim0.6-0.8$, once incompleteness effects are taken into account, with no strong dependence on redshift or luminosity. We derived the high-redshift AGN number counts down to $F_{\mathrm{0.5-2\,keV}}=7\times10^{-18}\,\mathrm{erg\,cm^{-2}\,s^{-1}}$, extending previous results to fainter fluxes, especially at $z>4$. We put the tightest constraints to date on the low-luminosity end of AGN luminosity function at high redshift. The space-density, in particular, declines at $z>3$ at all luminosities, with only a marginally steeper slope for low-luminosity AGN. By comparing the evolution of the AGN and galaxy densities, we suggest that such a decline at high luminosities is mainly driven by the underlying galaxy population, while at low luminosities there are hints of an intrinsic evolution of the parameters driving nuclear activity. Also, the black-hole accretion rate density and star-formation rate density, which are usually found to evolve similarly at $z\lesssim3$, appear to diverge at higher redshifts.
Context: In the last six years, the VVV survey mapped 562 sq. deg. across the bulge and southern disk of the Galaxy. However, a detailed study of these regions, which includes $\sim 36$ globular clusters (GCs) and thousands of open clusters is by no means an easy challenge. High differential reddening and severe crowding along the line of sight makes highly hamper to reliably distinguish stars belonging to different populations and/or systems. Aims: The aim of this study is to separate stars that likely belong to the Galactic GC NGC 6544 from its surrounding field by means of proper motion (PM) techniques. Methods: This work was based upon a new astrometric reduction method optimized for images of the VVV survey. Results: Photometry over the six years baseline of the survey allowed us to obtain a mean precision of $\sim0.51$ mas/yr, in each PM coordinate, for stars with Ks < 15 mag. In the area studied here, cluster stars separate very well from field stars, down to the main sequence turnoff and below, allowing us to derive for the first time the absolute PM of NGC 6544. Isochrone fitting on the clean and differential reddening corrected cluster color magnitude diagram yields an age of $\sim$ 11-13 Gyr, and metallicity [Fe/H] = -1.5 dex, in agreement with previous studies restricted to the cluster core. We were able to derive the cluster orbit assuming an axisymmetric model of the Galaxy and conclude that NGC 6544 is likely a halo GC. We have not detected tidal tail signatures associated to the cluster, but a remarkable elongation in the galactic center direction has been found. The precision achieved in the PM determination also allows us to separate bulge stars from foreground disk stars, enabling the kinematical selection of bona fide bulge stars across the whole survey area. Our results show that VVV data is perfectly suitable for this kind of analysis.
Cold hydrogen gas is the raw fuel for star formation in galaxies, and its partition into atomic and molecular phases is a key quantity for galaxy evolution. In this paper, we combine Atacama Large Millimeter Array and Arecibo single-dish observations to estimate the molecular-to-atomic hydrogen mass ratio for massive star-forming galaxies at $z\sim$ 0.2 extracted from the HIGHz survey, i.e., some of the most massive gas-rich systems currently known. We show that the balance between atomic and molecular hydrogen in these galaxies is similar to that of local main-sequence disks, implying that atomic hydrogen has been dominating the cold gas mass budget of star-forming galaxies for at least the last three billion years. In addition, despite harboring gas reservoirs that are more typical of objects at the cosmic noon, HIGHz galaxies host regular rotating disks with low gas velocity dispersions suggesting that high total gas fractions do not necessarily drive high turbulence in the interstellar medium.
Hot dust-obscured galaxies (Hot DOGs) are a new population recently discovered in the \wise All-Sky survey. Multiwavelength follow-up observations suggest that they are luminous, dust-obscured quasars at high redshift. Here we present the JCMT SCUBA-2 850 $\mu m$ follow-up observations of 10 Hot DOGs. Four out of ten Hot DOGs have been detected at $>3\sigma$ level. Based on the IR SED decomposition approach, we derive the IR luminosities of AGN torus and cold dust components. Hot DOGs in our sample are extremely luminous with most of them having $L_{\rm IR}^{\rm tot}>10^{14} L_\odot$. The torus emissions dominate the total IR energy output. However, the cold dust contribution is still non-negligible, with the fraction of the cold dust contribution to the total IR luminosity $(\sim 8-24\%)$ being dependent on the choice of torus model. The derived cold dust temperatures in Hot DOGs are comparable to those in UV bright quasars with similar IR luminosity, but much higher than those in SMGs. Higher dust temperatures in Hot DOGs may be due to the more intense radiation field caused by intense starburst and obscured AGN activities. Fourteen and five submillimeter serendipitous sources in the 10 SCUBA-2 fields around Hot DOGs have been detected at $>3\sigma$ and $>3.5\sigma$ levels, respectively. By estimating their cumulative number counts, we confirm the previous argument that Hot DOGs lie in dense environments. Our results support the scenario in which Hot DOGs are luminous, dust-obscured quasars lying in dense environments, and being in the transition phase between extreme starburst and UV-bright quasars.
We present direct images in the H$\alpha$ and [SII]$\lambda \lambda$6717,6731 $\text{\AA}$ lines of the Galactic Supernova Remnant G109.1-1.0 (CTB 109). We confirm that the filaments detected are the optical counterpart of the X-ray and radio supernova remnant due to their high [SII]/H$\alpha$ line-ratios. We study for the first time the kinematics of the optical counterpart of SNR CTB 109 using the UNAM scanning Fabry-Perot interferometer PUMA. We estimate a systemic velocity of V$_{LSR}$=-50$\pm$6 km s$^{-1}$ for this remnant and an expansion velocity of V$_{exp}$=230$\pm$5 km s$^{-1}$. From this velocity value and taking into account previous studies about the kinematics of objects at that Galactic longitude we derive a distance to the SNR CTB 109 of 3.1$\pm$0.2 kpc, locating it in the Perseus arm. Using the [SII]$\lambda$6717/[SII]$\lambda$6731 line-ratio we find an electronic density value around n$_e$= 580 cm$^{-3}$. Considering that this remnant is evolving in a low density medium with higher density cloudlets responsible of the optical emission, we determine the age and energy deposited in the ISM by the supernova explosion (E$_0$) in both the Sedov-Taylor phase and the radiative phase. For both cases the age is of thousands of years and the E$_0$ is rather typical of SNRs containing simple pulsars so that, the energy released to the ISM cannot be used to distinguish between supernova remnants hosting typical pulsars from those hosting powerful magnetars as in the case of CTB 109.
The planar arrangement of nearly half the satellite galaxies of M31 has been a source of mystery and speculation since it was discovered. With a growing number of other host galaxies showing these satellite galaxy planes, their stability and longevity have become central to the debate on whether the presence of satellite planes are a natural consequence of prevailing cosmological models, or represent a challenge. Given the dependence of their stability on host halo shape, we look into how a galaxy plane's dark matter environment influences its longevity. An increased number of dark matter subhalos results in increased interactions that hasten the deterioration of an already-formed plane of satellite galaxies in spherical dark halos. The role of total dark matter mass fraction held in subhalos in dispersing a plane of galaxies present non trivial effects on plane longevity as well. But any misalignments of plane inclines to major axes of flattened dark matter halos lead to their lifetimes being reduced to < 3 Gyrs. Distributing > 40% of total dark mass in subhalos in the overall dark matter distribution results in a plane of satellite galaxies that is prone to change through the 5 Gyr integration time period.
Spectral energy distribution (SED) fitting of stellar population synthesis models is an important and popular way to constrain the physical parameters ---e.g., the ages, metallicities, masses for stellar population analysis. The previous works suggest that both blue-bands and red-bands photometry works for the SED-fitting. Either blue-domained or red-domained SED-fitting usually lead to the unreliable or biased results. Meanwhile, it seems that extending the wavelength coverage could be helpful. Since the Galaxy Evolution Explorer ({\it GALEX}) and Wide-field Infrared Survey Explorer (WISE) provide the FUV/NUV and mid-infrared $W1$/$W2$ band data, we extend the SED-fitting to a wider wavelength coverage. In our work, we analyzed the effect of adding the FUV/NUV and $W1$/$W2$ band to the optical and near-infrared $UBVRIJHK$ bands for the fitting with the Bruzual \& Charlot 2003 (BC03) models and {\sc galev} models. It is found that the FUV/NUV bands data affect the fitting results of both ages and metallicities much more significantly than that of the WISE $W1$/$W2$ band with the BC03 models. While for the {\sc galev} models, the effect of the WISE $W1$/$W2$ band for the metallicity fitting seems comparable to that of {\it GALEX} FUV/NUV bands, but for age the effect of the $W1$/$W2$ band seems less crucial than that of the FUV/NUV bands. Thus we conclude that the {\it GALEX} FUV/NUV bands are more crucial for the SED-fitting of ages and metallicities, than the other bands, and the high-quality UV data (with high photometry precision) are required.
The results of the study of the maser emission source IRAS 18316$-$0602 in the H$_2$O line at $\lambda =1.35$ cm are reported. The observations have been carried out at the RT-22 radio telescope of the Pushchino Radio Astronomy Observatory (Russia) since June 2002 until March 2017. Three superflares have been detected, in 2002, 2010, and 2016, with peak flux densities of 3400, 19,000, and 46,000 Jy, respectively. The results of the analysis of the superflares are given. The flares took place during periods of high maser activity in a narrow interval of radial velocities (40.5-42.5 km/s) and could be associated with the passage of a strong shock wave. During our monitoring the emission of three groups of features at radial velocities of about 41, 42, and 43 km/s dominated. The flare of 2016 was accompanied by a considerable increase in the flux densities of several features in an interval of 35--56 km/s.
In an effort to search for Ly$\alpha$ emission from circum- and intergalactic gas on scales of hundreds of kpc around $z\sim3$ quasars, and thus characterise the physical properties of the gas in emission, we have initiated an extensive fast-survey with the Multi Unit Spectroscopic Explorer (MUSE): Quasar Snapshot Observations with MUse: Search for Extended Ultraviolet eMission (QSO MUSEUM). In this work, we report the discovery of an enormous Ly$\alpha$ nebula (ELAN) around the quasar SDSS~J102009.99+104002.7 at $z=3.164$, which we followed-up with deeper MUSE observations. This ELAN spans $\sim297$ projected kpc, has an average Ly$\alpha$ surface brightness ${\rm SB}_{\rm Ly\alpha}\sim 6.04\times10^{-18}$ erg s$^{-1}$ cm$^{-2}$ arcsec$^{-2}$ (within the $2\sigma$ isophote), and is associated with an additional four, previously unknown embedded sources: two Ly$\alpha$ emitters and two faint active galactic nuclei (one Type-1 and one Type-2 quasar). By mapping at high significance the line-of-sight velocity in the entirety of the observed structure, we unveiled a large-scale coherent rotation-like pattern spanning $\sim300$ km s$^{-1}$ with a velocity dispersion of $<270$ km s$^{-1}$, which we interpret as a signature of the inspiraling accretion of substructures within the quasar's host halo. Future multiwavelength data will complement our MUSE observations, and are definitely needed to fully characterise such a complex system. None the less, our observations reveal the potential of new sensitive integral-field spectrographs to characterise the dynamical state of diffuse gas on large scales in the young Universe, and thereby witness the assembly of galaxies.
We report the relative abundances of the three stable isotopes of silicon, $^{28}$Si, $^{29}$Si and $^{30}$Si, across the Galaxy using the $v = 0, J = 1 \to 0$ transition of silicon monoxide. The chosen sources represent a range in Galactocentric radii ($R_{\rm GC}$) from 0 to 9.8 kpc. The high spectral resolution and sensitivity afforded by the GBT permit isotope ratios to be corrected for optical depths. The optical-depth-corrected data indicate that the secondary-to-primary silicon isotope ratios $^{29}{\rm Si}/^{28}{\rm Si}$ and $^{30}{\rm Si}/^{28}{\rm Si}$ vary much less than predicted on the basis of other stable isotope ratio gradients across the Galaxy. Indeed, there is no detectable variation in Si isotope ratios with $R_{\rm GC}$. This lack of an isotope ratio gradient stands in stark contrast to the monotonically decreasing trend with $R_{\rm GC}$ exhibited by published secondary-to-primary oxygen isotope ratios. These results, when considered in the context of the expectations for chemical evolution, suggest that the reported oxygen isotope ratio trends, and perhaps that for carbon as well, require further investigation. The methods developed in this study for SiO isotopologue ratio measurements are equally applicable to Galactic oxygen, carbon and nitrogen isotope ratio measurements, and should prove useful for future observations of these isotope systems.
Core collapse supernova (CCSN) rates suffer from large uncertainties as many CCSNe exploding in regions of bright background emission and significant dust extinction remain unobserved. Such a shortfall is particularly prominent in luminous infrared galaxies (LIRGs), which have high star formation (and thus CCSN) rates and host bright and crowded nuclear regions, where large extinctions and reduced search detection efficiency likely lead to a significant fraction of CCSNe remaining undiscovered. We present the first results of project SUNBIRD (Supernovae UNmasked By Infra-Red Detection), where we aim to uncover CCSNe that otherwise would remain hidden in the complex nuclear regions of LIRGs, and in this way improve the constraints on the fraction that is missed by optical seeing-limited surveys. We observe in the near-infrared 2.15 {\mu}m $K_s$-band, which is less affected by dust extinction compared to the optical, using the multi-conjugate adaptive optics imager GeMS/GSAOI on Gemini South, allowing us to achieve a spatial resolution that lets us probe close in to the nuclear regions. During our pilot program and subsequent first full year we have discovered three CCSNe and one candidate with projected nuclear offsets as small as 200 pc. When compared to the total sample of LIRG CCSNe discovered in the near-IR and optical, we show that our method is singularly effective in uncovering CCSNe in nuclear regions and we conclude that the majority of CCSNe exploding in LIRGs are not detected as a result of dust obscuration and poor spatial resolution.
Recent work has suggested that mid-IR wavelengths are optimal for estimating the mass-to-light ratios of stellar populations and hence the stellar masses of galaxies. We compare stellar masses deduced from spectral energy distribution (SED) models, fitted to multi-wavelength optical-NIR photometry, to luminosities derived from {\it WISE} photometry in the $W1$ and $W2$ bands at 3.6 and 4.5$\mu$m for non-star forming galaxies. The SED derived masses for a carefully selected sample of low redshift ($z \le 0.15$) passive galaxies agree with the prediction from stellar population synthesis models that $M_*/L_{W1} \simeq 0.6$ for all such galaxies, independent of other stellar population parameters. The small scatter between masses predicted from the optical SED and from the {\it WISE} measurements implies that random errors (as opposed to systematic ones such as the use of different IMFs) are smaller than previous, deliberately conservative, estimates for the SED fits. This test is subtly different from simultaneously fitting at a wide range of optical and mid-IR wavelengths, which may just generate a compromise fit: we are directly checking that the best fit model to the optical data generates an SED whose $M_*/L_{W1}$ is also consistent with separate mid-IR data. We confirm that for passive low redshift galaxies a fixed $M_*/L_{W1} = 0.65$ can generate masses at least as accurate as those obtained from more complex methods. Going beyond the mean value, in agreement with expectations from the models, we see a modest change in $M_*/L_{W1}$ with SED fitted stellar population age but an insignificant one with metallicity.
Pulsars around the Massive Black Hole (MBH) in the Galactic Center (GC) are expected to be revealed by the incoming facilities (e.g., the Square Kilometre Array). Under a full relativistic framework with the pulsar approximated as a test particle, we investigate the constraints on the spinning of the MBH by monitoring the timing of surrounding pulsars. For GC pulsars orbiting closely around the MBH (e.g., $\lesssim1000$AU), we find that full relativistic treatment in modeling accurately their timing signals can be necessary, as the relativistic signals are orders of magnitude larger than the time of arrival measurement accuracies. Although usually there are near-degeneracies among MBH spin parameters, the constraints on the spinning of the MBH are still very tight. By continuously monitoring a normal pulsar in orbits with a period of $\sim2.6$yr and an eccentricity of $0.3-0.9$ under timing precision of $1-5$ms, within $\sim 8$yr the spin magnitude and the orientations of the GC MBH can be constrained with $2\sigma$ error of $10^{-3}-10^{-2}$ and $10^{-1}-10^\circ$, respectively. Even for pulsars in orbits similar to the detected star S2/S0-2 or S0-102, we find that the spinning of the MBH can still be constrained within $4-8$yr, with the most significant constraints provided near the pericenter passage. If the proper motion of the pulsars with astrometric accuracy of $10\mu$as can also be collected along with the timing measurement, then the position, velocity, mass and the distance to the Solar System of the MBH can be constrained about $\sim10\mu$as, $\sim1\mu$as$/$yr, $\sim 10 M_\odot$ and $\sim1$pc, respectively.
In this work, we perform a statistical investigation towards 50 high-mass clumps using the data from the Bolocam Galactic Plane Survey (BGPS) and the Millimetre Astronomy Legacy Team 90-GHz survey (MALT90). Eleven dense molecular lines (N$_2$H$^+$(1-0), HNC(1-0), HCO$^+$(1-0), HCN(1-0), HN$^{13}$C(1-0), H$^{13}$CO$^+$(1-0), C$_2$H(1-0), HC$_3$N(10-9), SiO(2-1), $^{13}$CS(2-1) and HNCO$(4_{4,0}-3_{0,3}))$ are detected. N$_2$H$^+$ and HNC are shown to be good tracers for clumps in virous evolutionary stages since they are detected in all the fields. And the detection rates of N-bearing molecules decrease as the clumps evolve, but those of O-bearing species increase with evolution. Furthermore, the abundance ratios [N$_2$H$^+$]/[HCO$^+$] and Log([HC$_3$N]/[HCO$^+$]) decline with Log([HCO$^+$]) as two linear functions, respectively. This suggests the transformation of N$_2$H$^+$ and HC$_3$N to HCO$^+$ as the clumps evolve. We also find that C$_2$H is the most abundant molecule with an order of $10^{-8}$. Besides, three new infall candidates G010.214-00.324, G011.121-00.128, and G012.215-00.118(a) are discovered to have large-scaled infall motions and infall rates in the magnitude of $10^{-3}$ M$_\odot$ yr$^{-1}$.
MHONGOOSE is a deep survey of the neutral hydrogen distribution in a representative sample of 30 nearby disk and dwarf galaxies with HI masses from 10^6 to ~10^{11} M_sun, and luminosities from M_R ~ -12 to M_R ~ -22. The sample is selected to uniformly cover the available range in log(M_HI). Our extremely deep observations, down to HI column density limits of well below 10^{18} cm^{-2} - or a few hundred times fainter than the typical HI disks in galaxies - will directly detect the effects of cold accretion from the intergalactic medium and the links with the cosmic web. These observations will be the first ever to probe the very low-column density neutral gas in galaxies at these high resolutions. Combination with data at other wavelengths, most of it already available, will enable accurate modelling of the properties and evolution of the mass components in these galaxies and link these with the effects of environment, dark matter distribution, and other fundamental properties such as halo mass and angular momentum. MHONGOOSE can already start addressing some of the SKA-1 science goals and will provide a comprehensive inventory of the processes driving the transformation and evolution of galaxies in the nearby universe at high resolution and over 5 orders of magnitude in column density. It will be a Nearby Galaxies Legacy Survey that will be unsurpassed until the advent of the SKA, and can serve as a highly visible, lasting statement of MeerKAT's capabilities.
In the framework of the STREGA (STRucture and Evolution of the GAlaxy) survey, two fields around the globular cluster Pal 12 were observed with the aim of detecting the possible presence of streams and/or an extended halo. The adopted stellar tracers are the Main Sequence, Turn-off and Red Giant Branch stars. We discuss the lumi- nosity function and the star counts in the observed region covering about 2 tidal radii, confirming that Pal 12 appears to be embedded in the Sagittarius Stream. Adopting an original approach to separate cluster and field stars, we do not find any evidence of sig- nificant extra-tidal Pal 12 stellar populations. The presence of the Sagittarius stream seems to have mimicked a larger tidal radius in previous studies. Indeed, adopting a King model, a redetermination of this value gives r_T = 0.22 +- 0.1 deg.
We analyze 3 mm emission of the ultraluminous infrared galaxy Arp 220 for spatially-resolved structure and spectral properties of the merger nuclei. ALMA archival data at ~0.05" resolution are used for extensive visibility fitting and deep imaging of continuum emission. The data are fitted well with two concentric components for each nucleus, such as two Gaussians or one Gaussian plus one exponential disk. The larger components in individual nuclei are similar in shape and extent, ~100-150 pc, to the cm-wave emission due to supernovae. They are therefore identified with the known starburst nuclear disks. The smaller components in both nuclei have about a few 10 pc sizes and peak brightness temperatures (Tb) more than twice higher than in previous single-Gaussian fitting. They correspond to the dust emission that we find centrally concentrated in both nuclei by subtracting the plasma emission measured at 33 GHz. The dust emission in the western nucleus is found to have a peak Tb ~ 530 K and a full width at half maximum of about 20 pc. This component is estimated to have a bolometric luminosity on the order of 10^{12.5} Lsun and a 20 pc-scale luminosity surface density 10^{15.5} Lsun/kpc^2. A luminous AGN is a plausible energy source for these high values while other explanations remain to be explored. Our continuum image also reveals a third structural component of the western nucleus --- a pair of faint spurs perpendicular to the disk major axis. We attribute it to a bipolar outflow from the highly inclined (i ~ 60 deg) western nuclear disk.
Both simulations and observations have confirmed that the spin of haloes/galaxies is correlated with the large scale structure (LSS) with a mass dependence such that the spin of low-mass haloes/galaxies tend to be parallel with the LSS, while that of massive haloes/galaxies tend to be perpendicular with the LSS. It is still unclear how this mass dependence is built up over time. We use N-body simulations to trace the evolution of the halo spin-LSS correlation and find that at early times the spin of all halo progenitors is parallel with the LSS. As time goes on, mass collapsing around massive halo is more isotropic, especially the recent mass accretion along the slowest collapsing direction is significant and it brings the halo spin to be perpendicular with the LSS. Adopting the $fractional$ $anisotropy$ (FA) parameter to describe the degree of anisotropy of the large-scale environment, we find that the spin-LSS correlation is a strong function of the environment such that a higher FA (more anisotropic environment) leads to an aligned signal, and a lower anisotropy leads to a misaligned signal. In general, our results show that the spin-LSS correlation is a combined consequence of mass flow and halo growth within the cosmic web. Our predicted environmental dependence between spin and large-scale structure can be further tested using galaxy surveys.
Polarized Galactic foregrounds are one of the primary sources of systematic error in measurements of the B-mode polarization of the Cosmic Microwave Background (CMB). Experiments are becoming increasingly sensitive to complexities in the foreground frequency spectra that are not captured by standard parametric models, potentially affecting our ability to efficiently separate out these components. Employing a suite of dust models encompassing a variety of physical effects, we simulate observations of a future seven-band CMB experiment to assess the impact of these complexities on parametric component separation. We identify configurations of frequency bands that minimize the `model errors' caused by fitting simple parametric models to more complex `true' foreground spectra, which bias the inferred CMB signal. We find that: (a) fits employing a simple two parameter modified blackbody (MBB) dust model tend to produce significant bias in the recovered polarized CMB signal in the presence of physically realistic dust foregrounds; (b) generalized MBB models with three additional parameters reduce this bias in most cases, but non-negligible biases can remain, and can be hard to detect; and (c) line of sight effects, which give rise to frequency decorrelation, and the presence of iron grains are the most problematic complexities in the dust emission for recovering the true CMB signal. More sophisticated simulations will be needed to demonstrate that future CMB experiments can successfully mitigate these more physically realistic dust foregrounds.
We present a statistical characterization of the $\gamma$-ray emission from the four \emph{Fermi}-LAT sources: FR I radio galaxy NGC 1275, BL Lac Mrk 421, FSRQs B2 1520+31 and PKS 1510-089 detected almost continuously over a time integration of 3-days between August 2008 - October 2015. The observed flux variation is large, spanning $\gtrsim 2$ orders of magnitude between the extremes except for Mrk~421. We compute the flux distributions and compare with Gaussian and lognormal ones. We find that the 3 blazars have distribution consistent with a lognormal, suggesting that the variability is of a non-linear, multiplicative nature. This is further supported by the computation of the flux-rms relation, which is observed to be linear for the 3 blazars. However, for NGC 1275, the distribution does not seem to be represented either by a lognormal or a Gaussian, while its flux-rms relation is still found to be linear. We also compute the power spectra, which suggest the presence of a break, but are consistent with typical scale-free power-law shot noise. The results are broadly consistent with the statistical properties of the magnetic reconnection powered minijets-in-a-jet model. We discuss other possible scenarios and implications of these observations on jet processes and connections with the central engine.
We present deep high-resolution (R~15,000) and high-quality UVES optical spectrophotometry of nine planetary nebulae with dual-dust chemistry. We compute physical conditions from several diagnostics. Ionic abundances for a large number of ions of N, O, Ne, S, Cl, Ar, K, Fe and Kr are derived from collisionally excited lines. Elemental abundances are computed using state-of-the-art ionization correction factors. We derive accurate C/O ratios from optical recombination lines. We have re-analyzed additional high-quality spectra of 14 PNe from the literature following the same methodology. Comparison with asymptotic giant branch models reveals that about half of the total sample objects are consistent with being descendants of low-mass progenitor stars (M < 1.5 Msun). Given the observed N/O, C/O, and He/H ratios, we cannot discard that some of the objects come from more massive progenitor stars (M > 3--4 Msun) that have suffered a mild HBB. None of the objects seem to be a descendant of very massive progenitors. We propose that in most of the planetary nebulae studied here, the PAHs have been formed through the dissociation of the CO molecule. The hypothesis of a last thermal pulse that turns O-rich PNe into C-rich PNe is discarded, except in three objects, that show C/O > 1. We also discuss the possibility of a He pre-enrichment to explain the most He-enriched objects. We cannot discard other scenarios like extra mixing, stellar rotation or binary interactions to explain the chemical abundances behaviour observed in our sample.
Many galaxy clusters have giant halos of non-thermal radio emission, indicating the presence of relativistic electrons in the clusters. Relativistic protons may also be accelerated by merger and/or accretion shocks in galaxy clusters. These cosmic-ray (CR) protons are expected to produce gamma-rays through inelastic collisions with thermal gas of the intra-cluster medium. Despite of intense efforts in searching for high-energy gamma-ray emission from galaxy clusters, conclusive evidence is still missing so far. Here we report the discovery of $\ge100$ MeV gamma-rays from the Coma cluster at a confidence level of $\sim6 \sigma$ with the {\it Fermi} Large Area Telescope (LAT). The result is obtained with the unbinned likelihood analysis of the nine years of {\it Fermi}-LAT Pass 8 data. The gamma-ray emission shows an extended spatial morphology roughly coincident with the giant radio halo, with an apparent excess at the southwest of the cluster. The integral energy flux in the energy range of 0.1-300 GeV is about $3\times 10^{-12}{\rm \ erg\ cm^{-2}\ s^{-1}}$. The spectrum of the gamma-ray emission is soft with a photon index of $\simeq-2.7$. Interpreting the gamma-ray emission as arising from CR proton interaction, we find that the volume-averaged value of the CR to thermal pressure ratio in the Coma cluster is about $\sim 2\%$. Our results show that galaxy clusters are a new type of GeV gamma-ray sources, and they are probably also giant reservoirs of CR protons.
Spatially resolved observations of protoplanetary discs are revealing that their inner regions can be warped or broken from the outer disc. A few mechanisms are known to lead to such 3D structures; among them, the interaction with a stellar companion. We perform a 3D SPH simulation of a circumbinary disc misaligned by $60^\circ$ with respect to the binary orbital plane. The inner disc breaks from the outer regions, precessing as a rigid body, and leading to a complex evolution. As the inner disc precesses, the misalignment angle between the inner and outer discs varies by more than $100^\circ$. Different snapshots of the evolution are post-processed with a radiative transfer code, in order to produce observational diagnostics of the process. Even though the simulation was produced for the specific case of a circumbinary disc, most of the observational predictions hold for any disc hosting a precessing inner rim. Synthetic scattered light observations show strong azimuthal asymmetries, where the pattern depends strongly on the misalignment angle between inner and outer disc. The asymmetric illumination of the outer disc leads to azimuthal variations of the temperature structure, in particular in the upper layers, where the cooling time is short. These variations are reflected in asymmetric surface brightness maps of optically thick lines, as CO $J$=3-2. The kinematical information obtained from the gas lines is unique in determining the disc structure. The combination of scattered light images and (sub-)mm lines can distinguish between radial inflow and misaligned inner disc scenarios.
According to the current paradigm of circumstellar disk evolution, gas-rich primordial disks evolve into gas-poor debris disks compose of second-generation dust. To explore the transition between these phases, we searched for $^{12}$CO, $^{13}$CO, and C$^{18}$O emission in seven dust-rich debris disks around young A-type stars, using ALMA in Band 6. We discovered molecular gas in three debris disks. In all these disks, the $^{12}$CO line was optically thick, highlighting the importance of less abundant molecules in reliable mass estimates. Supplementing our target list by literature data, we compiled a volume-limited sample of dust-rich debris disks around young A-type stars within 150 pc. We obtained a CO detection rate of 11/16 above a $^{12}$CO J=2$-$1 line luminosity threshold of $\sim 1.4 \times 10 ^4$ Jykms$^{-1}$pc$^2$ in the sample. This high incidence implies that the presence of CO gas in bright debris disks around young A-type stars is likely more the rule than the exception. Interestingly, dust-rich debris disks around young FG-type stars exhibit, with the same detectability threshold as for A-type stars, significantly lower gas incidence. While the transition from protoplanetary to debris phase is associated with a drop of dust content, our results exhibit a large spread in the CO mass in our debris sample, with peak values comparable to those in protoplanetary Herbig Ae disks. In the particularly CO-rich debris systems the gas may have primordial origin, characteristic of a hybrid disk.
Annihilating dark matter particles in nearby subhalos could generate potentially observable fluxes of gamma rays, unaccompanied by emission at other wavelengths. Furthermore, this gamma-ray emission is expected to be spatially extended, providing us with a powerful way to discriminate dark matter subhalos from other astrophysical gamma-ray sources. Fermi has detected two dark matter subhalo candidates which exhibit a statistically significant degree of spatial extension (3FGL J2212.5+0703 and 3FGL J1924.8-1034). It has been argued that the most likely non-dark matter interpretation of these observations is that they are each in fact multiple nearby point sources, too close to one another on the sky to be individually resolved. In this study, we consider the ability of next generation gamma-ray telescopes to spatially resolve the gamma-ray emission from subhalo candidates, focusing on the proposed e-ASTROGAM mission. We find that such an instrument could significantly clarify the nature of Fermi's dark matter subhalo candidates, and provide an unprecedented level of sensitivity to the presence of annihilating dark matter in nearby subhalos.
We have developed a model atom for Cu with which we perform statistical equilibrium computations that allow us to compute the line formation of Cu I lines in stellar atmospheres without assuming Local Thermodynamic Equilibrium (LTE). We validate this model atom by reproducing the observed line profiles of the Sun, Procyon and eleven metal-poor stars. Our sample of stars includes both dwarfs and giants. Over a wide range of stellar parameters we obtain excellent agreement among different Cu I lines. The eleven metal-poor stars have iron abundances in the range -4.2 <= [Fe/H] <= -1.4, the weighted mean of the [Cu/Fe] ratios is -0.22 dex, with a scatter of -0.15 dex. This is very different from the results from LTE analysis (the difference between NLTE and LTE abundances reaches 1 dex) and in spite of the small size of our sample it prompts for a revision of the Galactic evolution of Cu.
Links to: arXiv, form interface, find, astro-ph, recent, 1709, contact, help (Access key information)
We study the kinematic properties of the Eastern Banded Structure (EBS) and Hydra I overdensity using exquisite proper motions derived from the Sloan Digital Sky Survey (SDSS) and Gaia source catalog. Main sequence turn-off stars in the vicinity of the EBS are identified from SDSS photometry; we use the proper motions and, where applicable, spectroscopic measurements of these stars to probe the kinematics of this apparent stream. We find that the EBS and Hydra I share common kinematic and chemical properties with the nearby Monoceros Ring. In particular, the proper motions of the EBS, like Monoceros, are indicative of prograde rotation (V_phi ~ 180-220 km/s), which is similar to the Galactic thick disc. The kinematic structure of stars in the vicinity of the EBS suggest that it is not a distinct stellar stream, but rather marks the "edge" of the Monoceros Ring. The EBS and Hydra I are the latest substructures to be linked with Monoceros, leaving the Galactic anti-centre a mess of interlinked overdensities which likely share a unified, Galactic disc origin.
We present cosmological radiation-hydrodynamic simulations, performed with the code Ramses-RT, of radiatively-driven outflows in a massive quasar host halo at z = 6. Our simulations include both single- and multi-scattered radiation pressure on dust from a quasar and are compared against simulations performed with thermal feedback. For radiation pressure-driving, we show that there is a critical quasar luminosity above which a galactic outflow is launched, set by the equilibrium of gravitational and radiation forces. While this critical luminosity is unrealistically high in the single-scattering limit for plausible black hole masses, it is in line with a 3 x 10^9 M_SUN black hole accreting at its Eddington limit, if infrared (IR) multi-scattering radiation pressure is included. The outflows are fast (v > 1000 km/s) and strongly mass-loaded with peak mass outflow rates 10^3 - 10^4 M_SUN/yr, but short-lived (< 10 Myr). Outflowing material is multi-phase, though predominantly composed of cool gas, forming via a thermal instability in the shocked swept-up component. Radiation pressure- and thermally-driven outflows both affect their host galaxies profoundly, but in different, complementary ways. Thermally-driven outflows couple more efficiently to diffuse halo gas, generating more powerful, hotter and more volume-filling outflows. IR radiation, through its ability to penetrate dense gas via diffusion, is more efficient at ejecting gas from the bulge. The combination of gas ejection through outflows with internal pressurisation by trapped IR radiation leads to a complete shut down of star formation in the bulge. We hence argue that radiation pressure-driven feedback may be a crucial ingredient in regulating star formation in compact starbursts, especially during the quasar's "obscured"' phase.
In symmetric gravitating systems experiencing rapid mass loss, particle orbits change almost instantaneously, which can lead to the development of a sharply contoured density profile, including singular caustics for collisionless systems. This framework can be used to model a variety of dynamical systems, such as accretion disks following a massive black hole merger and dwarf galaxies following violent early star formation feedback. Particle interactions in the high-density peaks seem a promising source of observable signatures of these mass loss events (i.e. a possible EM counterpart for black hole mergers or strong gamma-ray emission from dark matter annihilation around young galaxies), because the interaction rate depends on the square of the density. We study post-mass-loss density profiles, both analytic and numerical, in idealised cases and present arguments and methods to extend to any general system. An analytic derivation is presented for particles on Keplerian orbits responding to a drop in the central mass. We argue that this case, with initially circular orbits, gives the most sharply contoured profile possible. We find that despite the presence of a set of singular caustics, the total particle interaction rate is reduced compared to the unperturbed system; this is a result of the overall expansion of the system dominating over the steep caustics. Finally we argue that this result holds more generally, and the loss of central mass decreases the particle interaction rate in any physical system.
The BLASTPol observations of Vela C have provided the most detailed characterization of the polarization fraction $p$ and dispersion in polarization angles $S$ for a molecular cloud. We compare the observed distributions of $p$ and $S$ with those obtained in synthetic observations of simulations of molecular clouds, assuming homogeneous grain alignment. We find that the orientation of the mean magnetic field relative to the observer has a significant effect on the $p$ and $S$ distributions. These distributions for Vela C are most consistent with synthetic observations where the mean magnetic field is close to the line-of-sight. Our results point to apparent magnetic disorder in the Vela C molecular cloud, although it can be due to either an inclination effect (i.e., observing close to the mean field direction) or significant field tangling from strong turbulence/low magnetization. The joint correlations of $p$ with column density and of $S$ with column density for the synthetic observations generally agree poorly with the Vela C joint correlations, suggesting that understanding these correlations require a more sophisticated treatment of grain alignment physics.
We present a photometric analysis of 65 galaxies in the rich cluster Abell 1689 at $z=0.183$, using the Hubble Space Telescope Advanced Camera for Surveys archive images in the rest-frame $V$-band. We perform two-dimensional multi-component photometric decomposition of each galaxy adopting different models of the surface-brightness distribution. We present an accurate morphological classification for each of the sample galaxies. For 50 early-type galaxies, we fit both a de Vaucouleurs and S\'ersic law; S0s are modelled by also including a disc component described by an exponential law. Bars of SB0s are described by the profile of a Ferrers ellipsoid. For the 15 spirals, we model a S\'ersic bulge, exponential disc, and, when required, a Ferrers bar component. We derive the Fundamental Plane by fitting 40 early-type galaxies in the sample, using different surface-brightness distributions. We find that the tightest plane is that derived by S\'ersic bulges. We find that bulges of spirals lie on the same relation. The Fundamental Plane is better defined by the bulges alone rather than the entire galaxies. Comparison with local samples shows both an offset and rotation in the Fundamental Plane of Abell 1689.
It is now possible for hydrodynamical simulations to reproduce a representative galaxy population. Accordingly, it is timely to assess critically some of the assumptions of traditional semi-analytic galaxy formation models. We use the Eagle simulations to assess assumptions built into the Galform semi-analytic model, focussing on those relating to baryon cycling, angular momentum and feedback. We show that the assumption in Galform that newly formed stars have the same specific angular momentum as the total disc leads to a significant overestimate of the total stellar specific angular momentum of disc galaxies. In Eagle, stars form preferentially out of low specific angular momentum gas in the interstellar medium (ISM) due to the assumed gas density threshold for stars to form, leading to more realistic galaxy sizes. We find that stellar mass assembly is similar between Galform and Eagle but that the evolution of gas properties is different, with various indications that the rate of baryon cycling in Eagle is slower than is assumed in Galform. Finally, by matching individual galaxies between Eagle and Galform, we find that an artificial dependence of AGN feedback and gas infall rates on halo mass doubling events in Galform drives most of the scatter in stellar mass between individual objects. Put together our results suggest that the Galform semi-analytic model can be significantly improved in light of recent advances.
The processes taking place in the outermost reaches of spiral disks (the 'proto-disk') are intimately connected to the build-up of mass and angular momentum in galaxies. The thinness of spiral disks suggests that the activity is mostly quiescent and presumably this region is fed by cool flows coming into the halo from the intergalactic medium. While there is abundant evidence for the presence of a circumgalactic medium (CGM) around disk galaxies as traced by quasar absorption lines, it has been very difficult to connect this material to the outer gas disk. This has been a very difficult transition region to explore because baryon tracers are hard to observe. In particular, HI disks have been argued to truncate at a critical column density N(H) $\approx 3\times 10^{19}$ cm$^{-2}$ at 30 kpc for an L* galaxy where the gas is vulnerable to the external ionizing background. But new deep observations of nearby L* spirals (e.g. Milky Way, NGC 2997) suggest that HI disks may extend much further than recognised to date, up to 60 kpc at N(H) $\approx 10^{18}$ cm$^{-2}$. Motivated by these observations, here we show that a clumpy outer disk of dense clouds or cloudlets is potentially detectable to much larger radii and lower HI column densities than previously discussed. This extended proto-disk component is likely to explain some of the MgII forest seen in quasar spectra as judged from absorption-line column densities and kinematics. We fully anticipate that the armada of new radio facilities and planned HI surveys coming online will detect this extreme outer disk (scree) material. We also propose a variant on the successful 'Dragonfly' technique to go after the very weak H$\alpha$ signals expected in the proto-disk region.
We describe the POLAMI program for the monitoring of all four Stokes parameters of a sample of bright radio-loud active galactic nuclei with the IRAM 30m telescope at 3.5 and 1.3mm. The program started in October 2006 and accumulated, until August 2014, 2300 observations at 3.5mm, achieving a median time sampling interval of 22 days for the sample of 37 sources. This first paper explains the source selection, mostly blazars, the observing strategy and data calibration, and gives the details of the instrumental polarisation corrections. The sensitivity (1sigma) reached at 3.5mm is 0.5% (linear polarisation degree), 4.7 deg. (polarisation angle), and 0.23% (circular polarisation), while the corresponding values at 1.3mm are 1.7%, 9.9 deg., and 0.72%, respectively. The data quality is demonstrated by the time sequences of our calibrators Mars and Uranus. For the quasar 3C286, widely used as a linear polarisation calibrator, we give improved estimates of its linear polarisation, and show for the first time occasional detections of its weak circular polarisation, which suggests a small level of variability of the source at millimeter wavelengths.
We analyse the circular polarisation data accumulated in the first 7 years of the POLAMI project introduced in an accompanying paper (Agudo et al.). In the 3mm wavelength band, we acquired more than 2600 observations, and all but one of our 37 sample sources were detected, most of them several times. For most sources, the observed distribution of the degree of circular polarisation is broader than that of unpolarised calibrators, indicating that weak (<0.5%) circular polarisation is present most of the time. Our detection rate and the maximum degree of polarisation found, 2.0%, are comparable to previous surveys, all made at much longer wavelengths. We argue that the process generating circular polarisation must not be strongly wavelength dependent, and we propose that the widespread presence of circular polarisation in our short wavelength sample dominated by blazars is mostly due to Faraday conversion of the linearly polarised synchrotron radiation in the helical magnetic field of the jet. Circular polarisation is variable, most notably on time scales comparable to or shorter than our median sampling interval <1 month. Longer time scales of about one year are occasionally detected, but severely limited by the weakness of the signal. At variance with some longer wavelength investigations we find that the sign of circular polarisation changes in most sources, while only 7 sources, including 3 already known, have a strong preference for one sign. The degrees of circular and linear polarisation do not show any systematic correlation. We do find however one particular event where the two polarisation degrees vary in synchronism during a time span of 0.9 years. The paper also describes a novel method for calibrating the sign of circular polarisation observations.
We report on the first results of the POLAMI program, a simultaneous 3.5 and 1.3mm full-Stokes-polarisation monitoring of a sample of 36 of the brightest active galactic nuclei in the northern sky with the IRAM 30m Telescope. Through a systematic statistical study of data taken from October 2006 (from December 2009 for the case of the 1.3mm observations) to August 2014, we characterise the variability of the total flux density and linear polarisation. We find that all sources in the sample are highly variable in total flux density at both 3.5 and 1.3mm, as well as in spectral index, that is found to be optically thin in general. The total flux-density variability at 1.3mm is found, in general, to be faster, and to have larger amplitude and flatter PSD slopes than 3.5mm. The polarisation degree is on average larger at 1.3mm than at 3.5mm, by a factor of 2.6. The variability of linear polarisation degree is faster and has higher fractional amplitude than for total flux density, with the typical time scales during prominent polarisation peaks being significantly faster at 1.3mm than at 3.5mm. The polarisation angle at both 3.5 and 1.3mm is highly variable. Most of the sources show one or two excursions of >180 deg. on time scales from a few weeks to about a year during the course of our observations. The 3.5 and 1.3mm polarisation angle evolution follow rather well each other, although the 1.3mm data show a clear preference to more prominent variability on the short time scales, i.e. weeks. The data are compatible with multi-zone models of conical jets involving smaller emission regions for the shortest-wavelength emitting sites. Such smaller emitting regions should also be more efficient in energising particle populations. The data also favours the integrated emission at 1.3mm to have better ordered magnetic fields than the one at 3.5mm.
Motivated by the existence of the relationship between the dynamical state of clusters and the shape of the velocity dispersion profiles (VDP), we study the VDPs for Gaussian (G) and Non-Gaussian (NG) systems for a subsample of clusters from the Yang catalog. The groups cover a redshift interval of $0.03\leq z\leq0.1$ with halo mass $\geq 10^{14}$M$_{\odot}$. We use a robust statistical method, Hellinger Distance, to classify the dynamical state of the systems according to their velocity distribution. The stacked VDP of each class, G and NG, is then determined using either Bright or Faint galaxies. The stacked VDP for G groups displays a central peak followed by a monotonically decreasing trend which indicates a predominance of radial orbits, with the Bright stacked VDP showing lower velocity dispersions in all radii. The distinct features we find in NG systems are manifested not only by the characteristic shape of VDP, with a depression in the central region, but also by a possible higher infall rate associated with galaxies in the Faint stacked VDP.
Based on current models of the Cosmic X-ray Background (CXB), heavily obscured Active Galactic Nuclei (AGN) are expected to make up ~10% of the peak emission of the CXB and ~20% of the total population of AGN, yet few of these sources have been recorded and characterized in current surveys. Here we present the Chandra follow-up observation of 14 AGN detected by Swift-BAT. For five sources in the sample, NuSTAR observations in the 3-80 keV band are also available. The X-ray spectral fitting over the 0.3-150 keV energy range allows us to determine the main X-ray spectral parameters, such as the photon index and the intrinsic absorption, of these objects, and to make hypotheses on the physical structures responsible for the observed spectra. We find that 13 of the 14 objects are absorbed AGN, and one is a candidate Compton thick AGN, having intrinsic absorption NH>1E24 cm^{-2}. Finally, we verified that the use of NuSTAR observations is strategic to strongly constrain the properties of obscured AGN, since the best-fit values we obtained for parameters such as the power-law photon index and the intrinsic absorption NH changed sometimes significantly fitting the spectra with and without the use of NuSTAR data.
We conducted an SiO maser survey towards 221 O-rich AGB stars with the aim of identifying maser emission associated with the Sagittarius stellar stream. In this survey, maser emission was detected in 44 targets, of which 35 were new detections. All of these masers are within 5 kpc of the Sun. We also compiled a Galactic SiO maser catalogue including ~2300 SiO masers from the literature. The distribution of these SiO masers give a scale height of 0.40 kpc, while 42 sources deviate from the Galactic plane by more than 1.2 kpc, half of which were found in this survey. Regarding SiO masers in the disc, we found both the rotational speeds and the velocity dispersions vary with the Galactic plane distance. Assuming Galactic rotational speed $\Theta$0 = 240 km/s , we derived the velocity lags are 15 km/s and 55 km/s for disc and off-plane SiO masers respectively. Moreover, we identified three groups with significant peculiar motions (with 70% confidence). The most significant group is in the thick disc that might trace stream/peculiar motion of the Perseus arm. The other two groups are mainly made up of off-plane sources. The northern and southern off-plane sources were found to be moving at ~33 km/s and 54 km/s away from the Galactic plane, respectively. Causes of these peculiar motions are still unclear. For the two off-plane groups, we suspect they are thick disc stars whose kinematics affected by the Sgr stellar stream or very old Sgr stream debris.
Aims. We aimed at investigating the structure and kinematics of the gaseous disk and outflows around the massive YSO S255 NIRS3 in the S255IR-SMA1 dense clump. Methods. Observations of the S255IR region were carried out with ALMA at two epochs in the compact and extended configurations. Results. We serendipitously detected a new, never predicted, bright maser line at about 349.1 GHz, which most probably represents the CH$_3$OH $14_{1} - 14_{0}$ A$^{- +}$ transition. The emission covers most of the 6.7 GHz methanol maser emission area of almost 1$^{\prime\prime}$ in size and shows a velocity gradient in the same sense as the disk rotation. No variability was found on the time interval of several months. It is classified as Class II maser and probably originates in a ring at a distance of several hundreds AU from the central star.
We investigate the triggering of star formation and the formation of stellar clusters in molecular clouds that form as the ISM passes through spiral shocks. The spiral shock compresses gas into $\sim$100 pc long main star formation ridge, where clusters forming every 5-10 pc along the merger ridge. We use a gravitational potential based cluster finding algorithm, which extracts individual clusters, calculates their physical properties and traces cluster evolution over multiple time steps. Final cluster masses at the end of simulation range between 1000 and 30000 M$_{\odot}$ with their characteristic half-mass radii between 0.1 pc and 2 pc. These clusters form by gathering material from 10-20 pc size scales. Clusters also show a mass - specific angular momentum relation, where more massive clusters have larger specific angular momentum due to the larger size scales, and hence angular momentum from which they gather their mass. The evolution shows that more massive clusters experiences hierarchical merging process, which increases stellar age spreads up to 2-3 Myr. Less massive clusters appear to grow by gathering nearby recently formed sinks, while more massive clusters with their large global gravitational potentials are increasing their mass growth from gas accretion.
We present stellar-dynamical measurements of the central supermassive black hole (SMBH) in the S0 galaxy NGC 307, using adaptive-optics IFU data from VLT-SINFONI. We investigate the effects of including dark-matter haloes as well as multiple stellar components with different mass-to-light (M/L) ratios in the dynamical modeling. Models with no halo and a single stellar component yield a relatively poor fit with a low value for the SMBH mass ($7.0 \pm 1.0 \times 10^{7} M_{\odot}$) and a high stellar M/L ratio (K-band M/L = $1.3 \pm 0.1$). Adding a halo produces a much better fit, with a significantly larger SMBH mass ($2.0 \pm 0.5 \times 10^{8} M_{\odot}$) and a lower M/L ratio ($1.1 \pm 0.1$). A model with no halo but with separate bulge and disc components produces a similarly good fit, with a slightly larger SMBH mass ($3.0 \pm 0.5 \times 10^{8} M_{\odot}$) and an identical M/L ratio for the bulge component, though the disc M/L ratio is biased high (disc M/L $ = 1.9 \pm 0.1$). Adding a halo to the two-stellar-component model results in a much more plausible disc M/L ratio of $1.0 \pm 0.1$, but has only a modest effect on the SMBH mass ($2.2 \pm 0.6 \times 10^{8} M_{\odot}$) and leaves the bulge M/L ratio unchanged. This suggests that measuring SMBH masses in disc galaxies using just a single stellar component and no halo has the same drawbacks as it does for elliptical galaxies, but also that reasonably accurate SMBH masses and bulge M/L ratios can be recovered (without the added computational expense of modeling haloes) by using separate bulge and disc components.
The Galactic Center is the closest galactic nucleus that can be studied with unprecedented angular resolution and sensitivity. We summarize recent basic observational results on Sagittarius A* and the conditions for star formation in the central stellar cluster. We cover results from the radio, infrared, and X-ray domain and include results from simulation as well. From (sub-)mm and near-infrared variability and near-infrared polarization data we find that the SgrA* system (supermassive black hole spin, a potential temporary accretion disk and/or outflow) is well ordered in its geometrical orientation and in its emission process that we assume to reflect the accretion process onto the supermassive black hole (SMBH).
We report the discovery, spectroscopic confirmation, and mass modelling of the gravitationally lensed quasar system PS J0630-1201. The lens was discovered by matching a photometric quasar catalogue compiled from Pan-STARRS and WISE photometry to the Gaia DR1 catalogue, exploiting the high spatial resolution of the latter (FWHM $\sim $0.1") to identify the three brightest components of the lens. Follow-up spectroscopic observations with the WHT confirm the multiple objects are quasars at redshift $z_{q}=3.34$. Further follow-up with Keck AO high-resolution imaging reveals that the system is composed of two lensing galaxies and the quasar is lensed into a $\sim$2.8" separation four-image cusp configuration with a fifth image clearly visible, and a 1.0" arc due to the lensed quasar host galaxy. The system is well-modelled with two singular isothermal ellipsoids, reproducing the position of the fifth image. We discuss future prospects for measuring time delays between the images and constraining any offset between mass and light using the faintly detected Einstein arcs associated with the quasar host galaxy.
Gaia's exceptional resolution (FWHM $\sim$ 0.1$^{\prime\prime}$) allows identification and cataloguing of the multiple images of gravitationally lensed quasars. We investigate a sample of 49 known lensed quasars in the SDSS footprint, with image separations less than 2$^{\prime\prime}$, and find that 8 are detected with multiple components in the first Gaia data release. In the case of the 41 single Gaia detections, we generally are able to distinguish these lensed quasars from single quasars when comparing Gaia flux and position measurements to those of Pan-STARRS and SDSS. This is because the multiple images of these lensed quasars are typically blended in ground-based imaging and therefore the total flux and a flux-weighted centroid are measured, which can differ significantly from the fluxes and centroids of the individual components detected by Gaia. We compare the fluxes through an empirical fit of Pan-STARRS griz photometry to the wide optical Gaia bandpass values using a sample of isolated quasars. The positional offsets are calculated from a recalibrated astrometric SDSS catalogue. Applying flux and centroid difference criteria to spectroscopically confirmed quasars, we discover 4 new sub-arcsecond-separation lensed quasar candidates which have two distinct components of similar colour in archival CFHT or HSC data. Our method based on single Gaia detections can be used to identify the $\sim$ 1400 lensed quasars with image separation above 0.5$^{\prime\prime}$, expected to have only one image bright enough to be detected by Gaia.
Sensitive, imaging observations of the 1.1 mm dust continuum emission from a 1 deg^2 area collected with the AzTEC bolometer camera on the Large Millimeter Telescope are presented. A catalog of 1545 compact sources is constructed based on a Wiener-optimization filter. These sources are linked to larger clump structures identified in the Bolocam Galactic Plane Survey. Hydrogen column densities are calculated for all sources and mass and mean volume densities are derived for the subset of sources for which kinematic distances can be assigned. The AzTEC sources are localized, high density peaks within the massive clumps of molecular clouds and comprise 5-15% of the clump mass. We examine the role of the gravitational instability in generating these fragments by comparing the mass of embedded AzTEC sources to the Jeans' mass of the parent BGPS object. For sources with distances less than 6 kpc the fragment masses are comparable to the clump Jeans' mass, despite having isothermal Mach numbers between 1.6 and 7.2. AzTEC sources linked to ultra-compact HII regions have mass surface densities greater than the critical value implied by the mass-size relationship of infrared dark clouds with high mass star formation while AzTEC sources associated with Class II methanol masers have mass surface densities greater than 0.7 g cm^{-2} that approaches the proposed threshold required to form massive stars.
A number of groups have recently been active in searching for gradients in the ob- served Faraday rotation measure (RM) across jets of Active Galactic Nuclei (AGNs) on various scales and estimating their reliability. Such RM structures provide direct evidence for the presence of an azimuthal magnetic field component, which may be associated with a helical jet magnetic field, as is expected based on the results of many theoretical studies. We present new parsec-scale RM maps of 4 AGNs here, and analyze their transverse RM structures together with those for 5 previously published RM maps. All these jets display transverse RM gradients with significances of at least 3{\sigma}. This is part of an ongoing effort to establish how common transverse RM gradients that may be associated with helical or toroidal magnetic fields are in AGNs on parsec scales.
The Galactic Centre is of fundamental astrophysical interest, but existing near-infrared surveys fall short to cover it adequately. Here we introduce the GALACTICNUCLEUS survey, a JHKs imaging survey of the centre of the Milky Way with a 0.2" angular resolution. We present the observations of Field 1 of our survey, centred approximately on SgrA* with an approximate size of 7.95 ' x 3.43 '. We describe the observational set-up and data reduction pipeline and discuss the quality of the data. Finally, we present some preliminary analysis of the data. The data were acquired with the near-infrared camera HAWK-I at the ESO VLT. Short readout times in combination with the speckle holography algorithm allowed us to produce final images with a stable, Gaussian PSF of 0.2" FWHM. Astrometric calibration is achieved via the VVV survey and photometric calibration is based on the SIRIUS/IRSF survey. The quality of the data is assessed by comparison between observations of the same field with different detectors of HAWK-I and at different times. We reach 5 sigma detection limits of approximately J = 22, H = 21, and Ks = 20. The photometric uncertainties are less than 0.05 at J < 20, H < 17 and Ks < 16. We distinguish five stellar populations in the colour-magnitude diagrams; three of them appear to belong to foreground spiral arms, and the other two correspond to a high- and a low-extinction star groups at the Galactic Centre. We use our data to analyse the near-infrared extinction curve and conclude that it can be described very well by a power-law with an index of JHKs = 2.31 +- 0.03. We do not find any evidence that this index depends on the position along the line-of-sight, or on the absolute value of the extinction. We produce extinction maps that show the clumpiness of the ISM at the Galactic Centre. Finally, we estimate that the majority of the stars have solar or super-solar metallicity.
We show that mid-infrared data from the all-sky WISE survey can be used as a robust photometric redshift indicator for powerful radio AGN, in the absence of other spectroscopic or multi-band photometric information. Our work is motivated by a desire to extend the well-known K-z relation for radio galaxies to the wavelength range covered by the all-sky WISE mid-infrared survey. Using the LARGESS radio spectroscopic sample as a training set, and the mid-infrared colour information to classify radio sources, we generate a set of redshift probability distributions for the hosts of high-excitation and low-excitation radio AGN. We test the method using spectroscopic data from several other radio AGN studies, and find good agreement between our WISE-based redshift estimates and published spectroscopic redshifts out to z ~ 1 for galaxies and z ~ 3-4 for radio-loud QSOs. Our chosen method is also compared against other classification methods and found to perform reliably. This technique is likely to be particularly useful in the analysis of upcoming large-area radio surveys with SKA pathfinder telescopes, and our code is publicly available. As a consistency check, we show that our WISE-based redshift estimates for sources in the 843 MHz SUMSS survey reproduce the redshift distribution seen in the CENSORS study up to z ~ 2. We also discuss two specific applications of our technique for current and upcoming radio surveys; an interpretation of large scale HI absorption surveys, and a determination of whether low-frequency peaked spectrum sources lie at high redshift.
We performed a multiwavelength study towards HII region Sh2-104. New maps of 12CO J=1-0 and 13CO J=1-0 were obtained from the Purple Mountain Observatory (PMO) 13.7 m radio telescope. Sh2-104 displays a double-ring structure. The outer ring with a radius of 4.4 pc is dominated by 12 um, 500 um, 12CO J=1-0, and 13CO J=1-0 emission, while the inner ring with a radius of 2.9 pc is dominated by 22 um and 21 cm emission. We did not detect CO emission inside the outer ring. The north-east portion of the outer ring is blueshifted, while the south-west portion is redshifted. The present observations have provided evidence that the collected outer ring around Sh2-104 is a two-dimensional structure. From the column density map constructed by the Hi-GAL survey data, we extract 21 clumps. About 90\% of all the clumps will form low-mass stars. A power-law fit to the clumps yields M=281Msun(r/pc)^1.31. The selected YSOs are associated with the collected material on the edge of Sh2-104. The derived dynamical age of Sh2-104 is 1.6*10^6 yr. Compared the Sh2-104 dynamical age with the YSOs timescale and the fragmentation time of the molecular ring, we further confirm that collect-and-collapse process operates in this region, indicating a positive feedback from a massive star for surrounding gas.
We present results from the analysis of 2997 fundamental mode RR Lyrae variables located in the Small Magellanic Cloud (SMC). For these objects near-infrared time-series photometry from the VISTA survey of the Magellanic Clouds system (VMC) and visual light curves from the OGLE IV survey are available. In this study the multi-epoch $K_{\rm s}$-band VMC photometry was used for the first time to derive intensity-averaged magnitudes of the SMC RR Lyrae stars. We determined individual distances to the RR Lyrae stars from the near-infrared period-absolute magnitude-metallicity ($PM_{K_{\rm s}}Z$) relation, which has a number of advantages in comparison with the visual absolute magnitude-metallicity ($M_{V}-{\rm [Fe/H]}$) relation, such as a smaller dependence of the luminosity on interstellar extinction, evolutionary effects and metallicity. The distances we have obtained were used to study the three-dimensional structure of the SMC. The distribution of the SMC RR Lyrae stars is found to be ellipsoidal. The actual line-of-sight depth of the SMC is in the range from 1 to 10 kpc, with an average depth of 4.3 $\pm$ 1.0 kpc. We found that RR Lyrae stars in the eastern part of the SMC are affected by interactions of the Magellanic Clouds. However, we do not see a clear bimodality in the distribution of RR Lyrae stars as observed for red clump (RC) stars.
Following the first two annual intensity mapping workshops at Stanford in March 2016 and Johns Hopkins in June 2017, we report on the recent advances in theory, instrumentation and observation that were presented in these meetings and some of the opportunities and challenges that were identified looking forward. With preliminary detections of CO, [CII], Lya and low-redshift 21cm, and a host of experiments set to go online in the next few years, the field is rapidly progressing on all fronts, with great anticipation for a flood of new exciting results. This current snapshot provides an efficient reference for experts in related fields and a useful resource for nonspecialists. We begin by introducing the concept of line-intensity mapping and then discuss the broad array of science goals that will be enabled, ranging from the history of star formation, reionization and galaxy evolution to measuring baryon acoustic oscillations at high redshift and constraining theories of dark matter, modified gravity and dark energy. After reviewing the first detections reported to date, we survey the experimental landscape, presenting the parameters and capabilities of relevant instruments such as COMAP, mmIMe, AIM-CO, CCAT-p, TIME, CONCERTO, CHIME, HIRAX, HERA, STARFIRE, MeerKAT/SKA and SPHEREx. Finally, we describe recent theoretical advances: different approaches to modeling line luminosity functions, several techniques to separate the desired signal from foregrounds, statistical methods to analyze the data, and frameworks to generate realistic intensity map simulations.
Links to: arXiv, form interface, find, astro-ph, recent, 1709, contact, help (Access key information)
We present dynamical measurements of the central mass-to-light ratio ($M/L$) of a sample of 27 low-mass early-type ATLAS$^{3D}$ galaxies. We consider all ATLAS$^{3D}$ galaxies with 9.7$<$log(M$_\star/$M$_\odot)$$<$10.5 in our analysis, selecting out galaxies with available high-resolution Hubble Space Telescope (HST) data, and eliminating galaxies with significant central color gradients or obvious dust features. We use the HST images to derive mass models for these galaxies and combine these with the central velocity dispersion values from ATLAS$^{3D}$ data to obtain a central dynamical $M/L$ estimate. These central dynamical $M/L$s are higher than dynamical $M/L$s derived at larger radii and stellar population estimates of the galaxy centers in $\sim$80\% of galaxies, with a median enhancement of $\sim$14\% and a statistical significance of 3.3$\sigma$. We show that the enhancement in the central $M/L$ is best described either by the presence of black holes in these galaxies or by radial IMF variations. Assuming a black hole model, we derive black hole masses for the sample of galaxies. In two galaxies, NGC 4458 and NGC 4660, the data suggests significantly over-massive BHs, while in most others only upper limits are obtained. We also show that the level of $M/L$ enhancements we see in these early-type galaxy nuclei are consistent with the larger enhancements seen in ultracompact dwarf galaxies (UCDs), supporting the scenario where massive UCDs are created by stripping galaxies of these masses.
Do the theorised different formation mechanisms of fast and slow rotators produce an observable difference in their star formation histories? To study this we identify quenching slow rotators in the MaNGA sample by selecting those which lie below the star forming sequence and identify a sample of quenching fast rotators which were matched in stellar mass. This results in a total sample of 194 kinematically classified galaxies, which is agnostic to visual morphology. We use u-r and NUV-u colours from SDSS and GALEX and an existing inference package, STARPY, to conduct a first look at the onset time and exponentially declining rate of quenching of these galaxies. An Anderson-Darling test on the distribution of the inferred quenching rates across the two kinematic populations reveals they are statistically distinguishable ($3.2\sigma$). We find that fast rotators quench at a much wider range of rates than slow rotators, consistent with a wide variety of physical processes such as secular evolution, minor mergers, gas accretion and environmentally driven mechanisms. Quenching is more likely to occur at rapid rates ($\tau \lesssim 1~\rm{Gyr}$) for slow rotators, in agreement with theories suggesting slow rotators are formed in dynamically fast processes, such as major mergers. Interestingly, we also find that a subset of the fast rotators quench at these same rapid rates as the bulk of the slow rotator sample. We therefore discuss how the total gas mass of a merger, rather than the merger mass ratio, may decide a galaxy's ultimate kinematic fate.
While 2% of active galactic nuclei (AGNs) exhibit narrow emission lines with line-of-sight velocities that are significantly offset from the velocity of the host galaxy's stars, the nature of these velocity offsets is unknown. We investigate this question with Chandra/ACIS and Hubble Space Telescope/Wide Field Camera 3 observations of seven velocity-offset AGNs at z<0.12, and all seven galaxies have a central AGN but a peak in emission that is spatially offset by < kpc from the host galaxy's stellar centroid. These spatial offsets are responsible for the observed velocity offsets and are due to shocks, either from AGN outflows (in four galaxies) or gas inflowing along a bar (in three galaxies). We compare our results to a velocity-offset AGN whose velocity offset originates from a spatially offset AGN in a galaxy merger. The optical line flux ratios of the offset AGN are consistent with pure photoionization, while the optical line flux ratios of our sample are consistent with contributions from photoionization and shocks. We conclude that these optical line flux ratios could be efficient for separating velocity-offset AGNs into subgroups of offset AGNs -- which are important for studies of AGN fueling in galaxy mergers -- and central AGNs with shocks -- where the outflows are biased towards the most energetic outflows that are the strongest drivers of feedback.
The gas surrounding galaxies outside their disks or interstellar medium and inside their virial radii is known as the circumgalactic medium (CGM). In recent years this component of galaxies has assumed an important role in our understanding of galaxy evolution owing to rapid advances in observational access to this diffuse, nearly invisible material. Observations and simulations of this component of galaxies suggest that it is a multiphase medium characterized by rich dynamics and complex ionization states. The CGM is a source for a galaxy's star-forming fuel, the venue for galactic feedback and recycling, and perhaps the key regulator of the galactic gas supply. We review our evolving knowledge of the CGM with emphasis on its mass, dynamical state, and coevolution with galaxies. Observations from all redshifts and from across the electromagnetic spectrum indicate that CGM gas has a key role in galaxy evolution. We summarize the state of this field and pose unanswered questions for future research.
We present a study of photometric redshift performance for galaxies and active galactic nuclei detected in deep radio continuum surveys. Using two multi-wavelength datasets, over the NOAO Deep Wide Field Survey Bo\"otes and COSMOS fields, we assess photometric redshift (photo-z) performance for a sample of $\sim 4,500$ radio continuum sources with spectroscopic redshifts relative to those of $\sim 63,000$ non radio-detected sources in the same fields. We investigate the performance of three photometric redshift template sets as a function of redshift, radio luminosity and infrared/X-ray properties. We find that no single template library is able to provide the best performance across all subsets of the radio detected population, with variation in the optimum template set both between subsets and between fields. Through a hierarchical Bayesian combination of the photo-z estimates from all three template sets, we are able to produce a consensus photo-z estimate which equals or improves upon the performance of any individual template set.
We analyze the diffuse ionized gas (DIG) in the first Galactic quadrant from l=18deg to 40deg using radio recombination line (RRL) data from the Green Bank Telescope. These data allow us to distinguish DIG emission from HII region emission and thus study the diffuse gas essentially unaffected by confusion from discrete sources. We find that the DIG has two dominant velocity components, one centered around 100km/s associated with the luminous HII region W43, and the other centered around 45km/s not associated with any large HII region. Our analysis suggests that the two velocity components near W43 may be caused by non-circular streaming motions originating near the end of the Galactic bar. At lower Galactic longitudes, the two velocities may instead arise from gas at two distinct distances from the Sun, with the most likely distances being ~6kpc for the 100km/s component and ~12kpc for the 45km/s component. We show that the intensity of diffuse Spitzer GLIMPSE 8.0um emission caused by excitation of polyaromatic hydrocarbons (PAHs) is correlated with both the locations of discrete HII regions and the intensity of the RRL emission from the DIG. This implies that the soft ultra-violet photons responsible for creating the infrared emission have a similar origin as the harder ultra-violet photons required for the RRL emission. The 8.0um emission increases with RRL intensity but flattens out for directions with the most intense RRL emission, suggesting that PAHs are partially destroyed by the energetic radiation field at these locations.
We present Gemini/GNIRS cross-dispersed near-infrared spectra of 12 nearby early-type galaxies, with the aim of testing commonly used stellar population synthesis models. We select a subset of galaxies from the atlas/ sample which span a wide range of ages (SSP-equivalent ages of 1--15~Gyr) at approximately solar metallicity. We derive star formation histories using four different stellar population synthesis models, namely those of [bruzual_stellar_2003], Conroy, Gunn \& White (2009; 2010), [maraston_stellar_2011] and [vazdekis_uv-extended_2016]. We compare star formation histories derived from near-infrared spectra with those derived from optical spectra using the same models. We find that while all models agree in the optical, the derived star formation histories vary dramatically from model to model in the near-infrared. We find that this variation is largely driven by the choice of stellar spectral library, such that models including high quality spectral libraries provide the best fits to the data, and are the most self-consistent when comparing optically-derived properties with near-infrared ones. We also find the impact of age variation in the near-infrared to be subtle, and largely encoded in the shape of the continuum, meaning that the common approach of removing continuum information with a high-order polynomial greatly reduces our ability to constrain ages in the near-infrared.
We present our observations of the optical intra-day variability (IDV) in $\gamma$-ray BL Lac object Mrk 501. The observations were run with the 1.02 m and 2.4 m optical telescopes at Yunnan Observatories from 2005 April to 2012 May. The light curve at the $R$ band on 2010 May 15 passes both variability tests (the $F$ test and the ANOVA test). A flare within the light curve on 2010 May 15 has a magnitude change $\Delta m = 0.03 \pm 0.005_{\rm{stat}} \pm 0.007_{\rm{sys}}$ mag, \textbf{a darkening timescale of $\tau_{\rm{d}}=$ 26.7 minutes}, and an amplitude of IDV $Amp=2.9\% \pm0.7\%$. A decline \textbf{described by 11 consecutive flux measurements} within the flare can be fitted linearly with a Pearson's correlation coefficient $r = 0.945$ at the confidence level of $> 99.99\%$. Under the assumptions that the IDV is tightly connected to the mass of the black hole, \textbf{and that the flare duration, being two times $\tau_{\rm{d}}$, is representative of the minimum characteristic timescale, we can derive upper bounds to the mass of the black hole}. In the case of the Kerr black hole, the timescale of $\Delta t_{\rm{min}}^{\rm{ob}}=$ 0.89 hours gives $M_{\bullet}\la 10^{9.20} M_{\odot}$, which is consistent with measurements reported in the literature. This agreement indicates that the hypothesis about $M_{\bullet}$ and $\Delta t_{\rm{min}}^{\rm{ob}}$ is consistent with the measurements/data.
Phylogenetics is a widely used concept in evolutionary biology. It is the
reconstruction of evolutionary history by building trees that represent
branching patterns and sequences. These trees represent shared history, and it
is our intention for this approach to be employed in the analysis of Galactic
history. In Galactic archaeology the shared environment is the interstellar
medium in which stars form and provides the basis for tree-building as a
methodological tool.
Using elemental abundances of solar-type stars as a proxy for DNA, we built
in Jofre et al 2017 such an evolutionary tree to study the chemical evolution
of the solar neighbourhood. In this proceeding we summarise these results and
discuss future prospects.
We describe methods designed to determine the astrophysical parameters of quasars based on spectra coming from the red and blue spectrophotometers of the Gaia satellite. These methods principally rely on two already published algorithms that are the weighted principal component analysis and the weighted phase correlation. The presented approach benefits from a fast implementation; an intuitive interpretation as well as strong diagnostic tools on the potential errors that may arise during predictions. The production of a semi-empirical library of spectra as they will be observed by Gaia is also covered and subsequently used for validation purpose. We detail the pre-processing that is necessary in order for these spectra to be fully exploitable by our algorithms along with the procedures that are used in order to predict the redshifts of the quasars; their continuum slopes; the total equivalent width of their emission lines and whether these are broad absorption line (BAL) quasars or not. Performances of these procedures were assessed in comparison with the Extremely Randomized Trees learning method and were proven to provide better results on the redshift predictions and on the ratio of correctly classified observations though the probability of detection of BAL quasars remains restricted by the low resolution of these spectra as well as by their limited signal-to-noise ratio. Finally, the triggering of some warning flags allows us to obtain an extremely pure subset of redshift predictions where approximately 99% of the observations come along with absolute errors that are below 0.1.
The Central Molecular Zone (CMZ) of the Galactic Center has to date only been fully mapped at mm wavelengths with singledish telescopes, with resolution about 30$^{\prime\prime}$ (1.2 pc). Using CARMA, we mapped the innermost 0.25 square degrees of the CMZ over the region between -0$.\kern-.25em ^{^\circ}$2$ \leq l \leq $0$.\kern-.25em ^{^\circ}$5 and -0$.\kern-.25em ^{^\circ}$2$ \leq b \leq $0$.\kern-.25em ^{^\circ}$2 (90$\times$50 pc) with spatial and spectral resolution of $\sim$10$^{\prime\prime}$ (0.4 pc) and $\sim$2.5 km/s, respectively. We provide a catalog of 3mm continuum sources as well as spectral line images of SiO(J=2-1), HCO$^{+}$(J=1-0), HCN(J=1-0), N$_{2}$H$^{+}$(J=1-0), and CS(J=2-1), with velocity coverage VLSR= -200 to 200 km/s. To recover the large scale structure resolved out by the interferometer, the continuum-subtracted spectral line images were combined with data from the Mopra 22-m telescope survey, thus providing maps containing all spatial frequencies down to the resolution limit. We find that integrated intensity ratio of I(HCN)/I(HCO$^{+}$) is anti-correlated with the intensity of the 6.4 keV Fe K$\alpha$, which is excited either by high energy photons or low energy cosmic rays, and the gas velocity dispersion as traced by HCO$^{+}$ is correlated with Fe K$\alpha$ intensity. The intensity ratio and velocity dispersion patterns are consistent with variation expected from the interaction of low energy cosmic rays with molecular gas.
We study the prototypical Seyfert 2 galaxy, Markarian 3, based on imaging and high-resolution spectroscopy observations taken by the Chandra X-ray Observatory. We construct a deconvolved X-ray image, which reveals the S-shaped morphology of the hot gas in the narrow line region (NLR). While this morphology is similar to the radio and [O III] emission, the distribution of the X-ray gas is broader than that obtained at these other wavelengths. By mapping the density and temperature distribution of the hot gas in the NLR, we demonstrate the presence of shocks towards the west ($M=2.5^{+1.0}_{-0.6}$) and east ($M=1.5^{+1.0}_{-0.5}$). Moreover, we compute the flux ratios between the [O III] and $0.5-2$ keV band X-ray luminosity and show that it is non-uniform in the NLR with the western side of the NLR being more highly ionized. In addition, based on the Chandra grating data we investigate the line ratios of the Si XIII triplet, which are not consistent with pure photoionization. Based on these results, we suggest that in the NLR of Mrk 3 both photoionization and collisional ionization act as excitation mechanisms. We conclude that the canonical picture, in which photoionization is solely responsible for exciting the interstellar medium in the NLR of Seyfert galaxies, may be overly simplistic. Given that weak and small-scale radio jets are commonly detected in Seyfert galaxies, it is possible that shock heating plays a non-negligible role in the NLR of these galaxies.
Blazar jets are renowned for their rapid violent variability and multiwavelength flares, however, the physical processes responsible for these flares are not well understood. In this paper we develop a time-dependent inhomogeneous fluid jet emission model for blazars. We model optically thick radio flares for the first time and show that they are delayed with respect to the prompt optically thin emission by ~ months to decades, with a lag that increases with the jet power and observed wavelength. This lag is caused by a combination of the travel time of the flaring plasma to the optically thin radio emitting sections of the jet and the slow rise time of the radio flare. We predict two types of flares: symmetric flares - with the same rise and decay time, which occur for flares whose duration is shorter than both the radiative lifetime and the geometric path-length delay timescale; extended flares - whose luminosity tracks the power of particle acceleration in the flare, which occur for flares with a duration longer than both the radiative lifetime and geometric delay. Our model naturally produces orphan X-ray and $\gamma$-ray flares. These are caused by flares which are only observable above the quiescent jet emission in a narrow band of frequencies. Our model is able to successfully fit to the observed multiwavelength flaring spectra and lightcurves of PKS1502+106 across all wavelengths, using a transient flaring front located within the broad-line region.
The short-spacing problem describes the inherent inability of
radio-interferometric arrays to measure the integrated flux and structure of
diffuse emission associated with extended sources. New interferometric arrays,
such as SKA, require solutions to efficiently combine interferometer and
single-dish data.
We present a new and open source approach for merging single-dish and cleaned
interferometric data sets requiring a minimum of data manipulation while
offering a rigid flux determination and full high angular resolution. Our
approach combines single-dish and cleaned interferometric data in the image
domain. This approach is tested for both Galactic and extragalactic HI data
sets. Furthermore, a quantitative comparison of our results to commonly used
methods is provided. Additionally, for the interferometric data sets of NGC4214
and NGC5055, we study the impact of different imaging parameters as well as
their influence on the combination for NGC4214.
The approach does not require the raw data (visibilities) or any additional
special information such as antenna patterns. This is advantageous especially
in the light of upcoming radio surveys with heterogeneous antenna designs.
The Gaia benchmark stars (GBS) are very bright stars of different late spectral types, luminosities and metallicities. They are well-known in the Galactic archaeology community because they are widely used to calibrate and validate the automatic pipelines delivering parameters of on-going and future spectroscopic surveys. The sample provides us with consistent fundamental parameters as well as a library of high resolution and high signal-to-noise spectra. This allows the community to study details of high resolution spectroscopy and to compare results between different survey pipelines, putting the GBS at the heart of this community. Here we discuss some results arising from using the GBS as main data source for spectral analyses.
Aims. We use accurate data on distances and radial velocities of galaxies
around the Local Group, as well as around 14 other massive nearby groups, to
estimate their radius of the zero-velocity surface, $R_0$, which separates any
group against the global cosmic expansion.
Methods. Our $R_0$ estimate was based on fitting the data to the velocity
field expected from the spherical infall model, including effects of the
cosmological constant. The reported uncertainties were derived by a Monte Carlo
simulation.
Results. Testing various assumptions about a location of the group
barycentre, we found the optimal estimates of the radius to be
$0.91\pm0.05$~Mpc for the Local Group, and $0.93\pm0.02$~Mpc for a synthetic
group stacked from 14 other groups in the Local Volume. Under the standard
Planck model parameters, these quantities correspond to the total mass of the
group $\sim (1.6\pm0.2) 10^{12} M_{\odot}$. Thus, we are faced with the
paradoxical result that the total mass estimate on the scale of $R_0 \approx
(3- 4) R_{vir}$ is only $~60$% of the virial mass estimate. Anyway, we conclude
that wide outskirts of the nearby groups do not contain a large amount of
hidden mass outside their virial radius.
The S-type asymptotic giant branch (AGB) star $\pi^{1}$ Gruis has a known companion at a separation of $\approx$400 AU. The envelope structure, including an equatorial torus and a fast bipolar outflow, is rarely seen in the AGB phase and is particularly unexpected in such a wide binary system. Therefore a second, closer companion has been suggested, but the evidence is not conclusive. The new ALMA $^{12}$CO and $^{13}$CO $J$=3-2 data, together with previously published $^{12}$CO $J$=2-1 data from the Submillimeter Array (SMA), and the $^{12}$CO $J$=5-4 and $J$=9-8 lines observed with Herschel/Heterodyne Instrument for the Far-Infrared (HIFI), is modeled with the 3D non-LTE radiative transfer code SHAPEMOL. The data analysis clearly confirms the torus-bipolar structure. The 3D model of the CSE that satisfactorily reproduces the data consists of three kinematic components: a radially expanding torus with velocity slowly increasing from 8 to 13 km s$^{-1}$ along the equator plane; a radially expanding component at the center with a constant velocity of 14 km s$^{-1}$; and a fast, bipolar outflow with velocity proportionally increasing from 14 km s$^{-1}$ at the base up to 100 km s$^{-1}$ at the tip, following a linear radial dependence. The results are used to estimate an average mass-loss rate during the creation of the torus of 7.7$\times$10$^{-7}$ M$_{\odot}$ yr$^{-1}$. The total mass and linear momentum of the fast outflow are estimated at 7.3$\times$10$^{-4}$ M$_{\odot}$ and 9.6$\times$10$^{37}$ g cm s$^{-1}$, respectively. The momentum of the outflow is in excess (by a factor of about 20) of what could be generated by radiation pressure alone, in agreement with recent findings for more evolved sources. The best-fit model also suggests a $^{12}$CO/$^{13}$CO abundance ratio of 50. Possible shaping scenarios for the gas envelope are discussed
High-resolution mid-infrared observations carried out by the Spitzer Space Telescope allowed one to resolve the fine structure of many astrospheres. In particular, they showed that the astrosphere around the B0.7 Ia star kappa Cas (HD 2905) has a clear-cut arc structure with numerous cirrus-like filaments beyond it. Previously, we suggested a physical mechanism for the formation of such filamentary structures. Namely, we showed theoretically that they might represent the non-monotonic spatial distribution of the interstellar dust in astrospheres (viewed as filaments) caused by interaction of the dust grains with the interstellar magnetic field disturbed in the astrosphere due to colliding of the stellar and interstellar winds. In this paper, we invoke this mechanism to explain the structure of the astrosphere around kappa Cas. We performed 3D magnetohydrodynamic modelling of the astrosphere for realistic parameters of the stellar wind and space velocity. The dust dynamics and the density distribution in the astrosphere were calculated in the framework of a kinetic model. It is found that the model results with the classical MRN size distribution of dust in the interstellar medium do not match the observations, and that the observed filamentary structure of the astrosphere can be reproduced only if the dust is composed mainly of big (micron-sized) grains. Comparison of the model results with observations allowed us to estimate parameters (number density and magnetic field strength) of the surrounding interstellar medium.
The large majority of the accreting supermassive black holes in the Universe are obscured by large columns of gas and dust. The location and evolution of this obscuring material have been the subject of intense research in the past decades, and are still highly debated. A decrease in the covering factor of the circumnuclear material with increasing accretion rates has been found by studies carried out across the electromagnetic spectrum. The origin of this trend has been suggested to be driven either by the increase in the inner radius of the obscuring material with incident luminosity due to the sublimation of dust; by the gravitational potential of the black hole; by radiative feedback; or by the interplay between outflows and inflows. However, the lack of a large, unbiased and complete sample of accreting black holes, with reliable information on gas column density, luminosity and mass, has left the main physical mechanism regulating obscuration unclear. Using a systematic multi-wavelength survey of hard X-ray-selected black holes, here we show that radiation pressure on dusty gas is indeed the main physical mechanism regulating the distribution of the circumnuclear material. Our results imply that the bulk of the obscuring dust and gas in these objects is located within the sphere of influence of the black hole (i.e., a few to tens of parsecs), and that it can be swept away even at low radiative output rates. The main physical driver of the differences between obscured and unobscured accreting black holes is therefore their mass-normalized accretion rate.
Links to: arXiv, form interface, find, astro-ph, recent, 1709, contact, help (Access key information)
We present ALMA Cycle 4 observations of CO(1-0), CO(3-2), and $^{13}$CO(3-2) line emission in the brightest cluster galaxy of RXJ0821+0752. This is one of the first detections of $^{13}$CO line emission in a galaxy cluster. Half of the CO(3-2) line emission originates from two clumps of molecular gas that are spatially offset from the galactic center. These clumps are surrounded by diffuse emission that extends $8~{\rm kpc}$ in length. The detected $^{13}$CO emission is confined entirely to the two bright clumps, with any emission outside of this region lying below our detection threshold. Two distinct velocity components with similar integrated fluxes are detected in the $^{12}$CO spectra. The narrower component ($60~{\rm km}~{\rm s}^{-1}$ FWHM) is consistent in both velocity centroid and linewidth with $^{13}$CO(3-2) emission, while the broader ($130-160~{\rm km}~{\rm s}^{-1}$), slightly blueshifted wing has no associated $^{13}$CO(3-2) emission. A simple local thermodynamic model indicates that the $^{13}$CO emission traces $2.1\times 10^{9}~{\rm M}_\odot$ of molecular gas. Isolating the $^{12}$CO velocity component that accompanies the $^{13}$CO emission yields a CO-to-H$_2$ conversion factor of $\alpha_{\rm CO}=2.3~{\rm M}_{\odot}~({\rm K~km~s^{-1}})^{-1}$, which is a factor of two lower than the Galactic value. Adopting the Galactic CO-to-H$_2$ conversion factor in brightest cluster galaxies may therefore overestimate their molecular gas masses by a factor of two. This is within the object-to-object scatter from extragalactic sources, so calibrations in a larger sample of clusters are necessary in order to confirm a sub-Galactic conversion factor.
We present a survey of 56 massive star-forming regions in the 44 GHz methanol maser transition made with the Karl G. Jansky Very Large Array (VLA); 24 of the 56 fields showed maser emission. The data allow us to demonstrate associations, at arcsecond precision, of the Class I maser emission with outflows, HII regions, and shocks traced by 4.5 micron emission. We find a total of 83 maser components with linewidths ranging from 0.17 to 3.3 km s$^{-1}$ with a nearly flat distribution and a median value of 1.1 km s$^{-1}$. The relative velocities of the masers with respect to the systemic velocity of the host clouds range from $-$2.5 to 3.1 km s$^{-1}$ with a distribution peaking near zero. We also study the correlation between the masers and the so-called extended green objects (EGOs) from the GLIMPSE survey. Multiple sources in each field are revealed from IR images as well as from centimeter continuum emission from VLA archival data; in the majority of cases the 44 GHz masers are positionally correlated with EGOs which seem to trace the younger sources in the fields. We report a possible instance of a 44 GHz maser associated with a low-mass protostar. If confirmed, this region will be the fifth known star-forming region that hosts Class I masers associated with low-mass protostars. We discuss three plausible cases of maser variability.
We address a question whether the observed light curves of X-ray flares originating deep in galactic cores can give us independent constraints on the mass of the central supermassive black hole. To this end we study four brightest flares that have been recorded from Sagittarius A*. They all exhibit an asymmetric shape consistent with a combination of two intrinsically separate peaks that occur at a certain time-delay with respect to each other, and are characterized by their mutual flux ratio and the profile of raising/declining parts. Such asymmetric shapes arise naturally in the scenario of a temporary flash from a source orbiting near a super- massive black hole, at radius of only 10-20 gravitational radii. An interplay of relativistic effects is responsible for the modulation of the observed light curves: Doppler boosting, gravitational redshift, light focusing, and light-travel time delays. We find the flare properties to be in agreement with the simulations (our ray-tracing code sim5lib). The inferred mass for each of the flares comes out in agreement with previous estimates based on orbits of stars; the latter have been observed at radii and over time-scales two orders of magnitude larger than those typical for the X-ray flares, so the two methods are genuinely different. We test the reliability of the method by applying it to another object, namely, the Seyfert I galaxy RE J1034+396.
Density profiles of isolated cores derived from thermal dust continuum emission rely on models of dust properties, such as mass opacity, which are poorly constrained. With complementary measures from near-infrared extinction maps, we can assess the reliability of commonly-used dust models. In this work, we compare Herschel-derived maps of the optical depth with equivalent maps derived from CFHT WIRCAM near-infrared observations for three isolated cores: CB68, L429, and L1552. We assess the dust opacities provided from four models: OH1a, OH5a, Orm1, and Orm4. Although the consistency of the models differs between the three sources, the results suggest that the optical properties of dust in the envelopes of the cores are best described by either silicate and bare graphite grains (e.g., Orm1) or carbonaceous grains with some coagulation and either thin or no ice mantles (e.g., OH5a). None of the models, however, individually produced the most consistent optical depth maps for every source. The results suggest that either the dust in the cores is not well described by any one dust property model, the application of the dust models cannot be extended beyond the very center of the cores, or more complex SED fitting functions are necessary.
We compare the optical properties of the host galaxies of radio-quiet (RQ) and radio-loud (RL) Type 2 active galactic nuclei (AGNs) to infer whether the jet production efficiency depends on the host properties or is determined just by intrinsic properties of the accretion flows. We carefully select galaxies from SDSS, FIRST, and NVSS catalogs. We confirm previous findings that the fraction of RL AGNs depends on the black-hole (BH) masses, and on the Eddington ratio. The comparison of the nature of the hosts of RL and RQ AGNs, therefore, requires pair-matching techniques. Matching in BH mass and Eddington ratio allows us to study the differences between galaxies hosting RL and RQ AGNs that have the same basic accretion parameters. We show that these two samples differ predominantly in the host-galaxy concentration index, morphological type (in the RL sample the frequency of elliptical galaxies becoming larger with increasing radio loudness), and nebular extinction (galaxies with highest radio loudness showing only low nebular extinction). Contrary to some previous studies, we find no significant difference between our radio-loud and radio-quiet samples regarding merger/interaction features.
We report the characterization of the first $62$ MaNGA Active Galactic Nuclei (AGN) hosts in the Fifth Product Launch (MPL-5) and the definition of a control sample of non-active galaxies. This control sample - comprising two galaxies for each AGN - was selected in order to match the AGN hosts in terms of stellar mass, redshift, visual morphology and inclination. The stellar masses are in the range $9.4<\log(M/M_\odot)<11.5$, and most objects have redshifts $\leq 0.08$. The AGN sample is mostly comprised of low-luminosity AGN, with only 17 nuclei with $L([OIII]\lambda 5007)\geq 3.8\times 10^{40}$ erg s$^{-1}$ (that we call "strong AGN"). The stellar population of the control sample galaxies within the inner $1$-$3$ kpc is dominated by the old ($\sim$ $4$ - $13$ Gyr) age component, with a small contribution of intermediate age ($\sim 640$-$940$ Myr) and young stars ($\leq 40$ Myr) to the total light at $5700\AA$. While the weaker AGN show a similar age distribution to that of the control galaxies, the strong AGN show an increased contribution of younger stars and a decreased contribution of older stars. Examining the relationship between the AGN stellar population properties and $L([OIII])$, we find that with increasing $L([OIII])$, the AGN exhibit a decreasing contribution from the oldest ($>4$ Gyr) stellar population relative to control galaxies, but have an increasing contribution from the younger components with ages $\sim 40$ Myr. We also find a correlation of the mean age differences (AGN - control) with $L([OIII])$, in the sense that more luminous AGN are younger than the control objects, while the low-luminosity AGN seem to be older. These results support a connection between the growth of the galaxy bulge via formation of new stars and the growth of the Supermassive Black Hole via matter accretion in the AGN phase.
In a recent work based on 3200 stacked H$\alpha$ maps of galaxies at $z \sim 1$, Nelson et al.~find evidence for `coherent star formation': the stacked SFR profiles of galaxies above (below) the 'star formation main sequence' (MS) are above (below) that of galaxies on the MS at all radii. One might interpret this result as inconsistent with highly bursty star formation and evidence that galaxies evolve smoothly along the MS rather than crossing it many times. We analyze six simulated galaxies at $z\sim1$ from the Feedback in Realistic Environments (FIRE) project in a manner analogous to the observations to test whether the above interpretations are correct. The trends in stacked SFR profiles are qualitatively consistent with those observed. However, SFR profiles of individual galaxies are much more complex than the stacked profiles: the former can be flat or even peak at large radii because of the highly clustered nature of star formation in the simulations. Moreover, the SFR profiles of individual galaxies above (below) the MS are not systematically above (below) those of MS galaxies at all radii. We conclude that the time-averaged coherent star formation evident stacks of observed galaxies is consistent with highly bursty, clumpy star formation of individual galaxies and is not evidence that galaxies evolve smoothly along the MS.
We consider a self-interacting dark matter model in which the massive dark photon mediating the self-interaction decays to light dark fermions to avoid over-closing the universe. We find that if the model is constrained to explain the stellar kinematic data for spiral galaxies and galaxy clusters, it implies the presence of dark radiation, late kinetic decoupling for dark matter, and a suppressed linear power spectrum due to dark acoustic damping. The matter power spectrum is essentially fixed by the kinetic decoupling temperature and independent of other combinations of the parameters. We find the minimum halo mass is in the range of $10^5-10^8M_{\odot}$, with the upper limit coming from the Lyman-$\alpha$ forest power spectrum measurements. The presence of dark radiation and the damping of the matter power spectrum, in tandem with the impact of self-interactions in galactic halos, makes it possible to measure the gauge coupling and masses of the dark sector particles even when signals in conventional dark matter searches are absent.
Constraints on inflationary $B$-modes using Cosmic Microwave Background polarization data commonly rely on either template cleaning or cross-spectra between maps at different frequencies to disentangle galactic foregrounds from the cosmological signal. Assumptions about how the foregrounds scale with frequency are therefore crucial to interpreting the data. Recent results from the Planck satellite collaboration claim significant evidence for a decorrelation in the polarization signal of the spatial pattern of galactic dust between 353 GHz and 217 GHz. Such a decorrelation would suppress power in the cross spectrum between high frequency maps, where the dust is strong, and lower frequency maps, where the sensitivity to cosmological $B$-modes is strongest. Alternatively, it would leave residuals in lower frequency maps cleaned with a template derived from the higher frequency maps. If not accounted for, both situations would result in an underestimate of the dust contribution and thus an upward bias on measurements of the tensor-to-scalar ratio, $r$. In this paper we revisit this measurement and find that the no-decorrelation hypothesis cannot be excluded with the Planck data. There are three main reasons for this: i) there is significant noise bias in cross spectra between Planck data splits that needs to be accounted for; ii) there is strong evidence for unknown instrumental systematics whose amplitude we estimate using alternative Planck data splits; iii) there are significant correlations between measurements in different sky patches that need to be taken into account when assessing the statistical significance. Between $\ell=55-90$ and over $72\%$ of the sky, the dust $BB$ correlation between 217 GHz and 353 GHz is $1.001^{+.004/.021}_{-.004/.000}$ ($68\%~stat./syst.$) and shows no significant trend with sky fraction.
The origin of super-massive black holes in the early universe remains poorly understood.Gravitational collapse of a massive primordial gas cloud is a promising initial process,but theoretical studies have difficulty growing the black hole fast enough.We report numerical simulations of early black hole formation starting from realistic cosmological conditions.Supersonic gas motions left over from the Big Bang prevent early gas cloud formation until rapid gas condensation is triggered in a proto-galactic halo. A protostar is formed in the dense, turbulent gas cloud, and it grows by sporadic mass accretion until it acquires 34,000 solar masses.The massive star ends its life with a catastrophic collapse to leave a black hole -- a promising seed for the formation of a monstrous black hole.
The global survey of star clusters in the Milky Way (MWSC) is a comprehensive list of 3061 objects that provides, among other parameters, distances to clusters based on isochrone fitting. The Tycho-Gaia Astrometric Solution (TGAS) catalogue, which is a part of Gaia data release 1 (Gaia DR1), delivers accurate trigonometric parallax measurements for more than 2 million stars, including those in star clusters. We compare the open cluster photometric distance scale with the measurements given by the trigonometric parallaxes from TGAS to evaluate the consistency between these values. The average parallaxes of probable cluster members available in TGAS provide the trigonometric distance scale of open clusters, while the photometric scale is given by the distances published in the MWSC. Sixty-four clusters are suited for comparison as they have more than 16 probable members with parallax measurements in TGAS. We computed the average parallaxes of the probable members and compared these to the photometric parallaxes derived within the MWSC. We find a good agreement between the trigonometric TGAS-based and the photometric MWSC-based distance scales of open clusters, which for distances less than 2.3 kpc coincide at a level of about 0.1 mas with no dependence on the distance. If at all, there is a slight systematic offset along the Galactic equator between $30^\circ$ and $160^\circ$ galactic longitude.
Links to: arXiv, form interface, find, astro-ph, recent, 1709, contact, help (Access key information)