Negative feedback from active galactic nuclei (AGN) is considered a key mechanism in shaping galaxy evolution. Fast, extended outflows are frequently detected in the AGN host galaxies at all redshifts and luminosities, both in ionised and molecular gas. However, these outflows are only "potentially" able to quench star formation and we are still missing a decisive evidence of negative feedback in action. Here we present Spectrograph for INtegral Field Observations in the Near Infrared (SINFONI) H- and K-band integral-field spectroscopic observations of two quasars at $z\sim$2.4 characterised by fast, extended outflows detected through the [OIII]$\lambda$5007 line (Carniani et al. 2015). The high signal-to-noise ratio of our observations allows us to identify faint narrow (FWHM $< 500$ km/s), and spatially extended components in [OIII]$\lambda$5007 and H$\alpha$ emission associated with star formation in the host galaxy. Such star-formation powered emission is spatially anti-correlated with the fast outflows. The ionised outflows therefore appear to be able to suppress star formation in the region where the outflow is expanding. However the detection of narrow, spatially-extended H$\alpha$ emission indicates star formation rates of at least $\sim 50-90$ M$_{\odot}$/yr, suggesting either that AGN feedback does not affect the whole galaxy or that many feedback episodes are required before star formation is completely quenched. On the other hand, the narrow H$\alpha$ emission extending along the edges of the outflow cone may lead also to a positive feedback interpretation. Our results highlight the possible twofold role of galaxy-wide outflows in host galaxy evolution.
We present radio Active Galactic Nuclei (AGN) luminosity functions over the redshift range 0.005 < z < 0.75. The sample from which the luminosity functions are constructed is an optical spectroscopic survey of radio galaxies, identified from matched Faint Images of the Radio Sky at Twenty-cm survey (FIRST) sources and Sloan Digital Sky Survey (SDSS) images.The radio AGN are separated into Low Excitation Radio Galaxies (LERGs) and High Excitation Radio Galaxies (HERGs) using the optical spectra. We derive radio luminosity functions for LERGs and HERGs separately in the three redshift bins (0.005 < z < 0.3, 0.3 < z < 0.5 and 0.5 < z <0.75). The radio luminosity functions can be well described by a double power-law. Assuming this double power-law shape the LERG population displays little or no evolution over this redshift range evolving as ~$(1+z)^{0.06}$ assuming pure density evolution or ~ $(1+z)^{0.46}$ assuming pure luminosity evolution. In contrast, the HERG population evolves more rapidly, best fitted by ~$(1+z)^{2.93}$ assuming a double power-law shape and pure density evolution. If a pure luminosity model is assumed the best fitting HERG evolution is parameterised by ~$(1+z)^{7.41}$. The characteristic break in the radio luminosity function occurs at a significantly higher power (~1 dex) for the HERG population in comparison to the LERGs. This is consistent with the two populations representing fundamentally different accretion modes.
Supernova (SN) blast waves inject energy and momentum into the interstellar medium (ISM), control its turbulent multiphase structure and the launching of galactic outflows. Accurate modelling of the blast wave evolution is therefore essential for ISM and galaxy formation simulations. We present an efficient method to compute the input of momentum, thermal energy, and the velocity distribution of the shock-accelerated gas for ambient media with uniform (and with stellar wind blown bubbles), power-law, and turbulent density distributions. Assuming solar metallicity cooling, the blast wave evolution is followed to the beginning of the momentum conserving snowplough phase. The model recovers previous results for uniform ambient media. The momentum injection in wind-blown bubbles depend on the swept-up mass and the efficiency of cooling, when the blast wave hits the wind shell. For power-law density distributions with $n(r) \sim$ $r^{-2}$ (for $n(r) > n_{_{\rm floor}}$) the amount of momentum injection is solely regulated by the background density $n_{_{\rm floor}}$ and compares to $n_{_{\rm uni}}$ = $n_{_{\rm floor}}$. However, in turbulent ambient media with log-normal density distributions the momentum input can increase by a factor of 2 (compared to the homogeneous case) for high Mach numbers. The average momentum boost can be approximated as $p_{_{\rm turb}}/\mathrm{p_{_{0}}}\ =23.07\, \left(\frac{n_{_{0,\rm turb}}}{1\,{\rm cm}^{-3}}\right)^{-0.12} + 0.82 (\ln(1+b^{2}\mathcal{M}^{2}))^{1.49}\left(\frac{n_{_{0,\rm turb}}}{1\,{\rm cm}^{-3}}\right)^{-1.6}$. The velocity distributions are broad as gas can be accelerated to high velocities in low-density channels. The model values agree with results from recent, computationally expensive, three-dimensional simulations of SN explosions in turbulent media.
Kinetic processes in fractal stellar media are analyzed in terms of the approach developed in our earlier paper (Chumak, Rastorguev, 2016) involving a generalization of the nearest neighbor and random force distributions to fractal media. Diffusion is investigated in the approximation of scale-dependent conditional density based on an analysis of the solutions of the corresponding Langevin equations. It is shown that kinetic parameters (time scales, coefficients of dynamic friction, diffusion, etc.) for fractal stellar media can differ significantly both qualitatively and quantitatively from the corresponding parameters for a quasi-uniform random media with limited fluctuations. The most important difference is that in the fractal case kinetic parameters depend on spatial scale length and fractal dimension of the medium studied. A generalized kinetic equation for stellar media (fundamental equation of stellar dynamics) is derived in the Fokker-Planck approximation with the allowance for the fractal properties of the spatial stellar density distribution. Also derived are its limit forms that can be used to describe small departures of fractal gravitating medium from equilibrium.
Elaborating on a formalism that was first expressed some 40 years ago, we consider the brightness of low-lying mm-wave rotational lines of strongly polar molecules at the threshold of detectability. We derive a simple expression relating the brightness to the line of sight integral of the product of the total gas and molecular number densities and a suitably-defined temperature-dependent excitation rate into the upper level of the transition. Detectability of a line is contingent only on the ability of a molecule to channel enough of the ambient thermal energy into the line and the excitation can be computed in bulk by summing over rates without solving the multi-level rate equations or computing optical depths and excitation temperatures. Results for \hcop, HNC and CS are compared with escape-probability solutions of the rate equations using closed-form expressions for the expected range of validity of our {\it ansatz}, with the result that gas number densities as high as $10^4 \pccc$ or optical depths as high as 100 can be accommodated in some cases. For densities below a well-defined upper bound, the range of validity of the discussion can be cast as an upper bound on the line brightness which is 0.3 K for the J=1-0 lines and 0.8 - 1.7 K for the J=2-1 lines of these species. The discussion casts new light on interpretation of line brightnesses under conditions of weak excitation, simplifies derivation of physical parameters and eliminates the need to construct grids of numerical solutions of the rate equations.
Many Seyfert galaxies display weak 'coronal' emission features corresponding to [Fe VII], [Fe XI] and [Fe XIV] in their optical spectra, whereas elsewhere these lines seem to be entirely absent. These lines appear to highlight zones in the nucleus irradiated by high-energy photons. The presence of these zones and the conditions therein as determined by the relative line strengths and profiles impose important constraints on the physical models of active galactic nuclei, and Seyferts in particular. In 2009 the discovery was announced of the highly unusual spectrum of SDSS J0952+2143, where the coronal lines are exceptionally strong. This paper presents a second object with abnormally strong coronal features, SDSS J1055+5637. The spectrum, line ratios and related parameters are compared to those of SDSS J0952+2143, three AGN with moderate coronal lines and one where the coronal lines are missing altogether. Possible mechanisms are discussed that may account for the stronger than usual coronal features.
We present numerical simulations of molecular clouds (MCs) with self-consistent CO gas-phase and isotope chemistry in various environments. The simulations are post-processed with a line radiative transfer code to obtain $^{12}$CO and $^{13}$CO emission maps for the $J=1\rightarrow0$ rotational transition. The emission maps are analysed with commonly used observational methods, i.e. the $^{13}$CO column density measurement, the virial mass estimate and the so-called $X_{\textrm{CO}}$ (also CO-to-H$_2$) conversion factor, and then the inferred quantities (i.e. mass and column density) are compared to the physical values. We generally find that most methods examined here recover the CO-emitting H$_{2}$ gas mass of MCs within a factor of two uncertainty if the metallicity is not too low. The exception is the $^{13}$CO column density method. It is affected by chemical and optical depth issues, and it measures both the true H$_{2}$ column density distribution and the molecular mass poorly. The virial mass estimate seems to work the best in the considered metallicity and radiation field strength range, even when the overall virial parameter of the cloud is above the equilibrium value. This is explained by a systematically lower virial parameter (i.e. closer to equilibrium) in the CO-emitting regions; in CO emission, clouds might seem (sub-)virial, even when, in fact, they are expanding or being dispersed. A single CO-to-H$_{2}$ conversion factor appears to be a robust choice over relatively wide ranges of cloud conditions, unless the metallicity is low. The methods which try to take the metallicity dependence of the conversion factor into account tend to systematically overestimate the true cloud masses.
We investigate the presence of extended ionized outflows in 18 luminous type 2 AGNs (11 quasars and 7 high luminosity Seyfert 2s) at 0.3<z<0.6 based on VLT-FORS2 spectroscopy. We infer typical lower limits on the radial sizes of the outflows R>=severalx100 pc and upper limits R<1-2 kpc. Our results are inconsistent with related studies which suggest that large scale R~several-15 kpc are ubiquitous in QSO2. We study the possible causes of discrepancy and propose that seeing smearing is the cause of the large inferred sizes. The implications in our understanding of the feedback phenomenon are important since the mass M (through the density), mass injection M(dot) and energy injection rates E(dot) of the outflows become highly uncertain. One conclusion seems unavoidable: M, M(dot), E(dot) are modest or low compared with previous estimations. No evidence is found supporting that typical outflows can affect the interstellar medium of the host galaxies accross spatial scales >~1-2 kpc.
We investigate in detail the escape dynamics in an analytical gravitational model which describes the motion of stars in a quasar galaxy with a disk and a massive nucleus. 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 configuration $(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. We hope our outcomes to be useful for a further understanding of the escape mechanism in active galaxy models.
The escape dynamics in an analytical gravitational model which describes the motion of stars in a binary system of interacting dwarf spheroidal galaxies is investigated in detail. We conduct a 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 configuration $(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 hope our outcomes to be useful for a further understanding of the escape mechanism in binary galaxy models.
The recent detection of the binary black hole merger GW150914 demonstrates the existence of black holes more massive than previously observed in X-ray binaries in our Galaxy. This article explores different scenarios of black hole formation in the context of self-consistent cosmic chemical evolution models that simultaneously match observations of the cosmic star formation rate, optical depth to reionization and metallicity of the interstellar medium. This framework is used to calculate the mass distribution of merging black hole binaries and its evolution with redshift. We also study the implications of the black hole mass distribution for the stochastic gravitational wave background from mergers and from core collapse events.
We present a systematic characterization of multi-wavelength emission from blazar PKS 1510-089 using well-sampled data at infrared(IR)-optical, X-ray and $\gamma$-ray energies. The resulting flux distributions, except at X-rays, show two distinct lognormal profiles corresponding to a high and a low flux level. The dispersions exhibit energy dependent behavior except for the LAT $\gamma$-ray and optical B-band. During the low level flux states, it is higher towards the peak of the spectral energy distribution, with $\gamma$-ray being intrinsically more variable followed by IR and then optical, consistent with mainly being a result of varying bulk Lorentz factor. On the other hand, the dispersions during the high state are similar in all bands expect optical B-band, where thermal emission still dominates. The centers of distributions are a factor of $\sim 4$ apart, consistent with anticipation from studies of extragalactic $\gamma$-ray background with the high state showing a relatively harder mean spectral index compared to the low state.
Links to: arXiv, form interface, find, astro-ph, recent, 1604, contact, help (Access key information)
We present a model for the classification of Coronal-Line Forest Active Galactic Nuclei (CLiF AGN). CLiF AGN are of special interest due to their remarkably large number of emission lines, especially forbidden high ionization lines (FHILs). Rose et al. (2015a) suggest that their emission is dominated by reflection from the inner wall of the obscuring region rather than direct emission from the accretion disk. This makes CLiF AGN laboratories to test AGN-torus models. Modeling AGN as an accreting supermassive black hole, surrounded by a cylinder of dust and gas, we show a relationship between viewing angle and the revealed area of the inner wall. From the revealed area, we can determine the amount of FHIL emission at various angles. We calculate the strength of [FeVII]{\lambda}6087 emission for a number of intermediate angles (30{\deg}, 40{\deg}, and 50{\deg}) and compare the results with the luminosity of the observed emission line from six known CLiF AGN. We find that there is good agreement between our model and the observational results. The model also enables us to determine the relationship between the type 2 : type 1 AGN fraction vs the ratio of torus height to radius, h/r.
Optical spectropolarimetry carried out by Kishimoto et al. (2004) has shown that several luminous type 1 quasars show a strong decrease of the polarized continuum flux in the rest-frame near-UV wavelengths of $\lambda<4000$\AA. In the literature, this spectral feature is interpreted as evidence of the broadened hydrogen Balmer absorption edge imprinted in the accretion disk thermal emission due to the disk atmospheric opacity effect. On the other hand, the quasar flux variability studies have shown that the variable continuum component in UV-optical spectra of quasars, which is considered to be a good indicator of the intrinsic spectral shape of the accretion disk emission, generally have significantly flat spectral shape throughout the near-UV to optical spectral range. To examine whether the disk continuum spectral shapes revealed as the polarized flux and as the variable component spectra are consistent with each other, we carry out multi-band photometric monitoring observations for a sample of four polarization-decreasing quasars of Kishimoto et al. (2004) (4C09.72, 3C323.1, Ton 202, and B2 1208+32) to derive the variable component spectra and compare the spectral shape of them with that of the polarized flux spectra. Contrary to expectation, we confirm that the two spectral components of these quasars have totally different spectral shape in that the variable component spectra are significantly bluer compared to the polarized flux spectra. This discrepancy in the spectral shape may imply either (1) the decrease of polarization degree in the rest-frame UV wavelengths is not indicating the Balmer absorption edge feature but is induced by some unknown (de)polarization mechanisms, or (2) the UV-optical flux variability is occurring preferentially at the hot inner radii of the accretion disk and thus the variable component spectra do not reflect the whole accretion disk emission.
Observation shows that cooling instabilities leading to nebular emission, molecular gas, and star formation in giant galaxies are formed behind buoyantly-rising X-ray bubbles inflated by radio jets launched from massive nuclear black holes. We propose a model where molecular clouds condense from hot but relatively low entropy gas lifted by X-ray bubbles to an altitude where its cooling time is shorter than the time required for it to fall to its equilibrium location in the galaxy i.e., t_c/t_I <~1$. Here the infall time can exceed the free-fall time, t_ff, by factors of a few. This mechanism, which we refer to as stimulated feedback, is motivated by recent ALMA observations of central galaxies in clusters and groups revealing molecular clouds apparently forming in the wakes of rising X-ray bubbles and with surprisingly low cloud velocities. Supported by recent numerical simulations, our model would naturally sustain a continual feedback-loop in galaxies fuelled by cooling gas stimulated by radio-mechanical feedback itself, that otherwise stabilizes cooling atmospheres on larger scales. The observed cooling time threshold for the onset of nebular emission and star formation of ~ 5x10^8 yr may result from the limited ability of radio bubbles to lift low entropy gas to an altitude where thermal instabilities can ensue. The molecular clouds condensing from the outflowing hot gas are unlikely to escape, but instead return to the central galaxy in a circulating flow.
Studies measuring the star formation rate density, luminosity function and properties of star-forming galaxies are numerous. However, it exists a gap at $0.5<z<0.8$ in H$\alpha$-based studies. Our main goal is to study the properties of a sample of faint H$\alpha$ emitters at $z\sim0.62$. We focus on their contribution to the faint-end of the luminosity function and derived star formation rate density, characterising their morphologies and basic photometric and spectroscopic properties. We use a narrow-band technique in the near-infrared, with a filter centered at 1.06 $\mu$m. The data come from ultra-deep VLT/HAWK-I observations in the GOODS-S field with a total of 31.9 h in the narrow-band filter. We perform a visual classification of the sample and study their morphologies from structural parameters available in CANDELS. Our 28 H$\alpha$-selected sample of faint star-forming galaxies reveals a robust faint-end slope of the luminosity function $\alpha=-1.46_{-0.08}^{+0.16}$. The derived star formation rate density at $z\sim0.62$ is $\rho_\mathrm{SFR} = 0.036_{-0.008}^{+0.012} M_{\odot}~\mathrm{yr^{-1}~Mpc^{-3}}$. The sample is mainly composed of disks, but an important contribution of compact galaxies with S\'ersic indexes $n\sim2$ display the highest specific star formation rates. The luminosity function at $z\sim0.62$ from our ultra-deep data points towards a steeper $\alpha$ when an individual extinction correction for each object is applied. Compact galaxies are low-mass, low-luminosity and starburst-dominated objects with a light profile in an intermediate stage from early to late-types.
We present near-IR (1.1-2.4 micron) position-position-velocity cubes of the 500-year-old Orion BN/KL explosive outflow with spatial resolution 1" and spectral resolution 86 km/s. We construct integrated intensity maps free of continuum sources of 15 H2 and [Fe II] lines while preserving kinematic information of individual outflow features. Included in the detected H2 lines are the 1-0 S(1) and 1-0 Q(3) transitions, allowing extinction measurements across the outflow. Additionally, we present dereddened flux ratios for over two dozen outflow features to allow for the characterization of the true excitation conditions of the BN/KL outflow. All ratios show the dominance of shock excitation of the H2 emission, although some features exhibit signs of fluorescent excitation from stellar radiation or J-type shocks. We also detect tracers of the PDR/ionization front north of the Trapezium stars in [O I] and [Fe II] and analyze other observed outflows not associated with the BN/KL outflow.
The Fundamental Plane (FP) describes the relation between stellar mass, size, and velocity dispersion of elliptical galaxies; the Faber-Jackson relation (FJR) is its projection onto {mass, velocity} space. In this work we redeploy and expand the framework of Desmond & Wechsler (2015) to ask whether abundance matching-based $\Lambda$CDM models that have shown success in matching the spatial distribution of galaxies are also capable of explaining key properties of the FJR and FP, including their scatter. Within our framework, agreement with the normalisation of the FJR requires haloes to expand in response to disc formation. We find that the tilt of the FP may be explained by a combination of the observed non-homology in galaxy structure and the variation in mass-to-light ratio produced by abundance matching with a universal IMF, provided that the anisotropy of stellar motions is taken into account. However, the predicted scatter around the FP is considerably increased by situating galaxies in cosmologically-motivated haloes due to variations in halo properties at fixed stellar mass, and appears to exceed that of the data. This implies that additional correlations between galaxy and halo variables may be required to fully reconcile these models with elliptical galaxy scaling relations.
I present an overview of the science goals and achievements of ongoing spectroscopic surveys of individual stars in the nearby Universe. I include a brief discussion of the development of the field of Galactic Archaeology - using the fossil record in old stars nearby to infer how our Galaxy evolved and place the Milky Way in cosmological context.
Nantais et al. used the Hubble Space Telescope to localize probable globular clusters (GCs) in M81, a spiral galaxy at a distance of 3.63 Mpc. Theory predicts that GCs can host intermediate-mass black holes (IMBHs) with masses M_BH \sim 100 - 100,000 M_sun. Finding IMBHs in GCs could validate a formation channel for seed BHs in the early universe, bolster gravitational-wave predictions for space missions, and test scaling relations between stellar systems and the central BHs they host. We used the NRAO Karl G. Jansky Very Large Array (VLA) to search for the radiative signatures of IMBH accretion from 206 probable GCs in a mosaic of M81. The observing wavelength was 5.5 cm and the spatial resolution was 1.5 arcsec (26.4 pc). None of the individual GCs are detected, nor are weighted-mean image stacks of the 206 GCs and the 49 massive GCs with stellar masses M_star \gtrsim 200,000 M_sun. We apply a semi-empirical model to predict the mass of an IMBH that, if undergoing accretion in the long-lived hard X-ray state, is consistent with a given radio luminosity. The 3$\sigma$ radio-luminosity upper limits correspond to mean IMBH masses of M_BH(all) < 42,000 M_sun for the all-cluster stack and M_BH(massive) < 51,000 M_sun for the massive-cluster stack. We also apply the empirical fundamental-plane relation to two X-ray-detected clusters, finding that their individual IMBH masses at 95% confidence are M_BH < 99,000 M_sun and M_BH < 15,000 M_sun. Finally, no analog of HLX-1, a strong IMBH candidate in an extragalactic star cluster, occurs in any individual GC in M81. This underscores the uniqueness or rarity of the HLX-1 phenomenon.
In this paper, by cross-correlating an archive sample of 542 extragalactic radio sources with the \emph{Fermi}-LAT Third Source Catalog(3FGL), we have compiled a sample of 80 $\gamma$-ray sources and 462 non-\emph{Fermi} sources with available core dominance parameter($R_{CD}$), core and extended radio luminosity; all the parameters are directly measured or derived from available data in the literature. We found that $R_{CD}$ have significant correlations with radio core luminosities, $\gamma$-ray luminosity and $\gamma$-ray flux respectively; the \emph{Fermi} sources have on average higher $R_{CD}$ than non-\emph{Fermi} sources. These results indicate that the \emph{Fermi} sources should be more compact, and beaming effect should play a crucial role for the detection of $\gamma$-ray emission. Moreover, our results also show \emph{Fermi} sources have systematically larger radio flux than non-\emph{Fermi} sources at fixed $R_{CD}$, indicating larger intrinsic radio flux in \emph{Fermi} sources. These results show a strong connection between radio and $\gamma$-ray flux for the present sample and indicate that the non-\emph{Fermi} sources is likely due to low beaming effect, and/or the low intrinsic $\gamma$-ray flux, support a scenario in the literature: a co-spatial origin of the activity for the radio and $\gamma$-ray emission, suggesting that the origin of the seed photons for the high-energy $\gamma$-ray emission is within the jet.
An abundance analysis is reported of 58 K giants identified by Famaey et al. (2005) as highly probable members of the Hercules stream selected from stars north of the celestial equator in the Hipparcos catalogue. The giants have compositions spanning the interval [Fe/H] from $-$0.17 to $+$0.42 with a mean value of $+$0.15 and relative elemental abundances [El/Fe] representative of the Galactic thin disc. Selection effects may have biassed the selection from the Hipparcos catalogue against the selection of metal-poor stars. Our reconsideration of the recent extensive survey of FG dwarfs which included metal-poor stars (Bensby et al. 2014) provides a [Fe/H] distribution for the Hercules stream which is similar to that from the 58 giants. It appears that the stream is dominated by metal-rich stars from the thin disc. Suggestions in the literature that the stream includes metal-poor stars from the thick disc are discussed.
The Santa Barbara cluster comparison project (Frenk et al. Frenk+1999) revealed that there is a systematic difference between entropy profiles of clusters of galaxies obtained by Eulerian mesh and Lagrangian smoothed particle hydrodynamics (SPH) codes: Mesh codes gave a core with a constant entropy whereas SPH codes did not. One possible reason for this difference is that mesh codes are not Galilean invariant. Another possible reason is the problem of the SPH method, which might give too much "protection" to cold clumps because of the unphysical surface tension induced at contact discontinuities. In this paper, we apply the density independent formulation of SPH (DISPH), which can handle contact discontinuities accurately, to simulations of a cluster of galaxies, and compare the results with those with the standard SPH. We obtained the entropy core when we adopt DISPH. The size of the core is, however, significantly smaller than those obtained with mesh simulations, and is comparable to those obtained with quasi-Lagrangian schemes such as "moving mesh" and "mesh free" schemes. We conclude that both the standard SPH without artificial conductivity and Eulerian mesh codes have serious problems even such an idealized simulation, while DISPH, SPH with artificial conductivity, and quasi-Lagrangian schemes have sufficient capability to deal with it.
We present the results of ALMA observations of dust continuum emission and molecular rotational lines, including the ALMA Compact Array, toward a dense core MC27 (a.k.a. L1521F) in Taurus, which is considered to be at very early stage of star formation. Detailed column density distribution with a size scale from a few tens AU to ~10000 AU scale are revealed by combining the ALMA data and the single-dish data. The high angular resolution observation at 0.87 mm reveals that a protostellar source, MMS-1, is still not spatially resolved without gas association and a starless high-density core, MMS-2, has substructures both in dust and molecular emission. The averaged radial column density distribution of the inner part (r < 3000 AU) is N_H2 ~r^-0.4, clearly flatter than that of the outer part, ~r^-1.0. We found the complex velocity/spatial structure obtained with previous ALMA observations is located inside the inner flatter region, which may reflect the dynamical status of the dense core.
We have performed systemmatic local linear stability analysis on a radially stratified infinite self-gravitating cylinder of rotating plasma under the influence of magnetic field. In order to render the system analytically tractable, we have focussed solely on the axisymmetric modes of perturbations. Using cylindrical coordinate system, we have derived the critical linear mass density of a non-rotating filament required for gravitational collapse to ensue in the presence of azimuthal magnetic field. Moreover, for such filaments threaded by axial magnetic field, we show that the growth rates of the modes having non-zero radial wavenumber are reduced more strongly by the magnetic field than that of the modes having zero radial wavenumber. More importantly, our study contributes to the understanding of the stability property of rotating astrophysical filaments that are more often than not influenced by magnetic fields. In addition to complementing many relevant numerical studies reported the literature, our results on filaments under the influence of magnetic field generalize some of the very recent analytical works (e.g.,~\citet{jog2014}, etc.). For example, here we prove that even a weak magnetic field can play a dominant role in determining stability of the filament when the rotation timescale is larger than the free fall timescale. A filamentary structure with faster rotation is, however, comparatively more stable for the same magnetic field. The results reported herein, due to strong locality assumption, are strictly valid for the modes for which one can ignore the radial variations in the density and the magnetic field profiles.
Aims. We aim to investigate whether a multilayer ice model can be as successful as a bulk ice model in reproducing the observed abundances of various deuterated gas-phase species toward starless cores. Methods. We calculate abundances for various deuterated species as functions of time adopting fixed physical conditions. We also estimate abundance gradients by adopting a modified Bonnor-Ebert sphere as a core model. In the multilayer ice scenario, we consider desorption from one or several monolayers on the surface. Results. We find that the multilayer model predicts abundances of $\rm DCO^+$ and $\rm N_2D^+$ that are about an order of magnitude lower than observed, caused by the trapping of CO and $\rm N_2$ into the grain mantle. As a result of the mantle trapping, deuteration efficiency in the gas phase increases and we find stronger deuterium fractionation in ammonia than what has been observed. Another distinguishing feature of the multilayer model is that $\rm D_3^+$ becomes the main deuterated ion at high density. The bulk ice model is generally easily reconciled with observations. Conclusions. Our results underline that more theoretical and experimental work is needed to understand the composition and morphology of interstellar ices, and the desorption processes that can act on them. With the current constraints, the bulk ice model appears to be better in reproducing observations than the multilayer ice model. According to our results, the $\rm H_2D^+$ to $\rm N_2D^+$ abundance ratio is higher than 100 in the multilayer model, while only a few $\times$ 10 in the bulk model, and so observations of this ratio could provide information on the ice morphology in starless cores. Observations of the abundance of $\rm D_3^+$ compared to $\rm H_2D^+$ and $\rm D_2H^+$ would provide additional constraints for the models.
Context. "Snow lines", marking regions where abundant volatiles freeze out onto the surface of dust grains, play an important role for planet growth and bulk composition in protoplanetary disks. They can already be observed in the envelopes of the much younger, low-mass Class 0 protostars that are still in their early phase of heavy accretion. Aims. We aim at using the information on the sublimation regions of different kinds of ices to understand the chemistry of the envelope, its temperature and density structure, and the history of the accretion process. Methods. As part of the CALYPSO IRAM Large Program, we have obtained observations of C$^{18}$O, N$_2$H$^+$ and CH$_3$OH towards nearby Class 0 protostars with the IRAM Plateau de Bure interferometer at sub-arcsecond resolution. For four of these sources we have modeled the emission using a chemical code coupled with a radiative transfer module. Results. We observe an anti-correlation of C$^{18}$O and N$_2$H$^+$ in NGC 1333-IRAS4A, NGC 1333-IRAS4B, L1157, and L1448C, with N$_2$H$^+$ forming a ring around the centrally peaked C$^{18}$O emission due to N$_2$H$^+$ being chemically destroyed by CO. The emission regions of models and observations match for a CO binding energy of 1200 K, which is higher than the binding energy of pure CO ices ($\sim$855 K). Furthermore, we find very low CO abundances inside the snow lines in our sources, about an order of magnitude lower than the total CO abundance observed in the gas on large scales in molecular clouds before depletion sets in. Conclusions. The high CO binding energy may hint at CO being frozen out in a polar ice environment like amorphous water ice or in non-polar CO$_2$-rich ice. The low CO abundances are comparable to values found in protoplanetary disks, which may indicate an evolutionary scenario where these low values are already established in the protostellar phase. (Abbr. Version)
This paper presents 1.4-GHz radio continuum observations of 15 very extended radio galaxies. These sources are so large that most interferometers lose partly their structure and total flux density. Therefore, single-dish detections are required to fill in the central (u,v) gap of interferometric data and obtain reliable spectral index patterns across the structures, and thus also an integrated radio continuum spectrum. We have obtained such 1.4-GHz maps with the 100-m Effelsberg telescope and combined them with the corresponding maps available from the NVSS. The aggregated data allow us to produce high-quality images, which can be used to obtain physical parameters of the mapped sources. The combined images reveal in many cases extended low surface-brightness cocoons.
Observationally, the quenching of star-forming galaxies appears to depend both on their mass and environment. The exact cause of the environmental dependence is still poorly understood, yet semi-analytic models (SAMs) of galaxy formation need to parameterise it to reproduce observations of galaxy properties. In this work, we use hydrodynamical simulations to investigate the quenching of disk galaxies through ram-pressure stripping (RPS) as they fall into galaxy clusters with the goal of characterising the importance of this effect for the reddening of disk galaxies. Our set-up employs a live model of a galaxy cluster that interacts with infalling disk galaxies on different orbits. We use the moving-mesh code AREPO, augmented with a special refinement strategy to yield high resolution around the galaxy on its way through the cluster in a computationally efficient way. Our direct simulations differ substantially from stripping models employed in current SAMs, which in most cases overpredict the mass loss from RPS. Furthermore, after pericentre passage, as soon as ram pressure becomes weaker, gas that remains bound to the galaxy is redistributed to the outer parts, an effect that is not captured in simplified treatements of RPS. Star formation in our model galaxies is quenched mainly because the hot gas halo is stripped, depriving the galaxy of its gas supply. The cold gas disk is only stripped completely in extreme cases, leading to full quenching and significant reddening on timescale of ~200 Myr. On the other hand, galaxies experiencing only mild ram pressure actually show an enhanced star formation rate that is consistent with observations and are quenched on timescales > 1 Gyr. Stripped gas in the wake is mixed efficiently with intracluster gas already a few tens of kpc behind the disk, and this gas is free of residual star formation.
We present an extended version of the 2-phase gas-grain code NAUTILUS to the 3-phase modelling of gas and grain chemistry of cold cores. In this model, both the mantle and the surface are considered as chemically active. We also take into account the competition among reaction, diffusion and evaporation. The model predictions are confronted to ice observations in the envelope of low-mass and massive young stellar objects as well as toward background stars. Modelled gas-phase abundances are compared to species observed toward TMC-1 (CP) and L134N dark clouds. We find that our model successfully reproduces the observed ice species. It is found that the reaction-diffusion competition strongly enhances reactions with barriers and more specifically reactions with H2, which is abundant on grains. This finding highlights the importance to have a good approach to determine the abundance of H2 on grains. Consequently, it is found that the major N-bearing species on grains go from NH3 to N2 and HCN when the reaction-diffusion competition is accounted. In the gas-phase and before few 10^5 yrs, we find that the 3-phase model does not have a strong impact on the observed species compared to the 2-phase model. After this time, the computed abundances dramatically decrease due to the strong accretion on dust, which is not counterbalanced by the desorption less efficient than in the 2-phase model. This strongly constrains the chemical-age of cold cores to be of the order of few 10^5 yrs.
A new approach to determining the solar galactocentric distance, $R_0$, from the geometry of spiral-arm segments is proposed. Geometric aspects of the problem are analyzed and a simplified three-point method for estimating $R_0$ from objects in a spiral segment is developed in order to test the proposed approach. An estimate of $R_0 = 8.44 \pm 0.45$ kpc is obtained by applying the method to masers with measured trigonometric parallaxes, and statistical properties of the $R_0$ estimation from spiral segments are analyzed.
We report the results of our follow-up campaign of the supernova impostor PSN J09132750+7627410, based on optical data covering $\sim250\,\rm{d}$. From the beginning, the transient shows prominent narrow Balmer lines with P-Cygni profiles, with a blue-shifted absorption component becoming more prominent with time. Along the $\sim3\,\rm{months}$ of the spectroscopic monitoring, broad components are never detected in the hydrogen lines, suggesting that these features are produced in slowly expanding material. The transient reaches an absolute magnitude $M_r=-13.60\pm0.19\,\rm{mag}$ at maximum, a typical luminosity for supernova impostors. Amateur astronomers provided $\sim4\,\rm{years}$ of archival observations of the host galaxy, NGC 2748. The detection of the quiescent progenitor star in archival images obtained with the Hubble Space Telescope suggests it to be an $18-20$\msun white-yellow supergiant.
We report the detection of near-IR H$_2$ lines emission from low-ionization structures (LISs) in planetary nebulae. The deepest, high-angular resolution H$_2$ 1-0 S(1) at 2.122 $\mu$m, and H$_2$ 2-1 S(1) at 2.248 $\mu$m images of K 4-47 and NGC 7662, obtained using NIRI@Gemini-North, are presented here. K 4-47 reveals a remarkable high-collimated bipolar structure. The H$_2$ emission emanates from the walls of the bipolar outflows and a pair of knots at the tips of these outflows. The H$_2$ 1-0 S(1)/2-1 S(1) line ratio is $\sim$7-10 which indicates shock interaction due to both the lateral expansion of the gas in the outflows and the high-velocity knots. The strongest line, H$_2$ v=1-0 S(1), is also detected in several LISs located at the periphery of the outer shell of the elliptical PN NGC 7662, whereas only four knots are detected in the H$_2$ v=2-1 S(1) line. These knots have H$_2$ v=1-0 S(1)/v=2-1 S(1) values between 2 and 3. These data confirm the presence of molecular gas in both highly (K 4-47) and slowly moving LISs (NGC 7662). The H$_2$ emission in K 4-47 is powered by shocks, whereas in NGC 7662 is due to photo-ionization by the central star. Moreover, a likely correlation is found between the H$_2$ v=1-0 S(1)/H$_2$ v=2-1 S(1) and [N II]/H$\alpha$ line ratios.
In order to understand the rate of merger of stellar-mass black hole binaries (BHBs) by gravitational wave (GW) emission it is important to determine the major pathways to merger. We use numerical simulations to explore the evolution of BHBs inside the radius of influence of supermassive black holes (SMBHs) in galactic centers. In this region the evolution of binaries is dominated by perturbations from the central SMBH. In particular, as first pointed out by Antonini and Perets, the Kozai-Lidov (KL) mechanism trades relative inclination of the BHB to the SMBH for eccentricity of the BHB, and for some orientations can bring the BHB to an eccentricity near unity. At very high eccentricities, GW emission from the BHB can become efficient, causing the members of the BHB to coalesce. We use a novel combination of two N-body codes to follow this evolution. We are forced to simulate small systems to follow the behavior accurately. We have completed 400 simulations that range from $\sim$ 300 stars around a $10^{3}$ M$_{\odot}$ black hole to $\sim$ 4500 stars around a $10^{4}$ M$_{\odot}$ black hole. These simulations are the first to follow the internal orbit of a binary near a SMBH while also following the changes to its external orbit self-consistently. We find that this mechanism could produce mergers at a maximum rate per volume of $\sim 100$ Gpc$^{-3}$ yr$^{-1}$ or considerably less if the inclination oscillations of the binary remain constant as the BHB inclination to the SMBH changes, or if the binary black hole fraction is small.
Galactic Globular clusters (GCs) are now known to harbour multiple stellar populations, which are chemically distinct in many light element abundances. It is becoming increasingly clear that asymptotic giant branch (AGB) stars in GCs show different abundance distributions in light elements compared to those in the red giant branch (RGB) and other phases, skewing toward more primordial, field-star-like abundances, which we refer to as subpopulation one (SP1). As part of a larger program targeting giants in GCs, we obtained high-resolution spectra for a sample of 106 RGB and 15 AGB stars in Messier 4 (NGC 6121) using the 2dF+HERMES facility on the Anglo-Australian Telescope. In this Letter we report an extreme paucity of AGB stars with [Na/O] > -0.17 in M4, which contrasts with the RGB that has abundances up to [Na/O] =0.55. The AGB abundance distribution is consistent with all AGB stars being from SP1. This result appears to imply that all subpopulation two stars (SP2; Na-rich, O-poor) avoid the AGB phase. This is an unexpected result given M4's horizontal branch morphology -- it does not have an extended blue horizontal branch. This is the first abundance study to be performed utilising the HERMES spectrograph.
Links to: arXiv, form interface, find, astro-ph, recent, 1604, contact, help (Access key information)
We examine the effect of an accretion disc on the orbits of stars in the central star cluster surrounding a central massive black hole by performing a suite of 39 high-accuracy direct N-body simulations using state-of-the art software and accelerator hardware, with particle numbers up to 128k. The primary focus is on the accretion rate of stars by the black hole (equivalent to their tidal disruption rate for black holes in the small to medium mass range) and the eccentricity distribution of these stars. Our simulations vary not only the particle number, but disc model (two models examined), spatial resolution at the centre (characterised by the numerical accretion radius) and softening length. The large parameter range and physically realistic modelling allow us for the first time to confidently extrapolate these results to real galactic centres. While in a real galactic centre both particle number and accretion radius differ by a few orders of magnitude from our models, which are constrained by numerical capability, we find that the stellar accretion rate converges for models with N > 32k. The eccentricity distribution of accreted stars, however, does not converge. We find that there are two competing effects at work when improving the resolution: larger particle number leads to a smaller fraction of stars accreted on nearly-circular orbits, while higher spatial resolution increases this fraction. We scale our simulations to some nearby galaxies and find that the expected boost in stellar accretion (or tidal disruption, which could be observed as X-ray flares) in the presence of a gas disc is about a factor of 10. Even with this boost, the accretion of mass from stars is still a factor of ~ 100 slower than the accretion of gas from the disc. Thus, it seems accretion of stars is not a major contributor to black hole mass growth.
This is the third paper in a series describing the spectroscopic properties of a sample of 39 AGN at $z \sim 1.5$, selected to cover a large range in black hole mass ($M_{BH}$) and Eddington ratio ($L/L_{Edd}$). In this paper, we continue the analysis of the VLT/X-shooter observations of our sample with the addition of 9 new sources. We use an improved Bayesian procedure, which takes into account intrinsic reddening, and improved $M_{BH}$ estimates, to fit thin accretion disc (AD) models to the observed spectra and constrain the spin parameter ($a_*$) of the central black holes. We can fit 37 out of 39 AGN with the thin AD model, and for those with satisfactory fits, we obtain constraints on the spin parameter of the BHs, with the constraints becoming generally less well defined with decreasing BH mass. Our spin parameter estimates range from $\sim$$-$0.6 to maximum spin for our sample, and our results are consistent with the "spin-up" scenario of BH spin evolution. We also discuss how the results of our analysis vary with the inclusion of non-simultaneous GALEX photometry in our thin AD fitting. Simultaneous spectra covering the rest-frame optical through far-UV are necessary to definitively test the thin AD theory and obtain the best constraints on the spin parameter.
We investigate the burstiness of star formation histories (SFHs) of galaxies at $0.4<z<1$ by using the ratio of star formation rates (SFRs) measured from FUV (1500 \AA) and H$\beta$ (FUV--to--H$\beta$ ratio). Our sample contains 164 galaxies down to stellar mass (M*) of $10^{8.5} M_\odot$ in the CANDELS GOODS-N region, where TKRS Keck/DEIMOS spectroscopy and HST/WFC3 F275W images from CANDELS and HDUV are available. When the ratio of FUV- and H$\beta$-derived SFRs is measured, dust extinction correction is negligible (except for very dusty galaxies) with the Calzetti attenuation curve. The FUV--to--H$\beta$ ratio of our sample increases with the decrease of M* and SFR. The median ratio is $\sim$1 at M* $\sim 10^{10} M_\odot$ (or SFR = 20 $M_\odot$/yr) and increases to $\sim$1.6 at M* $\sim 10^{8.5} M_\odot$ (or SFR $\sim 0.5 M_\odot$/yr). At M* $< 10^{9.5} M_\odot$, our median FUV--to--H$\beta$ ratio is higher than that of local galaxies at the same M*, implying a redshift evolution. Bursty SFH on a timescale of $\sim$10 Myr on galactic scales provides a plausible explanation of our results, and the importance of the burstiness increases as M* decreases. Due to sample selection effects, our FUV--to--H$\beta$ ratio may be a lower limit of the true value of a complete sample, which strengthens our conclusions. Other models, e.g., non-universal initial mass function or stochastic star formation on star cluster scales, are unable to plausibly explain our results.
We present the first study of evolution of galaxy groups in the Illustris simulation. We focus on dynamically relaxed and unrelaxed galaxy groups representing dynamically evolved and evolving galaxy systems, respectively. The evolutionary state of a group is probed from its luminosity gap and separation between the brightest group galaxy and the center of mass of the group members. We find that the Illustris simulation, over-produces large luminosity gap galaxy systems, known as fossil systems, in comparison to observations and the probed semi-analytical predictions. However, this simulation is equally successful in recovering the correlation between luminosity gap and luminosity centroid offset, in comparison to the probed semi-analytic model. We find evolutionary tracks based on luminosity gap which indicate that a large luminosity gap group is rooted in a small luminosity gap group, regardless of the position of the brightest group galaxy within the halo. This simulation helps, for the first time, to explore the black hole mass and its accretion rate in galaxy groups. For a given stellar mass of the brightest group galaxies, the black hole mass is larger in dynamically relaxed groups with a lower rate of mass accretion. We find this consistent with the latest observational studies of the radio activities in the brightest group galaxies in fossil groups. We also find that the IGM in dynamically evolved groups is hotter for a given halo mass than that in evolving groups, again consistent with earlier observational studies.
We present the first results from the EMPIRE survey, an IRAM large program that is mapping tracers of high density molecular gas across the disks of nine nearby star-forming galaxies. Here, we present new maps of the 3-mm transitions of HCN, HCO+, and HNC across the whole disk of our pilot target, M51. As expected, dense gas correlates with tracers of recent star formation, filling the "luminosity gap" between Galactic cores and whole galaxies. In detail, we show that both the fraction of gas that is dense, f_dense traced by HCN/CO, and the rate at which dense gas forms stars, SFE_dense traced by IR/HCN, depend on environment in the galaxy. The sense of the dependence is that high surface density, high molecular gas fraction regions of the galaxy show high dense gas fractions and low dense gas star formation efficiencies. This agrees with recent results for individual pointings by Usero et al. 2015 but using unbiased whole-galaxy maps. It also agrees qualitatively with the behavior observed contrasting our own Solar Neighborhood with the central regions of the Milky Way. The sense of the trends can be explained if the dense gas fraction tracks interstellar pressure but star formation occurs only in regions of high density contrast.
We show influences of the mass of seed black holes on black hole mass -- bulge mass relation at z ~ 0 by using a semi-analytic model of galaxy formation combined with large cosmological N-body simulations. We constrain our model to reproduce observed properties of galaxies at z ~ 0. Similar to other semi-analytic models, we place a seed black hole immediately after a galaxy forms. When we set the seed black hole mass to 10^5 M_sun, we find that the model result becomes inconsistent with recent observational results of black hole mass -- bulge mass relation for dwarf galaxies. Namely, the model predicts that bulges with ~ 10^9 M_sun harbor black holes more massive than observed. On the other hand, when we employ seed black holes with 10^3 M_sun or randomly choose their masses in the range of 10^{3-5} M_sun, the black hole mass -- bulge mass relation obtained from these models are consistent with observational results including dispersions. We find that to obtain more stringent restrictions of the mass of seed black holes, observations of less massive bulges at z ~ 0 are more powerful than that of black hole mass -- bulge mass relations at higher redshifts.
To study the strength and structure of the magnetic field in the Galactic centre (GC) we measured Faraday rotation of the radio emission of pulsars which are seen towards the GC. Three of these pulsars have the largest rotation measures (RMs) observed in any Galactic object with the exception of Sgr A*. Their large dispersion measures, RMs and the large RM variation between these pulsars and other known objects in the GC implies that the pulsars lie in the GC and are not merely seen in projection towards the GC. The large RMs of these pulsars indicate large line-of-sight magnetic field components between ~ 16-33 microgauss; combined with recent model predictions for the strength of the magnetic field in the GC this implies that the large-scale magnetic field has a very small inclination angle with respect to the plane of the sky (~ 12 degrees). Foreground objects like the Radio Arc or possibly an ablated, ionized halo around the molecular cloud G0.11-0.11 could contribute to the large RMs of two of the pulsars. If these pulsars lie behind the Radio Arc or G0.11-0.11 then this proves that low-scattering corridors with lengths >~ 100 pc must exist in the GC. This also suggests that future, sensitive observations will be able to detect additional pulsars in the GC. Finally, we show that the GC component in our most accurate electron density model oversimplifies structure in the GC.
[Abridged] We exploit the synergy between low-resolution spectroscopy and photometric redshifts to study environmental effects on galaxy evolution in slitless spectroscopic surveys from space. As a test case, we consider the future Euclid Deep survey (~40deg$^2$), which combines a slitless spectroscopic survey limited at H$\alpha$ flux $\leq5\times 10^{-17}$ erg cm$^{-2}$ s$^{-1}$ and a photometric survey limited in H-band ($H\leq26$). To test the power of the method, we use Euclid-like galaxy mock catalogues, in which we anchor the photometric redshifts to the 3D galaxy distribution of the available spectroscopic redshifts. We then estimate the local density contrast by counting objects in cylindrical cells with radius ranging from 1 to 10 h$^{-1}$Mpc over the redshift range 0.9<z<1.8. We compare this density field with the one computed in a mock catalogue with the same depth as the Euclid Deep survey (H=26) but without redshift measurement errors. We find that our method is successful in separating high from low density environments with an efficiency that increases at low redshift and large cells. The fraction of low density regions mistaken by high density peaks is below 1% for all scales and redshift explored, but for scales of 1 h$^{-1}$Mpc for which is a few percent. When small (1 h$^{-1}$Mpc) cells are used, our technique is successful, at z~1.5, at spotting the regions where the most massive galaxy clusters reside. These results show that we can efficiently study environment in photometric samples if spectroscopic information is available for a smaller sample of objects that sparsely samples the same volume. We demonstrate that these studies will be possible in the Euclid Deep survey, i.e. in a redshift range in which environmental effects are different from those observed in the local universe, hence providing new constraints for galaxy evolution models.
At redshifts beyond $z{\sim}1$ measuring the black hole galaxy relations proves to be a difficult task. The bright light of the AGN aggravates deconvolution of black hole and galaxy properties. On the other hand high redshift data on these relations is vital to understand in what ways galaxies and black holes co-evolve and in what ways they don't. In this work we use black hole (BHMDs) and stellar mass densities (SMDs) to constrain the possible co-evolution of black holes with their host galaxies since $z{\sim}5$. The BHMDs are calculated from quasar luminosity functions (QLF) using the Soltan argument, while we use integrals over stellar mass functions (SMFs) or the star formation rate density to obtain values for the stellar mass density. We find that both quantities grow in lock-step below redshifts of $z{\sim}3$ with a non-evolving BHMD to SMD ratio. A fit to the data assuming a power law relation between the BHMD and the SMD yields exponents around unity ($1.0{-}1.5$). Up to $z{\sim}5$ the BHMD to SMD ratio doesn't show a strong evolution given the larger uncertainty in the completeness of high-redshift datasets. Our results, always applying the same analysis technique, seem to be consistent across all adopted data sets.
By combining the data from surveys for HI 21-cm absorption at various impact parameters in near-by galaxies, we report an anti-correlation between the 21-cm absorption strength (velocity integrated optical depth) and the impact parameter. Also, by combining the 21-cm absorption strength with that of the emission, giving the neutral hydrogen column density, we find no evidence that the spin temperature of the gas (degenerate with the covering factor) varies significantly across the disk. This is consistent with the uniformity of spin temperature measured across the Galactic disk. Furthermore, comparison with the Galactic distribution suggests that intervening 21-cm absorption preferentially arises in disks of high inclinations (near face-on). We also investigate the hypothesis that 21-cm absorption is favourably detected towards compact radio sources. Although there is insufficient data to determine whether there is a higher detection rate towards quasar, rather than radio galaxy, sight-lines, the 21-cm detections intervene objects with a mean turnover frequency of 5 x 10^8 Hz, compared to 1 x 10^8 Hz for the non-detections. Since the turnover frequency is anti-correlated with radio source size, this does indicate a preferential bias for detection towards compact background radio sources.
In the density wave theory of spiral structure, the grand-design two-armed spiral pattern is taken to rotate rigidly in a galactic disc with a constant, definite pattern speed. The observational measurement of the pattern speed of the spiral arms, though difficult, has been achieved in a few galaxies such as NGC 6946, NGC 2997, and M 51 which we consider here. We examine whether the theoretical dispersion relation permits a real solution for wavenumber corresponding to a stable wave, for the observed rotation curve and the pattern speed values. We find that the disc when treated to consist of stars alone, as is usually done in literature, does not generally support a stable density wave for the observed pattern speed. Instead the inclusion of the low velocity dispersion component, namely, gas, is essential to obtain a stable density wave. Further, we obtain a theoretical range of allowed pattern speeds that correspond to a stable density wave at a certain radius, and check that for the three galaxies considered, the observed pattern speeds fall in the respective prescribed range. The inclusion of even a small amount (~ 15 %) of gas by mass fraction in the galactic disc is shown to have a significant dynamical effect on the dispersion relation and hence on the pattern speed that is likely to be seen in a real, gas-rich spiral galaxy.
The observed structure function (SF) of rotation measure (RM) varies as a broken power-law function of angular scales. The systematic shallowness of its spectral slope is inconsistent with the standard Kolmogorov scaling. This motivates us to examine the statistical analysis on RM fluctuations. The correlations of RM constructed by Lazarian & Pogosyan (2016) are demonstrated to be adequate in explaining the observed features of RM SFs through a direct comparison between the theoretically obtained and observationally measured SF results. By segregating the density and magnetic field fluctuations and adopting arbitrary indices for their respective power spectra, we find that when the SFs of RM and emission measure have a similar form over the same range of angular scales, the statistics of the RM fluctuations reflect the properties of density fluctuations. RM SFs can be used to evaluate the mean magnetic field along the line of sight, but cannot serve as an informative source on the properties of turbulent magnetic field in the interstellar medium. We identify the spectral break of RM SFs as the inner scale of a shallow spectrum of electron density fluctuations, which characterizes the typical size of discrete electron density structures in the observed region.
We analyze the effects of flattening on the annihilation (J) and decay (D) factors of dwarf spheroidal galaxies with both analytic and numerical methods. Flattening has two consequences: first, there is a geometric effect as the squeezing (or stretching) of the dark matter distribution enhances (or diminishes) the J-factor; second, the line of sight velocity dispersion of stars must hold up the flattened baryonic component in the flattened dark matter halo. We provide analytic formulae and a simple numerical approach to estimate the correction to the J- and D-factors required over simple spherical modeling. The formulae are validated with a series of equilibrium models of flattened stellar distributions embedded in flattened dark-matter distributions. We compute corrections to the J- and D-factors for the Milky Way dwarf spheroidal galaxies under the assumption that they are prolate or oblate and find that the hierarchy of J-factors for the dwarf spheroidals is slightly altered. We demonstrate that spherical estimates of the D-factors are very insensitive to the flattening and introduce uncertainties significantly less than the uncertainties in the D-factors from the other observables for all the dwarf spheroidals. We conclude by investigating the spread in correction factors produced by triaxial figures and provide uncertainties in the J-factors for the dwarf spheroidals using different physically-motivated assumptions for their intrinsic shape and axis alignments. We find that the uncertainty in the J-factors due to triaxiality increases with the observed ellipticity and, in general, introduces uncertainties of order 25 per cent in the J-factors. We discuss our results in light of the reported gamma-ray annihilation signal from the highly-flattened ultrafaint Reticulum II.
Based on the analysis of available published data and archival data along 24 sightlines (5 of which are new) we derive more accurate estimates of the column densities of OH and CH towards diffuse/translucent clouds and revisit the typically observed correlation between the abundances of these species. The increase in the sample size was possible because of the equivalence of the column densities of CH derived from a combination of the transitions at 3137 & 3143 Angstrom, and a combination of transitions at 3886 & 3890 Angstrom, which we have demonstrated here. We find that with the exception of four diffuse clouds, the entire source sample shows a clear correlation between the column densities of OH and CH similar to previous observations. The analysis presented also verifies the theoretically predicted oscillator strengths of the OH A--X (3078 & 3082 Angstrom), CH B--X (3886 & 3890 Angstrom) and C--X (3137 & 3143 Angstrom) transitions. We estimate N(H) and N(H2) from the observed E(B-V) and N(CH) respectively. The N(OH)/N(CH) ratio is not correlated with the molecular fraction of hydrogen in the diffuse/translucent clouds. We show that with the exception of HD 34078 for all the clouds the observed column density ratios of CH and OH can be reproduced by simple chemical models which include gas-grain interaction and gas-phase chemistry. The enhanced N(OH)/N(CH) ratio seen towards the 3 new sightlines can be reproduced primarily by considering different cosmic ray ionization rates.
Observational implications are derived for two standard models of supernovae-driven galactic winds: a freely expanding steady-state wind and a wind sourced by a self-similarly expanding superbubble including thermal heat conduction. It is shown that, for the steady-state wind, matching the measured correlation between the soft x-ray luminosity and star formation rate of starburst galaxies is equivalent to producing a scaled wind mass-loading factor relative to the star-formation rate of 0.5 - 3, in agreement with the amount inferred from metal absorption line measurements. The match requires the asymptotic wind velocity v_inf to scale with the star formation rate SFR (in solar masses per year) approximately as v_inf ~ (700 - 1000) km/s SFR^{1/6}. The corresponding mass injection rate is close to the amount naturally provided by thermal evaporation from the wall of a superbubble in a galactic disc, suggesting thermal evaporation may be a major source of mass-loading. The predicted mass-loading factors from thermal evaporation within the galactic disc alone, however, are somewhat smaller, 0.2 - 2, so that a further contribution from cloud ablation or evaporation may be required. Both models may account for the 1.4GHz luminosity of unresolved radio sources within starburst galaxies for plausible parameters describing the distribution of relativistic electrons. Further observational tests to distinguish the models are suggested.
The Planck Collaboration has recently released maps of the microwave sky in both temperature and polarization. Diffuse astrophysical components (including Galactic emissions, cosmic far infrared (IR) background, y-maps of the thermal Sunyaev-Zeldovich (SZ) effect) and catalogs of many thousands of Galactic and extragalactic radio and far-IR sources, and galaxy clusters detected through the SZ effect are the main astrophysical products of the mission. A concise overview of these results and of astrophysical studies based on Planck data is presented.
We report two new tidal debris nearby the Sagittarius (Sgr) tidal stream in the north Galactic cap identified from the M giant stars in LAMOST DR2 data. The M giant stars with sky area of $210^\circ<$\Lambda$<290^\circ$, distance of 10--20kpc, and [Fe/H]$<-0.75$ show clear bimodality in velocity distribution. We denote the two peaks as Vel-3+83 for the one within mean velocity of -3kms$^{-1}$ with respect to that of the well observed Sgr leading tail at the same $\Lambda$ and Vel+162+26 for the other one with mean velocity of 162kms$^{-1}$ with respect to the Sgr leading tail. Although the projected $\Lambda$--$V_{gar}$ relation of Vel-3+83 is very similar to the Sgr leading tail, the opposite trend in $\Lambda$--distance relation against the Sgr leading tail suggests Vel-3+83 has a different 3D direction of motion with any branch of the simulated Sgr tidal stream from Law & Majewski. Therefore, we propose it to be a new tidal debris not related to the Sgr stream. Similarly, the another substructure Vel+162+26, which is the same one as the NGC group discovered by Chou et al., also moves toward a different direction with the Sgr stream, implying that it may have different origin with the Sgr tidal stream.
J-factors (or D-factors) describe the distribution of dark matter in an astrophysical system and determine the strength of the signal provided by annihilating (or decaying) dark matter respectively. We provide simple analytic formulae to calculate the J-factors for spherical cusps obeying the empirical relationship between enclosed mass, velocity dispersion and half-light radius. We extend the calculation to the spherical Navarro-Frenk-White (NFW) model, and demonstrate that our new formulae give accurate results in comparison to more elaborate Jeans models driven by Markov Chain Monte Carlo methods. Of the known ultrafaint dwarf spheroidals, we show that Ursa Major II, Reticulum II, Tucana II and Horologium I have the largest J-factors and so provide the most promising candidates for indirect dark matter detection experiments. Amongst the classical dwarfs, Draco, Sculptor and Ursa Minor have the highest J-factors. We show that the behaviour of the J-factor as a function of integration angle can be inferred for general dark halo models with inner slope $\gamma$ and outer slope $\beta$. The central and asymptotic behaviour of the J-factor curves are derived as a function of the dark halo properties. Finally, we show that models obeying the empirical relation on enclosed mass and velocity dispersion have J-factors that are most robust at the integration angle equal to the projected half-light radius of the dSph divided by heliocentric distance. For most of our results, we give the extension to the D-factor which is appropriate for the decaying dark matter picture.
Ammonia and its deuterated isotopologues probe physical conditions in dense molecular cloud cores. With the aim of testing the current understanding of the spin-state chemistry of these molecules, we observed spectral lines of NH3, NH2D, NHD2 , ND3 , and N2D+ towards a dense, starless core in Ophiuchus with the APEX, GBT, and IRAM 30-m telescopes. The observations were interpreted using a gas-grain chemistry model combined with radiative transfer calculations. The chemistry model distinguishes between the different nuclear spin states of light hydrogen molecules, ammonia, and their deuterated forms. High deuterium fractionation ratios with NH2D/NH3=0.5, NHD2/NH2D=0.2, and ND3/NHD2=0.07 were found in the core. The observed ortho/para ratios of NH2D and NHD2 are close to the corresponding nuclear spin statistical weights. The chemistry model can approximately reproduce the observed abundances, but a clearly better agreement with the observations is found using the assumption of statistical fractionation and spin ratios, which implies NH2D/NH3 = 3 NHD2/NH2D = 9 ND3/NHD2. The longevity of N2H+ and NH3 in dense gas, which is prerequisite to their strong deuteration, can be attributed to the chemical inertia of N2 on grain surfaces. The discrepancies between the chemistry model and the observations are likely to be caused by the fact that the model assumes complete scrambling in principal gas-phase deuteration reactions of ammonia, which means that all the nuclei are mixed in reactive collisions. If, instead, these reactions occur through proton hop/hydrogen abstraction processes, statistical spin ratios are to be expected.
We have performed a spectral decomposition to search for recoiling supermassive black holes (rSMBH) in the SDSS QSOs with $z<0.25$. Out of 1271 QSOs, we have identified 26 rSMBH candidates that are recoiling toward us. The projected recoil velocities range from $-76\ \kms$ to $-307\ \kms$ with a mean of $-149\pm58\ \kms$. Most of the rSMBH candidates are hosted by gas-rich LIRGs/ULIRGs, but only 23\% of them shows signs of tidal features suggesting majority of them are advanced mergers. We find that the black hole masses $M_{BH}$ of the rSMBH candidates are on average $\sim$5 times smaller than that of their stationary counterparts and cause a scatter in $M_{BH}-\sigma_*$ relation. The Eddington ratios of all of the rSMBH candidates are larger than 0.1, with mean of 0.52$\pm$0.27, suggesting they are actively accreting mass. Velocity shifts in high-excitation coronal lines suggest that the rSMBH candidates are recoiling with an average velocity of about $-265\ \kms$. Electron density in the narrow line region of the H II rSMBH candidates is about 1/10 of that in AGN rSMBH candidates probably because AGN in the former was more spatially offset than that in the latter. The estimated spatial offsets between the rSMBH candidate and center of host galaxy range from 0.21\as \ to 1.97\as \ and need to be confirmed spatially with high-resolution adaptive optics imaging observations.
We develop a detailed chemical network relevant to the conditions characteristic of prestellar core collapse. We solve the system of time-dependent differential equations to calculate the equilibrium abundances of molecules and dust grains, with a size distribution given by size-bins for these latter. These abundances are used to compute the different non-ideal magneto-hydrodynamics resistivities (ambipolar, Ohmic and Hall), needed to carry out simulations of protostellar collapse. For the first time in this context, we take into account the evaporation of the grains, the thermal ionisation of Potassium, Sodium and Hydrogen at high temperature, and the thermionic emission of grains in the chemical network, and we explore the impact of various cosmic ray ionisation rates. All these processes significantly affect the non-ideal magneto-hydrodynamics resistivities, which will modify the dynamics of the collapse. Ambipolar diffusion and Hall effect dominate at low densities, up to n_H = 10^12 cm^-3, after which Ohmic diffusion takes over. We find that the time-scale needed to reach chemical equilibrium is always shorter than the typical dynamical (free fall) one. This allows us to build a large, multi-dimensional multi-species equilibrium abundance table over a large temperature, density and ionisation rate ranges. This table, which we make accessible to the community, is used during first and second prestellar core collapse calculations to compute the non-ideal magneto-hydrodynamics resistivities, yielding a consistent dynamical-chemical description of this process.
Recent studies show the importance of feedback in the evolution of the star formation rate in the Universe. However, the nature and physics of the feedback are still pressing questions. Radio continuum observations can provide unique dust-unbiased tracers of massive star formation and of the interstellar medium (ISM) and hence are ideal to address the regulation of star formation in galaxies. Our multi-frequency and multi-resolution radio surveys in nearby galaxies enable us to trace various phases of star formation and dissect the thermal and nonthermal ISM in galaxies. Mapping the cosmic ray electron energy index and magnetic field strength, we have found observational evidence that massive star formation significantly affects the energy balance in the ISM through the injection and acceleration of cosmic rays and the amplification of magnetic fields. How the next generation of stars could form in such a magnetized and turbulent ISM will be addressed in our 'EVLA cloud-scale survey of the local group galaxy M33' and in forthcoming surveys with the SKA.
We present the first attempt to detect outflows from galaxies approaching the Epoch of Reionization (EoR) using a sample of 9 star-forming (5 < SFR < 70 Msun/yr) z ~ 6 galaxies for which high-quality spectra of the [CII]158 micron line has been previously obtained with ALMA. We first fit each line with a Gaussian function and compute the residuals by subtracting the best fitting model from the data. We combine the residuals of all sample galaxies and find that the total signal is characterized by a flux excess that can be ascribed to broad wings of the [CII] line, which we interpret as a signature of starburst-driven outflows. The tentatively inferred outflow rate is dM/dt ~ 65 Msun/yr. Our interpretation is consistent with results from zoomed hydro- simulations of Dahlia, a z ~ 7 galaxy (SFR ~ 100 Msun/yr) whose feedback-regulated star formation results in an outflow rate dM/dt ~ 30 Msun/yr. These results suggest that starburst-driven outflows are in place in the EoR. Deeper observations of the [CII] line in the galaxies of this sample are required to better characterize stellar feedback at high-z and to understand the role of outflows in shaping early galaxy formation.
\textbf{Aims}: We present POLARIS (\textbf{POLA}rized \textbf{R}ad\textbf{I}ation \textbf{S}imulator), a newly developed three-dimensional Monte-Carlo radiative transfer code. POLARIS was designed to calculate dust temperature, polarization maps, and spectral energy distributions. It is optimized to handle data that results from sophisticated magneto-hydrodynamic simulations. The main purpose of the code is to prepare and analyze multi-wavelength continuum polarization measurements in the context of magnetic field studies in the interstellar medium. An exemplary application is the investigation of the role of magnetic fields in star formation processes.\\ \textbf{Methods}: We combine currently discussed state-of-the-art grain alignment theories with existing dust heating and polarization algorithms. We test the POLARIS code on multiple scales in complex astrophysical systems that are associated with different stages of star formation. POLARIS uses the full spectrum of dust polarization mechanisms to trace the underlying magnetic field morphology.\\ \textbf{Results}: Resulting temperature distributions are consistent with the density and position of radiation sources resulting from magneto-hydrodynamic (MHD) - collapse simulations. The calculated layers of aligned dust grains in the considered cirumstellar disk models are in excellent agreement with theoretical predictions. Finally, we compute unique patterns in synthetic multi-wavelength polarization maps that are dependent on applied dust-model and grain-alignment theory in analytical cloud models.
Radio-mode active galactic nucleus (AGN) feedback plays a key role in the evolution of galaxy groups and clusters. Its physical origin lies in the kpc-scale interaction of AGN jets with the hot halo gas, where jet properties may play an important role. Large-scale jet simulations often initiate light internally-supersonic jets with density contrast $0.01<\eta<1$. Here we argue for the importance of very-light ($\eta<0.01$) internally-subsonic jets in AGN feedback. We investigated the shapes of young X-ray cavities produced by AGN jets in a suite of hydrodynamic simulations, and found that bottom-wide cavities are always produced by internally-subsonic jets, while internally-supersonic jets produce cylindrical, center-wide, or top-wide cavities. We found examples of real cavities inflated by internally-subsonic and internally-supersonic jets, suggesting a dichotomy of AGN jets according to their internal Mach numbers. We further studied the long-term cavity evolution, and found that old cavities resulted from light jets tend to spread along the jet direction, while those produced by very light jets are significantly elongated along the perpendicular direction. The northwestern ghost cavity in Perseus is pancake-shaped, providing tentative evidence for the existence of very light jets. We speculate that very-light internally-subsonic jets exist in galaxy clusters. Our simulations show that they decelerate faster and rise much slower in the ICM than light internally-supersonic jets, possibly depositing a larger fraction of jet energy to cluster cores and alleviating the problem of low coupling efficiencies found in previous studies. The internal Mach number points to the energy content of AGN jets, and internally-subsonic jets are energetically dominated by non-kinetic energy, such as thermal energy, cosmic rays, or magnetic fields.
Angular momentum distribution and its redistribution are of crucial importance in formation and evolution of circumstellar disks. Many molecular line observations toward young stellar objects indicate that radial distributions of the specific angular momentum $j$ are more or less universal. In small scales, typically $R\lesssim$ 100 AU, the specific angular momenta distribute like $j\propto r^{1/2}$, indicating existence of rotationally supported disk. In outer regions, $R\gtrsim$ 5000 AU, $j$ increases as the radius increases and the slope is steeper than unity, which is supposed to reflect the original angular momentum distributions in the maternal molecular clouds. And lastly there is a connecting region, 100 AU $\lesssim R \lesssim$ 5000 AU, in which $j$-distribution looks almost flat. While this is often interpreted as a consequence of conservation of the specific angular momentum, it actually is insufficient and requires a stronger condition that the initial distribution of $j$ must be spatially uniform. However, this requirement is unrealistic and inconsistent with observations. In this work, we propose a simple alternative explanation; the flat $j$ profile is produced by strong prolongation due to the large velocity gradient in the accreting flow no matter what the initial $j$ distribution is. We provide a simple analytic model for gravitational collapse of molecular clouds. This model can be used to estimate ages of protostars based solely on the observed rotation profile. We demonstrate its validity in comparison with hydrodynamic simulations, and apply the model to young stellar objects such as L1527 IRS, TMC-1A and B335.
We report the results of VERA multi-epoch VLBI 22 GHz water maser observations of S255IR-SMA1, a massive young stellar object located in the S255 star forming region. By annual parallax the source distance was measured as D = 1.78 +-0.12 kpc and the source systemic motion was (u alpha cos d, u d) = (-0.13 +- 0.20, -0.06 +- 0.27) mas yr-1. Masers appear to trace a U-shaped bow shock whose morphology and proper motions are well reproduced by a jet-driven outflow model with a jet radius of about 6 AU. The maser data, in the context of other works in the literature, reveal ejections from S255IR-SMA1 to be episodic, operating on timescales of ~1000 years.
Classical Cepheid variable stars are high-sensitivity probes of stellar evolution and fundamental tracers of cosmic distances. While rotational mixing significantly affects the evolution of Cepheid progenitors (intermediate-mass stars), the impact of the resulting changes in stellar structure and composition on Cepheids on their pulsational properties is hitherto unknown. Here we present the first detailed pulsational instability analysis of stellar evolution models that include the effects of rotation, for both fundamental mode and first overtone pulsation. We employ Geneva evolution models spanning a three-dimensional grid in mass (1.7 - 15 $M_\odot$), metallicity (Z = 0.014, 0.006, 0.002), and rotation (non-rotating, average & fast rotation). We determine (1) hot and cool instability strip (IS) boundaries taking into account the coupling between convection and pulsation, (2) pulsation periods, and (3) rates of period change. We investigate relations between period and (a) luminosity, (b) age, (c) radius, (d) temperature, (e) rate of period change, (f) mass, (g) the flux-weighted gravity-luminosity relation (FWGLR). We confront all predictions aside from those for age with observations, finding generally excellent agreement. We tabulate period-luminosity relations (PLRs) for several photometric pass-bands and investigate how the finite IS width, different IS crossings, metallicity, and rotation affect PLRs. We show that a Wesenheit index based on H, V, and I photometry is expected to have the smallest intrinsic PLR dispersion. We confirm that rotation resolves the Cepheid mass discrepancy. Period-age relations depend significantly on rotation (rotation increases Cepheid ages), offering a straightforward explanation for evolved stars in binary systems that cannot be matched by conventional isochrones assuming a single age. Finally, we show that Cepheids obey a tight FWGLR.
Links to: arXiv, form interface, find, astro-ph, recent, 1604, contact, help (Access key information)
We present the first detection of molecular emission from a galaxy selected to be near a projected background quasar using the Atacama Large Millimeter/submillimeter Array (ALMA). The ALMA detection of CO(1$-$0) emission from the $z=0.101$ galaxy toward quasar PKS 0439-433 is coincident with its stellar disk and yields a molecular gas mass of $M_{\rm mol} \approx 4.2 \times 10^9 M_\odot$ (for a Galactic CO-to-H$_2$ conversion factor), larger than the upper limit on its atomic gas mass. We resolve the CO velocity field, obtaining a rotational velocity of $134 \pm 11$ km s$^{-1}$, and a resultant dynamical mass of $\geq 4 \times 10^{10} M_\odot$. Despite its high metallicity and large molecular mass, the $z=0.101$ galaxy has a low star formation rate, implying a large gas consumption timescale, larger than that typical of late-type galaxies. Most of the molecular gas is hence likely to be in a diffuse extended phase, rather than in dense molecular clouds. By combining the results of emission and absorption studies, we find that the strongest molecular absorption component toward the quasar cannot arise from the molecular disk, but is likely to arise from diffuse gas in the galaxy's circumgalactic medium. Our results emphasize the potential of combining molecular and stellar emission line studies with optical absorption line studies to achieve a more complete picture of the gas within and surrounding high-redshift galaxies.
The South Pole Telescope has discovered one hundred gravitationally lensed, high-redshift, dusty, star-forming galaxies (DSFGs). We present 0.5" resolution 870um Atacama Large Millimeter/submillimeter Array imaging of a sample of 47 DSFGs spanning z=1.9-5.7, and construct gravitational lens models of these sources. Our visibility-based lens modeling incorporates several sources of residual interferometric calibration uncertainty, allowing us to properly account for noise in the observations. At least 70% of the sources are strongly lensed by foreground galaxies (mu_870um > 2), with a median magnification mu_870um = 6.3, extending to mu_870um > 30. We compare the intrinsic size distribution of the strongly lensed sources to a similar number of unlensed DSFGs and find no significant differences in spite of a bias between the magnification and intrinsic source size. This may indicate that the true size distribution of DSFGs is relatively narrow. We use the source sizes to constrain the wavelength at which the dust optical depth is unity and find this wavelength to be correlated with the dust temperature. This correlation leads to discrepancies in dust mass estimates of a factor of 2 compared to estimates using a single value for this wavelength. We investigate the relationship between the [CII] line and the far-infrared luminosity and find that the same correlation between the [CII]L_FIR ratio and Sigma_FIR found for low-redshift star-forming galaxies applies to high-redshift galaxies and extends at least two orders of magnitude higher in Sigma_FIR. This lends further credence to the claim that the compactness of the IR-emitting region is the controlling parameter in establishing the "[CII] deficit."
We present Gemini/GMOS imaging of twelve candidate intergalactic globular clusters (IGCs) in the Local Group, identified in a recent survey of the SDSS footprint by di Tullio Zinn & Zinn (2015). Our image quality is sufficiently high, at $\sim 0.4^{\prime\prime} - 0.7^{\prime\prime}$, that we are able to unambiguously classify all twelve targets as distant galaxies. To reinforce this conclusion we use GMOS images of globular clusters in the M31 halo, taken under very similar conditions, to show that any genuine clusters in the putative IGC sample would be straightforward to distinguish. Based on the stated sensitivity of the di Tullio Zinn & Zinn (2015) search algorithm, we conclude that there cannot be a significant number of IGCs with $M_V \le -6$ lying unseen in the SDSS area if their properties mirror those of globular clusters in the outskirts of M31 -- even a population of $4$ would have only a $\approx 1\%$ chance of non-detection.
Rest-frame UV-optical (i.e., NUV-B) color index is sensitive to the low-level recent star formation and dust extinction, but is insensitive to the metallicity. In this Letter, we have measured the rest-frame NUV-B color gradients in ~1400 large ($\rm r_e>0.18^{\prime\prime}$), nearly face-on (b/a>0.5) main-sequence star-forming galaxies (SFGs) between redshift 0.5 and 1.5 in the CANDELS/GOODS-S and UDS fields. With this sample, we study the origin of UV-optical color gradients in the SFGs at z~1 and discuss their link with the buildup of stellar mass. We find that more massive, centrally compact and more dust extinguished SFGs tend to have statistically more negative raw color gradients (redder centers) than less massive, centrally diffuse and less dusty SFGs. After correcting for dust reddening based on optical-SED fitting, the color gradients in the low-mass ($M_{\ast} <10^{10}M_{\odot}$) SFGs generally become quite flat, while most of the high-mass ($M_{\ast} > 10^{10.5}M_{\odot}$) SFGs still retain shallow negative color gradients. These findings imply that dust reddening is likely the principal cause of negative color gradients in the low-mass SFGs, while both increased central dust reddening and buildup of compact, old bulges are likely the origins of negative color gradients in the high-mass SFGs. These findings also imply that at these redshifts the low-mass SFGs build up their stellar masses in a self-similar way, while the high-mass SFGs grow inside-out.
The massive molecular outflows erupting from the high-mass young stellar objects provide important clues to understand the mechanism of the high-mass star formation. Based on new CO J=3-2 and J=1-0 observations with ASTE and Mopra, we discovered a new and very young massive bipolar outflow associated with the dense dust core AGAL G337.916-00.477 (AGAL337.9-S), located at 3.2kpc. The outflow is very compact less than 1,pc and is not fully resolved in the angular resolutions of ASTE and Mopra. The maximum velocities of the outflow lobes are as high as 30-40km/s. The total source luminosity for AGAL337.9-S was estimated to be ~10^5 Lo. Compared with the other massive outflows, we concluded that the AGAL337.9-S outflow is certainly one of the most massive and youngest high-mass young stellar objects ever known in the Milky Way. The short dynamical timescale ~10^4yrs of the outflow and mid-infrared-faint spectral energy distribution of AGAL337.9-S support for its young age. We also found that another dust core AGAL G337.922-00.456 (AGAL337.9-N) at 0.2 degree north of AGAL337.9-S is also a high-mass young stellar object in an earlier evolutional stage than AGAL337.9-S, although it is less bright in infrared than AGAL337.9-S. These two dense cores are embedded at the western rim of the galactic HII region G337.90-00.50, and we discussed that the formation of AGAL337.9-S and AGAL337.9-N and their star formation were not triggered by the interaction with G337.90-00.50 and that they were formed before the birth of the exciting high-mass star(s) of G337.90-00.50.
Three dimensional (3D) distribution of the volume-density molecular fraction, defined by fmol_rho=rho_H2/(rho_HI + rho_H2), is studied in the Milky Way Galaxy. The molecular front appears at galacto-centric distance of R \sim 8 kpc, where the phase transition from atomic to molecular hydrogen occurs suddenly with fmol_rho dropping from \sim 0.8 to 0.2 within a radial interval as narrow as \sim 0.5 kpc. The front in fmol_rho is much sharper than that for surface density molecular fraction. The front also appears in the vertical direction with a full width of the high-fmol_rho disk to be \sim 100 pc. The radial and vertical fmol_rho profiles, particularly the front behaviors, are fitted by theoretical curves calculated using the observed density profile and assumed radiation field and metallicity with exponential gradients. The molecular fraction was found to be enhanced along the Perseus and some other spiral arms. The fmol_rho arms imply that the molecular clouds are produced from HI gas in the spiral arms and are dissociated in the interarm region. We also show that there is a threshold HI density over which the gas is transformed into molecules.
The CH$^+$ ion is a key species in the initial steps of interstellar carbon chemistry. Its formation in diverse environments where it is observed is not well understood, however, because the main production pathway is so endothermic (4280 K) that it is unlikely to proceed at the typical temperatures of molecular clouds. We investigate the formation of this highly reactive molecule with the first velocity-resolved spectral mapping of the CH$^+$ $J=1-0, 2-1$ rotational transitions, three sets of CH $\Lambda$-doubled triplet lines, $^{12}$C$^+$ and $^{13}$C$^+$ $^2P_{3/2} - ^2P_{1/2}$, and CH$_3$OH 835 GHz E-symmetry Q branch transitions, obtained with Herschel/HIFI over a $\approx$12 arcmin$^2$ region centered on the Orion BN/KL source. We present the spatial morphologies and kinematics, cloud boundary conditions, excitation temperatures, column densities, and $^{12}$C$^+$ optical depths. Emission from all of C$^+$, CH$^+$, and CH is indicated to arise in the diluted gas, outside of the explosive, dense BN/KL outflow. Our models show that UV-irradiation provides favorable conditions for steady-state production of CH$^+$ in this environment. Surprisingly, no spatial or kinematic correspondences of the observed species are found with H$_2$ S(1) emission tracing shocked gas in the outflow. We propose that C$^+$ is being consumed by rapid production of CO to explain the lack of both C$^+$ and CH$^+$ in the outflow. Hence, in star-forming environments containing sources of shocks and strong UV radiation, a description of the conditions leading to CH$^+$ formation and excitation is incomplete without including the important --- possibly dominant --- role of UV irradiation.
Recent studies adopting $\lambda_{\rm Re}$, a proxy for specific angular momentum, have highlighted how early-type galaxies (ETGs) are composed of two kinematic classes for which distinct formation mechanisms can be inferred. With upcoming surveys expected to obtain $\lambda_{\rm Re}$ from a broad range of environments (e.g. SAMI, MaNGA), we investigate in this numerical study how the $\lambda_{\rm Re}$-$\epsilon_{\rm e}$ distribution of fast-rotating dwarf satellite galaxies reflects their evolutionary state. By combining N-body/SPH simulations of progenitor disc galaxies (stellar mass $\simeq$10$^{\rm 9}$ M$_{\odot}$), their cosmologically-motivated sub-halo infall history and a characteristic group orbit/potential, we demonstrate the evolution of a satellite ETG population driven by tidal interactions (e.g. harassment). As a general result, these satellites remain intrinsically fast-rotating oblate stellar systems since their infall as early as $z=2$; mis-identifications as slow rotators often arise due to a bar/spiral lifecycle which plays an integral role in their evolution. Despite the idealistic nature of its construction, our mock $\lambda_{\rm Re}$-$\epsilon_{\rm e}$ distribution at $z<0.1$ reproduces its observational counterpart from the ATLAS$^{\rm 3D}$/SAURON projects. We predict therefore how the observed $\lambda_{\rm Re}$-$\epsilon_{\rm e}$ distribution of a group evolves according to these ensemble tidal interactions.
Discrepancies between reported structure function (SF) slopes and their overall flatness as compared to expectations from the damped random walk (DRW) model, generally well-describing the variability of active galactic nuclei (AGN), has triggered us to study this problem in detail. We review common AGN variability observables and identify their most common problems. Equipped with this knowledge, we study ~9000 r-band AGN light curves from Stripe 82 of the Sloan Digital Sky Survey, using SFs described by stochastic processes with the power exponential covariance matrix of the signal. We model the "sub-ensemble" SFs in the redshift-absolute magnitude bins with the full SF equation (including the turnover and the noise part) and a single power-law (SPL; in the "red noise regime" after subtracting the noise term). The distribution of full-equation SF (SPL) slopes peaks at gamma=0.55+/-0.08 (0.52+/-0.06) and is consistent with DRW. There is a hint of a weak correlation of gamma with the luminosity and a lack of correlation with the black hole mass. The typical decorrelation time scale in optical is tau=0.97+\-0.46 year. The SF amplitude at one year obtained from the SPL fitting is SF_0=0.22+/-0.06 mag and is overestimated because SF is already at the turnover part, hence the true value is SF_0=0.20+/-0.06 mag. The asymptotic variability is SF_\infty=0.25+/-0.06 mag. It is strongly anti-correlated with both the luminosity and Eddington ratio, and is correlated with the black hole mass. The reliability of these results is fortified with Monte Carlo simulations.
We investigate the UV continuum slope $\alpha$ of a large quasar sample from SDSS DR7. By using specific continuum windows, we build two samples at lower ($0.71<z<1.19$) and higher ($1.90<z<3.15$) redshifts, which correspond to the continuum slopes at longer (NUV) and shorter (FUV) rest wavelength ranges respectively. Overall, the average continuum slopes are $-0.36$ and $-0.51$ for $\alpha_{\rm NUV}$ and $\alpha_{\rm FUV}$ with similar dispersions $\sigma_{\alpha} \sim 0.5$. For both samples, we confirm the luminosity dependence of the continuum slope, i.e., fainter quasars have redder spectra. We further find that both $\alpha_{\rm NUV}$ and $\alpha_{\rm FUV}$ have a common upper limit ($\sim 1/3$) which is almost independent of the quasar luminosity $L_{\rm bol}$. This finding implies that the intrinsic quasar continuum (or the bluest quasar), at any luminosity, obey the standard thin disk model. We propose that the other quasars with redder $\alpha$ are caused by the reddening from the dust {\it locally}. With this assumption, we employ the dust extinction scenario to model the observed $L_{\rm bol}-\alpha$ relation. We find that, a typical value of $E(B-V)\sim0.1$ to $0.3$ mag (depending on the types of extinction curve) of the quasar {\it local} dust is enough to explain the discrepancy of $\alpha$ between the observation ($\sim-0.5$) and the standard accretion disk model prediction ($\sim 1/3$).
(Abridged) We present the first public release of high-quality data products (DR1) from Hi-GAL, the {\em Herschel} infrared Galactic Plane Survey. Hi-GAL is the keystone of a suite of continuum Galactic Plane surveys from the near-IR to the radio, and covers five wavebands at 70, 160, 250, 350 and 500 micron, encompassing the peak of the spectral energy distribution of cold dust for 8 < T < 50K. This first Hi-GAL data release covers the inner Milky Way in the longitude range 68{\deg} > l > -70{\deg} in a |b|<1{\deg} latitude strip. Photometric maps have been produced with the ROMAGAL pipeline, that optimally capitalizes on the excellent sensitivity and stability of the bolometer arrays of the {\em Herschel} PACS and SPIRE photometric cameras, to deliver images of exquisite quality and dynamical range, absolutely calibrated with {\em Planck} and {\em IRAS}, and recovering extended emission at all wavelengths and all spatial scales. The compact source catalogues have been generated with the CuTEx algorithm, specifically developed to optimize source detection and extraction in the extreme conditions of intense and spatially varying background that are found in the Galactic Plane in the thermal infrared. Hi-GAL DR1 images will be accessible via a dedicated web-based image cutout service. The DR1 Compact Source Catalogues are delivered as single-band photometric lists containing, in addition to source position, peak and integrated flux and source sizes, a variety of parameters useful to assess the quality and reliability of the extracted sources, caveats and hints to help this assessment are provided. Flux completeness limits in all bands are determined from extensive synthetic source experiments and depend on the specific line of sight along the Galactic Plane. Hi-GAL DR1 catalogues contain 123210, 308509, 280685, 160972 and 85460 compact sources in the five bands, respectively.
We investigate the 1.4 GHz radio properties of 92 nearby (z<0.05) ultra hard X-ray selected Active Galactic Nuclei (AGN) from the Swift Burst Alert Telescope (BAT) sample. Through the ultra hard X-ray selection we minimise the biases against obscured or Compton-thick AGN as well as confusion with emission derived from star formation that typically affect AGN samples selected from the UV, optical and infrared wavelengths. We find that all the objects in our sample of nearby, ultra-hard X-ray selected AGN are radio quiet; 83\% of the objects are classed as high-excitation galaxies (HEGs) and 17\% as low-excitation galaxies (LEGs). While these low-z BAT sources follow the radio--far-infrared correlation in a similar fashion to star forming galaxies, our analysis finds that there is still significant AGN contribution in the observed radio emission from these radio quiet AGN. In fact, the majority of our BAT sample occupy the same X-ray--radio fundamental plane as have been observed in other samples, which include radio loud AGN --- evidence that the observed radio emission (albeit weak) is connected to the AGN accretion mechanism, rather than star formation.
We present the ALMA Band 3 and Band 6 results of 12CO(2-1), 13$CO(2-1), H30alpha recombination line, free-free emission around 98 GHz, and the dust thermal emission around 230 GHz toward the N159 East Giant Molecular Cloud (N159E) in the Large Magellanic Cloud (LMC). LMC is the nearest active high-mass star forming face-on galaxy at a distance of 50 kpc and is the best target for studing high-mass star formation. ALMA observations show that N159E is the complex of filamentary clouds with the width and length of ~1 pc and 5 pc - 10 pc, respectively. The total molecular mass is 0.92 x 10^5 Msun from the 13CO(2-1) intensity. N159E harbors the well-known Papillon Nebula, a compact high-excitation HII region. We found that a YSO associated with the Papillon Nebula has the mass of 35 Msun and is located at the intersection of three filamentary clouds. It indicates that the formation of the high-mass YSO was induced by the collision of filamentary clouds. Fukui et al. 2015 reported a similar kinematic structure toward a YSO in the N159 West region which is another YSO that has the mass larger than 35 Msun in these two regions. This suggests that the collision of filamentary clouds is a primary mechanism of high-mass star formation. We found a small molecular hole around the YSO in Papillon Nebula with sub-pc scale. It is filled by free-free and H30alpha emission. Temperature of the molecular gas around the hole reaches ~ 80 K. It indicates that this YSO has just started the distruction of parental molecular cloud.
Ultralight scalar dark matter with mass at or below the eV scale and pressure from repulsive self-interaction could form a Bose-Einstein condensate in the early Universe and maybe in galaxies as well. It has been suggested to be a possible solution to the cusp/core problem or even to explain MOND phenomenology. In this paper, I initiate a study of possible self-interactions of ultralight scalar dark matter from the particle physics point of view. To protect its mass, the scalar dark matter is identified as a pseudo Nambu-Goldstone boson (pNGB). Quite a few pNGB models with different potentials such as the QCD axion and the dilaton lead to attractive self-interactions. Yet if an axion is a remnant of a 5D gauged U(1) symmetry, its self-interactions could be repulsive provided the masses and charges of the 5D matter contributing to its potential satisfy certain constraints. Collective symmetry breaking could also lead to a repulsive self-interaction yet with too large a strength that is ruled out by Bullet Cluster constraints. I also discuss cosmological and astrophysical constraints on ultralight repulsive dark matter in terms of a parametrization motivated by particle physics considerations.
Neutron-capture elements, those with Z > 35, are the least well-understood in terms of nucleosynthesis and formation environments. The rapid neutron-capture, or r-process, elements are formed in the environments and/or remnants of massive stars, while the slow neutron-capture, or s-process, elements are primarily formed in low-mass AGB stars. These elements can provide much information about Galactic star formation and enrichment, but observational data is limited. We have assembled a sample of 68 stars in 23 open clusters that we use to probe abundance trends for six neutron-capture elements (Eu, Gd, Dy, Mo, Pr, and Nd) with cluster age and location in the disk of the Galaxy. In order to keep our analysis as homogenous as possible, we use an automated synthesis fitting program, which also enables us to measure multiple (3-10) lines for each element. We find that the pure r-process elements (Eu, Gd, and Dy) have positive trends with increasing cluster age, while the mixed r- and s- process elements (Mo, Pr, and Nd) have insignificant trends consistent with zero. Pr, Nd, Eu, Gd, and Dy have similar, slight (though mostly statistically significant) gradients of ~0.04 dex/kpc. The mixed elements also appear to have nonlinear relationships with Galactocentric radius.
Supernova (SN) explosions play an important role in the development of galactic structures. The energy and momentum imparted on the interstellar medium (ISM) in so called "supernova feedback" drives turbulence, heats the gas, enriches it with heavy elements, can lead to the formation of new stars or even suppress star formation by disrupting stellar nurseries. In the numerical simulation at the sub-galactic level, not including the energy and momentum of supernovas in the physical description of the problem can also lead to several problems that might partially be resolved by including a description of supernovas. In this thesis such an implementation is attempted for the combined numerical hydrodynamics and N-body simulation software Arepo (Springel, 2010). In a stochastic process a large amount of thermal energy is imparted on a number of neighbouring cells, mimicking the effect of a supernova explosions. We test this approach by modelling the explosion of a single supernova in a uniform density medium and comparing the evolution of the resulting supernova remnant to the theoretically-predicted behaviour. We also run a simulation with our feedback code and a fixed supernova rate derived from the Kennicutt-Schmidt relation (Kennicutt, 1998) for a duration of about 20 Myrs. We describe our method in detail in this text and discuss the properties of our implementation.
Links to: arXiv, form interface, find, astro-ph, recent, 1604, contact, help (Access key information)
We extend the analysis of the fluctuations in the velocity slices of Position-Position- Velocity (PPV) spectroscopic data from Doppler broadened lines, i.e. Velocity Channel Analysis (VCA) introduced by Lazarian & Pogosyan (2000), to study anisotropy of the underlying velocity and density turbulence statistics that arises from the presence of magnetic field. In particular, we study analytically how the measurable anisotropy of the statistics of the channel map fluctuations changes with the thickness of velocity channels. In agreement with the earlier VCA studies we find that the anisotropy of the thick channels reflects the anisotropy of the density field, while the relative contribution of density and velocity fluctuations to the thin velocity channels depends on the density spectral slope. We show that the anisotropies arising from Alfven, slow and fast modes are different, in particular, the anisotropy in PPV created by fast modes is opposite to that created by Alfven and slow modes and this can be used to separate their contributions. We successfully compare our results with the earlier numerical study of the PPV anisotropies measured with synthetic observations in Esquivel et al. 2015. We also obtain the anisotropies for the media with absorption as well as for the absorption lines. In addition, we demonstrate how the studies of anisotropy can be performed using interferometers.
We present the stellar mass ($M_{*}$), and K-corrected $K$-band absolute magnitude ($M_{K}$) Tully-Fisher relations (TFRs) for sub-samples of the 584 galaxies spatially resolved in H$\alpha$ emission by the KMOS Redshift One Spectroscopic Survey (KROSS). We model the velocity field of each of the KROSS galaxies and extract a rotation velocity, $V_{80}$ at a radius equal to the major axis of an ellipse containing 80% of the total integrated H$\alpha$ flux. The large sample size of KROSS allowed us to select 210 galaxies with well measured rotation speeds. We extract from this sample a further 56 galaxies that are rotationally supported, using the stringent criterion $V_{80}/\sigma > 3$, where $\sigma$ is the flux weighted average velocity dispersion. We find the $M_{K}$ and $M_{*}$ TFRs for this sub-sample to be $M_{K} / \rm{mag}= (-7.3 \pm 0.9) \times [(\log(V_{80}/\rm{km\ s^{-1}})-2.25]- 23.4 \pm 0.2$ , and $\log(M_{*} / M_{\odot})= (4.7 \pm 0.4) \times [(\log(V_{80}/\rm{km\ s^{-1}}) - 2.25] + 10.0 \pm 0.3$, respectively. We find an evolution of the $M_{*}$ TFR zero-point of $-0.41 \pm 0.08$ dex over the last $\sim $8 billion years. However, we measure no evolution in the $M_{K}$ TFR zero-point over the same period. We conclude that rotationally supported galaxies of a given dynamical mass had less stellar mass at $z \sim 1$ than the present day, yet emitted the same amounts of $K$-band light. The ability of KROSS to differentiate, using integral field spectroscopy with KMOS, between those galaxies that are rotationally supported and those that are not explains why our findings are at odds with previous studies without the same capabilities.
We investigate the relative significance of radiation pressure and gas pressure in the extended narrow line regions (ENLRs) of four Seyfert galaxies from the integral field Siding Spring Southern Seyfert Spectroscopic Snapshot Survey (S7). We demonstrate that there exist two distinct types of starburst-AGN mixing curves on standard emission line diagnostic diagrams which reflect the balance between gas pressure and radiation pressure in the ENLR. In two of the galaxies the ENLR is radiation pressure dominated throughout and the ionization parameter remains constant (log U ~ 0). In the other two galaxies radiation pressure is initially important, but gas pressure becomes dominant as the ionization parameter in the ENLR decreases from log U ~ 0 to -3.4 <= log U <= -3.2. Where radiation pressure is dominant, the AGN regulates the density of the interstellar medium on kpc scales and may therefore have a direct impact on star formation activity and/or the incidence of outflows in the host galaxy to scales far beyond the zone of influence of the black hole. We find that both radiation pressure dominated and gas pressure dominated ENLRs are dynamically active with evidence for outflows, indicating that radiation pressure may be an important source of AGN feedback even when it is not dominant over the entire ENLR.
Recently, Li et al. (2016) claimed to have found evidence for multiple generations of stars in the intermediate age clusters NGC 1783, NGC 1696 and NGC 411 in the Large and Small Magellanic Clouds (LMC/SMC). Here we show that these young stellar populations are present in the field regions around these clusters and are not likely associated with the clusters themselves. Using the same datasets, we find that the background subtraction method adopted by the authors does not adequately remove contaminating stars in the small number Poisson limit. Hence, we conclude that their results do not provide evidence of young generations of stars within these clusters.
The evolutionary classification of massive clumps that are candidate
progenitors of high-mass young stars and clusters relies on a variety of
independent diagnostics based on observables from the near-infrared to the
radio. A promising evolutionary indicator for massive and dense
cluster-progenitor clumps is the L/M ratio between the bolometric luminosity
and the mass of the clumps. With the aim of providing a quantitative
calibration for this indicator we used SEPIA/APEX to obtain CH3C2H(12-11)
observations, that is an excellent thermometer molecule probing densities >
10^5 cm^-3 , toward 51 dense clumps with M>1000 solar masses, and uniformly
spanning -2 < Log(L/M) < 2.3.
We identify three distinct ranges of L/M that can be associated to three
distinct phases of star formation in massive clumps. For L/M <1 no clump is
detected in CH3C2H , suggesting an inner envelope temperature below 30K. For 1<
L/M < 10 we detect 58% of the clumps, with a temperature between 30 and 35 K
independently from the exact value of L/M; such clumps are building up
luminosity due to the formation of stars, but no star is yet able to
significantly heat the inner clump regions. For L/M> 10 we detect all the
clumps, with a gas temperature rising with Log(L/M), marking the appearance of
a qualitatively different heating source within the clumps; such values are
found towards clumps with UCHII counterparts, suggesting that the quantitative
difference in T - L/M behaviour above L/M >10 is due to the first appearance of
ZAMS stars in the clumps.
We present wide-field $g$ and $i$ band stellar photometry of the Sextans dwarf spheroidal galaxy and its surrounding area out to four times its half-light radius ($r_h=695\,$pc), based on images obtained with the Dark Energy Camera at the 4-m Blanco telescope at CTIO. We find clear evidence of stellar substructure associated with the galaxy, extending to a distance of $82\arcmin$ (2\,kpc) from its centre. We perform a statistical analysis of the over-densities and find three distinct features, as well as an extended halo-like structure, to be significant at the $99.7\%$ confidence level or higher. Unlike the extremely elongated and extended substructures surrounding the Hercules dwarf spheroidal galaxy, the over-densities seen around Sextans are distributed evenly about its centre, and do not appear to form noticeable tidal tails. Fitting a King model to the radial distribution of Sextans stars yields a tidal radius $r_t =83.2\arcmin\pm7.1\arcmin$ (2.08$\pm$0.18\,kpc), which implies the majority of detected substructure is gravitationally bound to the galaxy. This finding suggests that Sextans is not undergoing significant tidal disruption from the Milky Way, supporting the scenario in which the orbit of Sextans has a low eccentricity.
This paper reports results of the third-year campaign of monitoring super-Eddington accreting massive black holes (SEAMBHs) in active galactic nuclei (AGNs) between 2014-2015. Ten new targets were selected from quasar sample of Sloan Digital Sky Survey (SDSS), which are generally more luminous than the SEAMBH candidates in last two years. H$\beta$ lags ($\tau_{_{\rm H\beta}}$) in five of the 10 quasars have been successfully measured in this monitoring season. We find that the lags are generally shorter, by large factors, than those of objects with same optical luminosity, in light of the well-known $R_{_{\rm H\beta}}-L_{5100}$ relation. The five quasars have dimensionless accretion rates of $\dot{\mathscr{M}}=10-10^3$. Combining measurements of the previous SEAMBHs, we find that the reduction of H$\beta$ lags tightly depends on accretion rates, $\tau_{_{\rm H\beta}}/\tau_{_{R-L}}\propto\dot{\mathscr{M}}^{-0.42}$, where $\tau_{_{R-L}}$ is the H$\beta$ lag from the normal $R_{_{\rm H\beta}}-L_{5100}$ relation. Fitting 63 mapped AGNs, we present a new scaling relation for the broad-line region: $R_{_{\rm H\beta}}=\alpha_1\ell_{44}^{\beta_1}\,\min\left[1,\left(\dot{\mathscr{M}}/\dot{\mathscr{M}}_c\right)^{-\gamma_1}\right]$, where $\ell_{44}=L_{5100}/10^{44}\,\rm erg~s^{-1}$ is 5100 \AA\ continuum luminosity, and coefficients of $\alpha_1=(29.6_{-2.8}^{+2.7})$ lt-d, $\beta_1=0.56_{-0.03}^{+0.03}$, $\gamma_1=0.52_{-0.16}^{+0.33}$ and $\dot{\mathscr{M}}_c=11.19_{-6.22}^{+2.29}$. This relation is applicable to AGNs over a wide range of accretion rates, from $10^{-3}$ to $10^3$. Implications of this new relation are briefly discussed.
Understanding how galaxies evolved from the early Universe through cosmic
time is a fundamental part of modern astrophysics. In order to study this
evolution it is important to sample the galaxies at various times in a
consistent way through time. In regular luminosity selected samples, our
analyses are biased towards the brightest galaxies at all times (as these are
easier to observe and identify). A complementary method relies on the
absorption imprint from neutral gas in galaxies, the so-called damped Ly-alpha
absorbers (DLAs) seen towards distant bright objects. This thesis seeks to
understand how the absorption selected galaxies relate to the emission selected
galaxies by identifying the faint glow from the absorbing galaxies at redshift
z~2.
In the last Chapter, a study of the more evolved, massive galaxies is
presented. These galaxies are observed to be a factor of 2 to 6 times smaller
than local galaxies of similar masses. A new spectroscopically selected sample
is presented and the increased precision of the redshifts allows a more
detailed measurement of the scatter in the mass-size relation. The size
evolution of massive, quiescent galaxies is modelled by a `dilution' scenario,
in which progressively larger galaxies at later times are added to the
population of denser galaxies, causing an increase of the mean size of the
population. This model describes the evolution of both sizes and number
densities very well, however, the scatter in the model increases with time,
contrary to the data. It is thus concluded that a combination of `dilution' and
individual growth, e.g., through mergers, is needed.
For brevity, the individual chapters based on published peer-reviewed
articles are omitted and a link to the given article is given instead.
Principal supervisor: Johan P. U. Fynbo
We present the star formation history of the extremely metal-poor dwarf galaxy DDO~68, based on our $V-$ and $I-$ band photometry with the Advanced Camera for Surveys on board of the Hubble Space Telescope. With a metallicity of only $12+\log(O/H)=7.15$ and an isolated location in the periphery of the nearby Lynx-Cancer void, DDO~68 is one of the most metal poor galaxies known. It has been argued in the past that DDO~68 is a young system that started forming stars only $\sim 0.15$~Gyr ago. Our data provide a deep and uncontaminated optical color-magnitude diagram that now allows us to disprove this hypothesis, since we find a population of at least $\sim 1$~Gyr old stars. The star formation activity has been fairly continuous over all the look-back time. The current rate is quite low, and the highest activity occurred between 10 and 100 Myr ago. The average star formation rate over the whole Hubble time is \mbox{$\simeq 0.01$ M$_{\odot}$ yr$^{-1}$}, providing a total mass of formed stars of \mbox{$\simeq 1.3 \times 10^8$ M$_{\odot}$}. Our photometry allows us to infer the distance to this galaxy: based on the tip of the red giant branch we estimate $D = 12.08 \pm 0.67$~Mpc, or $(m-M)_0 = 30.41 \pm 0.12$~mag, while to let our synthetic color-magnitude diagram reproduce the observed ones we need a slightly higher distance, $D=12.65$~Mpc, or $(m-M)_0 = 30.51$, still inside the errors of the previous determination. DDO~68 shows a very interesting and complex history, with its quite disturbed shape and a long tail probably due to tidal interactions.
We investigate the impact of f(R) modified gravity on the internal properties of Milky Way sized dark matter halos in a set of cosmological zoom simulations of seven halos from the Aquarius suite, carried out with our code MG-GADGET in the Hu & Sawicki f(R) model. Also, we calculate the fifth forces in ideal NFW-halos as well as in our cosmological simulations and compare them against analytic model predictions for the fifth force inside spherical objects. We find that these theoretical predictions match the forces in the ideal halos very well, whereas their applicability is somewhat limited for realistic cosmological halos. Our simulations show that f(R) gravity significantly affects the dark matter density profile of Milky Way sized objects as well as their circular velocities. In unscreened regions, the velocity dispersions are increased by up to 40% with respect to LCDM for viable f(R) models. This difference is larger than reported in previous works. The Solar circle is fully screened in $f_{R0} = -10^{-6}$ models for Milky Way sized halos, while this location is unscreened for slightly less massive objects. Within the scope of our limited halo sample size, we do not find a clear dependence of the concentration parameter of dark matter halos on $f_{R0}$.
Multiplicity is common in field stars and among protostellar systems. Models suggest two paths of formation: turbulent fragmentation and protostellar disk fragmentation. We attempt to find whether or not the coevality frequency of multiple protostellar systems can help to better understand their formation mechanism. The coevality frequency is determined by constraining the relative evolutionary stages of the components in a multiple system. SEDs for known multiple protostars in Perseus were constructed from literature data. Herschel PACS photometric maps were used to sample the peak of the SED for systems with separations >7", a crucial aspect in determining the evolutionary stage of a protostellar system. Inclination effects and the surrounding envelope and outflows were considered to decouple source geometry from evolution. This together with the shape and derived properties from the SED was used to determine each system's coevality as accurately as possible. SED models were used to examine the frequency of non-coevality that is due to geometry. We find a non-coevality frequency of 33+/-10% from the comparison of SED shapes of resolved multiple systems. Other source parameters suggest a somewhat lower frequency of non-coevality. The frequency of apparent non-coevality that is due to random inclination angle pairings of model SEDs is 17+/-0.5%. Observations of the outflow of resolved multiple systems do not suggest significant misalignments within multiple systems. Effects of unresolved multiples on the SED shape are also investigated. We find that 1/3 of the multiple protostellar systems sampled here are non-coeval, which is more than expected from random geometric orientations. The other 2/3 are found to be coeval. Higher order multiples show a tendency to be non-coeval. The frequency of non-coevality found here is most likely due to formation and enhanced by dynamical evolution.
Host galaxy identification is a crucial step for modern supernova (SN) surveys such as the Dark Energy Survey (DES) and the Large Synoptic Survey Telescope (LSST), which will discover SNe by the thousands. Spectroscopic resources are limited, so in the absence of real-time SN spectra these surveys must rely on host galaxy spectra to obtain accurate redshifts for the Hubble diagram and to improve photometric classification of SNe. In addition, SN luminosities are known to correlate with host-galaxy properties. Therefore, reliable identification of host galaxies is essential for cosmology and SN science. We simulate SN events and their locations within their host galaxies to develop and test methods for matching SNe to their hosts. We use both real and simulated galaxy catalog data from the Advanced Camera for Surveys General Catalog and MICECATv2.0, respectively. We also incorporate "hostless" SNe residing in undetected faint hosts into our analysis, with an assumed hostless rate of 5%. Our fully automated algorithm is run on catalog data and matches SNe to their hosts with 91% accuracy. We find that including a machine learning component, run after the initial matching algorithm, improves the accuracy (purity) of the matching to 97% with a 2% cost in efficiency (true positive rate). Although the exact results are dependent on the details of the survey and the galaxy catalogs used, the method of identifying host galaxies we outline here can be applied to any transient survey.
In this work, we present a study of 207 quasars selected from the Sloan Digital Sky Survey quasar catalogs and the Herschel Stripe 82 survey. Quasars within this sample are high luminosity quasars with a mean bolometric luminosity of $10^{46.4}$ erg s$^{-1}$. The redshift range of this sample is within $z<4$, with a mean value of $1.5\pm0.78$. Because we only selected quasars that have been detected in all three Herschel-SPIRE bands, the quasar sample is complete yet highly biased. Based on the multi-wavelength photometric observation data, we conducted a spectral energy distribution (SED) fitting through UV to FIR. Parameters such as active galactic nucleus (AGN) luminosity, FIR luminosity, stellar mass, as well as many other AGN and galaxy properties are deduced from the SED fitting results. The mean star formation rate (SFR) of the sample is 419 $M_{\odot}$ yr$^{-1}$ and the mean gas mass is $\sim 10^{11.3}$ $M_{\odot}$. All these results point to an IR luminous quasar system. Comparing with star formation main sequence (MS) galaxies, at least 80 out of 207 quasars are hosted by starburst galaxies. It supports the statement that luminous AGNs are more likely to be associated with major mergers. The SFR increases with the redshift up to $z=2$. It is correlated with the AGN bolometric luminosity, where $L_{\rm FIR} \propto L_{\rm Bol}^{0.46\pm0.03}$. The AGN bolometric luminosity is also correlated with the host galaxy mass and gas mass. Yet the correlation between $L_{\rm FIR}$ and $L_{\rm Bol}$ has higher significant level, implies that the link between AGN accretion and the SFR is more primal. The $M_{\rm BH}/M_{\ast}$ ratio of our sample is 0.02, higher than the value 0.005 in the local Universe. It might indicate an evolutionary trend of the $M_{\rm BH} - M_{\ast}$ scaling relation.
We study the possible rotation of cluster galaxies, developing, testing and applying a novel algorithm which identifies rotation, if such does exits, as well as its rotational centre, its axis orientation, rotational velocity amplitude and, finally, the clockwise or counterclockwise direction of rotation on the plane of the sky. To validate our algorithms we construct realistic Monte-Carlo mock rotating clusters and confirm that our method provides robust indications of rotation. We then apply our methodology on a sample of Abell clusters with z<~0.1 with member galaxies selected from the SDSS DR10 spectroscopic database. We find that ~35% of our clusters are rotating when using a set of strict criteria, while loosening the criteria we find this fraction increasing to ~48%. We correlate our rotation indicators with the cluster dynamical state, provided either by their Bautz-Morgan type or by their X-ray isophotal shape and find for those clusters showing rotation that the significance and strength of their rotation is correlated to dynamical youth but only when limiting our analysis to within a 1.5h^{-1}_{70} Mpc radius. The lack of significant such correlations when we increase the radius of analysis to 2.5h^{-1}_{70} Mpc points towards a different mechanism being responsible for the rotation of the inner and outer cluster region. The outer cluster rotation could possibly be related to coherent motions of infalling substructures. Finally, finding rotational modes in galaxy clusters could lead to the necessity of correcting the dynamical cluster mass calculations.
Links to: arXiv, form interface, find, astro-ph, recent, 1604, contact, help (Access key information)