It is still on open question to what degree the cluster environment influences the sizes of protoplanetary discs surrounding young stars. Particularly so for the short-lived clusters typical for the solar neighbourhood in which the stellar density and therefore the influence of the cluster environment changes considerably over the first 10 Myr. In previous studies often the effect of the gas on the cluster dynamics has been neglected, this is remedied here. Using the code NBody6++ we study the stellar dynamics in different developmental phases - embedded, expulsion, expansion - including the gas and quantify the effect of fly-bys on the disc size. We concentrate on massive clusters ($M_{\text{cl}} \geq 10^3 - 6 \cdot 10^4 M_{\text{Sun}}$), which are representative for clusters like the Orion Nebula Cluster (ONC) or NGC 6611. We find that not only the stellar density but also the duration of the embedded phase matters. The densest clusters react fastest to the gas expulsion and drop quickly in density, here 98% of relevant encounters happen before gas expulsion. By contrast, discs in sparser clusters are initially less affected but as they expand slower 13% of discs are truncated after gas expulsion. For ONC-like clusters we find that usually discs larger than 500 AU are affected by the environment, which corresponds to the observation that 200 AU-sized discs are common. For NGC 6611-like clusters disc sizes are cut-down on average to roughly 100 AU. A testable hypothesis would be that the discs in the centre of NGC 6611 should be on average ~20 AU and therefore considerably smaller than in the ONC.
We present new results from the widest narrow band survey search for Lyman-alpha (Lya) emitters at z=5.7, just after reionization. We survey a total of 7 deg$^2$ spread over the COSMOS, UDS and SA22 fields. We find over 11,000 line emitters, out of which 514 are robust Lya candidates at z=5.7 within a volume of 6.3x10$^6$ Mpc$^3$. Our Lya emitters span a wide range in Lya luminosities, from faint to bright (L$_{\rm Ly\alpha}\sim10^{42.5-44}$ erg s$^{-1}$) and rest-frame equivalent widths (EW$_0$~25-1000 \AA) in a single, homogeneous data-set. By combining all our fields we find that the faint end slope of the z=5.7 Lya luminosity function is very steep, with $\alpha=-2.3^{+0.4}_{-0.3}$. We also present an updated z=6.6 Lya luminosity function, based on comparable volumes and obtained with the same methods, which we directly compare with that at z=5.7. We find a significant decline of the number density of faint Lya emitters from z=5.7 to z=6.6 (by $0.5\pm0.1$ dex), but no evolution at the bright end/no evolution in L*. Faint Lya emitters at z=6.6 show much more extended haloes than those at z=5.7, suggesting that neutral Hydrogen plays an important role, increasing the scattering and leading to observations missing faint Lya emission within the epoch of reionization. All together, our results indicate that we are observing patchy reionization which happens first around the brightest Lya emitters, allowing the number densities of those sources to remain unaffected by the increase of neutral Hydrogen fraction from z~5 to z~7.
The analytic "equilibrium model" for galaxy evolution using a mass balance equation is able to reproduce mean observed galaxy scaling relations between stellar mass, halo mass, star formation rate (SFR) and metallicity across the majority of cosmic time with a small number of parameters related to feedback. Here we aim to test this data-constrained model to quantify deviations from the mean relation between stellar mass and SFR, i.e. the star-forming galaxy main sequence (MS). We implement fluctuation in halo accretion rates parameterised from merger-based simulations, and quantify the intrinsic scatter introduced into the MS under the assumption that fluctuations in star formation follow baryonic inflow fluctuations. We predict the 1-sigma MS scatter to be ~ 0.2 - 0.25 dex over the stellar mass range 10^8 Mo to 10^11 Mo and a redshift range 0.5 < z < 3 for SFRs averaged over 100 Myr. The scatter increases modestly at z > 3, as well as by averaging over shorter timescales. The contribution from merger-induced star formation is generally small, around 5% today and 10 - 15% during the peak epoch of cosmic star formation. These results are generally consistent with available observations, suggesting that deviations from the MS primarily reflect stochasticity in the inflow rate owing to halo mergers.
We describe a weak lensing view of the downsizing of star forming galaxies based on cross correlating a weak lensing ($\kappa$) map with a predicted map constructed from a redshift survey. Moderately deep and high resolution images with Subaru/Hyper Suprime-Cam covering the 4 deg^2 DLS F2 field provide a $\kappa$ map with 1 arcmin resolution. A dense complete redshift survey of the F2 field including 12,705 galaxies with $R\leq20.6$ is the basis for construction of the predicted map. The zero-lag cross-correlation between the \kappa and predicted maps is significant at the $30\sigma$ level. The width of the cross-correlation peak is comparable with the angular scale of rich cluster at $z\sim0.3$, the median depth of the redshift survey. Slices of the predicted map in $\delta{z} = 0.05$ redshift bins enable exploration of the impact of structure as a function of redshift. The zero-lag normalised cross-correlation has significant local maxima at redshifts coinciding with known massive X-ray clusters. Even in slices where there are no known massive clusters, there is significant signal in the cross-correlation originating from lower mass groups that trace the large-scale of the universe. Spectroscopic $D_n4000$ measurements enable division of the sample into star-forming and quiescent populations. The significance of the cross-correlation with structure containing star-forming galaxies increases with redshift from $5\sigma$ at $z = 0.3$ to $7 \sigma$ at $z = 0.5$. The weak lensing results are consistent with the downsizing view of galaxy evolution established on the basis of many other independent studies.
In this paper we study the impact of cosmic variance and observational uncertainties in constraining the mass and occupation fraction, $f_{\rm occ}$, of dark matter halos hosting Ly-$\alpha$ Emitting Galaxies (LAEs) at high redshift. To this end, we construct mock catalogs from an N-body simulation to match the typical size of observed fields at $z=3.1$ ($\sim 1 {\rm deg^2}$). In our model a dark matter halo with mass in the range $M_{\rm min}<M_{\mathrm h}<M_{\rm max}$ can only host one detectable LAE at most. We proceed to explore the parameter space determined by $M_{\rm min}$,$M_{\rm max}$ and $f_{\rm occ}$ with a Markov Chain Monte-Carlo algorithm using the angular correlation function (ACF) and the LAEs number density as observational constraints. We find that the preferred minimum and maximum masses in our model span a wide range $10^{10.0}h^{-1}{\rm{M_{\odot}}}\leq M_{\rm min} \leq 10^{11.1}h^{-1}{\rm{M_{\odot}}}$ , $10^{11.0}h^{-1}{\rm{M_{\odot}}}\leq M_{\rm max} \leq 10^{13.0}h^{-1}{\rm{M_{\odot}}}$; followed by a wide range in the occupation fraction $0.02\leq f_{\rm occ} \leq 0.30$. As a consequence the median mass, $M_{50}$, of all the consistent models has a large uncertainty $M_{50} = 3.16^{+9.34}_{-2.37}\times 10^{10}$$h^{-1}{\rm{M_{\odot}}}$. However, we find that the same individual models have a relatively tight $1\sigma$ scatter around the median mass $\Delta M_{1\sigma} = 0.55^{+0.11}_{-0.31}$ dex. We are also able to show that \focc\ is uniquely determined by $M_{\rm min}$, regardless of $M_{\rm max}$. We argue that upcoming large surveys covering at least $25$ deg$^{2}$ should be able to put tighter constraints on $M_{\rm min}$ and $f_{\rm occ}$ through the LAE number density distribution width constructed over several fields of $\sim 1$ deg$^{2}$.
We obtained Keck/OSIRIS near-IR adaptive optics-assisted integral-field spectroscopy to probe the morphology and kinematics of the ionized gas in four velocity-offset active galactic nuclei (AGNs) from the Sloan Digital Sky Survey. These objects possess optical emission lines that are offset in velocity from systemic as measured from stellar absorption features. At a resolution of ~0.18", OSIRIS allows us to distinguish which velocity offset emission lines are produced by the motion of an AGN in a dual supermassive black hole system, and which are produced by outflows or other kinematic structures. In three galaxies, J1018+2941, J1055+1520 and J1346+5228, the spectral offset of the emission lines is caused by AGN-driven outflows. In the remaining galaxy, J1117+6140, a counterrotating nuclear disk is observed that contains the peak of Pa$\alpha$ emission 0.2" from the center of the galaxy. The most plausible explanation for the origin of this spatially and kinematically offset peak is that it is a region of enhanced Pa$\alpha$ emission located at the intersection zone between the nuclear disk and the bar of the galaxy. In all four objects, the peak of ionized gas emission is not spatially coincident with the center of the galaxy as traced by the peak of the near-IR continuum emission. The peaks of ionized gas emission are spatially offset from the galaxy centers by 0.1"-0.4" (0.1-0.7 kpc). We find that the velocity offset originates at the location of this peak of emission, and the value of the offset can be directly measured in the velocity maps. The emission-line ratios of these four velocity-offset AGNs can be reproduced only with a mixture of shocks and AGN photoionization. Shocks provide a natural explanation for the origin of the spatially and spectrally offset peaks of ionized gas emission in these galaxies.
In this letter we describe how we use stellar dynamics information to constrain the shape of the stellar IMF in a sample of 27 early-type galaxies from the CALIFA survey. We obtain dynamical and stellar mass-to-light ratios, $\Upsilon_\mathrm{dyn}$ and $\Upsilon_{\ast}$, over a homogenous aperture of 0.5~$R_{e}$. We use the constraint $\Upsilon_\mathrm{dyn} \ge \Upsilon_{\ast}$ to test two IMF shapes within the framework of the extended MILES stellar population models. We rule out a single power law IMF shape for 75% of the galaxies in our sample. Conversely, we find that a double power law IMF shape with a varying high-mass end slope is compatible (within 1$\sigma$) with 95% of the galaxies. We also show that dynamical and stellar IMF mismatch factors give consistent results for the systematic variation of the IMF in these galaxies.
We present an overview of the "KIFF" project, which provides ultra-deep Ks-band imaging of all six of the Hubble Frontier Fields clusters Abell 2744, MACS-0416, Abell S1063, Abell 370, MACS-0717 and MACS-1149. All of these fields have recently been observed with large allocations of Directors' Discretionary Time with the HST and Spitzer telescopes covering 0.4 < lambda < 1.6 microns and 3.6--4.5 microns, respectively. VLT/HAWK-I integrations of the first four fields reach 5-sigma limiting depths of Ks~26.0 (AB, point sources) and have excellent image quality (FWHM ~ 0."4). Shorter Keck/MOSFIRE integrations of the MACS-0717 (MACS-1149) field better observable in the north reach limiting depths Ks=25.5 (25.1) with seeing FWHM ~0."4 (0."5). In all cases the Ks-band mosaics cover the primary cluster and parallel HST/ACS+WFC3 fields. The total area of the Ks-band coverage is 490 arcmin^2. The Ks-band at 2.2 microns crucially fills the gap between the reddest HST filter (1.6 micron ~ H-band) and the IRAC 3.6 micron passband. While reaching the full depths of the space-based imaging is not currently feasible from the ground, the deep Ks-band images provide important constraints on both the redshifts and the stellar population properties of galaxies extending well below the characteristic stellar mass across most of the age of the universe, down to, and including, the redshifts of the targeted galaxy clusters (z < 0.5).
Radio galaxies are among the most massive galaxies in the high redshift universe and are known to often lie in protocluster environments. We have studied the fields of seven z = 2.2 radio galaxies with HAWK-I narrow-band and broad-band imaging in order to map out their environment using Halpha emitters (HAEs). The results are compared to the blank field HAE survey HiZELS. All of the radio galaxy fields are overdense in HAEs relative to a typical HiZELS field of the same area and four of the seven are richer than all except one of 65 essentially random HiZELS subfields of the same size. The star formation rates of the massive HAEs are lower than those necessary to have formed their stellar population in the preceding Gyr - indicating that these galaxies are likely to have formed the bulk of their stars at higher redshifts, and are starting to quench.
Life on Earth relies on chiral molecules, that is, species not superimposable on their mirror images. This manifests itself in the selection of a single molecular handedness, or homochirality, across the biosphere. We present the astronomical detection of a chiral molecule, propylene oxide (CH$_3$CHCH$_2$O), in absorption toward the Galactic Center. Propylene oxide is detected in the gas phase in a cold, extended molecular shell around the embedded, massive protostellar clusters in the Sagittarius B2 star-forming region. This material is representative of the earliest stage of solar system evolution in which a chiral molecule has been found.
We compute the isotropic gravitational wave (GW) background produced by binary supermassive black holes (SBHs) in galactic nuclei. In our model, massive binaries evolve at early times via gravitational-slingshot interaction with nearby stars, and at later times by the emission of GWs. Our expressions for the rate of binary hardening in the "stellar" regime are taken from the recent work of Vasiliev et al., who show that in the non-axisymmetric galaxies expected to form via mergers, stars are supplied to the center at high enough rates to ensure binary coalescence on Gyr timescales. We also include, for the first time, the extra degrees of freedom associated with evolution of the binary's orbital plane; in rotating nuclei, interaction with stars causes the orientation and the eccentricity of a massive binary to change in tandem, leading in some cases to very high eccentricities (e>0.9) before the binary enters the GW-dominated regime. We argue that previous studies have over-estimated the mean ratio of SBH mass to galaxy bulge mass by factors of 2 - 3. In the frequency regime currently accessible to pulsar timing arrays (PTAs), our assumptions imply a factor 2 - 3 reduction in the characteristic strain compared with the values computed in most recent studies, removing the tension that currently exists between model predictions and the non-detection of GWs.
Massive galaxies at higher redshifts ($\emph{z}$ $>$ 2) show different characteristics from their local counterparts: They are compact and most likely have a disk. In this study, we trace the evolution of local massive galaxies by performing a detailed morphological analysis, namely, fitting single S\'{e}rsic profiles and performing bulge+disk decompositions. We analyze $\sim$ 250 massive galaxies selected from all CANDELS fields (COSMOS, UDS, EGS, GOODS-South and GOODS-North). We confirm that both star-forming and quiescent galaxies increase their sizes significantly from $\emph{z}$ $\approx$ 2.5 to the present day. The global S\'{e}rsic index of quiescent galaxies increases over time (from $n$ $\approx$ 2.5 to $n$ $>$ 4), while that of star-forming galaxies remains roughly constant ($n$ $\approx$ 2.5). By decomposing galaxy profiles into bulge+disk components, we find that massive galaxies at high redshift have prominent stellar disks, which are also evident from visual inspection of the images. By $z$ $\approx$ 0.5, the majority of the disks disappear and massive quiescent galaxies begin to resemble the local elliptical galaxies. Star-forming galaxies have lower bulge-to-total ratios ($B/T$) than their quiescent counterparts at each redshift bin. The bulges of star-forming and quiescent galaxies follow different evolutionary histories, while their disks evolve similarly. We conclude that major mergers, along with minor mergers, have played a crucial role in the significant size increase of high-\emph{z} galaxies and the destruction of their massive and large-scale disks.
We have developed an iterative procedure to systematically combine the millimeter and submillimeter images of OB cluster-forming molecular clouds, which were taken by ground based (CSO, JCMT, APEX, IRAM-30m) and space telescopes (Herschel, Planck). For the seven luminous ($L$$>$10$^{6}$ $L_{\odot}$) Galactic OB cluster-forming molecular clouds selected for our analyses, namely W49A, W43-Main, W43-South, W33, G10.6-0.4, G10.2-0.3, G10.3-0.1, we have performed single-component, modified black-body fits to each pixel of the combined (sub)millimeter images, and the Herschel PACS and SPIRE images at shorter wavelengths. The $\sim$10$"$ resolution dust column density and temperature maps of these sources revealed dramatically different morphologies, indicating very different modes of OB cluster-formation, or parent molecular cloud structures in different evolutionary stages. The molecular clouds W49A, W33, and G10.6-0.4 show centrally concentrated massive molecular clumps that are connected with approximately radially orientated molecular gas filaments. The W43-Main and W43-South molecular cloud complexes, which are located at the intersection of the Galactic near 3-kpc (or Scutum) arm and the Galactic bar, show a widely scattered distribution of dense molecular clumps/cores over the observed $\sim$10 pc spatial scale. The relatively evolved sources G10.2-0.3 and G10.3-0.1 appear to be affected by stellar feedback, and show a complicated cloud morphology embedded with abundant dense molecular clumps/cores. We find that with the high angular resolution we achieved, our visual classification of cloud morphology can be linked to the systematically derived statistical quantities (i.e., the enclosed mass profile, the column density probability distribution function, the two-point correlation function of column density, and the probability distribution function of clump/core separations).
The results of long XMM-Newton X-ray observations of the NW radio relic of Abell 3667 are presented. A shock is detected at the sharp outer edge of the radio relic, both in the X-ray surface brightness and the temperature profiles. The Mach number is M = 2.54^+0.80_-0.43. The temperature jump at the shock is larger than expected from the density jump, which may indicate that a dynamically important magnetic field aligned primarily parallel to the shock front is present. The gas temperature rises gradually over several arc minutes within the shock region. This could indicate that the shock energy is initially dissipated into some mix of thermal and nonthermal (e.g., turbulence) components, and that the nonthermal energy decays into heat in the post-shock region. The observed radio relic can be powered if ~0.2% of the energy dissipated in the shock goes into the (re)acceleration of relativistic electrons. We show that the observed steepening of the radio spectrum with distance behind the shock is consistent with radiative losses by the radio-emitting electrons. However, the radio spectrum immediately behind the shock is flatter than expected for linear diffusive shock acceleration of thermal electrons. This suggests that the shock re-accelerates a pre-existing population of relativistic electrons. We also detect a bright, cool region (the "Mushroom") to the south of the radio relic, which we propose is the remnant cool core of a merging subcluster, and that this subcluster was the driver for the observed NW shock. In this model, the properties of Abell 3667 are mainly the result of an offset binary merger, and the cluster is being observed about 1 Gyr after first core passage. We predict that deeper X-ray or SZ observations of the SE radio relic will reveal a second merger shock at the outer edge.
We present a software package, IDA, which can easily handle two-dimensional spectroscopy data. IDA has been written in IDL and offers a window-based interface. The available tools can visualize a recovered image from spectra at any desired wavelength interval, obtain velocity fields, velocity dispersion distributions, etc.
We have analyzed the atomic and molecular gas using the 21 cm HI and 2.6/1.3 mm CO emissions toward the TeV $\gamma$-ray supernova remnant (SNR) RCW 86 in order to identify the interstellar medium interacting with the shock waves of the SNR. We have found intensity depression in HI at a velocity range of $-46$-$-28$ km s$^{-1}$ toward the SNR, suggesting a cavity of the interstellar medium. The HI cavity corresponds with the X-ray shell consisting of thermal and non-thermal emission. The thermal X-rays are coincident with the edge of the HI distribution which indicates strong density gradient, while the non-thermal X-rays are found toward the less dense inner part of the HI cavity. The most significant non-thermal X-rays are seen toward the southwestern part of the shell where HI gas includes the dense and cold component. We also identified CO clouds which are likely interacting with the SNR shock waves in the same velocity range with HI whereas the CO clouds are distributed only in a limited part of the SNR shell. The most massive cloud is located in the southeastern part of the shell, showing detailed correspondence with the thermal X-rays. These CO clouds show an enhanced CO $J$ = 2-1/1-0 intensity ratio, suggesting heating/compression by the shock front. We interpret the present results that the shock-cloud interaction enhances non-thermal X-rays in the southwest and the thermal X-rays are emitted from the shock-heated gas of density 10-100 cm$^{-3}$. It seems likely that the progenitor of RCW 86 had low-velocity stellar winds with the white dwarf(s), because the thermal X-rays indicate the remaining medium density HI gas which was not fully swept up by the progenitor winds.
The $\Lambda$CDM model, or concordance cosmology, as it is often called, is a
paradigm at its maturity. It has been checked against a large quantity of
observations, and it passed almost all tests.
The paradigm is clearly able to describe the universe at large scale, even if
some issues remain open, like the cosmological constant problem, or the
unexplained anomalies in the CMB. However, $\Lambda$CDM clearly shows
difficulty at small scales, that could be related to our scant understanding,
from the nature of dark matter to that of gravity, or to the role of baryon
physics, which is not well understood and implemented in simulation codes or in
semi-analytic models. At this stage, it is of fundamental importance to
understand if the problems encountered by the $\Lambda$DCM model are a sign of
its limits or a sign of our failures in getting the finer details right. In the
present paper, we will review the small scale problems of the $\Lambda$CDM
model, we will discuss the proposed solutions and to what extent they are able
to give us a theory accurately describing the phenomena in the complete range
of scale of the observed universe.
Links to: arXiv, form interface, find, astro-ph, recent, 1606, contact, help (Access key information)
We present new molecular gas maps of NGC5195 (alternatively known as M51b) from the Combined Array for Research in Millimeter Astronomy (CARMA), including 12CO(1-0), 13CO(1-0), CN(1-0), and CS(2-1). NGC5195 has also been detected in 3mm continuum. NGC5195 has a 12CO/13CO ratio consistent with normal star-forming galaxies. The CN(1-0) intensity is higher than is seen in an average star-forming galaxy, possibly enhanced in the diffuse photo-dissociation regions. Stellar template fitting of the nuclear spectrum of NGC5195 shows two stellar populations: an 80% mass fraction of old (>10Gyr) and a 20% mass fraction of intermediate-aged (~1Gyr) stellar populations, providing a constraint on the timescale over which NGC5195 experienced enhanced star formation during its interaction with M51a. The average molecular gas depletion timescale in NGC5195 is: tdep=3.08Gyr, a factor of ~2 larger than the depletion timescales in nearby star-forming galaxies, but consistent with the depletion seen in CO-detected early-type galaxies. While the radio continuum emission at centimeter and millimeter wavelengths is present in the vicinity of the nucleus of NGC5195, we find it is most likely associated with nuclear star formation rather than a radio-loud AGN. Thus, despite having a substantial interaction with M51a ~1/2Gyr ago, the molecular gas in NGC5195 has resettled and is forming stars at an efficiency consistent with settled early-type galaxies at the present time.
Recent studies of galaxies in the local Universe, including those in the Local Group, find that the efficiency of environmental (or satellite) quenching increases dramatically at satellite stellar masses below ~ $10^8\ {\rm M}_{\odot}$. This suggests a physical scale where quenching transitions from a slow "starvation" mode to a rapid "stripping" mode at low masses. We investigate the plausibility of this scenario using observed HI surface density profiles for a sample of 66 nearby galaxies as inputs to analytic calculations of ram-pressure and viscous stripping. Across a broad range of host properties, we find that stripping becomes increasingly effective at $M_{*} < 10^{8-9}\ {\rm M}_{\odot}$, reproducing the critical mass scale observed. However, for canonical values of the circumgalactic medium density ($n_{\rm halo} < 10^{-3.5}$ ${\rm cm}^{-3}$), we find that stripping is not fully effective; infalling satellites are, on average, stripped of < 40 - 70% of their cold gas reservoir, which is insufficient to match observations. By including a host halo gas distribution that is clumpy and therefore contains regions of higher density, we are able to reproduce the observed HI gas fractions (and thus the high quenched fraction and short quenching timescale) of Local Group satellites, suggesting that a host halo with clumpy gas may be crucial for quenching low-mass systems in Local Group-like (and more massive) host halos.
Galaxy populations at different cosmic epochs are often linked together by comoving cumulative number density in observational studies. Many theoretical works, however, have shown that the number densities of tracked galaxy populations evolve in bulk and spread out over time. We present a number density method for linking progenitor and descendant galaxy populations which takes both of these effects into account. We define probability distribution functions that capture the evolution and dispersion of galaxy populations in comoving number density space, and use these functions to assign galaxies at one redshift $z_f$ probabilities of being progenitors or descendants of a galaxy population at another redshift $z_0$. These probabilities are then used as weights for calculating distributions of physical properties such as stellar mass, star formation rate, or velocity dispersion within the progenitor/descendant population. We demonstrate that this probabilistic method provides more accurate predictions for the evolution of physical properties then either the assumption of a constant number density or the assumption of an evolving number density in a bin of fixed width by comparing the predictions against galaxy populations directly tracked through a cosmological simulation. We find that the constant number density method performs most poorly at recovering galaxy properties, the evolving number method density slightly better, and the probabilistic number density method best of all. The improvement is present for predictions of both stellar mass as well as inferred quantities such as star formation rate and velocity dispersion which were not included in the number density fits. We demonstrate that this method can also be applied robustly and easily to observational data, and provide a code package for doing so.
We constrain the newly-introduced decorated Halo Occupation Distribution (HOD) model using SDSS DR7 measurements of projected galaxy clustering or r-band luminosity threshold samples. The decorated HOD is a model for the galaxy-halo connection that augments the HOD by allowing for the possibility of galaxy assembly bias: galaxy luminosity may be correlated with dark matter halo properties besides mass, Mvir. We demonstrate that it is not possible to rule out galaxy assembly bias using DR7 measurements of galaxy clustering alone. Moreover, galaxy samples with Mr < -20 and Mr < -20.5 favor strong central galaxy assembly bias. These samples prefer scenarios in which high-concentration are more likely to host a central galaxy relative to low-concentration halos of the same mass. We exclude zero assembly bias with high significance for these samples. Satellite galaxy assembly bias is significant for the faintest sample, Mr < -19. We find no evidence for assembly bias in the Mr < -21 sample. Assembly bias should be accounted for in galaxy clustering analyses or attempts to exploit galaxy clustering to constrain cosmology. In addition to presenting the first constraints on HOD models that accommodate assembly bias, our analysis includes several improvements over previous analyses of these data. Therefore, our inferences supersede previously-published results even in the case of a standard HOD analysis.
Hyper-luminous infrared galaxies (HyLIRGs) lie at the extreme luminosity end of the IR galaxy population with L_IR>10^13L_sun. They are thought to be closer counterparts of the more distant sub-mm galaxies, and should therefore be optimal targets to study the most massive systems in formation. We present deep Chandra observations of IRAS F15307+3252 (100 ks), a classical HyLIRG located at z=0.93, hosting a radio-loud AGN (L_1.4GHz=3.5*10^25 W/Hz). The Chandra X-ray images reveal extended, asymmetric X-ray emission in the soft 0.3-2.0 keV band, extending to 160 kpc in the southern direction. VLA observations at 1.4 GHz and 8.4 GHz reveal no radio counterpart to this extended X-ray emission. The emission is therefore most likely of thermal origin originating from a hot intragroup or intracluster medium virializing in the potential. The temperature (2 keV) and bolometric X-ray luminosity (3*10^43 erg/s) of the gas follow the expected L_X-ray-T correlation for groups and clusters of galaxies. We also find that the gas has a remarkably short cooling time of 1.2 Gyrs and a central entropy of 23 keV cm^-2. In addition, the radio observations reveal that the galaxy hosts an unresolved compact steep-spectrum (CSS) source, most likely indicating the presence of a young radio source in its early stages of evolution. We also confirm that the spectrum of the nucleus is dominated by a redshifted 6.4 keV Fe K_alpha emission line, indicating the presence of a reflection-dominated Compton-thick AGN. Finally, HST images reveal an over-density of galaxies, as well as a detection of sub-structure in IRAS F15307+3252 that correlates with soft X-ray emission. This could be a snapshot view of on-going groupings expected in a growing cluster environment. IRAS F15307+3252 might therefore be a rare example of a group in the process of transforming into a cluster.
The study of molecular gas is crucial for understanding star formation, feedback, and the broader ecosystem of a galaxy as a whole. However, we have limited understanding of its physics and distribution in all but the nearest galaxies. We present a new technique for studying the composition and distribution of molecular gas in high-redshift galaxies inaccessible to existing methods. Our proposed approach is an extension of carbon monoxide intensity mapping methods, which have garnered significant experimental interest in recent years. These intensity mapping surveys target the 115 GHz $^{12}$CO (1-0) line, but also contain emission from the substantially fainter 110 GHz $^{13}$CO (1-0) transition. The method leverages the information contained in the $^{13}$CO line by cross-correlating pairs of frequency channels in an intensity mapping survey. Since $^{13}$CO is emitted from the same medium as the $^{12}$CO, but saturates at a much higher column density, this cross-correlation provides valuable information about both the gas density distribution and isotopologue ratio, inaccessible from the $^{12}$CO alone. Using a simple model of these molecular emission lines, we show that a future intensity mapping survey can constrain the abundance ratio of these two species and the fraction of emission from optically thick regions to order $\sim15\%$. These measurements cannot be made by traditional CO observations, and consequently the proposed method will provide unique insight into the physics of star formation, feedback, and galactic ecology at high redshifts.
We present the stellar mass functions (SMFs) of passive and star-forming galaxies with a limiting mass of 10$^{10.1}$ M$_{\odot}$ in four spectroscopically confirmed Spitzer Adaptation of the Red-sequence Cluster Survey (SpARCS) galaxy clusters at 1.37 $<$ z $<$ 1.63. The clusters have 113 spectroscopically confirmed members combined, with 8-45 confirmed members each. We construct $Ks$-band-selected photometric catalogs for each cluster with an average of 11 photometric bands ranging from $u$ to 8 $\mu$m. We compare our cluster galaxies to a field sample derived from a similar $Ks$-band-selected catalog in the UltraVISTA/COSMOS field. The SMFs resemble those of the field, but with signs of environmental quenching. We find that 30 $\pm$ 20\% of galaxies that would normally be forming stars in the field are quenched in the clusters. The environmental quenching efficiency shows little dependence on projected cluster-centric distance out to $\sim$ 4 Mpc, providing tentative evidence of pre-processing and/or galactic conformity in this redshift range. We also compile the available data on environmental quenching efficiencies from the literature, and find that the quenching efficiency in clusters and in groups appears to decline with increasing redshift in a manner consistent with previous results and expectations based on halo mass growth.
In this work we model galaxy-like structures as self-gravitating fluids, and analyse their properties in the Newtonian framework. For isotropic fluids, we show that this leads to a generalised Hernquist profile that admits flat rotation curves at large radial distances. For two-fluid component models, we show analytically that physicality of the solutions demand that one of the fluids is necessarily exotic, i.e has negative pressure, excepting for the case where the density profile is that of the isothermal sphere. We reconcile this result with a corresponding relativistic analysis. Our work can be applied to cases where the gravitating fluids are interpreted as dark fluids, whose microscopic constituents are dark matter particles, which may accompany or cause gravitational collapse giving birth to galaxy like structures. We elaborate on such collapse processes, which might lead to naked singularities.
We compute the expected X-ray diffuse background and radiative feedback on the intergalactic medium (IGM) from X-ray binaries prior and during the epoch of reionization. The cosmic evolution of compact binaries is followed using a population synthesis technique that treats separately neutron stars and black hole binaries in different spectral states and is calibrated to reproduce the observed X-ray properties of galaxies at z<4. Together with an updated empirical determination of the cosmic history of star formation, recent modeling of the stellar mass-metallicity relation, and a scheme for absorption by the IGM that accounts for the presence of ionized HII bubbles during the epoch of reionization, our detailed calculations provide refined predictions of the X-ray volume emissivity and filtered radiation background from "normal" galaxies at z>6. Radiative transfer effects modulate the background spectrum, which shows a characteristic peak between 1 and 2 keV. While the filtering of X-ray radiation through the IGM slightly increases the mean excess energy per photoionization, it also weakens the radiation intensity below 1 keV, lowering the mean photoionization and heating rates. Numerical integration of the rate and energy equations shows that the contribution of X-ray binaries to the ionization of the bulk IGM is negligible, with the electron fraction never exceeding 1%. Direct HeI photoionizations are the main source of IGM heating, and the temperature of the largely neutral medium in between HII cavities increases above the temperature of the cosmic microwave background (CMB) only at z<10, when the volume filling factor of HII bubbles is already >0.1. Therefore, in this scenario, it is only at relatively late epochs that the bulk of neutral intergalactic hydrogen may be observable in 21-cm emission against the CMB.
Short version: We apply a new synergetic radiative transfer method: the derived 1D density profiles are both consistent with a cut through the Herschel PACS/SPIRE and JCMT SCUBA-2 continuum maps of L1689B and with a derived local interstellar radiation field. Choosing an appropriate cut along the filament major axis, we minimize the impact of the filament emission on the modeling. For the bulk of the core (5000-20000 au) an isothermal sphere model with a temperature of around 10 K provides the best fits. We show that the power law index of the density profile, as well as the constant temperature can be derived directly from the radial surface brightness profiles. For the inner region (< 5000 au), we find a range of densities and temperatures that are consistent with the surface brightness profiles and the local interstellar radiation field. Based on our core models, we find that pixel-by-pixel single temperature spectral energy distribution fits are incapable of determining dense core properties. We conclude that, to derive physical core properties, it is important to avoid an azimuthal average of core and filament. Correspondingly, derived core masses are too high since they include some mass of the filament, and might introduce errors when determining core mass functions. The forward radiative transfer methods also avoids the loss of information owing to smearing of all maps to the coarsest spatial resolution. We find the central core region to be colder and denser than estimated in recent inverse radiative transfer modeling, possibly indicating the start of star formation in L1689B.
Recent large-scale infrared surveys have been revealing stellar populations in the inner Galaxy seen through strong interstellar extinction in the disk. In particular, classical Cepheids with their period-luminosity and period-age relations are useful tracers of Galactic structure and evolution. Interesting groups of Cepheids reported recently include four Cepheids in the Nuclear Stellar Disk (NSD), about 200 pc around the Galactic Centre, found by Matsunaga et al. and those spread across the inner part of the disk reported by Dekany and collaborators. We here report our discovery of nearly thirty classical Cepheids towards the bulge region, some of which are common with Dekany et al., and discuss the large impact of the reddening correction on distance estimates for these objects. Assuming that the four Cepheids in the NSD are located at the distance of the Galactic Centre and that the near-infrared extinction law, i.e. wavelength dependency of the interstellar extinction, is not systematically different between the NSD and other bulge lines-of-sight, most of the other Cepheids presented here are located significantly further than the Galactic Centre. This suggests a lack of Cepheids in the inner 2.5 kpc region of the Galactic disk except the NSD. Recent radio observations show a similar distribution of star-forming regions.
Due to their relation to massive stars, long-duration gamma-ray bursts (GRBs) allow pinpointing star formation in galaxies independently of redshift, dust obscuration, or galaxy mass/size, thus providing a unique tool to investigate the star-formation history over cosmic time. About half of the optical afterglows of long-duration GRBs are missed due to dust extinction, and are primarily located in the most massive GRB hosts. In order to understand this bias it is important to investigate the amount of obscured star-formation in these GRB host galaxies. Radio emission of galaxies correlates with star-formation, but does not suffer extinction as do the optical star-formation estimators. We selected 11 GRB host galaxies with either large stellar mass or large UV-/optical-based star-formation rates (SFRs) and obtained radio observations of these with the Australia Telescope Compact Array and the Karl Jansky Very Large Array. Despite intentionally selecting GRB hosts with expected high SFRs, we do not find any star-formation-related radio emission in any of our targets. Our upper limit for GRB 100621A implies that the earlier reported radio detection was due to afterglow emission. We do detect radio emission from the position of GRB 020819B, but argue that it is in large parts, if not all, due to afterglow contamination. Half of our sample has radio-derived SFR limits which are only a factor 2--3 above the optically measured SFRs. This supports other recent studies that the majority of star formation in GRB hosts is not obscured by dust.
We present the clustering properties of a complete sample of 968 radio sources detected at 1.4 GHz by the VLA-COSMOS survey with radio fluxes brighter than 0.15 mJy. Ninety-two per cent have redshift determinations from the Laigle et al. (2016) catalogue. Based on their radio-luminosity, these objects have been divided into two populations of 644 AGN and 247 star-forming galaxies. We find r_0=11.7^{+1.0}_{-1.1} Mpc for the clustering length of the whole sample, while r_0=11.2^{+2.5}_{-3.3} Mpc and r_0=7.8^{+1.6}_{-2.1} Mpc (r_0=6.8^{+1.4}_{-1.8} Mpc if we restrict our analysis to z<0.9) are respectively obtained for AGN and star-forming galaxies. These values correspond to minimum masses for dark matter haloes of M_min=10^[13.6^{+0.3}_{-0.6}] M_sun for radio-selected AGN and M_min=10^[13.1^{+0.4}_{-1.6}] M_sun for radio-emitting star-forming galaxies (M_min=10^[12.7^{+0.7}_{-2.2}] M_sun for z<0.9). Comparisons with previous works imply an independence of the clustering properties of the AGN population with respect to both radio luminosity and redshift. We also investigate the relationship between dark and luminous matter in both populations. We obtain <M*>/M_halo<~10^{-2.7} for AGN, and <M*>/M_halo<~10^{-2.4} in the case of star-forming galaxies. Furthermore, if we restrict to z<~0.9 star-forming galaxies, we derive <M*>/M_halo<~10^{-2.1}, result which clearly indicates the cosmic process of stellar build-up as one moves towards the more local universe. Comparisons between the observed space density of radio-selected AGN and that of dark matter haloes shows that about one in two haloes is associated with a black hole in its radio-active phase. This suggests that the radio-active phase is a recurrent phenomenon.
We present The Machine, an artificial neural network (ANN) capable of differentiating between the numbers of Gaussian components needed to describe the emission lines of Integral Field Spectroscopic (IFS) observations. Here we show the preliminary results of the S7 first data release (Siding Spring Southern Seyfert Spectro- scopic Snapshot Survey, Dopita et al. 2015) and SAMI Galaxy Survey (Sydney-AAO Multi-object Integral Field Unit, Croom et al. 2012) to classify whether the emission lines in each spatial pixel are composed of 1, 2, or 3 different Gaussian components. Previously this classification has been done by individual people, taking an hour per galaxy. This time investment is no longer feasible with the large spectroscopic surveys coming online.
A previous study of correlations between element abundance ratios, [X/Fe], and ages of solar twin stars is extended to include Sc, Mn, Cu, and Ba. HARPS spectra with S/N > 600 are used to derive very precise (+/- 0.01 dex) differential abundances, and stellar ages with internal errors less than 1 Gyr are obtained by interpolation in the logg - Teff diagram between isochrones calculated with the Aarhus Stellar Evolution Code. For stars younger than 6 Gyr, [X/Fe] is tightly correlated with stellar age for all elements. For ages between 6 and 9 Gyr, the [X/Fe] - age correlations break down and the stars split up into two groups having respectively high and low [X/Fe] for the odd-Z elements. It is concluded that while stars in the solar neighborhood younger than about 6 Gyr were formed from interstellar gas with a smooth chemical evolution, older stars have originated from regions enriched by supernovae with different neutron excesses. Furthermore, the correlations between abundance ratios and stellar age suggest that: i) Sc is made in Type II supernovae along with the alpha-capture elements, ii) the Type II to Ia SNe yield ratio is about the same for Mn and Fe, iii) Cu is mainly made by the weak s-process in massive stars, iv) the Ba/Y yield ratio for AGB stars increases with decreasing stellar mass, v) [Y/Mg] and [Y/Al] can be used as "chemical clocks" when determining ages of solar metallicity stars.
Links to: arXiv, form interface, find, astro-ph, recent, 1606, contact, help (Access key information)
Using the L-Galaxies semi-analytic model we simultaneously fit the HI mass function, stellar mass function and galaxy colours. We find good fits to all three observations at z = 0 and to the stellar mass function and galaxy colours at z = 2. Using Markov Chain Monte Carlo (MCMC) techniques we adjust the L-Galaxies parameters to best fit the constraining data. In order to fit the HI mass function we must greatly reduce the gas surface density threshold for star formation, thus lowering the number of low HI mass galaxies. A simultaneous reduction in the star formation efficiency prevents the over production of stellar content. A simplified model in which the surface density threshold is eliminated altogether also provides a good fit to the data. Unfortunately, these changes weaken the fit to the Kennicutt-Schmidt relation and raise the star-formation rate density at recent times, suggesting that a change to the model is required to prevent accumulation of gas onto dwarf galaxies in the local universe.
Feedback from active galactic nuclei (AGN) is widely considered to be the main driver in regulating the growth of massive galaxies through heating or driving gas out of the galaxy, preventing further increase in stellar mass. Observational proof for this scenario has, however, been scarce. We have assembled a sample of 132 radio-quiet type-2 and red AGN at 0.1<z<1. We measure the kinematics of the AGN-ionized gas, the host galaxies' stellar masses and star formation rates and investigate the relationships between AGN luminosities, specific star formation rates (sSFR) and outflow strengths W90 -- the 90\% velocity width of the [OIII]5007 line power and a proxy for the AGN-driven outflow speed. Outflow strength is independent of sSFR for AGN selected on their mid-IR luminosity, in agreement with previous work demonstrating that star formation is not sufficient to produce the observed ionized gas outflows which have to be powered by AGN activity. More importantly, we find a negative correlation between W90 and sSFR in the AGN hosts with the highest SFRs, i.e., with the highest gas content, where presumably the coupling of the AGN-driven wind to the gas is strongest. This implies that AGN with strong outflow signatures are hosted in galaxies that are more `quenched' than galaxies with weaker outflow signatures. Despite the galaxies' high SFRs, we demonstrate that the outflows are not star-formation driven but indeed due to AGN-powering. This observation is consistent with the AGN having a net suppression, `negative' impact, through feedback on the galaxies' star formation history.
With the aim of improving predictions on the far infrared (FIR) line emission from Giant Molecular Clouds (GMC), we have studied the effects of photoevaporation (PE) produced by external far-ultraviolet (FUV) and ionizing (extreme-ultraviolet, EUV) radiation on the GMC structure. Our model includes: (i) an observationally-based inhomogeneous GMC density field, and (ii) its time evolution during the PE process. We find that the photoevaporation timescale, $t_{pe}$, of a typical GMC ($M_{\rm GMC}=9 \times 10^{4}\,\rm{M_{\odot}}$) is in the range $1-30$ Myr, for gas metallicity $Z=0.05-1\,\rm Z_{\odot}$, respectively. At fixed metallicity, $t_{pe}$ decreases for higher FUV fluxes, $G_0$, due to the larger temperature of the photodissociation region (PDR). The presence of EUV radiation is important at $Z<0.2\,\rm{Z_{\odot}}$, when the size of the HII layer becomes comparable to the cloud radius. We then compute the time-dependent luminosity (for which we give fitting formulae) of several key FIR lines ([CII], [OI] [OIII], CO(1-0), and CO(6-5)) during the PE process as a function of $G_0$, and $Z$ until complete photoevaporation at $t_{pe}$. Stronger FUV fluxes produce higher [CII], [OIII], and [OI] luminosities, however lasting for progressively shorter times as $t_{pe}$ decreases. Dense clumps dominate [OI] and CO(6-5) line emission as long ($\approx 3\times 10^4$ yr) as they survive photoevaporation. Noticeably, and consistently with the recent detection by Inoue et al. (2016) of a galaxy at redshift $z\approx7.2$, at low $Z$ and early times, the [OIII] line outshines [CII] emission by $\approx 100$ times. Hence, we suggest that [OIII] is a key diagnostic of low metallicity ISM, especially in galaxies characterized by very young stellar populations.
In this Letter, we present a new theoretical method for solving the chemical evolution of galaxies, by assuming the instantaneous recycling approximation for chemical elements restored by massive stars and the Delay Time Distribution formalism for the delayed chemical enrichment by Type Ia Supernovae. The galaxy gas mass assembly history, together with the assumed stellar yields and initial mass function, represent the starting point of this method. We derive a very simple and general equation which closely relates the Laplace transforms of the galaxy gas accretion and star formation history, which can be used to simplify the problem of retrieving these quantities in most of current galaxy evolution models. We find that - once the galaxy star formation history has been reconstructed from our assumptions - the differential equation for the evolution of the chemical element $X$ can be suitably solved with classical methods. We apply our model to reproduce the $[\text{O/Fe}]$ and $[\text{Si/Fe}]$ vs. $[\text{Fe/H}]$ chemical abundance patterns as observed in the Milky Way halo and disc stars, by assuming a decaying exponential infall rate of gas and different delay time distributions for Type Ia Supernovae. Although approximate, we conclude that our model with the single degenerate scenario for Type Ia Supernovae provides the best agreement with the observed set of data. Our method will be very useful in cosmological simulations and other complementary stellar population synthesis models to predict the chemical evolution of galaxies.
Magnetic fields permeate the Universe on all scales and play a key role during star formation. We study the evolution of magnetic fields around a massive metal-free (Population III) star at $z \sim 15$ during the growth of its HII region and subsequent supernova explosion by conducting three cosmological magnetohydrodynamic simulations with radiation transport. Given the theoretical uncertainty and weak observational constraints of magnetic fields in the early universe, we initialize the simulations with identical initial conditions only varying the seed field strength. We find that magnetic fields grow as $\rho^{2/3}$ during the gravitational collapse preceding star formation, as expected from ideal spherical collapse models. Massive Population III stars can expel a majority of the gas from the host halo through radiative feedback, and we find that the magnetic fields are not amplified above the spherical collapse scaling relation during this phase. However, afterwards when its supernova remnant can radiatively cool and fragment, the turbulent velocity field in and around the shell causes the magnetic field to be significantly amplified on average by $\sim$100 in the shell and up to 6 orders of magnitude behind the reverse shock. Within the shell, field strengths are on the order of a few nG at a number density of 1 cm$^{-3}$. We show that this growth is primarily caused by small-scale dynamo action initiated by Rayleigh-Taylor instabilities that form in the remnant. These strengthened fields will propagate into the first generations of galaxies, possibly affecting the nature of their star formation.
We investigate angular momentum acquisition in Milky Way sized galaxies by comparing five high resolution zoom-in simulations, each implementing identical cosmological initial conditions, but utilizing different hydrodynamic codes: Enzo, Art, Ramses, Arepo, and Gizmo-PSPH. Each code implements a distinct set of feedback and star formation prescriptions. We find that while many galaxy and halo properties vary between the different codes (and feedback prescriptions), there is qualitative agreement on the process of angular momentum acquisition in the galaxy's halo. In all simulations, cold filamentary gas accretion to the halo results in ~4 times more specific angular momentum in cold halo gas ($\lambda_{cold} \simeq 0.15$) than in the dark matter halo. At z>1, this inflow frequently results in the formation of transient cold flow disks---large co-rotating gaseous structures in the halo of the galaxy that are fueled, aligned, and kinematically connected to filamentary gas infall along the cosmic web. Due to the qualitative agreement among disparate simulations, we conclude that the buildup of high angular momentum halo gas and the transitory formation of cold flow disks are robust predictions of LCDM galaxy formation. A growing body of observational evidence suggests that this process is borne out in the real universe.
Low surface brightness (LSB) galaxies are dominated by dark matter halo from the innermost radii; hence they are ideal candidates to investigate the influence of dark matter on different dynamical aspects of spiral galaxies. Here, we study the effect of dark matter halo on grand-design, m = 2, spiral modes in a galactic disc, treated as a collisionless system, by carrying out a global modal analysis within the WKB approximation. First, we study a superthin, LSB galaxy UGC 7321 and show that it does not support discrete global spiral modes when modelled as a disc-alone system or as a disc plus dark matter system. Even a moderate increase in the stellar central surface density does not yield any global spiral modes. This naturally explains the observed lack of large-scale spiral structure in LSBs. An earlier work (Ghosh, Saini, & Jog 2016) where the galactic disc was treated as a fluid system for simplicity had shown that the dominant halo could not arrest global modes. This difference arises due to the different dispersion relation used in the two cases and which plays a crucial role in the search for global spiral modes. Thus the correct treatment of stars as a collisionless system as done here results in the non-existence of global spiral modes, in agreement with the observations. We performed a similar modal analysis for the Galaxy, and found that the dark matter halo has a negligible effect on large-scale spiral structure.
Mergers of stellar-mass black holes were recently observed in the gravitational wave window opened by LIGO. This puts the spotlight on dense stellar systems and their ability to create intermediate-mass black holes (IMBHs) through repeated merging. Unfortunately, attempts at direct and indirect IMBH detection in star clusters in the nearby universe have proven inconclusive as of now. Indirect detection methods attempt to constrain IMBHs through their effect on star cluster photometric and kinematic observables. They are usually based on looking for a specific, physically motivated signature. While this approach is justified, it may be suboptimal in its usage of the available data. Here I present a new indirect detection method, based on machine learning, that is unaffected by these restrictions. I reduce the scientific question whether a star cluster hosts an IMBH to a classification problem in the machine learning framework. I present preliminary results to illustrate how machine learning models are trained on simulated datasets and measure their performance on previously unseen, simulated data.
We present the results of CCD $UBV$ photometric and spectroscopic observations of the open cluster NGC 225. In order to determine the structural parameters of NGC 225, we calculated the stellar density profile in the cluster's field. We estimated the probabilities of the stars being physical members of the cluster using the existing astrometric data. The most likely members of the cluster were used in the determination of the astrophysical parameters of the cluster. We calculated the mean radial velocity of the cluster as $V_{r}=-8.3\pm 5.0$ km s$^{-1}$ from the optical spectra of eight stars in the cluster's field. Using the U-B vs B-V two-colour diagram and UV excesses of the F-G type main-sequence stars, the reddening and metallicity of NGC 225 were inferred as $E(B-V)=0.151\pm 0.047$ mag and $[Fe/H]=-0.11\pm 0.01$ dex, respectively. We fitted the colour-magnitude diagrams of NGC 225 with the PARSEC isochrones and derived the distance modulus, distance and age of the cluster as $\mu_{V}=9.3\pm 0.07$ mag, d=585$\pm$20 pc and $t=900\pm 100$ Myr, respectively. We also estimated the galactic orbital parameters and space velocity components of the cluster and found that the cluster has a slightly eccentric orbit of $e=0.07\pm 0.01$ and an orbital period of $P_{orb}= 255\pm 5$ Myr.
The joint use of accurate near- and mid-infrared photometry from the 2MASS and WISE catalogues has allowed the variations of the extinction law and the dust grain size distribution in high Galactic latitudes (|b|>50) at distances up to 3 kpc from the Galactic midplane to be analyzed. The modified method of extrapolation of the extinction law applied to clump giants has turned out to be efficient for separating the spatial variations of the sample composition, metallicity, reddening, and properties of the medium. The detected spatial variations of the coefficients E(H-W1)/E(H-Ks), E(H-W2)/E(H-Ks), and E(H-W3)/E(H-Ks) are similar for all high latitudes and depend only on the distance from the Galactic midplane. The ratio of short-wavelength extinction to long-wavelength one everywhere outside the Galactic disk has been found to be smaller than that in the disk and, accordingly, the mean dust grain size is larger, while the grain size distribution in the range 0.5-11 microns is shifted toward coarse dust. Specifically, the mean grain size initially increases sharply with distance from the Galactic midplane, then decreases gradually, approaching a value typical of the disk at $|Z|\approx2.4$ kpc, and, further out, stabilizes or may increase again. The coefficients under consideration change with coordinate $Z$ with a period of about 1312\pm40 pc, coinciding every 656\pm20 pc to the south and the north and showing a significant anticorrelation between their values in the southern and northern hemispheres at intermediate Z. Thus, there exists a unified large-scale periodic structure of the interstellar medium at high latitudes within at least 5 kpc. The same periodic variations have also been found for the extinction coefficient Rv within 600 pc of the Galactic midplane through the reduction of different photometric data for stars of different classes.
Infrared photometry in the J (1.2 microns), H (1.7 microns), Ks (2.2 microns) bands from the 2MASS catalogue and in the W1 (3.4 microns), W2 (4.6 microns), W3 (12 microns), W4 (22 microns) bands from the WISE catalogue is used to reveal the spatial variations of the interstellar extinction law in the infrared near the midplane of the Galaxy by the method of extrapolation of the extinction law applied to clump giants. The variations of the coefficients E(H-W1)/E(H-Ks), E(H-W2)/E(H-Ks), E(H-W3)/E(H-Ks), and E(H-W4)/E(H-Ks) along the line of sight in 2 deg per 2 deg squares of the sky centered at b=0 and l=20, 30, ..., 330, 340 deg as well as in several 4 deg per 4 deg squares with |b|=10 are considered. The results obtained here agree with those obtained by Zasowski et al. in 2009 using 2MASS and Spitzer-IRAC photometry for the same longitudes and similar photometric bands, confirming their main result: in the inner (relative to the Sun) Galactic disk, the fraction of fine dust increases with Galactocentric distance (or the mean dust grain size decreases). However, in the outer Galactic disk that was not considered by Zasowski et al., this trend is reversed: at the disk edge, the fraction of coarse dust is larger than that in the solar neighborhood. This general Galactic trend seems to be explained by the influence of the spiral pattern: its processes sort the dust by size and fragment it so that coarse and fine dust tend to accumulate, respectively, at the outer and inner (relative to the Galactic center) edges of the spiral arms. As a result, fine dust may exist only in the part of the Galactic disk far from both the Galactic center and the edge, while coarse dust dominates at the Galactic center, at the disk edge, and outside the disk.
The present study aims to characterise the physical and chemical properties of the protostellar core Orion B9-SMM3. The APEX telescope was used to perform a follow-up molecular line survey of SMM3. The following species were identified from the frequency range 218.2-222.2 GHz: $^{13}$CO, C$^{18}$O, SO, para-H$_2$CO, and E$_1$-type CH$_3$OH. The on-the-fly mapping observations at 215.1-219.1 GHz revealed that SMM3 is associated with a dense gas core as traced by DCO$^+$ and p-H$_2$CO. Altogether three different p-H$_2$CO transitions were detected with clearly broadened linewidths (8.2-11 km s$^{-1}$ in FWHM). The derived p-H$_2$CO rotational temperature, $64\pm15$ K, indicates the presence of warm gas. We also detected a narrow p-H$_2$CO line (FWHM=0.42 km s$^{-1}$) at the systemic velocity. The p-H$_2$CO abundance for the broad component appears to be enhanced by two orders of magnitude with respect to the narrow line value ($\sim3\times10^{-9}$ versus $\sim2\times10^{-11}$). The detected methanol line shows a linewidth similar to those of the broad p-H$_2$CO lines, which indicates their coexistence. The CO isotopologue data suggest that the CO depletion factor decreases from $\sim27\pm2$ towards the core centre to a value of $\sim8\pm1$ towards the core edge. In the latter position, the N$_2$D$^+$/N$_2$H$^+$ ratio is revised down to $0.14\pm0.06$. The origin of the subfragments inside the SMM3 core we found previously can be understood in terms of the Jeans instability if non-thermal motions are taken into account. The estimated fragmentation timescale, and the derived chemical abundances suggest that SMM3 is a few times $10^5$ yr old, in good agreement with its Class 0 classification inferred from the spectral energy distribution analysis. The broad p-H$_2$CO and CH$_3$OH lines, and the associated warm gas provide the first clear evidence of a molecular outflow driven by SMM3.
We attempt to find a progenitor for the ultra-faint object Segue 1 under the assumption that it formed as a dark matter free star cluster in the past. We look for orbits, using the elongation of Segue 1 on the sky as a tracer of its path. Those orbits are followed backwards in time to find the starting points of our N-body simulations. The successful orbit, with which we can reproduce Segue 1 has a proper motion of mu_alpha = -0.19 mas/yr and mu_delta = -1.9 mas/yr, placing Segue 1 near its apo-galacticon today. Our best fitting model has an initial mass of 6224 Msun and an initial scale-length of 5.75 pc.
We model the time evolution of gaps in tidal streams caused by the impact of a dark matter subhalo, while both orbit a spherical gravitational potential. To this end, we make use of the simple behaviour of orbits in action-angle space. A gap effectively results from the divergence of two nearby orbits whose initial phase-space separation is, for very cold thin streams, largely given by the impulse induced by the subhalo. We find that in a spherical potential the size of a gap increases linearly with time, while its density decreases as 1/t^2 for sufficiently long timescales. We have derived an analytic expression that shows how the growth rate depends on the mass of the perturbing subhalo, its scale and its relative velocity with respect to the stream. We have verified these scalings using N-body simulations and find excellent agreement. For example, a subhalo of mass 10^8 Msun directly impacting a very cold thin stream on an inclined orbit can induce a gap that may reach a size of several tens of kpc after a few Gyr. The gap size fluctuates importantly with phase on the orbit, and it is largest close to pericentre. This indicates that it may not be fully straightforward to invert the spectrum of gaps present in a stream to recover the mass spectrum of the subhalos.
Accurate modeling of the spectrum of thermal dust emission at millimeter wavelengths is important for improving the accuracy of foreground subtraction for CMB measurements, for improving the accuracy with which the contributions of different foreground emission components can be determined, and for improving our understanding of dust composition and dust physics. We fit four models of dust emission to high Galactic latitude COBE/FIRAS and COBE/DIRBE observations from 3 millimeters to 100 microns and compare the quality of the fits. We consider the two-level systems model because it provides a physically motivated explanation for the observed long wavelength flattening of the dust spectrum and the anticorrelation between emissivity index and dust temperature. We consider the model of Finkbeiner, Davis, and Schlegel because it has been widely used for CMB studies, and the generalized version of this model recently applied to Planck data by Meisner and Finkbeiner. For comparison we have also fit a phenomenological model consisting of the sum of two graybody components. We find that the two-graybody model gives the best fit and the FDS model gives a significantly poorer fit than the other models. The Meisner and Finkbeiner model and the two-level systems model remain viable for use in Galactic foreground subtraction, but the FIRAS data do not have sufficient signal-to-noise ratio to provide a strong test of the predicted spectrum at millimeter wavelengths.
We use a large data-set of realistic synthetic observations (PaperI) to assess how observational techniques affect the measurement of physical properties of star-forming regions. In this paper (PaperII), we explore the reliability of the measured total gas mass, dust surface density and dust temperature maps derived from modified blackbody fitting of synthetic Herschel observations. We found from our pixel-by-pixel analysis of the measured dust surface density and dust temperature a worrisome error spread especially close to star-formation sites and low-density regions, where for those "contaminated" pixels the surface densities can be under/overestimated by up to three orders of magnitude. In light of this, we recommend to treat the pixel-based results from this technique with caution in regions with active star formation. In regions of high background typical in the inner Galactic plane, we are not able to recover reliable surface density maps of individual synthetic regions, since low-mass regions are lost in the FIR background. When measuring the total gas mass of regions in moderate background, we find that modified blackbody fitting works well (absolute error:+9%;-13%) up to 10kpc distance (errors increase with distance). Commonly, the initial images are convolved to the largest common beam-size, which smears contaminated pixels over large areas. The resulting information loss makes this commonly-used technique less verifiable as now chi^2-values cannot be used as a quality indicator of a fitted pixel. Our control measurements of the total gas mass (without the step of convolution to the largest common beam size) produce similar results (absolute error:+20%;-7%) while having much lower median errors especially for the high-mass stellar feedback phase. In upcoming papers (III&IV) we test the reliability of measured star-formation rate with direct and indirect techniques.
We present radial velocities and chemical abundances for red giant branch stars in the Galactic bulge globular clusters NGC 6342 and NGC 6366. The velocities and abundances are based on measurements of high resolution (R > 20,000) spectra obtained with the MMT-Hectochelle and WIYN-Hydra spectrographs. We find that NGC 6342 has a heliocentric radial velocity of +112.5 km/s (sigma = 8.6 km/s), NGC 6366 has a heliocentric radial velocity of -122.3 km/s (sigma = 1.5 km/s), and that both clusters have nearly identical metallicities ([Fe/H] ~ -0.55). NGC 6366 shows evidence of a moderately extended O-Na anti-correlation, but more data are needed for NGC 6342 to determine if this cluster also exhibits the typical O-Na relation likely found in all other Galactic globular clusters. The two clusters are distinguished from similar metallicity field stars as having larger [Na/Fe] spreads and enhanced [La/Fe] ratios, but we find that NGC 6342 and NGC 6366 display alpha and Fe-peak element abundance patterns that are typical of other metal-rich ([Fe/H] > -1) inner Galaxy clusters. However, the median [La/Fe] abundance may vary from cluster-to-cluster.
We investigate mass segregation in group and cluster environments by
identifying galaxy analogues in high-resolution dark matter simulations.
Subhalos identified by the AHF and ROCKSTAR halo finders have similar mass
functions, independent of resolution, but different radial distributions due to
significantly different subhalo hierarchies. We propose a simple way to
classify subhalos as galaxy analogues. The radial distributions of galaxy
analogues agree well at large halo-centric radii for both AHF and ROCKSTAR but
disagree near parent halo centres where the phase-space information used by
ROCKSTAR is essential.
We see clear mass segregation at small radii (within $0.5\,r_{vir}$) with
average galaxy analogue mass decreasing with radius. Beyond the virial radius,
we find a mild trend where the average galaxy analogue mass increases with
radius. These mass segregation trends are strongest in small groups and
dominated by the segregation of low mass analogues. The lack of mass
segregation in massive galaxy analogues suggests that the observed trends are
driven by the complex accretion histories of the parent halos rather than
dynamical friction.
The James Webb Space Telescope (JWST) will be an exquisite new near-infrared
observatory with imaging and multi-object spectroscopy through ESA's NIRspec
instrument with its unique Micro-Shutter Array (MSA), allowing for slits to be
positioned on astronomical targets by opening specific 0.002"-wide micro
shutter doors.
To ensure proper target acquisition, the on-sky position of the MSA needs to
be verified before spectroscopic observations start. An onboard centroiding
program registers the position of pre-identified guide stars in a Target
Acquisition (TA) image, a short pre-spectroscopy exposure without dispersion
(image mode) through the MSA with all shutters open.
The outstanding issue is the availability of Galactic stars in the right
luminosity range for TA relative to typical high redshift targets. We explore
this here using the stars and $z\sim8$ candidate galaxies identified in the
source extractor catalogs of Brightest of Reionizing Galaxies survey
(BoRG[z8]), a pure-parallel program with Hubble Space Telescope Wide-Field
Camera 3.
We find that (a) a single WFC3 field contains enough Galactic stars to
satisfy the NIRspec astrometry requirement (20 milli-arcseconds), provided its
and the NIRspec TA's are $m_{lim}>24.5$ AB in WFC3 F125W, (b) a single WFC3
image can therefore serve as the pre-image if need be, (c) a WFC3 mosaic and
accompanying TA image satisfy the astrometry requirement at $\sim23$ AB mag in
WFC3 F125W, (d) no specific Galactic latitude requires deeper TA imaging due to
a lack of Galactic stars, and (e) a depth of $\sim24$ AB mag in WFC3 F125W is
needed if a guide star in the same MSA quadrant as a target is required.
We take the example of a BoRG identified $z\sim8$ candidate galaxy and
require a Galactic star within 20" of it. In this case, a depth of 25.5 AB in
F125W is required (with $\sim$97% confidence).
We investigate the properties of the galaxies that reionized the Universe and the history of cosmic reionization using the "Evolution and Assembly of GaLaxies and their environments" (EAGLE) cosmological hydrodynamical simulations. We obtain the evolution of the escape fraction of ionizing photons in galaxies assuming that galactic winds create channels through which 20~percent of photons escape when the local surface density of star formation is greater than $0.1$ M$_\odot$ yr$^{-1}$ kpc$^{-2}$. Such threshold behaviour for the generation of winds is observed, and the rare local objects which have such high star formation surface densities exhibit high escape fractions. In our model the luminosity-weighted mean escape fraction increases with redshift as $\bar f_{\rm esc}=0.045~((1+z)/4)^{1.1}$ at $z>3$, and the galaxy number weighted mean as $\langle f_{\rm esc} \rangle=2.2\times10^{-3}~((1+z)/4)^4$, and becomes constant $\approx0.2$ at redshift $z>10$. The escape fraction evolves as an increasingly large fraction of stars forms above the critical surface density of star formation at earlier times. This evolution of the escape fraction, combined with that of the star formation rate density from \eagle, reproduces the inferred evolution of the filling factor of ionized regions during the reionization epoch ($6<z<8$), the evolution of the post-reionization ($0\leq z<6$) hydrogen photoionization rate, and the optical depth due to Thomson scattering of the cosmic microwave background photons measured by the Planck satellite.
Links to: arXiv, form interface, find, astro-ph, recent, 1606, contact, help (Access key information)
Galaxies are thought to grow through accretion; as less massive galaxies are disrupted and merge over time, their debris results in diffuse, clumpy stellar halos enveloping the central galaxy. Here we present a study of the variation in the stellar halos of galaxies, using data from the Dragonfly Nearby Galaxies Survey (DNGS). The survey consists of wide field, deep ($\mu_{g} > 31$ mag arcsec$^{-2}$) optical imaging of nearby galaxies using the Dragonfly Telephoto Array. Our sample includes eight spiral galaxies with stellar masses similar to that of the Milky Way, inclinations of $16-90$ degrees and distances between $7-18$ Mpc. We construct stellar mass surface density profiles from the observed $g$-band surface brightness in combination with the $g-r$ color as a function of radius, and compute the halo fractions from the excess stellar mass (relative to a disk$+$bulge fit) beyond $5$ half-mass radii. We find a mean halo fraction of $0.009 \pm 0.005$ and a large RMS scatter of $1.01^{+0.9}_{-0.26}$ dex. The peak-to-peak scatter is a factor of $>100$ -- while some galaxies feature strongly structured halos resembling that of M31, three of the eight have halos that are completely undetected in our data. We conclude that spiral galaxies as a class exhibit a rich variety in stellar halo properties, implying that their assembly histories have been highly non-uniform. We find no convincing evidence for an environmental or stellar mass dependence of the halo fraction in the sample.
There is evidence that the well-established mass-metallicity relation in galaxies is correlated with a third parameter: star formation rate (SFR). The strength of this correlation may be used to disentangle the relative importance of different physical processes (e.g., infall of pristine gas, metal-enriched outflows) in governing chemical evolution. However, all three parameters are susceptible to biases that might affect the observed strength of the relation between them. We analyze possible sources of systematic error, including sample bias, application of S/N cuts on emission lines, choice of metallicity calibration, uncertainty in stellar mass determination, aperture effects, and dust. We present the first analysis of the relation between stellar mass, gas phase metallicity, and SFR using strong line abundance diagnostics from Dopita et al. (2013) for ~130,000 star-forming galaxies in the Sloan Digital Sky Survey and provide a detailed comparison of these diagnostics in an appendix. Using these abundance diagnostics yields a 30-55% weaker anti-correlation between metallicity and SFR at fixed stellar mass than that reported by Mannucci et al. (2010). We find that, for all abundance diagnostics, the anti-correlation with SFR is stronger for the relatively few galaxies whose current SFRs are elevated above their past average SFRs. This is also true for the new abundance diagnostic of Dopita et al. (2016), which gives anti-correlation between metallicity and SFR only in the high sSFR regime, in contrast to the recent results of Kashino et al. (2016). The poorly constrained strength of the relation between stellar mass, metallicity, and SFR must be carefully accounted for in theoretical studies of chemical evolution.
We present a 3 mm spectral line and continuum survey of L1451 in the Perseus Molecular Cloud. These observations are from the CARMA Large Area Star Formation Survey (CLASSy), which also imaged Barnard 1, NGC 1333, Serpens Main and Serpens South. L1451 is the survey region with the lowest level of star formation activity---it contains no confirmed protostars. HCO+, HCN, and N2H+ (J=1-0) are all detected throughout the region, with HCO+ the most spatially widespread, and molecular emission seen toward 90% of the area above N(H_2) column densities of 1.9x10^21 cm^-2. HCO+ has the broadest velocity dispersion, near 0.3 km/s on average, compared to ~0.15 km/s for the other molecules, thus representing a range from supersonic to subsonic gas motions. Our non-binary dendrogram analysis reveals that the dense gas traced by each molecule has similar hierarchical structure, and that gas surrounding the candidate first hydrostatic core (FHSC), L1451-mm, and other previously detected single-dish continuum clumps have similar hierarchical structure; this suggests that different sub-regions of L1451 are fragmenting on the pathway to forming young stars. We determined the three-dimensional morphology of the largest detectable dense gas structures to be relatively ellipsoidal compared to other CLASSy regions, which appeared more flattened at largest scales. A virial analysis shows the most centrally condensed dust structures are likely unstable against collapse. Additionally, we identify a new spherical, centrally condensed N2H+ feature that could be a new FHSC candidate. The overall results suggest L1451 is a young region starting to form its generation of stars within turbulent, hierarchical structures.
The origin of striations aligned along the local magnetic field direction in the translucent envelope of the Taurus molecular cloud is examined with new observations of 12CO and 13CO J=2-1 emission obtained with the 10~m submillimeter telescope of the Arizona Radio Observatory. These data identify a periodic pattern of excess blue and redshifted emission that is responsible for the striations. For both 12CO and 13CO, spatial variations of the J=2-1 to J=1-0 line ratio are small and are not spatially correlated with the striation locations. A medium comprised of unresolved CO emitting substructures (cells) with a beam area filling factor less than unity at any velocity is required to explain the average line ratios and brightness temperatures. We propose that the striations result from the modulation of velocities and the beam filling factor of the cells as a result of either the Kelvin-Helmholtz instability or magnetosonic waves propagating through the envelope of the Taurus molecular cloud. Both processes are likely common features in molecular clouds that are sub-Alfvenic and may explain low column density, cirrus-like features similarly aligned with the magnetic field observed throughout the interstellar medium in far-infrared surveys of dust emission.
The growth of the most massive black holes in the early universe, consistent with the detection of highly luminous quasars at $z> 6$ implies sustained, critical accretion of material to grow and power them. Given a black hole seed scenario, it is still uncertain which conditions in the early Universe allow the fastest black hole growth. Large scale hydrodynamical cosmological simulations of structure formation allow us to explore the conditions conducive to the growth of the earliest supermassive black holes. We use the cosmological hydrodynamic simulation BlueTides, which incorporates a variety of baryon physics in a (400 Mpc/h)^3 volume with 0.7 trillion particles to follow the earliest phases of black hole critical growth. At z=8 the most massive black holes (a handful) approach masses of 10^8 Msun with the most massive (with M_BH = 4 x 10^8 Msun ) being found in an extremely compact spheroid-dominated host galaxy. Examining the large-scale environment of hosts, we find that the initial tidal field is more important than overdensity in setting the conditions for early BH growth. In regions of low tidal fields gas accretes 'cold' onto the black hole and falls along thin, radial filaments straight into the center forming the most compact galaxies and most massive black holes at earliest times. Regions of high tidal fields instead induce larger coherent angular momenta and influence the formation of the first population of massive compact disks. The extreme early growth depends on the early interplay of high gas densities and the tidal field that shapes the mode of accretion. Mergers play a minor role in the formation of the first generation, rare massive BHs.
We use the Balmer decrements of the broad-line regions (BLRs) and narrow-line regions (NLRs) of active galactic nuclei (AGNs) as reddening indicators to investigate the location of the dust for four samples of AGNs with reliable estimates of the NLR contribution to the Balmer lines. Intercomparison of the NLR and BLR Balmer decrements indicates that the reddening of the NLR sets a lower limit to the reddening of the BLR. Almost no objects have high NLR reddening but low BLR reddening. The reddening of the BLR is often substantially greater than the reddening of the NLR. The BLR reddening is correlated with the equivalent widths of [O III] lines and the intensity of the [O III] lines relative to broad H\beta. We find these relationships to be consistent with the predictions of a simple model where the additional dust reddening the BLR is interior to the NLR. We thus conclude that the dust causing the additional reddening of the accretion disc and BLR is mostly located at a smaller radius than the NLR.
High-redshift radio galaxies present signs of both star formation and AGN activity, making them ideal candidates to investigate the connection and coevolution of AGN and star formation in the progenitors of present-day massive galaxies. We make use of a sample of 11 powerful radio galaxies spanning 1<z<4 which have complete coverage of their spectral energy distribution (SED) from UV to FIR wavelengths. Using Herschel data, we disentangle the relative contribution of the AGN and star formation by combining the galaxyevolutioncodePEGASE.3 with an AGN torus model. We find that three components are necessary to reproduce the observed SEDs: an evolved and massive stellar component, a submm bright young starburst, and an AGN torus. We find that powerful radio galaxies form at very high-redshift, but experience episodic and important growth at 1<z<4 as the mass of the associated starburst varies from 5 to 50% of the total mass of the system. The properties of star formation differ from source to source, indicating no general trend of the star formation properties in the most infrared luminous high-redshift radio galaxies and no correlation with the AGN bolometric luminosity. Moreover, we find that AGN scattered light have a very limited impact on broad-band SED fitting on our sample. Finally, our analysis also suggests a wide range in origins for the observed star formation,which we partially constrain for some sources.
It has long been known that environment has a large effect on star formation
in galaxies. There are several known plausible mechanisms to remove the cool
gas needed for star formation, such as strangulation, harassment and
ram-pressure stripping. It is unclear which process is dominant, and over what
range of stellar mass. In this paper, we find evidence for suppression of the
cross-correlation function between massive galaxies and less massive
star-forming galaxies, giving a measure of how less likely a galaxy is to be
star-forming in the vicinity of a more massive galaxy. We develop a formalism
for modelling environmental quenching mechanisms within the Halo Occupation
Distribution formalism.
We find that at $z \sim 2$ environment is not a significant factor in
determining quenching of star-forming galaxies, and that galaxies are quenched
with similar probabilities in group environments as they are globally. However,
by $z \sim 0.5$ galaxies are much less likely to be star forming when in a
group environment than when not. This increased probability of being quenched
does not appear to have significant radial dependence within the halo,
supportive of the quenching being caused by the halting of fresh inflows of
pristine gas, as opposed to by tidal stripping. Furthermore, by separating the
massive sample into passive and star-forming, we see that this effect is
further enhanced when the central galaxy is passive. This effect is present
only in the 1-halo term (within a halo) at high redshifts ($z>1$), but is
apparent in the 2-halo term at lower redshifts ($z<1$), a manifestation of
galactic conformity.
Searches for extra-terrestrial intelligence (SETI) using large survey data often look for possible signatures of astroengineering. We propose to search for physically impossible effects caused by highly advanced technology, by carrying out a search for disappearing galaxies and Milky Way stars. We select $\sim$ 10 million objects from USNO-B1.0 with low proper motion ($\mu$ $<$ 20 milli arcseconds / year) imaged on the sky in two epochs. We search for objects not found at the expected positions in the Sloan Digital Sky Survey (SDSS) by visually examining images of $\sim$ 290 000 USNO-B1.0 objects with no counterpart in the SDSS. We identify some spurious targets in the USNO-B1.0. We find one candidate of interest for follow-up photometry, although it is very uncertain. If the candidate eventually is found, it defines the probability of observing a disappearing-object event the last decade to less than one in one million in the given samples. Nevertheless, since the complete USNO-B1.0 dataset is 100 times larger than any of our samples, we propose an easily accessible citizen science project in search of USNO-B1.0 objects which have disappeared from the SDSS.
Observations suggest that dwarf spheroidal (dSph) galaxies exhibit large constant-density cores in the centers, which can hardly be explained by dissipationless cold dark matter simulations. Wave dark matter (${\psi {\rm DM}}$), characterized by a single parameter, the dark matter particle mass $m_{\psi}$, predicts a central soliton core in every galaxy arising from quantum pressure against gravity. Here we apply Jeans analysis to the kinematic data of eight classical dSphs so as to constrain $m_{\psi}$, and obtain $m_{\psi}=1.18_{-0.24}^{+0.28}\times10^{-22}{\,\rm eV}$ and $m_{\psi}=1.79_{-0.33}^{+0.35}\times10^{-22}{\,\rm eV}~(2\sigma)$ using the observational data sets of Walker et al. (2007) and Walker et al. (2009b), respectively. We show that the estimate of $m_{\psi}$ is sensitive to the dSphs kinematic data sets and is robust to various models of stellar density profile. We also consider multiple stellar subpopulations in dSphs and find consistent results. This mass range of $m_{\psi}$ is in good agreement with other independent estimates, such as the high-redshift luminosity functions, the reionization history, and the Thomson optical depth to the cosmic microwave background.
In most of Seyfert-1 active galactic nucei (AGN) the optical linear continuum polarization degree is usually small (less than 1%) and the polarization position angle is nearly parallel to the AGN radio-axis. However, there are many types-1 AGNs with unexplained intermediate values for both positional angles and polarization degrees. Our explanation of polarization degree and positional angle of Seyfert-1 AGNs focuses on the reflection of non-polarized radiation from sub-parsec jets in optically thick accretion discs. The presence of a magnetic field surrounding the scattering media will induce Faraday rotation of the polarization plane that may explain the intermediate values of positional angles if there is a magnetic field component normal to the accretion disc. The Faraday rotation depolarization effect in disc diminishes the competition between polarization of the reflected radiation with the parallel component of polarization and the perpendicular polarization from internal radiation of disc (the Milne problem) in favor of polarization of reflected radiation. This effect allows us to explain the observed polarization of Seyfert-1 AGN radiation even though the jet optical luminosity is much lower than the luminosity of disc. We present the calculation of polarization degrees for a number of Seyfert-1 AGNs.
Energetic winds and radiation from massive star clusters push the surrounding gas and blow superbubbles in the interstellar medium (ISM). Using 1-D hydrodynamic simulations, we study the role of radiation in the dynamics of superbubbles driven by a young star cluster of mass $10^{6}$ M$_{\odot}$. We have considered a realistic time evolution of the mechanical power as well as radiation power of the star cluster, and detailed heating and cooling processes. We find that the ratio of the radiation pressure on the shell (shocked ISM) to the thermal pressure ($\sim10^{7}$ K) of the shocked wind region is almost independent of the ambient density, and it is greater than unity before $\lesssim 1$ Myr. We explore the parameter space of density and dust opacity of the ambient medium, and find that the size of the hot gas ($\sim$ 10$^{7}$ K) cavity is insensitive to the dust opacity ($\sigma_{d}\approx(0.1-1.5)\times 10^{-21}$ cm$^{2}$), but the structure of the photoionized ($\sim10^4$ K) gas depends on it. Most of the radiative losses occur at $\sim10^{4}$ K, with sub-dominant losses at $\lesssim 10^3$ K and $\sim10^{6}-10^{8}$ K. The superbubbles can retain as high as $\sim 10\%$ of its input energy, for an ambient density of $10^{3}\,m{\rm_{H}\,cm^{-3}}$. We discuss the role of ionization parameter and recombination-averaged density in understanding the dominant feedback mechanism. Finally, we compare our results with the observations of 30 Doradus.
Photometry from the Tycho-2, 2MASS, and WISE catalogues for clump and branch giants at a distance up to 25 kpc toward the Galactic poles has allowed the variations of various characteristics of the infrared interstellar extinction law with distance to be analyzed. The results obtained by the extinction law extrapolation method are consistent for different classes of stars and different characteristics as well as with previous studies. The conventional extinction law with a low infrared extinction is characteristic of only a thin layer no farther than 100 pc from the Galactic plane and of two thin layers near $Z=-600$ and $+500$ pc. Far from the Galactic plane, in the Galactic halo, the infrared extinction law is different: the extinction in the $Ks$, $W1$, $W2$, $W3$, and $W4$ bands is, respectively, $0.17$, $0.16$, $0.16$, $0.07$, and $0.03$ of the extinction in the $V$ band. The accuracy of these coefficients is $0.03$. If the extinction law reflects primarily the grain size distribution, then the fraction of large dust grains far from the Galactic plane is greater than that in the circumsolar interstellar medium.
Through an extensive set of realistic synthetic observations (produced in Paper I), we assess in this part of the paper series (Paper III) how the choice of observational techniques affects the measurement of star-formation rates (SFRs) in star-forming regions. We test the accuracy of commonly used techniques and construct new methods to extract the SFR, so that these findings can be applied to measure the SFR in real regions throughout the Milky Way. We investigate diffuse infrared SFR tracers such as those using 24 {\mu}m, 70 {\mu}m and total infrared emission, which have been previously calibrated for global galaxy scales. We set up a toy model of a galaxy and show that the infrared emission is consistent with the intrinsic SFR using extra-galactic calibrated laws (although the consistency does not prove their reliability). For local scales, we show that these techniques produce completely unreliable results for single star-forming regions, which are governed by different characteristic timescales. We show how calibration of these techniques can be improved for single star-forming regions by adjusting the characteristic timescale and the scaling factor and give suggestions of new calibrations of the diffuse star-formation tracers. We show that star-forming regions that are dominated by high-mass stellar feedback experience a rapid drop in infrared emission once high-mass stellar feedback is turned on, which implies different characteristic timescales. Moreover, we explore the measured SFRs calculated directly from the observed young stellar population. We find that the measured point sources follow the evolutionary pace of star formation more directly than diffuse star-formation tracers.
We present the JCMT Gould Belt Survey's first look results of the southern extent of the Orion A Molecular Cloud ($\delta \leq -5\mathrm{:}31\mathrm{:}27.5$). Employing a two-step structure identification process, we construct individual catalogues for large-scale regions of significant emission labelled as islands and smaller-scale subregions called fragments using the 850 $\mu$m continuum maps obtained using SCUBA-2. We calculate object masses, sizes, column densities, and concentrations. We discuss fragmentation in terms of a Jeans instability analysis and highlight interesting structures as candidates for follow-up studies. Furthermore, we associate the detected emission with young stellar objects (YSOs) identified by Spitzer and Herschel. We find that although the population of active star-forming regions contains a wide variety of sizes and morphologies, there is a strong positive correlation between the concentration of an emission region and its calculated Jeans instability. There are, however, a number of highly unstable subregions in dense areas of the map that show no evidence of star formation. We find that only $\sim$72\% of the YSOs defined as Class 0+I and flat-spectrum protostars coincide with dense 850 $\mu$m emission structures (column densities $>3.7\times10^{21}\mathrm{\:cm}^{-2}$). The remaining 28\% of these objects, which are expected to be embedded in dust and gas, may be misclassified. Finally, we suggest that there is an evolution in the velocity dispersion of young stellar objects such that sources which are more evolved are associated with higher velocities.
We present Intercut, a Python-based program that applies secondary line identification and photometric cuts to mock galaxy surveys, in order to simulate interloper identification. This program can be used to optimize the removal of interloper contamination in upcoming surveys. Intercut reads a mock galaxy survey and an emission line sensitivity and simulates interloper removal through secondary line identification and broad-band photometry. This program is designed to use the COSMOS mock catalog, although the program can be modified for an alternative mock catalog. The output of the program returns an interloper fraction for each emission line, as well as the total fraction over all lines, as a function of redshift. We test Intercut by predicting interloper rates for the WFIRST emission line sensitivity, finding agreement with previous work. This program is publically available on Github
Is life most likely to emerge at the present cosmic time near a star like the Sun? We consider the habitability of the Universe throughout cosmic history, and conservatively restrict our attention to the context of "life as we know it" and the standard cosmological model, LCDM. The habitable cosmic epoch started shortly after the first stars formed, about 30 Myr after the Big Bang, and will end about 10 Tyr from now, when all stars will die. We review the formation history of habitable planets and find that unless habitability around low mass stars is suppressed, life is most likely to exist near 0.1 solar mass stars ten trillion years from now. Spectroscopic searches for biosignatures in the atmospheres of transiting Earth-mass planets around low mass stars will determine whether present-day life is indeed premature or typical from a cosmic perspective.
In this work, we present the new catalog of carbon stars from the LAMOST DR2 catalog. In total, 894 carbon stars are identified from multiple line indices measured from the stellar spectra. Combining the CN bands in the red end with \ctwo\ and other lines, we are able to identify the carbon stars. Moreover, we also classify the carbon stars into spectral sub-types of \ch, \CR, and \cn. These sub-types approximately show distinct features in the multi-dimensional line indices, implying that in the future we can use them to identify carbon stars from larger spectroscopic datasets. Meanwhile, from the line indices space, while the \cn\ stars are clearly separated from the others, we find no clear separation between \CR\ and \ch\ sub-types. The \CR\ and \ch\ stars seem to smoothly transition from one to another. This may hint that the \CR\ and \ch\ stars may not be different in their origins but look different in their spectra because of different metallicity. Due to the relatively low spectral resolution and lower signal-to-noise ratio, the ratio of $^{12}$C/$^{13}$C is not measured and thus the \cj\ stars are not identified.
The sample of 37 485 suspected OB stars selected by Gontcharov (2008) from the Tycho-2 catalogue has been cleaned of the stars that are not of spectral types OV--A0V. For this purpose, the apparent magnitude $V_T$ from Tycho-2, the absolute magnitude $M_{V_T}$ calibrated as a function of the dereddened color index $(B_T-V_T)_0$, the interstellar extinction $A_{V_T}$ calculated from the 3D analytical model by Gontcharov (2009) as a function of the Galactic coordinates, and the photometric distance $r_{ph}$ calculated as a function of $V_T$, $M_{V_T}$, and $A_{V_T}$ have been reconciled in an iterative process. The 20 514 stars that passed the iterations have $(B_T-V_T)_0<0$ and $M_{V_T}>-5$ and are considered as a sample of OV--A0V stars complete within 350 pc of the Sun. Based on the theoretical relation between the dereddened color and age of the stars, the derived sample has been divided into three subsamples: $(B_T-V_T)_0<-0.2^m$, $-0.2^m<(B_T-V_T)_0<-0.1^m$, and $-0.1^m<(B_T-V_T)_0<0^m$, younger than 100, $100-200$, and $200-400$ Myr, respectively. The spatial distribution of all 20 514 stars and the kinematics analyzed for more than 1500 stars with radial velocities from the PCRV and RAVE catalogues are different for the subsamples, showing smooth rotations, shears, and deformations of the layer of gas producing stars with the formation of the Gould Belt, the Great Tunnel, the Local Bubble, and other structures within the last 200 Myr. The detected temporal variations of the velocity dispersions, solar motion components, Ogorodnikov--Milne model parameters, and Oort constants are significant, agree with the results of other authors, and show that it is meaningless to calculate the kinematic parameters for samples of stars with uncertain ages or with a wide range of ages.
A major hurdle for modified gravity theories is to explain the dynamics of galaxy clusters. This paper makes the case for a generalized gravitational theory called Scalar-Tensor-Vector-Gravity (STVG) to explain merging cluster dynamics, and it will be the first of a series of papers intended to investigate this issue. The paper presents the results of a re-analysis of the Bullet Cluster as well as an analysis of the Train Wreck Cluster (using data from Jee et al. and Harvey et al.) in the weak gravitational field limit without dark matter. The King-$\beta$ model is used to fit the X-ray data of both clusters, and the $\kappa$-maps are computed using the parameters of this fit. The amount of galaxies in the clusters is estimated by subtracting the predicted $\kappa$-map from the $\kappa$-map data. The estimate suggests that $3.2\%$ of the Bullet Cluster is composed of galaxies. For the Train Wreck Cluster, if the Jee et al. data is used, it is found that $3.5\%$ of it is galaxies, and $22\%$ if the Harvey et al. data is used. The matter in galaxies and the enhanced gravity shift the lensing peaks making the peaks offset from the X-ray gas. The work demonstrates that this generalized gravitational theory has the potential to explain merging cluster dynamics without dark matter.
A significant fraction of nearby young moving group members harbor circumstellar debris dust disks. Due to their proximity and youth, these disks are attractive targets for studying the early evolution of debris dust and planetesimal belts. Here we present 70 and 160$\mu$m observations of 31 systems in the $\beta$ Pic moving group, and in the Tucana-Horologium, Columba, Carina and Argus associations, using the Herschel Space Observatory. None of these stars were observed at far-infrared wavelengths before. Our Herschel measurements were complemented by photometry from the WISE satellite for the whole sample, and by submillimeter/millimeter continuum data for one source, HD 48370. We identified six stars with infrared excess, four of them are new discoveries. By combining our new findings with results from the literature, we examined the incidence and general characteristics of debris disks around Sun-like members of the selected groups. With their dust temperatures of <45 K the newly identified disks around HD 38397, HD 48370, HD 160305, and BD-20 951 represent the coldest population within this sample. For HD 38397 and HD 48370, the emission is resolved in the 70$\mu$m PACS images, the estimated radius of these disks is ~90 au. Together with the well-known disk around HD 61005, these three systems represent the highest mass end of the known debris disk population around young G-type members of the selected groups. In terms of dust content, they resemble the hypothesized debris disk of the ancient Solar System.
Recently Stanchev et al. (2015) proposed a technique to derive density scaling relations in a star-forming region from analysis of the probability distribution function of column density. We address the possible dependence of the outcome on the selection of probe zones, applying the method to Planck dust-opacity data on Orion A. The derived steep scaling relation of mean density with index -1.6 in the molecular cloud (so called `Central filament') points to its self-gravitating nature. The result is reproduced also for large parts of the clouds' vicinity which indicates major role of gravity in the energy balance of the entire star-forming region.
Links to: arXiv, form interface, find, astro-ph, recent, 1606, contact, help (Access key information)
Hypervelocity stars (HVSs) that are observed in the Galactic halo, are believed to be accelerated to large velocities by a process of tidal disruption of binary stars passing close to a supermassive black hole (SMBH) which resides in the center of the Galaxy. It is, however, still unclear, where these relatively young stars were born and which dynamical process pushed them to nearly radial orbits around the SMBH. In this paper we investigate the possibility that the young binaries originated from a thin eccentric disc, similar to the one observed in the Galactic center nowadays. By means of direct N-body simulations, we follow the dynamical evolution of an initially thin and eccentric disc of stars with a 100% binary fraction orbiting around the SMBH. Such a configuration leads to Kozai-Lidov oscillations of orbital elements, bringing considerable amount of binaries to close vicinity of the black hole. Subsequent tidal disruption of these binaries accelerates one of their component to velocities well above the escape velocity from the SMBH while the second component becomes tightly bound to the SMBH. We describe the main kinematic properties of the escaping and tightly bound stars within our model and compare them qualitatively to the properties of the observed HVSs and S-stars, respectively. The most prominent feature is a strong anisotropy in the directions of the escaping stars which is observed for the Galactic HVSs but not explained yet.
We introduce SPARC (Spitzer Photometry & Accurate Rotation Curves): a sample of 175 nearby galaxies with new surface photometry at 3.6 um and high-quality rotation curves from previous HI/Halpha studies. SPARC spans a broad range of morphologies (S0 to Irr), luminosities (~5 dex), and surface brightnesses (~4 dex). We derive [3.6] surface photometry and study structural relations of stellar and gas disks. We find that both the stellar mass-HI mass relation and the stellar radius-HI radius relation have significant intrinsic scatter, while the HI mass-radius relation is extremely tight. We build detailed mass models and quantify the ratio of baryonic-to-observed velocity (Vbar/Vobs) for different characteristic radii and values of the stellar mass-to-light ratio (M/L) at [3.6]. Assuming M/L=0.5 Msun/Lsun (as suggested by stellar population models) we find that (i) the gas fraction linearly correlates with total luminosity, (ii) the transition from star-dominated to gas-dominated galaxies roughly corresponds to the transition from spiral galaxies to dwarf irregulars in line with density wave theory; and (iii) Vbar/Vobs varies with luminosity and surface brightness: high-mass, high-surface-brightness galaxies are nearly maximal, while low-mass, low-surface-brightness galaxies are submaximal. These basic properties are lost for low values of M/L=0.2 Msun/Lsun as suggested by the DiskMass survey. The mean maximum-disk limit in bright galaxies is M/L=0.7 Msun/Lsun at [3.6]. The SPARC data are publicly available and represent an ideal test-bed for models of galaxy formation.
We present an analysis of the flow of metals through the circumgalactic medium in the Feedback in Realistic Environments (FIRE) simulations of galaxy formation, ranging from isolated dwarfs to $L*$. We find that nearly all metals produced in high-redshift galaxies are carried out in winds that escape the galaxy and reach $0.25 R_{\rm vir}$. When measured at $0.25 R_{\rm vir}$ the metallicity of outflows is greater than, but similar to the ISM metallicity. Many metals thus reside in a reservoir in the CGM. We find that the outflowing flux through $R_{\rm vir}$ is lower than that at $0.25 R_{\rm vir}$ by a factor of $\sim2-5$. Cooling and recycling from this reservoir determine the metal budget in the ISM. The inflow metallicity at $R_{\rm vir}$ is generally very low, but outflow and inflow metallicities are similar in the inner halo. At low redshift, massive galaxies no longer generate outflows that reach the CGM, causing a divergence in CGM and ISM metallicity. Dwarf galaxies continue to generate outflows, which become increasingly dominated by metal-poor gas, while the galaxies themselves preferentially retains metal ejecta. In all but the least massive galaxy considered, a majority of the metals are within the halo at $z=0$. We measure the fraction of metals in CGM, ISM, stars, and roughly quantify the thermal state of CGM metals in each halo. The total amount of metals in the low-redshift CGM of two simulated L* galaxies is consistent with estimates from the COS halos survey, while for the other two it appears to be lower.
This paper is the second in a pair of articles presenting data release 1 (DR1) of the Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS), the largest single open-time key project carried out with the Herschel Space Observatory. The H-ATLAS is a wide-area imaging survey carried out in five photometric bands at 100, 160, 250, 350 and 500$\mu$m covering a total area of 600deg$^2$. In this paper we describe the identification of optical counterparts to submillimetre sources in DR1, comprising an area of 161 deg$^2$ over three equatorial fields of roughly 12$^\circ$x4.5$^\circ$ centred at 9$^h$, 12$^h$ and 14.5$^h$ respectively. Of all the H-ATLAS fields, the equatorial regions benefit from the greatest overlap with current multi-wavelength surveys spanning ultraviolet (UV) to mid-infrared regimes, as well as extensive spectroscopic coverage. We use a likelihood-ratio technique to identify SDSS counterparts at r<22.4 for 250-$\mu$m-selected sources detected at $\geq$ 4$\sigma$ ($\approx$28mJy). We find `reliable' counterparts (reliability R$\geq$0.8) for 44,835 sources (39 per cent), with an estimated completeness of 73.0 per cent and contamination rate of 4.7 per cent. Using redshifts and multi-wavelength photometry from GAMA and other public catalogues, we show that H-ATLAS-selected galaxies at $z<0.5$ span a wide range of optical colours, total infrared (IR) luminosities, and IR/UV ratios, with no strong disposition towards mid-IR-classified AGN in comparison with optical selection. The data described herein, together with all maps and catalogues described in the companion paper (Valiante et al. 2016), are available from the H-ATLAS website at www.h-atlas.org.
Andromeda II (And II) has been known for a few decades but only recently observations have unveiled new properties of this dwarf spheroidal galaxy. The presence of two stellar populations, the bimodal star formation history (SFH) and an unusual rotation velocity of And II put strong constrains on its formation and evolution. Following Lokas et al. (2014), we propose a detailed model to explain the main properties of And II involving (1) a gas-rich major merger between two dwarf galaxies at high redshift in the field and (2) a close interaction with M31 about 5 Gyr ago. The model is based on N-body/hydrodynamical simulations including gas dynamics, star formation and feedback. One simulation is designed to reproduce the gas-rich major merger explaining the origin of stellar populations and the SFH. Other simulations are used to study the effects of tidal forces and the ram pressure stripping during the interaction between And II and M31. The model successfully reproduces the SFH of And II including the properties of stellar populations, its morphology, kinematics and the lack of gas. Further improvements to the model are possible via joint modelling of all processes and better treatment of baryonic physics.
(Abridged) We present Herschel/HIFI spectra of the H2O 1113 GHz and H2O+ 1115
GHz lines toward five nearby prototypical starburst/AGN systems, and OH+ 971
GHz spectra toward three of these. The beam size of 20" corresponds to
resolutions between 0.35 and 7 kpc. The observed line profiles range from pure
absorption (NGC 4945, M82) to P-Cygni indicating outflow (NGC 253, Arp 220) and
inverse P-Cygni indicating infall (Cen A). The similarity of the H2O, OH+, and
H2O+ profiles to each other and to HI indicates that diffuse and dense gas
phases are well mixed. We estimate column densities assuming negligible
excitation (for absorption features) and using a non-LTE model (for emission
features), adopting calculated collision data for H2O and OH+, and rough
estimates for H2O+. Column densities range from ~10^13 to ~10^15 cm^-2 for each
species, and are similar between absorption and emission components, indicating
that the nuclear region does not contribute much to the emission in these
ground-state lines. The N(H2O)/N(H2O+) ratios of 1.4-5.6 indicate an origin of
the lines in diffuse gas, and the N(OH+)/N(H2O+) ratios of 1.6-3.1 indicate a
low H2 fraction (~11%) in the gas.
Adopting recent Galactic values for the average gas density and the
ionization efficiency, we find ionization rates for our sample galaxies of
~3x10^-16 s^-1 which are similar to the value for the Galactic disk, but ~10x
below that of the Galactic Center and ~100x below estimates for AGN from
excited-state H3O+ lines. We conclude that the ground-state lines of water and
its associated ions probe primarily non-nuclear gas in the disks of these
centrally active galaxies. Our data thus provide evidence for a decrease in
ionization rate by a factor of ~10 from the nuclei to the disks of galaxies, as
found before for the Milky Way.
We present LOFAR High-Band Array (HBA) observations of the Herschel-ATLAS North Galactic Pole survey area. The survey we have carried out, consisting of four pointings covering around 142 square degrees of sky in the frequency range 126--173 MHz, does not provide uniform noise coverage but otherwise is representative of the quality of data to be expected in the planned LOFAR wide-area surveys, and has been reduced using recently developed `facet calibration' methods at a resolution approaching the full resolution of the datasets ($\sim 10 \times 6$ arcsec) and an rms off-source noise that ranges from 100 $\mu$Jy beam$^{-1}$ in the centre of the best fields to around 2 mJy beam$^{-1}$ at the furthest extent of our imaging. We describe the imaging, cataloguing and source identification processes, and present some initial science results based on a 5-$\sigma$ source catalogue. These include (i) an initial look at the radio/far-infrared correlation at 150 MHz, showing that many Herschel sources are not yet detected by LOFAR; (ii) number counts at 150 MHz, including, for the first time, observational constraints on the numbers of star-forming galaxies; (iii) the 150-MHz luminosity functions for active and star-forming galaxies, which agree well with determinations at higher frequencies at low redshift, and show strong redshift evolution of the star-forming population; and (iv) some discussion of the implications of our observations for studies of radio galaxy life cycles.
This review summarises the main properties of Young Massive Clusters (YMCs), including their population properties, particularly focusing on extragalactic cluster samples. We discuss potential biases and caveats that can affect the construction of cluster samples and how incompleteness effects can result in erroneous conclusions regarding the long term survival of clusters. In addition to the luminosity, mass and age distributions of the clusters, we discuss the size distribution and profile evolution of the clusters. We also briefly discuss the stellar populations within YMCs. The final part of the review focusses on the connections between YMCs and the ancient globular clusters, whether or not they are related objects and how we can use what we know about YMC formation and evolution to understand how GCs formed in the early universe and how they relate to galaxy formation/evolution.
A population of quasars at z ~ 2 is determined based on dust luminosities vLv(7.8 um) that includes unobscured, partially obscured, and obscured quasars. Quasars are classified by the ratio vLv(0.25 um)/vLv(7.8 um) = UV/IR, assumed to measure obscuration of UV luminosity by the dust which produces IR luminosity. Quasar counts at rest frame 7.8 um are determined for quasars in the Bootes field of the NOAO Deep Wide Field Survey using 24 um sources with optical redshifts from the AGN and Galaxy Evolution Survey (AGES) or infrared redshifts from the Spitzer Infrared Spectrograph. Spectral energy distributions are extended to far infrared wavelengths using observations from the Herschel Space Observatory Spectral and Photometric Imaging Receiver (SPIRE), and new SPIRE photometry is presented for 77 high redshift quasars from the Sloan Digital Sky Survey. It is found that unobscured and obscured quasars have similar space densities at rest frame 7.8 um, but the ratio Lv(100 um)/Lv(7.8 um) is about three times higher for obscured quasars compared to unobscured, so that far infrared or submm discoveries are dominated by obscured quasars. Quasar source counts for these samples are determined for comparison to the number of submm sources that have been discovered with the SCUBA-2 camera at z ~ 2 using the Lv(100 um)/Lv(7.8 um) results together with the Bootes 7.8 um counts, and we find that only ~ 5% of high redshift submm sources are quasars, including even the most obscured quasars. Illustrative source counts are predicted to z = 10, and we show that existing SCUBA-2 850 um surveys or 2 mm surveys with the Goddard-IRAM Superconducting 2 Millimeter Observer (GISMO) survey camera should already have detected sources at z ~ 10 if quasar and starburst luminosity functions remain the same from z = 2 until z = 10.
Nuclear Clusters (NCs) are common stellar systems in the centres of galaxies. Yet, the physical mechanisms involved in their formation are still debated. Using a parsec-resolution hydrodynamical simulation of a dwarf galaxy, we propose an updated formation scenario for NCs. In this 'wet migration scenario', a massive star cluster forms in the gas-rich disc, keeping a gas reservoir, and growing further while it migrates to the centre via a combination of interactions with other substructures and dynamical friction. A wet merger with another dense cluster and its own gas reservoir can occur, although this is not a pre-requisite for the actual formation of the NC. The merging process does significantly alter the properties of the NC (mass, morphology, star formation history), also quenching the on-going local star formation activity, thus leading to interesting observational diagnostics for the physical origin of NCs. A population of lower mass clusters co-exist during the simulation, but these are either destroyed via tidal forces, or have high angular momentum preventing them to interact with the NC and contribute to its growth. The proposed updated scenario emphasises the role of gas reservoirs associated with the densest star clusters formed in a gas-rich low-mass galaxy.
The recent detections of GW150914 and GW151226 imply an abundance of stellar-mass binary-black-hole mergers in the local universe. While ground-based gravitational-wave detectors are limited to observing the final moments before a binary merges, space-based detectors, such as the Laser Interferometer Space Antenna (LISA), can observe binaries at lower orbital frequencies where such systems may still encode information about their formation histories. In particular, the orbital eccentricity of binary black holes in the LISA frequency band can be used discriminate between binaries formed in isolation in galactic fields, and those formed in dense stellar environments such as globular clusters. In this letter, we explore the differences in orbital eccentricities of binary black hole populations as they evolve through the LISA frequency band. Overall we find that there are three distinct populations of orbital eccentricities discernible by LISA. We show that, depending on gravitational-wave frequency, anywhere from 6% to 100% of binaries from globular clusters will be detectable with a measurable eccentricity. We note a bimodality in the eccentricities of binary black holes evolved in isolation that could provide detailed constraints on the physics of black-hole natal kicks and common-envelope evolution. Finally, we show how the measurement of both chirp mass and eccentricity can be used to constrain formation channels for binary black holes.
We present the first major data release of the largest single key-project in area carried out in open time with the Herschel Space Observatory. The Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS) is a survey of 600 deg^2 in five photometric bands - 100, 160, 250, 350 and 500 um - with the PACS and SPIRE cameras. In this paper and a companion paper (Bourne et al. 2016) we present the survey of three fields on the celestial equator, covering a total area of 161.6 deg^2 and previously observed in the Galaxy and Mass Assembly (GAMA) spectroscopic survey. This paper describes the Herschel images and catalogues of the sources detected on the SPIRE 250 um images. The 1-sigma noise for source detection, including both confusion and instrumental noise, is 7.4, 9.4 and 10.2 mJy at 250, 350 and 500 um. Our catalogue includes 120230 sources in total, with 113995, 46209 and 11011 sources detected at >4-sigma at 250, 350 and 500 um. The catalogue contains detections at >3-sigma at 100 and 160 um for 4650 and 5685 sources, and the typical noise at these wavelengths is 44 and 49 mJy. We include estimates of the completeness of the survey and of the effects of flux bias and also describe a novel method for determining the true source counts. The H-ATLAS source counts are very similar to the source counts from the deeper HerMES survey at 250 and 350 um, with a small difference at 500 um. Appendix A provides a quick start in using the released datasets, including instructions and cautions on how to use them.
CONTEXT. A number of recent papers have claimed the discovery of an X-shape
structure in the bulge of our Galaxy in the population of the red clumps.
AIMS. We endeavor to analyze the stellar density of bulge stars in the same
regions using a different stellar population that is characteristic of the
young bulge ($\lesssim 5$ Gyr). Particularly, we use F0-F5 main-sequence stars
with distances derived through photometric parallax.
METHODS. We extract these stars from extinction-corrected color-magnitude
diagrams in the near-infrared of VISTA-VVV data in some bulge regions and
calculate the densities along the line of sight. We take the uncertaintity in
the photometric parallax and the contamination of other sources into account,
and we see that these errors do not avoid the detection of a possible double
peak along some lines of sight as expected for a X-shape bulge if it existed.
RESULTS. Only a single peak in the density distribution along the line of
sight is observed, so apparently there is no X-shape structure for this
population of stars. Nonetheless, the effects of the dispersion of absolute
magnitudes in the selected population might be an alternative explanation,
although in principle these effects are insufficient to explain this lack of
double peak according to our calculations.
CONCLUSIONS. The results of the present paper do not demonstrate that
previous claims of X-shaped bulge using only red clump stars are incorrect, but
there are apparently some puzzling questions if we want to maintain the
validity of both the red-clump results and the results of this paper.
We report new IRAM/PdBI, JCMT/SCUBA-2, and VLA observations of the ultraluminous quasar SDSSJ010013.02+280225.8 (hereafter, J0100+2802) at z=6.3, which hosts the most massive supermassive black hole (SMBH) of 1.24x10^10 Msun known at z>6. We detect the [C II] 158 $\mu$m fine structure line and molecular CO(6-5) line and continuum emission at 353 GHz, 260 GHz, and 3 GHz from this quasar. The CO(2-1) line and the underlying continuum at 32 GHz are also marginally detected. The [C II] and CO detections suggest active star formation and highly excited molecular gas in the quasar host galaxy. The redshift determined with the [C II] and CO lines shows a velocity offset of ~1000 km/s from that measured with the quasar Mg II line. The CO (2-1) line luminosity provides direct constraint on the molecular gas mass which is about (1.0+/-0.3)x10^10 Msun. We estimate the FIR luminosity to be (3.5+/-0.7)x10^12 Lsun, and the UV-to-FIR spectral energy distribution of J0100+2802 is consistent with the templates of the local optically luminous quasars. The derived [C II]-to-FIR luminosity ratio of J0100+2802 is 0.0010+/-0.0002, which is slightly higher than the values of the most FIR luminous quasars at z~6. We investigate the constraint on the host galaxy dynamical mass of J0100+2802 based on the [C II] line spectrum. It is likely that this ultraluminous quasar lies above the local SMBH-galaxy mass relationship, unless we are viewing the system at a small inclination angle.
It has been suggested that the GeV excess, observed from the region surrounding the Galactic Center, might originate from a population of millisecond pulsars that formed in globular clusters. With this in mind, we employ the publicly available Fermi data to study the gamma-ray emission from 157 globular clusters, identifying a statistically significant signal from 25 of these sources (ten of which are not found in existing gamma-ray catalogs). We combine these observations with the predicted pulsar formation rate based on the stellar encounter rate of each globular cluster to constrain the gamma-ray luminosity function of millisecond pulsars in the Milky Way's globular cluster system. We find that this pulsar population exhibits a luminosity function that is quite similar to those millisecond pulsars observed in the field of the Milky Way (i.e. the thick disk). After pulsars are expelled from a globular cluster, however, they continue to lose rotational kinetic energy and become less luminous, causing their luminosity function to depart from the steady-state distribution. Using this luminosity function and a model for the globular cluster disruption rate, we show that millisecond pulsars born in globular clusters can account for only a few percent or less of the observed GeV excess. Among other challenges, scenarios in which the entire GeV excess is generated from such pulsars are in conflict with the observed mass of the Milky Way's Central Stellar Cluster.
We derived chemical abundances for C, N, O, Na, Mg and Al in 20 asymptotic giant branch (AGB) stars in the globular cluster NGC 6752. All these elements (but Mg) show intrinsic star-to-star variations and statistically significant correlations or anticorrelations analogous to those commonly observed in red giant stars of globular clusters hosting multiple populations. This demonstrates that, at odds with previous findings, both first and second generation stars populate the AGB of NGC 6752. The comparison with the Na abundances of red giant branch stars in the same cluster reveals that second generation stars (with mild Na and He enrichment) do reach the AGB phase. The only objects that are not observed along the AGB of NGC 6752 are stars with extreme Na enhancement. This is also consistent with standard stellar evolution models, showing that highly Na and He enriched stars populate the bluest portion of the horizontal branch and, because of their low stellar masses, evolve directly to the white dwarf cooling sequence, skipping the AGB phase.
Due to their heavily obscured central engines, the growth rate of Compton-thick (CT) active galactic nuclei (AGN) is difficult to measure. A statistically significant correlation between the Eddington ratio, {\lambda}$_{Edd}$, and the X-ray power-law index, {\Gamma}, observed in unobscured AGN offers an estimate of their growth rate from X-ray spectroscopy (albeit with large scatter). However, since X-rays undergo reprocessing by Compton scattering and photoelectric absorption when the line-of-sight to the central engine is heavily obscured, the recovery of the intrinsic {\Gamma} is challenging. Here we study a sample of local, predominantly Compton-thick megamaser AGN, where the black hole mass, and thus Eddington luminosity, are well known. We compile results on X-ray spectral fitting of these sources with sensitive high-energy (E> 10 keV) NuSTAR data, where X-ray torus models which take into account the reprocessing effects have been used to recover the intrinsic {\Gamma} values and X-ray luminosities, L$_X$. With a simple bolometric correction to L$_X$ to calculate {\lambda}$_{Edd}$, we find a statistically significant correlation between {\Gamma} and {\lambda}$_{Edd}$ (p = 0.007). A linear fit to the data yields {\Gamma} = (0.41$\pm$0.18)log$_{10}${\lambda}$_{Edd}$+(2.38$\pm$ 0.20), which is statistically consistent with results for unobscured AGN. This result implies that torus modeling successfully recovers the intrinsic AGN parameters. Since the megamasers have low-mass black holes (M$_{BH}\approx10^6-10^7$ M$_{sol}$) and are highly inclined, our results extend the {\Gamma}-{\lambda}$_{Edd}$ relationship to lower masses and argue against strong orientation effects in the corona, in support of AGN unification. Finally this result supports the use of {\Gamma} as a growth-rate indicator for accreting black holes, even for Compton-thick AGN.
We present a new approach aimed at constraining the typical size and optical properties of carbon dust grains in Circumstellar envelopes (CSEs) of carbon-rich stars (C-stars) in the Small Magellanic Cloud (SMC). To achieve this goal, we apply our recent dust growth description, coupled with a radiative transfer code to the CSEs of C-stars evolving along the TP-AGB, for which we compute spectra and colors. Then we compare our modeled colors in the near- and mid-infrared (NIR and MIR) bands with the observed ones, testing different assumptions in our dust scheme and employing several data sets of optical constants for carbon dust available in the literature. Different assumptions adopted in our dust scheme change the typical size of the carbon grains produced. We constrain carbon dust properties by selecting the combination of grain size and optical constants which best reproduces several colors in the NIR and MIR at the same time. The different choices of optical properties and grain size lead to differences in the NIR and MIR colors greater than two magnitudes in some cases. We conclude that the complete set of observed NIR and MIR colors are best reproduced by small grains, with sizes between $\sim$0.035 and $\sim$0.12~$\mu$m, rather than by large grains between $\sim0.2$ and $0.7$~$\mu$m. The inability of large grains to reproduce NIR and MIR colors seems independent of the adopted optical data set. We also find a possible trend of the grain size with mass-loss and/or carbon excess in the CSEs of these stars.
We report the serendipitous discovery of HSC J142449-005322, a double source plane lens system in the Hyper Suprime-Cam Subaru Strategic Program. We dub the system Eye of Horus. The lens galaxy is a very massive early-type galaxy with stellar mass of ~7x10^11 Msun located at z_L=0.795. The system exhibits two arcs/rings with clearly different colors, including several knots. We have performed spectroscopic follow-up observations of the system with FIRE on Magellan. The outer ring is confirmed at z_S2=1.988 with multiple emission lines, while the inner arc and counterimage is confirmed at z_S1=1.302. This makes it the first double source plane system with spectroscopic redshifts of both sources. Interestingly, redshifts of two of the knots embedded in the outer ring are found to be offset by delta_z=0.002 from the other knots, suggesting that the outer ring consists of at least two distinct components in the source plane. We perform lens modeling with two independent codes and successfully reproduce the main features of the system. However, two of the lensed sources separated by ~0.7 arcsec cannot be reproduced by a smooth potential, and the addition of substructure to the lens potential is required to reproduce them. Higher-resolution imaging of the system will help decipher the origin of this lensing feature and potentially detect the substructure.
Weak lensing statistics is typically measured as weighted sum of shear estimators or their products (shear-shear correlation). The weighting schemes are designed in the hope of minimizing the statistical error without introducing systematic errors. It would be ideal to approach the Cramer-Rao bound (the lower bound of the statistical uncertainty) in shear statistics, though it is generally difficult to do so in practice. The reasons may include: difficulties in galaxy shape measurement, inaccurate knowledge of the probability-distribution-function (PDF) of the shear estimator, misidentification of point sources as galaxies, etc.. Using the shear estimators defined in Zhang et al. (2015), we show that one can overcome all these problems, and allow shear measurement accuracy to approach the Cramer-Rao bound. This can be achieved by symmetrizing the PDF of the shear estimator, or the joint PDF of shear estimator pairs (for shear-shear correlation), without any prior knowledge of the PDF. Using simulated galaxy images, we demonstrate that under general observing conditions, this idea works as expected: it minimizes the statistical uncertainty without introducing systematic error.
Multicolor photometry from the Tycho-2 and 2MASS catalogues for 11 990 OB and 30 671 K-type red giant branch stars is used to detect systematic large-scale variations of the interstellar extinction law within the nearest kiloparsec. The characteristic of the extinction law, the total-to-selective extinction ratio $R_V$, which also characterizes the size and other properties of interstellar dust grains, has been calculated for various regions of space by the extinction law extrapolation method. The results for the two classes of stars agree: the standard deviation of the "red giants minus OB" $R_V$ differences within 500 pc of the Sun is 0.2. The detected $R_V$ variations between 2.2 and 4.4 not only manifest themselves in individual clouds but also span the entire space near the Sun, following Galactic structures. In the Local Bubble within about 100 pc of the Sun, $R_V$ has a minimum. In the inner part of the Gould Belt and at high Galactic latitudes, at a distance of about 150 pc from the Sun, $R_V$ reaches a maximum and then decreases to its minimum in the outer part of the Belt and other directions at a distance of about 500 pc from the Sun, returning to its mean values far from the Sun. The detected maximum of $R_V$ at high Galactic latitudes is important when allowance is made for the interstellar extinction toward extragalactic objects. In addition, a monotonic increase in $R_V$ by 0.3 per kpc toward the Galactic center has been found near the Galactic equator. It is consistent with the result obtained by Zasowski et al. (2009) for much of the Galaxy. Ignoring the $R_V$ variations and traditionally using a single value for the entire space must lead to systematic errors in the calculated distances reaching 10\%.
The product of the previously constructed 3D maps of stellar reddening (Gontcharov 2010) and $R_V$ variations (Gontcharov 2012) has allowed us to produce a 3D interstellar extinction map within the nearest kiloparsec from the Sun with a spatial resolution of 50 pc and an accuracy of $0.2^m$. This map is compared with the 2D reddening map by Schlegel et al. (1998), the 3D extinction map at high latitudes by Jones et al. (2011), and the analytical extinction models by Arenou et al. (1992) and Gontcharov (2009). In all cases, we have found good agreement and show that there are no systematic errors in the new map everywhere except the direction toward the Galactic center. We have found that the map by Schlegel et al. (1998) reaches saturation near the Galactic equator at $E_{(B-V)}>0.8^m$, has a zero-point error and systematic errors gradually increasing with reddening, and among the analytical models those that take into account the extinction in the Gould Belt are more accurate. Our extinction map shows that it is determined by reddening variations at low latitudes and $R_V$ variations at high ones. This naturally explains the contradictory data on the correlation or anticorrelation between reddening and $R_V$ available in the literature. There is a correlation in a thin layer near the Galactic equator, because both reddening and $R_V$ here increase toward the Galactic center. There is an anticorrelation outside this layer, because higher values of $R_V$ correspond to lower reddening at high and middle latitudes. Systematic differences in sizes and other properties of the dust grains in different parts of the Galaxy manifest themselves in this way. The largest structures within the nearest kiloparsec, including the Local Bubble, the Gould Belt, the Great Tunnel, the Scorpius, Perseus, Orion, and other complexes, have manifested themselves in the constructed map.
A new analytical 3D model of interstellar extinction within 500 pc of the Sun as a function of the Galactic spherical coordinates is suggested. This model is physically more justified than the widely used Arenou model, since it takes into account the presence of absorbing matter both in the layer along the equatorial Galactic plane and in the Gould Belt. The extinction in the equatorial layer varies as the sine of the Galactic longitude and in the Gould Belt as the sine of twice the longitude in the Belt plane. The extinction across the layers varies according to a barometric law. It has been found that the absorbing layers intersect at an angle of 17 deg and that the Sun is located near the axial plane of the absorbing layer of the Gould Belt and is probably several parsecs below the axial plane of the equatorial absorbing layer but above the Galactic plane. The model has been tested using the extinction of real stars from three catalogs.
Links to: arXiv, form interface, find, astro-ph, recent, 1606, contact, help (Access key information)