We present a multiwavelength investigation of a region of a nearby giant molecular cloud that is distinguished by a minimal level of star formation activity. With our new 12CO(J=2-1) and 13CO(J=2-1) observations of a remote region within the middle of the California molecular cloud, we aim to investigate the relationship between filaments, cores, and a molecular outflow in a relatively pristine environment. An extinction map of the region from Herschel Space Observatory observations reveals the presence of two 2-pc-long filaments radiating from a high-extinction clump. Using the 13CO observations, we show that the filaments have coherent velocity gradients and that their mass-per-unit-lengths may exceed the critical value above which filaments are gravitationally unstable. The region exhibits structure with eight cores, at least one of which is a starless, prestellar core. We identify a low-velocity, low-mass molecular outflow that may be driven by a flat spectrum protostar. The outflow does not appear to be responsible for driving the turbulence in the core with which it is associated, nor does it provide significant support against gravitational collapse.
We present a systematic analysis of X-ray archival data of all the 29 quasars (QSOs) at $z$ > 5.5 observed so far with Chandra, XMM-Newton and Swift-XRT, including the most-distant quasar ever discovered, ULAS J1120+0641 ($z$ = 7.08). This study allows us to place constraints on the mean spectral properties of the primordial population of luminous Type 1 (unobscured) quasars. Eighteen quasars are detected in the X-ray band, and we provide spectral-fitting results for their X-ray properties, while for the others we provide upper limits to their soft (0.5-2.0 keV) X-ray flux. We measured the power-law photon index and derived an upper limit to the column density for the five quasars (J1306+0356, J0100+2802, J1030+0524, J1148+5251, J1120+0641) with the best spectra (> 30 net counts in the 0.5-7.0 keV energy range) and find that they are consistent with values from the literature and lower-redshift quasars. By stacking the spectra of ten quasars detected by Chandra in the redshift range 5.7 $\le$ $z$ $\le$ 6.1 we find a mean X-ray power-law photon index of $\Gamma = 1.92_{-0.27}^{+0.28}$ and a neutral intrinsic absorption column density of $N_H \le 10^{23}$ cm$^{-2}$. These results suggest that the X-ray spectral properties of luminous quasars have not evolved up to $z$ $\approx$ 6. We also derived the optical-X-ray spectral slopes ($\alpha_{ox}$) of our sample and combined them with those of previous works, confirming that $\alpha_{ox}$ strongly correlates with UV monochromatic luminosity at 2500 \AA . These results strengthen the non-evolutionary scenario for the spectral properties of luminous active galactic nuclei (AGN).
We identify an object previously thought to be a star in the disk of M31, J0045+41, as a background $z\approx0.215$ AGN seen through a low-absorption region of M31. We present moderate resolution spectroscopy of J0045+41 obtained using GMOS at Gemini-North. The spectrum contains features attributable to the host galaxy. We model the spectrum to estimate the AGN contribution, from which we estimate the luminosity and virial mass of the central engine. We also detect \ion{Na}{1} absorption in the Milky Way restframe. We search for evidence of periodicity using $g$-band photometry from the Palomar Transient Factory and find evidence for a 355 day period. If this signal arises due to the orbit of a binary supermassive black hole system, the system is well within the gravitational wave regime. We calculate the time until inspiral due to gravitational radiation, assuming reasonable values of the mass ratio of the two black holes. We discuss the implications of our findings and forthcoming work to identify other such interlopers.
We present a proper motion investigation of a sample of Be star candidates towards the Magellanic Clouds, which has resulted in the identification of separate populations, in the Galactic foreground and in the Magellanic background. Be stars are broadly speaking B-type stars that have shown emission lines in their spectra. In this work, we studied a sample of 2446 and 1019 Be star candidates towards the LMC and SMC respectively, taken from the literature and proposed as possible Be stars due to their variability behaviour in the OGLE-II I band. JHKs magnitudes from the IRSF catalog and proper motions from the SPM4 catalog, were obtained for 1188 and 619 LMC and SMC Be stars candidates, respectively. Color-color and vector-point diagrams were used to identify different populations among the Be star candidates. In the LMC sample, two populations with distinctive infrared colours and kinematics were found, the bluer sample is consistent with being in the LMC and the redder one with belonging to the Milky Way disk. This settles the nature of the redder sample which had been described in previous publications as a possible unknown subclass of stars among the Be candidates in the LMC. In the SMC sample, a similar but less evident result was obtained, since this apparent unknown subclass was not seen in this galaxy. We confirm that in the selection of Be stars by their variability, although generally successful, there is a higher risk of contamination by Milky Way objects towards redder B$-$V and V$-$I colors.
In this paper we investigate the statistical properties of the Tully-Fisher
relation for a sample of 32 galaxies with measured distances from the Cepheid
period-luminosity relation and/or TRGB stars.
We take advantage of panchromatic photometry in 12 bands (from FUV to 4.5
{\mu}m) and of spatially resolved HI kinematics. We use these data together
with three kinematic measures ($W^{i}_{50}$, $V_{max}$ and $V_{flat}$)
extracted from the global HI profiles or HI rotation curves, so as to construct
36 correlations allowing us to select the one with the least scatter. We
introduce a tightness parameter $\sigma_{\perp}$ of the TFr, in order to obtain
a slope-independent measure of the goodness of fit. We find that the tightest
correlation occurs when we select the 3.6 {\mu}m photometric band together with
the $V_{flat}$ parameter extracted from the HI rotation curve.
Recently, some very strong correlations between the distribution of dark matter and baryons (the dark matter-baryon relations) in galaxies with very different morphologies, masses, sizes, and gas fractions have been obtained. Some models have been suggested to explain why the dark matter contribution is fully specified by that of the baryons. In this article, we derive two analytic expressions to explain the observed dark matter-baryon relations based on the cold dark matter (CDM) model. The resultant expressions give excellent agreement with the observational data. The parameters involved in the analytic expressions are closely related to the amount of the baryon content. This model can provide a theoretical understanding of the strong correlations observed. We suggest that the observed relation represents the end product of galaxy formation.
Characterizing the evolution of the faint end of the cluster red sequence
(RS) galaxy luminosity function (GLF) with redshift is a milestone in
understanding galaxy evolution. However, the community is still divided in that
respect, hesitating between an enrichment of the RS due to efficient quenching
of blue galaxies from $z\sim1$ to present-day or a scenario in which the RS is
built at a higher redshift and does not evolve afterwards. Recently, it has
been proposed that surface brightness (SB) selection effects could possibly
solve the literature disagreement, accounting for the diminishing of the RS
faint population in ground based observations. We investigate this hypothesis
by comparing the RS GLFs of 16 CLASH clusters computed independently from
ground-based Subaru/Suprime-Cam and HST/ACS images in the redshift range
$0.187\leq z\leq0.686$. We stack individual cluster GLFs in redshift and mass
bins.
We find similar RS GLFs for space and ground based data, with a difference of
0.2$\sigma$ in the faint end parameter $\alpha$ when stacking all clusters
together and a maximum difference of 0.9$\sigma$ in the case of the high
redshift stack, demonstrating a weak dependence on the type of observations in
the probed range of redshift and mass. When considering the full sample, we
estimate $\alpha = -0.76 \pm 0.07$ and $\alpha = -0.78 \pm 0.06$ with HST and
Subaru respectively. We note a mild variation of the faint end with redshift at
a 1.7$\sigma$ and 2.6$\sigma$ significance. We investigate the effect of SB
dimming by simulating our low redshift galaxies at high redshift. We measure an
evolution in the faint end slope of less than 1$\sigma$ in this case, implying
that the observed signature is moderately larger than one would expect from SB
dimming alone, and indicating a true evolution in the faint end slope.
(Abridged...)
This work analyses the effect of the Helium content on synthetic Period-Luminosity Relations (PLRs) and Period-Wesenheit Relations (PWRs) of Cepheids and the systematic uncertainties on the derived distances that a hidden population of He-enhanced Cepheids may generate. We use new stellar and pulsation models to build a homogeneous and consistent framework to derive the Cepheid features. The Cepheid populations expected in synthetic color-magnitude diagrams of young stellar systems (from 20 Myr to 250 Myr) are computed in several photometric bands for Y = 0.25 and Y = 0.35, at a fixed metallicity (Z = 0.008). The PLRs appear to be very similar in the two cases, with negligible effects (few %) on distances, while PWRs differ somewhat, with systematic uncertainties in deriving distances as high as about 7% at log P < 1.5. Statistical effects due to the number of variables used to determine the relations contribute to a distance systematic error of the order of few percent, with values decreasing from optical to near-infrared bands. The empirical PWRs derived from multi-wavelength datasets for the Large Magellanic Cloud (LMC) is in a very good agreement with our theoretical PWRs obtained with a standard He content, supporting the evidence that LMC Cepheids do not show any He effect.
In the earliest (so-called "Class 0") phase of sunlike (low-mass) star formation, circumstellar disks are expected to form, feeding the protostars. However, such disks are difficult to resolve spatially because of their small sizes. Moreover, there are theoretical difficulties in producing such disks in the earliest phase, due to the retarding effects of magnetic fields on the rotating, collapsing material (so-called "magnetic braking"). With the Atacama Large Millimeter/submillimeter Array (ALMA), it becomes possible to uncover such disks and study them in detail. HH 212 is a very young protostellar system. With ALMA, we not only detect but also spatially resolve its disk in dust emission at submillimeter wavelength. The disk is nearly edge-on and has a radius of ~ 60 AU. Interestingly, it shows a prominent equatorial dark lane sandwiched between two brighter features, due to relatively low temperature and high optical depth near the disk midplane. For the first time, this dark lane is seen at submillimeter wavelength, producing a "hamburger"-shaped appearance that is reminiscent of the scattered-light image of an edge-on disk in optical and near infrared. Our observations open up an exciting possibility of directly detecting and characterizing small disks around the youngest protostars through high-resolution imaging with ALMA, which provides strong constraints on theories of disk formation.
A number of earth's tremor spectral peaks show a persistent narrow bandwidth incompatible with any geophysical or instrumental origin. These peaks, located at frequencies lower than a few mHz, are in principle consistent with the earth strain waves induced by monochromatic gravitational waves. Exploring this hypothesis under the current cosmological constraints yields that the tremor peaks below 2 mHz are in apparently significant coincidence with the theoretical emission of two binary systems each consisting of a small main sequence star with mass $\sim 10^{-1} M_{\odot}$, captured by Sgr A* in a close orbit.
We present hydrodynamic simulations of gas clouds inflowing from the disk to a few hundred parsec region of the Milky Way. A gravitational potential is generated to include realistic Galactic structures by using thousands of multipole expansions that describe 6.4 million stellar particles of a self-consistent Galaxy simulation. We find that a hybrid multipole expansion model, with two different basis sets and a thick disk correction, accurately reproduces the overall structures of the Milky Way. Through non-axisymmetric Galactic structures of an elongated bar and spiral arms, gas clouds in the disk inflow to the nuclear region and form a central molecular zone (CMZ)-like nuclear ring. We find that the size of the nuclear ring evolves into ~240 pc at T~1500 Myr, regardless of the initial size. For most simulation runs, the rate of gas inflow to the nuclear region is equilibrated to ~0.02 M_sun/yr. The nuclear ring is off-centered, relative to the Galactic center, by the lopsided central mass distribution of the Galaxy model, and thus an asymmetric mass distribution of the nuclear ring arises accordingly. The vertical asymmetry of the the Galaxy model also causes the nuclear ring to be tilted along the Galactic plane. During the first ~100 Myr, the vertical frequency of the gas motion is twice that of the orbital frequency, thus the projected nuclear ring shows a twisted, infinity-like shape.
We study the radio emission of the most massive galaxies in a sample of dynamically relaxed and un-relaxed galaxy groups from Galaxy and Mass Assembly (GAMA). The dynamical state of the group is defined by the stellar dominance of the brightest group galaxy, e.g. the luminosity gap between the two most luminous members, and the offset between the position of the brightest group galaxy and the luminosity centroid of the group. We find that the radio luminosity of the most massive galaxy in the group strongly depends on its environment, such that the brightest group galaxies in dynamically young (evolving) groups are an order of magnitude more luminous in the radio than those with a similar stellar mass but residing in dynamically old (relaxed) groups. This observation has been successfully reproduced by a newly developed semi-analytic model which allows us to explore the various causes of these findings. We find that the fraction of radio loud brightest group galaxies in the observed dynamically young groups is ~2 times that in the dynamically old groups. We discuss the implications of this observational constraint on the central galaxy properties in the context of galaxy mergers and the super-massive blackhole accretion rate.
CONTEXT: For more than two decades, stellar atmosphere codes have been used to derive the stellar and wind parameters of massive stars. Although they have become a powerful tool and sufficiently reproduce the observed spectral appearance, they can hardly be used for more than measuring parameters. One major obstacle is their inconsistency between the calculated radiation field and the wind stratification due to the usage of prescribed mass-loss rates and wind-velocity fields. AIMS: We present the concepts for a new generation of hydrodynamically consistent non-local thermodynamical equilibrium (non-LTE) stellar atmosphere models that allow for detailed studies of radiation-driven stellar winds. As a first demonstration, this new kind of model is applied to a massive O star. METHODS: Based on earlier works, the PoWR code has been extended with the option to consistently solve the hydrodynamic equation together with the statistical equations and the radiative transfer in order to obtain a hydrodynamically consistent atmosphere stratification. In these models, the whole velocity field is iteratively updated together with an adjustment of the mass-loss rate. RESULTS: The concepts for obtaining hydrodynamically consistent models using a comoving-frame radiative transfer are outlined. To provide a useful benchmark, we present a demonstration model, which was motivated to describe the well-studied O4 supergiant Zeta Pup. The obtained stellar and wind parameters are within the current range of literature values. CONCLUSIONS: For the first time, the PoWR code has been used to obtain a hydrodynamically consistent model for a massive O star. This has been achieved by a profound revision of earlier concepts used for Wolf-Rayet stars. The velocity field is shaped by various elements contributing to the radiative acceleration, especially in the outer wind. (...)
The survey volume of a proper motion-limited sample is typically much smaller than a magnitude-limited sample. This is because of the noisy astrometric measurements from detectors that are not dedicated for astrometric missions. In order to apply an empirical completeness correction, existing works limit the survey depth to the shallower parts of the sky that hamper the maximum potential of a survey. The number of epoch of measurement is a discrete quantity that cannot be interpolated across the projected plane of observation, so that the survey properties change in discrete steps across the sky. This work proposes a method to dissect the survey into small parts with Voronoi tessellation using candidate objects as generating points, such that each part defines a `mini-survey' that has its own properties. Coupling with a maximum volume density estimator, the new method is demonstrated to be unbiased and recovered {\sim}20% more objects than the existing method in a mock catalogue of a white dwarf-only solar neighbourhood with Pan--STARRS 1-like characteristics. Towards the end of this work, we demonstrate one way to increase the tessellation resolution with artificial generating points, which would be useful for analysis of rare objects with small number counts.
An attempt to predict the new atomic dark matter lines is done on the example of a dark lepton atom - positronium. Its Layman-alpha line with the energy near 3 GeV may be observable if the appropriate conditions realize. For this we have studied a {\gamma}-ray excess in the Center of our Galaxy. In principle, this excess may be produced by the L{\alpha} line of a dark positronium in the medium with Compton scattering. The possibility of observations of an annihilation line (E~300 TeV) of dark positronium is also predicted. Other proposals to observe the atomic dark matter are shortly described. Besides, H{\alpha} line (1.3{\mu}) of usual positronum must be observable in the direction on the Center of our Galaxy.
Spectroscopic studies have demonstrated that nearly all Galactic globular clusters (GCs) harbour multiple stellar populations with different chemical compositions. Moreover, colour-magnitude diagrams based exclusively on Str\"omgrem photometry have allowed us to identify and characterise multiple populations along the RGB of a large number of clusters. In this paper we show for the first time that Str\"omgren photometry is also very effcient at identifying multiple populations along the AGB, and demonstrate that the AGB of M3, M92, NGC362, NGC1851, and NGC6752 are not consistent with a single stellar population. We also provide a catalogue of RGB and AGB stars photometrically identified in these clusters for further spectroscopic follow-up studies.We combined photometry and elemental abundances from the literature for RGB and AGB stars in NGC6752 where the presence of multiple populations along the AGB has been widely debated. We find that, while the MS, SGB, and RGB host three stellar populations with different helium and light element abundances, only two populations of AGB stars are present in the cluster. These results are consistent with standard evolutionary theory.
The interaction of the cosmic microwave background (CMB) with the hot gas in
clusters of galaxies, the so-called Sunyaev--Zel'dovich (SZ) effect, is a very
useful tool that allows us to determine the physical conditions in such
clusters and fundamental parameters of the cosmological models. In this work,
we determine the dependence of the the SZ surface brightness amplitude with
redshift and mass of the clusters. We have used PLANCK+IRAS data in the
microwave-far infrared and a catalog with >10^5 clusters of galaxies extracted
from the SDSS by Wen et al. (2012). We estimate and subtract the dust emission
from those clusters. From the residual flux, we extract its SZ flux densities.
The absolute value of the SZ amplitude indicates that the gas mass is around
10% of the total mass for cluster masses of M~10^{14} M_sun. This amplitude is
compatible with no evolution with redshift and proportional to M^{2.70+/-0.37}
(using X-ray derived masses) or M^{2.51+/-0.38} (using weak-lensing derived
masses), with some tension regarding the expectations of the self-similar
dependence (amplitude proportional to M^{5/3}).
Other secondary products of our analysis include that clusters have a dust
emission with emissivity index beta~2 and temperature T~25 K; we confirm that
the CMB temperature agrees with a dependence of T_0(1+z) with clusters of much
lower mass than those explored previously; and we find that the cluster masses
derived by Wen et al. (2012) from a richness-mass relationship are biased by a
factor of (1+z)^{-1.8} with respect to the X-ray and weak-lensing measurements.
Links to: arXiv, form interface, find, astro-ph, recent, 1704, contact, help (Access key information)
We use multi-wavelength SDSS images and Galaxy Zoo morphologies to identify a sample of $\sim$$270$ late-type galaxies with an off-centre bar. We measure offsets in the range 0.2-2.5 kpc between the photometric centres of the stellar disc and stellar bar. The measured offsets correlate with global asymmetries of the galaxies, with those with largest offsets showing higher lopsidedness. These findings are in good agreement with predictions from simulations of dwarf-dwarf tidal interactions producing off-centre bars. We find that the majority of galaxies with off-centre bars are of Magellanic type, with a median mass of $10^{9.6} M_{\odot}$, and 91% of them having $M_{\star}<3\times10^{10} M_{\odot}$, the characteristic mass at which galaxies start having higher central concentrations attributed to the presence of bulges. We conduct a search for companions to test the hypothesis of tidal interactions, but find that a similar fraction of galaxies with offset bars have companions within 100 kpc as galaxies with centred bars. Although this may be due to the incompleteness of the SDSS spectroscopic survey at the faint end, alternative scenarios that give rise to offset bars such as interactions with dark companions or the effect of lopsided halo potentials should be considered. Future observations are needed to confirm possible low mass companion candidates and to determine the shape of the dark matter halo, in order to find the explanation for the off-centre bars in these galaxies.
We analyze the spectra of 300,000 luminous red galaxies (LRGs) with stellar masses $M_* \gtrsim 10^{11} M_{\odot}$ from the SDSS-III Baryon Oscillation Spectroscopic Survey (BOSS). By studying their star-formation histories, we find two main evolutionary paths converging into the same quiescent galaxy population at $z\sim0.55$. Fast-growing LRGs assemble $80\%$ of their stellar mass very early on ($z\sim5$), whereas slow-growing LRGs reach the same evolutionary state at $z\sim1.5$. Further investigation reveals that their clustering properties on scales of $\sim$1-30 Mpc are, at a high level of significance, also different. Fast-growing LRGs are found to be more strongly clustered and reside in overall denser large-scale structure environments than slow-growing systems, for a given stellar-mass threshold. Our results imply a dependence of clustering on stellar-mass assembly history (naturally connected to the mass-formation history of the corresponding halos) for a homogeneous population of similar mass and color, which constitutes a strong observational evidence of galaxy assembly bias.
We present the results from a broadband (1 to 3 GHz), spectro-polarimetry study of the integrated emission from 100 extragalactic radio sources with the ATCA, selected to be highly linearly polarized at 1.4 GHz. We use a general purpose, polarization model-fitting procedure that describes the Faraday rotation measure (RM) and intrinsic polarization structure of up to three distinct polarized emission regions or 'RM components' of a source. Overall, 37%/52%/11% of sources are best fit by one/two/three RM components. However, these fractions are dependent on the signal-to-noise ratio (S/N) in polarization (more RM components more likely at higher S/N). In general, our analysis shows that sources with high integrated degrees of polarization at 1.4 GHz have low Faraday depolarization, are typically dominated by a single RM component, have a steep spectral index, and a high intrinsic degree of polarization. After classifying our sample into radiative-mode and jet-mode AGN, we find no significant difference between the Faraday rotation or Faraday depolarization properties of jet-mode and radiative-mode AGN. However, there is a statistically significant difference in the intrinsic degree of polarization between the two types, with the jet-mode sources having more intrinsically ordered magnetic field structures than the radiative-mode sources. We also find a preferred perpendicular orientation of the intrinsic magnetic field structure of jet-mode AGN with respect to the jet direction, while no clear preference is found for the radiative-mode sources.
We carried out a pilot campaign of radio and optical band intra-day variability (IDV) observations of five blazars (3C66A, S5 0716+714, OJ287, B0925+504, and BL Lacertae) on December 18--21, 2015 by using the radio telescope in Effelsberg (Germany) and several optical telescopes in Asia, Europe, and America. After calibration, the light curves from both 5 GHz radio band and the optical R band were obtained, although the data were not smoothly sampled over the sampling period of about four days. We tentatively analyse the amplitudes and time scales of the variabilities, and any possible periodicity. The blazars vary significantly in the radio (except 3C66A and BL Lacertae with only marginal variations) and optical bands on intra- and inter-day time scales, and the source B0925+504 exhibits a strong quasi-periodic radio variability. No significant correlation between the radio- and optical-band variability appears in the five sources, which we attribute to the radio IDV being dominated by interstellar scintillation whereas the optical variability comes from the source itself. However, the radio- and optical-band variations appear to be weakly correlated in some sources and should be investigated based on well-sampled data from future observations.
Recent state-of-the-art calculations of A-values and electron impact excitation rates for Fe III are used in conjunction with the Cloudy modeling code to derive emission line intensity ratios for optical transitions among the fine-structure levels of the 3d$^6$ configuration. A comparison of these with high resolution, high signal-to-noise spectra of gaseous nebulae reveals that previous discrepancies found between theory and observation are not fully resolved by the latest atomic data. Blending is ruled out as a likely cause of the discrepancies, because temperature- and density-independent ratios (arising from lines with common upper levels) match well with those predicted by theory. For a typical nebular plasma with electron temperature $T_{\rm e} = 9000$ K and electron density $\rm N_{e}=10^4 \, cm^{-3}$, cascading of electrons from the levels $\rm ^3G_5$, $\rm ^3G_4$ and $\rm ^3G_3$ plays an important role in determining the populations of lower levels, such as $\rm ^3F_4$, which provide the density diagnostic emission lines of Fe III, such as $\rm ^5D_4$ - $\rm ^3F_4$ at 4658 \AA. Hence further work on the A-values for these transitions is recommended, ideally including measurements if possible. However, some Fe III ratios do provide reliable $N_{\rm e}$-diagnostics, such as 4986/4658. The Fe III cooling function calculated with Cloudy using the most recent atomic data is found to be significantly greater at $T_e$ $\simeq$ 30000 K than predicted with the existing Cloudy model. This is due to the presence of additional emission lines with the new data, particularly in the 1000--4000 \AA\ wavelength region.
Recent studies have suggested that star formation (SF) in dense gas may be governed by essentially the same law in Galactic clouds and external galaxies. This conclusion remains controversial, however, because different tracers have been used to probe the dense gas mass in Galactic and extragalactic studies. We conducted observations of Aquila, Oph., and Orion B in HCN(1-0), HCO$^+$(1-0), and their isotopomers to calibrate the HCN and HCO$^+$ lines used as dense gas tracers in extragalactic studies and to test the possible universality of the star formation efficiency in dense gas, SFE$_{dense}$. H$^{13}$CO$^+$ and H$^{13}$CN were observed to be good tracers of the filaments detected with Herschel. Comparing the luminosities $L_{HCN}$ and $L_{HCO^+}$ with the reference masses $M_{Herschel}^{Av>8}$, the empirical conversion factors $\alpha_{Herschel-HCN}$ and $\alpha_{Herschel-HCO^+}$ were found to be anti-correlated with the local FUV strength. In agreement with Pety et al. (2017), HCN and HCO$^+$ were also found to trace gas down to Av>2. As a result, published extragalactic HCN studies must be tracing all of the moderate density gas. Estimating the contribution of this moderate density gas from the typical column density PDFs in nearby clouds, we obtained the following $G_0$-dependent HCN conversion factor for external galaxies:$\alpha_{Herschel-HCN}^{fit'}=64\times G_0^{-0.34}$. Re-estimating the dense gas masses in external galaxies with $\alpha_{Herschel-HCN}^{fit'}$, we found that SFE$_{dense}$ is remarkably constant over 8 orders of magnitude in dense gas mass. Our results confirm that SFE$_{dense}$ of galaxies is quasi-universal on a wide range of scales from ~1-10pc to >10kpc. Based on the tight link between SF and filamentary structure found in Herschel studies of nearby clouds, we argue that SFE$_{dense}$ is primarily set by the microphysics of SF in filaments.
In this paper, we extend the study initiated in PaperI by modelling grain ensemble evolution in a dynamical model of an expanding HII region and checking the effects of momentum transfer from dust to gas. The radiation pressure on the dust, the dust drift, and the lug on the gas by the dust are all important process that should be considered simultaneously to describe the dynamics of HII regions. With accounting for the momentum transfer from the dust to the gas, the expansion time of the HII region is notably reduced (for our model of RCW120, the time to reach the observed radius of the HII region is reduced by a factor of 1.5). Under the common approximation of frozen dust, where there is no relative drift between the dust and gas, the radiation pressure from the ionizing star drives the formation of the very deep gas cavity near the star. Such a cavity is much less pronounced when the dust drift is taken into account. The dust drift leads to the two-peak morphology of the dust density distribution and significantly reduces the dust-to-gas ratio in the ionized region (by a factor of 2 to 10). The dust-to-gas ratio is larger for higher temperatures of the ionizing star since the dust grains have a larger electric charge and are more strongly coupled to the gas.
Recent availability of high quality infrared (IR) data for diffuse regions in the Galaxy and external galaxies have added to our understanding of interstellar dust. A comparison of ultraviolet (UV) and IR observations may be used to estimate absorption, scattering and thermal emission from interstellar dust. In this paper, we report IR and UV observations for selective diffuse sources in the Galaxy. Using archival mid-infrared (MIR) and far-infrared (FIR) observations from Spitzer Space Telescope, we look for counterparts of diffuse far-ultraviolet (FUV) sources observed by the Voyager, Far Ultraviolet Spectroscopic Explorer (FUSE) and Galaxy Evolution Explorer (GALEX) telescopes in the Galaxy. IR emission features at 8micron are generally attributed to Polycyclic Aromatic Hydrocarbon (PAH) molecules, while emission at 24micron are attributed to Very Small Grains (VSGs). The data presented here is unique and our study tries to establish a relation between various dust populations. By studying the FUV-IR correlations separately at low and high latitude locations, we have identified the grain component responsible for the diffuse FUV emission.
With a sample of 48,161 K giant stars selected from the LAMOST DR 2 catalogue, we construct torus models in a large volume extending, for the first time, from the solar vicinity to a Galactocentric distance of $\sim 20$ kpc, reaching the outskirts of the Galactic disc. We show that the kinematics of the K giant stars match conventional models, e.g. as created by Binney in 2012, in the Solar vicinity. However such two-disc models fail if they are extended to the outer regions, even if an additional disc component is utilised. If we loosen constraints in the Sun's vicinity, we find that an effective thick disc model could explain the anti-centre of the MW. The LAMOST data imply that the sizes of the Galactic discs are much larger, and that the outer disc is much thicker, than previously thought, or alternatively that the outer structure is not a conventional disc at all. However, the velocity dispersion $\sigma_{0z}$ of the kinematically thick disc in the best-fitting model is about 80 km s$^{-1}$ and has a scale parameter $R_{\sigma}$ for an exponential distribution function of $\sim 19$ kpc. Such a height $\sigma_{0z}$ is strongly rejected by current measurements in the solar neighbourhood, and thus a model beyond quasi-thermal, two or three thin or thick discs is required.
We investigate the conditions of ideal magnetohydrodynamic (MHD) turbulence responsible for the relative orientation between density structures, characterized by their gradient, $\vec{\nabla}\rho$, and the magnetic field, $\vec{B}$, in molecular clouds (MCs). For that purpose, we construct an expression for the time evolution of the angle, $\phi$, between $\vec{\nabla}\rho$ and $\vec{B}$ based on the transport equations of MHD turbulence. Using this expression, we find that the configuration where $\vec{\nabla}\rho$ and $\vec{B}$ are mostly parallel, $\cos\phi=1$, and where $\vec{\nabla}\rho$ and $\vec{B}$ are mostly perpendicular, $\cos\phi=0$, constitute attractors, that is, the system tends to evolve towards either of these configurations and they are more represented than others. This fact would explain the predominant alignment or anti-alignment between column density, $N_H$, structures and the projected magnetic field orientation, $\hat{B}_\perp$, reported in observations. Additionally, we find that departures from the $\cos\phi=0$ configurations are related to convergent flows, quantified by the divergence of the velocity field, $\vec{\nabla}\cdot\vec{v}$, in the presence of a relatively strong magnetic field. This would explain the observed change in relative orientation between $N_H$-structures and $\hat{B}_\perp$ towards MCs, from mostly parallel at low $N_H$ to mostly perpendicular at the highest $N_H$, as the result of the gravitational collapse and/or convergence of flows. Finally, we show that the density threshold that marks the observed change in relative orientation towards MCs, from $N_H$ and $\hat{B}_\perp$ being mostly parallel at low $N_H$ to mostly perpendicular at the highest $N_H$, is related to the magnetic field strength and constitutes a crucial piece of information for determining the role of the magnetic field in the dynamics of MCs.
We present an analysis of the properties of neutral hydrogen (HI) in 248 nearby (0.02<z<0.25) radio galaxies with $S_{\rm 1.4\, GHz}>30$ mJy and for which optical spectroscopy is available. The observations were carried out with the Westerbork Synthesis Radio Telescope as the last large project before the upgrade of the telescope with phased array feed receivers (Apertif). The sample covers almost four orders of magnitude in radio power from $\log\,P_{\rm 1.4 \,GHz}=22.5$ W Hz$^{-1}$ and $26.2$ W Hz$^{-1}$. We detect HI in absorption in $27\% \pm 5.5\%$ of the objects. The distribution and kinematics of the absorbing HI gas appear to depend on radio power, the properties of the radio continuum emission, and the dust content of the sources. Among the sources where HI is detected, gas with kinematics deviating from regular rotation is more likely found as the radio power increases. In these cases, the HI profile is often asymmetric with a significant blue-shifted component. This is particularly common for sources with $\log\,P_{\rm 1.4 \, GHz}>24$ W Hz$^{-1}$, where the radio emission is small, possibly because these radio sources are young. The same is found for sources that are bright in the mid-infrared, i.e. sources rich in heated dust. In these sources, the HI is outflowing likely under the effect of the interaction with the radio emission. Conversely, in dust-poor galaxies, and in sources with extended radio emission, at all radio powers we only detect HI distributed in a rotating disk. Stacking experiments show that in sources for which we do not detect HI in absorption directly, the HI has a column density that is lower than $3.5\times 10^{17} (T_{ \rm spin}/c_f)$ cm$^{-2}$. We use our results to predict the number and type of HI absorption lines that will be detected by the upcoming surveys of the Square Kilometre Array precursors and pathfinders (Apertif, MeerKAT, and ASKAP).
Deep observations of nearby galaxy clusters with Chandra have revealed concave 'bay' structures in a number of systems (Perseus, Centaurus and Abell 1795), which have similar X-ray and radio properties. These bays have all the properties of cold fronts, where the temperature rises and density falls sharply, but are concave rather than convex. By comparing to simulations of gas sloshing, we find that the bay in the Perseus cluster bears a striking resemblance in its size, location and thermal structure, to a giant ($\approx$50 kpc) roll resulting from Kelvin-Helmholtz instabilities. If true, the morphology of this structure can be compared to simulations to put constraints on the initial average ratio of the thermal and magnetic pressure, $\beta= p_{\rm th} / p_{\rm B}$, throughout the overall cluster before the sloshing occurs, for which we find $\beta=200$ to best match the observations. Simulations with a stronger magnetic field ($\beta=100$) are disfavoured, as in these the large Kelvin-Helmholtz rolls do not form, while in simulations with a lower magnetic field ($\beta=500$) the level of instabilities is much larger than is observed. We find that the bay structures in Centaurus and Abell 1795 may also be explained by such features of gas sloshing.
Previously we identified a new class of early galaxy that we estimate contributes up to 30% of the ionizing photons responsible for reionization. These are low mass halos in the range $M_{\rm halo} =10^{6.5}-10^{8} M_{\odot}$ that have been chemically enriched by supernova ejecta from prior Pop III star formation. Despite their low star formation rates, these Metal Cooling halos (MCs) are significant sources of ionizing radiation, especially at the onset of reionization, due to their high number density and ionizing escape fractions. Here we present a fully-coupled radiation hydrodynamic simulation of reionization that includes these MCs as well the more massive hydrogen atomic line cooling halos. Our method is novel: we perform halo finding inline with the radiation hydrodynamical simulation, and assign escaping ionizing fluxes to halos using a probability distribution function (PDF) measured from the galaxy-resolving Renaissance Simulations. The PDF captures the mass dependence of the ionizing escape fraction as well as the probability that a halo is actively forming stars. With MCs, reionization starts earlier than if only halos of $10^8 M_{\odot}$ and above are included, however the redshift when reionization completes is only marginally affected as this is driven by more massive galaxies. Because star formation is intermittent in MCs, the earliest phase of reionization exhibits a stochastic nature, with small HII regions forming and recombining. Only later, once halos of mass $\sim 10^9 M_{\odot}$ and above begin to dominate the ionizing emissivity, does reionization proceed smoothly in the usual manner deduced from previous studies. This occurs at $z\approx 10$ in our simulation.
We study the mass-richness relation of 116 spectroscopically-confirmed massive clusters at $0.4 < z < 2$ by mining the $Spitzer$ archive. We homogeneously measure the richness at 4.5$\mu$m for our cluster sample within a fixed aperture of $2^{\prime}$ radius and above a fixed brightness threshold, making appropriate corrections for both background galaxies and foreground stars. We have two subsamples, those which have a) literature X-ray luminosities and b) literature Sunyaev-Zeldovich effect masses. For the X-ray subsample we re-derive masses adopting the most recent calibrations. We then calibrate an empirical mass-richness relation for the combined sample spanning more than one decade in cluster mass and find the associated uncertainties in mass at fixed richness to be $\pm 0.25$ dex. We study the dependance of the scatter of this relation with galaxy concentration, defined as the ratio between richness measured within an aperture radius of 1 and 2 arcminutes. We find that at fixed aperture radius the scatter increases for clusters with higher concentrations. We study the dependance of our richness estimates with depth of the [4.5]$\mu$m imaging data and find that reaching a depth of at least [4.5]= 21 AB mag is sufficient to derive reasonable mass estimates. We discuss the possible extension of our method to the mid-infrared $WISE$ all-sky survey data, and the application of our results to the $Euclid$ mission. This technique makes richness-based cluster mass estimates available for large samples of clusters at very low observational cost.
We explore the relationship between young, embedded binaries and their parent cores, using observations within the Perseus Molecular Cloud. We combine recently published VLA observations of young stars with core properties obtained from SCUBA-2 observations at 850 um. Most embedded binary systems are found toward the centres of their parent cores, although several systems have components closer to the core edge. Wide binaries, defined as those systems with physical separations greater than 500 au, show a tendency to be aligned with the long axes of their parent cores, whereas tight binaries show no preferred orientation. We test a number of simple, evolutionary models to account for the observed populations of Class 0 and I sources, both single and binary. In the model that best explains the observations, all stars form initially as wide binaries. These binaries either break up into separate stars or else shrink into tighter orbits. Under the assumption that both stars remain embedded following binary breakup, we find a total star formation rate of 168 Myr^-1. Alternatively, one star may be ejected from the dense core due to binary breakup. This latter assumption results in a star formation rate of 247 Myr^-1. Both production rates are in satisfactory agreement with current estimates from other studies of Perseus. Future observations should be able to distinguish between these two possibilities. If our model continues to provide a good fit to other star-forming regions, then the mass fraction of dense cores that becomes stars is double what is currently believed.
Planetary nebulae are traditionally considered to represent the final evolutionary stage of all intermediate-mass stars ($\sim$0.7-8Msol). Recent evidence seems to contradict this picture. In particular, since the launch of the Hubble Space Telescope it has become clear that planetary nebulae display a wide range of striking morphologies which cannot be understood in a single star scenario, instead pointing towards a binary evolution in a majority of systems. Here, we summarise our current understanding of the importance of binarity in the formation and shaping of planetary nebulae, as well as the surprises that recent observational studies have revealed with respect to our understanding of binary evolution in general. These advances have critical implications, including for the understanding of mass transfer processes in binary stars - particularly the all-important but ever-so poorly understood `common envelope phase' - as well as the formation of cosmologically important type Ia supernovae.
We obtain a new set of analytical solutions for the evolution of a self-gravitating accretion disc by holding the Toomre parameter close to its threshold, and obtaining the stress parameter from the cooling rate. In agreement with the previous numerical solutions, furthermore, the accretion rate is assumed to be independent of the disc radius. Extreme situations where the entire disc is either optically thick or optically thin are studied independently and the obtained solutions can be used for exploring the early or the final phases of a protoplanetary disc evolution. Our solutions exhibit decay of the accretion rate as a power-law function of the age of the system with the exponent -0.75 and -1.04 for optically thick and thin cases, respectively. Our calculations permit us to explore evolution of the snow line analytically. Location of the snow line in the optically thick regime evolves as a power-law function of time with the exponent -0.16, however, when the disc is optically thin, location of the snow line as a function of time with the exponent -0.7 has a stronger dependence on time. It means that in an optically thin disc inward migration of the snow line is faster than an optically thick disc.
Triangulum-II, a newly discovered dwarf spheroidal galaxy, is a strong candidate for the indirect search of dark matter through the detection of {\gamma}-ray emission that could originate from the pair- annihilation of weakly interacting massive particles (WIMPs). We here report on the analysis of almost seven years Fermi Gamma-Ray Space Telescope data of Triangulum-II which was taken during its all sky survey operation mode. No excess {\gamma}-ray emission has been detected above 100 MeV from Triangulum-II. We derive the upper limits on {\gamma}-ray flux assuming both the power-law spectra and the spectra related to WIMP annihilation. In this work, we have considered several theoretical WIMP (neutralinos here) models envisioning both thermal and non-thermal production of WIMPs, and put limits on pair-annihilation cross-section of WIMPs to constrain the parameter space related to those theoretical models.
Links to: arXiv, form interface, find, astro-ph, recent, 1705, contact, help (Access key information)
Several studies of nearby active galactic nuclei (AGN) have shown that the soft X-ray emission presents a size and morphology that resembles that of the narrow-line region (NLR) traced by [O III]. Since the NLR is mainly constituted by gas photoionised by the AGN, it seems logical to assume that this is also the primary source of the soft X-ray emission. However, these results are based on individual sources or small samples, particularly focused on type-2 Seyfert galaxies. Very little has been said concerning other types of AGN. The purpose of this work is to compare the circumnuclear morphologies of soft X-ray and [O III] images to test whether they match in different optical classes of AGN. Our sample is composed of 27 AGN: nine type-1 Seyferts, 10 type-2 Seyferts, and eight low ionisation nuclear emission-line regions (LINERs). We find a good match in 100% of the type-2 Seyferts in our sample. This correspondence is less frequent in type-1 Seyferts (22%) and it is not seen in LINERs. The good resemblance in type-2 Seyferts constitutes an evidence for a common physical origin. We argue that the lack of correspondence in type-1 Seyferts might be due to the line of sight perpendicular to the accretion disk. Based on the morphologies of the eight LINERs in our sample, we discard a common origin for the soft X-ray and [O III] emissions in these objects. Regarding the X-ray properties, both high column density and hard X-ray luminosity are associated with matched morphologies.
Understanding the nature of bulges in disc galaxies can provide important insights into the formation and evolution of galaxies. For instance, the presence of a classical bulge suggests a relatively violent history, in contrast, the presence of simply an inner disc (also referred to as a "pseudobulge") indicates the occurrence of secular evolution processes in the main disc. However, we still lack criteria to effectively categorise bulges, limiting our ability to study their impact on the evolution of the host galaxies. Here we present a recipe to separate inner discs from classical bulges by combining four different parameters from photometric and kinematic analyses: The bulge S\'ersic index $n_\mathrm{b}$, the concentration index $C_{20,50}$, the Kormendy (1977) relation and the inner slope of the radial velocity dispersion profile $\nabla\sigma$. With that recipe we provide a detailed bulge classification for a sample of 45 galaxies from the integral-field spectroscopic survey CALIFA. To aid in categorising bulges within these galaxies, we perform 2D image decomposition to determine bulge S\'ersic index, bulge-to-total light ratio, surface brightness and effective radius of the bulge and use growth curve analysis to derive a new concentration index, $C_{20,50}$. We further extract the stellar kinematics from CALIFA data cubes and analyse the radial velocity dispersion profile. The results of the different approaches are in good agreement and allow a safe classification for approximately $95\%$ of the galaxies. In particular, we show that our new "inner" concentration index performs considerably better than the traditionally used $C_{50,90}$ when yielding the nature of bulges. We also found that a combined use of this index and the Kormendy (1977) relation gives a very robust indication of the physical nature of the bulge.
We present Herschel spectroscopy of atomic lines arising in photodissociation regions as well as ionization regions of nearby early-type galaxies (ETGs), focusing on the volume-limited Atlas3D sample. Our data include the [CII], [OI], and [NII] 122 and 205 micron lines, along with ancillary data including CO and HI maps. We find ETGs have [CII]/FIR ratios slightly lower than spiral galaxies in the KINGFISH sample, and several ETGs have unusually large [NII] 122/[CII] ratios. The [NII] 122/[CII] ratio is correlated with UV colors and there is a strong anti-correlation of [CII]/FIR with NUV-K seen in both spirals and ETGs, likely due to a softer radiation field with fewer photons available to ionize carbon and heat the gas. The correlation thus makes a [CII] deficit in galaxies with redder stellar populations. The high [NII] 122/[CII] (and low [CII]/FIR) line ratios could also be affected by the removal of much of the diffuse, low density gas, which is consistent with the low HI/H2 ratios. [CII] is now being used as a star formation indicator, and we find it is just as good for ETGs as in spirals. The [CII]/CO ratios found are also similar to those found in spiral galaxies. Through use of the [NII] 205 micron line, estimates of the percentage of [CII] emission arising from ionized gas indicate a significant portion could arise in ionized regions.
Post-starbursts (PSBs) are candidate for rapidly transitioning from star-bursting to quiescent galaxies. We study the molecular gas evolution of PSBs at z ~ 0.03 - 0.2. We undertook new CO (2-1) observations of 22 Seyfert PSBs candidates using the ARO Submillimeter Telescope. This sample complements previous samples of PSBs by including green valley PSBs with Seyfert-like emission, allowing us to analyze for the first time the molecular gas properties of 116 PSBs with a variety of AGN properties. The distribution of molecular gas to stellar mass fractions in PSBs is significantly different than normal star-forming galaxies in the COLD GASS survey. The combined samples of PSBs with Seyfert-like emission line ratios have a gas fraction distribution which is even more significantly different and is broader (~ 0.03-0.3). Most of them have lower gas fractions than normal star-forming galaxies. We find a highly significant correlation between the WISE 12 micron to 4.6 micron flux ratios and molecular gas fractions in both PSBs and normal galaxies. We detect molecular gas in 27% of our Seyfert PSBs. Taking into account the upper limits, the mean and the dispersion of the distribution of the gas fraction in our Seyfert PSB sample are much smaller (mean = 0.025, std dev. = 0.018) than previous samples of Seyfert PSBs or PSBs in general (mean ~ 0.1 - 0.2, std dev. ~ 0.1 - 0.2).
Recent observations show a population of active galaxies with milliarcseconds offsets between optical and radio emission. Such offsets can be an indication of extreme phenomena associated with supermassive black holes including relativistic jets, binary supermassive black holes, or even recoiling supermassive black holes. However, the multi-wavelength structure of active galaxies at a few milliarcseconds cannot be fathomed with direct observations. We propose using strong gravitational lensing to elucidate the multi-wavelength structure of sources. When sources are located close to the caustic of lensing galaxy, even small offset in the position of the sources results in a drastic difference in the position and magnification of mirage images. We show that the angular offset in the position of the sources can be amplified more than 50 times in the observed position of mirage images. We find that at least 8% of the observed gravitationally lensed quasars will be in the caustic configuration. The synergy between SKA and Euclid will provide an ideal set of observations for thousands of gravitationally lensed sources in the caustic configuration, which will allow us to elucidate the multi-wavelength structure for a large ensemble of sources, and study the physical origin of radio emissions, their connection to supermassive black holes, and their cosmic evolution.
We present 1-second cadence observations of M32 (NGC221) with the CHIMERA instrument at the Hale 200-inch telescope of the Palomar Observatory. Using field stars as a baseline for relative photometry, we are able to construct a light curve of the nucleus in the g-prime and r-prime band with 1sigma=36 milli-mag photometric stability. We derive a temporal power spectrum for the nucleus and find no evidence for a time-variable signal above the noise as would be expected if the nuclear black hole were accreting gas. Thus, we are unable to constrain the spin of the black hole although future work will use this powerful instrument to target more actively accreting black holes. Given the black hole mass of (2.5+/-0.5)*10^6 Msun inferred from stellar kinematics, the absence of a contribution from a nuclear time-variable signal places an upper limit on the accretion rate which is 4.6*10^{-8} of the Eddington rate, a factor of two more stringent than past upper limits from HST. The low mass of the black hole despite the high stellar density suggests that the gas liberated by stellar interactions was primarily at early cosmic times when the low-mass black hole had a small Eddington luminosity. This is at least partly driven by a top-heavy stellar initial mass function at early cosmic times which is an efficient producer of stellar mass black holes. The implication is that supermassive black holes likely arise from seeds formed through the coalescence of 3-100 Msun mass black holes that then accrete gas produced through stellar interaction processes.
We present Lya and UV-nebular emission line properties of bright Lya emitters (LAEs) at z=6-7 with a luminosity of log L_Lya/[erg s-1] = 43-44 identified in the 21-deg2 area of the SILVERRUSH early sample developed with the Subaru Hyper Suprime-Cam (HSC) survey data. Our optical spectroscopy newly confirm 21 bright LAEs with clear Lya emission, and contribute to make a spectroscopic sample of 97 LAEs at z=6-7 in SILVERRUSH. From the spectroscopic sample, we select 7 remarkable LAEs as bright as Himiko and CR7 objects, and perform deep Keck/MOSFIRE and Subaru/nuMOIRCS near-infrared spectroscopy reaching the 3sigma-flux limit of ~ 2x10^{-18} erg s-1 for the UV-nebular emission lines of He II1640, C IV1548,1550, and O III]1661,1666. Except for one tentative detection of C IV, we find no strong UV-nebular lines down to the flux limit, placing the upper limits of the rest-frame equivalent widths (EW_0) of ~2-4 A for He II, C IV, and O III] lines. Here we also investigate the VLT/X-SHOOTER spectrum of CR7 whose 6 sigma detection of He II is claimed by Sobral et al. Although two individuals and the ESO-archive service carefully re-analyze the X-SHOOTER data that are used in the study of Sobral et al., no He II signal of CR7 is detected, supportive of weak UV-nebular lines of the bright LAEs even for CR7. Spectral properties of these bright LAEs are thus clearly different from those of faint dropouts at z~7 that have strong UV-nebular lines shown in the various studies. Comparing these bright LAEs and the faint dropouts, we find anti-correlations between the UV-nebular line EW_0 and UV-continuum luminosity, which are similar to those found at z~2-3.
We present stacked average far-infrared spectra of a sample of 197 dusty, star-forming galaxies (DSFGs) at $0.005 < z < 4$ using close to 90% of the SPIRE Fourier Transform Spectrometer (FTS) extragalactic data archive from the Herschel Space Observatory based on 3.5 years of science operations. These spectra explore an observed-frame $\rm 447\,GHz-1568\,GHz$ ($\rm 191\,\mu m-671\,\mu m$) frequency (wavelength) range allowing us to observe the main atomic and molecular lines emitted by gas in the interstellar medium. The sample is sub-divided into five redshift bins at $0.005 < z < 0.05$, $0.05 < z < 0.2$, $0.2 < z < 0.5$, $0.8 < z <2$, and $2 < z < 4$. To study the dependence of observed spectral lines on total infrared luminosity, the sources in a subset of the redshift bins are stacked in luminosity bins. These stacked spectra are used to determine the average properties of the interstellar medium and dense molecular gas properties of DSFGs, in particular, the fine-structure line ([CII] 158 $\mu$m and [OI] 63 $\mu$m) luminosity ratios, and the line to far-IR luminosity ratios are used to model the gas density and radiation field strength in the photodissociation regions (PDRs). For the low-redshift sample, we additionally present the average spectral line energy distributions (SLED) of CO and $\rm{H_2O}$ rotational transitions and also consider PDR conditions based on a combination of [CI] 370 $\mu$m and 609 $\mu$m and $\rm CO (7-6)$ lines. For the high-z ($0.8 < z < 4$) sample PDR models suggest a molecular gas distribution in the presence of a radiation field that is at least a factor of 10$^3$ larger than the Milky-Way and with a neutral gas density of roughly 10$^3$ to 10$^5$ cm$^{-3}$. The corresponding PDR models for the low-z sample suggest a UV radiation field and gas density comparable to those at high-z.
A radio-far infrared correlation is constructed for a sample of blue early-type galaxies in which optical-line emissions are known to be originating in star-forming activities without any significant contribution from active galactic nuclei. A tight radio-far infrared correlation similar to that in star-forming late-type galaxies is seen in these blue early-type galaxies. The average value of the 'q' parameter is estimated as $2.35\pm0.15$ with a scatter of 0.18 dex in the radio far-infrared correlation. The tight radio-far infrared correlation implies that the radio emission is primarily from the synchrotron process, which allowed us to estimate magnetic field strengths in early-type galaxies assuming equipartition between cosmic rays and magnetic field energy densities. The average value of the magnetic field strengths in early-type galaxies is estimated as $\sim10\mu$G, which is similar to magnetic field strengths in spiral galaxies. The magnetic fields in early-type galaxies are most likely captured during a tidal encounter with a gas-rich galaxy. The subsequent evolution in terms of amplification, regularization and sustainment of fields over long periods may take place in some dynamical components such as rings.
We present results based on the observations of the blazar 3C 66A from 2005 November 06 to 2016 February 14 in the BVR and I broadbands using 1.2m telescope of the Mt. Abu InfraRed Observatory (MIRO). The source was observed on 160 nights out of which on 89 nights it was monitored for more than 1 hr to check for the presence of any intra-day variability (IDV). The blazar 3C 66A exhibited significant variations in the optical flux on short and long term time scales. However, unlike highly variable S5 0716+71, it showed IDV duty cycle of about 8% only. Our statistical studies suggest the IDV time scales ranging from $\sim$ 37 min to about 3.12 hours and, in one case, a possibility of the quasi-periodic variations with characteristic timescale of $\sim$ 1.4 hrs. The IDV amplitudes in R$-$band were found to vary from 0.03 mag to as much as 0.6 mag, with larger amplitude of variation when the source was relatively fainter. The typical rate of the flux variation was estimated to be $\sim$0.07 mag hr$^{-1}$ in both, the rising and the falling phases. However, the rates of the brightness variation as high as 1.38 mag/hr were also detected. The shortest timescale of the variation, 37 min, sets an upper limit of $6.92 \times 10^{14}$ cm on the size of the emission region and about $3.7 \times 10^8\ \mathrm{M}_{\odot}$ as an estimate of the mass of the black hole, assuming the origin of the rapid optical variability is in close vicinity of the central SMBH. The long-term study suggests a mild bluer-when-brighter behavior, typical for BL Lacs.
Nuclear spirals are ubiquitous in galaxy centers. They exist not only in strong barred galaxies but also in galaxies without noticeable bars. We use high-resolution hydrodynamic simulations to study the properties of nuclear gas spirals driven by weak bar-like and oval potentials. The amplitude of the spirals increases toward the center by a geometric effect, readily developing into shocks at small radii even for very weak potentials. The shape of the spirals and shocks depends rather sensitively on the background shear. When shear is low, the nuclear spirals are loosely wound and the shocks are almost straight, resulting in large mass inflows toward the center. When shear is high, on the other hand, the spirals are tightly wound and the shocks are oblique, forming a circumnuclear disk through which gas flows inward at a relatively lower rate. The induced mass inflow rates are enough to power black hole accretion in various types of Seyfert galaxies as well as to drive supersonic turbulence at small radii.
It is not well known what drives the chemistry of a protostellar envelope, in particular the role of the stellar mass and the outflows on its chemical enrichment. We study the chemical structure of NGC 1333 IRAS 4A in order to (i) investigate the influence of the outflows on the chemistry, (ii) constrain the age of our object, (iii) compare it with a typical high-mass protostellar envelope. In our analysis we use JCMT line mapping and HIFI pointed spectra. To study the influence of the outflow on the degree of deuteration, we compare JCMT maps of HCO+ and DCO+ with non-LTE (RADEX) models in a region that spatially covers the outflow activity of IRAS 4A. To study the envelope chemistry, we derive empirical molecular abundance profiles for the observed species using the radiative transfer code (RATRAN) and adopting a 1D dust density/temperature profile from the literature. We compare our best-fit observed abundance profiles with the predictions from the time dependent gas grain chemical code (ALCHEMIC). The CO, HCN, HNC and CN abundance require an enhanced UV field which points towards an outflow cavity. The abundances (wrt H2) are 1 to 2 orders of magnitude lower than those observed in the high mass protostellar envelope (AFGL 2591), while they are found to be similar within factors of a few with respect to CO. Differences in UV radiation may be responsible for such chemical differentiation, but temperature differences seem a more plausible explanation. The CH3OH modeled abundance profile points towards an age of > 4x10^4 yrs for IRAS 4A. The spatial distribution of H2D+ differs from that of other deuterated species, indicating an origin from a foreground colder layer (<20 K). The observed abundances can be explained by passive heating towards the high mass protostellar envelope, while the presence of UV cavity channels become more important toward the low mass protostellar envelope.
Low mass galaxies are thought to provide the bulk of the ionizing radiation necessary to reionize the Universe. The amount of photons escaping the galaxies is poorly constrained theoretically, and difficult to measure observationally. Yet it is an essential parameter of reionization models. We study in detail how ionizing radiation can leak from high redshift galaxies. For this purpose, we use a series of high resolution radiation hydrodynamics simulations, zooming on three dwarf galaxies in a cosmological context. We find that the energy and momentum input from the supernova explosions has a pivotal role in regulating the escape fraction, by disrupting dense star forming clumps, and clearing sight lines in the halo. In the absence of supernovae, photons are absorbed very locally, within the birth clouds of massive stars. We follow the time evolution of the escape fraction, and find that it can vary by more than six orders of magnitude. This explains the large scatter in the value of the escape fraction found by previous studies. This fast variability also impacts the observability of the sources of reionization: a survey even as deep as $M_{\rm UV} = -14$ would miss about half of the underlying population of Lyman-continuum emitters.
We use data on extreme radio scintillation to demonstrate that this phenomenon is associated with hot stars in the solar neighbourhood. The ionized gas responsible for the scattering is found at distances up to 1.75pc from the host star, and on average must comprise 1.E5 distinct structures per star. We detect azimuthal velocities of the plasma, relative to the host star, up to 9.7 km/s, consistent with warm gas expanding at the sound speed. The circumstellar plasma structures that we infer are similar in several respects to the cometary knots seen in the Helix, and in other planetary nebulae. There the ionized gas appears as a skin around tiny molecular clumps. Our analysis suggests that molecular clumps are ubiquitous circumstellar features, unrelated to the evolutionary state of the star. The total mass in such clumps is comparable to the stellar mass.
The galaxy NGC 4418 harbours a compact ($<20$ pc) core with a very high bolometric luminosity ($\sim10^{11}$L$_\odot$). As most of the galaxy's energy output comes from this small region, it is of interest to determine what fuels this intense activity. An interaction with VV 655 has been proposed, where gas aquired by NGC 4418 could trigger intense star formation and/or black hole accretion in the centre. We aim to constrain the interaction hypothesis by studying neutral hydrogen structures around the two galaxies. We present observations at 1.4 GHz with the Very Large Array of radio continuum as well as emission and absorption from atomic hydrogen. Gaussian distributions are fitted to observed HI emission and absorption spectra. An atomic HI bridge is seen in emission, connecting NGC 4418 to VV 655. While NGC 4418 is bright in continuum emission and seen in HI absorption, VV 655 is barely detected in the continuum but show bright HI emission (M$_\mathrm{HI}\sim10^9$ M$_\odot$). We estimate SFRs from 1.4 GHz of 3.2 M$_\odot$ yr$^{-1}$ and 0.13 M$_\odot$ yr$^{-1}$ for NGC 4418 and VV 655 respectively. Systemic HI velocities of 2202$\pm$20 km s$^{-1}$ (emission) and 2105.4$\pm$10 km s$^{-1}$ (absorption) are measured for VV 655 and NGC 4418 respectively. Redshifted HI absorption is seen towards NGC 4418, suggesting gas infall. Blueshifted HI-emission is seen north-west of NGC 4418, which we interpret as a continuation of the outflow previously discussed by Sakamoto et al. (2013). The morphology and velocity structure seen in HI is consistent with an interaction scenario, where gas was transferred from VV 655 to NGC 4418, and may fuel the activity in the centre.
We present a framework to simultaneously constrain the values and uncertainties of the strength of convective core overshooting, metallicity, extinction, distance, and age in stellar populations. We then apply the framework to archival Hubble Space Telescope observations of six stellar clusters in the Large Magellanic Cloud that have reported ages between ~1-2.5 Gyr. Assuming a canonical value of the strength of core convective overshooting, we recover the well-known age-metallicity correlation, and additional correlations between metallicity and extinction and metallicity and distance. If we allow the strength of core overshooting to vary, we find that for intermediate-aged stellar clusters, the measured values of distance and extinction are negligibly effected by uncertainties of core overshooting strength. However, cluster age and metallicity may have disconcertingly large systematic shifts when core overshooting strength is allowed to vary by more than +/- 0.05 Hp. Using the six stellar clusters, we combine their posterior distribution functions to obtain the most probable core overshooting value, 0.500 +0.016 -0.134 Hp, which is in line with canonical values.
Parallaxes for 331 classical Cepheids, 31 Type II Cepheids and 364 RR Lyrae stars in common between Gaia and the Hipparcos and Tycho-2 catalogues are published in Gaia Data Release 1 (DR1) as part of the Tycho-Gaia Astrometric Solution (TGAS). In order to test these first parallax measurements of the primary standard candles of the cosmological distance ladder, that involve astrometry collected by Gaia during the initial 14 months of science operation, we compared them with literature estimates and derived new period-luminosity ($PL$), period-Wesenheit ($PW$) relations for classical and Type II Cepheids and infrared $PL$, $PL$-metallicity ($PLZ$) and optical luminosity-metallicity ($M_V$-[Fe/H]) relations for the RR Lyrae stars, with zero points based on TGAS. The new relations were computed using multi-band ($V,I,J,K_{\mathrm{s}},W_{1}$) photometry and spectroscopic metal abundances available in the literature, and applying three alternative approaches: (i) by linear least squares fitting the absolute magnitudes inferred from direct transformation of the TGAS parallaxes, (ii) by adopting astrometric-based luminosities, and (iii) using a Bayesian fitting approach. TGAS parallaxes bring a significant added value to the previous Hipparcos estimates. The relations presented in this paper represent first Gaia-calibrated relations and form a "work-in-progress" milestone report in the wait for Gaia-only parallaxes of which a first solution will become available with Gaia's Data Release 2 (DR2) in 2018.
We performed very deep searches for 2 ground-state water transitions in 13 protoplanetary disks with the HIFI instrument on-board the Herschel Space Observatory, with integration times up to 12 hours per line. Two other water transitions that sample warmer gas were also searched for with shallower integrations. The detection rate is low, and the upper limits provided by the observations are generally much lower than predictions of thermo-chemical models with canonical inputs. One ground-state transition is newly detected in the stacked spectrum of AA Tau, DM Tau, LkCa 15, and MWC 480. We run a grid of models to show that the abundance of gas-phase oxygen needs to be reduced by a factor of at least ~100 to be consistent with the observational upper limits (and positive detections) if a dust-to-gas mass ratio of 0.01 were to be assumed. As a continuation of previous ideas, we propose that the underlying reason for the depletion of oxygen (hence the low detection rate) is the freeze-out of volatiles such as water and CO onto dust grains followed by grain growth and settling/migration, which permanently removes these gas-phase molecules from the emissive upper layers of the outer disk. Such depletion of volatiles is likely ubiquitous among different disks, though not necessarily to the same degree. The volatiles might be returned back to the gas phase in the inner disk (within about 15 AU), which is consistent with current constraints. Comparison with studies on disk dispersal due to photoevaporation indicates that the timescale for volatile depletion is shorter than that of photoevaporation.
Fermi Large Area Telescope data reveal an excess of GeV gamma rays from the direction of the Galactic Center and bulge. Several explanations have been proposed for this excess including an unresolved population of millisecond pulsars (MSPs) and self-annihilating dark matter. It has been claimed that a key discriminant for or against the MSP explanation can be extracted from the properties of the luminosity function describing this source population. Specifically, is the luminosity function of the putative MSPs in the Galactic Center consistent with that characterizing the resolved MSPs in the Galactic disk? To investigate this we have used a Bayesian Markov Chain Monte Carlo to evaluate the posterior distribution of the parameters of the MSP luminosity function describing both resolved MSPs and the Galactic Center excess. At variance with some other claims, our analysis reveals that, within current uncertainties, both data sets can be well fit with the same luminosity function.
We use the most homogeneous $Geneva$ seven-colour photometric system to derive new metallicity calibrations for early A- to K-type stars that cover both, dwarf stars and giants. The calibrations are based on several spectroscopic data sets that were merged to a common scale, and we applied them to open cluster data to obtain an additional proof of the metallicity scale and accuracy. In total, metallicities of 54 open clusters are presented. The accuracy of the calibrations for single stars is in general below 0.1 dex, but for the open cluster sample with mean values based on several stars we find a much better precision, a scatter as low as about 0.03 dex. Furthermore, we combine the new results with another comprehensive photometric data set to present a catalogue of mean metallicities for more than 3000 F and G-type dwarf stars with $\sigma \sim 0.06$ dex. The list was extended by more than 1200 hotter stars up to about 8500 K (or spectral type A3) by taking advantage of their almost reddening free characteristic in the new $Geneva$ metallicity calibrations. These two large samples are well suited as primary or secondary calibrators of other data, and we already identified about 20 spectroscopic data sets that show offsets up to about 0.4 dex.
Galactic archaeology - the study of the formation and evolution of the Milky Way by reconstructing its past from its current constituents - requires precise and accurate knowledge of stellar parameters for as many stars as possible. To achieve this a number of large spectroscopic surveys have been undertaken and are still ongoing. So far consortia carrying out the different spectroscopic surveys have used different tools to determine stellar parameters of stars from their derived effective temperatures (Teff), surface gravities (log g) and metallicities ([Fe/H]) possibly combined with photometric, astrometric, interferometric or asteroseismic information. Here we aim to homogenise the stellar characterisation by applying a unified tool to a large set of publicly available spectrophotometric data. We use spectroscopic data from a variety of large surveys combined with infra-red photometry from 2MASS and AllWISE and compare these in a Bayesian manner with PARSEC isochrones to derive probability density functions (PDFs) for stellar masses, ages and distances. We treat PDFs of pre-helium-core burning, helium-core burning and post helium-core burning solutions as well as different peaks in multi-modal PDFs (i.e. each unimodal sub-PDF) of the different evolutionary phases separately. For over 2.5 million stars we report mass, age and distance estimate for each evolutionary phase and unimodal sub-PDF. We report Gaussian, skewed Gaussian, truncated Gaussian, modified truncated exponential distribution or truncated Student's t-distribution functions to represent each sub-PDF, allowing to reconstruct detailed PDFs. Comparisons with stellar parameter estimates from the literature show good agreement within uncertainties. We present UniDAM - the unified tool applicable to spectrophotometric data of different surveys to obtain a homogenised set of stellar parameters.
(Sub-)Millimeter observations of the polarized emission of aligned aspherical dust grains enable us to study the magnetic fields within protoplanetary disk. However, the interpretation of these observations is complex. One must consider the various effects that alter the measured polarized signal, such as the shape of dust grains, the efficiency of grain alignment, the magnetic field properties, and the projection of the signal along the line of sight. We aim at analyzing observations of the polarized dust emission by disentangling the effects on the polarization signal in the context of 3D radiative transfer simulations. For this purpose, we developed a code capable of simulating dust grain alignment of aspherical grains and intrinsical polarization of thermal dust emission. We find that the influence of thermal polarization and dust grain alignment on the polarized emission displayed as spatially resolved polarization map or as spectral energy distribution trace disk properties which are not traced in total (unpolarized) emission such as the magnetic field topology. The radiative transfer simulations presented in this work enable the 3D analysis of intrinsically polarized dust emission - observed with, e.g., ALMA - which is essential to constrain magnetic field properties.
Links to: arXiv, form interface, find, astro-ph, recent, 1705, contact, help (Access key information)
We investigate the clustering properties of $\sim 7000$ H$\beta$+[OIII] and [OII] narrowband-selected emitters at $z \sim 0.8 - 4.7$ from the High-$z$ Emission Line Survey. We find clustering lengths, $r_0$, of $1.5 - 4.0h^{-1}$ Mpc and minimum dark matter halo masses of $10^{10.7 - 12.1}\rm{M}_\odot$ for our $z = 0.8 - 3.2$ H$\beta$+[OIII] emitters and $r_0 \sim 2.0 - 8.3h^{-1}$ Mpc and halo masses of $10^{11.5 - 12.6}\rm{M}_\odot$ for our $z = 1.5 - 4.7$ [OII] emitters. We find $r_0$ to strongly increase both with increasing line luminosity and redshift. By taking into account the evolution of the characteristic line luminosity, $L^\star(z)$, and using our model predictions of halo mass given $r_0$, we find a strong, redshift-independent increasing trend between $L/L^\star(z)$ and minimum halo mass. The faintest H$\beta$+[OIII] emitters are found to reside in $10^{9.5}\rm{M}_\odot$ halos and the brightest emitters in $10^{13.0}\rm{M}_\odot$ halos. For [OII] emitters, the faintest emitters are found in $10^{10.5} \rm{M}_\odot$ halos and the brightest emitters in $10^{12.6}\rm{M}_\odot$ halos. A redshift-independent stellar mass dependency is also observed where the halo mass increases from $10^{11}\rm{M}_\odot$ to $10^{12.5} \rm{M}_\odot$ for stellar masses of $10^{8.5}\rm{M}_\odot$ to $10^{11.5}\rm{M}_\odot$, respectively. We investigate the interdependencies of these trends by repeating our analysis in a $L_\textrm{line} - \rm{M}_\textrm{star}$ grid space for our most populated samples (H$\beta$+[OIII] $z = 0.84$ and [OII] $z = 1.47$) and find that the line luminosity dependency is stronger than the stellar mass dependency on halo mass. For $L > L^\star$ emitters at all epochs, we find a relatively flat trend with halo masses of $10^{12.5 - 13}\rm{M}_\odot$ which may be due to quenching mechanisms in massive halos which is consistent with a transitional halo mass predicted by models.
M87, the active galaxy at the center of the Virgo cluster, is ideal for studying the interaction of a supermassive black hole (SMBH) with a hot, gas-rich environment. A deep Chandra observation of M87 exhibits an approximately circular shock front (13 kpc radius, in projection) driven by the expansion of the central cavity (filled by the SMBH with relativistic radio-emitting plasma) with projected radius $\sim$1.9 kpc. We combine constraints from X-ray and radio observations of M87 with a shock model to derive the properties of the outburst that created the 13 kpc shock. Principal constraints for the model are 1) the measured Mach number ($M$$\sim$1.2), 2) the radius of the 13 kpc shock, and 3) the observed size of the central cavity/bubble (the radio-bright cocoon) that serves as the piston to drive the shock. We find an outburst of $\sim$5$\times$$10^{57}$ ergs that began about 12 Myr ago and lasted $\sim$2 Myr matches all the constraints. In this model, $\sim$22% of the energy is carried by the shock as it expands. The remaining $\sim$80% of the outburst energy is available to heat the core gas. More than half the total outburst energy initially goes into the enthalpy of the central bubble, the radio cocoon. As the buoyant bubble rises, much of its energy is transferred to the ambient thermal gas. For an outburst repetition rate of about 12 Myrs (the age of the outburst), 80% of the outburst energy is sufficient to balance the radiative cooling.
We present the Pristine survey, a new narrow-band photometric survey focused on the metallicity-sensitive Ca H & K lines and conducted in the northern hemisphere with the wide-field imager MegaCam on the Canada-France-Hawaii Telescope (CFHT). This paper reviews our overall survey strategy and discusses the data processing and metallicity calibration. Additionally we review the application of these data to the main aims of the survey, which are to gather a large sample of the most metal-poor stars in the Galaxy, to further characterise the faintest Milky Way satellites, and to map the (metal-poor) substructure in the Galactic halo. The current Pristine footprint comprises over 1,000 deg2 in the Galactic halo ranging from b~30 to 78 and covers many known stellar substructures. We demonstrate that, for SDSS stellar objects, we can calibrate the photometry at the 0.02-magnitude level. The comparison with existing spectroscopic metallicities from SDSS/SEGUE and LAMOST shows that, when combined with SDSS broad-band g and i photometry, we can use the CaHK photometry to infer photometric metallicities with an accuracy of ~0.2 dex from [Fe/H]=-0.5 down to the extremely metal-poor regime ([Fe/H]<-3.0). After the removal of various contaminants, we can efficiently select metal-poor stars and build a very complete sample with high purity. The success rate of uncovering [Fe/H]SEGUE<-3.0 stars among [Fe/H]Pristine<-3.0 selected stars is 24% and 85% of the remaining candidates are still very metal poor ([Fe/H]<-2.0). We further demonstrate that Pristine is well suited to identify the very rare and pristine Galactic stars with [Fe/H]<-4.0, which can teach us valuable lessons about the early Universe.
We present constraints on variations in the initial mass function (IMF) of nine local early-type galaxies based on their low mass X-ray binary (LMXB) populations. Comprised of accreting black holes and neutron stars, these LMXBs can be used to constrain the important high mass end of the IMF. We consider the LMXB populations beyond the cores of the galaxies ($>0.2R_{e}$; covering $75-90\%$ of their stellar light) and find no evidence for systematic variations of the IMF with velocity dispersion ($\sigma$). We reject IMFs which become increasingly bottom heavy with $\sigma$, up to steep power-laws (exponent, $\alpha>2.8$) in massive galaxies ($\sigma>300$km/s), for galactocentric radii $>1/4\ R_{e}$. Previously proposed IMFs that become increasingly bottom heavy with $\sigma$ are consistent with these data if only the number of low mass stars $(<0.5M_{\odot}$) varies. We note that our results are consistent with some recent work which proposes that extreme IMFs are only present in the central regions of these galaxies. We also consider IMFs that become increasingly top-heavy with $\sigma$, resulting in significantly more LMXBs. Such a model is consistent with these observations, but additional data are required to significantly distinguish between this and an invariant IMF. For six of these galaxies, we directly compare with published IMF mismatch parameters from the Atlas3D survey, $\alpha_{dyn}$. We find good agreement with the LMXB population if galaxies with higher $\alpha_{dyn}$ have more top-heavy IMFs -- although we caution that our sample is quite small. Future LMXB observations can provide further insights into the origin of $\alpha_{dyn}$ variations.
Because the timescale of H$\alpha$ emission (several tens of Myr) following star formation is significantly shorter than that of UV radiation (a few hundred Myr), the H$\alpha$/UV flux ratio provides insight on the star formation histories (SFHs) of galaxies on timescales shorter than $\sim100$ Myr. We present H$\alpha$/UV ratios for galaxies at $z=$ 0.02--0.1 on the familiar star-forming main sequence based on the AKARI-GALEX-SDSS archive dataset. The data provide us with robust measurements of dust-corrected SFRs in both H$\alpha$ and UV for 1,050 galaxies. The results show a correlation between the H$\alpha$/UV ratio and the deviation from the main sequence in the sense that galaxies above/below the main sequence tend to have higher/lower H$\alpha$/UV ratios. This trend increases the dispersion of the main sequence by 0.04 dex (a small fraction of the total scatter of 0.36 dex), suggesting that diversity of recent SFHs of galaxies has a direct impact on the observed main sequence scatter. We caution that the results suffer from incompleteness and a selection bias which may lead us to miss many sources with high H$\alpha$/UV ratios; this could further increase the scatter from SFHs in the star-forming main sequence.
The drift and diffusion coefficients of the inhomogeneous multi-mass degenerate Landau equation are computed to describe the self-induced resonant relaxation of a discrete self-gravitating quasi-Keplerian razor-thin axisymmetric disc orbiting a massive black hole while relying on Gauss' method. For a disc-like configuration in our Galactic centre, secular diffusion induces an adiabatic distortion of orbits. When considering a disc composed of multiple masses similarly distributed, the population of lighter stars will gain eccentricity, driving it closer to the central black hole provided the distribution function increases with angular momentum. The quenching of the diffusion of a test star in the vicinity of the black hole due to the divergence of the relativistic precessions (the "Schwarzschild barrier") is correctly recovered by the kinetic equation. The dual stochastic Langevin formulation yields consistent results and provides a versatile framework in which to incorporate other stochastic processes.
Deep narrow-band imaging of the iconic spiral galaxy M101 has revealed over a thousand Wolf Rayet (WR) candidates. We report spectrographic confirmation of 10 HeII emission line sources hosting 15 WR stars. We find WR stars present at both sub-- and super--solar metalicities with WC stars favouring more metal-rich regions compared to WN stars. We investigate the association of WR stars with HII regions using archival HST imaging and conclude that the majority of WR stars are in or associated with HII regions with only 10% truly isolated. Our spectroscopic survey provides confidence that our narrow-band photometric candidates are in fact bonafide WR stars, which will allow us to characterise the progenitors of core-collapse supernovae that erupt in the future in M101.
We test the impact of using variable star forming histories (SFHs) and the use of the IR luminosity (LIR) as a constrain on the physical parameters of high redshift dusty star-forming galaxies. We explore in particular the stellar properties of galaxies in relation with their location on the SFR-M* diagram. We perform SED fitting of the UV-NIR and FIR emissions of a large sample of GOODS-Herschel galaxies, for which rich multi-wavelength observations are available. We test different SFHs and imposing energy conservation in the SED fitting process, to face issues like the age-extinction degeneracy and produce SEDs consistent with observations. Our models work well for the majority of the sample, with the notable exception of the high LIR end, for which we have indications that our simple energy conservation approach cannot hold true. We find trends in the SFHs fitting our sources depending on stellar mass M* and z. Trends also emerge in the characteristic timescales of the SED models depending on the location on the SFR-M* diagram. We show that whilst using the same available observational data, we can produce galaxies less star-forming than usually inferred, if we allow declining SFHs, while properly reproducing their observables. These sources can be post-starbursts undergoing quenching, and their SFRs are potentially overestimated if inferred from their LIR. Fitting without the IR constrain leads to a strong preference for declining SFHs, while its inclusion increases the preference of rising SFHs, more so at high z, in tentative agreement with the cosmic star formation history. Keeping in mind that the sample is biased towards high LIR, the evolution shaped by our model appears as both bursty (initially) and steady-lasting (later on). The global SFH of the sample follows the cosmic SFH with a small scatter, and is compatible with the "downsizing" scenario of galaxy evolution.
We present the Ly$\alpha$ luminosity functions (LFs) at $z=$5.7 and 6.6 derived from a new large sample of 2,354 Ly$\alpha$ emitters (LAEs) identified in total areas of 14 and 21 deg$^2$, respectively, based on the early narrowband data of the Subaru/Hyper Suprime-Cam (HSC) survey. Together with careful Monte-Carlo simulations that account for the incompleteness of the LAE selection and the flux estimate systematics in the narrowband imaging, we have determined the Ly$\alpha$ LFs with the unprecedentedly small statistical and systematic uncertainties in a wide Ly$\alpha$ luminosity range of $10^{42.8-43.8}$ erg s$^{-1}$. We obtain the best-fit Schechter parameters of $L^{*}_{{\rm Lya}}=1.1^{+0.7}_{-0.4}\times10^{43}\ (8.2^{+8.6}_{-3.0}\times10^{42})$ erg s$^{-1}$, $\phi^{*}_{{\rm Lya}}=2.5^{+3.4}_{-1.9}\ (2.8^{+3.6}_{-2.3})\times10^{-4}$ Mpc$^{-3}$, and $\alpha=-2.3^{+0.8}_{-0.4}\ (-1.9^{+0.8}_{-0.6})$ at $z=5.7$ ($6.6$). We confirm that our best-estimate Ly$\alpha$ LFs are consistent with the majority of the previous studies, but find that our Ly$\alpha$ LFs do not agree with the large Ly$\alpha$ LFs recently claimed by Matthee/Santos et al.'s studies that may overcorrect the incompleteness and the flux systematics. Our Ly$\alpha$ LFs at $z=$5.7 and 6.6 show a very steep slope ($\alpha\simeq-2$) or a significant ($\simeq3\sigma$) excess beyond the best-fit Schechter function at the bright end ($\gtrsim10^{43.5}$ erg s$^{-1}$) that is claimed by previous studies. This bright-end LF excess at $z\sim6-7$ may be explained by the contribution from AGNs, blended merging galaxies, and/or large ionized bubbles around bright LAEs. Comparing our Ly$\alpha$ LF measurements with four independent reionization models, we estimate the neutral hydrogen fraction of the IGM to be $x_{\rm HI}=0.3\pm0.2$ at $z=$6.6 that is consistent with the small Thomson scattering optical depth obtained by Planck 2016.
The integrated properties of nuclear rings are correlated with their host
galaxy's secular evolution and its dynamics, as well as with the formation and
evolution of the ring's star cluster population(s). Here we present a new
method to accurately measure the spectral energy distribution and current
star-formation rate (SFR) of the nuclear ring in the barred spiral galaxy NGC
1512 based on high-resolution {\sl Hubble} and {\sl Spitzer Space Telescope}
images. Image degradation does not have a significant negative effect on the
robustness of the results. To obtain the ring's SFR for the period spanning
$\sim$3--10 Myr, we apply our method to the continuum-subtracted H$\alpha$ and
8 $\mu$m images. The resulting SFR surface density, $\Sigma_{\rm
SFR}$=$0.09\, {M}_{\odot} \,{\rm yr}^{-1}$ ${\rm kpc}^{-2}$, which is much
higher than the disk-averaged SFR densities in normal galaxies. We also
estimate the ring's total stellar mass, log (${M}/{M}_{\odot}$) = 7.1 $\pm$
0.11 for an average age of $\sim$40 Myr.
We consider a sample of 412 galaxies with radial velocities $V_{\rm LG} < 2500$ km s$^{-1}$ situated in the sky region of ${\rm RA}=13^h\hspace{-0.4em}.\,0$ ... $19^h\hspace{-0.4em}.\,0$, ${\rm Dec}=+10^{\circ}$ ... $+40^{\circ}$ between the Local Void and the Supergalactic plane. One hundred and eighty-one of them have individual distance estimates. Peculiar velocities of the galaxies as a function of Supergalactic latitude SGB show signs of Virgocentric infall at $SGB < 10^{\circ}$ and motion from the Local Void at $SGB > 60^{\circ}$. A half of the Hercules-Bootes galaxies belong to 17 groups and 29 pairs, with the richest group around NGC5353. A typical group is characterized by the velocity dispersion of $67$ km s$^{-1}$, the harmonic radius of $182$ kpc, the stellar mass of $4.3 \times10^{10} M_{\odot}$ and the virial-to-stellar mass ratio of $32$. The binary galaxies have the mean radial velocity difference of $37$ km s$^{-1}$, the projected separation of $96$ kpc, the mean integral stellar mass of $2.6\times 10^9 M_{\odot}$ and the mean virial-to-stellar mass ratio of about $8$. The total dark-matter-to-stellar mass ratio in the considered sky region amounts to $37$ being almost the same as that in the Local Volume.
The Galactic habitable zone is defined as the region with highly enough metallicity to form planetary systems in which Earth-like planets could be born and might be capable of sustaining life surviving to the destructive effects of nearby supernova explosion events. Galactic chemical evolution models can be useful tools for studying the galactic habitable zones in different systems. Our aim here is to find the Galactic habitable zone using chemical evolution models for the Milky Way disc, adopting the most recent prescriptions for the evolution of dust and for the probability of finding planetary systems around M and FGK stars. Moreover, for the first time, we will express those probabilities in terms of the dust-to-gas ratio of the ISM in the solar neighborhood as computed by detailed chemical evolution models. At a fixed Galactic time and Galactocentric distance we determine the number of M and FGK stars having Earths (but no gas giant planets) which survived supernova explosions, using the formalism of our Paper I. The probabilities of finding terrestrial planets but not gas giant planets around M stars deviate substantially from the ones around FGK stars for supersolar values of [Fe/H]. For both FGK and M stars the maximum number of stars hosting habitable planets is at 8 kpc from the Galactic Centre, if destructive effects by supernova explosions are taken into account. At the present time the total number of M stars with habitable planets are $\simeq$ 10 times the number of FGK stars. Moreover, we provide a sixth order polynomial fit (and a linear one but more approximated) for the relation found with chemical evolution models in the solar neighborhood between the [Fe/H] abundances and the dust-to-gas ratio.
Most old globular clusters (GCs) in the Galaxy are observed to have internal chemical abundance spreads in light elements. We discuss a new GC formation scenario based on hierarchical star formation within fractal molecular clouds. In the new scenario, a cluster of bound and unbound star clusters (`star cluster complex', SCC) that have a power-law cluster mass function with a slope (beta) of 2 is first formed from a massive gas clump developed in a dwarf galaxy. Such cluster complexes and beta=2 are observed and expected from hierarchical star formation. The most massive star cluster (`main cluster'), which is the progenitor of a GC, can accrete gas ejected from asymptotic giant branch (AGB) stars initially in the cluster and other low-mass clusters before the clusters are tidally stripped or destroyed to become field stars in the dwarf. The SCC is initially embedded in a giant gas hole created by numerous supernovae of the SCC so that cold gas outside the hole can be accreted onto the main cluster later. New stars formed from the accreted gas have chemical abundances that are different from those of the original SCC. Using hydrodynamical simulations of GC formation based on this scenario, we show that the main cluster with the initial mass as large as [2-5]x10^5 Msun can accrete more than 10^5 Msun gas from AGB stars of the SCC. We suggest that merging of hierarchical star cluster complexes can play key roles in stellar halo formation around GCs and self-enrichment processes of GCs.
This paper estimates the specific accretion-rate distribution of AGN using a sample of 4821 X-ray sources from both deep and shallow surveys. The specific accretion-rate distribution is defined as the probability of a galaxy with a given stellar mass and redshift hosting an active nucleus with a certain specific accretion rate. We find that the probability of a galaxy hosting an AGN increases with decreasing specific accretion rate. There is evidence that this trend reverses at low specific accretion rates, $\lambda<10^{-4}-10^{-3}$ (in Eddington units). There is also a break close to the Eddington limit, above which the probability of an accretion event decreases steeply. The specific accretion-rate distribution evolves such that the fraction of AGN among galaxies drops toward lower redshifts. This decrease in the AGN duty cycle is responsible for the strong evolution of the accretion density of the Universe from redshift $z\approx1-1.5$ to the present day. Our analysis also suggests that this evolution is accompanied by a decoupling of accretion events onto black holes from the formation of stars in galaxies. There is also evidence that at earlier times the relative probability of high vs low specific accretion-rate events among galaxies increases. We argue that this differential redshift evolution of the AGN duty cycle with respect to $\lambda$ produces the AGN downsizing trend, whereby luminous sources peak at earlier epochs compared to less luminous ones. Finally, we also find a stellar-mass dependence of the specific accretion-rate distribution, with more massive galaxies avoiding high specific accretion-rate events.
We report the discovery of a new Herbig-Haro jet, HH 1165, in SOAR narrow-band imaging of the vicinity of the sigma Orionis cluster. HH 1165 shows a spectacular extended and collimated spatial structure, with a projected length of 0.26 pc, a bent C-shaped morphology, multiple knots, and fragmented bow-shocks at the apparent ends of the flow. The Halpha image shows a bright halo with a clumpy distribution of material seen around the driving source, and curved reflection nebulosity tracing the outflow cavities. The driving source of HH 1165 is a Class I proto-brown dwarf, Mayrit 1701117 (M1701117), with a total (dust+gas) mass of ~36 MJup and a bolometric luminosity of ~0.1 Lsun. High-resolution VLT/UVES spectra of M1701117 show a wealth of emission lines indicative of strong outflow and accretion activity. SOAR/Goodman low-resolution spectra along the jet axis show an asymmetrical morphology for HH 1165. We find a puzzling picture wherein the north-west part exhibits a classical HH jet running into a pre-dominantly neutral medium, while the southern part resembles an externally irradiated jet. The C-shaped bending in HH 1165 may be produced by the combined effects from the massive stars in the ionization front to the east, the sigma Orionis core to the west, and the close proximity to the B2-type star HR 1950. HH 1165 shows all of the signatures to be considered as a scaled-down version of parsec-length HH jets, and can be termed as the first sub-stellar analog of a protostellar HH jet system.
We present the results of a search for extremely metal-poor (EMP), carbon-enhanced metal-poor (CEMP), and cataclysmic variable (CV) stars using a new exploration tool based on linked scatter plots (LSPs). Our approach is especially designed to work with very large spectrum data sets such as the SDSS, LAMOST, RAVE, and Gaia data sets, and it can be applied to stellar, galaxy, and quasar spectra. As a demonstration, we conduct our search using the SDSS DR10 data set. We first created a 3326-dimensional phase space containing nearly 2 billion measures of the strengths of over 1600 spectral features in 569,738 SDSS stars. These measures capture essentially all the stellar atomic and molecular species visible at the resolution of SDSS spectra. We show how LSPs can be used to quickly isolate and examine interesting portions of this phase space. To illustrate, we use LSPs coupled with cuts in selected portions of phase space to extract EMP stars, CEMP stars, and CV stars. We present identifications for 59 previously unrecognized candidate EMP stars and 11 previously unrecognized candidate CEMP stars. We also call attention to 2 candidate He~II emission CV stars found by the LSP approach that have not yet been discussed in the literature.
Galactic rotation curves have proven to be the testing ground for dark matter bounds in spiral galaxies of all morphologies. Dwarf galaxies serve as an increasingly interesting case of rotation curve dynamics due to their typically rising rotation curve as opposed to the flattening curve of large spirals. These galaxies usually vary in galactic structure and mostly terminate at small radial distances. This, coupled with the fact that Cold Dark Matter theories struggle with the universality of galactic rotation curves, allow for exclusive features of alternative gravitational models to be analyzed. Recently, the THINGS (The HI Nearby Galactic Survey) has been extended to include a sample of 25 dwarf galaxies now known as the LITTLE THINGS Survey. Here, we present a thorough application of alternative gravitational models to the LITTLE THINGS survey, specifically focusing on conformal gravity and Modified Newtonian Dynamics. An analysis and discussion of the results of the fitting procedure of the two alternative gravitational models are explored, as well as the resulting rotation curve predictions of each. Further, we show how these two alternative gravitational models account for the recently observed universal trends in centripetal accelerations in spiral galaxies. We posit here that both conformal gravity and MOND can provide an accurate description of the galactic dynamics without the need for dark matter.
High-resolution Very-Long-Baseline Interferometry observations of active galactic nuclei have revealed asymmetric structures in the jets of radio galaxies. These asymmetric structures may be due to internal asymmetries in the jet, could be induced by the different conditions in the surrounding ambient medium including the obscuring torus, or a combination of the two. In this paper we investigate the influence of the ambient medium (including the obscuring torus) on the observed properties of jets from radio galaxies. We performed special-relativistic hydrodynamic (RHD) simulations of over-pressured and pressure-matched jets using the special-relativistic hydrodynamics code \texttt{Ratpenat}, which is based on a second-order accurate finite-volume method and an approximate Riemann solver. Using a newly developed emission code to compute the electromagnetic emission, we have investigated the influence of different ambient medium and torus configurations on the jet structure and subsequently computed the non-thermal emission produced by the jet and the thermal absorption due to the torus. To better compare the emission simulations with observations we produced synthetic radio maps, taking into account the properties of the observatory. The detailed analysis of our simulations shows that the observed asymmetries can be produced by the interaction of the jet with the ambient medium and by the absorption properties of the obscuring torus.
We have obtained low-resolution optical (0.7-0.98 micron) and near-infrared (1.11-1.34 micron and 0.8-2.5 micron) spectra of twelve isolated planetary-mass candidates (J = 18.2-19.9 mag) of the 3-Myr sigma Orionis star cluster with a view to determining the spectroscopic properties of very young, substellar dwarfs and assembling a complete cluster mass function. We have classified our targets by visual comparison with high- and low-gravity standards and by measuring newly defined spectroscopic indices. We derived L0-L4.5 and M9-L2.5 using high- and low-gravity standards, respectively. Our targets reveal clear signposts of youth, thus corroborating their cluster membership and planetary masses (6-13 Mjup). These observations complete the sigma Orionis mass function by spectroscopically confirming the planetary-mass domain to a confidence level of $\sim$75 percent. The comparison of our spectra with BT-Settl solar metallicity model atmospheres yields a temperature scale of 2350-1800 K and a low surface gravity of log g ~ 4.0 [cm/s2], as would be expected for young planetary-mass objects. We discuss the properties of the cluster least-massive population as a function of spectral type. We have also obtained the first optical spectrum of S Ori 70, a T dwarf in the direction of sigma Orionis. Our data provide reference optical and near-infrared spectra of very young L dwarfs and a mass function that may be used as templates for future studies of low-mass substellar objects and exoplanets. The extrapolation of the sigma Orionis mass function to the solar neighborhood may indicate that isolated planetary-mass objects with temperatures of 200-300 K and masses in the interval 6-13-Mjup may be as numerous as very low-mass stars.
Links to: arXiv, form interface, find, astro-ph, recent, 1705, contact, help (Access key information)
We reveal multiple components of an interacting galaxy system at $z\approx 3.35$ through a detailed analysis of the exquisite high-resolution Keck/HIRES spectrum of the afterglow of a gamma-ray burst (GRB). Through Voigt-profile fitting of absorption lines from the Lyman-series, we constrain the neutral hydrogen column density to $N_{\mathrm{HI}} \leq 10^{18.35}$ cm$^{-2}$ for the densest of four distinct systems at the host redshift of GRB 080810, among the lowest $N_{\mathrm{HI}}$ ever observed in a GRB host, despite the line of sight passing within a projected 5 kpc of the galaxy centres. By detailed analysis of the corresponding metal absorption lines, we derive chemical, ionic and kinematic properties of the individual absorbing systems, and thus build a picture of the host as a whole. Striking differences between the systems imply that the line of sight passes through several phases of gas: the star-forming regions of the GRB host; enriched material in the form of a galactic outflow; the hot and ionised halo of a second, interacting galaxy falling towards the host at a line-of-sight velocity of 700 km s$^{-1}$; and a cool, metal-poor cloud which may represent one of the best candidates yet for the inflow of metal-poor gas from the intergalactic medium.
We present an imaging and spectral analysis of the nuclear region of the ULIRG merger Arp 220, using deep \textit{Chandra}-ACIS observations summing up to \(\sim 300\mbox{ ks}\). Narrow-band imaging with sub-pixel resolution of the innermost nuclear region reveals two distinct Fe-K emitting sources, coincident with the infrared and radio nuclear clusters. These sources are separated by 1' (\(\sim 380\) pc). The X-ray emission is extended and elongated in the eastern nucleus, like the disk emission observed in millimeter radio images, suggesting starburst dominance in this region. We estimate Fe-K equivalent width \(\gtrsim 1\) keV for both sources, and observed 2-10 keV luminosities \(\sim 2\times{10}^{40}\mbox{ erg}\mbox{ s}^{-1}\) (W) and \(\sim 3 \times {10}^{40}\mbox{ erg}\mbox{ s}^{-1}\) (E). In the 6-7 keV band the emission from these regions is dominated by the 6.7 keV Fe \textsc{xxv} line, suggesting contribution from collisionally ionized gas. The thermal energy content of this gas is consistent with kinetic energy injection in the interstellar medium by Type II SNe. However, nuclear winds from hidden AGN (\(\varv\sim 2000 \mbox{ km}\mbox{ s}^{-1}\)) cannot be excluded. The \(3\sigma\) upper limits on the neutral Fe-K\(\alpha\) flux of the nuclear regions correspond to intrinsic AGN 2-10 keV luminosities \(< 1\times {10}^{42}\mbox{ erg}\mbox{ s}^{-1}\) (W) and \(< 0.4\times {10}^{42}\mbox{ erg}\mbox{ s}^{-1}\) (E). For typical AGN SEDs the bolometric luminosities are \(< 3\times {10}^{43}\mbox{ erg}\mbox{ s}^{-1}\) (W) and \(< 8\times {10}^{43}\mbox{ erg}\mbox{ s}^{-1}\) (E), and black hole masses \(<1\times{10}^5 M_{\astrosun}\) (W) and \(< 5\times{10}^5 M_{\astrosun}\) (E) for Eddington limited AGNs with a standard 10\% efficiency.
I present a simple phenomenological model for the observed linear scaling of the stellar mass in old globular clusters (GCs) with $z=0$ halo mass in which the stellar mass in GCs scales linearly with progenitor halo mass at $z=6$ above a minimum halo mass for GC formation. This model reproduces the observed $M_{\rm GCs}-M_{\rm halo}$ relation at $z=0$ and results in a prediction for the minimum halo mass at $z=6$ required for hosting one GC: $M_{\rm min}(z=6)=1.07 \times 10^9\,\msun$. Translated to $z=0$, the mean threshold mass is $M_{\rm halo}(z=0) \approx 2\times 10^{10}\,\msun$. I explore the observability of GCs in the reionization era and their contribution to cosmic reionization, both of which depend sensitively on the (unknown) ratio of GC birth mass to present-day stellar mass, $\xi$. Based on current detections of $z \ga 6$ objects with $M_{1500} < -17$, values of $\xi > 10$ are strongly disfavored; this, in turn, has potentially important implications for GC formation scenarios. Even for low values of $\xi$, some observed high-$z$ galaxies may actually be GCs, complicating estimates of reionization-era galaxy ultraviolet luminosity functions and constraints on dark matter models. GCs are likely important reionization sources if $5 \la \xi \la 10$. I also explore predictions for the fraction of accreted versus in situ GCs in the local Universe and for descendants of systems at the halo mass threshold of GC formation (dwarf galaxies). An appealing feature of the model presented here is the ability to make predictions for GC properties based solely on dark matter halo merger trees.
We study the observed relation between accretion rate (in terms of L/L_Edd) and shape of the hard X-ray spectral energy distribution (namely the photon index Gamma_X) for a large sample of 228 hard X-ray selected, low-redshift active galactic nuclei (AGN), drawn from the Swift/BAT AGN Spectroscopic Survey (BASS). This includes 30 AGN for which black hole mass (and therefore L/L_Edd) is measured directly through masers, spatially resolved gas or stellar dynamics, or reverberation mapping. The high quality and broad energy coverage of the data provided through BASS allow us to examine several alternative determinations of both Gamma_X and L/L_Edd. For the BASS sample as a whole, we find a statistically significant, albeit very weak correlation between Gamma_X and L/L_Edd. The best-fitting relations we find, Gamma_X=0.15 log(L/L_Edd)+const., are considerably shallower than those reported in previous studies. Moreover, we find no corresponding correlations among the subsets of AGN with different M_BH determination methodology. In particular, we find no robust evidence for a correlation when considering only those AGN with direct or single-epoch M_BH estimates. This latter finding is in contrast to several previous studies which focused on z>0.5 broad-line AGN. We discuss this tension and conclude that it can be partially accounted for if one adopts a simplified, power-law X-ray spectral model, combined with L/L_Edd estimates that are based on the continuum emission and on single-epoch broad line spectroscopy in the optical regime. We finally highlight the limitations on using Gamma_X as a probe of supermassive black hole evolution in deep extragalactic X-ray surveys.
Tidal streams of disrupting dwarf galaxies orbiting around their host galaxy offer a unique way to constrain the shape of galactic gravitational potentials. Such streams can be used as leaning tower gravitational experiments on galactic scales. The most well motivated modification of gravity proposed as an alternative to dark matter on galactic scales is Milgromian dynamics (MOND), and we present here the first ever N-body simulations of the dynamical evolution of the disrupting Sagittarius dwarf galaxy in this framework. Using a realistic baryonic mass model for the Milky Way, we attempt to reproduce the present-day spatial and kinematic structure of the Sagittarius dwarf and its immense tidal stream that wraps around the Milky Way. With very little freedom on the original structure of the progenitor, constrained by the total luminosity of the Sagittarius structure and by the observed stellar mass-size relation for isolated dwarf galaxies, we find reasonable agreement between our simulations and observations of this system. The observed stellar velocities in the leading arm can be reproduced if we include a massive hot gas corona around the Milky Way that is flattened in the direction of the principal plane of its satellites. This is the first time that tidal dissolution in MOND has been tested rigorously at these mass and acceleration scales.
We present results on the dust attenuation of galaxies at redshift $\sim3-6$ by studying the relationship between the UV spectral slope ($\beta_{\rm UV}$) and the infrared excess (IRX; $L_{\rm IR}$/$L_{\rm UV}$) using ALMA far-infrared continuum observations. Our study is based on a sample of 67 massive, star-forming galaxies with a median mass of $M_{\ast}\sim 10^{10.7}\,M_{\rm \odot}$ spanning a redshift range $z=2.6-3.7$ (median $z=3.2$) that were observed with ALMA band-6. Both the individual ALMA detections (41 sources) and stacks including all galaxies show the IRX-$\beta_{\rm UV}$ relationship at $z\sim3$ is mostly consistent with that of local starburst galaxies on average. However, we find evidence for a large dispersion around the mean relationship by up to $\pm0.5$ dex. Nevertheless, the locally calibrated dust correction factors based on the IRX-$\beta_{\rm UV}$ relation are on average applicable to main-sequence $z\sim3$ galaxies. This does not appear to be the case at even higher redshifts, however. Using public ALMA observations of $z\sim4-6$ galaxies we find evidence for a significant evolution in the IRX-$\beta_{\rm UV}$ and the IRX-$M_{\ast}$ relations beyond $z\sim3$ toward lower IRX values. We discuss several caveats that could affect these results, including the assumed dust temperature. ALMA observations of larger $z>3$ galaxy samples will be required to confirm this intriguing redshift evolution.
We report the large effort which is producing comprehensive high-level young star cluster (YSC) catalogues for a significant fraction of galaxies observed with the Legacy ExtraGalactic UV Survey (LEGUS) Hubble treasury program. We present the methodology developed to extract cluster positions, verify their genuine nature, produce multiband photometry (from NUV to NIR), and derive their physical properties via spectral energy distribution fitting analyses. We use the nearby spiral galaxy NGC628 as a test case for demonstrating the impact that LEGUS will have on our understanding of the formation and evolution of YSCs and compact stellar associations within their host galaxy. Our analysis of the cluster luminosity function from the UV to the NIR finds a steepening at the bright end and at all wavelengths suggesting a dearth of luminous clusters. The cluster mass function of NGC628 is consistent with a power-law distribution of slopes $\sim -2$ and a truncation of a few times $10^5$ M$_\odot$. After their formation YSCs and compact associations follow different evolutionary paths. YSCs survive for a longer timeframe, confirming their being potentially bound systems. Associations disappear on time scales comparable to hierarchically organized star-forming regions, suggesting that they are expanding systems. We find mass-independent cluster disruption in the inner region of NGC628, while in the outer part of the galaxy there is little or no disruption. We observe faster disruption rates for low mass ($\leq$ $10^4$ M$_\odot$) clusters suggesting that a mass-dependent component is necessary to fully describe the YSC disruption process in NGC628.
We describe a new Large Program in progress on the Gemini North and South telescopes: Gemini Observations of Galaxies in Rich Early Environments (GOGREEN). This is an imaging and deep spectroscopic survey of 21 galaxy systems at $1<z<1.5$, selected to span a factor $>10$ in halo mass. The scientific objectives include measuring the role of environment in the evolution of low-mass galaxies, and measuring the dynamics and stellar contents of their host haloes. The targets are selected from the SpARCS, SPT, COSMOS and SXDS surveys, to be the evolutionary counterparts of today's clusters and groups. The new red-sensitive Hamamatsu detectors on GMOS, coupled with the nod-and-shuffle sky subtraction, allow simultaneous wavelength coverage over $\lambda\sim 0.6$--$1.05\mu$m, and this enables a homogeneous and statistically complete redshift survey of galaxies of all types. The spectroscopic sample targets galaxies with AB magnitudes $z^{\prime}<24.25$ and [3.6]$\mu$m$<22.5$, and is therefore statistically complete for stellar masses $M_\ast\gtrsim10^{10.3}M_\odot$, for all galaxy types and over the entire redshift range. Deep, multiwavelength imaging has been acquired over larger fields for most systems, spanning $u$ through $K$, in addition to deep IRAC imaging at 3.6$\mu$m. The spectroscopy is $\sim 50$ per cent complete as of semester 17A, and we anticipate a final sample of $\sim 500$ new cluster members. Combined with existing spectroscopy on the brighter galaxies from GCLASS, SPT and other sources, GOGREEN will be a large legacy cluster and field galaxy sample at this redshift that spectroscopically covers a wide range in stellar mass, halo mass, and clustercentric radius.
Present-day clusters are massive halos containing mostly quiescent galaxies, while distant protoclusters are extended structures containing numerous star-forming galaxies. We investigate the implications of this fundamental change in a cosmological context using a set of N-body simulations and semi-analytic models. We find that the fraction of the cosmic volume occupied by all (proto)clusters increases by nearly three orders of magnitude from z=0 to z=7. We show that (proto)cluster galaxies are an important, and even dominant population at high redshift, as their expected contribution to the cosmic star-formation rate density rises (from 1% at z=0) to 20% at z=2 and 40% at z=7. Protoclusters thus provide a significant fraction of the cosmic ionizing photons, and may have been crucial in driving the timing and topology of cosmic reionization. Internally, the average history of cluster formation can be described by three distinct phases: at z~10-5, galaxy growth in protoclusters proceeded in an inside-out manner, with centrally dominant, compact cores that are excellent targets for JWST; at z~5-1.5, rapid star formation occurred within the entire 10-20 Mpc structures, forming half of their total stellar mass by z~2 (2 Gyr earlier than that of the Universe as a whole); at z~1.5, violent gravitational collapse drove these stellar contents into single cluster halos, largely erasing the details of cluster galaxy formation due to relaxation and virialization. Our results motivate observations of distant protoclusters in order to understand the rapid, extended stellar growth during Cosmic Noon, and their connection to reionization during Cosmic Dawn.
GalRotpy is an educational python-based visual tool, which is useful to understand how is the contribution of each mass components to the gravitational potential of disk-like galaxies seen by its rotation curve. The standard gravitational potential of disk-like galaxies can be built by the contribution of a Miyamoto-Nagai potential model for both a Bulge/Core and Thin Disk together with a NFW (Navarro-Frenk-White) potential model for the Dark matter Halo. The GalRotpy tool can give the first approximation of the galaxy rotation curve, using models for the stellar component: (i) a rotation curve composed by both Thin a Thick Disks and (ii) exponential disk model (iii) Bulge and (iv) Dark Halo. GalRotpy routines are based on galpy package. This add-on routine is useful to modeling rotation curves of the Milky Way galaxy and other extragalactic sources. In this work we present two study cases, the rotation curve of NGC6361 CALIFA galaxy and M33 galaxy based in the data of L\'opez et al. 2017. This development began in Galactic Dynamics lecture held at Observatorio Astron\'omico Nacional.
In order to understand the origin of observed molecular cloud properties, it is critical to understand how clouds interact with their environments during their formation, growth, and collapse. It has been suggested that accretion-driven turbulence can maintain clouds in a highly turbulent state, preventing runaway collapse, and explaining the observed non-thermal velocity dispersions. We present 3D, AMR, MHD simulations of a kiloparsec-scale, stratified, supernova-driven, self-gravitating, interstellar medium, including diffuse heating and radiative cooling. These simulations model the formation and evolution of a molecular cloud population in the turbulent interstellar medium. We use zoom-in techniques to focus on the dynamics of the mass accretion and its history for individual molecular clouds. We find that mass accretion onto molecular clouds proceeds as a combination of turbulent and near free-fall accretion of a gravitationally bound envelope. Nearby supernova explosions have a dual role, compressing the envelope, boosting accreted mass, but also disrupting parts of the envelope and eroding mass from the cloud's surface. It appears that the inflow rate of kinetic energy onto clouds from supernova explosions is insufficient to explain the net rate of charge of the cloud kinetic energy. In the absence of self-consistent star formation, conversion of gravitational potential into kinetic energy during contraction seems to be the main driver of non-thermal motions within clouds. We conclude that although clouds interact strongly with their environments, bound clouds are always in a state of gravitational contraction, close to runaway, and their properties are a natural result of this collapse.
We report on new JVLA observations performed at 3 GHz and 5.5 GHz of Abell 2626. The cluster has been the object of several studies in the recent years due to its peculiar radio emission, which shows a complex system of symmetric radio arcs characterized by a steep spectrum. The origin of these radio sources is still unclear. Due to their mirror symmetry toward the center, it has been proposed that they may be created by pairs of precessing jets powered by the inner AGN. The new JVLA observations were requested with the specific aim of detecting extended emission on frequencies higher than 1.4 GHz, in order to constrain the jet-precession model by analyzing the spectral index and radiative age patterns alongs the arcs. We performed a standard data reduction of the JVLA datasets with the software CASA. By combining the new 3 GHz data with the archival 1.4 GHz VLA dataset we produced a spectral index maps of the extended emission, and then we estimated the radiative age of the arcs by assuming that the plasma was accelerated in moving hot-spots tracing the arcs. Thanks to the high sensitivity of the JVLA, we achieve the detection of the arcs at 3 GHz and extended emission at 5.5 GHz. We measure a mean spectral index <-2.5 for the arcs up to 3 GHz. No clear spectral index, or radiative age, trend is detected across the arcs which may challenge the interpretation based on precession or put strong constraints on the jet-precession period. In particular, by analyzing the radiative age distribution along the arcs, we were able to provide for the first time a time-scale < 26 Myr of the jet-precession period.
The shock L1157-B1 driven by the low-mass protostar L1157-mm is an unique environment to investigate the chemical enrichment due to molecules released from dust grains. IRAM-30m and Plateau de Bure Interferometer observations allow a census of Si-bearing molecules in L1157-B1. We detect SiO and its isotopologues and, for the first time in a shock, SiS. The strong gradient of the [SiO/SiS] abundance ratio across the shock (from >=180 to ~25) points to a different chemical origin of the two species. SiO peaks where the jet impacts the cavity walls ([SiO/H2] ~ 1e-6), indicating that SiO is directly released from grains or rapidly formed from released Si in the strong shock occurring at this location. In contrast, SiS is only detected at the head of the cavity opened by previous ejection events ([SiS/H2] ~ 2e-8). This suggests that SiS is not directly released from the grain cores but instead should be formed through slow gas-phase processes using part of the released silicon. This finding shows that Si-bearing molecules can be useful to distinguish regions where grains or gas-phase chemistry dominates.
The Purple Mountain Observatory 13.7 m radio telescope has been used to search for 95 GHz (8$_0$--7$_1$A$^+$) class I methanol masers towards 1020 Bolocam Galactic Plane Survey (BGPS) sources, leading to 213 detections. We have compared the line width of the methanol and HCO$^+$ thermal emission in all of the methanol detections and on that basis we find 205 of the 213 detections are very likely to be masers. This corresponds to an overall detection rate of 95 GHz methanol masers towards our BGPS sample of 20%. Of the 205 detected masers 144 (70%) are new discoveries. Combining our results with those of previous 95 GHz methanol masers searches, a total of four hundred and eighty-one 95 GHz methanol masers are now known, we have compiled a catalog listing the locations and properties of all known 95 GHz methanol masers.
We present a new technique to probe the central dark matter (DM) density profile of galaxies that harnesses both the survival and observed properties of star clusters. As a first application, we apply our method to the `ultra-faint' dwarf Eridanus II (Eri II) that has a lone star cluster ~45 pc from its centre. Using a grid of collisional $N$-body simulations, incorporating the effects of stellar evolution, external tides and dynamical friction, we show that a DM core for Eri II naturally reproduces the size, radial light profile and projected position of its star cluster. By contrast, a dense cusped galaxy requires the cluster to lie implausibly far from the centre of Eri II (>1 kpc), with a high inclination orbit that must be observed at a particular orbital phase. Our results imply that either a cold DM cusp was `heated up' at the centre of Eri II by bursty star formation, or we are seeing the first unambiguous evidence for physics beyond cold DM.
We use integrated-light spectroscopic observations to measure metallicities and chemical abundances for two extragalactic young massive star clusters (NGC1313-379 and NGC1705-1). The spectra were obtained with the X-Shooter spectrograph on the ESO Very Large Telescope. We compute synthetic integrated-light spectra, based on colour-magnitude diagrams for the brightest stars in the clusters from Hubble Space Telescope photometry and theoretical isochrones. Furthermore, we test the uncertainties arising from the use of Colour Magnitude Diagram (CMD) +Isochrone method compared to an Isochrone-Only method. The abundances of the model spectra are iteratively adjusted until the best fit to the observations is obtained. In this work we mainly focus on the optical part of the spectra. We find metallicities of [Fe/H] = $-$0.84 $\pm$ 0.07 and [Fe/H] = $-$0.78 $\pm$ 0.10 for NGC1313-379 and NGC1705-1, respectively. We measure [$\alpha$/Fe]=$+$0.06 $\pm$ 0.11 for NGC1313-379 and a super-solar [$\alpha$/Fe]=$+$0.32 $\pm$ 0.12 for NGC1705-1. The roughly solar [$\alpha$/Fe] ratio in NGC1313-379 resembles those for young stellar populations in the Milky Way (MW) and the Magellanic Clouds, whereas the enhanced [$\alpha$/Fe] ratio in NGC1705-1 is similar to that found for the cluster NGC1569-B by previous studies. Such super-solar [$\alpha$/Fe] ratios are also predicted by chemical evolution models that incorporate the bursty star formation histories of these dwarf galaxies. Furthermore, our $\alpha$-element abundances agree with abundance measurements from H II regions in both galaxies. In general we derive Fe-peak abundances similar to those observed in the MW and Large Magellanic Cloud (LMC) for both young massive clusters. For these elements, however, we recommend higher-resolution observations to improve the Fe-peak abundance measurements.
We explore the structure and statistics of multiphase, magnetized ISM turbulence in the local Milky Way by means of driven periodic box numerical MHD simulations. Using the higher order-accurate PPML solver, we carry out a small parameter survey varying the mean magnetic field strength and density while fixing the rms velocity to observed values. We quantify numerous characteristics of the transient and steady-state turbulence, including its thermodynamics and phase structure, kinetic and magnetic energy power spectra, structure functions, and distribution functions of density, column density, pressure, and magnetic field strength. The simulations reproduce many observables of the local ISM, including molecular clouds, such as the ratio of turbulent to mean magnetic field at 100 pc scale, the mass and volume fractions of thermally stable HI, the lognormal distribution of column densities, the mass-weighted distribution of thermal pressure, and the linewidth-size relationship for molecular clouds. Our models predict the shape of magnetic field PDFs, which are strongly non-Gaussian, and the relative alignment of magnetic field and density structures. Finally, our models show how the observed low rates of star formation per free-fall time are controlled by the multiphase thermodynamics and large-scale turbulence.
Any measurement made using an N-body simulation is subject to noise due to the finite number of particles used to sample the dark matter distribution function, and the lack of structure below the simulation resolution. This noise can be particularly significant when attempting to measure intrinsically small quantities, such as halo spin. In this work we develop a model to describe the effects of particle noise on halo spin parameters. This model is calibrated using N-body simulations in which the particle noise can be treated as a Poisson process on the underlying dark matter distribution function, and we demonstrate that this calibrated model reproduces measurements of halo spin parameter error distributions previously measured in N-body convergence studies. Utilizing this model, along with previous measurements of the distribution of halo spin parameters in N-body simulations, we place constraints on the noise-free distribution of halo spins. We find that the noise-free median spin is 3% lower than that measured directly from the N-body simulation, corresponding to a shift of approximately 40 times the statistical uncertainty in this measurement arising purely from halo counting statistics. We also show that measurement of the spin of an individual halo to 10% precision requires at least $4\times 10^4$ particles in the halo - for halos containing 200 particles the fractional error on spins measured for individual halos is of order unity.
Context. We are creating the AKARI mid-infrared all-sky diffuse maps. Through a foreground removal of the zodiacal emission, we serendipitously detected a bright residual component whose angular size is about 50 x 20 deg. at a wavelength of 9 micron. Aims. We investigate the origin and the physical properties of the residual component. Methods. We measured the surface brightness of the residual component in the AKARI mid-infrared all-sky maps. Results. The residual component was significantly detected only in 2007 January, even though the same region was observed in 2006 July and 2007 July, which shows that it is not due to the Galactic emission. We suggest that this may be a small cloud passing near the Earth. By comparing the observed intensity ratio of I_9um/I_18um with the expected intensity ratio assuming thermal equilibrium of dust grains at 1 AU for various dust compositions and sizes, we find that dust grains in the moving cloud are likely to be much smaller than typical grains that produce the bulk of the zodiacal light. Conclusions. Considering the observed date and position, it is likely that it originates in the solar coronal mass ejection (CME) which took place on 2007 January 25.
PSR B1820$-$30A is located in the globular cluster NGC 6624 and is the closest known pulsar to the centre of any globular cluster. We present more than 25 years of high-precision timing observations of this millisecond pulsar and obtain four rotational frequency time derivative measurements. Modelling these higher-order derivatives as being due to orbital motion, we find solutions which indicate that the pulsar is in either a low-eccentricity ($0.33\lesssim e\lesssim0.4$) smaller orbit with a low mass companion (such as a main sequence star, white dwarf, neutron star, or stellar mass black hole) or a high-eccentricity ($e\gtrsim0.9$) larger orbit with a massive companion. The cluster mass properties and the observed properties of 4U 1820$-$30 and the other pulsars in the cluster argue against the low-eccentricity possibility. The high-eccentricity solution reveals that the pulsar is most likely orbiting around an intermediate-mass black hole (IMBH) of mass $> 7,500$~M$_\odot$ located at the cluster centre. A gravitational model for the globular cluster, which includes such a central black hole (BH), predicts an acceleration that is commensurate with that measured for the pulsar. It further predicts that the model-dependent minimum mass of the IMBH is $\sim60,000$~M$_\odot$. Accounting for the associated contribution to the observed period derivative indicates that the $\gamma$-ray efficiency of the pulsar should be between 0.08 and 0.2. Our results suggest that other globular clusters may also contain central black holes and they may be revealed by the study of new pulsars found sufficiently close to their centres.
The recent detection of fullerene (C$_{60}$) in space and the positive assignment of five diffuse interstellar bands to C$_{60}^+$ reinforce the notion that fullerene-related compounds can be efficiently formed in circumstellar envelopes and be present in significant quantities in the interstellar medium. Experimental studies have shown that C$_{60}$ can be readily hydrogenated, raising the possibility that hydrogenated fullerenes (or fulleranes, C$_{60}$H$_m$, $m=1-60$) may be abundant in space. In this paper, we present theoretical studies of the vibrational modes of isomers of C$_{60}$H$_m$. Our results show that the four mid-infrared bands from the C$_{60}$ skeletal vibrations remain prominent in slightly hydrogenated C$_{60}$, but their strengths diminish in different degrees with increasing hydrogenation. It is therefore possible that the observed infrared bands assigned to C$_{60}$ could be due to a mixture of fullerenes and fulleranes. This provides a potential explanation for the observed scatter of the C$_{60}$ band ratios. Our calculations suggest that a feature around 15$\mu$m due to the breathing mode of heavily hydrogenated C$_{60}$ may be detectable astronomically. A preliminary search for this feature in 35 C$_{60}$ sources is reported.
The GW Ori system is a pre-main sequence triple system (GW Ori A/B/C) with companions (GW Ori B/C) at $\sim$1 AU and $\sim$8 AU, respectively, from the primary (GW Ori A). The primary of the system has a mass of 3.9 $M_{\odot}$, but shows a spectral type of G8. Thus, GW Ori A could be a precursor of a B star, but it is still at an earlier evolutionary stage than Herbig Be stars. GW Ori provides us an ideal target for experiments and observations (being a "blown-up" upscaled Solar System with a very massive "sun" and at least two "upscaled planets"). We present the first spatially-resolved millimeter interferometric observations of the disk around the triple pre-main-sequence system GW Ori, obtained with the the Submillimeter Array, both in continuum and in the $^{12}{\rm CO} J=2-1$, $^{13}{\rm CO} J=2-1$, and ${\rm C^{18}O} J=2-1$ lines. These new data reveal a huge, massive, and bright disk in the GW Ori system. The dust continuum emission suggests a disk radius around 400 AU. But, the $^{12}{\rm CO} J=2-1$ emission shows much more extended disk with a size around 1300 AU. Due to the spatial resolution ($\sim$1$"$), we cannot detect the gap in the disk which is inferred from spectral energy distribution (SED) modeling. We characterize the dust and gas properties in the disk by comparing the observations with the predictions from the disk models with various parameters calculated with a Monte Carlo radiative transfer code RADMC-3D. The disk mass is around 0.12 $M_{\odot}$, and the disk inclination with respect to the line of sight is around $\sim$ 35$^\circ$. The kinematics in the disk traced by the CO line emission strongly suggest that the circumstellar material in the disk is in Keplerian rotation around GW Ori. Tentatively substantial ${\rm C^{18}O}$ depletion in gas phase is required to explain the characteristics of the line emission from the disk.
Links to: arXiv, form interface, find, astro-ph, recent, 1705, contact, help (Access key information)