We constrain the evolution of the galaxy stellar mass function from 2 < z < 5 for galaxies with stellar masses as low as 10^5 Msun by combining star formation histories of Milky Way satellite galaxies derived from deep Hubble Space Telescope observations with merger trees from the ELVIS suite of N-body simulations. This approach extends our understanding more than two orders of magnitude lower in stellar mass than is currently possible by direct imaging. We find the faint end slopes of the mass functions to be alpha= -1.42(+0.07/-0.05) at z = 2 and alpha = -1.57^(+0.06/-0.06) at z = 5, and show the slope only weakly evolves from z = 5 to z = 0. Our findings are in stark contrast to a number of direct detection studies that suggest slopes as steep as alpha = -1.9 at these epochs. Such a steep slope would result in an order of magnitude too many luminous Milky Way satellites in a mass regime that is observationally complete (Mstar > 2*10^5 Msun at z = 0). The most recent studies from ZFOURGE and CANDELS also suggest flatter faint end slopes that are consistent with our results, but with a lower degree of precision. This work illustrates the strong connections between low and high-z observations when viewed through the lens of LCDM numerical simulations.
We present J-band spectroscopy of passive galaxies focusing on the Na I doublet at 1.14 {\mu}m. Like the Na I 0.82 {\mu}m doublet, this feature is strong in low-mass stars and hence may provide a useful probe of the initial mass function (IMF). From high signal-to-noise composite spectra, we find that Na I 1.14 {\mu}m increases steeply with increasing velocity dispersion, {\sigma}, and for the most massive galaxies (\sigma > 300 km/s) is much stronger than predicted from synthetic spectra with Milky-Way-like IMFs and solar abundances. Reproducing Na I 1.14 {\mu}m at high {\sigma} likely requires either a very high [Na/H], or a bottom-heavy IMF, or a combination of both. Using the Na D line to break the degeneracy between IMF and abundance, we infer [Na/H] $\approx$ +0.5 and a steep IMF (single-slope-equivalent x $\approx$ 3.2, where x = 2.35 for Salpeter), for the high-\sigma galaxies. At lower mass ({\sigma} = 50-100 km/s), the line strengths are compatible with MW-like IMFs and near-solar [Na/H]. We highlight two galaxies in our sample where strong gravitational lensing masses favour MW-like IMFs. Like the high-{\sigma} sample on average, these galaxies have strong Na I 1.14 \mu m; taken in isolation their sodium indices imply bottom-heavy IMFs which are hard to reconcile with the lensing masses. An alternative full-spectrum-fitting approach, applied to the high-\sigma sample, recovers an IMF less heavy than Salpeter, but under-predicts the Na I 1.14 \mu m line at the 5{\sigma} level. We conclude that current models struggle to reproduce this feature in the most massive galaxies without breaking other constraints, and caution against over-reliance on the sodium lines in spectroscopic IMF studies.
We use the GALFORM semi-analytical model to study high density regions traced by radio galaxies and quasars at high redshifts. We explore the impact that baryonic physics has upon the properties of galaxies in these environments. Star-forming emission-line galaxies (Ly{\alpha} and H{\alpha} emitters) are used to probe the environments at high redshifts. Radio galaxies are predicted to be hosted by more massive haloes than quasars, and this is imprinted on the amplitude of galaxy overdensities and cross-correlation functions. We find that Ly{\alpha} radiative transfer and AGN feedback indirectly affect the clustering on small scales and also the stellar masses, star- formation rates and gas metallicities of galaxies in dense environments. We also investigate the relation between protoclusters associated with radio galaxies and quasars, and their present- day cluster descendants. The progenitors of massive clusters associated with radio galaxies and quasars allow us to determine an average protocluster size in a simple way. Overdensities within the protoclusters are found to correlate with the halo descendant masses. We present scaling relations that can be applied to observational data. By computing projection effects due to the wavelength resolution of modern spectrographs and narrow-band filters we show that the former have enough spectral resolution to map the structure of protoclusters, whereas the latter can be used to measure the clustering around radio galaxies and quasars over larger scales to determine the mass of dark matter haloes hosting them.
We present the largest number of Milky Way sized dark matter halos simulated at very high mass ($\sim$$10^4$ M$_\odot$/particle) and temporal resolution ($\sim$5 Myrs/snapshot) done to date, quadrupling what is currently available in the literature. This initial suite consists of the first 24 halos of the $Caterpillar$ $Project$ (www.caterpillarproject.org) whose project goal of 60 - 70 halos will be made public when complete. We resolve $\sim$20,000 gravitationally bound subhalos within the virial radius of each host halo. Over the ranges set by our spatial resolution our convergence is excellent and improvements were made upon current state-of-the-art halo finders to better identify substructure at such high resolutions (e.g., on average we recover $\sim$4 subhalos in each host halo above 10$^8$ M$_\odot$ which would have otherwise not been found using conventional methods). For our relaxed halos, the inner profiles are reasonably fit by Einasto profiles ($\alpha$ = 0.169 $\pm$ 0.023) though this depends on the relaxed nature and assembly history of a given halo. Averaging over all halos, the substructure mass fraction is $f_{m,subs} = 0.121 \pm 0.041$, and mass function slope is d$N$/d$M\propto M^{-1.88 \pm 0.10}$ though we find scatter in the normalizations for fixed halo mass due to more concentrated hosts having less subhalos at fixed subhalo mass. There are no biases stemming from Lagrangian volume selection as all Lagrangian volume types are included in our sample. Our detailed contamination study of 264 low resolution halos has resulted in obtaining very large and unprecedented, high-resolution regions around our host halos for our target resolution (sphere of radius $\sim$$1.4 \pm 0.4$ Mpc) allowing for accurate studies of low mass dwarf galaxies at large galactocentric radii and the very first stellar systems at high redshift ($z \geq$ 10).
This is the fourth paper in a series that reports on our investigation of the clustering properties of active galactic nuclei (AGN) identified in the ROSAT All-Sky Survey (RASS) and Sloan Digital Sky Survey (SDSS). In this paper we investigate the cause of the X-ray luminosity dependence of the clustering of broad-line, luminous AGN at 0.16<z<0.36. We fit the H-alpha line profile in the SDSS spectra for all X-ray and optically-selected broad-line AGN, determine the mass of the super-massive black hole (SMBH), M_BH, and infer the accretion rate relative to Eddington (L/L_EDD). Since M_BH and L/L_EDD are correlated, we create AGN subsamples in one parameter while maintaining the same distribution in the other parameter. In both the X-ray and optically-selected AGN samples we detect a weak clustering dependence with M_BH and no statistically significant dependence on L/L_EDD. We find a difference of up to 2.7sigma when comparing the objects that belong to the 30% least and 30% most massive M_BH subsamples, in that luminous broad-line AGN with more massive black holes reside in more massive parent dark matter halos at these redshifts. These results provide evidence that higher accretion rates in AGN do not necessarily require dense galaxy environments in which more galaxy mergers and interactions are expected to channel large amounts of gas onto the SMBH. We also present semi-analytic models which predict a positive M_DMH dependence on M_BH, which is most prominent at M_BH ~ 10^{8-9} M_SUN.
This paper is the third in a series presenting the results of direct numerical integrations of the Fokker-Planck equation for stars orbiting a supermassive black hole (SBH) at the center of a galaxy. The algorithm of Paper II included diffusion coefficients that described the effects of random ("classical") and correlated ("resonant") relaxation. In this paper, the diffusion coefficients of Paper II have been generalized to account for the effects of "anomalous relaxation," the qualitatively different way in which eccentric orbits evolve in the regime of rapid relativistic precession. Two functional forms for the anomalous diffusion coefficients are investigated, based on power-law or exponential modifications of the resonant diffusion coefficients. The parameters defining the modified coefficients are first constrained by comparing the results of Fokker-Planck integrations with previously-published N-body integrations. Steady-state solutions are then obtained via the Fokker-Planck equation for models with properties similar to those of the Milky Way nucleus. Inclusion of anomalous relaxation leads to the formation of less prominent cores than in the case of resonant relaxation alone, due to the lengthening of diffusion timescales for eccentric orbits. Steady-state capture rates of stars by the SBH are found to always be less, by as much as an order of magnitude, than capture rates in the presence of resonant relaxation alone.
The Stripe 82 Massive Galaxy Catalog (S82-MGC) is the largest-volume stellar mass-limited sample of galaxies beyond z~1 constructed to date. Spanning 139.4 deg2, the S82-MGC includes a mass-limited sample of 41,770 galaxies with log Mstar > 11.2 to z~0.7, sampling a volume of 0.3 Gpc3, roughly equivalent to the volume of the Sloan Digital Sky Survey-I/II (SDSS-I/II) z < 0.15 MAIN sample. The catalog is built on three pillars of survey data: the SDSS Stripe 82 Coadd photometry which reaches r-band magnitudes of 23.5 AB, YJHK photometry at depths of 20th magnitude (AB) from the UK Infrared Deep Sky Survey (UKIDSS) Large Area Survey, and over 70,000 spectroscopic galaxy redshifts from SDSS-I/II and the Baryon Oscillation Spectroscopic Survey (BOSS). We describe the catalog construction and verification, the production of 9-band matched aperture photometry, tests of existing and newly estimated photometric redshifts required to supplement spectroscopic redshifts for 55% of the log Mstar > 11.2 sample, and geometric masking. We provide near-IR based stellar mass estimates and compare these to previous estimates. All catalog products are made publicly available. The S82-MGC not only addresses previous statistical limitations in high-mass galaxy evolution studies but begins tackling inherent data challenges in the coming era of wide-field imaging surveys.
We discuss spatially resolved emission line spectroscopy secured for a total sample of 15 gravitationally lensed star-forming galaxies at a mean redshift of $z\simeq2$ based on Keck laser-assisted adaptive optics observations undertaken with the recently-improved OSIRIS integral field unit (IFU) spectrograph. By exploiting gravitationally lensed sources drawn primarily from the CASSOWARY survey, we sample these sub-L$^{\ast}$ galaxies with source-plane resolutions of a few hundred parsecs ensuring well-sampled 2-D velocity data and resolved variations in the gas-phase metallicity. Such high spatial resolution data offers a critical check on the structural properties of larger samples derived with coarser sampling using multiple-IFU instruments. We demonstrate how serious errors of interpretation can only be revealed through better sampling. Although we include four sources from our earlier work, the present study provides a more representative sample unbiased with respect to emission line strength. Contrary to earlier suggestions, our data indicates a more diverse range of kinematic and metal gradient behavior inconsistent with a simple picture of well-ordered rotation developing concurrently with established steep metal gradients in all but merging systems. Comparing our observations with the predictions of hydrodynamical simulations suggests that strong feedback plays a key role in flattening metal gradients in early star-forming galaxies.
We present a new image of the 5.5 GHz radio emission from the extended Chandra Deep Field South. Deep radio observations at 5.5 GHz were obtained in 2010 and presented in the first data release. A further 76 hours of integration has since been obtained, nearly doubling the integration time. This paper presents a new analysis of all the data. The new image reaches 8.6 microJy rms, an improvement of about 40% in sensitivity. We present a new catalogue of 5.5 GHz sources, identifying 212 source components, roughly 50% more than were detected in the first data release. Source counts derived from this sample are consistent with those reported in the literature for S_{5.5GHz} > 0.1 mJy but significantly lower than published values in the lowest flux density bins (S_{5.5GHz} < 0.1 mJy), where we have more detected sources and improved statistical reliability. The 5.5 GHz radio sources were matched to 1.4 GHz sources in the literature and we find a mean spectral index of -0.35 +- 0.10 for S_{5.5GHz} > 0.5 mJy, consistent with the flattening of the spectral index observed in 5 GHz sub-mJy samples. The median spectral index of the whole sample is \alpha_{med} = -0.58, indicating that these observations may be starting to probe the star forming population. However, even at the faintest levels (0.05 < S_{5.5GHz} < 0.1 mJy), 39% of the 5.5 GHz sources have flat or inverted radio spectra. Four flux density measurements from our data, across the full 4.5 to 6.5 GHz bandwidth, are combined with those from literature and we find 10% of sources (S_{5.5GHz} >~ 0.1 mJy) show significant curvature in their radio spectral energy distribution spanning 1.4 to 9 GHz.
We compared optical spectroscopic and photometric data for 18 AGN galaxies over 2 to 3 epochs, with time intervals of typically 5 to 10 years. We used the Multi-Object Double Spectrograph (MODS) at the Large Binocular Telescope (LBT) and compared the spectra to data taken from the SDSS database and the literature. We find variations in the forbidden oxygen lines as well as in the hydrogen recombination lines of these sources. For 4 of the sources we find that, within the calibration uncertainties, the variations in continuum and line spectra of the sources are very small. We argue that it is mainly the difference in black hole mass between the samples that is responsible for the different degree of continuum variability. In addition we find that for an otherwise constant accretion rate the total line variability (dominated by the narrow line contributions) reverberates the continuum variability with a dependency $\Delta L_{line} \propto (\Delta L_{cont.})^{\frac{3}{2}}$. Since this dependency is prominently expressed in the narrow line emission it implies that the luminosity dominating part of the narrow line region must be very compact with a size of the order of at least 10 light years. A comparison to literature data shows that these findings describe the variability characteristics of a total of 61 broad and narrow line sources.
Double-lobe radio galaxies in the local Universe have traditionally been found to be hosted in elliptical or lenticular galaxies. We report the discovery of four spiral-host double-lobe radio galaxies (J0836+0532, J1159+5820, J1352+3126 and J1649+2635) that are discovered by cross-matching a large sample of 187005 spiral galaxies from SDSS DR7 to the full catalogues of FIRST and NVSS. J0836+0532 is reported for the first time. The host galaxies are forming stars at an average rate of 1.7 $-$ 10 M$_{\odot}$ yr$^{-1}$ and possess Super Massive Black Holes (SMBHs) with masses of a few times 10$^{8}$ M$_{\odot}$. Their radio morphologies are similar to FR-II radio galaxies with total projected linear sizes ranging from 86 kpc to 420 kpc, but their total 1.4 GHz radio luminosities are only in the range 10$^{24}$ $-$ 10$^{25}$ W Hz$^{-1}$. We propose that the formation of spiral-host double-lobe radio galaxies can be attributed to more than one factor, such as the occurrence of strong interactions, mergers, and the presence of unusually massive SMBHs, such that the spiral structures are not destroyed. Only one of our sources (J1649+2635) is found in a cluster environment, indicating that processes other than accretion through cooling flows e.g., galaxy-galaxy mergers or interactions could be plausible scenarios for triggering radio-loud AGN activity in spiral galaxies.
Metal-poor stars in the Milky Way are local relics of the epoch of the first stars and the first galaxies. However, a low metallicity does not prove that a star formed in this ancient era, as metal-poor stars form over a range of redshift in different environments. Theoretical models of Milky Way formation have shown that at constant metallicity, the oldest stars are those closest to the center of the Galaxy on the most tightly-bound orbits. For that reason, the most metal-poor stars in the bulge of the Milky Way provide excellent tracers of the chemistry of the high-redshift universe. We report the dynamics and detailed chemical abundances of three stars in the bulge with [Fe/H] $\lesssim-2.7$, two of which are the most metal-poor stars in the bulge in the literature. We find that with the exception of scandium, all three stars follow the abundance trends identified previously for metal-poor halo stars. These three stars have the lowest [Sc II/Fe] abundances yet seen in $\alpha$-enhanced giant stars in the Galaxy. Moreover, all three stars are outliers in the otherwise tight [Sc II/Fe]-[Ti II/Fe] relation observed among metal-poor halo stars. Theoretical models predict that there is a 30% chance that at least one of these stars formed at $z\gtrsim15$, while there is a 70% chance that at least one formed at $10 \lesssim z \lesssim 15$. These observations imply that by $z\sim10$, the progenitor galaxies of the Milky Way had both reached [Fe/H] $\sim-3.0$ and established the abundance pattern observed in extremely metal-poor stars.
The number of galaxies at a given flux as a function of the redshift, $z$, is derived when the $z$-distance relation is non-standard. In order to compare different models, the same formalism is also applied to the standard cosmology. The observed luminosity function for galaxies of the zCOSMOS catalog at different redshifts is modelled by a new luminosity function for galaxies, which is derived by the truncated beta probability density function. Three astronomical tests, which are the photometric maximum as a function of the redshift for a fixed flux, the mean value of the redshift for a fixed flux, and the luminosity function for galaxies as a function of the redshift, compare the theoretical values of the standard and non-standard model with the observed value. The tests are performed on the FORS Deep Field (FDF) catalog up to redshift $z=1.5$ and on the zCOSMOS catalog extending beyond $z=4$. These three tests show minimal differences between the standard and the non-standard models.
We used ultra-deep UV observations obtained with the Hubble Space Telescope to search for optical companions to binary millisecond pulsars (MSPs) in the globular cluster 47 Tucanae. We identified four new counterparts (to MSPs 47TucQ, 47TucS, 47TucT and 47TucY) and confirmed those already known (to MSPs 47TucU and 47TucW). In the color magnitude diagram, the detected companions are located in a region between the main sequence and the CO white dwarf cooling sequences, consistent with the cooling tracks of He white dwarfs of mass between 0.15 Msun and 0.20 Msun. For each identified companion, mass, cooling age, temperature and pulsar mass (as a function of the inclination angle) have been derived and discussed. For 47TucU we also found that the past accretion history likely proceeded in a sub-Eddington rate. The companion to the redback 47TucW is confirmed to be a non degenerate star, with properties particularly similar to those observed for black widow systems. Two stars have been identified within the 2-sigma astrometric uncertainty from the radio positions of 47TucH and 47TucI, but the available data prevent us from firmly assessing whether they are the true companions of these two MSPs.
We present wide field near-infrared photometry of 12 Galactic globular clusters, typically extending from the tip of the cluster red giant branch (RGB) to the main sequence turnoff. Using recent homogenous values of cluster distance, reddening and metallicity, the resulting photometry is directly compared to the predictions of several recent libraries of stellar evolutionary models. Of the sets of models investigated, Dartmouth and Victoria-Regina models best reproduce the observed RGB morphology, albeit with offsets in J-Ks color which vary in their significance in light of all sources of observational uncertainty. Therefore, we also present newly recalibrated relations between near-IR photometric indices describing the upper RGB versus cluster iron abundance as well as global metallicity. The influence of enhancements in alpha elements and helium are analyzed, finding that the former affect the morphology of the upper RGB in accord with model predictions. Meanwhile, the empirical relations we derive are in good agreement with previous results, and minor discrepancies can likely be attributed to differences in the assumed cluster distances and reddenings. In addition, we present measurements of the horizontal branch (HB) and RGB bump magnitudes, finding a non-negligible dependence of the near-IR HB magnitude on cluster metallicity. Lastly, we discuss the influence of assumed cluster distances, reddenings and metallicities on our results, finding that our empirical relations are generally insensitive to these factors to within their uncertainties.
Links to: arXiv, form interface, find, astro-ph, recent, 1509, contact, help (Access key information)
The Galactic center is the closest region in which we can study star formation under extreme physical conditions like those in high-redshift galaxies. We measure the temperature of the dense gas in the central molecular zone (CMZ) and examine what drives it. We mapped the inner 300 pc of the CMZ in the temperature-sensitive J = 3-2 para-formaldehyde (p-H$_2$CO) transitions. We used the $3_{2,1} - 2_{2,0} / 3_{0,3} - 2_{0,2}$ line ratio to determine the gas temperature in $n \sim 10^4 - 10^5 $cm$^{-3}$ gas. We have produced temperature maps and cubes with 30" and 1 km/s resolution and published all data in FITS form. Dense gas temperatures in the Galactic center range from ~60 K to > 100 K in selected regions. The highest gas temperatures T_G > 100 K are observed around the Sgr B2 cores, in the extended Sgr B2 cloud, the 20 km/s and 50 km/s clouds, and in "The Brick" (G0.253+0.016). We infer an upper limit on the cosmic ray ionization rate ${\zeta}_{CR} < 10^{-14}$ 1/s. The dense molecular gas temperature of the region around our Galactic center is similar to values found in the central regions of other galaxies, in particular starburst systems. The gas temperature is uniformly higher than the dust temperature, confirming that dust is a coolant in the dense gas. Turbulent heating can readily explain the observed temperatures given the observed line widths. Cosmic rays cannot explain the observed variation in gas temperatures, so CMZ dense gas temperatures are not dominated by cosmic ray heating. The gas temperatures previously observed to be high in the inner ~75 pc are confirmed to be high in the entire CMZ.
We report the detection of a strong Milky Way-type 2175 \AA$ $ extinction bump at $z$ = 2.1166 in the quasar spectrum towards SDSS J121143.42+083349.7 from the Sloan Digital Sky Survey (SDSS) Data Release 10. We conduct follow up observations with the Echelle Spectrograph and Imager (ESI) onboard the Keck-II telescope and the Ultraviolet and Visual Echelle Spectrograph (UVES) on the VLT. This 2175 \AA$ $ absorber is remarkable in that we simultaneously detect neutral carbon (C I), neutral chlorine (Cl I), and carbon monoxide (CO). It also qualifies as a damped Lyman alpha system. The J1211+0833 absorber is found to be metal-rich and has a dust depletion pattern resembling that of the Milky Way disk clouds. We use the column densities of the C I fine structure states and the C II/C I ratio (under the assumption of ionization equilibrium) to derive the temperature and volume density in the absorbing gas. A Cloudy photoionization model is constructed, which utilizes additional atoms/ions to constrain the physical conditions. The inferred physical conditions are consistent with a canonical cold (T $\sim$ 100 K) neutral medium with a high density ($n$(H I) $\sim$ 100 cm$^{-3}$) and a slightly higher pressure than the local interstellar medium. Given the simultaneous presence of C I, CO, and the 2175 \AA$ $ bump, combined with the high metallicity, high dust depletion level and overall low ionization state of the gas, the absorber towards J1211+0833 supports the scenario that the presence of the bump requires an evolved stellar population.
We report the first identification of a compact dense galaxy overdensity at $z\sim 8$ dubbed A2744z8OD. A2744z8OD consists of eight $Y$-dropout galaxies behind Abell 2744 that is originally pinpointed by Hubble Frontier Fields studies. However, no studies have, so far, derived basic physical quantities of structure formation or made comparisons with theoretical models. We obtain a homogeneous sample of dropout galaxies at $z\sim 8$ from eight field data of Hubble legacy images that are as deep as the A2744z8OD data. Using the sample, we find that a galaxy surface overdensity value of A2744z8OD is very high $\delta\simeq 130$ that is defined in a small circle of $6"$ ($\simeq 30$ physical kpc) radius. Because there is no such a large $\delta$ value reported for high-$z$ overdensities to date, A2744z8OD is a system clearly different from those found in previous high-$z$ overdensity studies. The total stellar mass of A2744z8OD is estimated to be $3.8\times 10^9 M_\odot$ that is as small as today's Large Magellanic Cloud. In the galaxy+structure formation models of Henriques et al. (2015), there exist a very similar overdensity, Modelz8OD, that is made of eight model dropout galaxies at $z\sim 8$ in a $6"$-radius circle. Eight out of seven galaxies in Modelz8OD form a filament within $\Delta z \sim 0.03$ that is elongated in the line of sight, which enhance $\delta$ of Modelz8OD. Modelz8OD is a progenitor of a today's cluster with $10^{14} M_\odot$, and more than a half of the seven Modelz8OD galaxies are merged into the brightest cluster galaxy (BCG) in the today's cluster. If Modelz8OD is a counterpart of A2744z8OD, the models suggest that A2744z8OD would be a forming cluster core of a today's $10^{14} M_\odot$ cluster that started formation earlier than the most of the other BCG progenitors at $z>12$.
The Extremely Large Telescopes currently under construction have a collecting area that is an order of magnitude larger than the present largest optical telescopes. For seeing-limited observations the performance will scale as the collecting area but, with the successful use of adaptive optics, for many applications it will scale as $D^4$ (where $D$ is the diameter of the primary mirror). Central to the success of the ELTs, therefore, is the successful use of multi-conjugate adaptive optics (MCAO) that applies a high degree correction over a field of view larger than the few arcseconds that limits classical adaptive optics systems. In this letter, we report on the analysis of crowded field images taken on the central region of the Galactic globular cluster NGC 1851 in $K_s$ band using GeMS at the Gemini South telescope, the only science-grade MCAO system in operation. We use this cluster as a benchmark to verify the ability to achieve precise near-infrared photometry by presenting the deepest $K_s$ photometry in crowded fields ever obtained from the ground. We construct a colour-magnitude diagram in combination with the F606W band from HST/ACS. As well as detecting the "knee" in the lower main sequence at $K_s\simeq20.5$, we also detect the double subgiant branch of NGC 1851, that demonstrates the high photometric accuracy of GeMS in crowded fields.
Two major features of the prestellar CMF are: 1) a broad peak below 1 Msun, presumably corresponding to a mean gravitational fragmentation scale, and 2) a characteristic power-law slope, very similar to the Salpeter slope of the stellar initial mass function (IMF) at the high-mass end. While recent Herschel observations have shown that the peak of the prestellar CMF is close to the thermal Jeans mass in marginally supercritical filaments, the origin of the power-law tail of the CMF/IMF at the high-mass end is less clear. Inutsuka (2001) proposed a theoretical scenario in which the origin of the power-law tail can be understood as resulting from the growth of an initial spectrum of density perturbations seeded along the long axis of filaments by interstellar turbulence. Here, we report the statistical properties of the line-mass fluctuations of filaments in nearby molecular clouds observed with Herschel using a 1-D power spectrum analysis. The observed filament power spectra were fitted by a power-law function $(P_{true}(s) \propto s^{\alpha})$ after removing the effect of beam convolution at small scales. A Gaussian-like distribution of power-spectrum slopes was found centered at -1.6, close to that of the one-dimensional velocity power spectrum generated by subsonic Kolomogorov turbulence (-1.67). An empirical correlation, $P^{0.5}(s_0) \propto <N_{\rm H_2}>^{1.4 \pm 0.1} $, was also found between the amplitude of each filament power spectrum $P(s_0)$ and the mean column density along the filament $<N_{\rm H_2}>$. Finally, the dispersion of line-mass fluctuations along each filament $\sigma_{\rm M_{line}}$ was found to scale with the physical length $L$ of the filament, roughly as $\sigma_{M_{line}} \propto L^{0.7}$. Overall, our results are consistent with the suggestion that the bulk of the CMF/IMF results from the gravitational fragmentation of filaments.
We present the detection of compact radio structures of fourteen radio-loud narrow line Seyfert 1 (NLS1) galaxies from Very Long Baseline Array observations at 5 GHz, which were performed in 2013. While 50\% of the sources of our sample show a compact core only, the remaining 50\% exhibit a core-jet structure. The measured brightness temperatures of the cores range from $10^{8.4}$ to $10^{11.4}$ K with a median value of $10^{10.1}$ K, indicating that the radio emission is from non-thermal jets, and that, likely, most sources are not strongly beamed, then implying a lower jet speed in these radio-loud NLS1 galaxies. In combination with archival data taken at multiple frequencies, we find that seven sources show flat or even inverted radio spectra, while steep spectra are revealed in the remaining seven objects. Although all these sources are very radio-loud with $R > 100$, their jet properties are diverse, in terms of their milli-arcsecond (mas) scale (pc scale) morphology and their overall radio spectral shape. The evidence for slow jet speeds (i.e., less relativistic jets), in combination with the low kinetic/radio power, may offer an explanation for the compact VLBA radio structure in most sources. The mildly relativistic jets in these high accretion rate systems are consistent with a scenario, where jets are accelerated from the hot corona above the disk by the magnetic field and the radiation force of the accretion disk. Alternatively, a low jet bulk velocity can be explained by low spin in the Blandford-Znajek mechanism.
Carbon monoxide (CO) is widely used as a tracer of molecular hydrogen (H2) in metal-rich galaxies, but is known to become ineffective in low metallicity dwarf galaxies. Atomic carbon has been suggested as a superior tracer of H2 in these metal-poor systems, but its suitability remains unproven. To help us to assess how well atomic carbon traces H2 at low metallicity, we have performed a series of numerical simulations of turbulent molecular clouds that cover a wide range of different metallicities. Our simulations demonstrate that in star-forming clouds, the conversion factor between [CI] emission and H2 mass, $X_{\rm CI}$, scales approximately as $X_{\rm CI} \propto Z^{-1}$. We recover a similar scaling for the CO-to-H2 conversion factor, $X_{\rm CO}$, but find that at this point in the evolution of the clouds, $X_{\rm CO}$ is consistently smaller than $X_{\rm CI}$, by a factor of a few or more. We have also examined how $X_{\rm CI}$ and $X_{\rm CO}$ evolve with time. We find that $X_{\rm CI}$ does not vary strongly with time, demonstrating that atomic carbon remains a good tracer of H2 in metal-poor systems even at times significantly before the onset of star formation. On the other hand, $X_{\rm CO}$ varies very strongly with time in metal-poor clouds, showing that CO does not trace H2 well in starless clouds at low metallicity.
Near-infrared observations of stellar orbits at the Galactic Center provide conclusive evidence for a massive black hole associated with the compact radio source Sgr A*. The astrometric reference frame for these observations is tied to a set of red giant stars, which are also detectable at radio wavelengths through SiO maser emission in their envelopes. We have improved the precision and long-term stability of this reference frame, in which Sgr A* is localized to within a factor 5 better than previously: ~0.17 mas in position (in 2009) and ~0.07 mas/yr in velocity. This improvement is the result of modeling and correcting optical distortion in the VLT/NACO imager to a sub-mas level and including new infrared and radio measurements, which now both span more than a decade in time. A further improvement will follow future observations and facilitate the detection of relativistic orbital effects.
Gas inflow feeds galaxies with low metallicity gas from the cosmic web, sustaining star formation across the Hubble time. We make a connection between these inflows and metallicity inhomogeneities in star-forming galaxies, by using synthetic narrow-band images of the Halpha emission line from zoom-in AMR cosmological simulations of galaxies with stellar masses of $M \simeq 10^9 $Msun at redshifts z=2-7. In $\sim$50\% of the cases at redshifts lower than 4, the gas inflow gives rise to star-forming, Halpha-bright, off-centre clumps. Most of these clumps have gas metallicities, weighted by Halpha luminosity, lower than the metallicity in the surrounding interstellar medium by $\sim$0.3 dex, consistent with observations of chemical inhomogeneities at high and low redshifts. Due to metal mixing by shear and turbulence, these metallicity drops are dissolved in a few disc dynamical times. Therefore, they can be considered as evidence for rapid gas accretion coming from cosmological inflow of pristine gas.
We propose an observational test for gravitationally recoiling supermassive black holes (BHs) in active galactic nuclei, based on a correlation between the velocities of BHs relative to their host galaxies, |\Delta v|, and their obscuring dust column densities, \Sigma_{dust} (both measured along the line of sight). Proxies for both quantities can be derived from spectral features of individual quasars. We use toy models for the distribution of recoil velocities, BH trajectories, and the geometry of obscuring dust tori in galactic centres, to simulate 2.5x10^5 random observations of recoiling quasars. BHs with recoil velocities comparable to the escape velocity from the galactic centre remain bound to the nucleus, and do not fully settle back to the centre of the torus due to dynamical friction in a typical quasar lifetime. We find that |\Delta v| and \Sigma_ {dust} for these BHs are positively correlated. For obscured (\Sigma_{dust}>0) and for partially obscured (0<\Sigma_{dust}<~2.3 g/m^2) quasars with |\Delta v|>=45 km/s, the sample correlation coefficient between log10(|\Delta v|) and \Sigma_{dust} is r_{45} = 0.28+/-0.02 and r_{45} = 0.13+/-0.02, respectively. Allowing for random +/-100 km/s errors in |\Delta v| unrelated to the recoil dilutes the correlation for the partially obscured quasars to r_{45} = 0.026+/-0.004 measured between |\Delta v| and \Sigma_{dust}. A random sample of >~3,500 obscured quasars with |\Delta v|>=45 km/s would allow rejection of the no-correlation hypothesis with 3 sigma significance 95% of the time. Finally, we find that the fraction of obscured quasars, F_{obs}(|\Delta v|), decreases with |\Delta v| from F_{obs}(<10 km/s)>~0.8 to F_{obs}(>10^3 km/s)<~0.4. This predicted trend can be compared to the observed fraction of type II quasars, and can further test combinations of recoil, trajectory, and dust torus models.
We extend our effort to calibrate stellar isochrones in the Johnson-Cousins ($BVI_C$) and the 2MASS ($JHK_s$) filter systems based on observations of well-studied open clusters. Using cool main-sequence (MS) stars in Praesepe, we define empirical corrections to the Lejeune et al. color-effective temperature ($T_{\rm eff}$) relations down to $T_{\rm eff} \sim 3600$ K, complementing our previous work based on the Hyades and the Pleiades. We apply empirically corrected isochrones to existing optical and near-infrared photometry of cool ($T_{\rm eff} \leq 5500$ K) and metal-rich ([Fe/H]$=+0.37$) MS stars in NGC 6791. The current methodology relies on an assumption that color-$T_{\rm eff}$ corrections are independent of metallicity, but we find that estimates of color-excess and distance from color-magnitude diagrams with different color indices converge on each other at the precisely known metallicity of the cluster. Along with a satisfactory agreement with eclipsing binary data in the cluster, we view the improved internal consistency as a validation of our calibrated isochrones at super-solar metallicities. For very cool stars ($T_{\rm eff} \leq 4800$ K), however, we find that $B - V$ colors of our models are systematically redder than the cluster photometry by $\sim0.02$ mag. We use color-$T_{\rm eff}$ transformations from the infrared flux method (IRFM) and alternative photometry to examine a potential color-scale error in the input cluster photometry. After excluding $B - V$ photometry of these cool MS stars, we derive $E(B - V)=0.105\pm0.014$, [M/H]$=+0.42\pm0.07$, $(m - M)_0 = 13.04\pm0.08$, and the age of $9.5\pm0.3$ Gyr for NGC 6791.
Light and spectroscopy are among the most important and frequently taught topics in introductory, college-level, general education astronomy courses. This is due to the fact that the vast majority of observational data studied by astronomers arrives at Earth in the form of light. While there are many processes by which matter can emit and absorb light, Astro 101 courses typically limit their instruction to the Bohr model of the atom and electron energy level transitions. In this paper, we report on the development of a new Lecture-Tutorial to help students learn about other processes that are responsible for the emission and absorption of light, namely molecular rotations, molecular vibrations, and the acceleration of charged particles by magnetic fields.
We compare the absolute gain photometric calibration of the Planck/HFI and Herschel/SPIRE instruments on diffuse emission. The absolute calibration of each of HFI and SPIRE relies on planet flux measurements and comparison with theoretical far-infrared emission models of planetary atmospheres. We measure the photometric cross calibration between the instruments at two overlapping bands, 545 GHz / 500 $\mu$m and 857 GHz / 350 $\mu$m. The SPIRE maps used have been processed in the Herschel Interactive Processing Environment (Version 12) and the HFI data are from the 2015 Public Data Release 2. For our study we used 15 large fields observed with SPIRE, which cover a total of about 120 deg^2. We have selected these fields carefully to provide a high signal-to-noise ratio, avoid residual systematics in the SPIRE maps, and span a wide range of surface brightness. The HFI maps are bandpass-corrected to match the emission observed by the SPIRE bandpasses. SPIRE maps are convolved to match the HFI beam and put on a common pixel grid. We measure the relative gain between the instruments using two methods in each field, pixel-to-pixel correlation and angular power spectrum measurements. Adopting the current SPIRE point-source to extended-emission conversion, we find that the two calibrations are in very good agreement. The SPIRE / HFI relative gains are 1.047 ($\pm$ 0.0069) and 1.003 ($\pm$ 0.0080) at 545 and 857 GHz, respectively. These relative gains deviate from unity by much less than the current uncertainty of the absolute extended emission calibration, which is about 6.4% and 9.5% for HFI and SPIRE, respectively, but the deviations are comparable to the values 1.4% and 5.5% for HFI and SPIRE if the uncertainty from models of the common calibrator can be discounted. Of the 5.5% for SPIRE, 4% arises from the beam area, highlighting that as focus for refinement.
We present a photometric catalog of 8,735,004 proper motion selected low-mass stars (KML-spectral types) within the Sloan Digital Sky Survey (SDSS) footprint, from the combined SDSS Data Release 10 (DR10), Two-Micron All-Sky Survey (2MASS) Point Source Catalog (PSC), and Wide-field Infrared Survey Explorer (WISE) AllWISE catalog. Stars were selected using $r-i$, $i-z$, $r-z$, $z-J$, and $z-W1$ colors, and SDSS, WISE, and 2MASS astrometry was combined to compute proper motions. The resulting 3,518,150 stars were augmented with proper motions for 5,216,854 earlier type stars from the combined SDSS and United States Naval Observatory B1.0 catalog (USNO-B). We used SDSS+USNO-B proper motions to determine the best criteria for selecting a clean sample of stars. Only stars whose proper motions were greater than their $2$$\sigma$ uncertainty were included. Our Motion Verified Red Stars (MoVeRS) catalog is available through SDSS CasJobs and VizieR.
We present confirmation of the cluster MOO J1142+1527, a massive galaxy cluster discovered as part of the Massive and Distant Clusters of WISE Survey. The cluster is confirmed to lie at $z=1.19$, and using the Combined Array for Research in Millimeter-wave Astronomy we robustly detect the Sunyaev-Zel'dovich (SZ) decrement at 13.2$\sigma$. The SZ data imply a mass of $\mathrm{M}_{200m}=(1.1\pm0.2)\times10^{15}$ $\mathrm{M}_\odot$, making MOO J1142+1527 the most massive galaxy cluster known at $z>1.15$ and the second most massive cluster known at $z>1$. For a standard $\Lambda$CDM cosmology it is further expected to be one of the $\sim 5$ most massive clusters expected to exist at $z\ge1.19$ over the entire sky. Our ongoing Spitzer program targeting $\sim1750$ additional candidate clusters will identify comparably rich galaxy clusters over the full extragalactic sky.
Links to: arXiv, form interface, find, astro-ph, recent, 1509, contact, help (Access key information)
We measure the projected 2-point correlation function of galaxies in the 180 deg$^2$ equatorial regions of the GAMA II survey, for four different redshift slices between z = 0.0 and z=0.5. To do this we further develop the Cole (2011) method of producing suitable random catalogues for the calculation of correlation functions. We find that more r-band luminous, more massive and redder galaxies are more clustered. We also find that red galaxies have stronger clustering on scales less than ~3 $h^{-1}$ Mpc. We compare to two different versions of the GALFORM galaxy formation model, Lacey et al (in prep.) and Gonzalez-Perez et al. (2014), and find that the models reproduce the trend of stronger clustering for more massive galaxies. However, the models under predict the clustering of blue galaxies, can incorrectly predict the correlation function on small scales and under predict the clustering in our sample of galaxies with ~3$L_r$ . We suggest possible avenues to explore to improve these cluster- ing predictions. The measurements presented in this paper can be used to test other galaxy formation models, and we make the measurements available online to facilitate this.
We present an end-to-end, two-phase model for the origin of globular clusters (GCs). In the model, populations of stellar clusters form in the high-pressure discs of high-redshift ($z>2$) galaxies (a rapid-disruption phase due to tidal perturbations from the dense interstellar medium), after which the galaxy mergers associated with hierarchical galaxy formation redistribute the surviving, massive clusters into the galaxy haloes, where they remain until the present day (a slow-disruption phase due to tidal evaporation). The high galaxy merger rates of $z>2$ galaxies allow these clusters to be `liberated' into the galaxy haloes before they are disrupted within the high-density discs. This physically-motivated toy model is the first to include the rapid-disruption phase, which is shown to be essential for simultaneously reproducing the wide variety of properties of observed GC systems, such as their universal characteristic mass-scale, the dependence of the specific frequency on metallicity and galaxy mass, the GC system mass-halo mass relation, the constant number of GCs per unit supermassive black hole mass, and the colour bimodality of GC systems. The model predicts that most of these observables were already in place at $z=1$-$2$, although under rare circumstances GCs may still form in present-day galaxies. In addition, the model provides important constraints on models for multiple stellar populations in GCs by putting limits on initial GC masses and the amount of pristine gas accretion. The paper is concluded with a discussion of these and several other predictions and implications, as well as the main open questions in the field.
We study high-resolution hydrodynamic simulations of Milky Way type galaxies obtained within the "Evolution and Assembly of GaLaxies and their Environments" (EAGLE) project, and identify the those that best satisfy observational constraints on the Milky Way total stellar mass, rotation curve, and galaxy shape. Contrary to mock galaxies selected on the basis of their total virial mass, the Milky Way analogues so identified consistently exhibit very similar dark matter profiles inside the solar circle, therefore enabling more accurate predictions for indirect dark matter searches. We find in particular that high resolution simulated haloes satisfying observational constraints exhibit, within the inner few kiloparsecs, dark matter profiles shallower than those required to explain the so-called Fermi GeV excess via dark matter annihilation.
Recent Spitzer/InfraRed Array Camera (IRAC) photometric observations have revealed that rest-frame optical emission lines contribute signficantly to the broadband fluxes of high-redshift galaxies. Specifically, in the narrow redshift range z~5.1-5.4 the [3.6]-[4.5] color is expected to be very red, due to contamination of the 4.5-micron band by the dominant Halpha line, while the 3.6-micron filter is free of nebular emission lines. We take advantage of new reductions of deep Spitzer/IRAC imaging over the GOODS-North+South fields (Labbe+2015) to obtain a clean measurement of the mean Halpha equivalent width from the [3.6]-[4.5] color in the redshift range z=5.1-5.4. The selected sources either have measured spectroscopic redshifts (13 sources) or lie very confidently in the redshift range z=5.1-5.4 based on the photometric redshift likelihood intervals (11 sources). Our z_{phot}=5.1-5.4 sample and z_{spec}=5.10-5.40 spectroscopic sample have a mean [3.6]-[4.5] color of 0.31+/-0.05 mag and 0.35+/-0.07 mag, implying a rest-frame equivalent width EW(Halpha+[NII]+[SII]) of 665+/-53 Angstroms and 707+/-74 Angstroms, respectively, for sources in these samples. These values are consistent albeit slightly higher than derived by Stark+2013 at z~4, suggesting an evolution to higher values of the Halpha+[NII]+[SII] EW at z>2. Using the 3.6micron band, which is free of emission line contamination, we perform robust SED fitting and find a median specific star formation rate of sSFR = 17_{-5}^{+2} Gyr^{-1}, 7_{-2}^{+1}x higher than at z~2. We find no strong correlation (<2sigma) between the Halpha+[NII]+[SII] EW and the stellar mass of sources. Before the advent of JWST, improvements in these results will come through an expansion of current spectroscopic samples and deeper Spitzer/IRAC measurements.
We analysed 385 galactic spectra from the Sloan Digital Sky Survey Data Release 7 (SDSS-DR7) that belong to the catalog of isolated pairs of galaxies by Karachentsev. The spectra corresponds to physical pairs of galaxies as defined by V $\leq$ 1200 Km/s and a pair separation $\leq$ 100 kpc. We search for the incidence of nuclear activity, both thermal (star-forming) and non-thermal -Active Galactic Nuclei (AGN). After a careful extraction of the nuclear spectra, we use diagnostic diagrams and find that the incidence of AGN activity is 48 \% in the paired galaxies with emission lines and 40\% for the total sample (as compared to $\sim$ 43 \% and 41\% respectively in a sample of isolated galaxies). These results remain after dissecting the effects of morphological type and galactic stellar mass (with only a small, non significant, enhancement of the AGN fraction in pairs of objects). These results suggest that weak interactions are not necessary or sufficient to trigger low-luminosity AGN. Since the fraction of AGN is predominant in early type spiral galaxies, we conclude that the role of a bulge, and a large gas reservoir are both essential for the triggering of nuclear activity. The most striking result is that type 1 galaxies are almost absent from the AGN sample. This result is in conflict with the Unified Model, and suggests that high accretion rates are essential to form the Broad Line Region in active galaxies.
Variations in the photometric parameters of stellar systems as a function of their evolution and the stellar populations comprising them are investigated. A set of seven evolutionary models with an exponential decrease in the star-formation rate and 672 models with a secondary burst of star formation are considered. The occurrence of a secondary burst of star formation can shift the position of a stellar system on two-color diagrams to the right or left of the normal color sequence for galaxies and the extinction line. This makes it possible to estimate the composition of the stellar population of a galaxy with a nonmonotonic star formation history from its position on two-color diagrams. Surface photometry in both the optical (UBVRI) and near-IR (JHK) is used to study the stellar populations and star-formation histories in the structural components (nucleus, bulge, disk, spiral arms, bar, ring) of 26 galaxies of various morphological types (from S0 to Sd). Components (nucleus, bulge, bar) with color characteristics corresponding to stellar systems with secondary bursts of star formation are indicated in 10 of the 26 galaxies. The parameters of these secondary bursts are estimated. Five of the 10 galaxies with complex star-formation histories display clear structural perturbations. Appreciable differences in the photometric characteristics of relatively red early-type galaxies (S0-Sb) and relatively blue later-type galaxies (Sb-Sd) have been found. Galaxies of both early and late types are encountered among the Sb-galaxies. Lenticular galaxies do not display different photometric characteristics from early-type spiral galaxies.
Estimates of the mass distribution and dark-matter (DM) content of dwarf spheroidal galaxies (dSphs) are usually derived under the assumption that the effect of the tidal field of the host galaxy is negligible over the radial extent probed by kinematic data-sets. We assess the implications of this assumption in the specific case of the Fornax dSph by means of N-body simulations of a satellite orbiting around the Milky Way. We consider observationally-motivated orbits and we tailor the initial distributions of the satellite's stars and DM to match, at the end of the simulations, the observed structure and kinematics of Fornax. In all our simulations the present-day observable properties of Fornax are not significantly influenced by tidal effects. The DM component is altered by the interaction with the Galactic field (up to 20% of the DM mass within 1.6 kpc is lost), but the structure and kinematics of the stellar component are only mildly affected even in the more eccentric orbit (more than 99% of the stellar particles remain bound to the dwarf). In the simulations that successfully reproduce Fornax's observables, the dark-to-luminous mass ratio within 1.6 kpc is in the range 5-6, and up to 16-18 if measured within 3 kpc.
We present a theory in 3 parts, of the long-term (or secular) evolution of stellar systems orbiting within the sphere of influence of massive black holes in galactic nuclei. Here we describe the secular collisionless dynamics of a (Keplerian) stellar system of mass $M$ orbiting a black hole of mass $M_\bullet \gg M$. The stellar distribution function (DF) $f$ obeys the collisionless Boltzmann equation (CBE) in 6-dim phase space. The small mass ratio, $\varepsilon = M/M_\bullet \ll 1$, implies a separation of time scales in the motions of stars: the fast Kepler orbital periods and the secular time scale which is longer by a factor $\varepsilon^{-1}$. We orbit-average the CBE over the fast Keplerian orbital phase using the Method of Multiple Scales. Then $f$ is expressed as the sum of a secular DF $F$ in a 5-dim (Gaussian Ring) space, and small fluctuations that remain of $O(\varepsilon)$ over secular times. $F$ obeys a secular CBE that includes stellar self-gravity, general relativistic corrections up to 1.5 post-Newtonian order, and external sources. Secular dynamics conserves the semi-major axis of every star. This additional integral of motion promotes extra regularity of the stellar orbits, and enables the construction of secular equilibrium DFs ($F_0$) through a Secular Jeans theorem. Secular equilibria allow for varied spatial geometries including figure rotation. A linearized secular CBE determines the linear response and stability of $F_0$. Spherical, non-rotating equilibria may support small-amplitude, long-lived, warp-like distortions. We also prove that an axisymmetric, zero-thickness, flat disc is secularly stable to all in-plane perturbations, when its DF $F_0$ is a monotonic function of the angular momentum at fixed energy.
We present a first-principles theory of Resonant Relaxation (RR) of stellar systems orbiting within the sphere of influence of massive black holes in galactic nuclei. We extend the rigorous kinetic theory of Gilbert (1968) to include the Keplerian field of a black hole of mass $M_\bullet$, and specialize to a (Keplerian) stellar system of mass $M \ll M_\bullet$. Using the results of the secular collisionless theory of Paper I, we orbit-average the kinetic equation through perturbative development in the small parameter $\varepsilon = M/M_\bullet$. This is supplemented with contributions from general relativistic corrections up to 1.5 post-Newtonian order and external gravitational sources. The result is a kinetic equation for a secular distribution function (DF) in 5-dim (Gaussian Ring) space, with explicit forms for the fluctuation and dissipation components of the collision integral. For general DFs, both apsidal and nodal precessions contribute to RR; so the traditional, physically-motivated distinction between scalar-RR and vector-RR disappears. Irreversible 2-particle correlations, that build up through secular gravitational interactions, are the driving agents of RR. The correlation function can be written in terms of the wake function, which is the linear response of the system to the perturbation offered by any chosen stellar orbit. The relationship includes direct interactions between pairs of stars as well as collective effects (gravitational polarization). We discuss the interplay of secular dynamics and RR in the evolution toward secular thermodynamic equilibria. In Paper III we apply our RR theory to axisymmetric discs, and provide explicit formulae for the loss cone rates at which mass, energy and angular momentum are fed to the black hole.
Accurate star formation histories (SFHs) of galaxies are fundamental for
understanding the build-up of their stellar content. However, the most accurate
SFHs - those obtained from colour-magnitude diagrams (CMDs) of resolved stars
reaching the oldest main sequence turnoffs (oMSTO) - are presently limited to a
few systems in the Local Group. It is therefore crucial to determine the
reliability and range of applicability of SFHs derived from integrated light
spectroscopy, as this affects our understanding of unresolved galaxies from low
to high redshift.
To evaluate the reliability of current full spectral fitting techniques in
deriving SFHs from integrated light spectroscopy by comparing SFHs from
integrated spectra to those obtained from deep CMDs of resolved stars.
We have obtained a high signal--to--noise (S/N $\sim$ 36.3 per \AA)
integrated spectrum of a field in the bar of the Large Magellanic Cloud (LMC)
using EFOSC2 at the 3.6 meter telescope at La Silla Observatory. For this same
field, resolved stellar data reaching the oMSTO are available. We have compared
the star formation rate (SFR) as a function of time and the age-metallicity
relation (AMR) obtained from the integrated spectrum using {\tt STECKMAP}, and
the CMD using the IAC-star/MinnIAC/IAC-pop set of routines. For the sake of
completeness we also use and discuss other synthesis codes ({\tt STARLIGHT} and
{\tt ULySS}) to derive the SFR and AMR from the integrated LMC spectrum.
We find very good agreement (average differences $\sim$ 4.1 $\%$) between the
SFR(t) and the AMR obtained using {\tt STECKMAP} on the integrated light
spectrum, and the CMD analysis. {\tt STECKMAP} minimizes the impact of the
age-metallicity degeneracy and has the advantage of preferring smooth solutions
to recover complex SFHs by means of a penalized $\chi^2$. [abridged]
We investigate how the Milky Way tidal field can affect the spatial mixing of multiple stellar populations in the globular cluster NGC 6362. We use $N$-body simulations of multiple population clusters on the orbit of this cluster around the Milky Way. Models of the formation of multiple populations in globular clusters predict that the second population should initially be more centrally concentrated than the first. However, NGC 6362 is comprised of two chemically distinct stellar populations having the same radial distribution. We show that the high mass loss rate experienced on this cluster's orbit significantly accelerates the spatial mixing of the two populations expected from two body relaxation. We also find that for a range of initial second population concentrations, cluster masses, tidal filling factors and fraction of first population stars, a cluster with two populations should be mixed when it has lost 70-80 per cent of its initial mass. These results fully account for the complete spatial mixing of NGC 6362, since, based on its shallow present day mass function, independent studies estimate that the cluster has lost 85 per cent of its initial mass.
The Decadal IRAC Bootes Survey (DIBS) is a mid-IR variability survey of the ~9 sq. deg. of the NDWFS Bootes Field and extnds the time baseline of its predecessor, the Spitzer Deep, Wide-Field Survey (SDWFS), from 4 to 10 years. The Spitzer Space Telescope visited the field five times between 2004 and 2014 at 3.6 and 4.5 microns. We provide the difference image analysis photometry for a half a million mostly extragalactic sources. In the mid-IR color-color plane, sources with quasar colors constitute the largest variability class (75%), 16% of the variable objects have stellar colors and the remaining 9% have the colors of galaxies. Adding the fifth epoch doubles the number of variable AGNs for the same false positive rates as in SDWFS, or increases the number of sources by 20% while decreasing the false positive rates by factors of 2-3 for the same variability amplitude. We quantify the ensemble mid-IR variability of ~1500 spectroscopically confirmed AGNs using single power-law structure functions, which we find to be steeper (index $\gamma=0.45$) than in the optical ($\gamma=0.3$), leading to much lower amplitudes at short time-lags. This provides evidence for large emission regions, smoothing out any fast UV/optical variations, as the origin of infrared quasar variability. The mid-IR AGN structure function slope $\gamma$ seems to be uncorrelated with both the luminosity and rest-frame wavelength, while the amplitude shows an anti-correlation with the luminosity and a correlation with the rest-frame wavelength.
We have searched for continuous gravitational wave (CGW) signals produced by individually resolvable, circular supermassive black hole binaries (SMBHBs) in the latest EPTA dataset, which consists of ultra-precise timing data on 41 millisecond pulsars. We develop frequentist and Bayesian detection algorithms to search both for monochromatic and frequency-evolving systems. None of the adopted algorithms show evidence for the presence of such a CGW signal, indicating that the data are best described by pulsar and radiometer noise only. Depending on the adopted detection algorithm, the 95\% upper limit on the sky-averaged strain amplitude lies in the range $6\times 10^{-15}<A<1.5\times10^{-14}$ at $5{\rm nHz}<f<7{\rm nHz}$. This limit varies by a factor of five, depending on the assumed source position, and the most constraining limit is achieved towards the positions of the most sensitive pulsars in the timing array. The most robust upper limit -- obtained via a full Bayesian analysis searching simultaneously over the signal and pulsar noise on the subset of ours six best pulsars -- is $A\approx10^{-14}$. These limits, the most stringent to date at $f<10{\rm nHz}$, exclude the presence of sub-centiparsec binaries with chirp mass $\cal{M}_c>10^9$M$_\odot$ out to a distance of about 25Mpc, and with $\cal{M}_c>10^{10}$M$_\odot$ out to a distance of about 1Gpc ($z\approx0.2$). We show that state-of-the-art SMBHB population models predict $<1\%$ probability of detecting a CGW with the current EPTA dataset, consistent with the reported non-detection. We stress, however, that PTA limits on individual CGW have improved by almost an order of magnitude in the last five years. The continuing advances in pulsar timing data acquisition and analysis techniques will allow for strong astrophysical constraints on the population of nearby SMBHBs in the coming years.
We obtain predictions for the properties of cold dark matter annihilation radiation using high resolution hydrodynamic zoom-in cosmological simulations of Milky Way-like galaxies carried out as part of the "Evolution and Assembly of GaLaxies and their Environments" (EAGLE) programme. Galactic halos in the simulation have significantly different properties from those assumed by the "standard halo model" often used in dark matter detection studies. The formation of the galaxy causes a contraction of the dark matter halo, whose density profile develops a steeper slope than the Navarro-Frenk-White profile between $r\approx1.5~\rm{kpc}$ and $r\approx10~\rm{kpc}$, and a flatter slope at smaller radii. The inner regions of the halos are almost perfectly spherical (axis ratios $b/a > 0.96$ within $r=500~\rm{pc}$) and there is no offset larger than $45~\rm{pc}$ between the centre of the stellar distribution and the centre of the dark halo. The morphology of the predicted dark matter annihilation radiation signal is in broad agreement with $\gamma$-ray observations at large Galactic latitudes ($b\gtrsim3^\circ$). At smaller angles, the inferred signal in one of our four galaxies is similar to that which is observed but it is significantly weaker in the other three.
N-body simulations suggest that the substructures that survive inside dark matter haloes follow universal distributions in mass and radial number density. We demonstrate that a simple analytical model can explain these subhalo distributions as resulting from tidal stripping which increasingly reduces the mass of subhaloes with decreasing halo-centric distance. As a starting point, the spatial distribution of subhaloes of any given infall mass is shown to be largely indistinguishable from the overall mass distribution of the host halo. Using a physically motivated statistical description of the amount of mass stripped off individual subhaloes, the model fully describes the joint distribution of subhaloes in final mass, infall mass and radius. As a result, it can be used to predict several derived distributions involving combinations of these quantities including, but not limited to, the universal subhalo mass function, the subhalo spatial distribution, the lensing profile, the dark matter annihilation radiation profile and boost factor. This model clarifies a common confusion when comparing the spatial distributions of galaxies and subhaloes, the so called "anti-bias", as a simple selection effect. We provide a Python code SubGen for populating haloes with subhaloes at this http URL
(abridged) Context: The envelopes of molecular gas around embedded low-mass
protostars show different chemistries, which can be used to trace their
formation history and physical conditions. The excitation of some molecular
species can also be used to trace these physical conditions, making it possible
to constrain e.g. sources of heating and excitation.
Aims: To study the range of influence of an intermediate-mass Herbig Be
protostar, and to find what chemical and physical impact feedback effects from
the environment may have on embedded protostars.
Methods: We follow up on an earlier line survey of the Class 0/I source R CrA
IRS7B in the 0.8 mm window with an unbiased line survey of the same source in
the 1.3 mm window using the APEX telescope. We also study the excitation of the
key species H2CO, CH3OH, and c-C3H2 in a complete sample of the 18 embedded
protostars in the Corona Australis star-forming region. Radiative transfer
models are used to establish abundances of the molecular species.
Results: We detect line emission from 20 molecular species (32 including
isotopologues) in the two surveys. The most complex species detected are CH3OH,
CH3CCH, CH3CHO, and CH3CN. Several complex organics are significantly
under-abundant in comparison with "hot corino" protostars. The H2CO
temperatures of the sources in the region decrease with the distance to the
Herbig Be star R CrA, whereas the c-C3H2 temperatures remain constant across
the star-forming region.
Conclusions: The high H2CO temperatures observed towards objects close to R
CrA suggest that this star has a sphere of influence of several 10000 AU in
which it increases the temperature of the molecular gas to 30-50 K through
irradiation. The chemistry in the IRS7B envelope differs significantly from
many other embedded protostars, which could be an effect of the external
irradiation from R CrA.
Links to: arXiv, form interface, find, astro-ph, recent, 1509, contact, help (Access key information)
Ultraviolet and optical spectra of interstellar gas along the lines of sight to nearby stars have been interpreted by Redfield & Linsky (2008) and previous studies as a set of discrete warm, partially ionized clouds each with a different flow vector, temperature, and metal depletion. Recently, Gry & Jenkins (2014) have proposed a fundamentally different model consisting of a single cloud with nonrigid flows filling space out to 9 parsecs from the Sun that they propose better describes the local ISM. Here we test these fundamentally different morphological models against the spatially unbiased Malamut et al. (2014) spectroscopic data set, and find that the multiple cloud morphology model provides a better fit to both the new and old data sets. The detection of three or more velocity components along the lines of sight to many nearby stars, the presence of nearby scattering screens, the observed thin elongated structures of warm interstellar gas, and the likely presence of strong interstellar magnetic fields also support the multiple cloud model. The detection and identification of intercloud gas and the measurement of neutral hydrogen density in clouds beyond the Local Interstellar Cloud could provide future morphological tests.
Using cosmological simulations, we address the properties of high-redshift star-forming galaxies (SFGs) across their main sequence (MS) in the plane of star-formation rate (SFR) versus stellar mass. We relate them to the evolution of galaxies through phases of gas compaction, depletion, possible replenishment, and eventual quenching. We find that the high-SFR galaxies in the upper envelope of the MS are compact, with high gas fractions and short depletion times ("blue nuggets"), while the lower-SFR galaxies in the lower envelope have lower central gas densities, lower gas fractions and longer depletion times, consistent with observed gradients across the MS. Stellar-structure gradients are negligible. The SFGs oscillate about the MS ridge on timescales $\sim0.4~t_{\mathrm{Hubble}}$ ($\sim1$ Gyr at $z\sim3$). The propagation upwards is due to gas compaction, triggered, e.g., by mergers, counter-rotating streams, and/or violent disc instabilities. The downturn at the upper envelope is due to central gas depletion by peak star formation and outflows while inflow from the shrunken gas disc is suppressed. An upturn at the lower envelope can occur once the extended disc has been replenished by fresh gas and a new compaction can be triggered, namely as long as the replenishment time is shorter than the depletion time. The mechanisms of gas compaction, depletion and replenishment confine the SFGs to the narrow ($\pm0.3$ dex) MS. Full quenching occurs in massive haloes ($M_{\mathrm{vir}}>10^{11.5}~M_\odot$) and/or at low redshifts ($z<3$), where the replenishment time is long compared to the depletion time, explaining the observed bending down of the MS at the massive end.
We conduct a comprehensive numerical study of the orbital dependence of harassment on early-type dwarfs consisting of 168 different orbits within a realistic, Virgo-like cluster, varying in eccentricity and pericentre distance. We find harassment is only effective at stripping stars or truncating their stellar disks for orbits that enter deep into the cluster core. Comparing to the orbital distribution in cosmological simulations, we find that the majority of the orbits (more than three quarters) result in no stellar mass loss. We also study the effects on the radial profiles of the globular cluster systems of early-type dwarfs. We find these are significantly altered only if harassment is very strong. This suggests that perhaps most early-type dwarfs in clusters such as Virgo have not suffered any tidal stripping of stars or globular clusters due to harassment, as these components are safely embedded deep within their dark matter halo. We demonstrate that this result is actually consistent with an earlier study of harassment of dwarf galaxies, despite the apparent contradiction. Those few dwarf models that do suffer stellar stripping are found out to the virial radius of the cluster at redshift=0, which mixes them in with less strongly harassed galaxies. However when placed on phase-space diagrams, strongly harassed galaxies are found offset to lower velocities compared to weakly harassed galaxies. This remains true in a cosmological simulation, even when halos have a wide range of masses and concentrations. Thus phase-space diagrams may be a useful tool for determining the relative likelihood that galaxies have been strongly or weakly harassed.
We present new spectroscopic observations that are part of our continuing monitoring campaign of 88 quasars at z<0.7 whose broad H-beta lines are offset from their systemic redshifts by a few thousand km/s. These quasars have been considered candidates for hosting supermassive black hole binaries (SBHBs) by analogy with single-lined spectroscopic binary stars. We present the data and describe our improved analysis techniques, which include an extensive evaluation of uncertainties. We also present a variety of measurements from the spectra that are of general interest and will be useful in later stages of our analysis. Additionally, we take this opportunity to study the variability of the optical continuum and integrated flux of the broad H-beta line. We compare the variability properties of the SBHB candidates to those of a sample of typical quasars with similar redshifts and luminosities observed multiple times during the Sloan Digital Sky Survey. We find that the variability properties of the two samples are similar (variability amplitudes of 10-30% on time scales of approximately 1-7 years) and that their structure functions can be described by a common model with parameters characteristic of typical quasars. These results suggest that the broad-line regions of SBHB candidates have a similar extent as those of typical quasars. We discuss the implications of this result for the SBHB scenario and ensuing constraints on the orbital parameters.
Context: Pre-merger interactions between galaxies can induce significant changes in the morphologies and kinematics of the stellar and ISM components. Large amounts of gas and stars are often found to be disturbed or displaced as tidal debris. This debris then evolves, sometimes forming stars and occasionally tidal dwarf galaxies. Here we present results from our HI study of Arp 65, an interacting pair hosting extended HI tidal debris. Aims: In an effort to understand the evolution of tidal debris produced by interacting pairs of galaxies, including in situ star and tidal dwarf galaxy formation, we are mapping HI in a sample of interacting galaxy pairs. The Arp 65 pair is one of them. Methods: Our resolved HI 21 cm line survey is being carried out using the Giant Metrewave Radio Telescope (GMRT). We used our HI survey data as well as available SDSS optical, Spitzer infra-red and GALEX UV data to study the evolution of the tidal debris and the correlation of HI with the star-forming regions within it. Results: In Arp 65 we see a high impact pre-merger interaction involving a pair of massive galaxies (NGC 90 and NGC 93) that have a stellar mass ratio of ~ 1:3. The interaction, which probably occurred ~ 1.0 -- 2.5 $\times$ 10$^8$ yr ago, appears to have displaced a large fraction of the HI in NGC 90 (including the highest column density HI) beyond its optical disk. We also find extended ongoing star formation in the outer disk of NGC 90. In the major star-forming regions, we find the HI column densities to be ~ 4.7 $\times$ 10$^{20}$ cm$^{-2}$ or lower. But no signature of star formation was found in the highest column density HI debris, SE of NGC 90. This indicates conditions within the highest column density HI debris remain hostile to star formation and it reaffirms that high HI column densities may be a necessary but not sufficient criterion for star formation.
In previous studies, it has been shown that the long term time average jet power, $\overline{Q}$, is correlated with the spectral index in the extreme ultraviolet (EUV), $\alpha_{EUV}$ (defined by $F_{\nu} \sim \nu^{-\alpha_{EUV}}$ computed between 700\AA\, and 1100\AA\,). Larger $\overline{Q}$ tends to decrease the EUV emission. This is a curious relationship because it connects a long term average over $\sim 10^{6}$ years with an instantaneous measurement of the EUV. The EUV appears to be emitted adjacent to the central supermassive black hole and the most straightforward explanation of the correlation is that the EUV emitting region interacts in real time with the jet launching mechanism. Alternatively stated, the $\overline{Q}$ - $\alpha_{EUV}$ correlation is a manifestation of a contemporaneous (real time) jet power, $Q(t)$, correlation with $\alpha_{EUV}$. In order to explore this possibility, this paper considers the time variability of the strong radio jet of the quasar 1442+101 that is not aberrated by strong Doppler enhancement. This high redshift (z = 3.55) quasar is uniquely suited for this endeavor as the EUV is redshifted into the optical observing window allowing for convenient monitoring. More importantly, it is bright enough to be seen through the Lyman forest and its radio flux is strong enough that it has been monitored frequently. Quasi-simultaneous monitoring (five epochs spanning $\sim 40$ years) show that increases in $Q(t)$ correspond to decreases in the EUV as expected.
We present the SAMI Pilot Survey, consisting of integral field spectroscopy
of 106 galaxies across three galaxy clusters, Abell 85, Abell 168 and Abell
2399. The galaxies were selected by absolute magnitude to have $M_r<-20.25$
mag. The survey, using the Sydney-AAO Multi-object Integral field spectrograph
(SAMI), comprises observations of galaxies of all morphological types with 75\%
of the sample being early-type galaxies (ETGs) and 25\% being late-type
galaxies (LTGs). Stellar velocity and velocity dispersion maps are derived for
all 106 galaxies in the sample.
The $\lambda_{R}$ parameter, a proxy for the specific stellar angular
momentum, is calculated for each galaxy in the sample. We find a trend between
$\lambda_{R}$ and galaxy concentration such that LTGs are less concentrated
higher angular momentum systems, with the fast-rotating ETGs (FRs) more
concentrated and lower in angular momentum. This suggests that some dynamical
processes are involved in transforming LTGs to FRs, though a significant
overlap between the $\lambda_{R}$ distributions of these classes of galaxies
implies that this is just one piece of a more complicated picture.
We measure the kinematic misalignment angle, $\Psi$, for the ETGs in the
sample, to probe the intrinsic shapes of the galaxies. We find the majority of
FRs (83\%) to be aligned, consistent with them being oblate spheroids (i.e.
disks). The slow rotating ETGs (SRs), on the other hand, are significantly more
likely to show kinematic misalignment (only 38\% are aligned). This confirms
previous results that SRs are likely to be mildly triaxial systems.
The South Galactic Cap u-band Sky Survey (SCUSS) is a deep u-band imaging survey in the Southern Galactic Cap, using the 90Prime wide-field imager on the 2.3m Bok telescope at Kitt Peak. The survey observations started in 2010 and ended in 2013. The final survey area is about 5000 deg2 with a median 5-sigma point source limiting magnitude of about 23.2. This paper describes the survey data reduction process, which includes basic imaging processing, astrometric and photometric calibrations, image stacking, and photometric measurements. Survey photometry is performed on objects detected both on SCUSS u-band images and in the SDSS database. Automatic, aperture, point-spread function (PSF), and model magnitudes are measured on stacked images. Co-added aperture, PSF, and model magnitudes are derived from measurements on single-epoch images. We also present comparisons of the SCUSS photometric catalog with those of the SDSS and CFHTLS.
Dwarf spheroidal galaxies (dSphs) are mostly investigated in the Local Group. DSphs are difficult targets for observations because of their small size and very low surface brightness. Here we measure spectroscopic and photometric parameters of three candidates for isolated dSphs, KKH65=BTS23, KK180, and KK227, outside the Local Group. The galaxies are found to be of low metallicity and low velocity dispersion. They are among the lowest surface brightness objects in the Local Universe. According to the measured radial velocities, metallicities, and structural and photometric parameters, KKH65 and KK227 are representatives of the ultra-diffuse quenched galaxies. KKH65 and KK227 belong to the outer parts of the groups NGC3414 and NGC5371, respectively. KK180 is located in the Virgo cluster infall region.
We investigated the habitability of the Milky Way, making use of recent observational analysis on the prevalence of Earth-sized planets, in order to estimate where and when potentially habitable star systems may have formed over the course of the Galaxy's history. We were then able to estimate the age distribution of potential intelligent life in our Galaxy using our own evolution and the age of the Sun as a proxy. To do this we created a galactic chemical evolution model and applied the following habitability constraints to the Sun-like (G-type) stars formed in our model: an environment free from life-extinguishing supernovae, a high enough metallicity for Earth-sized planet formation and sufficient time for the evolution of complex life. We determined a galactic habitable zone as the region containing all the potentially habitable star systems in our model. Our galactic habitable zone contains stars formed between 11 and 3.8 billion years ago at radial distances of between 7 and 14 kiloparsecs. We found that most potentially habitable star systems are much older than the Sun and located farther from the galactic centre. By comparing the ages of these systems we estimated that 77% of potentially habitable star systems are on average 3.13 billion years older than the Sun. This suggests that any intelligent life in the Galaxy is likely to be incredibly more advanced than we are assuming that they have evolved under similar timescales than we have. Implications and limitations of our study are discussed.
To understand cosmic mass assembly in the Universe at early epochs, we primarily rely on measurements of stellar mass and star formation rate of distant galaxies. In this paper, we present stellar masses and star formation rates of six high-redshift ($2.8\leq z \leq 5.7$) dusty, star-forming galaxies (DSFGs) that are strongly gravitationally lensed by foreground galaxies. These sources were first discovered by the South Pole Telescope (SPT) at millimeter wavelengths and all have spectroscopic redshifts and robust lens models derived from ALMA observations. We have conducted follow-up observations, obtaining multi-wavelength imaging data, using {\it HST}, {\it Spitzer}, {\it Herschel} and the Atacama Pathfinder EXperiment (APEX). We use the high-resolution {\it HST}/WFC3 images to disentangle the background source from the foreground lens in {\it Spitzer}/IRAC data. The detections and upper limits provide important constraints on the spectral energy distributions (SEDs) for these DSFGs, yielding stellar masses, IR luminosities, and star formation rates (SFRs). The SED fits of six SPT sources show that the intrinsic stellar masses span a range more than one order of magnitude with a median value $\sim$ 5 $\times 10^{10}M_{\Sun}$. The intrinsic IR luminosities range from 4$\times 10^{12}L_{\Sun}$ to 4$\times 10^{13}L_{\Sun}$. They all have prodigious intrinsic star formation rates of 510 to 4800 $M_{\Sun} {\rm yr}^{-1}$. Compared to the star-forming main sequence (MS), these six DSFGs have specific SFRs that all lie above the MS, including two galaxies that are a factor of 10 higher than the MS. Our results suggest that we are witnessing the ongoing strong starburst events which may be driven by major mergers.
We present the discovery of the highest velocity CIV broad absorption line to date in the z=2.47 quasar SDSS J023011.28+005913.6, hereafter J0230. In comparing the public DR7 and DR9 spectra of J0230, we discovered an emerging broad absorption trough outflowing at~60,000 km/s. In pursuing follow up observations we discovered a second emergent broad absorption trough outflowing at ~40,000 km/s. We collected seven spectral epochs of J0230 that demonstrate emergent and rapidly (~10 days in the rest-frame) varying broad absorption. We investigate two possible scenarios that could cause these rapid changes: bulk motion and ionization variability. Given our multi-epoch data, a transverse motion scenario would likely be a flow-tube feature travelling across the emitting region at 8,000 < v (km/s) < 56,000. If ionization variability is the cause for the changes observed, the absorber either has n_e >=1540 cm^-3 and is at r_{eq} >= 1.37 kpc, or is at r < 1.37 kpc with no constraint on the density.
The shape of the probability distribution function (PDF) of molecular clouds is an important ingredient for modern theories of star formation and turbulence. Recently, several studies have pointed out observational difficulties with constraining the low column density (i.e. Av <1) PDF using dust tracers. In order to constrain the shape and properties of the low column density probability distribution function, we investigate the PDF of multiphase atomic gas in the Perseus molecular cloud using opacity-corrected GALFA-HI data and compare the PDF shape and properties to the total gas PDF and the N(H2) PDF. We find that the shape of the PDF in the atomic medium of Perseus is well described by a lognormal distribution, and not by a power-law or bimodal distribution. The peak of the atomic gas PDF in and around Perseus lies at the HI-H2 transition column density for this cloud, past which the N(H2) PDF takes on a powerlaw form. We find that the PDF of the atomic gas is narrow and at column densities larger than the HI-H2 transition the HI rapidly depletes, suggesting that the HI PDF may be used to find the HI-H2 transition column density. We also calculate the sonic Mach number of the atomic gas by using HI absorption line data, which yields a median value of Ms=4.0 for the CNM, while the HI emission PDF, which traces both the WNM and CNM, has a width more consistent with transonic turbulence.
Recent observations of the Lyman-alpha forest show large-scale spatial variations in the intergalactic Lyman-alpha opacity that grow rapidly with redshift at z>5, far in excess of expectations from empirically motivated models. Previous studies have attempted to explain this excess with spatial fluctuations in the ionizing background, but found that this required either extremely rare sources or problematically low values for the mean free path of ionizing photons. Here we report that much -- or potentially all -- of the observed excess likely arises from residual spatial variations in temperature that are an inevitable byproduct of a patchy and extended reionization process. The amplitude of opacity fluctuations generated in this way depends on the timing and duration of reionization. If the entire excess is due to temperature variations alone, the observed fluctuation amplitude favors a late-ending but extended reionization process that was roughly half complete by z~9 and that ended at z~6. In this scenario, the highest opacities occur in regions that reionized earliest, since they have had the most time to cool, while the lowest opacities occur in the warmer regions that reionized most recently. This correspondence potentially opens a new observational window into patchy reionization.
We present a realistic modeling of the NVSS number count dipole across the sky. The modeling relies on mock catalogues generated within the context of $\Lambda$CDM cosmology in the linear regime of structure formation. After removal of the solar motion dipole from the observation, the mocks show that the remaining signal is mostly (70\%) due to contribution from large scale structure within $\sim 500$Mpc ($z\sim0.1$). The amplitude of this contribution depends on the bias factor of the NVSS galaxies. An effective bias factor $b(z<0.1)<2.0$ is ruled out at the $\sim 2.8\sigma$ significance by the comparison between the model and observed NVSS dipole. The mismatch is toned down to $\sim 2.3\sigma$ for $b(z)=3$.
We developed a dual-linear-polarization HEMT (High Electron Mobility Transistor) amplifier receiver system of the 45-GHz band (hereafter Z45), and installed it in the Nobeyama 45-m radio telescope. The receiver system is designed to conduct polarization observations by taking the cross correlation of two linearly-polarized components, from which we process full-Stokes spectroscopy. We aim to measure the magnetic field strength through the Zeeman effect of the emission line of CCS ($J_N=4_3-3_2$) toward pre-protostellar cores. A linear-polarization receiver system has a smaller contribution of instrumental polarization components to the Stokes $V$ spectra than that of the circular polarization system, so that it is easier to obtain the Stokes $V$ spectra. The receiver has an RF frequency of 42 $-$ 46 GHz and an intermediate frequency (IF) band of 4$-$8 GHz. The typical noise temperature is about 50 K, and the system noise temperature ranges from 100 K to 150K over the frequency of 42 $-$ 46 GHz. The receiver system is connected to two spectrometers, SAM45 and PolariS. SAM45 is a highly flexible FX-type digital spectrometer with a finest frequency resolution of 3.81 kHz. PolariS is a newly-developed digital spectrometer with a finest frequency resolution of 60 Hz, having a capability to process the full-Stokes spectroscopy. The Half Power Beam Width (HPBW) of the beam was measured to be 37$"$ at 43 GHz. The main beam efficiency of the Gaussian main beam was derived to be 0.72 at 43 GHz. The SiO maser observations show that the beam pattern is reasonably round at about 10 \% of the peak intensity and the side-lobe level was less than 3 \% of the peak intensity. Finally, we present some examples of astronomical observations using Z45.
We present new multi-band (UBVI) time-series data of helium burning variables in the Carina dwarf spheroidal galaxy. The current sample includes 92 RR Lyrae-six of them are new identifications-and 20 Anomalous Cepheids, one of which is new identification. The analysis of the Bailey diagram shows that the luminosity amplitude of the first overtone component in double-mode variables is located along the long-period tail of regular first overtone variables, while the fundamental component is located along the short-period tale of regular fundamental variables. This evidence further supports the transitional nature of these objects. Moreover, the distribution of Carina double-mode variables in the Petersen diagram (P_1/P_0 vs P_0) is similar to metal-poor globulars (M15, M68), to the dwarf spheroidal Draco and to the Galactic Halo. This suggests that the Carina old stellar population is metal-poor and affected by a small spread in metallicity. We use trigonometric parallaxes for five field RR Lyrae stars to provide an independent estimate of the Carina distance using the observed reddening free Period--Wesenheit [PW, (BV)] relation. Theory and observations indicate that this diagnostic is independent of metallicity. We found a true distance modulus of \mu=20.01\pm0.02 (standard error of the mean) \pm0.05 (standard deviation) mag. We also provided independent estimates of the Carina true distance modulus using four predicted PW relations (BV, BI, VI, BVI) and we found: \mu=(20.08\pm0.007\pm0.07) mag, \mu=(20.06\pm0.006\pm0.06) mag, \mu=(20.07\pm0.008\pm0.08) mag and \mu=(20.06\pm0.006\pm0.06) mag. Finally, we identified more than 100 new SX Phoenicis stars that together with those already known in the literature (340) make Carina a fundamental laboratory to constrain the evolutionary and pulsation properties of these transitional variables.
A multifrequency campaign on the BL Lac object PG 1553+113 was organized by the Whole Earth Blazar Telescope (WEBT) in 2013 April-August, involving 19 optical, two near-IR, and three radio telescopes. The aim was to study the source behaviour at low energies during and around the high-energy observations by the Major Atmospheric Gamma-ray Imaging Cherenkov (MAGIC) telescopes in April-July. We also analyse the UV and X-ray data acquired by the Swift and XMM-Newton satellites in the same period. The WEBT and satellite observations allow us to detail the synchrotron emission bump in the source spectral energy distribution (SED). In the optical we found a general bluer-when-brighter trend. The X-ray spectrum remained stable during 2013, but a comparison with previous observations suggests that it becomes harder when the X-ray flux increases. The long XMM-Newton exposure reveals a curved X-ray spectrum. In the SED, the XMM-Newton data show a hard near-UV spectrum, while Swift data display a softer shape that is confirmed by previous HST-COS and IUE observations. Polynomial fits to the optical-X-ray SED show that the synchrotron peak likely lies in the 4-30 eV energy range, with a general shift towards higher frequencies for increasing X-ray brightness. However, the UV and X-ray spectra do not connect smoothly. Possible interpretations include: i) orientation effects, ii) additional absorption, iii) multiple emission components, and iv) a peculiar energy distribution of relativistic electrons. We discuss the first possibility in terms of an inhomogeneous helical jet model.
Using moderate-resolution optical spectra from 58 background Lyman-break galaxies and quasars at $z\sim 2.3-3$ within a $11.5'\times13.5'$ area of the COSMOS field ($\sim 1200\,\mathrm{deg}^2$ projected area density or $\sim 2.4\,h^{-1}\,\mathrm{Mpc}$ mean transverse separation), we reconstruct a 3D tomographic map of the foreground Ly$\alpha$ forest absorption at $2.2<z<2.5$ with an effective smoothing scale of $\sigma_{3d}\approx3.5\,h^{-1}\,\mathrm{Mpc}$ comoving. Comparing with 61 coeval galaxies with spectroscopic redshifts in the same volume, we find that the galaxy positions are clearly biased towards regions with enhanced IGM absorption in the tomographic map. We find an extended IGM overdensity with deep absorption troughs at $z=2.45$ associated with a recently-discovered galaxy protocluster at the same redshift. Based on simulations matched to our data, we estimate the enclosed dark matter mass within this IGM overdensity to be $M_{\rm dm} (z=2.45) = (9\pm4)\times 10^{13}\,h^{-1}\,\mathrm{M_\odot}$, and argue based on this mass and absorption strength that it will form at least one $z\sim0$ galaxy cluster with $M(z=0) = (3\pm 2) \times 10^{14}\,h^{-1}\mathrm{M_\odot}$, although its elongated nature suggests that it will likely collapse into two separate clusters. We also point out a compact overdensity of six MOSDEF galaxies at $z=2.30$ within a $r\sim 1\,h^{-1}\,\mathrm{Mpc}$ radius and $\Delta z\sim 0.006$, which does not appear to have a large associated IGM overdensity. These results demonstrate the potential of Ly$\alpha$ forest tomography on larger volumes to study galaxy properties as a function of environment, as well as revealing the large-scale IGM overdensities associated with protoclusters and other features of large-scale structure.
Up to ages of ~100 Myr, massive clusters are still swamped in large amounts of gas and dust, with considerable and uneven levels of extinction. At the same time, large grains (ices?) produced by type II supernovae profoundly alter the interstellar medium (ISM), thus resulting in extinction properties very different from those of the diffuse ISM. To obtain physically meaningful parameters of stars, from basic luminosities and effective temperatures to masses and ages, we must understand and measure the local extinction law. This problem affects all the massive young clusters discussed in his volume. We have developed a powerful method to unambiguously determine the extinction law in an uniform way across a cluster field, using multi-band photometry of red giant stars belonging to the red clump (RC). In the Large Magellanic Cloud, with about 20 RC stars per arcmin^2, we can easily derive a solid and self-consistent absolute extinction curve over the entire wavelength range of the photometry. Here, we present the extinction law of the Tarantula nebula (30 Dor) based on thousands of stars observed as part of the Hubble Tarantula Treasury Project.
Links to: arXiv, form interface, find, astro-ph, recent, 1509, contact, help (Access key information)
We summarise some of the recent progress in understanding the formation and evolution of globular clusters (GCs) in the context of galaxy formation and evolution. It is discussed that an end-to-end model for GC formation and evolution should capture four different phases: (1) star and cluster formation in the high-pressure interstellar medium of high-redshift galaxies, (2) cluster disruption by tidal shocks in the gas-rich host galaxy disc, (3) cluster migration into the galaxy halo, and (4) the final evaporation-dominated evolution of GCs until the present day. Previous models have mainly focussed on phase 4. We present and discuss a simple model that includes each of these four steps - its key difference with respect to previous work is the simultaneous addition of the high-redshift formation and early evolution of young GCs, as well as their migration into galaxy haloes. The new model provides an excellent match to the observed GC mass spectrum and specific frequency, as well as the relations of GCs to the host dark matter halo mass and supermassive black hole mass. These results show (1) that the properties of present-day GCs are reproduced by assuming that they are the natural outcome of regular high-redshift star formation (i.e. they form according to same physical processes that govern massive cluster formation in the local Universe), and (2) that models only including GC evaporation strongly underestimate their integrated mass loss over a Hubble time.
Mergers of gas-rich galaxies are key events in the hierarchical built-up of cosmic structures, and can lead to the formation of massive black hole binaries. By means of high-resolution hydrodynamical simulations we consider the late stages of a gas-rich major merger, detailing the dynamics of two circumnuclear discs, and of the hosted massive black holes during their pairing phase. During the merger gas clumps with masses of a fraction of the black hole mass form because of fragmentation. Such high-density gas is very effective in forming stars, and the most massive clumps can substantially perturb the black hole orbits. After $\sim 10$ Myr from the start of the merger a gravitationally bound black hole binary forms at a separation of a few parsecs, and soon after, the separation falls below our resolution limit of $0.39$ pc. At the time of binary formation the original discs are almost completely disrupted because of SNa feedback, while on pc scales the residual gas settles in a circumbinary disc with mass $\sim 10^5 M_\odot$. We also test that binary dynamics is robust against the details of the SNa feedback employed in the simulations, while gas dynamics is not. We finally highlight the importance of the SNa time-scale on our results.
We present a wide area (~ 8 x 8 kpc), sensitive map of CO (2-1) emission around the nearby starburst galaxy M82. Molecular gas extends far beyond the stellar disk, including emission associated with the well-known outflow as far as 3 kpc from M82's midplane. Kinematic signatures of the outflow are visible in both the CO and HI emission: both tracers show a minor axis velocity gradient and together they show double peaked profiles, consistent with a hot outflow bounded by a cone made of a mix of atomic and molecular gas. Combining our CO and HI data with observations of the dust continuum, we study the changing properties of the cold outflow as it leaves the disk. While H_2 dominates the ISM near the disk, the dominant phase of the cool medium changes as it leaves the galaxy and becomes mostly atomic after about a kpc. Several arguments suggest that regardless of phase, the mass in the cold outflow does not make it far from the disk; the mass flux through surfaces above the disk appears to decline with a projected scale length of ~ 1-2 kpc. The cool material must also end up distributed over a much wider angle than the hot outflow based on the nearly circular isophotes of dust and CO at low intensity and the declining rotation velocities as a function of height from the plane. The minor axis of M82 appears so striking at many wavelengths because the interface between the hot wind cavity and the cool gas produces Halpha, hot dust, PAH emission, and scattered UV light. We also show the level at which a face-on version of M82 would be detectable as an outflow based on unresolved spectroscopy. Finally, we consider multiple constraints on the CO-to-H$_2$ conversion factor, which must change across the galaxy but appears to be only a factor of ~ 2 lower than the Galactic value in the outflow.
What main mechanisms set the star formation rate (SFR) of galaxies? This PhD thesis is a quest into the influences of gas and Active Galactic Nuclei (AGNs) on the SFR, with particular focus on massive galaxies at z~2. First, a new code if presented; SImulator of GAlaxy Millimeter/submillimeter Emission (S\'IGAME) which can predict the atomic/molecular line emission from galaxies. By post-processing the outputs of cosmological simulations of galaxy formation with sub-grid physics recipes, S\'IGAME divides the Interstellar Medium (ISM) into different gas phases and derives density and temperature structure, with locally resolved radiation fields. This method is used to predict the strengths of CO rotational transitions as well as the [CII] emission line in normal star-forming galaxies at z~2. A CO ladder close to that of our own Galaxy is found, but with CO-H2 conversion factors about 3 times smaller. For a set of 7 simulated galaxies at z~2, the relation between [CII] luminosity and SFR displays a slope significantly steeper than that found for observed galaxies at z<0.5. A corresponding relation on kpc-scales is established for the first time theoretically. Finally, a separate study entails the analysis of CHANDRA CDF-S X-ray data with the aim of uncovering AGNs among massive galaxies at z~2. It is found that about every fifth massive galaxy, quenched or not, contains an X-ray luminous AGN. Interestingly, an even higher fraction of low-luminosity AGNs emerges in the X-ray undetected galaxies when performing a stacking analysis, and preferentially in the quenched ones, lending support to the importance of AGNs in impeding star formation during galaxy evolution.
Observations show the presence, in the halo of our Galaxy, of stars moving at velocities so high to require an acceleration mechanism involving the presence of a massive central black hole. Thus, in the frame of a galaxy hosting a supermassive black hole ($10^8$ $M_{\odot}$) we investigated a mechanism for the production of high velocity stars, which was suggested by the results of N-body simulations of the close interaction between a massive, orbitally decayed, globular cluster and the super massive black hole. The high velocity acquired by some stars of the cluster comes from the transfer of gravitational binding energy into kinetic energy of the escaping star originally orbiting around the cluster. After the close interaction with the massive black hole, stars could reach a velocity sufficient to travel in the halo and even overcome the galactic gravitational well, while some of them are just stripped from the globular cluster and start orbiting on precessing loops around the galactic centre.
The intermittent dissipation of interstellar turbulence is an important energy source in the diffuse ISM. Though on average smaller than the heating rates due to cosmic rays and the photoelectric effect on dust grains, the turbulent cascade can channel large amounts of energy into a relatively small fraction of the gas that consequently undergoes significant heating and chemical enrichment. In particular, this mechanism has been proposed as a solution to the long-standing problem of the high abundance of CH+ along diffuse molecular sight lines, which steady-state, low temperature models under-produce by over an order of magnitude. While much work has been done on the structure and chemistry of these small-scale dissipation zones, comparatively little attention has been paid to relating these zones to the properties of the large-scale turbulence. In this paper, we attempt to bridge this gap by estimating the temperature and CH+ column density along diffuse molecular sight-lines by post-processing 3-dimensional MHD turbulence simulations. Assuming reasonable values for the cloud density (30 / cm^3), size (20 pc), and velocity dispersion (2.3 km / s), we find that our computed abundances compare well with CH+ column density observations, as well as with observations of emission lines from rotationally excited H2 molecules.
Observations using the Fermi Large Area Telescope (Fermi-LAT) have found a significant gamma-ray excess surrounding the center of the Milky Way (GC). One possible interpretation of this excess invokes gamma-ray emission from an undiscovered population of either young or recycled pulsars densely clustered throughout the inner kiloparsec of the Milky Way. While these systems, by construction, have individual fluxes that lie below the point source sensitivity of the Fermi-LAT, they may already be observed in multiwavelength observations. Notably the Australia Telescope National Facility (ATNF) catalog of radio pulsars includes 270 sources observed in the inner 10 degrees around the GC. We calculate the gamma-ray emission observed from these 270 sources and obtain three key results: (1) point source searches in the GC region produce a plethora of highly significant gamma-ray "hotspots", compared to searches far from the Galactic plane, (2) there is no statistical correlation between the positions of these gamma-ray hotspots and the locations of ATNF pulsars, and (3) the spectrum of the most statistically significant gamma-ray hotspots is substantially softer than the spectrum of the GC gamma-ray excess. These results place strong constraints on models where young pulsars produce the majority of the gamma-ray excess, and disfavor some models where millisecond pulsars produce the gamma-ray excess.
Small-scale dark matter structure within the Milky Way is expected to affect pulsar timing. The change in gravitational potential induced by a dark matter halo passing near the line of sight to a pulsar would produce a varying delay in the light travel time of photons from the pulsar. Individual transits produce an effect that would either be too rare or too weak to be detected in 30-year pulsar observations. However, a population of dark matter subhalos would be expected to produce a detectable effect on the measured properties of pulsars if the subhalos constitute a significant fraction of the total halo mass. The effect is to increase the dispersion of measured period derivatives across the pulsar population. By statistical analysis of the ATNF pulsar catalogue, we place an upper limit on this dispersion of $\log \sigma_{\dot{P}} \leq -17.05$. We use this to place strong upper limits on the number density of ultracompact minihalos within the Milky Way. These limits are completely independent of the particle nature of dark matter.
The total population of Wolf-Rayet (WR) stars in the Galaxy is predicted by models to be as many as $\sim$6000 stars, and yet the number of catalogued WR stars as a result of optical surveys was far lower than this ($\sim$200) at the turn of this century. When beginning our WR searches using infrared techniques it was not clear whether WR number predictions were too optimistic or whether there was more hidden behind interstellar and circumstellar extinction. During the last decade we pioneered a technique of exploiting the near- and mid-infrared continuum colours for individual point sources provided by large-format surveys of the Galaxy, including 2MASS and Spitzer/GLIMPSE, to pierce through the dust and reveal newly discovered WR stars throughout the Galactic Plane. The key item to the colour discrimination is via the characteristic infrared spectral index produced by the strong winds of the WR stars, combined with dust extinction, which place WR stars in a relatively depopulated area of infrared colour-colour diagrams. The use of the Spitzer/GLIMPSE 8$\mu$m and, more recently, WISE 22$\mu$m fluxes together with cross-referencing with X-ray measurements in selected Galactic regions have enabled improved candidate lists that increased our confirmation success rate, achieved via follow-up infrared and optical spectroscopy. To date a total of 102 new WR stars have been found with many more candidates still available for follow-up. This constitutes an addition of $\sim$16\% to the current inventory of 642 Galactic WR stars. In this talk we review our methods and provide some new results and a preliminary review of their stellar and interstellar medium environments. We provide a roadmap for the future of this search, including statistical modeling, and what we can add to star formation and high mass star evolution studies.
We report on annual parallax and proper motion observations of H2O masers in S235AB-MIR, which is a massive young stellar object in the Perseus Arm. Using multi-epoch VLBI astrometry we measured a parallax of pi = 0.63 +- 0.03 mas, corresponding to a trigonometric distance of D = 1.56+-0.09 kpc, and source proper motion of ( u alpha cos d , u d) = (0.79 +- 0.12, -2.41 +- 0.14) mas/yr. Water masers trace a jet of diameter 15 au which exhibits a definite radial velocity gradient perpendicular to its axis. 3D maser kinematics were well modelled by a rotating cylinder with physical parameters: v_out = 45+-2 km/s, v_rot = 22+-3 km/s, i = 12+-2 degrees, which are the outflow velocity, tangential rotation velocity and line-of-sight inclination, respectively. One maser feature exhibited steady acceleration which may be related to the jet rotation. During our 15 month VLBI programme there were three `maser burst' events caught `in the act' which were caused by the overlapping of masers along the line of sight.
Links to: arXiv, form interface, find, astro-ph, recent, 1509, contact, help (Access key information)