We analyse the output of the hi-res cosmological zoom-in simulation ErisBH to study self-consistently the formation of a strong stellar bar in a Milky Way-type galaxy and its effect on the galactic structure, on the central gas distribution and on star formation. The simulation includes radiative cooling, star formation, SN feedback and a central massive black hole which is undergoing gas accretion and is heating the surroundings via thermal AGN feedback. A large central region in the ErisBH disk becomes bar-unstable after z~1.4, but a clear bar-like structure starts to grow significantly only after z~0.4, possibly triggered by the interaction with a massive satellite. At z~0.1 the bar reaches its maximum radial extent of l~2.2 kpc. As the bar grows, it becomes prone to buckling instability, which we quantify based on the anisotropy of the stellar velocity dispersion. The actual buckling event is observable at z~0.1, resulting in the formation of a boxy-peanut bulge clearly discernible in the edge-on view of the galaxy at z=0. The bar in ErisBH does not dissolve during the formation of the bulge but remains strongly non-axisymmetric down to the resolution limit of ~100 pc at z=0. During its early growth, the bar exerts a strong torque on the gas within its extent and drives gas inflows that enhance the nuclear star formation on sub-kpc scales. Later on the infalling gas is nearly all consumed into stars and, to a lesser extent, accreted onto the central black hole, leaving behind a gas-depleted region within the central ~2 kpc. Observations would more likely identify a prominent, large-scale bar at the stage when the galactic central region has already been quenched. Bar-driven quenching may play an important role in disk-dominated galaxies at all redshift. [Abridged]
We present coarse but robust star formation histories (SFHs) derived from spectro-photometric data of the Carnegie-Spitzer-IMACS Survey, for 22,494 galaxies at 0.3<z<0.9 with stellar masses of 10^9 Msun to 10^12 Msun. Our study moves beyond "average" SFHs and distribution functions of specific star formation rates (sSFRs) to individually measured SFHs for tens of thousands of galaxies. By comparing star formation rates (SFRs) with timescales of 10^10, 10^9, and 10^8 years, we find a wide diversity of SFHs: 'old galaxies' that formed most or all of their stars early; galaxies that formed stars with declining or constant SFRs over a Hubble time, and genuinely 'young galaxies' that formed most of their stars since z=1. This sequence is one of decreasing stellar mass, but, remarkably, each type is found over a mass range of a factor of 10. Conversely, galaxies at any given mass follow a wide range of SFHs, leading us to conclude that: (1) halo mass does not uniquely determine SFHs; (2) there is no 'typical' evolutionary track; and (3) "abundance matching" has limitations as a tool for inferring physics. Our observations imply that SFHs are set at an early epoch, and that--for most galaxies--the decline and cessation of star formation occurs over a Hubble-time, without distinct "quenching" events. SFH diversity is inconsistent with models where galaxy mass, at any given epoch, grows simply along relations between SFR and stellar mass, but is consistent with a 2-parameter lognormal form, lending credence to this model from a new and independent perspective.
We use the SPARC (Spitzer Photometry & Accurate Rotation Curves) database to study the relation between the central surface density of stars Sstar and dynamical mass Sdyn in 135 disk galaxies (S0 to dIrr). We find that Sdyn correlates tightly with Sstar over 4 dex. This central density relation can be described by a double power law. High surface brightness galaxies are consistent with a 1:1 relation, suggesting that they are self-gravitating and baryon dominated in the inner parts. Low surface brightness galaxies systematically deviate from the 1:1 line, indicating that the dark matter contribution progressively increases but remains tightly coupled to the stellar one. The observed scatter is small (~0.2 dex) and largely driven by observational uncertainties. The residuals show no correlations with other galaxy properties like stellar mass, size, or gas fraction.
We present observations with the Karl G. Jansky Very Large Array (VLA) at 3 GHz toward a sub-field of the XXL-North 25 deg$^{2}$ field targeting the first supercluster discovered in the XXL Survey. The structure has been found at a spectroscopic redshift of 0.43 and extending over $0.35\times0.1$ deg$^{2}$ on the sky. We present the 3 GHz VLA radio continuum observations, the radio mosaic and radio source catalogue, and, perform an analysis of the supercluster in the optical and radio regimes using photometric redshifts from the CFHTLS survey and our new VLA data. Our 3 GHz radio mosaic has a resolution of $3.2\times1.9$ arcsec$^2$, and covers an area of $41\times41$ arcmin$^2$ with rms noise level below $\sim20 \mu$Jy beam$^{-1}$ ($\sim11\mu$Jy beam$^{-1}$ in central $15\times15$ arcmin$^2$). We extract a catalogue of 155 radio sources $\geq6\sigma$, 8 of which are multicomponent sources, and 123 ($79\%$) can be associated with optical sources in the CFHTLS W1 catalogue. Applying Voronoi tessellation analysis (VTA) in the area around the X-ray identified supercluster we identify a total of 17 overdensities at $z=0.35-0.50$, 7 of which are associated with clusters detected in the $XMM-Newton$ XXL data. We find a mean photometric redshift of 0.43 for our overdensities, consistent with the spectroscopic redshifts of 7 X-ray detected clusters. The full VTA-identified structure extends over $\sim0.6\times0.2$ deg$^2$ on the sky, ($\sim12\times4$ Mpc$^{2}$ at $z=0.43$). We associate 8 radio sources with potential group/cluster member galaxies none of which are large radio galaxies. The spatial distribution of the red and blue potential group member galaxies suggests that the clusters are dynamically young and not virialised yet, as expected for hierarchical structure growth in a $\Lambda$CDM universe. Further spectroscopic data are required to analyze the dynamical state of the groups.
We provide evidence that UGC1382, long believed to be a passive elliptical galaxy, is actually a giant low surface brightness (GLSB) galaxy which rivals the archetypical GLSB Malin 1 in size. Like other GLSB galaxies, it has two components: a high surface brightness disk galaxy surrounded by an extended low surface brightness (LSB) disk. For UGC1382, the central component is a lenticular system with an effective radius of 6 kpc. Beyond this, the LSB disk has an effective radius of ~38 kpc and an extrapolated central surface brightness of ~26 mag/arcsec^2. Both components have a combined stellar mass of ~8x10^10 M_sun, and are embedded in a massive (10^10 M_sun) low-density (<3 M_sun/pc^2) HI disk with a radius of 110 kpc, making this one of the largest isolated disk galaxies known. The system resides in a massive dark matter halo of at least 2x10^12 M_sun. Although possibly part of a small group, its low density environment likely plays a role in the formation and retention of the giant LSB and HI disks. We model the spectral energy distributions and find that the LSB disk is likely older than the lenticular component. UGC1382 has UV-optical colors typical of galaxies transitioning through the green valley. Within the LSB disk are spiral arms forming stars at extremely low efficiencies. The gas depletion time scale of ~10^11 yr suggests that UGC1382 may be a very long term resident of the green valley. We find that the formation and evolution of the LSB disk is best explained by the accretion of gas-rich LSB dwarf galaxies.
We present a novel implementation of supermassive black hole (SMBH) formation, dynamics, and accretion in the massively parallel tree+SPH code, ChaNGa. This approach improves the modeling of SMBHs in fully cosmological simulations, allowing for a more detailed analysis of SMBH-galaxy co-evolution throughout cosmic time. Our scheme includes novel, physically motivated models for SMBH formation, dynamics and sinking timescales within galaxies, and SMBH accretion of rotationally supported gas. The sub-grid parameters that regulate star formation (SF) and feedback from SMBHs and SNe are optimized against a comprehensive set of z = 0 galaxy scaling relations using a novel, multi-dimensional parameter search. We have incorporated our new SMBH implementation and parameter optimization onto a new set of high resolution, large-scale cosmological simulations called Romulus. We present initial results from our flagship simulation, Romulus25, showing that our SMBH model results in SF efficiency, SMBH masses, and global cosmic SF and SMBH accretion histories that are consistent with observations. We discuss the importance of SMBH physics in shaping the evolution of massive galaxies and show how SMBH feedback is much more effective at regulating star formation compared to SNe feedback in this regime. Further, we show how each aspect of our SMBH model impacts this evolution compared to more common approaches. Finally, we present a science application of this scheme studying the properties and time evolution of an example dual AGN system, highlighting how our approach allows simulations to better study galaxy interactions and SMBH mergers in the context of galaxy-BH co-evolution.
The stellar halos of large galaxies represent a vital probe of the processes of galaxy evolution. They are the remnants of the initial bouts of star formation during the collapse of the proto-galactic cloud, coupled with imprint of ancient and on-going accretion events. Previously, we have reported the tentative detection of a possible, faint, extended stellar halo in the Local Group spiral, the Triangulum Galaxy (M33). However, the presence of substructure surrounding M33 made interpretation of this feature difficult. Here, we employ the final data set from the Pan-Andromeda Archaeological Survey (PAndAS), combined with an improved calibration and a newly derived contamination model for the region to revisit this claim. With an array of new fitting algorithms, fully accounting for contamination and the substantial substructure beyond the prominent stellar disk in M33, we reanalyse the surrounds to separate the signal of the stellar halo and the outer halo substructure. Using more robust search algorithms, we do not detect a large scale smooth stellar halo and place a limit on the maximum surface brightness of such a feature of ${\mu}_V$ = 35.5 mags per square arcsec, or a total halo luminosity of $L < 10^6L_{\odot}$.
We provide an analytic framework for interpreting observations of multiphase circumgalactic gas that is heavily informed by recent numerical simulations of thermal instability and precipitation in cool-core galaxy clusters. We start by considering the local conditions required for the formation of multiphase gas via two different modes: (1) uplift of ambient gas by galactic outflows, and (2) condensation in a stratified stationary medium in which thermal balance is explicitly maintained. Analytic exploration of these two modes provides insights into the relationships between the local ratio of the cooling and freefall time scales (i.e., t_cool / t_ff), the large-scale gradient of specific entropy, and development of precipitation and multiphase media in circumgalactic gas. We then use these analytic findings to interpret recent simulations of circumgalactic gas in which global thermal balance is maintained. We show that long-lasting configurations of gas with 5 < t_cool / t_ff < 20 and radial entropy profiles similar to observations of local cool-core galaxy cluster cores are a natural outcome of precipitation-regulated feedback. We conclude with some observational predictions that follow from these models. This work focuses primarily on precipitation and AGN feedback in galaxy cluster cores, because that is where the observations of multiphase gas around galaxies are most complete. However, many of the physical principles that govern condensation in those environments apply to circumgalactic gas around galaxies of all masses.
We present emission line templates for passively evolving ("retired") galaxies, useful for investigation of the evolution of the ISM in these galaxies, and characterization of their high-temperature source populations. The templates are based on high signal-to-noise ($>800$) co-added spectra ($3700-6800$\AA) of $\sim11500$ gas-rich Sloan Digital Sky Survey galaxies devoid of star-formation and active galactic nuclei. Stacked spectra are provided for the entire sample and sub-samples binned by mean stellar age. In Johansson~et al (2014), these spectra provided the first measurements of the He II 4686\AA\ line in passively-evolving galaxies, and the observed He II/H$\beta$ ratio constrained the contribution of accreting white dwarfs (the "single-degenerate" scenario) to the type Ia supernova rate. In this paper, the full range of unambiguously detected emission lines are presented. Comparison of the observed [O I] 6300\AA/H$\alpha$ ratio with photoionization models further constrains any high-temperature single-degenerate scenario for type Ia supernovae (with 1.5 $\lesssim$ T/$10^{5}K$ $\lesssim$ 10) to $\lesssim$3-6\% of the observed rate in the youngest age bin (i.e. highest SN Ia rate). Hence, for the same temperatures, in the presence of an ambient population of post-AGB stars, we exclude additional high-temperature sources with a combined ionizing luminosity of $\approx 1.35\times 10^{30} L_{\odot}/M_{\odot,*}$ for stellar populations with mean ages of 1 - 4 Gyrs. Furthermore, we investigate the extinction affecting both the stellar and nebular continuum. The latter shows about five times higher values. This contradicts isotropically distributed dust and gas that renders similar extinction values for both cases.
We present Submillimeter Array (SMA) 1.35 mm subarcsecond angular resolution observations toward the LkH{\alpha} 234 intermediate-mass star-forming region. The dust emission arises from a filamentary structure of $\sim$5 arcsec ($\sim$4500 au) enclosing VLA 1-3 and MM 1, perpendicular to the different outflows detected in the region. The most evolved objects are located at the southeastern edge of the dust filamentary structure and the youngest ones at the northeastern edge. The circumstellar structures around VLA 1, VLA 3, and MM 1 have radii between $\sim$200 and $\sim$375 au and masses in the $\sim$0.08-0.3 M$_{\odot}$ range. The 1.35 mm emission of VLA 2 arises from an unresolved (r$< 135$ au) circumstellar disk with a mass of $\sim$0.02 M$_{\odot}$. This source is powering a compact ($\sim$4000 au), low radial velocity ($\sim$7 km s$^{-1}$) SiO bipolar outflow, close to the plane of the sky. We conclude that this outflow is the "large-scale" counterpart of the short-lived, episodic, bipolar outflow observed through H$_2$O masers at much smaller scales ($\sim $180 au), and that has been created by the accumulation of the ejection of several episodic collimated events of material. The circumstellar gas around VLA 2 and VLA 3 is hot ($\sim$130 K) and exhibits velocity gradients that could trace rotation. There is a bridge of warm and dense molecular gas connecting VLA 2 and VLA 3. We discuss the possibility that this bridge could trace a stream of gas between VLA 3 and VLA 2, increasing the accretion rate onto VLA 2 to explain why this source has an important outflow activity.
In this paper, we present results of a photometric and kinematic study for a sample of 13 edge-on spiral galaxies with pronounced integral-shape warps of their stellar discs. The global structure of the galaxies is analyzed on the basis of the Sloan Digital Sky Survey (SDSS) imaging, in the g, r and i passbands. Spectroscopic observations are obtained with the 6-m Special Astrophysical Observatory telescope. In general, galaxies of the sample are typical bright spiral galaxies satisfying the Tully-Fisher relation. Most of the galaxies reside in dense spatial environments and, therefore, tidal encounters are the most probable mechanism for generating their stellar warps. We carried out a detailed analysis of the galaxies and their warps and obtained the following main results: (i) maximum angles of stellar warps in our sample are about 20{\deg}; (ii) warps start, on average, between 2 and 3 exponential scale lengths of a disc; (iii) stronger warps start closer to the center, weak warps start farther; (iv) warps are asymmetric, with the typical degree of asymmetry of about several degrees (warp angle); (v) massive dark halo is likely to preclude the formation of strong and asymmetric warps.
We present observations of the 85.69 GHz continuum emission and H42alpha line emission from the central 30 arcsec within NGC 4945. Both sources of emission originate from nearly identical structures that can be modelled as exponential discs with a scale length of ~2.1 arcsec (or ~40 pc). An analysis of the spectral energy distribution based on combining these data with archival data imply that 84% +/- 10% of the 85.69 GHz continuum emission originates from free-free emission. The electron temperature is 5400 +/- 600 K, which is comparable to what has been measured near the centre of the Milky Way Galaxy. The star formation rate (SFR) based on the H42alpha and 85.69 GHz free-free emission (and using a distance of 3.8 Mpc) is 4.35 +/- 0.25 M/yr. This is consistent with the SFR from the total infrared flux and with previous measurements based on recombination line emission, and it is within a factor of ~2 of SFRs derived from radio data. The Spitzer Space Telescope 24 micron data and Wide-field Infrared Survey Explorer 22 micron data yield SFRs ~10x lower than the ALMA measurements, most likely because the mid-infrared data are strongly affected by dust attenuation equivalent to A_V=150. These results indicate that SFRs based on mid-infrared emission may be highly inaccurate for dusty, compact circumnuclear starbursts.
We present a multi-wavelength analysis of five submillimeter sources (S_1.1mm = 0.54-2.02 mJy) that were detected during our 1.1-mm-deep continuum survey in the SXDF-UDS-CANDELS field (2 arcmin^2, 1sigma = 0.055 mJy beam^-1) using the Atacama Large Millimeter/submillimeter Array (ALMA). The two brightest sources correspond to a known single-dish (AzTEC) selected bright submillimeter galaxy (SMG), whereas the remaining three are faint SMGs newly uncovered by ALMA. If we exclude the two brightest sources, the contribution of the ALMA-detected faint SMGs to the infrared extragalactic background light is estimated to be ~ 4.1^{+5.4}_{-3.0} Jy deg^{-2}, which corresponds to ~ 16^{+22}_{-12}% of the infrared extragalactic background light. This suggests that their contribution to the infrared extragalactic background light is as large as that of bright SMGs. We identified multi-wavelength counterparts of the five ALMA sources. One of the sources (SXDF-ALMA3) is extremely faint in the optical to near-infrared region despite its infrared luminosity (L_IR ~ 1e12 L_sun or SFR ~ 100 M_sun yr^{-1}). By fitting the spectral energy distributions (SEDs) at the optical-to-near-infrared wavelengths of the remaining four ALMA sources, we obtained the photometric redshifts (z_photo) and stellar masses (M_*): z_photo ~ 1.3-2.5, M_* ~ (3.5-9.5)e10 M_sun. We also derived their star formation rates (SFRs) and specific SFRs (sSFRs) as ~ 30-200 M_sun yr^{-1} and ~ 0.8-2 Gyr^{-1}, respectively. These values imply that they are main-sequence star-forming galaxies.
We present a general method to identify infalling substructure in discrete datasets with position and line-of-sight velocity data. We exploit the fact that galaxies falling onto a brightest cluster galaxy (BCG) in a virialised cluster, or dwarf satellites falling onto a central galaxy like the Milky Way, follow nearly radial orbits. If the orbits are exactly radial, we show how to find the probability distribution for a satellite's energy, given a tracer density for the satellite population, by solving an Abel integral equation. This is an extension of Eddington (2016)'s classical formula for the isotropic distribution function. When applied to a system of galaxies, clustering in energy space can then be quantified using the Kullback-Leibler divergence, and groups of objects can be identified which, though separated in the sky, may be falling in on the same orbit. This method is tested using mock data and applied to the satellite galaxy population around M87, the BCG in Virgo, and a number of associations are found which may represent infalling galaxy groups.
We present an updated constrained hyperbolic/parabolic divergence cleaning algorithm for smoothed particle magnetohydrodynamics (SPMHD) that remains conservative with wave cleaning speeds which vary in space and time. This is accomplished by evolving the quantity $\psi / c_h$ instead of $\psi$. Doing so allows each particle to carry an individual wave cleaning speed, $c_h$, that can evolve in time without needing an explicit prescription for how it should evolve, preventing circumstances which we demonstrate could lead to runaway energy growth related to variable wave cleaning speeds. This modification requires only a minor adjustment to the cleaning equations and is trivial to adopt in existing codes. Finally, we demonstrate that our constrained hyperbolic/parabolic divergence cleaning algorithm, run for a large number of iterations, can reduce the divergence of the field to an arbitrarily small value, achieving $\nabla \cdot B=0$ to machine precision.
We define Baryon Acoustic Oscillation (BAO) distances $\hat{d}_\alpha(z, z_c)$, $\hat{d}_z(z, z_c)$, and $\hat{d}_/(z, z_c)$ that do not depend on cosmological parameters. These BAO distances are measured as a function of redshift $z$ with the Sloan Digital Sky Survey (SDSS) data release DR12. From these BAO distances alone, or together with the correlation angle $\theta_\textrm{MC}$ of the Cosmic Microwave Background (CMB), we constrain the cosmological parameters in several scenarios. We find $4.3 \sigma$ tension between the BAO plus $\theta_\textrm{MC}$ data and a cosmology with flat space and constant dark energy density $\Omega_\textrm{DE}(a)$. Releasing one and/or the other of these constraints obtains agreement with the data. We measure $\Omega_\textrm{DE}(a)$ as a function of $a$.
We describe the redmonster automated redshift measurement and spectral classification software designed for the extended Baryon Oscillation Spectroscopic Survey (eBOSS) of the Sloan Digital Sky Survey IV (SDSS-IV). We describe the algorithms, the template standard and requirements, and the newly developed galaxy templates to be used on eBOSS spectra. We present results from testing on early data from eBOSS, where we have found a 90.5% automated redshift and spectral classification success rate for the luminous red galaxy sample (redshifts 0.6 $\lesssim$ $z$ $\lesssim$ 1.0). The \texttt{redmonster} performance meets the eBOSS cosmology requirements for redshift classification and catastrophic failures, and represents a significant improvement over the previous pipeline. We describe the empirical processes used to determine the optimum number of additive polynomial terms in our models and an acceptable $\Delta\chi_r^2$ threshold for declaring statistical confidence. Statistical errors on redshift measurement due to photon shot noise are assessed, and we find typical values of a few tens of km s$^{-1}$. An investigation of redshift differences in repeat observations scaled by error estimates yields a distribution with a Gaussian mean and standard deviation of $\mu\sim$ 0.01 and $\sigma\sim$ 0.65, respectively, suggesting the accuracy of the statistical redshift errors estimates are over-estimated by $\sim$ 54%. We assess the effects of object magnitude, signal-to-noise ratio, fiber number, and fiber head location on the pipeline's redshift success rate. Finally, we describe directions of ongoing development.
Links to: arXiv, form interface, find, astro-ph, recent, 1607, contact, help (Access key information)
We show how the interplay between active galactic nuclei (AGN) and merger history determines whether a galaxy quenches star formation at high redshift. We first simulate, in a full cosmological context, a galaxy of total dynamical mass $10^{12}\,M_{\odot}$ at $z=2$. Then we systematically alter the accretion history of the galaxy by minimally changing the linear overdensity in the initial conditions. This "genetic modification" approach allows the generation of three sets of $\Lambda$CDM initial conditions leading to maximum merger ratios of 1:10, 1:5 and 2:3 respectively. The changes leave the final halo mass, large scale structure and local environment unchanged, providing a controlled numerical experiment. Interaction between the AGN physics and mergers in the three cases lead respectively to a star-forming, temporarily-quenched and permanently-quenched galaxy. However the differences do not primarily lie in the black hole accretion rates, but in the kinetic effects of the merger: the galaxy is resilient against AGN feedback unless its gaseous disk is first disrupted. Typical accretion rates are comparable in the three cases, falling below $0.1\,M_{\odot}$ yr$^{-1}$, equivalent to around $2\%$ of the Eddington rate or $10^{-3}$ times the pre-quenching star formation rate, in agreement with observations. This low level of black hole accretion can be sustained even when there is insufficient dense cold gas for star formation. Conversely, supernova feedback is too distributed to generate outflows in high-mass systems, and cannot maintain quenching over periods longer than the halo gas cooling time.
We present a coherent database of spectroscopic observations of far-IR fine-structure lines from the Herschel/PACS archive for a sample of 170 local AGN, plus a comparison sample of 20 starburst galaxies and 43 dwarf galaxies. Published Spitzer/IRS and Herschel/SPIRE line fluxes are included to extend our database to the full 10-600 $\mu m$ spectral range. The observations are compared to a set of CLOUDY photoionisation models to estimate the above physical quantities through different diagnostic diagrams. We confirm the presence of a stratification of gas density in the emission regions of the galaxies, which increases with the ionisation potential of the emission lines. The new [OIV]25.9$\mu m$/[OIII]88$\mu m$ vs [NeIII]15.6$\mu m$/[NeII]12.8$\mu m$ diagram is proposed as the best diagnostic to separate: $i)$ AGN activity from any kind of star formation; and $ii)$ low-metallicity dwarf galaxies from starburst galaxies. Current stellar atmosphere models fail to reproduce the observed [OIV]25.9$\mu m$/[OIII]88$\mu m$ ratios, which are much higher when compared to the predicted values. Finally, the ([NeIII]15.6$\mu m$ + [NeII]12.8$\mu m$)/([SIV]10.5$\mu m$ + [SIII]18.7$\mu m$) ratio is proposed as a promising metallicity tracer to be used in obscured objects, where optical lines fail to accurately measure the metallicity. The diagnostic power of mid- to far-infrared spectroscopy shown here for local galaxies will be of crucial importance to study galaxy evolution during the dust-obscured phase at the peak of the star formation and black-hole accretion activity ($1 < z < 4$). This study will be addressed by future deep spectroscopic surveys with present and forthcoming facilities such as JWST, ALMA, and SPICA.
We report interferometric measurements of [NII] 205 um fine-structure line emission from a representative sample of three galaxies at z=5-6 using the Atacama Large (sub)Millimeter Array (ALMA). These galaxies were previously detected in [CII] and far-infrared continuum emission and span almost two orders of magnitude in star formation rate (SFR). Our results show at least two different regimes of ionized inter-stellar medium properties for galaxies in the first billion years of cosmic time, separated by their L_[CII]/L_[NII] ratio. We find extremely low [NII] emission compared to [CII] (L_ [CII]/L_[NII]=68 [+200/-28]) from a "typical" L*_UV star-forming galaxy, likely directly or indirectly (by its effect on the radiation field) related to low dust abundance and low metallicity. The infrared-luminous modestly star-forming Lyman Break Galaxy (LBG) in our sample is characterized by an ionized-gas fraction (L_[CII]/L_[NII]<=20) typical of local star-forming galaxies and shows evidence for spatial variations in its ionized-gas fraction across an extended gas reservoir. The extreme SFR, warm and compact dusty starburst AzTEC-3 shows an ionized fraction lower than expected given its star-formation rate surface density (L_[CII]/L_[NII]=22+/-8) suggesting that [NII] dominantly traces a diffuse ionized medium rather than star-forming HII regions in this type of galaxy. This highest redshift sample of [NII] detections provides some of the first constraints on ionized and neutral gas modeling attempts and on the structure of the inter-stellar medium at z=5-6 in "normal" galaxies and starbursts.
Aims. Galaxy mass models based on simple and analytical functions for the density and potential pairs have been widely proposed in the literature. Disk models constrained by kinematic data alone give information on the global disk structure only very near the Galactic plane. We attempt to circumvent this issue by constructing disk mass models whose three-dimensional structures are constrained by a recent Galactic star counts model in the near-infrared and also by observations of the hydrogen distribution in the disk. Our main aim is to provide models for the gravitational potential of the Galaxy that are fully analytical but also with a more realistic description of the density distribution in the disk component. Methods. From the disk model directly based on the observations (here divided into the thin and thick stellar disks and the HI and H$_2$ disks subcomponents), we produce fitted mass models by combining three Miyamoto-Nagai disk profiles of any "model order" (1, 2, or 3) for each disk subcomponent. The Miyamoto-Nagai disks are combined with models for the bulge and "dark halo" components and the total set of parameters is adjusted by observational kinematic constraints. A model which includes a ring density structure in the disk, beyond the solar Galactic radius, is also investigated. Results. The Galactic mass models return very good matches to the imposed observational constraints. In particular, the model with the ring density structure provides a greater contribution of the disk to the rotational support inside the solar circle. The gravitational potential models and their associated force-fields are described in analytically closed forms, and in addition, they are also compatible with our best knowledge of the stellar and gas distributions in the disk component. The gravitational potential models are suited for investigations of orbits in the Galactic disk.
We present an analysis of the impact of bars on the radial distributions of the different types of supernovae (SNe) in the stellar discs of host galaxies with various morphologies. We find that in Sa-Sbc galaxies, the radial distribution of core-collapse (CC) SNe in barred hosts is inconsistent with that in unbarred ones, while the distributions of SNe Ia are not significantly different. At the same time, the radial distributions of both types of SNe in Sc-Sm galaxies are not affected by bars. We propose that the additional mechanism shaping the distributions of Type Ia and CC SNe can be explained within the framework of substantial suppression of massive star formation in the radial range swept by strong bars, particularly in early-type spirals. The radial distribution of CC SNe in unbarred Sa-Sbc galaxies is more centrally peaked and inconsistent with that in unbarred Sc-Sm hosts, while the distribution of SNe Ia in unbarred galaxies is not affected by host morphology. These results can be explained by the distinct distributions of massive stars in the discs of early-and late-type spirals.
We present new photometry of the dwarf irregular galaxy DDO 68, one of the most metal-poor and least massive dwarfs, located in the Lynx-Cancer Void. The images were acquired with the Large Binocular Telescope in the g and r passbands, and show unequivocally that DDO 68 has previously unknown stellar streams related to the accretion of at least two smaller companions: a flea with smaller fleas biting it, to put it in Jonathan Swift's words. Our data provide direct observational evidence of multiple merging occurring at very low galactic mass scales. We present the results of an N-body simulation of the interaction of three dwarf galaxies which reproduces well the main morphological features of DDO 68.
Most investigations of the X-ray variability of active galactic nuclei (AGN) have been concentrated on the detailed analyses of individual, nearby sources. A relatively small number of studies have treated the ensemble behaviour of the more general AGN population in wider regions of the luminosity-redshift plane. We want to determine the ensemble variability properties of a rich AGN sample, called Multi-Epoch XMM Serendipitous AGN Sample (MEXSAS), extracted from the latest release of the XMM-Newton Serendipitous Source Catalogue, with redshift between 0.1 and 5, and X-ray luminosities, in the 0.5-4.5 keV band, between 10^{42} and 10^{47} erg/s. We caution on the use of the normalised excess variance (NXS), noting that it may lead to underestimate variability if used improperly. We use the structure function (SF), updating our previous analysis for a smaller sample. We propose a correction to the NXS variability estimator, taking account of the light-curve duration in the rest-frame, on the basis of the knowledge of the variability behaviour gained by SF studies. We find an ensemble increase of the X-ray variability with the rest-frame time lag tau, given by tau^{0.12}. We confirm an inverse dependence on the X-ray luminosity, approximately as L_X^{-0.19}. We analyse the SF in different X-ray bands, finding a dependence of the variability on the frequency as nu^{-0.15}, corresponding to a softer when brighter trend. In turn, this dependence allows us to parametrically correct the variability estimated in observer-frame bands to that in the rest-frame, resulting in a moderate shift upwards (V-correction). Ensemble X-ray variability of AGNs is best described by the structure function. An improper use of the normalised excess variance may lead to an underestimate of the intrinsic variability, so that appropriate corrections to the data or the models must be applied to prevent these effects.
Due to their long mean free paths, X-rays are expected to have many significant impacts globally on the properties of the intergalactic medium (IGM) by their heating and ionizing processes on large scales. At high redshifts, X-rays from Population (Pop) III binaries might have important effects on cosmic reionization and the Lyman alpha forest. As a continuation of our previous work on Pop III binary X-rays (Xu et al. 2014), we use the Pop III distribution and evolution from the Renaissance Simulations, a suite of self-consistent cosmological radiation hydrodynamics simulations of the formation of the first galaxies, to calculate the X-ray luminosity density and background over the redshift range 20 > z > 7.6. As we find that Pop III star formation continues at a low, nearly constant rate to the end of reionization, X-rays are being continuously produced at significant rates compared to other possible X-ray sources, such as AGNs and normal X-ray binaries during the same period of time. We estimate that Pop III binaries produce approximately 6 eV of energy in the X-rays per hydrogen atom. We calculate the X-ray background for different monochromatic photon energies. KeV X-rays redshift and accumulate to produce a strong X-ray background spectrum extending to roughly 500 eV. The X-ray background is strong enough to heat the IGM to ~ 1000 K and to ionize a few percent of the neutral hydrogen. These effects are important for an understanding of the neutral hydrogen hyperfine transition 21-cm line signatures, the Ly alpha forest, and optical depth of the CMB to Thomson scattering.
Infrared-faint radio sources (IFRS) are a class of radio-loud (RL) active galactic nuclei (AGN) at high redshifts (z > 1.7) that are characterised by their relative infrared faintness, resulting in enormous radio-to-infrared flux density ratios of up to several thousand. We aim to test the hypothesis that IFRS are young AGN, particularly GHz peaked-spectrum (GPS) and compact steep-spectrum (CSS) sources that have a low frequency turnover. We use the rich radio data set available for the Australia Telescope Large Area Survey fields, covering the frequency range between 150 MHz and 34 GHz with up to 19 wavebands from different telescopes, and build radio spectral energy distributions (SEDs) for 34 IFRS. We then study the radio properties of this class of object with respect to turnover, spectral index, and behaviour towards higher frequencies. We also present the highest-frequency radio observations of an IFRS, observed with the Plateau de Bure Interferometer at 105 GHz, and model the multi-wavelength and radio-far-infrared SED of this source. We find IFRS usually follow single power laws down to observed frequencies of around 150 MHz. Mostly, the radio SEDs are steep, but we also find ultra-steep SEDs. In particular, IFRS show statistically significantly steeper radio SEDs than the broader RL AGN population. Our analysis reveals that the fractions of GPS and CSS sources in the population of IFRS are consistent with the fractions in the broader RL AGN population. We find that at least 18% of IFRS contain young AGN, although the fraction might be significantly higher as suggested by the steep SEDs and the compact morphology of IFRS. The detailed multi-wavelength SED modelling of one IFRS shows that it is different from ordinary AGN, although it is consistent with a composite starburst-AGN model with a star formation rate of 170 solar masses per year.
Almost 200 different species have been detected in the interstellar medium (ISM) during the last decades, revealing not only simple species but complex molecules with more than 6 atoms. Other exotic compounds, like the weakly-bound dimer (H2)2, have also been detected in astronomical sources like Jupiter. We aim at detecting for the first time the CO-H2 van der Waals complex in the ISM, which if detected can be a sensitive indicator for low temperatures. We use the IRAM30m telescope, located in Pico Veleta (Spain), to search for the CO-H2 complex in a cold, dense core in TMC-1C (with a temperature of 10 K). All the brightest CO-H2 transitions in the 3 mm (80-110 GHz) band have been observed with a spectral resolution of 0.5-0.7 km/s, reaching a rms noise level of 2 mK. The simultaneous observation of a broad frequency band, 16 GHz, has allowed us to conduct a serendipitous spectral line survey. No lines belonging to the CO-H2 complex have been detected. We have set up a new, more stringent upper limit for its abundance to be [CO-H2]/[CO] = 5x10^{-6}, while we expect the abundance of the complex to be in the range 10^{-8}-10^{-3}. The spectral line survey has allowed us to detect 75 lines associated with 41 different species (including isotopologues). We detect a number of complex organic species, e.g. methyl cyanide (CH3CN), methanol (CH3OH), propyne (CH3CCH) and ketene (CH2CO), associated with cold gas (excitation temperatures about 7 K), confirming the presence of these complex species not only in warm objects but also in cold regimes.
We analyse new observations with the International Low Frequency Array (LOFAR) telescope, and archival data from the Multi-Element Radio Linked Interferometer Network (MERLIN) and the Karl G. Jansky Very Large Array (VLA). We model the spatially resolved radio spectrum of Arp 220 from 150 MHz to 33 GHz. We present an image of Arp 220 at 150 MHz with resolution $0.65''\times0.35''$, sensitivity 0.15 mJy beam$^{-1}$, and integrated flux density $394\pm59$ mJy. More than 80% of the detected flux comes from extended ($6''\approx$2.2 kpc) steep spectrum ($\alpha=-0.7$) emission, likely from star formation in the molecular disk surrounding the two nuclei. We find elongated features extending $0.3''$ (110 pc) and $0.9''$ (330 pc) from the eastern and western nucleus respectively, which we interpret as evidence for outflows. The extent of radio emission requires acceleration of cosmic rays far outside the nuclei. We find that a simple three component model can explain most of the observed radio spectrum of the galaxy. When accounting for absorption at 1.4 GHz, Arp 220 follows the FIR/radio correlation with $q=2.36$, and we estimate a star formation rate of 220 M$_\odot\text{yr}^{-1}$. We derive thermal fractions at 1 GHz of less than 1% for the nuclei, which indicates that a major part of the UV-photons are absorbed by dust. International LOFAR observations shows great promise to detect steep spectrum outflows and probe regions of thermal absorption. However, in LIRGs the emission detected at 150 MHz does not necessarily come from the main regions of star formation. This implies that high spatial resolution is crucial for accurate estimates of star formation rates for such galaxies at 150 MHz.
Since their discovery, the Quintuplet proper members (QPMs) have been somewhat mysterious in nature. Originally dubbed the "cocoon stars" due to their cool featureless spectra, high-resolution near-infrared imaging observations have shown that at least two of the objects exhibit "pinwheel" nebulae consistent with binary systems with a carbon-rich Wolf-Rayet star and O/B companion. In this paper, we present 19.7, 25.2, 31.5, and 37.1 $\mu$m observations of the QPMs (with an angular resolution of 3.2-3.8") taken with the Faint Object Infrared Camera for the SOFIA Telescope (FORCAST) in conjunction with high-resolution ($\sim$0.1-0.2") images at 8.8 and 11.7 $\mu$m from the Thermal-Region Camera Spectrograph (TReCS). DUSTY models of the thermal dust emission of two of the four detected QPMs, Q2 and Q3, are fitted by radial density profiles which are consistent with constant mass loss rates ($\rho_d \propto r^{-2}$). For the two remaining sources, Q1 and Q9, extended structures ($\sim$ 1") are detected around these objects in high-resolution imaging data. Based on the fitted dust masses, Q9 has an unusually large dust reservoir ($\mathrm{M_d}=1.3^{+0.8}_{-0.4}\times 10^{-3} \mathrm{M_{\odot}}$) compared to typical dusty Wolf-Rayet stars which suggests that it may have recently undergone an episode of enhanced mass loss.
To investigate galaxy properties as a function of their total stellar mass, we obtained 21cm HI line observations at the 100-m class Nan\c{c}ay Radio Telescope of 2839 galaxies from the Sloan Digital Sky Survey (SDSS) in the Local Volume (900<cz<12,000 km/s), dubbed the Nan\c{c}ay Interstellar Baryons Legacy Extragalactic Survey (NIBLES) sample. They were selected evenly over their entire range of absolute SDSS z-band magnitudes (-13.5 to -24 mag), which were used as a proxy for their stellar masses. Here, a first, global presentation of the observations and basic results is given, their further analysis will be presented in other papers in this series. The galaxies were selected based on their properties, as listed in SDSS DR5. Comparing this photometry to their total HI masses, we noted that, for a few percent, the SDSS magnitudes appeared severely misunderestimated, as confirmed by our re-measurements for selected objects. Although using the later DR9 results eliminated this problem in most cases, 384 still required manual photometric source selection. Usable HI spectra were obtained for 2600 galaxies, of which 1733 (67%) were clearly detected and 174 (7%) marginally. The spectra for 241 other observed galaxies could not be used for further analysis owing to problems with either the HI or the SDSS data. We reached the target number of about 150 sources per half-magnitude bin over the Mz range -16.5 to -23 mag. Down to -21 mag the overall detection rate is rather constant at the ~75% level but it starts to decline steadily towards the 30% level at -23 mag. Making regression fits by comparing total HI and stellar masses for our sample, including our conservatively estimated HI upper limits for non-detections, we find the relationship log(M_HI/M*) = -0.59 log(M*) + 5.05, which lies significantly below the relationship found in the M_HI/M* - M* plane when only using HI detections.
We investigate the correlations in galaxy shapes between optical and radio wavelengths using archival observations of the COSMOS field. Cross-correlation studies between different wavebands will become increasingly important for precision cosmology as future large surveys may be dominated by systematic rather than statistical errors. In the case of weak lensing, galaxy shapes must be measured to extraordinary accuracy (shear systematics of $< 0.01\%$) in order to achieve good constraints on dark energy parameters. By using shape information from overlapping surveys in optical and radio bands, robustness to systematics may be significantly improved without loss of constraining power. Here we use HST-ACS optical data, VLA radio data, and extensive simulations to investigate both our ability to make precision measurements of source shapes from realistic radio data, and to constrain the intrinsic astrophysical scatter between the shapes of galaxies as measured in the optical and radio wavebands. By producing a new image from the VLA-COSMOS L-band radio visibility data that is well suited to galaxy shape measurements, we are able to extract precise measurements of galaxy position angles. Comparing to corresponding measurements from the HST optical image, we set a lower limit on the intrinsic astrophysical scatter in position angles, between the optical and radio bands, of $\sigma_\alpha > 0.212\pi$ radians (or $38.2^{\circ}$) at a $95\%$ confidence level.
The optical morphology of galaxies holds the cumulative record of their assembly history, and techniques for its quantitative characterization offer a promising avenue toward understanding galaxy formation and evolution. However, the morphology of star-forming galaxies is generally dictated by the youngest stellar component, which can readily overshine faint structural/morphological features in the older underlying stellar background (e.g., relics from recent minor mergers) that could hold important insights into the galaxy build-up process. Stripping off galaxy images from the emission from stellar populations younger than an adjustable age cutoff tcut can, therefore, provide a valuable tool in extragalactic research. RemoveYoung (RY), a publicly available tool that is presented here, exploits the combined power of integral field spectroscopy (IFS) and spectral population synthesis (SPS) toward this goal. Two-dimensional (2D) post-processing of SPS models to IFS data cubes with RY permits computation of the spectral energy, surface brightness, and stellar surface density distribution of stellar populations older than a user-defined tcut. This suggests a variety of applications of star-forming galaxies, such as interacting or merging galaxy pairs and lower mass starburst galaxies near and far; these include blue compact and tidal dwarf galaxies.
In this paper we present simulated observations of galaxies at z=2 and z=3 to probe the capabilities of next-generation telescopes (E-ELT and JWST) to measure the structural and photometrical properties of high redshift galaxies. We carry out an extensive set of simulations of high redshift galaxies adopting the specifications of the E-ELT first light instrument MICADO. The main parameters (sizes, Sersic index, magnitudes) of the galaxies are measured using GALFIT and the accuracy of the determinations is assessed by comparing the input values to the measurements from many runs with different statistical noise. We also address the effects on the measurements accuracy of possible spatial variation of the PSF in the field. We find that from 3h exposure E-ELT near-IR images of galaxies at z=2 and z=3 it will be possible to measure the size, the total magnitude and the galaxy morphology with an accuracy of 5% for objects as faint as H=25 and half light size of 0.2 arcsec. The effective radius of compact, early type galaxies will also be recovered with 5% accuracy, provided that their half light size exceeds 20 mas. These results are compared with those expected from simulated observations obtained with NIRCam on board of the JWST.
We compare common star-formation rate (SFR) indicators in the local Universe in the GAMA equatorial fields (around 160 sq. deg.), using ultraviolet (UV) photometry from GALEX, far-infrared (FIR) and sub-millimetre (sub-mm) photometry from H-ATLAS, and Halpha spectroscopy from the GAMA survey. With a high-quality sample of 745 galaxies (median redshift 0.08), we consider three SFR tracers: UV luminosity corrected for dust attenuation using the UV spectral slope beta (SFRUV,corr), Halpha line luminosity corrected for dust using the Balmer decrement (BD) (SFRHalpha,corr), and the combination of UV and IR emission (SFRUV+IR). We demonstrate that SFRUV,corr can be reconciled with the other two tracers after applying attenuation corrections by calibrating IRX (i.e. the IR to UV luminosity ratio) and attenuation in the Halpha (derived from BD) against beta. However, beta on its own is very unlikely to be a reliable attenuation indicator. We find that attenuation correction factors depend on parameters such as stellar mass, z and dust temperature (Tdust), but not on Halpha equivalent width (EW) or Sersic index. Due to the large scatter in the IRX vs beta correlation, when compared to SFRUV+IR, the beta-corrected SFRUV,corr exhibits systematic deviations as a function of IRX, BD and Tdust.
Aims. We aim to study the 250 micron luminosity function (LF) down to much
fainter luminosities than achieved by previous efforts.
Methods. We developed a modified stacking method to reconstruct the 250
micron LF using optically selected galaxies from the SDSS survey and Herschel
maps of the GAMA equatorial fields and Stripe 82. Our stacking method not only
recovers the mean 250 micron luminosities of galaxies that are too faint to be
individually detected, but also their underlying distribution functions.
Results. We find very good agreement with previous measurements in the
overlapping luminosity range. More importantly, we are able to derive the LF
down to much fainter luminosities (around 25 times fainter) than achieved by
previous studies. We find strong positive luminosity evolution \propto (1 +
z)^4.89\pm1.07 and moderate negative density evolution \propto (1 +
z)^-1.02\pm0.54 over the redshift range z=[0.02, 0.5].
Stellar evolution theory has been extraordinarily successful at explaining the different phases under which stars form, evolve and die. While the strongest constraints have traditionally come from binary stars, the advent of asteroseismology is bringing unique measures in well-characterised stars. For stellar populations in general, however, only photometric measures are usually available, and the comparison with the predictions of stellar evolution theory have mostly been qualitative. For instance, the geometrical shapes of isochrones have been used to infer ages of coeval populations, but without any proper statistical basis. In this chapter we provide a pedagogical review on a Bayesian formalism to make quantitative inferences on the properties of single, binary and small ensembles of stars, including unresolved populations. As an example, we show how stellar evolution theory can be used in a rigorous way as a prior information to measure the ages of stars between the ZAMS and the Helium flash, and their uncertainties, using photometric data only.
The characteristic and distinctive features of the visibility amplitude of interferometric observations for compact objects like stars in the immediate vicinity of the central black hole in our Galaxy are considered. These features are associated with the specifics of strong gravitational scattering of point sources by black holes, wormholes, or black_white holes. The revealed features will help to determine the most important topological characteristics of the central object in our Galaxy: whether this object possesses the properties of only a black hole or also has characteristics unique to wormholes or black_white holes. These studies can be used to interpret the results of optical, infrared, and radio interferometric observations.
The Tully-Fisher relation (TFR) links the stellar mass of a disk galaxy, $M_{\rm str}$, to its rotation speed: it is well approximated by a power law, shows little scatter, and evolves weakly with redshift. The relation has been interpreted as reflecting the mass-velocity scaling ($M\propto V^3$) of dark matter halos, but this interpretation has been called into question by abundance-matching (AM) models, which predict the galaxy-halo mass relation to be non-monotonic and rapidy evolving. We study the TFR of luminous spirals and its relation to AM using the EAGLE set of $\Lambda$CDM cosmological simulations. Matching both relations requires disk sizes to satisfy constraints given by the concentration of halos and their response to galaxy assembly. EAGLE galaxies approximately match these constraints and show a tight mass-velocity scaling that compares favourably with the observed TFR. The TFR is degenerate to changes in galaxy formation efficiency and the mass-size relation; simulations that fail to match the galaxy stellar mass function may fit the observed TFR if galaxies follow a different mass-size relation. The small scatter in the simulated TFR results because, at fixed halo mass, galaxy mass and rotation speed correlate strongly, scattering galaxies along the main relation. EAGLE galaxies evolve with lookback time following approximately the prescriptions of AM models and the observed mass-size relation of bright spirals, leading to a weak TFR evolution consistent with observation out to $z=1$. $\Lambda$CDM models that match both the abundance and size of galaxies as a function of stellar mass have no difficulty reproducing the observed TFR and its evolution.
In an intensive observational campaign in the 9 month duration of Chandra X-ray Visionary Project that was conducted in the year 2012, 39 large X-ray flares of Sgr A* were recorded. An analysis of the times of the observed flares reveals that the 39 flares are separated in time by intervals that are grouped around integer numbers times 0.10333 days. This time interval is thus the period of a uniform grid of equally spaced points on the time axis. The grouping of the flares around tic marks of this grid is derived from the data with at least a 3.2 {\sigma} level of statistical significance. No signal of any period can be found among 22 flares recorded by Chandra in the years 2013-2014. If the 0.10333 d period is that of a nearly circular Keplerian orbit around the blackhole at the center of the Galaxy, its radius is at 7.6 Schwarzschild radii. Large flares were more likely to be triggered when the agent responsible for their outbursts was near the peri-center phase of its slightly eccentric orbit.
Cosmic far-infrared background (CFIRB) is a powerful probe of the history of star formation rate and the connection between baryons and dark matter. In this work, we explore to which extent the CFIRB anisotropies can be reproduced by a simple physical framework for galaxy evolution, the gas regulator (bathtub) model. The model is based on continuity equations for gas, stars, and metals, taking into account cosmic gas accretion, star formation, and gas ejection. Our model not only provides a good fit to the CFIRB power spectra measured by Planck, but also agrees well with the correlation between CFIRB and gravitational lensing, far-infrared galaxy number counts, and bolometric infrared luminosity functions. The strong clustering of CFIRB indicates a large galaxy bias, which corresponds to haloes of mass 10^12.5 Msun at z=2; thus, CFIRB favors strong infrared emission in massive haloes, which is higher than the expectation from the star formation rate. We provide constraints and fitting functions for the cosmic star formation history and the infrared luminosity-halo mass relation.
Gamma-ray Bursts (GRBs) are the most powerful explosions known, capable of outshining the rest of gamma-ray sky during their short-lived prompt emission. Their cosmological nature makes them the best tool to explore the final stages in the lives of very massive stars up to the highest redshifts. Furthermore, studying the emission from their low-energy counterparts (optical and infrared) via rapid spectroscopy, we have been able to pin down the exact location of the most distant galaxies as well as placing stringent constraints on their host galaxies and intervening systems at low and high-redshift (e.g. metallicity and neutral hydrogen fraction). In fact, each GRB spectrum contains absorption features imprinted by metals in the host interstellar medium (ISM) as well as the intervening intergalactic medium (IGM) along the line of sight. In this chapter we summarize the progress made using a large dataset of GRB spectra in understanding the nature of both these absorbers and how GRBs can be used to study the early Universe, in particular to measure the neutral hydrogen fraction and the escape fraction of UV photons before and during the epoch of re-ionization.
Originally developed to image the shadow region of the central black hole in Sagittarius A* and in the nearby galaxy M87, the Event Horizon Telescope (EHT) provides deep, very high angular resolution data on other AGN sources too. The challenges of working with EHT data have spurred the development of new image reconstruction algorithms. This work briefly reviews the status of the EHT and its utility for observing AGN sources, with emphasis on novel imaging techniques that offer the promise of better reconstructions at 1.3 mm and other wavelengths.
Links to: arXiv, form interface, find, astro-ph, recent, 1607, contact, help (Access key information)
In this paper, we investigate the relationship between star formation and structure, using a mass-complete sample of 27,893 galaxies at 0.5<z<2.5 selected from 3D-HST. We confirm that, at fixed stellar mass, star-forming galaxies are larger than quiescent galaxies at all redshifts. However, within the star-forming population, there is no relation between star formation rate and size: when dividing star-forming galaxies into quartiles based on their residual offsets in star formation rate, we find that the sizes of galaxies in the highest quartile are 0.27$\pm$0.47 dex larger than galaxies in the lowest quartile, a difference that is not significant. We also show that 50\% of all star formation in galaxies at fixed stellar mass takes place within a narrow range of sizes (0.26 dex). Taken together, these results suggest that there is an abrupt cessation of star formation after galaxies attain particular structural properties. Confirming earlier results, we find that the central stellar density within a fixed physical radius of 1 kpc is the key parameter connecting galaxy morphology and star formation histories: galaxies with high central densities are red and have increasingly lower specific star formation rates, whereas galaxies with low central densities are blue and have a roughly constant (higher) specific star formation rate at a given redshift interval. We find remarkably little scatter in the average trends and a strong evolution in the central density (or equivalently central circular velocity) threshold correlated with quiescence. This quenching threshold decreases by >0.5 dex from z~2 to z~0.7. Neither a compact galaxy size nor a high n are sufficient to assess the likelihood of quiescence for the average galaxy; rather, it is the combination of these two parameters together with stellar mass that results in a unique quenching threshold in central density or velocity.
We report on the discovery of extended Ly-alpha nebulae at z~3.3 in the Hubble Ultra Deep Field (HUDF, ~ 40 kpc X 80 kpc) and behind the Hubble Frontier Fields galaxy cluster MACSJ0416 (~ 40kpc), spatially associated with groups of star-forming galaxies. VLT/MUSE integral field spectroscopy reveals a complex structure with a spatially-varying double peaked Ly-alpha emission. Overall, the spectral profiles of the two Ly-alpha nebulae are remarkably similar, both showing a prominent blue emission, more intense and slightly broader than the red peak. From the first nebula, located in the HUDF, no X-ray emission has been detected, disfavoring the possible presence of AGNs. Spectroscopic redshifts have been derived for 11 galaxies within two arcsec from the nebula and spanning the redshift range 1.037<z<5.97. The second nebula, behind MACSJ0416, shows three aligned star-forming galaxies plausibly associated to the emitting gas. In both systems, the associated galaxies reveal possible intense rest-frame-optical nebular emissions lines [OIII]4959-5007+Hbeta with equivalent widths as high as 1500A rest-frame and star formation rates ranging from a few to tens of solar masses per year. A possible scenario is that of a group of young, star-forming galaxies sources of escaping ionising radiation that induce Ly-alpha fluorescence, therefore revealing the kinematics of the surrounding gas. Also Ly-alpha powered by star-formation and/or cooling radiation may resemble the double peaked spectral properties and the morphology observed here. If the intense blue emission is associated with inflowing gas, then we may be witnessing an early phase of galaxy or a proto-cluster (or group) formation.
We fit the rotation curves of isolated dwarf galaxies to directly measure the
stellar mass-halo mass relation ($M_*-M_{200}$) over the mass range $5 \times
10^5 < M_{*}/{\rm M}_\odot < 10^{8}$. By accounting for cusp-core
transformations due to stellar feedback, we find a monotonic relation with
remarkably little scatter. Such monotonicity implies that abundance matching
should yield a similar $M_*-M_{200}$ if the cosmological model is correct.
Using the 'field galaxy' stellar mass function from the Sloan Digital Sky
Survey (SDSS) and the halo mass function from the $\Lambda$ Cold Dark Matter
Bolshoi simulation, we find remarkable agreement between the two. This holds
down to $M_{200} \sim 5 \times 10^9$ M$_\odot$, and to $M_{200} \sim 5 \times
10^8$ M$_\odot$ if we assume a power law extrapolation of the SDSS stellar mass
function below $M_* \sim 10^7$ M$_\odot$.
However, if instead of SDSS we use the stellar mass function of nearby galaxy
groups, then the agreement is poor. This occurs because the group stellar mass
function is shallower than that of the field below $M_* \sim 10^9$ M$_\odot$,
recovering the familiar 'missing satellites' and 'too big to fail' problems.
Our result demonstrates that both problems are confined to group environments
and must, therefore, owe to 'galaxy formation physics' rather than exotic
cosmology.
Finally, we repeat our analysis for a $\Lambda$ Warm Dark Matter cosmology,
finding that it fails at 68% confidence for a thermal relic mass of $m_{\rm
WDM} < 1.25$ keV, and $m_{\rm WDM} < 2$ keV if we use the power law
extrapolation of SDSS. We conclude by making a number of predictions for future
surveys based on these results.
We present new radial velocity (RV) measurements of old (horizontal branch) and intermediate-age (red clump) stellar tracers in the Carina dwarf spheroidal. They are based on more than 2,200 low-resolution spectra collected with VIMOS at VLT. The targets are faint (20<V<21.5 mag), but the accuracy at the faintest limit is <9 kms-1. These data were complemented with RV measurements either based on spectra collected with FORS2 and FLAMES/GIRAFFE at VLT or available in the literature. We ended up with a sample of 2748 stars and among them 1389 are candidate Carina stars. We found that the intermediate-age stellar component shows a well defined rotational pattern around the minor axis. The western and the eastern side of the galaxy differ by +5 and -4 km s-1 when compared with the main RV peak. The old stellar component is characterized by a larger RV dispersion and does not show evidence of RV pattern. We compared the observed RV distribution with N-body simulations for a former disky dwarf galaxy orbiting a giant MilkyWay-like galaxy (Lokas et al. 2015). We rotated the simulated galaxy by 60 degrees with respect to the major axis, we kept the observer on orbital plane of the dwarf and extracted a sample of stars similar to the observed one. Observed and predicted Vrot/{\sigma} ratios across the central regions are in remarkable agreement. This evidence indicates that Carina was a disky dwarf galaxy that experienced several strong tidal interactions with the Milky Way. Owing to these interactions, Carina transformed from a disky to a prolate spheroid and the rotational velocity transformed into random motions.
Dense cloud cores present chemical differentiation due to the different distribution of C-bearing and N-bearing molecules, the latter being less affected by freeze-out onto dust grains. In this letter we show that two C-bearing molecules, CH$_3$OH and $c$-C$_3$H$_2$, present a strikingly different (complementary) morphology while showing the same kinematics toward the prestellar core L1544. After comparing their distribution with large scale H$_2$ column density N(H$_2$) map from the Herschel satellite, we find that these two molecules trace different environmental conditions in the surrounding of L1544: the $c$-C$_3$H$_2$ distribution peaks close to the southern part of the core, where the surrounding molecular cloud has a N(H$_2$) sharp edge, while CH$_3$OH mainly traces the northern part of the core, where N(H$_2$) presents a shallower tail. We conclude that this is evidence of chemical differentiation driven by different amount of illumination from the interstellar radiation field: in the South, photochemistry maintains more C atoms in the gas phase allowing carbon chain (such as $c$-C$_3$H$_2$) production; in the North, C is mainly locked in CO and methanol traces the zone where CO starts to freeze out significantly. During the process of cloud contraction, different gas and ice compositions are thus expected to mix toward the central regions of the core, where a potential Solar-type system will form. An alternative view on carbon-chain chemistry in star-forming regions is also provided.
We quantify the impact that a variety of galactic and environmental properties have on the quenching of star formation. We collate a sample of $\sim$ 400,000 central and $\sim$ 100,000 satellite galaxies from the Sloan Digital Sky Survey Data Release 7 (SDSS DR7). Specifically, we consider central velocity dispersion ($\sigma_{c}$), stellar, halo, bulge and disk mass, local density, bulge-to-total ratio, group-centric distance and galaxy-halo mass ratio. We develop and apply a new statistical technique to quantify the impact on the quenched fraction ($f_{\rm Quench}$) of varying one parameter, while keeping the remaining parameters fixed. For centrals, we find that the $f_{\rm Quench} - \sigma_{c}$ relationship is tighter and steeper than for any other variable considered. We compare to the Illustris hydrodynamical simulation and the Munich semi-analytic model (L-Galaxies), finding that our results for centrals are qualitatively consistent with their predictions for quenching via radio-mode AGN feedback, hinting at the viability of this process in explaining our observational trends. However, we also find evidence that quenching in L-Galaxies is too efficient and quenching in Illustris is not efficient enough, compared to observations. For satellites, we find strong evidence that environment affects their quenched fraction at fixed central velocity dispersion, particularly at lower masses. At higher masses, satellites behave identically to centrals in their quenching. Of the environmental parameters considered, local density affects the quenched fraction of satellites the most at fixed central velocity dispersion.
We present new KAT-7 observations of the neutral hydrogen (HI) spectral line, and polarized radio continuum emission, in the grand design spiral M83. These observations provide a sensitive probe of the outer disk structure and kinematics, revealing a vast and massive neutral gas distribution that appears to be tightly coupled to the interaction of the galaxy with the environment. We present a new rotation curve extending out to a radius of 50 kpc. Based on our new HI dataset and comparison with multiwavelength data from the literature we consider the impact of mergers on the outer disk and discuss the evolution of M83. We also study the periphery of the HI distribution and reveal a sharp edge to the gaseous disk that is consistent with photoionization or ram pressure from the intergalactic medium (IGM). The radio continuum emission is not nearly as extended as the HI and is restricted to the main optical disk. Despite the relatively low angular resolution we are able to draw broad conclusions about the large-scale magnetic field topology. We show that the magnetic field of M83 is similar in form to other nearby star forming galaxies, and suggest that the disk-halo interface may host a large-scale regular magnetic field.
We present new MUSE observations of a galaxy group probed by a background quasar. The quasar sightline passes between multiple $z=0.28$ galaxies, whilst showing at the same redshift low ionised metal line species, including Ca II, Mg I, Mg II and Fe II. Based on the galaxy redshifts measured from the MUSE data, we estimate the galaxies to be part of a small galaxy group with a halo mass of $\approx6\times10^{12}$ M$_{\odot}$. We use the MUSE data to reveal the two dimensional dynamical properties of the gas and stars in the group galaxies, and relate these to the absorber kinematics. With these data we consider four scenarios for the nature of the gas probed by the sightline absorbers: a co-rotating gas halo associated with a single galaxy within the group; outflowing material from a single group member powered by recent star-formation; and cool dense gas embedded within the intra-group medium. We find that the dynamics, alongside the galaxy impact parameters and star-formation rates, favour the latter, in which the cool gas belongs to the intra-group medium and is likely pressure confined within a warmer gas halo.
We revisited the spectroscopic characteristics of narrow-line Seyfert 1 galaxies (NLS1s) by analysing a homogeneous sample of 296 NLS1s at redshift between 0.028 and 0.345, extracted from the Sloan Digital Sky Survey (SDSS-DR7) public archive. We confirm that NLS1s are mostly characterized by Balmer lines with Lorentzian profiles, lower black hole masses and higher Eddington ratios than classic broad-line Seyfert 1 (BLS1s), but they also appear to be active galactic nuclei (AGNs) contiguous with BLS1s and sharing with them common properties. Strong Fe II emission does not seem to be a distinctive property of NLS1s, as low values of Fe II/H$\beta$ are equally observed in these AGNs. Our data indicate that Fe II and Ca II kinematics are consistent with the one of H$\beta$. On the contrary, O I $\lambda$8446 seems to be systematically narrower and it is likely emitted by gas of the broad-line region more distant from the ionizing source and showing different physical properties. Finally, almost all NLS1s of our sample show radial motions of the narrow-line region highly-ionised gas. The mechanism responsible for this effect is not yet clear, but there are hints that very fast outflows require high continuum luminosities (> $10^{44}$ erg/s) or high Eddington ratios (log(L$_{\rm bol}$/L$_{\rm Edd}$) > -0.1).
CONTEXT. Stellar populations are the building blocks of galaxies including the Milky Way. The majority, if not all extragalactic studies are entangled with the use of stellar population models given the unresolved nature of their observation. Extragalactic systems contain multiple stellar populations with complex star formation histories. However, their study is mainly based upon the principles of simple stellar populations (SSP). Hence, it is critical to examine the validity of SSP models. AIMS. This work aims to empirically test the validity of SSP models. This is done by comparing SSP models against observations of spatially resolved young stellar population in the determination of its physical properties, i.e. age and metallicity. METHODS. Integral field spectroscopy of a young stellar cluster in the Milky Way, NGC 3603, is used to study the properties of the cluster both as a resolved and unresolved stellar population. The unresolved stellar population is analysed using the H$\alpha$ equivalent width as an age indicator, and the ratio of strong emission lines to infer metallicity. In addition, spectral energy distribution (SED) fitting using STARLIGHT, is used to infer these properties from the integrated spectrum. Independently, the resolved stellar population is analysed using the color-magnitude diagram (CMD) for age and metallicity determination. As the SSP model represents the unresolved stellar population, the derived age and metallicity are put to test whether they agree with those derived from resolved stars. RESULTS. The age and metallicity estimate of NGC 3603 derived from integrated spectroscopy are confirmed to be within the range of those derived from the CMD of the resolved stellar population, including other estimates found in the literature. The result from this pilot study supports the reliability of SSP models for studying unresolved young stellar populations.
ID11 is an actively star forming extremely compact galaxy and Ly alpha
emitter at z=3.117 that is gravitationally magnified by a factor of ~17 by the
cluster of galaxies Hubble Frontier Fields AS1063. Its observed properties
resemble those of low luminosity HII galaxies or Giant HII regions like
30-Doradus in the LMC.
Using the tight correlation correlation between the Balmer-line luminosities
and the width of the emission lines (typically L(Hbeta)-sigma(Hbeta)) valid for
HII galaxies and Giant HII regions to estimate its total luminosity, we are
able to measure the lensing amplification of ID11. We obtain an amplification
of 23 +- 11 similar within errors to the value of ~17 estimated or predicted by
the best lensing models of the massive cluster Abell S1063.
We also compiled from the literature luminosities and velocity dispersions
for a set of lensed compact starforming regions. There is more scatter in the
L-sigma correlation for these lensed systems but on the whole the results tend
to support the lensing models estimates of the magnification.
Our result indicates that the amplification can be independently measured
using the L-sigma relation in lensed Giant HII regions or HII galaxies. It also
supports the suggestion, even if lensing model dependent, that the L-sigma
relation is valid for low luminosity high-z objects. Ad-hoc observations of
lensed starforming systems are required to accurately determine the lensing
amplification.
We report the spectroscopic confirmation of 22 new multiply lensed sources behind the {\it Hubble Frontier Field} (HFF) galaxy cluster MACS J0416.1$-$2403 (MACS 0416), using archival data from the Multi Unit Spectroscopic Explorer (MUSE) on the VLT. Combining with previous spectroscopic measurements of other 15 multiply imaged sources, we obtain a sample of 102 secure multiple images with measured redshifts, the largest to date in a single strong lensing system. The newly identified sources are largely low-luminosity Lyman-$\alpha$ emitters with redshift in the range [3.1-6.15]. With such a large number of secure constraints, and a significantly improved sample of galaxy members in the cluster core, we have improved our previous strong lensing model and obtained a robust determination of the projected total mass distribution of MACS 0416. We find evidence of three cored dark-matter halos, adding to the known complexity of this merging system. The total mass density profile, as well as the sub-halo population, are found in good agreement with previous works. We make public an updated redshift catalog for MACS~0416 from our previous spectroscopic campaign, as well as lensing maps (convergence, shear, magnification) in the standard HFF format.
The pressure exerted by massive stars' radiation fields is an important mechanism regulating their formation. Detailed simulation of massive star formation therefore requires an accurate treatment of radiation. However, all published simulations have either used a diffusion approximation of limited validity; have only been able to simulate a single star fixed in space, thereby suppressing potentially-important instabilities; or did not provide adequate resolution at locations where instabilities may develop. To remedy this we have developed a new, highly accurate radiation algorithm that properly treats the absorption of the direct radiation field from stars and the re-emission and processing by interstellar dust. We use our new tool to perform three-dimensional radiation-hydrodynamic simulations of the collapse of massive pre-stellar cores with laminar and turbulent initial conditions and properly resolve regions where we expect instabilities to grow. We find that mass is channeled to the stellar system via gravitational and Rayleigh-Taylor (RT) instabilities, in agreement with previous results using stars capable of moving, but in disagreement with methods where the star is held fixed or with simulations that do not adequately resolve the development of RT instabilities. For laminar initial conditions, proper treatment of the direct radiation field produces later onset of instability, but does not suppress it entirely provided the edges of radiation-dominated bubbles are adequately resolved. Instabilities arise immediately for turbulent pre-stellar cores because the initial turbulence seeds the instabilities. Our results suggest that RT features are significant and should be present around accreting massive stars throughout their formation.
We present the detection of an absorpton feature at $E=8.77^{+0.05}_{-0.06}$ keV in the combined X-ray spectrum of the ultraluminous X-ray source NGC 1313 X-1 observed with XMM-Newton and NuSTAR, significant at the 3$\sigma$ level. If associated with blueshifted ionized iron, the implied outflow velocity is ~0.2$c$ for Fe XXVI, or ~0.25$c$ for Fe XXV. These velocities are similar to the ultrafast outflow seen in absorption recently discovered in this source at lower energies by XMM-Newton, and we therefore conclude that this is an iron component to the same outflow. Photoionization modeling marginally prefers the Fe XXV solution, but in either case the outflow properties appear to be extreme, potentially supporting a super-Eddington hypothesis for NGC 1313 X-1.
Links to: arXiv, form interface, find, astro-ph, recent, 1607, contact, help (Access key information)
The inefficiency of star formation in massive elliptical galaxies is widely believed to be caused by the interactions of an active galactic nucleus (AGN) with the surrounding gas. Achieving a sufficiently rapid reddening of moderately massive galaxies without expelling too many baryons has however proven difficult for hydrodynamical simulations of galaxy formation, prompting us to explore a new model for the accretion and feedback effects of supermassive black holes. For high accretion rates relative to the Eddington limit, we assume that a fraction of the accreted rest mass energy heats the surrounding gas thermally, similar to the `quasar mode' in previous work. For low accretion rates, we invoke a new, pure kinetic feedback model which imparts momentum into the surrounding gas in a stochastic manner. These two modes of feedback are motivated both by theoretical conjectures for the existence of different types of accretion flows as well as recent observational evidence for the importance of kinetic AGN winds in quenching galaxies. We find that a large fraction of the injected kinetic energy in this mode thermalises via shocks in the surrounding gas, thereby providing a distributed heating channel. In cosmological simulations, the resulting model produces red, non star-forming massive elliptical galaxies, and achieves realistic gas fractions, black hole growth histories and thermodynamic profiles in large haloes.
We present observations of the Pisces A and B galaxies with the Advanced Camera for Surveys on the Hubble Space Telescope. Photometry from these images clearly resolve a Red Giant Branch for both objects, demonstrating that they are nearby dwarf galaxies. We describe a Bayesian inferential approach to determining the distance to these galaxies using the magnitude of the tip of the RGB, and then apply this approach to these galaxies. We also provide the full probability distributions for parameters derived using this approach. This reveals the distance to these galaxies as $5.64^{+0.13}_{-0.15} \, {\rm Mpc}$ and $8.89^{+0.75}_{-0.85} \, {\rm Mpc}$ for Pisces A and B, respectively, placing both within the Local Volume but not the Local Group. We estimate the star formation histories of these galaxies, which suggests that they have recently undergone an increase in their star formation rates. Together these yield luminosities for Pisces A and B of $M_V=-11.57^{+0.06}_{-0.05}$ and $-12.9 \pm 0.2$, respectively, and estimated stellar masses of $\log(M_*/M_{\odot})= 7.0^{+0.4}_{-1.7}$ and $7.5^{+0.3}_{-1.8}$. We further show that these galaxies are likely at the boundary between nearby voids and higher-density filamentary structure. This suggests that they are entering a higher-density region from voids, where they would have experienced delayed evolution, consistent with their recent increased star formation rates. If this is indeed the case, they are useful for study as proxies of the galaxies that later evolved into typical Local Group satellite galaxies.
The [CII] 158$\mu$m emission line can arise in all phases of the ISM, therefore being able to disentangle the different contributions is an important yet unresolved problem when undertaking galaxy-wide, integrated [CII] observations. We present a new multi-phase 3D radiative transfer interface that couples Starburst99, a stellar spectrophotometric code, with the photoionisation and astrochemistry codes Mocassin and 3D-PDR. We model entire star forming regions, including the ionised, atomic and molecular phases of the ISM, and apply a Bayesian inference methodology to parametrise how the fraction of the [CII] emission originating from molecular regions, $f_{[CII],mol}$, varies as a function of typical integrated properties of galaxies in the local Universe. The main parameters responsible for the variations of $f_{[CII],mol}$ are specific star formation rate (sSFR), gas phase metallicity, HII region electron number density ($n_e$), and dust mass fraction. For example, $f_{[CII],mol}$ can increase from 60% to 80% when either $n_e$ increases from 10$^{1.5}$ to 10$^{2.5}$cm$^{-3}$, or SSFR decreases from $10^{-9.6}$ to $10^{-10.6}$ yr$^{-1}$. Our model predicts for the Milky Way that $f_{[CII],mol}$$=75.8\pm5.9$%, in agreement with the measured value of 75%. When applying the new prescription to a complete sample of galaxies from the Herschel Reference Survey (HRS), we find that anywhere from 60 to 80% of the total integrated [CII] emission arises from molecular regions.
We investigate the origin of the evolution of the population-averaged size of quenched galaxies (QGs) through a spectroscopic analysis of their stellar ages. The two most favoured scenarios for this evolution are either the size growth of individual galaxies through a sequence of dry minor merger events, or the addition of larger, newly quenched galaxies to the pre-existing population (i.e., a progenitor bias effect). We use the 20k zCOSMOS-bright spectroscopic survey to select bona fide quiescent galaxies at 0.2<z<0.8. We stack their spectra in bins of redshift, stellar mass and size to compute stellar population parameters in these bins through fits to the rest-frame optical spectra and through Lick spectral indices. We confirm a change of behaviour in the size-age relation below and above the ~10^11 MSun stellar mass scale: In our 10.5 < log M*/MSun < 11 mass bin, over the entire redshift window, the stellar populations of the largest galaxies are systematically younger than those of the smaller counterparts, pointing at progenitor bias as the main driver of the observed average size evolution at sub-10^11 MSun masses. In contrast, at higher masses, there is no clear trend in age as a function of galaxy size, supporting a substantial role of dry mergers in increasing the sizes of these most massive QGs with cosmic time. Within the errors, the [alpha/Fe] abundance ratios of QGs are (i) above-solar over the entire redshift range of our analysis, hinting at universally short timescales for the buildup of the stellar populations of QGs, and (ii) similar at all masses and sizes, suggesting similar (short) timescales for the whole QG population and strengthening the role of mergers in the buildup of the most massive QGs in the Universe.
The Central Molecular Zone (CMZ) at the center of our Galaxy is the best template to study star formation processes under extreme conditions, similar to those in high-redshift galaxies. We observed on-the-fly maps of para-H$_{2}$CO transitions at 218 GHz and 291 GHz towards seven Galactic Center clouds. From the temperature-sensitive integrated intensity line ratios of H$_{2}$CO(3$_{2,1}-$2$_{2,0}$)/H$_{2}$CO(3$_{0,3}-$2$_{0,2}$) and H$_{2}$CO(4$_{2,2}-$3$_{2,1}$)/H$_{2}$CO(4$_{0,4}-$3$_{0,3}$) in combination with radiative transfer models, we produce gas temperature maps of our targets. These transitions are sensitive to gas with densities of $\sim$10$^{5}$ cm$^{-3}$ and temperatures <150 K. The measured gas temperatures in our sources are all higher (>40 K) than their dust temperatures ($\sim$25 K). Our targets have a complex velocity structure that requires a careful disentanglement of the different components. We produce temperature maps for each of the velocity components and show that the temperatures of the components differ, revealing temperature gradients in the clouds. Combining the temperature measurements with the integrated intensity line ratio of H$_{2}$CO(4$_{0,4}-$3$_{0,3}$)/H$_{2}$CO(3$_{0,3}-$2$_{0,2}$), we constrain the density of this warm gas to 10$^{4}-$10$^{6}$ cm$^{-3}$. We find a positive correlation of the line width of the main H$_{2}$CO lines with the temperature of the gas, direct evidence for gas heating via turbulence. Our data is consistent with a turbulence heating model with a density of n = 10$^5$ cm$^{-3}$.
We present the first study of GALEX far ultra-violet (FUV) luminosity functions of individual star-forming regions within a sample of 258 nearby galaxies spanning a large range in total stellar mass and star formation properties. We identify ~65,000 star-forming regions (i.e., FUV sources), measure each galaxy's luminosity function, and characterize the relationships between the luminosity function slope (alpha) and several global galaxy properties. A final sample of 82 galaxies with reliable luminosity functions are used to define these relationships and represent the largest sample of galaxies with the largest range of galaxy properties used to study the connection between luminosity function properties and galaxy environment. We find that alpha correlates with global star formation properties, where galaxies with higher star formation rates and star formation rate densities (Sigma_SFR) tend to have flatter luminosity function slopes. In addition, we find that neither stochastic sampling of the luminosity function in galaxies with low-number statistics nor the effects of blending due to distance can fully account for these trends. We hypothesize that the flatter slopes in high Sigma_SFR galaxies is due to higher gas densities and higher star formation efficiencies which result in proportionally greater numbers of bright star-forming regions. Finally, we create a composite luminosity function composed of star-forming regions from many galaxies and find a break in the luminosity function at brighter luminosities. However, we find that this break is an artifact of varying detection limits for galaxies at different distances.
Stellar streams have become central to studies of the interaction histories of nearby galaxies. To characterize the most prominent parts of the stellar stream around the well-known nearby (d = 17 Mpc) edge-on disk galaxy NGC 5907, we have obtained and analyzed new, deep gri Subaru/Suprime-Cam and 3.6 micron Spitzer/Infrared Array Camera (IRAC) observations. Combining the near-infrared 3.6 micron data with visible-light images allows us to use a long wavelength baseline to estimate the metallicity and age of the stellar population along a ~60 kpc long segment of the stream. We have fitted the stellar spectral energy distribution (SED) with a single-burst stellar population synthesis model and we use it to distinguish between the proposed satellite accretion and minor/major merger formation models of the stellar stream around this galaxy. We conclude that a massive minor merger (stellar mass ratio of at least 1:8) can best account for the metallicity of -0.3 inferred along the brightest parts of the stream.
Nuclear rings at centers of barred galaxies exhibit strong star formation activities. They are thought to undergo gravitational instability when sufficiently massive. We approximate them as rigidly-rotating isothermal objects and investigate their gravitational instability. Using a self-consistent field method, we first construct their equilibrium sequences specified by two parameters: alpha corresponding to the thermal energy relative to gravitational potential energy, and R_B measuring the ellipticity or ring thickness. Unlike in the incompressible case, not all values of R_B yield an isothermal equilibrium, and the range of R_B for such equilibria shrinks with decreasing alpha. The density distributions in the meridional plane are steeper for smaller alpha, and well approximated by those of infinite cylinders for slender rings. We also calculate the dispersion relations of nonaxisymmetric modes in rigidly-rotating slender rings with angular frequency Omega_0 and central density rho_max. Rings with smaller alpha are found more unstable with a larger unstable range of the azimuthal mode number. The instability is completely suppressed by rotation when Omega_0 exceeds the critical value. The critical angular frequency is found to be almost constant at ~ 0.7 sqrt(G*rho_c) for alpha > 0.01 and increases rapidly for smaller alpha. We apply our results to a sample of observed star-forming rings and confirm that rings without a noticeable azimuthal age gradient of young star clusters are indeed gravitationally unstable.
We use the Advanced Camera for Surveys on the Hubble Space Telescope to study the rich population of young massive star clusters in the main body of NGC 3256, a merging pair of galaxies with a high star formation rate (SFR) and SFR per unit area ($\Sigma_{\rm{SFR}}$). These clusters have luminosity and mass functions that follow power laws, $dN/dL \propto L^{\alpha}$ with $\alpha = -2.23 \pm 0.07$, and $dN/dM \propto M^{\beta}$ with $\beta = -1.86 \pm 0.34$ for $\tau < 10$ Myr clusters, similar to those found in more quiescent galaxies. The age distribution can be described by $dN/d\tau \propto \tau ^ \gamma$, with $\gamma \approx -0.67 \pm 0.08$ for clusters younger than about a few hundred million years, with no obvious dependence on cluster mass. This is consistent with a picture where $\sim 80 \%$ of the clusters are disrupted each decade in time. We investigate the claim that galaxies with high $\Sigma_{\rm{SFR}}$ form clusters more efficiently than quiescent systems by determining the fraction of stars in bound clusters ($\Gamma$) and the CMF/SFR statistic (CMF is the cluster mass function) for NGC 3256 and comparing the results with those for other galaxies. We find that the CMF/SFR statistic for NGC 3256 agrees well with that found for galaxies with $\Sigma_{\rm{SFR}}$ and SFRs that are lower by $1-3$ orders of magnitude, but that estimates for $\Gamma$ are only robust when the same sets of assumptions are applied. Currently, $\Gamma$ values available in the literature have used different sets of assumptions, making it more difficult to compare the results between galaxies.
In simulations of the standard cosmological model (LCDM), dark matter halos are aspherical. However, so far the asphericity of an individual galaxy's halo has never been robustly established. We use the Jeans equations to define a quantity which robustly characterizes a deviation from rotational symmetry. This quantity is essentially the gravitational torque and it roughly provides the ellipticity projected along the line of sight. We show that the Thirty Meter Telescope (TMT), with a single epoch of observations combined with those of the Gaia space telescope, can distinguish the LCDM value of the torque from zero for each Sculptor-like dwarf galaxy with a confidence between 0 and 5 sigma, depending on the orientation of each halo. With two epochs of observations, TMT will achieve a 5 sigma discovery of torque and so asphericity for most such galaxies, and so will provide a new and powerful test of the LCDM model.
We present $Spitzer$/IRS mid-infrared spectral maps of the Galactic star-forming region M17 as well as IRSF/SIRIUS Br$\gamma$ and Nobeyama 45-m/FOREST $^{13}$CO ($J$=1--0) maps. The spectra show prominent features due to polycyclic aromatic hydrocarbons (PAHs) at wavelengths of 6.2, 7.7, 8.6, 11.3, 12.0, 12.7, 13.5, and 14.2 $\mu$m. We find that the PAH emission features are bright in the region between the HII region traced by Br$\gamma$ and the molecular cloud traced by $^{13}$CO, supporting that the PAH emission originates mostly from photo-dissociation regions. Based on the spatially-resolved maps, we examine spatial variations of the PAH properties in detail. As a result, we find that the interband ratio of PAH 7.7 $\mu$m/PAH 11.3 $\mu$m varies locally near M17SW, but rather independently of the distance from the OB stars in M17, suggesting that the ionization degree of PAHs is mainly controlled by local conditions rather than the global UV environments determined by the OB stars in M17. We also find that the interband ratios of the PAH 12.0 $\mu$m, 12.7 $\mu$m, 13.5 $\mu$m, and 14.2 $\mu$m features to the PAH 11.3 $\mu$m feature are high near the M17 center, which suggests structural changes of PAHs through processing due to intense UV radiation, producing abundant edgy irregular PAHs near the M17 center.
We present new observations of the galaxy cluster 3C 129 obtained with the Sardinia Radio Telescope in the frequency range 6000-7200 MHz, with the aim to image the large-angular-scale emission at high-frequency of the radio sources located in this cluster of galaxies. The data were acquired using the recently-commissioned ROACH2-based backend to produce full-Stokes image cubes of an area of 1 deg x 1 deg centered on the radio source 3C 129. We modeled and deconvolved the telescope beam pattern from the data. We also measured the instrumental polarization beam patterns to correct the polarization images for off-axis instrumental polarization. Total intensity images at an angular resolution of 2.9 arcmin were obtained for the tailed radio galaxy 3C 129 and for 13 more sources in the field, including 3C 129.1 at the galaxy cluster center. These data were used, in combination with literature data at lower frequencies, to derive the variation of the synchrotron spectrum of 3C 129 along the tail of the radio source. If the magnetic field is at the equipartition value, we showed that the lifetimes of radiating electrons result in a radiative age for 3C 129 of t_syn = 267 +/- 26 Myrs. Assuming a linear projected length of 488 kpc for the tail, we deduced that 3C 129 is moving supersonically with a Mach number of M=v_gal/c_s=1.47. Linearly polarized emission was clearly detected for both 3C 129 and 3C 129.1. The linear polarization measured for 3C 129 reaches levels as high as 70% in the faintest region of the source where the magnetic field is aligned with the direction of the tail.
We report a new 1-pc (30") resolution CS($J=2-1$) line map of the central 30 pc of the Galactic Center (GC), made with the Nobeyama 45m telescope. We revisit our previous study of the extraplanar feature called polar arc (PA), which is a molecular cloud located above SgrA* with a velocity gradient perpendicular to the Galactic plane. We find that the PA can be traced back to the Galactic disk. This provides clues of the launching point of the PA , roughly $6\times10^{6}$ years ago. Implications of the dynamical time scale of the PA might be related to the Galactic Center Lobe (GCL) at parsec scale. Our results suggest that in the central 30 pc of the GC, the feedback from past explosions could alter the orbital path of the molecular gas down to the central tenth of parsec. In the follow-up work of our new CS($J=2-1$) map, we also find that near the systemic velocity, the molecular gas shows an extraplanar hourglass-shaped feature (HG-feature) with a size of $\sim$13 pc. The latitude-velocity diagrams show that the eastern edge of the HG-feature is associated with an expanding bubble B1, $\sim$7 pc away from SgrA*. The dynamical time scale of this bubble is $\sim3\times10^{5}$ years. This bubble is interacting with the 50 km s$^{-1}$ cloud. Part of the molecular gas from the 50 km s$^{-1}$ cloud was swept away by the bubble to $b=-0.2deg$. The western edge of the HG-feature seems to be the molecular gas entrained from the 20 km s$^{-1}$ cloud towards the north of the Galactic disk. Our results suggest a fossil explosion in the central 30 pc of the GC a few 10$^{5}$ years ago.
We analyze new high spatial resolution (~60 pc) ALMA CO(2-1) observations of the isolated luminous infrared galaxy ESO 320-G030 (d=48 Mpc) in combination with ancillary HST optical and near-IR imaging as well as VLT/SINFONI near-IR integral field spectroscopy. We detect a high-velocity (~450 km/s) spatially resolved (size~2.5 kpc; dynamical time ~3 Myr) massive (~10^7 Msun; mass rate~2-8 Msun/yr) molecular outflow originated in the central ~250 pc. We observe a clumpy structure in the outflowing cold molecular gas with clump sizes between 60 and 150 pc and masses between 10^5.5 and 10^6.4 Msun. The mass of the clumps decreases with increasing distance, while the velocity is approximately constant. Therefore, both the momentum and kinetic energy of the clumps decrease outwards. In the innermost (~100 pc) part of the outflow, we measure a hot-to-cold molecular gas ratio of 7x10^-5, which is similar to that measured in other resolved molecular outflows. We do not find evidence of an ionized phase in this outflow. The nuclear IR and radio properties are compatible with strong and highly obscured star-formation (A_k ~ 4.6 mag; SFR~15 Msun/yr). We do not find any evidence for the presence of an active galactic nucleus. We estimate that supernova explosions in the nuclear starburst (\nu(SN) ~ 0.2 yr^-1) can power the observed molecular outflow. The kinetic energy and radial momentum of the cold molecular phase of the outflow correspond to about 2% and 20%, respectively, of the supernovae output. The cold molecular outflow velocity is lower than the escape velocity, so the gas will likely return to the galaxy disk. The mass loading factor is ~0.1-0.5, so the negative feedback due to this star-formation powered molecular outflow is probably limited.
Spiral galaxies are thought to acquire their gas through a protracted infall phase resulting in the inside-out growth of their associated discs. For field spirals, this infall occurs in the lower density environments of the cosmic web. The overall infall rate, as well as the galactocentric radius at which this infall is incorporated into the star-forming disc, plays a pivotal role in shaping the characteristics observed today. Indeed, characterising the functional form of this spatio-temporal infall in-situ is exceedingly difficult, and one is forced to constrain these forms using the present day state of galaxies with model or simulation predictions. We present the infall rates used as input to a grid of chemical evolution models spanning the mass spectrum of discs observed today. We provide a systematic comparison with alternate analytical infall schemes in the literature, including a first comparison with cosmological simulations. Identifying the degeneracies associated with the adopted infall rate prescriptions in galaxy models is an important step in the development of a consistent picture of disc galaxy formation and evolution.
A complete sample of 96 faint ($S > 0.5$ mJy) radio galaxies is selected from
the Tenth Cambridge (10C) survey at 15.7~GHz. Optical spectra are used to
classify 17 of the sources as high-excitation or low-excitation radio galaxies
(HERGs and LERGs respectively), for the remaining sources three other methods
are used; these are optical compactness, X-ray observations and mid-infrared
colour--colour diagrams. 32 sources are HERGs and 35 are LERGs while the
remaining 29 sources could not be classified. We find that the 10C HERGs tend
to have higher 15.7-GHz flux densities, flatter spectra, smaller linear sizes
and be found at higher redshifts than the LERGs. This suggests that the 10C
HERGs are more core dominated than the LERGs.
Lower-frequency radio images, linear sizes and spectral indices are used to
classify the sources according to their radio morphology; 18 are Fanaroff and
Riley type I or II sources, a further 13 show some extended emission, and the
remaining 65 sources are compact and are referred to as FR0 sources. The FR0
sources are sub-divided into compact, steep-spectrum (CSS) sources (13 sources)
or GHz-peaked spectrum (GPS) sources (10 sources) with the remaining 42 in an
unclassified class. FR0 sources are more dominant in the subset of sources with
15.7-GHz flux densities $<$1 mJy, consistent with the previous result that the
fainter 10C sources have flatter radio spectra.
The properties of the 10C sources are compared to the higher-flux density
Australia Telescope 20 GHz (AT20G) survey. The 10C sources are found at similar
redshifts to the AT20G sources but have lower luminosities. The nature of the
high-frequency selected objects change as flux density decreases; at high flux
densities the objects are primarily quasars, while at low flux densities radio
galaxies dominate.
We introduce a parametric family of models to characterize the properties of astrophysical systems in a quasi-stationary evolution under the incidence evaporation. We start from an one-particle distribution $f_{\gamma}\left(\mathbf{q},\mathbf{p}|\beta,\varepsilon_{s}\right)$ that considers an appropriate deformation of Maxwell-Boltzmann form with inverse temperature $\beta$, in particular, a power-law truncation at the scape energy $\varepsilon_{s}$ with exponent $\gamma>0$. This deformation is implemented using a generalized $\gamma$-exponential function obtained from the \emph{fractional integration} of ordinary exponential. As shown in this work, this proposal generalizes models of tidal stellar systems that predict particles distributions with \emph{isothermal cores and polytropic haloes}, e.g.: Michie-King models. We perform the analysis of thermodynamic features of these models and their associated distribution profiles. A nontrivial consequence of this study is that profiles with isothermal cores and polytropic haloes are only obtained for low energies whenever deformation parameter $\gamma<\gamma_{c}\simeq 2.13$.
The rotational spectra of thioisocyanic acid (HNCS), and its three energetic isomers (HSCN, HCNS, and HSNC) have been observed at high spectral resolution by a combination of chirped-pulse and Fabry-P\'{e}rot Fourier-transform microwave spectroscopy between 6 and 40~GHz in a pulsed-jet discharge expansion. Two isomers, thiofulminic acid (HCNS) and isothiofulminic acid (HSNC), calculated here to be 35-37~kcal/mol less stable than the ground state isomer HNCS, have been detected for the first time. Precise rotational, centrifugal distortion, and nitrogen hyperfine coupling constants have been determined for the normal and rare isotopic species of both molecules; all are in good agreement with theoretical predictions obtained at the coupled cluster level of theory. On the basis of isotopic spectroscopy, precise molecular structures have been derived for all four isomers by correcting experimental rotational constants for the effects of rotation-vibration calculated theoretically. Formation and isomerization pathways have also been investigated; the high abundance of HSCN relative to ground state HNCS, and the detection of strong lines of SH using CH$_3$CN and H$_2$S, suggest that HSCN is preferentially produced by the radical-radical reaction HS + CN. A radio astronomical search for HSCN and its isomers has been undertaken toward the high-mass star-forming region Sgr B2(N) in the Galactic Center with the 100 m Green Bank Telescope. While we find clear evidence for HSCN, only a tentative detection of HNCS is proposed, and there is no indication of HCNS or HSNC at the same rms noise level. HSCN, and tentatively HNCS, displays clear deviations from a single-excitation temperature model, suggesting weak masing may be occurring in some transitions in this source.
Dusty star-forming galaxies are among the most prodigious systems at high
redshift (z>1), characterized by high star formation rates and huge dust
reservoirs. The bright end of this population has been well characterized in
recent years, but considerable uncertainties remain for fainter dusty
star-forming galaxies, which are responsible for the bulk of star formation at
high redshift and thus play a key role in galaxy growth and evolution.
In this first paper of our series, we describe our methods for finding high
redshift faint dusty galaxies using millimeter observations with ALMA. We
obtained ALMA 1.1 mm mosaic images for three strong-lensing galaxy clusters
from the Frontier Fields survey. The 2'x2' mosaics overlap with the deep HST
WFC3/IR footprints and encompass the high magnification regions of each
cluster. The combination of extremely high ALMA sensitivity and the
magnification power of these clusters allows us to systematically probe the
sub-mJy population of dusty star-forming galaxies over a large surveyed area.
We present a description of the reduction and analysis of the ALMA continuum
observations for the galaxy clusters Abell 2744 (z=0.308), MACSJ0416.1-2403
(z=0.396) and MACSJ1149.5+2223 (z=0.543), for which we reach observed rms
sensitivities of 55, 59 and 71 $\mu$Jy beam$^{-1}$ respectively. We detect 12
dusty star-forming galaxies at S/N>5.0 across the three clusters, all of them
presenting coincidence with near-infrared detected counterparts in the HST
images. None of the sources fall close to the lensing caustics, thus they are
not strongly lensed. The observed 1.1 mm flux densities for the total sample of
galaxies range from 0.41 to 2.82 mJy, with observed effective radii spanning
0.05" to 0.37". The lensing-corrected sizes of the detected sources appear to
be in the same range as those measured in brighter samples, albeit with
possibly larger dispersion.
We analyze extragalactic extinction pro?les derived through gamma-ray burst afterglows, using a dust model speci?cally constructed on the assumption that dust grains are not immutable but respond time-dependently to the local physics. Such a model includes core-mantle spherical particles of mixed chemical composition (silicate core, sp2 and sp3 carbonaceous layers), and an additional molecular component, in the form of free-flying polycyclic aromatic hydrocarbons. We fit most of the observed extinction pro?les. Failures occur for lines of sight presenting remarkable rises blueward the bump. We find a tendency in the carbon chemical structure to become more aliphatic with the galactic activity, and to some extent with increasing redshifts. Moreover, the contribution of the moleclar component to the total extinction is more important in younger objects. The results of the ?tting procedure (either successes and failures) may be naturally interpreted through an evolutionary prescription based on the carbon cycle in the interstellar medium of galaxies.
Observations by the INTEGRAL satellite reveal that the Galaxy glows with the radiation from the annihilation of $(5.0_{-1.5}^{+1.0}) \times 10^{43}$ electron-positron pairs every second. Constrained to be injected into the interstellar medium (ISM) at only mildly relativistic energies, it is highly plausible most positrons originate from the $\beta^+$ decay of radionuclides synthesised in stars or supernovae. However, none of the initially most likely candidates -- massive stars, core-collapse (CC) supernovae (SNe) or ordinary thermonuclear supernovae (SNe Ia) -- have Galactic distributions that match the spatial distribution of positron injection across the Milky Way. Here we show that a class of transient positron source occurring in stars of age >5 Gyr can explain the global distribution of positron annihilation in the Galaxy. Such sources, occurring at a present Galactic rate $\sim$ 0.002 year$^{-1}$ and typically synthesising $\sim$ 0.03 solar masses of the $\beta^+$-unstable radionuclide $^{44}$Ti, can simultaneously explain the absolute positron luminosity of the Galaxy and the abundance of $^{44}$Ca in mainstream solar system material. Our binary evolution models show that mergers of helium-white dwarf (WD) and carbon-oxygen (CO) WD binary systems satisfy these age and rate requirements and plausibly match the $^{44}$Ti yield requirements. The $^{56}$Ni yield, delay time, and rates of these mergers strongly suggests they are the sub-luminous, thermonuclear SNe known as SN1991bg-like (SNe 91bg). These supernovae are, therefore, likely the main source of Galactic positrons. ONeMg WDs from the same WD population plausibly birth (via accretion induced collapse) the millisecond pulsars plausibly responsible for the 'Galactic Centre Excess'.
We compile an updated list of 28 independent measurements of the Hubble parameter $H(z)$ between redshifts $0.1 \leq z \leq 2.36$ and use them to place constraints on model parameters of constant and time-varying dark energy cosmological models, both spatially flat and curved. We use five models to measure the redshift of the cosmological deceleration-acceleration transition, $z_{\rm da}$, from these $H(z)$ data. Within the error bars, the measured $z_{\rm da}$ are insensitive to the model used, depending only on the value assumed for the Hubble constant $H_0$. The weighted mean of our measurements is $z_{\rm da} = 0.74 \pm 0.06\ (0.86 \pm 0.04)$ for $H_0 = 68 \pm 2.8\ (73.8 \pm 2.4)$ km s$^{-1}$ Mpc$^{-1}$ and should provide a reasonably model-independent estimate of this cosmological parameter. The $H(z)$ data are consistent with the standard spatially-flat $\Lambda$CDM cosmological model but do not rule out non-flat models or dynamical dark energy models.
Ionization-recombination balance in dense interstellar and circumstellar environments is a key factor for a variety of important physical processes, such as chemical reactions, dust charging and coagulation, coupling of the gas with magnetic field and the development of magnetorotational instability in protoplanetary disks. We present a self-consistent analytical model which allows us to exactly calculate abundances of charged species in dusty gas, in the regime where the dust-phase recombination dominates over the gas-phase recombination. The model is employed to verify applicability of a conventional approximation of low dust charges in protoplanetary disks, and to discuss the implications for the dust coagulation and the development of the "dead zone" in the disk. Furthermore, the importance of mutually consistent models for the ionization and dust evolution is addressed: These processes are coupled via several mechanisms operating in the disk, and therefore their interplay can be crucial for the ultimate dust evolution and the formation of macroscopic bodies.
Links to: arXiv, form interface, find, astro-ph, recent, 1607, contact, help (Access key information)
We present here 21 cm predictions from high dynamic range simulations for a range of reionization histories that have been tested against available Lyman-alpha and CMB data. We assess the observability of the predicted spatial 21 cm fluctuations by ongoing and upcoming experiments in the late stages of reionization in the limit in which the hydrogen spin temperature is significantly larger than the CMB temperature. Models consistent with the available Lyman-alpha data and CMB measurement of the Thomson optical depth predict typical values of 10--20 mK^2 for the variance of the 21 cm brightness temperature at redshifts z=7--10 at scales accessible to ongoing and upcoming experiments (k < 1 h/cMpc). This is only a factor of a few below the sensitivity claimed to have been already reached by ongoing experiments. Our different models for the reionization history make markedly different predictions for the redshift evolution and thus frequency dependence of the 21 cm power spectrum and should be easily discernible by LOFAR (and later HERA and SKA1) at their design sensitivity. Our simulations have sufficient resolution to assess the effect of high-density Lyman limit systems that can self-shield against ionizing radiation and stay 21 cm bright even if the hydrogen in their surroundings is highly ionized. Our simulations predict that including the effect of the self-shielded gas in highly ionized regions reduces the large scale 21 cm power by about 30%.
We report on the search for galaxies in the proximity of two very metal-poor gas clouds at z~3 towards the quasar Q0956+122. With a 5-hour MUSE integration in a ~500x500 kpc^2 region centred at the quasar position, we achieve a >80% complete spectroscopic survey of continuum-detected galaxies with m<25 mag and Ly{\alpha} emitters with luminosity L>3e41 erg/s. We do not identify galaxies at the redshift of a z~3.2 Lyman limit system (LLS) with log Z/Zsun = -3.35 +/- 0.05, placing this gas cloud in the intergalactic medium or circumgalactic medium of a galaxy below our sensitivity limits. Conversely, we detect five Ly{\alpha} emitters at the redshift of a pristine z~3.1 LLS with log Z/Zsun < -3.8, while ~0.4 sources were expected given the z~3 Ly{\alpha} luminosity function. Both this high detection rate and the fact that at least three emitters appear aligned in projection with the LLS suggest that this pristine cloud is tracing a gas filament that is feeding one or multiple galaxies. Our observations uncover two different environments for metal-poor LLSs, implying a complex link between these absorbers and galaxy halos, which ongoing MUSE surveys will soon explore in detail. Moreover, in agreement with recent MUSE observations, we detected a ~90 kpc Ly{\alpha} nebula at the quasar redshift and three Ly{\alpha} emitters reminiscent of a "dark galaxy" population.
We present a catalogue of nearly 3,000 submillimetre sources detected at 850um over ~5 square degrees surveyed as part of the James Clerk Maxwell Telescope (JCMT) SCUBA-2 Cosmology Legacy Survey (S2CLS). This is the largest survey of its kind at 850um, probing a meaningful cosmic volume at the peak of star formation activity and increasing the sample size of submillimetre galaxies selected at 850um by an order of magnitude. We describe the wide 850um survey component of S2CLS, which covers the key extragalactic survey fields: UKIDSS-UDS, COSMOS, Akari-NEP, Extended Groth Strip, Lockman Hole North, SSA22 and GOODS-North. The average 1-sigma depth of S2CLS is 1.2 mJy/beam, approaching the SCUBA-2 850um confusion limit, which we determine to be ~0.8 mJy/beam. We measure the single dish 850um number counts to unprecedented accuracy, reducing the Poisson errors on the differential counts to approximately 4% at S_850~3mJy. With several independent fields, we investigate field-to-field variance, finding that the number counts on 0.5-1 degree scales are generally within 50% of the S2CLS mean for S_850>3mJy, with scatter consistent with the Poisson and estimated cosmic variance uncertainties, although there is a marginal (2-sigma) density enhancement in the GOODS-North field. The observed number counts are in reasonable agreement with recent phenomenological and semi-analytic models. Finally, the large solid angle of S2CLS allows us to measure the bright-end counts: at S_850>10mJy there are approximately ten sources per square degree, and we detect the distinctive up-turn in the number counts indicative of the detection of local sources of 850um emission and strongly lensed high-redshift galaxies. Here we describe the data collection and reduction procedures and present calibrated maps and a catalogue of sources; these are made publicly available.
The evolution of the number density of galaxies in the universe, and thus also the total number of galaxies, is a fundamental question with implications for a host of astrophysical problems including galaxy evolution and cosmology. However there has never been a detailed study of this important measurement, nor a clear path to answer it. To address this we use observed galaxy stellar mass functions up to $z\sim8$ to determine how the number densities of galaxies changes as a function of time and mass limit. We show that the increase in the total number density of galaxies ($\phi_{\rm T}$), more massive than M$_{*} = 10^{6}$ M_0, decreases as $\phi_{\rm T} \sim t^{-1}$, where $t$ is the age of the universe. We further show that this evolution turns-over and rather increases with time at higher mass lower limits of M$_{*}>10^{7}$ M_0. By using the M$_{*}=10^{6}$ M_0 lower limit we further show that the total number of galaxies in the universe up to $z = 10$ is $2.0^{+0.7}_{-0.6} \times 10^{12}$ (two trillion), almost a factor of ten higher than would be seen in an all sky survey at Hubble Ultra-Deep Field depth. We discuss the implications for these results for galaxy evolution, as well as compare our results with the latest models of galaxy formation. These results also reveal that the cosmic background light in the optical and near-infrared likely arise from these unobserved faint galaxies. We also show how these results solve the question of why the sky at night is dark, otherwise known as Olbers' paradox.
This work presents the properties of 42 objects in the group of the most luminous, highest star formation rate LINERs at z = 0.04 - 0.11. We obtained long-slit spectroscopy of the nuclear regions for all sources, and FIR data (Herschel and IRAS) for 13 of them. We measured emission line intensities, extinction, stellar populations, stellar masses, ages, AGN luminosities, and star-formation rates. We find considerable differences from other low-redshift LINERs, in terms of extinction, and general similarity to star forming (SF) galaxies. We confirm the existence of such luminous LINERs in the local universe, after being previously detected at z ~ 0.3 by Tommasin et al. (2012). The median stellar mass of these LINERs corresponds to 6 - 7 $\times$ 10$^{10}$M$_{\odot}$ which was found in previous work to correspond to the peak of relative growth rate of stellar populations and therefore for the highest SFRs. Other LINERs although showing similar AGN luminosities have lower SFR. We find that most of these sources have LAGN ~ LSF suggesting co-evolution of black hole and stellar mass. In general among local LINERs being on the main-sequence of SF galaxies is related to their AGN luminosity.
We present 3D magnetohydrodynamic numerical simulations of the adiabatic interaction of a shock with a dense, filamentary cloud. We investigate the effects of various filament lengths and orientations on the interaction using different orientations of the magnetic field, and vary the Mach number of the shock, the density contrast of the filament, and the plasma beta, in order to determine their effect on the evolution and lifetime of the filament. We find that in a parallel magnetic field filaments have longer lifetimes if they are orientated more 'broadside' to the shock front, and that an increase in the density contrast hastens the destruction of the cloud, in terms of the modified cloud-crushing time-scale, tcs. The combination of a mild shock and a perpendicular or oblique field provides the best condition for extending the life of the filament, with some filaments able to survive almost indefinitely since they are cocooned by the magnetic field. A high value for the density contrast does not initiate large turbulent instabilities in either the perpendicular or oblique field cases but rather draws the filament out into long tendrils which may eventually fragment. In addition, flux ropes are only formed in parallel magnetic fields. The length of the filament is, however, not as important for the evolution and destruction of a filament.
We present the Green Bank Telescope absorption survey of cold atomic hydrogen ($\lesssim 300$K) in the inner halo of low-redshift galaxies. The survey aims to characterize the cold gas distribution and to address where condensation - the process where ionized gas accreted by galaxies condenses into cold gas within the disks of galaxies - occurs. Our sample consists of 16 galaxy-quasar pairs with impact parameters of $\le$ 20kpc. We detected an HI absorber associated with J0958+3222 (NGC 3067) and HI emission from six galaxies. We also found two \ion{Ca}{2} absorption system in the archival SDSS data associated with galaxies J0958+3222 and J1228+3706, although the sample was not selected based on the presence of metals in absorption. Our detection rate of HI absorbers with optical depths of $\ge 0.06$ is $\sim$7\%. We also find that cold HI phase ($\lesssim$300K) is 44($\pm$18)\% of the total atomic gas in the sightline probing J0958+3222. We find no correlation between the peak optical depth and impact parameter or stellar and HI radii normalized impact parameters, $\rho/\rm R_{90}$ and $\rho/\rm R_{HI}$. We conclude that the process of condensation of inflowing gas into cold ($\lesssim$ 300K) HI occurs at the $\rho << 20$kpc. However, the warmer phase of neutral gas (T $\sim$ 1000K) can exists out to much larger distances as seen in emission maps. Therefore, the process of condensation of warm to cold HI is likely occurring in stages from ionized to warm HI in the inner halo and then to cold HI very close to the galaxy disk.
A galaxy group catalog is constructed from the 2MASS Redshift Survey (2MRS) with the use of a halo-based group finder. The halo mass associated with a group is estimated using a `GAP' method based on the luminosity of the central galaxy and its gap with other member galaxies. Tests using mock samples shows that this method is reliable, particularly for poor systems containing only a few members. On average 80% of all the groups have completeness >0.8, and about 65% of the groups have zero contamination. Halo masses are estimated with a typical uncertainty $\sim 0.35\,{\rm dex}$. The application of the group finder to the 2MRS gives 29,904 groups from a total of 43,246 galaxies at $z \leq 0.08$, with 5,286 groups having two or more members. Some basic properties of this group catalog is presented, and comparisons are made with other groups catalogs in overlap regions. With a depth to $z\sim 0.08$ and uniformly covering about 91% of the whole sky, this group catalog provides a useful data base to study galaxies in the local cosmic web, and to reconstruct the mass distribution in the local Universe.
We investigate a correlation between star-formation rate (SFR) and stellar mass for Halpha emission line galaxies (HAEs) in one of the richest proto-clusters ever known at z~2.5, USS 1558-003 proto-cluster. This study is based on a 9.7-hour narrow-band imaging data with MOIRCS on the Subaru telescope. We are able to construct a sample, in combination with additional H-band data taken with WFC3 on Hubble Space Telescope (HST), of 100 HAEs reaching the dust-corrected SFRs down to 3 Msun/yr and the stellar masses down to $10^{8.0}$ Msun. We find that while the star-forming galaxies with >$10^{9.3}$ Msun are located on the universal SFR-mass main sequence irrespective of the environment, less massive star-forming galaxies with <$10^{9.3}$ Msun show a significant upward scatter from the main sequence in this proto-cluster. This suggests that some less massive galaxies are in a starburst phase, although we do not know yet if this is due to environmental effects.
We present the results of 325 MHz GMRT observations of a super-cluster field, known to contain five Abell clusters at redshift $z \sim 0.2$. We achieve a nominal sensitivity of $34\,\mu$Jy beam$^{-1}$ toward the phase centre. We compile a catalogue of 3257 sources with flux densities in the range $183\,\mu\rm{Jy}\,-\,1.5\,\rm{Jy}$ within the entire $\sim 6.5$ square degree field of view. Subsequently, we use available survey data at other frequencies to derive the spectral index distribution for a sub-sample of these sources, recovering two distinct populations -- a dominant population which exhibit spectral index trends typical of steep-spectrum synchrotron emission, and a smaller population of sources with typically flat or rising spectra. We identify a number of sources with ultra-steep spectra or rising spectra for further analysis, finding two candidate high-redshift radio galaxies and three gigahertz-peaked-spectrum radio sources. Finally, we derive the Euclidean-normalised differential source counts using the catalogue compiled in this work, for sources with flux densities in excess of $223 \, \mu$Jy. Our differential source counts are consistent with both previous observations at this frequency and models of the low-frequency source population. These represent the deepest source counts yet derived at 325 MHz. Our source counts exhibit the well-known flattening at mJy flux densities, consistent with an emerging population of star-forming galaxies; we also find marginal evidence of a downturn at flux densities below $308 \, \mu$Jy, a feature so far only seen at 1.4 GHz.
Galaxies which strongly deviate from the radio-far IR correlation are of great importance for studies of galaxy evolution as they may be tracing early, short-lived stages of starbursts and active galactic nuclei (AGNs). The most extreme FIR-excess galaxy NGC1377 has long been interpreted as a young dusty starburst, but millimeter observations of CO lines revealed a powerful collimated molecular outflow which cannot be explained by star formation alone. We present new radio observations at 1.5 and 10 GHz obtained with the Jansky Very Large Array (JVLA) and Chandra X-ray observations towards NGC1377. The observations are compared to synthetic starburst models to constrain the properties of the central energy source. We obtained the first detection of the cm radio continuum and X-ray emission in NGC1377. We find that the radio emission is distributed in two components, one on the nucleus and another offset by 4$"$.5 to the South-West. We confirm the extreme FIR-excess of the galaxy, with a $q_\mathrm{FIR}\simeq$4.2, which deviates by more than 7-$\sigma$ from the radio-FIR correlation. Soft X-ray emission is detected on the off-nucleus component. From the radio emission we estimate for a young ($<10$ Myr) starburst a star formation rate SFR$<$0.1 M$_\odot$ yr$^{-1}$. Such a SFR is not sufficient to power the observed IR luminosity and to drive the CO outflow. We find that a young starburst cannot reproduce all the observed properties of the nucleus of NGC1377. We suggest that the galaxy may be harboring a radio-quiet, obscured AGN of 10$^6$M$_\odot$, accreting at near-Eddington rates. We speculate that the off-nucleus component may be tracing an hot-spot in the AGN jet.
Physical conditions that could render a core starless(in the local Universe) is the subject of investigation in this work. To this end we studied the evolution of four starless cores, B68, L694-2, L1517B, L1689, and L1521F, a VeLLO. The density profile of a typical core extracted from an earlier simulation developed to study core-formation in a molecular cloud was used for the purpose. We demonstrate - (i) cores contracted in quasistatic manner over a timescale on the order of $\sim 10^{5}$ years. Those that remained starless did briefly acquire a centrally concentrated density configuration that mimicked the density profile of a unstable Bonnor Ebert sphere before rebounding, (ii) three of our test cores viz. L694-2, L1689-SMM16 and L1521F remained starless despite becoming thermally super-critical. On the contrary B68 and L1517B remained sub-critical; L1521F collapsed to become a VeLLO only when gas-cooling was enhanced by increasing the size of dust-grains. This result is robust, for other cores viz. B68, L694-2, L1517B and L1689 that previously remained starless could also be similarly induced to collapse. Our principle conclusions are : (a) acquiring the thermally super-critical state does not ensure that a core will necessarily become protostellar, (b) potentially star-forming cores (VeLLO L1521F here), could be experiencing coagulation of dust-grains that must enhance the gas-dust coupling and in turn lower the gas temperature, thereby assisting collapse. This hypothesis appears to have some observational support, and (c) depending on its dynamic state at any given epoch, a core could appear to be pressure-confined, gravitationally/virially bound, suggesting that gravitational/virial boundedness of a core is insufficient to ensure it will form stars, though it is crucial for gas in a contracting core to cool efficiently so it can collapse further to become protostellar.
We develop a galaxy cluster finding algorithm based on spectral clustering technique to identify optical counterparts and estimate optical redshifts for X-ray selected cluster candidates. As an application, we run our algorithm on a sample of X-ray cluster candidates selected from the third XMM-Newton serendipitous source catalog (3XMM-DR5) that are located in the Stripe 82 of the Sloan Digital Sky Survey (SDSS). Our method works on galaxies described in the color-magnitude feature space. We begin by examining 45 galaxy clusters with published spectroscopic redshifts in the range of 0.1 to 0.8 with a median of 0.36. As a result, we are able to identify their optical counterparts and estimate their photometric redshifts, which have a typical accuracy of 0.025 and agree with the published ones. Then, we investigate another 40 X-ray cluster candidates (from the same cluster survey) with no redshift information in the literature and found that 12 candidates are considered as galaxy clusters in the redshift range from 0.29 to 0.76 with a median of 0.57. These systems are newly discovered clusters in X-rays and optical data. Among them 7 clusters have spectroscopic redshifts for at least one member galaxy.
We report the discovery and analysis of the most metal-poor damped Lyman-alpha (DLA) system currently known, which also displays the Lyman series absorption lines of neutral deuterium. The average [O/H] abundance of this system is [O/H] = -2.804 +/- 0.015, which includes an absorption component with [O/H] = -3.07 +/- 0.03. Despite the unfortunate blending of many weak D I absorption lines, we report a precise measurement of the deuterium abundance of this system. Using the six highest quality and self-consistently analyzed measures of D/H in DLAs, we report tentative evidence for a subtle decrease of D/H with increasing metallicity. This trend must be confirmed with future high precision D/H measurements spanning a range of metallicity. A weighted mean of these six independent measures provides our best estimate of the primordial abundance of deuterium, 10^5 (D/H)_P = 2.547 +/- 0.033 (log_10 (D/H)_P = -4.5940 +/- 0.0056). We perform a series of detailed Monte Carlo calculations of Big Bang nucleosynthesis (BBN) that incorporate the latest determinations of several key nuclear cross sections, and propagate their associated uncertainty. Combining our measurement of (D/H)_P with these BBN calculations yields an estimate of the cosmic baryon density, 100 Omega_B,0 h^2(BBN) = 2.156 +/- 0.020, if we adopt the most recent theoretical determination of the d(p,gamma)3He reaction rate. This measure of Omega_B,0 h^2 differs by ~2.3 sigma from the Standard Model value estimated from the Planck observations of the cosmic microwave background. Using instead a d(p,gamma)3He reaction rate that is based on the best available experimental cross section data, we estimate 100 Omega_B,0 h^2(BBN) = 2.260 +/- 0.034, which is in somewhat better agreement with the Planck value. Forthcoming measurements of the crucial d(p,gamma)3He cross section may shed further light on this discrepancy.
The existence of satellite galaxy planes poses a major challenge for the standard picture of structure formation with non-baryonic dark matter. Recently Tully et al. (2015) reported the discovery of two almost parallel planes in the nearby Cen A group using mostly high-mass galaxies (M$_B$ $<$ -10 mag) in their analysis. Our team detected a large number of new group member candidates in the Cen A group (M\"uller et al. 2016). This dwarf galaxy sample combined with other recent results from the literature enables us to test the galaxy distribution in the direction of the Cen A group and to determine the statistical significance of the geometric alignment. Taking advantage of the fact that the two galaxy planes lie almost edge-on along the line of sight, the newly found 13 group members by Crnojevic et al. (2014, 2016) and our 16 new Cen A group candidates (M\"uller et al. 2016) can be assigned relative to the two planes. We use various statistical methods to test whether the distribution of galaxies follows a single normal distribution or shows evidence of bimodality as it has been reported earlier. We confirm that the data used for the Tully et al. (2015) study support the picture of a bimodal structure. However, when the new galaxy samples are included, the gap between the two galaxy planes is closing and the significance level of the bimodality is reduced. Instead, the plane that contains Cen A becomes more prominent. We found evidence that the galaxy system around Cen A is made up of only one plane of satellites. This plane is almost orthogonal to the dust plane of Cen A. Accurate distances to the new dwarf galaxies will be required to measure the precise 3D distribution of the galaxies around Cen A.
We present a detailed analysis of the pre-main-sequence (PMS) population of the young star cluster Westerlund~2 (Wd2), the central ionizing cluster of the HII region RCW 49, using data from a high resolution multi-band survey with the Hubble Space Telescope. The data were acquired with the Advanced Camera for Surveys in the F555W, F814W, and F658N filters and with the Wide Field Camera 3 in the F125W, F160W, and F128N filters. We find a mean age of the region of 1.04+-0.72 Myr. The combination of dereddened F555W and F814W photometry in combination with F658N photometry allows us to study and identify stars with H_alpha excess emission. With a careful selection of 240 bona-fide PMS H_alpha excess emitters we were able to determine their H_alpha luminosity, which has a mean value L(H_alpha)=1.67 x 10^{-31} erg s^{-1}. Using the PARSEC 1.2S isochrones to obtain the stellar parameters of the PMS stars we determined a mean mass accretion rate \dot M_acc=4.43 x 10^{-8} M_sun yr^{-1} per star. A careful analysis of the spatial dependence of the mass-accretion rate suggests that this rate is ~25% lower in center of the two density peaks of Wd2 in close proximity to the luminous OB stars, compared to the Wd2 average. This rate is higher with increasing distance from the OB stars, indicating that the PMS accretion disks are being rapidly destroyed by the far-ultra-violet radiation emitted by the OB population.
We apply a particular form of the inverse scattering theory to turbulent magnetic fluctuations in a plasma. In the present note we develop the theory, formulate the magnetic fluctuation problem in terms of its electrodynamic turbulent response function, and reduce it to the solution of a special form of the famous Gel$'$fand-Levitan-Marchenko equation of quantum mechanical scattering theory.
Timing results for the black-widow pulsar J2051-0827 are presented, using a 21-year dataset from four European Pulsar Timing Array telescopes and the Parkes radio telescope. This dataset, which is the longest published to date for a black-widow system, allows for an improved analysis that addresses previously unknown biases. While secular variations, as identified in previous analyses, are recovered, short-term variations are detected for the first time. Concurrently, a significant decrease of approx. 2.5x10-3 cm-3 pc in the dispersion measure associated with PSR J2051-0827 is measured for the first time and improvements are also made to estimates of the proper motion. Finally, PSR J2051-0827 is shown to have entered a relatively stable state suggesting the possibility of its eventual inclusion in pulsar timing arrays.
The ultraviolet background (UVB) emitted by quasars and galaxies governs the ionization and thermal state of the intergalactic medium (IGM), regulates the formation of high-redshift galaxies, and is thus a key quantity for modeling cosmic reionization. The vast majority of cosmological hydrodynamical simulations implement the UVB via a set of spatially uniform photoionization and photoheating rates derived from UVB synthesis models. We show that simulations using canonical UVB rates reionize, and perhaps more importantly, spuriously heat the IGM, much earlier z ~ 15 than they should. This problem arises because at z > 6, where observational constraints are non-existent, the UVB amplitude is far too high. We introduce a new methodology to remedy this issue, and generate self-consistent photoionization and photoheating rates to model any chosen reionization history. Following this approach, we run a suite of hydrodynamical simulations of different reionization scenarios, and explore the impact of the timing of reionization and its concomitant heat injection on the the thermal state of the IGM. We present a comprehensive study of the pressure smoothing scale of IGM gas, illustrating its dependence on the details of both hydrogen and helium reionization, and argue that it plays a fundamental role in interpreting Lyman-alpha forest statistics and the thermal evolution of the IGM. The premature IGM heating we have uncovered implies previous work has likely dramatically overestimated the impact of photoionization feedback on galaxy formation, which sets the minimum halo mass able to form stars at high redshifts. We make our new UVB photoionization and photoheating rates publicly available for use in future simulations.
Links to: arXiv, form interface, find, astro-ph, recent, 1607, contact, help (Access key information)