We measure the intrinsic Faber-Jackson (F-J) relation between velocity dispersion $\sigma$ and luminosity $L$ for massive, luminous red galaxies (LRGs) at redshift z~0.55. We achieve unprecedented precision by using a sample of 600,000 galaxies with spectra from the Baryon Oscillation Spectroscopic Survey (BOSS) of the third Sloan Digital Sky Survey (SDSS-III). We deconvolve the effects of photometric errors, limited spectroscopic signal-to-noise ratio, and red--blue galaxy confusion using a novel hierarchical Bayesian formalism that is generally applicable to any combination of photometric and spectroscopic observables. For a F-J relation of the form $L \propto \sigma^{\beta}$, we find $\beta = 7.8 \pm 1.1$ for $\sigma$ corrected to the effective radius. We find a very small intrinsic scatter of $s = 0.047 \pm 0.004$ in $\log_{10} \sigma$ at fixed $L$. Assuming plausible stellar population models, our measurements are consistent with no evolution in the parameters of the F-J relation over the range 0.5 < z < 0.7 covered by the sample. The steep F-J slope indicates that the scaling relations for the most massive LRGs are systematically different than the relations defined at lower masses, and the small scatter suggests that these galaxies more closely approximate a one-parameter family than their less massive counterparts. The curvature of the F-J relation has been observed previously in lower-mass and/or smaller galaxy samples; this new work provides a definitive measurement of the high-mass limit of the relation. Our results reinforce a picture in which the formation of LRGs is primarily driven by major dissipationless mergers.
We have observed two kinematically offset active galactic nuclei (AGN), whose ionised gas is at a different line-of-sight velocity to their host galaxies, with the SAMI integral field spectrograph (IFS). One of the galaxies shows gas kinematics very different to the stellar kinematics, indicating a recent merger or accretion event. We demonstrate that the star formation associated with this event was triggered within the last 100 Myr. The other galaxy shows simple disc rotation in both gas and stellar kinematics, aligned with each other, but in the central region has signatures of an outflow driven by the AGN. Other than the outflow, neither galaxy shows any discontinuity in the ionised gas kinematics at the galaxy's centre. We conclude that in these two cases there is no direct evidence of the AGN being in a supermassive black hole binary system. Our study demonstrates that selecting kinematically offset AGN from single-fibre spectroscopy provides, by definition, samples of kinematically peculiar objects, but IFS or other data are required to determine their true nature.
Molecular oxygen, O_2, has been the target of ground-based and space-borne searches for decades. Of the thousands of lines of sight surveyed, only those toward Rho Oph and Orion H_2 Peak 1 have yielded detections of any statistical significance. The detection of the O_2 N_J =3_3 -1_2 and 5_4 - 3_4 lines at 487.249 GHz and 773.840 GHz, respectively, toward Rho Ophiuchus has been attributed to a short-lived peak in the time-dependent, cold-cloud O_2 abundance, while the detection of the O_2 N_J =3_3 - 1_2, 5_4 - 3_4 lines, plus the 7_6 - 5_6 line at 1120.715 GHz, toward Orion has been ascribed to time-dependent preshock physical and chemical evolution and low-velocity (12 km/s) non-dissociative C-type shocks, both of which are fully shielded from far-ultraviolet (FUV) radiation, plus a postshock region that is exposed to a FUV field. We report a re-interpretation of the Orion O_2 detection based on new C-type shock models that fully incorporate the significant effects the presence of even a weak FUV field can have on the preshock gas, shock structure and postshock chemistry. In particular, we show that a family of solutions exists, depending on the FUV intensity, that reproduces both the observed O_2 intensities and O_2 line ratios. The solution in closest agreement with the shock parameters inferred for H_2 Peak 1 from other gas tracers assumes a 23 km/s shock impacting gas with a preshock density of 8x10^4 cm^-3 and G_0 =1, substantially different from that inferred for the fully-shielded shock case. As pointed out previously, the similarity between the LSR velocity of all three O_2 lines (~11 km/s) and recently measured H_2O 5_32 - 4_41 maser emission at 620.701 GHz toward H_2 Peak 1 suggests that the O_2 emission arises behind the same shocks responsible for the maser emission, though the O_2 emission is almost certainly more extended than the localized high density maser spots
The escape dynamics in a simple analytical gravitational model which describes the motion of stars in a Seyfert galaxy is investigated in detail. We conduct a thorough numerical analysis distinguishing between regular and chaotic orbits as well as between trapped and escaping orbits, considering only unbounded motion for several energy levels. In order to distinguish safely and with certainty between ordered and chaotic motion, we apply the Smaller ALingment Index (SALI) method. It is of particular interest to locate the escape basins through the openings around the collinear Lagrangian points $L_1$ and $L_2$ and relate them with the corresponding spatial distribution of the escape times of the orbits. Our exploration takes place both in the physical $(x,y)$ and in the phase $(x,\dot{x})$ space in order to elucidate the escape process as well as the overall orbital properties of the galactic system. Our numerical analysis reveals the strong dependence of the properties of the considered escape basins with the total orbital energy, with a remarkable presence of fractal basin boundaries along all the escape regimes. It was also observed, that for energy levels close to the critical escape energy the escape rates of orbits are large, while for much higher values of energy most of the orbits have low escape periods or they escape immediately to infinity. We also present evidence obtained through numerical simulations that our model can describe the formation and the evolution of the observed spiral structure in Seyfert galaxies. We hope our outcomes to be useful for a further understanding of the escape mechanism in galaxies with active nuclei.
Here we introduce moments of visibility function and discuss how those can be used to estimate the power spectrum of the turbulent velocity of external spiral galaxies. We perform numerical simulation to confirm the credibility of this method and found that for galaxies with lower inclination angles it works fine. This is the only method to estimate the power spectrum of the turbulent velocity fluctuation in the ISM of the external galaxies.
Context. The inner disc, linking the thin disc with the bulge, has been somehow neglected in the past because of intrinsic difficulties in its study, due, e.g., to crowding and high extinction. Open clusters located in the inner disc are among the best tracers of its chemistry at different ages and distances. Aims. We analyse the chemical patterns of four open clusters located within 7 kpc of the Galactic Centre and of field stars to infer the properties of the inner disc with the Gaia-ESO survey idr2/3 data release. Methods. We derive the parameters of the newly observed cluster, Berkeley 81, finding an age of about 1 Gyr and a Galactocentric distance of 5.4 kpc. We construct the chemical patterns of clusters and we compare them with those of field stars in the Solar neighbourhood and in the inner-disc samples. Results. Comparing the three populations we observe that inner-disc clusters and field stars are both, on average, enhanced in [O/Fe], [Mg/Fe] and [Si/Fe]. Using the idr2/3 results of M67, we estimate the non-local thermodynamic equilibrium (NLTE) effect on the abundances of Mg and Si in giant stars. After empirically correcting for NLTE effects, we note that NGC 6705 and Be 81 still have a high [{\alpha}/Fe]. Conclusions. The location of the four open clusters and of the field population reveals that the evolution of the metallicity [Fe/H] and of [alpha/Fe] can be explained within the framework of a simple chemical evolution model: both [Fe/H] and [{\alpha}/Fe] of Trumpler 20 and of NGC 4815 are in agreement with expectations from a simple chemical evolution model. On the other hand, NGC 6705, and at a lower level Berkeley 81, have higher [{\alpha}/Fe] than expected for their ages, location in the disc, and metallicity. These differences might originate from local enrichment processes as explained in the inhomogeneous evolution framework.
We present near-infrared emission-line flux distributions, excitation and kinematics, as well as stellar kinematics, of the inner 520x520 pc2$ of the Seyfert 2 galaxy NGC5929. The observations were performed with the Gemini's Near-Infrared Integral Field Spectrograph (NIFS) at a spatial resolution of 20 pc and spectral resolution of 40km/s in the J- and K-bands. The flux distributions of H2, [FeII], [PII] and recombination lines are extended over most of the field of view, with the highest intensity levels observed along PA=60/240deg, and well correlated with the radio emission. The H2 and [FeII] line emission are originated in thermal processes, mainly due to heating of the gas by X-rays from the central Active Galactic Nucleus (AGN). Contribution of shocks due to the radio jet is observed at locations co-spatial with the radio hotspots at 0.5" northeast and 0.6" southwest of the nucleus, as evidenced by the emission-line ratio and gas kinematics. The stellar kinematics shows rotation with an amplitude at 250pc from the nucleus of 200 km/s after corrected for the inferred inclination of 18.3deg. The stellar velocity dispersion obtained from the integrated K-band spectrum is sigma*=133+/-8 km/s, which implying on a mass for the supermassive black hole of M=5.2E7 Msun, using the M-sigma* relation. The gas kinematics present three components: (1) gas in the plane of the galaxy in counter-rotation relative to the stars; (2) an outflow perpendicular to the radio jet that seems to be due to an equatorial AGN outflow; (3) turbulence of the gas observed in association with the radio hot spots, supporting an interaction of the radio jet with the gas of the disk. We estimated the mass of ionized and warm molecular gas of ~1.3E6 Msun and ~470 Msun, respectively.
We have mapped the superwind/halo region of the nearby starburst galaxy M82 in the mid-infrared with $Spitzer-IRS$. The spectral regions covered include the H$_2 S(1)-S(3)$, [NeII], [NeIII] emission lines and PAH features. We estimate the total warm H$_2$ mass and the kinetic energy of the outflowing warm molecular gas to be between $M_{warm}\sim5-17\times10^6$ M$_{\odot}$ and $E_{K}\sim6-20\times10^{53}$ erg. Using the ratios of the 6.2, 7.7 and 11.3 micron PAH features in the IRS spectra, we are able to estimate the average size and ionization state of the small grains in the superwind. There are large variations in the PAH flux ratios throughout the outflow. The 11.3/7.7 and the 6.2/7.7 PAH ratios both vary by more than a factor of five across the wind region. The Northern part of the wind has a significant population of PAH's with smaller 6.2/7.7 ratios than either the starburst disk or the Southern wind, indicating that on average, PAH emitters are larger and more ionized. The warm molecular gas to PAH flux ratios (H$_2/PAH$) are enhanced in the outflow by factors of 10-100 as compared to the starburst disk. This enhancement in the H$_2/PAH$ ratio does not seem to follow the ionization of the atomic gas (as measured with the [NeIII]/[NeII] line flux ratio) in the outflow. This suggests that much of the warm H$_2$ in the outflow is excited by shocks. The observed H$_2$ line intensities can be reproduced with low velocity shocks ($v < 40$ km s$^{-1}$) driven into moderately dense molecular gas ($10^2 <n_H < 10^4$ cm$^{-3}$) entrained in the outflow.
Elemental abundance patterns in the Galactic disk constrain theories of the formation and evolution of the Milky Way. HII region abundances are the result of billions of years of chemical evolution. We made radio recombination line and continuum measurements of 21 HII regions located between Galactic azimuth Az = 90-130 degree, a previously unexplored region. We derive the plasma electron temperatures using the line-to-continuum ratios and use them as proxies for the nebular [O/H] abundances, because in thermal equilibrium the abundance of the coolants (O, N, and other heavy elements) in the ionized gas sets the electron temperature, with high abundances producing low temperatures. Combining these data with our previous work produces a sample of 90 HII regions with high quality electron temperature determinations. We derive kinematic distances in a self-consistent way for the entire sample. The radial gradient in [O/H] is -0.082 +/- 0.014 dex/kpc for Az = 90-130 degree, about a factor of two higher than the average value between Az = 0-60 degree. Monte Carlo simulations show that the azimuthal structure we reported for Az = 0-60 degree is not significant because kinematic distance uncertainties can be as high as 50% in this region. Nonetheless, the flatter radial gradients between Az = 0-60 degree compared with Az = 90-130 degree, are significant within the uncertainty. We suggest that this may be due to radial mixing from the Galactic Bar whose major axis is aligned toward Az ~30 degree.
In our recent paper, we argued that the broad hydrogen lines, as well as the low-ionization lines in quasars, are significantly contributed by the Cerenkov quasi-line emission of the fast electrons in the dense clouds/filaments/sheets ($N_{\rm H}\geq 10^{14}~{\rm cm^{-3}}$); whereas this line-like radiation mechanism is invalid for producing the high ionization lines. Therefore, the observed broad hydrogen lines or low-ionization lines should be blended by both the real line emission via the bound-bound transition in atoms/ions and the Cerenkov quasi-line emission. In this paper, we provide an evidence quantitatively supporting above conclusions. Until now the observed different redshifts of different broad hydrogen lines is still a problem at issue. The Cerenkov line-like radiation mechanism provides a plausible resolution for this difficulty: it is the `Cerenkov line redshift', which is different from line to line, causes the peculiar redshift-differences among Ly$\alpha$, H$\alpha$ and H$\beta$ lines. The good fitting to the observed redshifts of quasars confirms the existence of Cerenkov component in the broad hydrogen lines. Furthermore, we conclude that, in the blended Ly$\alpha$ line, the line-intensity of the Cerenkov component approximately equals that of the accompanying `normal line' (an approximate equipartition of intensity between the two components in the broad Ly$\alpha$ line). This result quantitatively illustrates the importance of the Cerenkov component in the broad lines of quasars.
Galaxy proto-clusters at z >~ 2 provide a direct probe of the rapid mass assembly and galaxy growth of present day massive clusters. Because of the need of precise galaxy redshifts for density mapping and the prevalence of star formation before quenching, nearly all the proto-clusters known to date were confirmed by spectroscopy of galaxies with strong emission lines. Therefore, large emission-line galaxy surveys provide an efficient way to identify proto-clusters directly. Here we report the discovery of a large-scale structure at z = 2.44 in the HETDEX Pilot Survey. On a scale of a few tens of Mpc comoving, this structure shows a complex overdensity of Lya emitters (LAE), which coincides with broad-band selected galaxies in the COSMOS/UltraVISTA photometric and zCOSMOS spectroscopic catalogs, as well as overdensities of intergalactic gas revealed in the Lya absorption maps of Lee et al. (2014). We construct mock LAE catalogs to predict the cosmic evolution of this structure. We find that such an overdensity should have already broken away from the Hubble flow, and part of the structure will collapse to form a galaxy cluster with 10^14.5 +- 0.4 M_sun by z = 0. The structure contains a higher median stellar mass of broad-band selected galaxies, a boost of extended Lya nebulae, and a marginal excess of active galactic nuclei relative to the field, supporting a scenario of accelerated galaxy evolution in cluster progenitors. Based on the correlation between galaxy overdensity and the z = 0 descendant halo mass calibrated in the simulation, we predict that several hundred 1.9 < z < 3.5 proto-clusters with z = 0 mass of > 10^14.5 M_sun will be discovered in the 8.5 Gpc^3 of space surveyed by the Hobby Eberly Telescope Dark Energy Experiment.
The barycenter of a massive black hole binary will lie outside the event horizon of the primary black hole for modest values of mass ratio and binary separation. Analagous to radial velocity shifts in stellar emission lines caused by the tug of planets, the radial velocity of the primary black hole around the barycenter can leave a tell-tale oscillation in the broad component of Fe K$\alpha$ emission from accreting gas. Near-future X-ray telescopes such as Astro-H and Athena will have the energy resolution ($\delta E/E \lesssim 10^{-3}$) to search nearby active galactic nuclei (AGN) for the presence of binaries with mass ratios $q \gtrsim 0.01$, separated by several hundred gravitational radii. The general-relativistic and Lense-Thirring precession of the periapse of the secondary orbit imprints a detectable modulation on the oscillations. The lowest mass binaries in AGN will oscillate many times within typical X-ray exposures, leading to a broadening of the line wings and an over-estimate of black hole spin in these sources. Detection of periodic oscillations in the AGN line centroid energy will reveal a massive black hole binary close to merger and will provide an early warning of gravitational radiation emission.
Links to: arXiv, form interface, find, astro-ph, recent, 1505, contact, help (Access key information)
Thermal mid-infrared emission of quasars requires an obscuring structure that can be modeled as a magneto-hydrodynamic wind in which radiation pressure on dust shapes the outflow. We have taken the dusty wind models presented by Keating and collaborators that generated quasar mid-infrared spectral energy distributions (SEDs), and explored their properties (such as geometry, opening angle, and ionic column densities) as a function of Eddington ratio and X-ray weakness. In addition, we present new models with a range of magnetic field strengths and column densities of the dust-free shielding gas interior to the dusty wind. We find this family of models -- with input parameters tuned to accurately match the observed mid-IR power in quasar SEDs -- provides reasonable values of the Type 1 fraction of quasars and the column densities of warm absorber gas, though it does not explain a purely luminosity-dependent covering fraction for either. Furthermore, we provide predictions of the cumulative distribution of E(B-V) values of quasars from extinction by the wind and the shape of the wind as imaged in the mid-infrared. Within the framework of this model, we predict that the strength of the near-infrared bump from hot dust emission will be correlated primarily with L/L_Edd rather than luminosity alone, with scatter induced by the distribution of magnetic field strengths. The empirical successes and shortcomings of these models warrant further investigations into the composition and behaviour of dust and the nature of magnetic fields in the vicinity of actively accreting supermassive black holes.
We extract the resonant orbits from an N-body bar that is a good representation of the Milky Way, using the method recently introduced by Molloy et al. (2015). By decomposing the bar into its constituent orbit families, we show that they are intimately connected to the boxy-peanut shape of the density. We highlight the imprint due solely to resonant orbits on the kinematic landscape towards the Galactic centre. The resonant orbits are shown to have distinct kinematic features and may be used to explain the cold velocity peak seen in the APOGEE commissioning data (Nidever et al. 2012). We show that high velocity peaks are a natural consequence of the motions of stars in the 2:1 orbit family. The locations of the peaks vary with bar angle and, with the tacit assumption that the observed peaks are due to the 2:1 family, we find that the locations of the high velocity peaks correspond to bar angles in the range 10 < theta_bar < 25 (deg). However, some important questions about the nature of the peaks remain, such as their apparent absence in other surveys of the Bulge and the deviations from symmetry between equivalent fields in the north and south.
With the use of a background Milky-Way-like potential model, we performed stellar orbital and magnetohydrodynamic (MHD) simulations. As a first experiment, we studied the gaseous response to a bisymmetric spiral arm potential: the widely employed cosine potential model and a self-gravitating tridimensional density distribution based model called PERLAS. Important differences are noticeable in these simulations, while the simplified cosine potential produces two spiral arms for all cases, the more realistic density based model produces a response of four spiral arms on the gaseous disk, except for weak arms -i.e. close to the linear regime- where a two-armed structure is formed. In order to compare the stellar and gas response to the spiral arms, we have also included a detailed periodic orbit study and explored different structural parameters within observational uncertainties. The four armed response has been explained as the result of ultra harmonic resonances, or as shocks with the massive bisymmetric spiral structure, among other. From the results of this work, and comparing the stellar and gaseous responses, we tracked down an alternative explanation to the formation of branches, based only on the orbital response to a self-gravitating spiral arms model. The presence of features such as branches, might be an indication of transiency of the arms.
Using new integral field observations of 106 galaxies in three nearby clusters we investigate how the intrinsic scatter of the Fundamental Plane depends on the way in which the velocity dispersion and effective radius are measured. Our spatially resolved spectroscopy, combined with a cluster sample with negligible relative distance errors allows us to derive a Fundamental Plane with minimal systematic uncertainties. From the apertures we tested, we find that velocity dispersions measured within a circular aperture with radius equal to one effective radius minimises the intrinsic scatter of the Fundamental Plane. Using simple yet powerful Jeans dynamical models we determine dynamical masses for our galaxies. Replacing luminosity in the Fundamental Plane with dynamical mass, we demonstrate that the resulting Mass Plane has further reduced scatter, consistent with zero intrinsic scatter. Using these dynamical models we also find evidence for a possibly non-linear relationship between dynamical mass-to-light ratio and velocity dispersion.
With the use of a detailed Milky Way nonaxisymmetric potential, observationally and dynamically constrained, the e?ects of the bar and the spiral arms in the Galaxy are studied in the disc and in the stellar halo. Especially the trapping of stars in the disc and Galactic halo by resonances on the Galactic plane, induced by the Galactic bar, has been analysed in detail. To this purpose, a new method is presented to delineate the trapping regions using empirical diagrams of some orbital properties obtained in the Galactic potential. In these diagrams we plot in the inertial Galactic frame a characteristic orbital energy versus a characteristic orbital angular momentum, or versus the orbital Jacobi constant in the reference frame of the bar, when this is the only nonaxisymmetric component in the Galactic potential. With these diagrams some trapping regions are obtained in the disc and halo using a sample of disc stars and halo stars in the solar neighbourhood. We compute several families of periodic orbits on the Galactic plane, some associated with this resonant trapping. In particular, we ?nd that the trapping e?ect of these resonances on the Galactic plane can extend several kpc from this plane, trapping stars in the Galactic halo. The purpose of our analysis is to investigate if the trapping regions contain some known moving groups in our Galaxy. We have applied our method to the Kapteyn group, a moving group in the halo, and we have found that this group appears not to be associated with a particular resonance on the Galactic plane.
Over the last decade there has been immense progress in the follow-up of short and long GRBs, resulting in a significant rise in the detection rate of X-ray and optical afterglows, in the determination of GRB redshifts, and of the identification of the underlying host galaxies. Nevertheless, our theoretical understanding on the progenitors and central engines powering these vast explosions is lagging behind, and a newly identified class of `ultra-long' GRBs has fuelled speculation on the existence of a new channel of GRB formation. In this paper we present high signal-to-noise X-shooter observations of the host galaxy of GRB130925A, which is the fourth unambiguously identified ultra-long GRB, with prompt gamma-ray emission detected for ~20ks. The GRB line of sight was close to the host galaxy nucleus, and our spectroscopic observations cover both this region along the bulge/disk of the galaxy, in addition to a bright star-forming region within the outskirts of the galaxy. From our broad wavelength coverage we obtain accurate metallicity and dust-extinction measurements at both the galaxy nucleus, and outer star-forming region, and measure a super-solar metallicity at both locations, placing this galaxy within the 10-20% most metal-rich GRB host galaxies. Such a high metal enrichment has implications on the progenitor models of both long and ultra-long GRBs, although the edge-on orientation of the host galaxy does not allow us to rule out a large metallicity variation along our line of sight. The spatially resolved spectroscopic data presented in this paper offer important insight into variations in the metal and dust abundance within GRB host galaxies. They also illustrate the need for IFU observations on a larger sample of GRB host galaxies at varies metallicities to provide a more quantitative view on the relation between the GRB circumburst and the galaxy-whole properties.
We use $N$-body/smoothed particle hydrodynamics simulations of encounters between an early-type galaxy (ETG) and a late-type galaxy (LTG) to study the effects of hot halo gas on the evolution for a case with the mass ratio of the ETG to LTG of 2:1 and the closest approach distance of $\sim$100 kpc. We find that the dynamics of the cold disk gas in the tidal bridge and the amount of the newly formed stars depend strongly on the existence of a gas halo. In the run of interacting galaxies not having a hot gas halo, the gas and stars accreted into the ETG do not include newly formed stars. However, in the run using the ETG with a gas halo and the LTG without a gas halo, a shock forms along the disk gas tidal bridge and induces star formation near the closest approach. The shock front is parallel to a channel along which the cold gas flows toward the center of the ETG. As a result, the ETG can accrete star-forming cold gas and newly born stars at and near its center. When both galaxies have hot gas halos, a shock is formed between the two gas halos somewhat before the closest approach. The shock hinders the growth of the cold gas bridge to the ETG and also ionizes it. Only some of the disk stars transfer through the stellar bridge. We conclude that the hot halo gas can give significant hydrodynamic effects during distant encounters.
The wavelength displacement of the Diffuse Interstellar Bands at 4502, 5705, 5780, 6284, and 7224 \AA\ with respect to the well known, narrow atomic/molecular interstellar lines (of Ca{\sc ii} and Na{\sc i}) have been measured in the spectra of the 2 Orion Trapezium stars HD 37022 and HD 37020, using the HARPS\textendash N spectrograph, fed with the 3.5 m Telescopio Nazionale Galileo, and the BOES spectrograph, fed with the 1.8m Korean telescope. The red shift is $\sim$25 km/s for all these DIBs. We discuss the various possible origins of this very peculiar wavelength shift in the light of the particular physical conditions in the Orion Trapezium. The above mentioned shift is seemingly absent in the DIBs at 6196 and 6993 \AA.
We use self-consistent N-body simulations to study the interaction of massive classical bulges (ClBs) with a bar that forms self-consistently in the disc. We show that the ClB gains significant angular momentum from the bar which, surprisingly, scales approximately linearly with the ClB mass. It is also tightly correlated with the ratio of the bulge size to the bar size. Most of the angular momentum gain occurs via low-order resonances, particularly 5:2 resonant orbits. A density wake forms in the ClB which corotates and aligns with the bar at the end of the evolution. The spin-up process creates a characteristic linear rotation profile and mild tangential anisotropy in the ClB. The induced rotation is small in the centre but significant beyond ~bulge half mass radii, where it leads to mass-weighted V/\sigma ~ 0.2, and reaches a local V/\sigma_in ~ 0.5 at around the scale of the bar. In all models a box/peanut bulge also forms suggesting that composite bulges may be common.
Interstellar dust in galaxies can be traced either through its extinction effects on the star light, or through its thermal emission at infrared wavelengths. Recent radiative transfer studies of several nearby edge-on galaxies have found an apparent inconsistency in the dust energy balance: the radiative transfer models that successfully explain the optical extinction underestimate the observed fluxes by an average factor of three. We investigate the dust energy balance for IC4225 and NGC5166, two edge-on spiral galaxies observed by the Herschel Space Observatory in the frame of the H-ATLAS survey. We start from models which were constrained from optical data and extend them to construct the entire spectral energy distribution of our galaxies. These predicted values are subsequently compared to the observed far-infrared fluxes. We find that including a young stellar population in the modelling is necessary as it plays a non-negligible part in the heating of the dust grains. While the modelling approach for both galaxies is nearly identical, we find two very different results. As is often seen in other edge-on spiral galaxies, the far-infrared emission of our radiative transfer model of IC4225 underestimates the observed fluxes by a factor of about three. For NGC5166 on the other hand, we find that both the predicted spectral energy distribution as well as the simulated images match the observations particularly well. We explore possible reasons for this difference and conclude that it is unlikely that one single mechanism is the cause of the dust energy balance problem in spiral galaxies. We discuss the different approaches that can be considered in order to get a conclusive answer on the origin this discrepancy.
Deep HI observations of the outer parts of disc galaxies demonstrate the frequent presence of extended, well-developed spiral arms far beyond the optical radius. To understand the nature and the origin of such outer spiral structure, we investigate the propagation in the outer gaseous disc of large-scale spiral waves excited in the bright optical disc. Using hydrodynamical simulations, we show that non-axisymmetric density waves, penetrating in the gas through the outer Lindblad resonance, can exhibit relatively regular spiral structures outside the bright optical stellar disc. For low-amplitude structures, the results of numerical simulations match the predictions of a simple WKB linear theory. The amplitude of spiral structure increases rapidly with radius. Beyond $\approx 2$ optical radii, spirals become nonlinear (the linear theory becomes quantitatively and qualitatively inadequate) and unstable to Kelvin-Helmholtz instability. In numerical simulations, in models for which gas is available very far out, spiral arms can extend out to 25 disc scale-lengths. A comparison between the properties of the models we have investigated and the observed properties of individual galaxies may shed light into the problem of the amount and distribution of dark matter in the outer halo.
We present a comparison of 14 galaxy formation models: 12 different semi-analytical models and 2 halo-occupation distribution models for galaxy formation based upon the same cosmological simulation and merger tree information derived from it. The participating codes have proven to be very successful in their own right but they have all been calibrated independently using various observational data sets, stellar models, and merger trees. In this paper we apply them without recalibration and this leads to a wide variety of predictions for the stellar mass function, specific star formation rates, stellar-to- halo mass ratios, and the abundance of orphan galaxies. The scatter is much larger than seen in previous comparison studies primarily because the codes have been used outside of their native environment within which they are well tested and calibrated. The purpose of the `nIFTy comparison of galaxy formation models' is to bring together as many different galaxy formation modellers as possible and to investigate a common approach to model calibration. This paper provides a unified description for all participating models and presents the initial, uncalibrated comparison as a baseline for our future studies where we will develop a common calibration framework and address the extent to which that reduces the scatter in the model predictions seen here.
We describe an overall picture of galactic-scale star formation. Recent high-resolution magneto-hydrodynamical simulations of two-fluid dynamics with cooling/heating and thermal conduction have shown that the formation of molecular clouds requires multiple episodes of supersonic compression. This finding enables us to create a scenario in which molecular clouds form in interacting shells or bubbles on a galactic scale. First we estimate the ensemble-averaged growth rate of molecular clouds over a timescale larger than a million years. Next we perform radiation hydrodynamics simulations to evaluate the destruction rate of magnetized molecular clouds by the stellar FUV radiation. We also investigate the resultant star formation efficiency within a cloud which amounts to a low value (a few percent) if we adopt the power-law exponent -2.5 for the mass distribution of stars in the cloud. We finally describe the time evolution of the mass function of molecular clouds over a long timescale (>1Myr) and discuss the steady state exponent of the power-law slope in various environments.
Galaxy evolution scenarios predict that the feedback of star formation and nuclear activity (AGN) can drive the transformation of gas-rich spiral mergers into ULIRGs, and, eventually, lead to the build-up of QSO/elliptical hosts. We study the role that star formation and AGN feedback have in launching and maintaining the molecular outflows in two starburst-dominated advanced mergers, NGC1614 and IRAS17208-0014, by analyzing the distribution and kinematics of their molecular gas reservoirs. We have used the PdBI array to image with high spatial resolution (0.5"-1.2") the CO(1-0) and CO(2-1) line emissions in NGC1614 and IRAS17208-0014, respectively. The velocity fields of the gas are analyzed and modeled to find the evidence of molecular outflows in these sources and characterize the mass, momentum and energy of these components. While most (>95%) of the CO emission stems from spatially-resolved (~2-3kpc-diameter) rotating disks, we also detect in both mergers the emission from high-velocity line wings that extend up to +-500-700km/s, well beyond the estimated virial range associated with rotation and turbulence. The kinematic major axis of the line wing emission is tilted by ~90deg in NGC1614 and by ~180deg in IRAS17208-0014 relative to their respective rotating disk major axes. These results can be explained by the existence of non-coplanar molecular outflows in both systems. In stark contrast with NGC1614, where star formation alone can drive its molecular outflow, the mass, energy and momentum budget requirements of the molecular outflow in IRAS17208-0014 can be best accounted for by the existence of a so far undetected (hidden) AGN of L_AGN~7x10^11 L_sun. The geometry of the molecular outflow in IRAS17208-0014 suggests that the outflow is launched by a non-coplanar disk that may be associated with a buried AGN in the western nucleus.
Observations suggest that high-redshift galaxies are either very dusty or essentially dust free. The evolution from one regime to the other must also be very fast, since evolved and dusty galaxies show up at redshifts corresponding to a Universe which is only about 500 Myr old. In the present paper models which predicts the existence of an apparent dichotomy between dusty and dust-free galaxies at high redshift are considered. Galaxies become dusty as soon as they reach an evolved state and the transition is very rapid. A special case suggests that while stellar dust production is overall relatively insignificant -- contrary to what has been argued recently -- it can at the same time be consistent with efficient dust production in supernovae in the local Universe. Special attention will be given to the recent discovery of a dusty normal galaxy (A1689-zD1) at a very high redshift z = 7.5 +/- 0.2.
MaNGA (Mapping Nearby Galaxies at Apache Point Observatory) is an integral-field spectroscopic survey of 10,000 nearby galaxies that is one of three core programs in the fourth-generation Sloan Digital Sky Survey (SDSS-IV). MaNGA's 17 pluggable optical fiber-bundle integral field units (IFUs) are deployed across a 3 deg field, they yield spectral coverage 3600-10,300 Ang at a typical resolution R ~ 2000, and sample the sky with 2" diameter fiber apertures with a total bundle fill factor of 56%. Observing over such a large field and range of wavelengths is particularly challenging for obtaining uniform and integral spatial coverage and resolution at all wavelengths and across each entire fiber array. Data quality is affected by the IFU construction technique, chromatic and field differential refraction, the adopted dithering strategy, and many other effects. We use numerical simulations to constrain the hardware design and observing strategy for the survey with the aim of ensuring consistent data quality that meets the survey science requirements while permitting maximum observational flexibility. We find that MaNGA science goals are best achieved with IFUs composed of a regular hexagonal grid of optical fibers with rms displacement of 5 microns or less from their nominal packing position, this goal is met by the MaNGA hardware, which achieves 3 microns rms fiber placement. We further show that MaNGA observations are best obtained in sets of three 15-minute exposures dithered along the vertices of a 1.44 arcsec equilateral triangle, these sets form the minimum observational unit, and are repeated as needed to achieve a combined signal-to-noise ratio of 5 per Angstrom per fiber in the r-band continuum at a surface brightness of 23 AB/arcsec^2. (abbrev.)
We use the Arepo moving mesh code to simulate the evolution of molecular clouds exposed to a harsh environment similar to that found in the galactic center (GC), in an effort to understand why the star formation efficiency (SFE) of clouds in this environment is so small. Our simulations include a simplified treatment of time-dependent chemistry and account for the highly non-isothermal nature of the gas and the dust. We model clouds with a total mass of 1.3x10^5 M_{sun} and explore the effects of varying the mean cloud density and the virial parameter, alpha = E_{kin}/|E_{pot}|. We vary the latter from alpha = 0.5 to alpha = 8.0, and so many of the clouds that we simulate are gravitationally unbound. We expose our model clouds to an interstellar radiation field (ISRF) and cosmic ray flux (CRF) that are both a factor of 1000 higher than the values found in the solar neighbourhood. As a reference, we also run simulations with local solar neighbourhood values of the ISRF and the CRF in order to better constrain the effects of the extreme conditions in the GC on the SFE. Despite the harsh environment and the large turbulent velocity dispersions adopted, we find that all of the simulated clouds form stars within less than a gravitational free-fall time. Increasing the virial parameter from alpha = 0.5 to alpha = 8.0 decreases the SFE by a factor ~4-10, while increasing the ISRF/CRF by a factor of 1000 decreases the SFE again by a factor ~2-6. However, even in our most unbound clouds, the SFE remains higher than that inferred for real GC clouds. We therefore conclude that high levels of turbulence and strong external heating are not enough by themselves to explain the low SFE at the center of the Galaxy.
We present results from 43 GHz (VLBA, six epochs from 2012.2 to 2013.2) and 86 GHz (GMVA, one epoch in 2012.4) observations toward the basis of the jet in the TeV Blazar Mrk 501. The 43-GHz data analysis reveals a new feature located northeast of the radio core, with a flux density of several tens of mJy, perpendicularly to the jet axis. The 86-GHz image shows the jet feature located 0.75 mas southeast of the radio core, which is consistent with the previous result. The location of Gaussian model for 0.75 mas feature does not coincide with those for the jets in the 43-GHz image, however, a distribution of emission is found. We also discuss the spectral indices of the core, the northeast feature, and the jet feature between 43 GHz and 86 GHz, which show flat-to-steep, steep, and flat-to-invert, respectively.
Numerical methods to improve the treatment of magnetic fields in smoothed field magnetohydrodynamics (SPMHD) are developed and tested. Chapter 2 is a review of SPMHD. In Chapter 3, a mixed hyperbolic/parabolic scheme is developed which cleans divergence error from the magnetic field. Average divergence error is an order of magnitude lower for all test cases considered, and allows for the stable simulation of the gravitational collapse of magnetised molecular cloud cores. The effectiveness of the cleaning may be improved by explicitly increasing the hyperbolic wave speed or by cycling the cleaning equations between timesteps. In the latter, it is possible to achieve DivB=0. Chapter 4 develops a switch to reduce dissipation of the magnetic field from artificial resistivity. Compared to the existing switch in the literature, this leads to sharper shock profiles in shocktube tests, lower overall dissipation of magnetic energy, and importantly, is able to capture magnetic shocks in the highly super-Alfvenic regime. Chapter 5 compares these numerical methods against grid-based MHD methods (using the Flash code) in simulations of the small-scale dynamo amplification of a magnetic field in driven, isothermal, supersonic turbulence. Both codes exponentially amplify the magnetic energy at a constant rate, though SPMHD shows a resolution dependence that arises from the scaling of the numerical dissipation terms. The time-averaged saturated magnetic spectra have similar shape, and both codes have PDFs of magnetic field strength that are log-normal, which become lopsided as the magnetic field saturates. We conclude that SPMHD is able to reliably simulate the small-scale dynamo amplification of magnetic fields. Chapter 6 concludes the thesis and presents some preliminary work demonstrating that SPMHD can activate the magneto-rotational instability in 2D shearing box tests.
Links to: arXiv, form interface, find, astro-ph, recent, 1505, contact, help (Access key information)
In light of recent findings from the kinematic morphology-density relation, we investigate whether the same trends exist in the original morphology density relation, using the same data as Dressler. In addition to Dressler's canonical relations, we find that further refinements are possible when considering the average local projected density of galaxies in a cluster. Firstly, the distribution of ellipticals in a cluster depends on the relative local density of galaxies in that cluster: equivalent rises in the elliptical fraction occur at higher local densities for clusters with higher average local densities. This is not true for the late-type fraction, where the variation with local density within a cluster is independent of the average local density of galaxies in that cluster, and is as Dressler originally found. Furthermore, the overall ratio of ellipticals to early-types in a cluster does not depend on the average density of galaxies in that cluster (unlike the ratio of lenticulars to disk systems), and is fixed at around 30%. In the paradigm of fast and slow rotators, we show that such an elliptical fraction in the early-type population is consistent with a slow rotator fraction of 15% in the early-type population, using the statistics of the ATLAS3D survey. We also find the scatter in the overall ratio of ellipticals to early-types is greatest for clusters with lower average densities, such that clusters with the highest elliptical fractions have the lowest average local densities. Finally, we show that average local projected density correlates well with global projected density, but the latter has difficulty in accurately characterising the density of irregular cluster morphologies.
Stars with a different vertical motion relative to the galactic disk have a different average acceleration. According to Modified Newtonian Dynamics (MOND) theories they should therefore have a different average orbital velocity while revolving around the Milky Way. We show that this property can be used to constrain MOND theories by studying stars in the local neighborhood. With the Hipparcos dataset we can only place marginal constraints. However, the forthcoming GAIA catalogue with its significantly fainter cutoff should allow placing a stringent constraint. The method cannot be used to prove MOND, since halo stars can contribute a similar signal which would be hard to discern.
In our recent investigation of the Oosterhoff dichotomy in the multiple population paradigm (Jang et al. 2014), we have suggested that the RR Lyrae variables in the Oosterhoff groups I, II, and III globular clusters (GCs) are produced mostly by the first, second, and third generation stars (G1, G2, and G3), respectively. Here we show, for the first time, that the observed dichotomies in the inner and outer halo GCs can be naturally reproduced when these models are extended to all metallicity regimes, while maintaining reasonable agreements in the horizontal-branch type versus [Fe/H] correlations. In order to achieve this, however, specific star formation histories are required for the inner and outer halos. In the inner halo GCs, the star formation commenced and ceased earlier with relatively short formation timescale between the subpopulations (~0.5 Gyr), while in the outer halo, the formation of G1 was delayed by ~0.8 Gyr with more extended timescale between G1 and G2 (~1.4 Gyr). This is consistent with the dual origin of the Milky Way halo. Despite the difference in detail, our models show that the Oosterhoff period groups observed in both outer and inner halo GCs are all manifestations of the "population-shift" effect within the instability strip, for which the origin can be traced back to the two or three discrete episodes of star formation in GCs.
We introduce a dust model for cosmological simulations implemented in the moving-mesh code AREPO and present a suite of cosmological hydrodynamical zoom-in simulations to study dust formation within galactic haloes. Our model accounts for the stellar production of dust, accretion of gas-phase metals onto existing grains, destruction of dust through local supernova activity, and dust driven by winds from star-forming regions. We find that accurate stellar and active galactic nuclei feedback is needed to reproduce the observed dust-metallicity relation and that dust growth largely dominates dust destruction. Our simulations predict a dust content of the interstellar medium which is consistent with observed scaling relations at $z = 0$, including scalings between dust-to-gas ratio and metallicity, dust mass and gas mass, dust-to-gas ratio and stellar mass, and dust-to-stellar mass ratio and gas fraction. We find that roughly two-thirds of dust at $z = 0$ originated from Type II supernovae, with the contribution from asymptotic giant branch stars below 20 per cent for $z \gtrsim 5$. While our suite of Milky Way-sized galaxies forms dust in good agreement with a number of key observables, it predicts a high dust-to-metal ratio in the circumgalactic medium, which motivates a more realistic treatment of thermal sputtering of grains and dust cooling channels.
The outer halo regions of early-type galaxies carry key information about their past accretion history. However, spectroscopically probing the stellar component at such galactocentric radii is still challenging. Using Keck/DEIMOS, we have been able to measure the metallicities of the stellar and globular cluster components in 12 early-type galaxies out to more than $10~\rm{R_{e}}$. We find similar metallicity gradients for the metal-poor and metal-rich globular cluster subpopulations, suggesting a common formation process for the two subpopulations. This is in conflict with most current theoretical predictions, where the metal-poor globular clusters are thought to be purely accreted and metal-rich globular clusters mostly formed in-situ. Moreover, we find that the globular cluster metallicity gradients show a trend with galaxy mass, being steeper in lower-mass galaxies than in higher-mass galaxies. This is similar to what we find for the outermost galaxy stars and suggests a more active accretion history, with a larger role played by major mergers, in the most massive galaxies. This conclusion is qualitatively consistent with expectations from two-phase galaxy assembly models. Finally, we find that the small difference in metallicity between galaxy stars and metal-rich globular clusters at $1~\rm{R_{e}}$ may correlate with galaxy mass. The origin of this difference is not currently clear.
It has been recently found that the characteristic photometric parameters of antitruncated discs in S0 galaxies follow tight scaling relations. We investigate if similar scaling relations are satisfied by galaxies of other morphological types. We have analysed the trends in several photometric planes relating the characteristic surface brightness and scalelengths of the breaks and the inner and outer discs of local antitruncated S0-Scd galaxies, using published data and fits performed to the surface brightness profiles of two samples of Type-III galaxies in the R and Spitzer 3.6 microns bands. We have performed linear fits to the correlations followed by different galaxy types in each plane, as well as several statistical tests to determine their significance. We have found that: 1) the antitruncated discs of all galaxy types from Sa to Scd obey tight scaling relations both in R and 3.6 microns, as observed in S0s; 2) the majority of these correlations are significant accounting for the numbers of the available data samples; 3) the trends are clearly linear when the characteristic scalelengths are plotted on a logarithmic scale; and 4) the correlations relating the characteristic surface brightnesses of the inner and outer discs and the breaks with the various characteristic scalelengths significantly improve when the latter are normalized to the optical radius of the galaxy. The results suggest that the scaling relations of Type-III discs are independent of the morphological type and the presence (or absence) of bars within the observational uncertainties of the available datasets, although larger and deeper samples are required to confirm this. The tight structural coupling implied by these scaling relations impose strong constraints on the mechanisms proposed for explaining the formation of antitruncated stellar discs in the galaxies across the whole Hubble Sequence (Abridged).
We investigate the ability of current implementations of galaxy group finders to recover colour-dependent halo occupation statistics. To test the fidelity of group catalogue inferred statistics, we run three different group finders used in the literature over a mock that includes galaxy colours in a realistic manner. Overall, the resulting mock group catalogues are remarkably similar, and most colour-dependent statistics are recovered with reasonable accuracy. However, it is also clear that certain systematic errors arise as a consequence of correlated errors in group membership determination, central/satellite designation, and halo mass assignment. We introduce a new statistic, the halo transition probability (HTP), which captures the combined impact of all these errors. As a rule of thumb, errors tend to equalize the properties of distinct galaxy populations (i.e. red vs. blue galaxies or centrals vs. satellites), and to result in inferred occupation statistics that are more accurate for red galaxies than for blue galaxies. A statistic that is particularly poorly recovered from the group catalogues is the red fraction of central galaxies as function of halo mass. Group finders do a good job in recovering galactic conformity, but also have a tendency to introduce weak conformity when none is present. We conclude that proper inference of colour-dependent statistics from group catalogues is best achieved using forward modelling (i.e., running group finders over mock data), or by implementing a correction scheme based on the HTP, as long as the latter is not too strongly model-dependent.
Reverberation mapping is now a well-established technique for investigating spatially-unresolved structures in the nuclei of distant galaxies with actively-accreting supermassive black holes. Structural parameters for the broad emission-line region, with angular sizes of microarcseconds, can be constrained through the substitution of time resolution for spatial resolution. Many reverberation experiments over the last 30 years have led to a practical understanding of the requirements necessary for a successful program. With reverberation measurements now in hand for 60 active galaxies, and more on the horizon, we are able to directly constrain black hole masses, derive scaling relationships that allow large numbers of black hole mass estimates throughout the observable Universe, and begin investigating the detailed geometry and kinematics of the broad line region. Reverberation mapping is therefore one of the few techniques available that will allow a deeper understanding of the physical mechanisms involved in AGN feeding and feedback at very small scales, as well as constraints on the growth and evolution of black holes across cosmic time. In this contribution, I will briefly review the background, implementation, and major results derived from this high angular resolution technique.
We study the baryon content of low mass galaxies selected from the Sloan Digital Sky Survey (SDSS DR8), focusing on galaxies in isolated environments where the complicating physics of galaxy-galaxy interactions are minimized. We measure neutral hydrogen (HI) gas masses and line-widths for 148 isolated galaxies with stellar mass between $10^7$ and $10^{9.5} M_{\odot}$. We compare isolated low mass galaxies to more massive galaxies and galaxies in denser environments by remeasuring HI emission lines from the Arecibo Legacy Fast ALFA (ALFALFA) survey 40% data release. All isolated low mass galaxies either have large atomic gas fractions or large atomic gas fractions cannot be ruled out via their upper limits. We measure a median atomic gas fraction of $f_{\rm gas} = 0.82 \pm 0.13$ for our isolated low mass sample with no systems below 0.30. At all stellar masses, the correlations between galaxy radius, baryonic mass and velocity width are not significantly affected by environment. Finally, we estimate a median baryon to total dynamical mass fraction of $f_{\rm baryon, disk} = 0.15 \pm 0.18$. We also estimate two different median baryon to halo mass fractions using the results of semi-analytic models $(f_{\rm baryon, halo} = 0.04 \pm 0.06)$ and abundance matching $(f_{\rm baryon, halo} = 0.04 \pm 0.02)$. Baryon fractions estimated directly using HI observations appear independent of environment and maximum circular velocity, while baryon fractions estimated using abundance matching show a significant depletion of baryons at low maximum circular velocities.
I briefly review advances in the understanding and modeling of relativistic stellar dynamics around massive black holes (MBHs) in galactic nuclei, following the inclusion of coherent relaxation and of secular processes in a new formal analytic description of the dynamics.
Using self-consistent three-dimensional (3D) N-body simulations, we investigate the physical properties of non-axisymmetric features in a disk galaxy created by a tidal interaction with its companion. The primary galaxy consists of a stellar disk, a bugle, and a live halo, corresponding to Milky-Way type galaxies, while the companion is represented by a halo alone. We vary the companion mass and the pericenter distance to explore situations with differing tidal strength parameterized by either the relative tidal force P or the relative imparted momentum S. We find that the formation of a tidal tail in the outer parts requires P > 0.05 or S > 0.07. A stronger interaction results in a stronger, less wound tail that forms earlier. Similarly, a stronger tidal forcing produces stronger, more loosely wound spiral arms in the inner parts. The arms are approximately logarithmic in shape, with both amplitude and pitch angle decaying with time. The derived pattern speed decreases with radius and is close to the Omega-kappa/2 curve at late time, with Omega and kappa denoting the angular and epicycle frequencies, respectively. This suggests that the tidally-induced spiral arms are most likely kinematic density waves weakly modified by self-gravity. Compared to the razor-thin counterparts, arms in the 3D models are weaker, have a smaller pitch angle, and wind and decay more rapidly. The 3D density structure of the arms is well described by the concentrated and sinusoidal models when the arms are in the nonlinear and linear regimes, respectively. We demonstrate that dynamical friction between interacting galaxies transfers the orbital angular momentum of one galaxy to the spin angular momentum of the companion halo.
We present a sample of about 120,000 red clump candidates selected from the LAMOST DR2 catalog based on the empirical distribution model in the effective temperature vs. surface gravity plane. Although, in general, red clump stars are considered as the standard candle, they do not exactly stay in a narrow range of absolute magnitude, but may extend to more than 1 magnitude depending on their initial mass. Consequently, conventional oversimplified distance estimations with assumption of fixed luminosity may lead to systematic bias related to the initial mass or the age, which may potentially affect the study of the evolution of the Galaxy with red clump stars. We therefore employ an isochrone-based method to estimate the absolute magnitude of red clump stars from their observed surface gravities, effective temperatures, and metallicities. We verify that the estimation well removes the systematics and provide an initial mass/age independent distance estimates with accuracy less than 10%.
We report the discovery of a new ultra-faint Milky Way satellite candidate, Horologium II, detected in the Dark Energy Survey Y1A1 public data. Horologium II features a half light radius of $r_{h}=47\pm10$ pc and a total luminosity of $M_{V}=-2.6^{+0.2}_{-0.3}$ that place it in the realm of ultra-faint dwarf galaxies on the size-luminosity plane. The stellar population of the new satellite is consistent with an old ($\sim13.5$ Gyr) and metal-poor ([Fe/H]$\sim-2.1$) isochrone at a distance modulus of $(m-M)=19.46$, or a heliocentric distance of 78 kpc, in the color-magnitude diagram. Horologium II has a distance similar to the Sculptor dwarf spheroidal galaxy (79 kpc) and the recently reported ultra-faint satellites Eridanus III (87 kpc) and Horologium I (79 kpc). All four satellites are well aligned on the sky, which suggests a possible common origin. As Sculptor is moving on a retrograde orbit within the Vast Polar Structure when compared to the other classical MW satellite galaxies including the Magellanic Clouds, this hypothesis can be tested once proper motion measurements become available.
Local starbursts have a higher efficiency of converting gas into stars, as compared to typical star-forming galaxies at a given stellar mass, possibly indicative of different modes of star formation. With the peak epoch of galaxy formation occurring at z > 1, it remains to be established whether such an efficient mode of star formation is occurring at high-redshift. To address this issue, we measure the CO molecular gas content of seven high-redshift starburst galaxies with ALMA and IRAM/PdBI. Our sample is selected from the FMOS-COSMOS near-infrared spectroscopic survey of star-forming galaxies at z ~ 1.6 with Subaru. All galaxies have star formation rates (~300-800 Msolar/yr) elevated, by at least four times, above the star-forming main sequence. We detect CO emission in all cases at high significance, indicative of plentiful gas supplies (f_gas ~ 30-50%). Even more compelling, we firmly establish for the first time that starbursts at high redshift systematically have a lower ratio of CO to total infrared luminosity as compared to typical 'main-sequence' star-forming galaxies, although with an offset less than expected based on past studies of local starbursts. We put forward a hypothesis that there exists a continuous increase in star formation efficiency with elevation from the main sequence with galaxy mergers as a possible physical driver. In support of this scenario, our high-redshift sample is similar in other respects to local starbursts such as being metal rich and having a higher ionization state of the ISM.
We present images from the Solar Blind Channel on HST that resolve hundreds of far ultraviolet (FUV) emitting stars in two ~1 kpc$^2$ interarm regions of the grand-design spiral M101. The luminosity functions of these stars are compared with predicted distributions from simple star formation histories, and are best reproduced when the star formation rate has declined recently (past 10-50 Myr). This pattern is consistent with stars forming within spiral arms and then streaming into the interarm regions. We measure the diffuse FUV surface brightness after subtracting all of the detected stars, clusters and background galaxies. A residual flux is found for both regions which can be explained by a mix of stars below our detection limit and scattered FUV light. The amount of scattered light required is much larger for the region immediately adjacent to a spiral arm, a bright source of FUV photons.
We present a study of the recent star formation of 30 Doradus in the Large Magellanic Cloud (LMC) using the panchromatic imaging survey Hubble Tarantula Treasury Project (HTTP). In this paper we focus on the stars within 20 pc of the center of the massive ionizing cluster of 30 Doradus, NGC 2070. We recovered the star formation history by comparing deep optical and NIR color-magnitude diagrams (CMDs) with state-of-the-art synthetic CMDs generated with the latest PARSEC models, which include all stellar phases from pre-main sequence to post- main sequence. For the first time in this region we are able to measure the star formation using intermediate and low mass stars simultaneously. Our results suggest that NGC2070 experienced a prolonged activity. In particular, we find that the star formation in the region: i) exceeded the average LMC rate ~ 20 Myr ago; ii) accelerated dramatically ~ 7 Myr ago; and iii) reached a peak value 1-3 Myr ago. We did not find significant deviations from a Kroupa initial mass function down to 0.5 Msun. The average internal reddening E(B-V) is found to be between 0.3 and 0.4 mag.
A search for high-energy neutrinos coming from the direction of the Galactic Centre is performed using the data recorded by the ANTARES neutrino telescope from 2007 to 2012. The event selection criteria are chosen to maximise the sensitivity to possible signals produced by the self-annihilation of weakly interacting massive particles accumulated around the centre of the Milky Way with respect to the atmospheric background. After data unblinding, the number of neutrinos observed in the line of sight of the Galactic Centre is found to be compatible with background expectations. The 90% C.L. upper limits in terms of the neutrino+anti-neutrino flux, $\rm \Phi_{\nu_{\mu}+\bar{\nu}_\mu}$, and the velocity averaged annihilation cross-section, $\rm <\sigma_{A}v>$, are derived for the WIMP self-annihilation channels into $\rm b\bar{b},W^{+}W^{-},\tau^{+}\tau^{-},\mu^{+}\mu^{-},\nu\bar{\nu}$. The ANTARES limits for $\rm <\sigma_{A}v>$ are shown to be the most stringent for a neutrino telescope over the WIMP masses $\rm 25\,GeV < M_{WIMP} < 10\,TeV$.
We study the H2 molecular content in high redshift damped Lyman-alpha systems (DLAs) as a function of the HI column density. We find a significant increase of the H2 molecular content around log N(HI) (cm^-2)~21.5-22, a regime unprobed until now in intervening DLAs, beyond which the majority of systems have log N(H2) > 17. This is in contrast with lines of sight towards nearby stars, where such H2 column densities are always detected as soon as log N(HI)>20.7. This can qualitatively be explained by the lower average metallicity and possibly higher surrounding UV radiation in DLAs. However, unlike in the Milky Way, the overall molecular fractions remain modest, showing that even at a large N(HI) only a small fraction of overall HI is actually associated with the self-shielded H2 gas. Damped Lyman-alpha systems with very high-N(HI) probably arise along quasar lines of sight passing closer to the centre of the host galaxy where the gas pressure is higher. We show that the colour changes induced on the background quasar by continuum (dust) and line absorption (HI Lyman and H2 Lyman & Werner bands) in DLAs with log N(HI)~22 and metallicity ~1/10 solar is significant, but not responsible for the long-discussed lack of such systems in optically selected samples. Instead, these systems are likely to be found towards intrinsically fainter quasars that dominate the quasar luminosity function. Colour biasing should in turn be severe at higher metallicities.
Links to: arXiv, form interface, find, astro-ph, recent, 1505, contact, help (Access key information)
We exploit long-baseline ALMA sub-mm observations of the lensed star-forming galaxy SDP 81 at z=3.042 to investigate the properties of inter-stellar medium on scales of 50-100pc. The kinematics of the CO gas within this system are well described by a rotationally-supported disk with an inclination-corrected rotation speed, v=320+/-20km/s and a dynamical mass of M=(3.5+/-1.0)x10^10Mo within a radius of 1.5 kpc. The disk is gas rich and unstable, with a Toomre parameter, Q=0.30+/-0.10 and so should collapse in to star-forming regions with Jeans length L_J~130pc. We identify five star-forming regions within the ISM on these scales and show that their scaling relations between luminosity, line-widths and sizes are significantly offset from those typical of molecular clouds in local Galaxies (Larson's relations). These offsets are likely to be caused by the high external hydrostatic pressure for the interstellar medium (ISM), P/kB=(40+/-20)x10^7K/cm3, which is ~10,000x higher than the typical ISM pressure in the Milky Way. The physical conditions of the star-forming ISM and giant molecular clouds appears to be similar to the those found in the densest environments in the local Universe, such as those in the Galactic center.
We report Hubble Space Telescope/Cosmic Origins Spectrograph observations of the Ly$\alpha$ emission and interstellar absorption lines in a sample of ten star-forming galaxies at z~0.2. Selected on the basis of high equivalent width optical emission lines, the sample, dubbed "Green Peas," make some of the best analogs for young galaxies in an early Universe. We detect Ly$\alpha$ emission in all ten galaxies, and 9/10 show double-peaked line profiles suggestive of low H I column density. We measure Ly$\alpha$/H$\alpha$ flux ratios of 0.5-5.6, implying that 5% to 60% of Ly$\alpha$ photons escape the galaxies. These data confirm previous findings that low-ionization metal absorption (LIS) lines are weaker when Ly$\alpha$ escape fraction and equivalent width are higher. However, contrary to previously favored interpretations of this trend, increased Ly$\alpha$ output cannot be the result of a varying H I covering: the Lyman absorption lines (Ly$\beta$ and higher) show a covering fraction near unity for gas with N_{H I} >~ 10^{16} cm^{-2}. Moreover, we detect no correlation between Ly$\alpha$ escape and the outflow velocity of the LIS lines, suggesting that kinematic effects do not explain the range of Ly$\alpha$/H$\alpha$ flux ratios in these galaxies. In contrast, we detect a strong anti-correlation between the Ly$\alpha$ escape fraction and the velocity separation of the Ly$\alpha$ emission peaks, driven primarily by the velocity of the blue peak. As this velocity separation is sensitive to H I column density, we conclude that Ly$\alpha$ escape in these Green Peas is likely regulated by the H I column density rather than outflow velocity or H I covering fraction.
We present high-resolution 870-um ALMA continuum maps of 30 bright sub-millimeter sources in the UKIDSS UDS field. These sources are selected from deep, 1-square degrees 850-um maps from the SCUBA--2 Cosmology Legacy Survey, and are representative of the brightest sources in the field (median SCUBA2 flux S_850=8.7+/-0.4 mJy). We detect 52 sub-millimeter galaxies (SMGs) at >4-sigma significance in our 30 ALMA maps. In 61+/-17% of the ALMA maps the single-dish source comprises a blend of >=2 SMGs, where the secondary SMGs are Ultra--Luminous Infrared Galaxies (ULIRGs) with L_IR>10^12 Lo. The brightest SMG contributes on average 80+/-4% of the single-dish flux density, and in the ALMA maps containing >=2 SMGs the secondary SMG contributes 25+/-3% of the integrated ALMA flux. We construct source counts and show that multiplicity boosts the apparent single-dish cumulative counts by 20% at S_870>7.5mJy, and by 60% at S_870>12mJy. We combine our sample with previous ALMA studies of fainter SMGs and show that the counts are well-described by a double power-law with a break at 8.5+/-0.6mJy. The break corresponds to a luminosity of ~6x10^12Lsol or a star-formation rate of ~1000Mo/yr. For the typical sizes of these SMGs, which are resolved in our ALMA data with r=1.2+/-0.1kpc, this yields a limiting SFR density of ~100Msol/yr/kpc2. Finally, the number density of S_870>2mJy SMGs is 80+/-30 times higher than that derived from blank-field counts. An over-abundance of faint SMGs is inconsistent with line-of-sight projections dominating multiplicity in the brightest SMGs, and indicates that a significant proportion of these high-redshift ULIRGs must be physically associated.
We study the environments of 49 WISE/NVSS-selected dusty, hyper-luminous, z~2 quasars using the Atacama Large Millimeter/Sub-millimeter Array (ALMA) 345GHz images. We find that 17 of the 49 WISE/NVSS sources show additional sub-mm galaxies within the ALMA primary beam, probing scales within ~150 kpc. We find a total of 23 additional sub-mm sources, four of which in the field of a single WISE/NVSS source. The measured 870 um source counts are ~10 times expectations for unbiased regions, suggesting such hyper-luminous dusty quasars are excellent at probing high-density peaks.
We present radio and infrared observations indicating on-going star formation activity inside the $\sim2-5$ pc circumnuclear ring at the Galactic center. Collectively these measurements suggest a continued disk-based mode of on-going star formation has taken place near Sgr A* over the last few million years. First, VLA observations with spatial resolution 2.17$"\times0.81"$ reveal 13 water masers, several of which have multiple velocity components. The presence of interstellar water masers suggests gas densities that are sufficient for self-gravity to overcome the tidal shear of the 4$\times10^6$ \msol\, black hole. Second, SED modeling of stellar sources indicate massive YSO candidates interior to the molecular ring, supporting in-situ star formation near Sgr A* and appear to show a distribution similar to that of the counter-rotating disks of $\sim$100 OB stars orbiting Sgr A*. Some YSO candidates (e.g., IRS~5) have bow shock structures suggesting that they have have gaseous disks that are phototoevaporated and photoionized by the strong radiation field. Third, we detect clumps of SiO (2-1) and (5-4) line emission in the ring based on CARMA and SMA observations. The FWHM and luminosity of the SiO emission is consistent with shocked protostellar outflows. Fourth, two linear ionized features with an extent of $\sim0.8$ pc show blue and redshifted velocities between $+50$ and $-40$ \kms, suggesting protostellar jet driven outflows with mass loss rates of $\sim5\times10^{-5}$ solar mass yr$^{-1}$. Finally, we present the imprint of radio dark clouds at 44 GHz, representing a reservoir of molecular gas that feeds star formation activity close to Sgr A*.
Modified Newtonian dynamics (MoND) is an empirical theory originally proposed to explain the rotation curves of spiral galaxies by modifying the gravitational acceleration, rather than by invoking dark matter. Here, we set constraints on MoND using an up-to-date compilation of kinematic tracers of the Milky Way and a comprehensive collection of morphologies of the baryonic component in the Galaxy. In particular, we find that the so-called "standard" interpolating function cannot explain at the same time the rotation curve of the Milky Way and that of external galaxies for any of the baryonic models studied, while the so-called "simple" interpolating function remains viable for a subset of models. Upcoming astronomical observations will refine our knowledge on the morphology of baryons and will ultimately confirm or rule out the validity of MoND in the Milky Way. We also present constraints on MoND-like theories without making any assumptions on the interpolating function.
We present observations and analysis of an unusual [C II] emission line in the very luminous QSO SDSS J155426.16+193703.0 at z~4.6. The line is extremely broad (FWHM 735 km/s) and seems to have a flat-topped or double-peaked line profile. A velocity map of the line shows a gradient across the source that indicates large-scale rotation of star-forming gas. Together, the velocity map and line profile suggest the presence of a massive rotating disc with a dynamical mass M_dyn > 5x10^10 M_sun. Using the assumption of a rotating disc origin, we employ an empirical relation between galaxy disc circular velocity and bulge velocity dispersion (sigma) to estimate that sigma > 310 km/s, subject to a correction for the unknown disc inclination. This result implies that this source is consistent with the local M--sigma relation, or offset at most by an order of magnitude in black hole mass. In contrast, the assumption of a bulge origin for the [C II] emission line would lead to a conclusion that the black hole is nearly two orders of magnitude more massive than predicted by the M--sigma relation, similar to previous findings for other high-redshift QSOs. As disc rotation may be a common origin for [C II] emission at high redshifts, these results stress that careful consideration of dynamical origins is required when using observations of this line to derive properties of high-redshift galaxies.
We investigate the global cold dust properties of 85 nearby (z < 0.5) QSOs, chosen from the Palomar-Green sample of optically luminous quasars. We determine their infrared spectral energy distributions and estimate their rest-frame luminosities by combining Herschel data from 70 to 500 microns with near-infrared and mid-infrared measurements from the Two Micron All Sky Survey (2MASS) and the Wide-Field Infrared Survey Explorer (WISE). In most sources the far-infrared (FIR) emission can be attributed to thermally heated dust. Single temperature modified black body fits to the FIR photometry give an average dust temperature for the sample of 33~K, with a standard deviation of 8~K, and an average dust mass of 7E6 Solar Masses with a standard deviation of 9E6 Solar Masses. Estimates of star-formation that are based on the FIR continuum emission correlate with those based on the 11.3 microns PAH feature, however, the star-formation rates estimated from the FIR continuum are higher than those estimated from the 11.3 microns PAH emission. We attribute this result to a variety of factors including the possible destruction of the PAHs and that, in some sources, a fraction of the FIR originates from dust heated by the active galactic nucleus and by old stars.
We present results from Washington CT1 photometry for eleven star fields located in the western outskirts of the Small Magellanic Cloud (SMC), which cover angular distances to its centre from 2 up to 13 degrees (~ 2.2 - 13.8 kpc). The colour- magnitude diagrams, cleaned from the unavoidable Milky Way (MW) and background galaxy signatures, reveal that the most distant dominant main sequence (MS) stellar populations from the SMC centre are located at an angular distance of ~ 5.7 deg (6.1 kpc); no sign of farther clear SMC MS is visible other than the residuals from the MW/background field contamination. The derived ages and metallicities for the dominant stellar populations of the western SMC periphery show a constant metallicity level ([Fe/H] = -1.0 dex) and an approximately constant age value (~ 7-8 Gyr). Their age-metallicity relationship (AMR) do not clearly differ from the most comprehensive AMRs derived for almost the entire SMC main body. Finally, the range of ages of the dominant stellar populations in the western SMC periphery confirms that the major stellar mass formation activity at the very early galaxy epoch peaked ~ 7-8 Gyr ago.
Motivated by the discovery of prolate rotation of stars in Andromeda II, a dwarf spheroidal companion of M31, we study the origin of this type of streaming motion via mergers of disky dwarf galaxies. We simulate merger events between two identical dwarfs changing the initial inclination of their disks with respect to the orbit and the amount of orbital angular momentum. On radial orbits the amount of prolate rotation in the merger remnants correlates strongly with the inclination of the disks and is well understood as due to the conservation of the angular momentum component of the disks along the merger axis. For non-radial orbits prolate rotation may still be produced if the orbital angular momentum is initially not much larger than the intrinsic angular momentum of the disks. The orbital structure of the remnants with significant rotation is dominated by box orbits in the center and long-axis tubes in the outer parts. We also detect significant figure rotation resulting from the tidal distortion of the disks before the merger. The frequency analysis of stellar orbits in the plane perpendicular to the major axis reveals the presence of two families roughly corresponding to inner and outer long-axis tubes. The fraction of inner tubes is largest in the remnant forming from disks oriented most vertically initially and is responsible for the boxy shape of the galaxy. We conclude that prolate rotation may result from mergers with a variety of initial conditions and no fine tuning is necessary to reproduce this feature.
We present a kinematic analysis of the globular cluster systems and diffuse stellar light of four intermediate luminosity (sub-$L^{\ast}$) early-type galaxies in the Virgo cluster based on Gemini/GMOS data. Our galaxy sample is fainter ($-23.8<M_K<-22.7$) than most previous studies, nearly doubling the number of galaxies in this magnitude range that now have GC kinematics. The data for the diffuse light extends to $4R_e$, and the data for the globular clusters reaches 8--$12R_e$. We find that the kinematics in these outer regions are all different despite the fact that these four galaxies have similar photometric properties, and are uniformly classified as "fast rotators" from their stellar kinematics within $1R_e$. The globular cluster systems exhibit a wide range of kinematic morphology. The rotation axis and amplitude can change between the inner and outer regions, including a case of counter-rotation. This difference shows the importance of wide-field kinematic studies, and shows that stellar and GC kinematics can change significantly as one moves beyond the inner regions of galaxies. Moreover, the kinematics of the globular cluster systems can differ from that of the stars, suggesting that the formation of the two populations are also distinct.
The Central Molecular Zone (CMZ) of the Milky Way shows several peculiar properties: a large star formation rate, some of the most massive young star clusters and molecular clouds in the Galaxy, and a twisted ring morphology in molecular gas. In this paper, I use SPH simulations to show that most of these properties can be explained as due to a recent outburst of AGN activity in Sgr A*, the central supermassive black hole of the Milky Way. In particular, the narrow ring of dense gas, massive gas clouds, young star clusters and an elevated SFR can all be caused by the passage of an AGN outflow through the system, which compresses the gas and triggers fragmentation. Furthermore, I show that the asymmetric distribution of gas, as observed in the CMZ, can be produced by outflow-induced instabilities from an initially axisymmetric gas disc. Angular momentum mixing in the disc produces some low angular momentum material, which can subsequently feed Sgr A*. These processes can occur in any galaxy that experiences an AGN episode, leading to bursts of nuclear star formation much stronger than pure bar-driven mass inflows would predict.
The centre of our Galaxy is one of the most studied and yet enigmatic places in the Universe. At a distance of about 8 kpc from our Sun, the Galactic centre (GC) is the ideal environment to study the extreme processes that take place in the vicinity of a supermassive black hole (SMBH). Despite the hostile environment, several tens of early-type stars populate the central parsec of our Galaxy. A fraction of them lie in a thin ring with mild eccentricity and inner radius ~0.04 pc, while the S-stars, i.e. the ~30 stars closest to the SMBH (<0.04 pc), have randomly oriented and highly eccentric orbits. The formation of such early-type stars has been a puzzle for a long time: molecular clouds should be tidally disrupted by the SMBH before they can fragment into stars. We review the main scenarios proposed to explain the formation and the dynamical evolution of the early-type stars in the GC. In particular, we discuss the most popular in situ scenarios (accretion disc fragmentation and molecular cloud disruption) and migration scenarios (star cluster inspiral and Hills mechanism). We focus on the most pressing challenges that must be faced to shed light on the process of star formation in the vicinity of a SMBH.
A binary supermassive black hole loses energy via ejection of stars in a galactic nucleus, until emission of gravitational waves becomes strong enough to induce rapid coalescence. Evolution via the gravitational slingshot requires that stars be continuously supplied to the binary, and it is known that in spherical galaxies the reservoir of such stars is quickly depleted, leading to stalling of the binary at parsec-scale separations. Recent N-body simulations of galaxy mergers and isolated nonspherical galaxies suggest that this stalling may not occur in less idealized systems. However, it remains unclear to what degree these conclusions are affected by collisional relaxation, which is much stronger in the numerical simulations than in real galaxies. In this study, we present a novel Monte Carlo method that can efficiently deal with both collisional and collisionless dynamics, and with galaxy models having arbitrary shapes. We show that without relaxation, the final-parsec problem may be overcome only in triaxial galaxies. Axisymmetry is not enough, but even a moderate departure from axisymmetry is sufficient to keep the binary shrinking. We find that the binary hardening rate is always substantially lower than the maximum possible, "full-loss-cone" rate, and that it decreases with time, but that stellar-dynamical interactions are nevertheless able to drive the binary to coalescence on a timescale <=1 Gyr in any triaxial galaxy.
High-mass X-ray binaries (HMXBs) might have contributed a non-negligible fraction of the energy feedback to the interstellar and intergalactic media at high redshift, becoming important sources for the heating and ionization history of the Universe. However, the importance of this contribution depends on the hypothesized increase in the number of HMXBs formed in low-metallicity galaxies and in their luminosities. In this work we test the aforementioned hypothesis, and quantify the metallicity dependence of HMXB population properties. We compile from the literature a large set of data on the sizes and X-ray luminosities of HMXB populations in nearby galaxies with known metallicities and star formation rates. We use Bayesian inference to fit simple Monte Carlo models that describe the metallicity dependence of the size and luminosity of the HMXB populations. We find that HMXBs are typically ten times more numerous per unit star formation rate in low-metallicity galaxies (12 + log(O/H) < 8, namely < 20% solar) than in solar-metallicity galaxies. The metallicity dependence of the luminosity of HMXBs is small compared to that of the population size. Our results support the hypothesis that HMXBs are more numerous in low-metallicity galaxies, implying the need to investigate the feedback in the form of X-rays and energetic mass outflows of these high-energy sources during cosmic dawn.
We report the multi-wavelength identification of the X-ray sources found in the Subaru-XMM-Newton Deep Survey (SXDS) using deep imaging data covering the wavelength range between the far-UV to the mid-IR. We select a primary counterpart of each X-ray source by applying the likelihood ratio method to R-band, 3.6micron, near-UV, and 24micron source catalogs as well as matching catalogs of AGN candidates selected in 1.4GHz radio and i'-band variability surveys. Once candidates of Galactic stars, ultra-luminous X-ray sources in a nearby galaxy, and clusters of galaxies are removed there are 896 AGN candidates in the sample. We conduct spectroscopic observations of the primary counterparts with multi-object spectrographs in the optical and NIR; 65\% of the X-ray AGN candidates are spectroscopically-identified. For the remaining X-ray AGN candidates, we evaluate their photometric redshift with photometric data in 15 bands. Utilising the multi-wavelength photometric data of the large sample of X-ray selected AGNs, we evaluate the stellar masses, M*, of the host galaxies of the narrow-line AGNs. The distribution of the stellar mass is remarkably constant from z=0.1 to 4.0. The relation between M* and 2--10 keV luminosity can be explained with strong cosmological evolution of the relationship between the black hole mass and M*. We also evaluate the scatter of the UV-MIR spectral energy distribution (SED) of the X-ray AGNs as a function of X-ray luminosity and absorption to the nucleus. The scatter is compared with galaxies which have redshift and stellar mass distribution matched with the X-ray AGN. The UV-NIR SEDs of obscured X-ray AGNs are similar to those of the galaxies in the matched sample. In the NIR-MIR range, the median SEDs of X-ray AGNs are redder, but the scatter of the SEDs of the X-ray AGN broadly overlaps that of the galaxies in the matched sample.
High--velocity features (HVF), especially of Ca II, are frequently seen in Type Ia supernovae observed prior to B-band maximum (Bmax). These HVF start at more than 25,000 km/s in the days after first light, and slow to about 18,000 km/s near Bmax. To recreate the Ca II near-infrared triplet (CaNIR) HVF in SN 2011fe, we consider the interaction between a Type Ia supernova and a compact circumstellar shell, employing a hydrodynamic 1-D simulation using FLASH. We generate synthetic spectra from the hydrodynamic results using syn++. We show that the CaNIR HVF and its velocity evolution is better explained by a supernova model interacting with a shell than a model without a shell, and briefly discuss the implications for progenitor models.
To study the atomic, molecular and ionized emission of Giant Molecular Clouds (GMCs), we have initiated a Large Program with the VLA: 'THOR - The HI, OH, Recombination Line survey of the Milky Way'. We map the 21cm HI line, 4 OH lines, 19 H_alpha recombination lines and the continuum from 1 to 2 GHz of a significant fraction of the Milky Way (l=15-67deg, |b|<1deg) at ~20" resolution. In this paper, we focus on the HI emission from the W43 star-formation complex. Classically, the HI 21cm line is treated as optically thin with column densities calculated under this assumption. This might give reasonable results for regions of low-mass star-formation, however, it is not sufficient to describe GMCs. We analyzed strong continuum sources to measure the optical depth, and thus correct the HI 21cm emission for optical depth effects and weak diffuse continuum emission. Hence, we are able to measure the HI mass of W43 more accurately and our analysis reveals a lower limit of M=6.6x10^6 M_sun, which is a factor of 2.4 larger than the mass estimated with the assumption of optically thin emission. The HI column densities are as high as N(HI)~150 M_sun/pc^2 ~ 1.9x10^22 cm^-2, which is an order of magnitude higher than for low mass star formation regions. This result challenges theoretical models that predict a threshold for the HI column density of ~10 M_sun/pc^2, at which the formation of molecular hydrogen should set in. By assuming an elliptical layered structure for W43, we estimate the particle density profiles. While at the cloud edge atomic and molecular hydrogen are well mixed, the center of the cloud is strongly dominated by molecular hydrogen. We do not identify a sharp transition between hydrogen in atomic and molecular form. Our results are an important characterization of the atomic to molecular hydrogen transition in an extreme environment and challenges current theoretical models.
In the present paper the results obtained in the investigation of possible diurnal effects for low-energy single-hit scintillation events of DAMA/LIBRA-phase1 (1.04 ton $\times$ yr exposure) have been analysed in terms of an effect expected in case of Dark Matter (DM) candidates inducing nuclear recoils and having high cross-section with ordinary matter, which implies low DM local density in order to fulfill the DAMA/LIBRA DM annual modulation results. This effect is due to the different Earth depths crossed by those DM candidates during the sidereal day.
We use the "Evolution and Assembly of GaLaxies and their Environments" ( EAGLE ) suite of hydrodynamical cosmological simulations to measure offsets between the centres of stellar and dark matter components of galaxies. We find that the vast majority (>95%) of the simulated galaxies display an offset smaller than the gravitational softening length of the simulations ($\epsilon = 700$ pc), both for field galaxies and satellites in clusters and groups. We also find no systematic trailing or leading of the dark matter along a galaxy's direction of motion. The offsets are consistent with being randomly drawn from a Maxwellian distribution with $\sigma = 196$ pc. Since astrophysical effects produce no feasible analogues for the $1.62^{+0.47}_{-0.49}$ kpc offset recently observed in Abell 3827, this observational result is in tension with the collisionless cold dark matter model assumed in the simulations.
Accurate photometric redshift are a lynchpin for many future experiments to pin down the cosmological model and for studies of galaxy evolution. In this study, a novel sparse regression framework for photometric redshift estimation is presented. Data from a simulated survey was used to train and test the proposed models. We show that approaches which include careful data preparation and model design offer a significant improvement in comparison with several competing machine learning algorithms. Standard implementation of most regression algorithms has as the objective the minimization of the sum of squared errors. For redshift inference, however, this induces a bias in the posterior mean of the output distribution, which can be problematic. In this paper we optimize to directly target minimizing $\Delta z = (z_\textrm{s} - z_\textrm{p})/(1+z_\textrm{s})$ and address the bias problem via a distribution-based weighting scheme, incorporated as part of the optimization objective. The results are compared with other machine learning algorithms in the field such as Artificial Neural Networks (ANN), Gaussian Processes (GPs) and sparse GPs. The proposed framework reaches a mean absolute $\Delta z = 0.002(1+z_\textrm{s})$, with a maximum absolute error of 0.0432, over the redshift range of $0.2 \le z_\textrm{s} \le 2$, a factor of three improvement over standard ANNs used in the literature. We also investigate how the relative size of the training affects the photometric redshift accuracy. We find that a training set of $>$30 per cent of total sample size, provides little additional constraint on the photometric redshifts, and note that our GP formalism strongly outperforms ANN in the sparse data regime.
Links to: arXiv, form interface, find, astro-ph, recent, 1505, contact, help (Access key information)
The Cheshire Cat is a relatively poor group of galaxies dominated by two luminous elliptical galaxies surrounded by at least four arcs from gravitationally lensed background galaxies that give the system a humorous appearance. Our combined optical/X-ray study of this system reveals that it is experiencing a line of sight merger between two groups with a roughly equal mass ratio with a relative velocity of ~1350 km/s. One group was most likely a low-mass fossil group, while the other group would have almost fit the classical definition of a fossil group. The collision manifests itself in a bimodal galaxy velocity distribution, an elevated central X-ray temperature and luminosity indicative of a shock, and gravitational arc centers that do not coincide with either large elliptical galaxy. One of the luminous elliptical galaxies has a double nucleus embedded off-center in the stellar halo. The luminous ellipticals should merge in less than a Gyr, after which observers will see a massive 1.2-1.5 x 10^14 solar mass fossil group with an M_r = -24.0 brightest group galaxy at its center. Thus, the Cheshire Cat offers us the first opportunity to study a fossil group progenitor. We discuss the limitations of the classical definition of a fossil group in terms of magnitude gaps between the member galaxies. We also suggest that if the merging of fossil (or near-fossil) groups is a common avenue for creating present-day fossil groups, the time lag between the final galactic merging of the system and the onset of cooling in the shock-heated core could account for the observed lack of well-developed cool cores in some fossil groups.
Using a cosmological N-body numerical simulation of the formation of a galaxy cluster- sized halo, we analyze the temporal evolution of its globular cluster population. We follow the dynamical evolution of 38 galactic dark matter halos orbiting in a galaxy cluster that at redshift z=0 has a virial mass of 1.71 * 10 ^14 Msol h^-1. In order to mimic both "blue" and "red" populations of globular clusters, for each galactic halo we select two different sets of particles at high redshift (z ~ 1), constrained by the condition that, at redshift z=0, their average radial density profiles are similar to the observed profiles. As expected, the general galaxy cluster tidal field removes a significant fraction of the globular cluster populations to feed the intracluster population. On average, halos lost approximately 16% and 29% of their initial red and blue globular cluster populations, respectively. Our results suggest that these fractions strongly depend on the orbital trajectory of the galactic halo, specifically on the number of orbits and on the minimum pericentric distance to the galaxy cluster center that the halo has had. At a given time, these fractions also depend on the current clustercentric distance, just as observations show that the specific frequencyof globular clusters S_N depends on their clustercentric distance.
We revisit the relation between magnetic-field strength ($B$) and gas density ($\rho$) for contracting interstellar clouds and fragments (or, cores), which is central in observationally determining the dynamical importance of magnetic fields in cloud evolution and star formation. Recently, it has been claimed that a relation $B \propto \rho^{2/3} $ is statistically preferred over $B \propto \rho^{1/2}$ in molecular clouds, when magnetic field detections and nondetections from Zeeman observations are combined. This finding has unique observational implications on cloud and core geometry: The relation $B \propto \rho^{2/3} $ can only be realized under spherical contraction. However, no indication of spherical geometry can be found for the objects used in the original statistical analysis of the $B-\rho$ relation. We trace the origin of the inconsistency to simplifying assumptions in the statistical model used to arrive at the $B\propto \rho^{2/3}$ conclusion and to an underestimate of observational uncertainties in the determination of cloud and core densities. We show that, when these restrictive assumptions are relaxed, $B \propto \rho^{1/2}$ is the preferred relation for the (self-gravitating) molecular-cloud data, as theoretically predicted four decades ago.
We explore trends in galaxy properties with Mpc-scale structures using catalogues of environment and large scale structure from the Galaxy And Mass Assembly (GAMA) survey. Existing GAMA catalogues of large scale structure, group and pair membership allow us to construct galaxy stellar mass functions for different environmental types. To avoid simply extracting the known underlying correlations between galaxy properties and stellar mass, we create a mass matched sample of galaxies with stellar masses between $9.5 \leq \log{M_*/h^{-2} M_{\odot}} \leq 11$ for each environmental population. Using these samples, we show that mass normalised galaxies in different large scale environments have similar energy outputs, $u-r$ colours, luminosities, and morphologies. Extending our analysis to group and pair environments, we show galaxies that are not in groups or pairs exhibit similar characteristics to each other regardless of broader environment. For our mass controlled sample, we fail to see a strong dependence of S\'{e}rsic index or galaxy luminosity on halo mass, but do find that it correlates very strongly with colour. Repeating our analysis for galaxies that have not been mass controlled introduces and amplifies trends in the properties of galaxies in pairs, groups, and large scale structure, indicating that stellar mass is the most important predictor of the galaxy properties we examine, as opposed to environmental classifications.
We present a method to estimate distances to stars with spectroscopically derived stellar parameters. The technique is a Bayesian approach with likelihood estimated via comparison of measured parameters to a grid of stellar isochrones, and returns a posterior probability density function for each star's absolute magnitude. This technique is tailored specifically to data from the Large Sky Area Multi-object Fiber Spectroscopic Telescope (LAMOST) survey. Because LAMOST obtains roughly 3000 stellar spectra simultaneously within each ~5-degree diameter "plate" that is observed, we can use the stellar parameters of the observed stars to account for the stellar luminosity function and target selection effects. This removes biasing assumptions about the underlying populations, both due to predictions of the luminosity function from stellar evolution modeling, and from Galactic models of stellar populations along each line of sight. Using calibration data of stars with known distances and stellar parameters, we show that our method recovers distances for most stars within ~20%, but with some systematic overestimation of distances to halo giants. We apply our code to the LAMOST database, and show that the current precision of LAMOST stellar parameters permits measurements of distances with ~40% error bars. This precision should improve as the LAMOST data pipelines continue to be refined.
This article explores the agreement between the predictions of Modified Newtonian Dynamics (MOND) and the rotation curves and stellar velocity dispersion profiles measured by the DiskMass Survey. A bulge-disk decomposition was made for each of the thirty published galaxies, and a MOND Poisson solver was used to simultaneously compute, from the baryonic mass distributions, model rotation curves and vertical velocity dispersion profiles, which were compared to the measured values. The two main free parameters, the stellar disk's mass-to-light ratio ($M/L$) and its exponential scale-height ($h_z$), were estimated by Markov Chain Monte Carlo modelling. The average best-fit K-band stellar mass-to-light ratio was $M/L \simeq 0.55 \pm 0.15$. However, to match the DiskMass Survey data, the vertical scale-heights would have to be in the range $h_z=200$ to $400$ pc which is a factor of two lower than those derived from observations of edge-on galaxies with a similar scale-length. The reason is that modified gravity versions of MOND characteristically require a larger $M/L$ to fit the rotation curve in the absence of dark matter and therefore predict a stronger vertical gravitational field than Newtonian models. It was found that changing the MOND acceleration parameter, the shape of the velocity dispersion ellipsoid, the adopted vertical distribution of stars, as well as the galaxy inclination, within any realistic range, all had little impact on these results.
We obtained [O III] narrow-band imaging and multi-slit MXU spectroscopy of the blue compact dwarf (BCD) galaxy NGC 1705 with FORS2@VLT to derive chemical abundances of PNe and H II regions and, more in general, to characterize the properties of the ionized gas. The auroral [O III]\lambda4363 line was detected in all but one of the eleven analyzed regions, allowing for a direct estimate of their electron temperature. The only object for which the [O III]\lambda4363 line was not detected is a possible low-ionization PN, the only one detected in our data. For all the other regions, we derived the abundances of Nitrogen, Oxygen, Neon, Sulfur and Argon out to ~1 kpc from the galaxy center. We detect for the first time in NGC 1705 a negative radial gradient in the oxygen metallicity of -0.24 \pm 0.08 dex kpc^{-1}. The element abundances are all consistent with values reported in the literature for other samples of dwarf irregular and blue compact dwarf galaxies. However, the average (central) oxygen abundance, 12 + log(O/H)=7.98 \pm 0.05, is ~0.2 dex lower than previous literature estimates for NGC 1705 based on the [O III]\lambda4363 line. From classical emission-line diagnostic diagrams, we exclude a major contribution from shock excitation. On the other hand, the radial behavior of the emission line ratios is consistent with the progressive dilution of radiation with increasing distance from the center of NGC~1705. This suggests that the strongest starburst located within the central ~150 pc is responsible for the ionization of the gas out to at least ~1 kpc. The gradual dilution of the radiation with increasing distance from the center reflects the gradual and continuous transition from the highly ionized H II regions in the proximity of the major starburst into the diffuse ionized gas.
We present extensive calculations of radiative transition rates and electron impact collision strengths for Fe II. The data sets involve 52 levels from the $3d\,^7$, $3d\,^64s$, and $3d\,^54s^2$ configurations. Computations of $A$-values are carried out with a combination of state-of-the-art multiconfiguration approaches, namely the relativistic Hartree--Fock, Thomas--Fermi--Dirac potential, and Dirac--Fock methods; while the $R$-matrix plus intermediate coupling frame transformation, Breit--Pauli $R$-matrix and Dirac $R$-matrix packages are used to obtain collision strengths. We examine the advantages and shortcomings of each of these methods, and estimate rate uncertainties from the resulting data dispersion. We proceed to construct excitation balance spectral models, and compare the predictions from each data set with observed spectra from various astronomical objects. We are thus able to establish benchmarks in the spectral modeling of [Fe II] emission in the IR and optical regions as well as in the UV Fe II absorption spectra. Finally, we provide diagnostic line ratios and line emissivities for emission spectroscopy as well as column densities for absorption spectroscopy. All atomic data and models are available online and through the AtomPy atomic data curation environment.
We apply a cross-correlation technique to infer the $S>3$mJy radio luminosity function (RLF) from the NRAO VLA sky survey (NVSS) to $z\sim3.5$. We measure $\Sigma$ the over density of radio sources around spectroscopically confirmed quasars. $\Sigma$ is related to the space density of radio sources at the distance of the quasars and the clustering strength between the two samples, hence knowledge of one constrains the other. Under simple assumptions we find $\Phi\propto (1+z)^{3.7\pm0.7}$ out to $z\sim2$. Above this redshift the evolution slows and we constrain the evolution exponent to $<1.01$ ($2\sigma$). This behaviour is almost identical to that found by previous authors for the bright end of the RLF potentially indicating that we are looking at the same population. This suggests that the NVSS is dominated by a single population; most likely radio sources associated with high-excitation cold-mode accretion. Inversely, by adopting a previously modelled RLF we can constrain the clustering of high-redshift radio sources and find a clustering strength consistent with $r_0=15.0\pm 2.5$ Mpc up to $z\sim3.5$. This is inconsistent with quasars at low redshift and some measurements of the clustering of bright FRII sources. This behaviour is more consistent with the clustering of lower luminosity radio galaxies in the local universe. Our results indicate that the high-excitation systems dominating our sample are hosted in the most massive galaxies at all redshifts sampled.
We continue the study of collisionless systems governed by additive $r^{-\alpha}$ interparticle forces by focusing on the influence of the force exponent $\alpha$ on radial orbital anisotropy. In this preparatory work we construct the radially anisotropic Osipkov-Merritt phase-space distribution functions for self-consistent spherical Hernquist models with $r^{-\alpha}$ forces and $1\leq\alpha<3$. The resulting systems are isotropic at the center and increasingly dominated by radial orbits at radii larger than the anisotropy radius $r_a$. For radially anisotropic models we determine the minimum value of the anisotropy radius $r_{ac}$ as a function of $\alpha$ for phase-space consistency (such that the phase-space distribution function is nowhere negative for $r_a\geq r_{ac}$). We find that $r_{ac}$ decreases for decreasing $\alpha$, and that the amount of kinetic energy that can be stored in the radial direction relative to that stored in the tangential directions for marginally consistent models increases for decreasing $\alpha$. In particular, we find that isotropic systems are consistent in the explored range of $\alpha$. By means of direct $N$-body simulations we finally verify that the isotropic systems are also stable.
We have analysed a sample of 574 Spitzer 4.5 micron-selected galaxies with [4.5]<23 and Ks>24 (AB) over the UltraVISTA ultra-deep COSMOS field. Our aim is to investigate whether these mid-IR bright, near-IR faint sources contribute significantly to the overall population of massive galaxies at redshifts z>=3. By performing a spectral energy distribution (SED) analysis using up to 30 photometric bands, we have determined that the redshift distribution of our sample peaks at redshifts z~2.5-3.0, and ~32% of the galaxies lie at z>=3. We have studied the contribution of these sources to the galaxy stellar mass function (GSMF) at high redshifts. We found that the [4.5]<23, Ks>24 galaxies produce a negligible change to the GSMF previously determined for Ks<24 sources at 3=<z<4, but their contribution is more important at 4=<z<5, accounting for >~50% of the galaxies with stellar masses Mst>~6 x 10^10 Msun. We also constrained the GSMF at the highest-mass end (Mst>~2 x 10^11 Msun) at z>=5. From their presence at 5=<z<6, and virtual absence at higher redshifts, we can pinpoint quite precisely the moment of appearance of the first most massive galaxies as taking place in the ~0.2 Gyr of elapsed time between z~6 and z~5. Alternatively, if very massive galaxies existed earlier in cosmic time, they should have been significantly dust-obscured to lie beyond the detection limits of current, large-area, deep near-IR surveys.
We study relations between global characteristics of low-redshift (0 < z < 1) compact star-forming galaxies, including absolute optical magnitudes, Hbeta emission-line luminosities (or equivalently star-formation rates), stellar masses, and oxygen abundances. The sample consists of 5182 galaxies with high-excitation HII regions selected from the SDSS DR7 and SDSS/BOSS DR10 surveys adopting a criterion [OIII]4959/Hbeta > 1. These data were combined with the corresponding data for high-redshift (2 < z < 3) star-forming galaxies. We find that in all diagrams low-z and high-z star-forming galaxies are closely related indicating a very weak dependence of metallicity on stellar mass, redshift, and star-formation rate. This finding argues in favour of the universal character of the global relations for compact star-forming galaxies with high-excitation HII regions over redshifts 0 < z < 3.
Many of the transformative processes in the Universe have taken place in
regions obscured by dust, and are best studied with far-IR spectroscopy. We
present the Cryogenic-Aperture Large Infrared-Submillimeter Telescope
Observatory (CALISTO), a 5-meter class, space-borne telescope actively cooled
to 4 K, emphasizing moderate-resolution spectroscopy in the crucial 35 to 600
micron band. CALISTO will enable NASA and the world to study the rise of heavy
elements in the Universe's first billion years, chart star formation and black
hole growth in dust-obscured galaxies through cosmic time, and conduct a census
of forming planetary systems in our region of the Galaxy. CALISTO will
capitalize on rapid progress in both format and sensitivity of far-IR
detectors. Arrays with a total count of a few 100,000 detector pixels will form
the heart of a suite of imaging spectrometers in which each detector reaches
the photon background limit.
This document contains a large overview paper on CALISTO, as well as six 2-3
page scientific white papers, all prepared in response to NASA's Cosmic Origins
Program Analysis Group (COPAG's) request for input on future mission concepts.
The Far-IR Science Interest Group will meet from 3-5 June 2015 with the
intention of reaching consensus on the architecture for the Far-IR Surveyor
mission. This white paper describes one of the architectures to be considered
by the community. One or more companion papers will describe alternative
architectures.
Merging galaxy clusters leave long-lasting signatures on the baryonic and non-baryonic cluster constituents, including shock fronts, cold fronts, X-ray substructure, radio halos, and offsets between the dark matter and the gas components. Using observations from Chandra, the Jansky Very Large Array, the Giant Metrewave Radio Telescope, and the Hubble Space Telescope, we present a multiwavelength analysis of the merging Frontier Fields cluster MACS J0416.1-2403 (z=0.396), which consists of a NE and a SW subclusters whose cores are separated on the sky by ~250 kpc. We find that the NE subcluster has a compact core and hosts an X-ray cavity, yet it is not a cool core. Approximately 450 kpc south-south west of the SW subcluster, we detect a density discontinuity that corresponds to a compression factor of ~1.5. The discontinuity was most likely caused by the interaction of the SW subcluster with a less massive structure detected in the lensing maps SW of the subcluster's center. For both the NE and the SW subclusters, the dark matter and the gas components are well-aligned, suggesting that MACS J0416.1-2403 is a pre-merging system. The cluster also hosts a radio halo, which is unusual for a pre-merging system. The halo has a 1.4 GHz power of (1.06 +/- 0.09) x 10^{24} W Hz^{-1}, which is somewhat lower than expected based on the X-ray luminosity of the cluster. We suggest that we are either witnessing the birth of a radio halo, or have discovered a rare ultra-steep spectrum halo.
An excess in $\gamma$-rays emanating from the galactic centre has recently been observed in the Fermi-LAT data. We investigate the new exciting possibility of fitting the signal spectrum by dark matter annihilating dominantly to a Higgs-pseudoscalar pair. We show that the fit to the $\gamma$-ray excess for the Higgs-pseudoscalar channel can be just as good as for annihilation into bottom-quark pairs. This channel arises naturally in a full model such as the next-to-minimal supersymmetric Standard Model (NMSSM) and we find regions where dark matter relic density, the $\gamma$-ray signal and other experimental constraints, can all be satisfied simultaneously. Annihilation into scalar pairs allows for the possibility of detecting the Higgs or pseudoscalar decay into two photons, providing a smoking-gun signal of the model.
Links to: arXiv, form interface, find, astro-ph, recent, 1505, contact, help (Access key information)