High critical density molecular lines like HCN(1-0) or HCO+(1-0) represent our best tool to study currently star-forming, dense molecular gas at extragalactic distances. The optical depth of these lines is a key ingredient to estimate the effective density required to excite emission. However, constraints on this quantity are even scarcer in the literature than measurements of the high density tracers themselves. Here, we combine new observations of HCN, HCO+ and HNC(1-0) and their optically thin isotopologues H13CN, H13CO+ and HN13C(1-0) to measure isotopologue line ratios. We use IRAM 30-m observations from the large program EMPIRE and new ALMA observations, which together target 6 nearby star-forming galaxies. Using spectral stacking techniques, we calculate or place strong upper limits on the HCN/H13CN, HCO+/H13CO+ and HNC/HN13C line ratios in the inner parts of these galaxies. Under simple assumptions, we use these to estimate the optical depths of HCN(1-0) and HCO+(1-0) to be \tau ~2-11 in the active, inner regions of our targets. The critical densities are consequently lowered to values between 5-20$\times 10^5$, 1-3$\times 10^5$ and 9$\times 10^4$ cm-3 for HCN, HCO+ and HNC, respectively. We study the impact of having different beam-filling factors, $\eta$, on these estimates and find that the effective critical densities decrease by a factor of $\frac{\eta_{12}}{\eta_{13}}\,\tau_{12}$. A comparison to existing work in NGC 5194 and NGC 253 shows HCN/H13CN and HCO+/H13CO+ ratios in agreement with our measurements within the uncertainties. The same is true for studies in other environments such as the Galactic Centre or nuclear regions of AGN-dominated nearby galaxies.
For the first time, we present the size evolution of a mass-complete (log(M*/Msol)>10) sample of star-forming galaxies over redshifts z=1-7, selected from the FourStar Galaxy Evolution Survey (ZFOURGE). Observed H-band sizes are measured from the Cosmic Assembly Near-Infrared Deep Extragalactic Legacy Survey (CANDELS) Hubble Space Telescope (HST)/F160W imaging. Distributions of individual galaxy masses and sizes illustrate that a clear mass-size relation exists up to z~7. At z~7, we find that the average galaxy size from the mass-size relation is more compact at a fixed mass of log(M*/Msol)=10.1, with r_1/2,maj=1.02+/-0.29 kpc, than at lower redshifts. This is consistent with our results from stacking the same CANDELS HST/F160W imaging, when we correct for galaxy position angle alignment. We find that the size evolution of star-forming galaxies is well fit by a power law of the form r_e = 7.07(1 + z)^-0.89 kpc, which is consistent with previous works for normal star-formers at 1<z<4. In order to compare our slope with those derived Lyman break galaxy studies, we correct for different IMFs and methodology and find a slope of -0.97+/-0.02, which is shallower than that reported for the evolution of Lyman break galaxies at z>4 (r_e\propto(1 +z)^-1.2+/-0.06). Therefore, we conclude the Lyman break galaxies likely represent a subset of highly star-forming galaxies that exhibit rapid size growth at z>4.
We study the dynamics of charged dust grains in turbulent molecular clouds (GMCs). Massive grains behave as aerodynamic particles in primarily neutral gas, and thus are able to produce dramatic small-scale fluctuations in the dust-to-gas ratio. Hopkins & Lee 2016 directly simulated the dynamics of neutral dust grains in supersonic MHD turbulence and showed that dust-to-gas fluctuations can exceed factor ~1000 on small scales, with important implications for star formation, stellar abundances, and dust growth. However, even in primarily neutral gas in GMCs, dust grains are negatively charged and Lorentz forces are non-negligible. Therefore, we extend our previous study by including the Lorentz forces on charged grains (in addition to drag). For small charged grains (<<0.1 micron), Lorentz forces suppress dust-to-gas ratio fluctuations, while for large grains (~micron), Lorentz forces have essentially no effect, trends that are well explained with a simple theory of dust magnetization. In some special intermediate cases, Lorentz forces can enhance dust-gas segregation. For physically expected scalings of dust charge with grain size, we find the most important effects depend on grain size with Lorentz forces/charge as a second-order correction. We show that the dynamics we consider are determined by three dimensionless numbers in the limit of weak background magnetic fields: the turbulent Mach number, a dust drag parameter (proportional to grain size) and a dust Lorentz parameter (proportional to grain charge); these allow us to generalize our simulations to a wide range of conditions.
The HI in disk galaxies frequently extends beyond the optical image, and can trace the dark matter there. I briefly highlight the history of high spatial resolution HI imaging, the contribution it made to the dark matter problem, and the current tension between several dynamical methods to break the disk-halo degeneracy. I then turn to the flaring problem, which could in principle probe the shape of the dark halo. Instead, however, a lot of attention is now devoted to understanding the role of gas accretion via galactic fountains. The current $\rm \Lambda$ cold dark matter theory has problems on galactic scales, such as the core-cusp problem, which can be addressed with HI observations of dwarf galaxies. For a similar range in rotation velocities, galaxies of type Sd have thin disks, while those of type Im are much thicker. After a few comments on modified Newtonian dynamics and on irregular galaxies, I close with statistics on the HI extent of galaxies.
The outskirts of galaxies offer extreme environments where we can test our understanding of the formation, evolution, and destruction of molecules and their relationship with star formation and galaxy evolution. We review the basic equations that are used in normal environments to estimate physical parameters like the molecular gas mass from CO line emission and dust continuum emission. Then we discuss how those estimates may be affected when applied to the outskirts, where the average gas density, metallicity, stellar radiation field, and temperature may be lower. We focus on observations of molecular gas in the outskirts of the Milky Way, extragalactic disk galaxies, early-type galaxies, groups, and clusters. The scientific results show the versatility of molecular gas, as it has been used to trace Milky Way spiral arms out to a galactocentric radius of 15 kpc, to study star formation in extended ultraviolet disk galaxies, to probe galaxy interactions in polar ring S0 galaxies, and to investigate ram pressure stripping in clusters. We highlight the physical stimuli that accelerate the formation of molecular gas, including internal processes such as spiral arm compression and external processes such as interactions.
Many massive star forming disc galaxies in the redshift range 3 to 0.5 are observed to have a clumpy morphology showing giant clumps of size $\sim$1 kpc and masses of about $10^{7}M_{\odot}$ to $10^{10} M_{\odot}$. The nature and fate of these giant clumps is still under debate. In this work we use 19 high-resolution simulations of disc galaxies from the NIHAO sample to study the formation and the evolution of clumps in the discs of high redshift galaxies. For the first time we use mock HST - CANDELS observations created with the radiative transfer code GRASIL-3D to carry out an observationally motivated comparison between the clumpy fraction in the light distribution of observations and simulations. We can well reproduce the observed fraction of clumpy galaxies and its evolution with redshift. We find that dust attenuation can suppress intrinsically bright clumps and enhance less luminous ones. In our galaxy sample we only find clumps in light (u-band) from young stars but not in stellar mass surface density maps. This means that the NIHAO sample does not show clumpy stellar discs but rather a clumpy light distribution originating from clumpy star formation events. The clumps found in the NIHAO sample match observed age/color gradients as a function of distance from the galaxy center but they show no sign of inward migration. Clumps in our simulations disperse on timescales of a about a hundred Myr and their contribution to bulge growth is negligible.
The analysis of the chemical composition of galaxies provides fundamental insights into their evolution. This holds true also in the case of the outer regions of spiral galaxies. This Chapter presents the observational data, accumulated in the past few years mostly from the analysis of HII region spectra, concerning the metallicity of the outer disks of spirals that are characterized by extended HI envelopes and low star formation rates. I present evidence from the literature that the metal radial distribution flattens at large galactocentric distances, with levels of enrichment that exceed those expected given the large gas mass fractions and the weak star formation activity. The interpretation of these results leads to speculations regarding mechanisms of metal mixing in galactic disks and the possibility that metal-enriched gas infall plays a role in determining the chemical evolution of the outskirts of spirals.
The Star Formation in Nearby Clouds (SFiNCs) project is aimed at providing detailed study of the young stellar populations and star cluster formation in nearby 22 star forming regions (SFRs) for comparison with our earlier MYStIX survey of richer, more distant clusters. As a foundation for the SFiNCs science studies, here, homogeneous data analyses of the Chandra X-ray and Spitzer mid-infrared archival SFiNCs data are described, and the resulting catalogs of over 15300 X-ray and over 1630000 mid-infrared point sources are presented. On the basis of their X-ray/infrared properties and spatial distributions, nearly 8500 point sources have been identified as probable young stellar members of the SFiNCs regions. Compared to the existing X-ray/mid-infrared publications, the SFiNCs member list increases the census of YSO members by 6-200% for individual SFRs and by 40% for the merged sample of all 22 SFiNCs SFRs.
Appealing to an analytical result from mean-field theory, we show, using a generic galaxy model, that galactic dynamo action can be suppressed by small-scale magnetic fluctuations. This is caused by the magnetic analogue of the R\"{a}dler or $\Omega\times J$ effect, where rotation-induced corrections to the mean-field turbulent transport result in what we interpret to be an effective reduction of the standard $\alpha$ effect in the presence of small-scale magnetic fields.
QSO absorption spectroscopy provides a sensitive probe of both the neutral medium and diffuse ionized gas in the distant Universe. It extends 21cm maps of gaseous structures around low-redshift galaxies both to lower gas column densities and to higher redshifts. Combining galaxy surveys with absorption-line observations of gas around galaxies enables comprehensive studies of baryon cycles in galaxy outskirts over cosmic time. This Chapter presents a review of the empirical understanding of the cosmic neutral gas reservoir from studies of damped Lya absorbers (DLAs). It describes the constraints on the star formation relation and chemical enrichment history in the outskirts of distant galaxies from DLA studies. A brief discussion of available constraints on the ionized circumgalactic gas from studies of lower column density Lya absorbers and associated ionic absorption transitions is presented at the end.
A systematic study is performed on radio properties of H$_{2}$O megamaser host Seyfert 2 galaxies, through multi-band radio continuum observations (at 11cm, 6.0 cm, 3.6 cm, 2.0 cm and 1.3 cm) with the Effelsberg 100-m radio telescope within a total time duration of four days. For comparison, a control Seyfert 2 galaxy sample without detected maser emission was also observed. Spectral indices were determined for those sources for which measurements exist at two adjacent bands assuming a power-law dependence S$_\nu \propto \nu^{-\alpha}$, where S is the flux density and $\nu$ is the frequency. Comparisons of the radio continuum properties between megamaser and non-masing Seyfert 2s show no difference in spectral indices. However, a difference in radio luminosity is statistically significant, i.e. the maser galaxies tend to have higher radio luminosities by a factor of 2 to 3 than the non-masing ones, commonly reaching values above a critical threshold of 10$^{29}$ ergs$^{-1}$Hz$^{-1}$. This result confirms an earlier conclusion by Zhang et al. (2012), but is based on superior data with respect to the time interval within which the data were obtained, with respect to the observational facility (only one telescope used), the number of frequency bands.
Deep imaging is a fundamental tool in the study of the outermost structures of galaxies. We review recent developments in ultra-deep imaging of galaxy disks and haloes, highlighting the technical advances as well as the challenges, and summarizing observational results in the context of modern theory and simulations. The deepest modern galaxy imaging comes from three main sources: (1) surveys such as the Sloan Digital Sky Survey's Stripe 82 project, (2) very long exposures on small telescopes, including by amateurs, and (3) long exposures on the largest professional telescopes. The technical challenges faced are common in all these approaches, and include the treatment of light scattered by atmosphere and telescope/instrument, correct flat fielding, and the subtraction of non-galaxy light in the images. We review scientific results on galaxy disks and haloes obtained with deep imaging, including the detection and characterization of stellar haloes, tidal features and stellar streams, disk truncations, and thick disks. The area of ultra-deep imaging is still very much unexplored territory, and future work in this area promises significant advances in our understanding of galaxy formation and evolution.
Galaxy haloes contain fundamental clues about the galaxy formation and evolution process: hierarchical cosmological models predict haloes to be ubiquitous, and to be (at least in part) the product of past merger and/or accretion events. The advent of wide-field surveys in the last two decades has revolutionized our view of our own Galaxy and its closest "sister", Andromeda, revealing copious tidal streams from past and ongoing accretion episodes, as well as doubling the number of their known faint satellites. The focus shall now be shifted to galaxy haloes beyond the Local Group: resolving individual stars over significant areas of galaxy haloes will enable estimates of their ages, metallicities and gradients. The valuable information collected for galaxies with a range of masses, morphologies and within diverse environments will ultimately test and quantitatively inform theoretical models of galaxy formation, and shed light onto the many challenges faced by simulations on galactic scales.
In a first systematic effort to characterize the intra-night optical variability (INOV) of different classes of narrow line Seyfert 1 (NLSy1) galaxies, we have carried out observations on a sample of radio-loud (RL) and radio-quiet (RQ) NLSy1 galaxies. The RL-NLSy1 galaxies are further divided into {\gamma}-ray loud (GL) and {\gamma}-ray quiet (GQ) NLSy1 galaxies. Our sample consists of four sets, each set consisting of a RQ-NLSy1, a GQ-NLSy1 and a GL-NLSy1 galaxy, closely matched in redshift and optical luminosity. Our observations on both RQ and GQ-NLSy1 galaxies consist of a total of 19 nights, whereas the data for GL-NLSy1 galaxies (18 nights) were taken from literature published earlier by us. This enabled us to do a comparison of the duty cycle (DC) of different classes of NLSy1 galaxies. Using power-enhanced F-test, with a variability threshold of 1%, we find DCs of about 55%, 39% and 0% for GL, GQ and RQ-NLSy1 galaxies respectively. The high DC and large amplitude of INOV (24.0 +/- 13.7%) shown by GL-NLSy1 galaxies relative to the other two classes might be due to their inner aligned relativistic jets having large bulk Lorentz factors. The null DC of RQ-NLSy1 galaxies could mean the presence of low power and/or largely misaligned jets in them. However, dividing RL-NLSy1 galaxies into low and high optical polarization sources, we find that sources with large polarization show somewhat higher DCs (69%) and amplitudes (29%) compared to those with low polarization. This points to a possible link between INOV and optical polarization.
The properties and star formation processes in the far-outer disks of nearby spiral and dwarf irregular galaxies are reviewed. The origin and structure of the generally exponential profiles in stellar disks is considered to result from cosmological infall combined with a non-linear star formation law and a history of stellar migration and scattering from spirals, bars, and random collisions with interstellar clouds. In both spirals and dwarfs, the far-outer disks tend to be older, redder and thicker than the inner disks, with the overall radial profiles suggesting inside-out star formation plus stellar scattering in spirals, and outside-in star formation with a possible contribution from scattering in dwarfs. Dwarf irregulars and the far-outer parts of spirals both tend to be gas dominated, and the gas radial profile is often non-exponential although still decreasing with radius. The ratio of H-alpha to far-UV flux tends to decrease with lower surface brightness in these regions, suggesting either a change in the initial stellar mass function or the sampling of that function, or a possible loss of H-alpha photons.
We analyze the inspiral dynamics of equal-mass precessing black-hole binaries using multi-timescale techniques. The orbit-averaged post-Newtonian evolutionary equations admit two constants of motion in the equal-mass limit, namely the magnitude of the total spin $S$ and the effective spin $\xi$. This feature makes the entire dynamics qualitatively different compared to the generic unequal-mass case, where only $\xi$ is constant while the variable $S$ parametrizes the precession dynamics. For fixed individual masses and spin magnitudes, an equal-mass black-hole inspiral is uniquely characterized by the two parameters $(S,\xi)$: these two numbers completely determine the entire evolution under the effect of radiation reaction. In particular, for equal-mass binaries we find that (i) the black-hole binary spin morphology is constant throughout the inspiral, and that (ii) the precessional motion of the two black-hole spins about the total spin takes place on a longer timescale than the precession of the total spin and the orbital plane about the total angular momentum.
The existence of multiple radial solutions to the elliptic equation modeling fermionic cloud of interacting particles is proved for the limiting Planck constant and intermediate values of mass parameters. It is achieved by considering the related nonautonomous dynamical system for which the passage to the limit can be established due to the continuity of the solutions with respect to the parameter going to zero.
Pan-STARRS1 has carried out a set of distinct synoptic imaging sky surveys including the $3\pi$ Steradian Survey and Medium Deep Survey in 5 bands $grizy_{P1}$. The Pan-STARRS system andthe design of the PS1 surveys is described and an overview of the resulting image and catalog dataproducts and their basic characteristics are described. The images, reduced data products, and derived data products from the Pan-STARRS1 surveys are available to the community from the STScI MAST archive.
Abridged. The 12CO/13CO ratio in the circumstellar envelope (CSE) of asymptotic giant branch (AGB) stars has been extensively used as the tracer of the photospheric 12C/13C ratio. However, spatially-resolved ALMA observations of R Scl, a carbon rich AGB star, have shown that the 12CO/13CO ratio is not consistent over the entire CSE. Hence, it can not necessarily be used as a tracer of the 12C/13C ratio. The most likely hypothesis to explain the observed discrepancy between the 12CO/13CO and 12C/13C ratios is CO isotopologue selective photodissociation by UV radiation. Unlike the CO isotopologue ratio, the HCN isotopologue ratio is not affected by UV radiation. Therefore, HCN isotopologue ratios can be used as the tracer of the atomic C ratio in UV irradiated regions. We have performed a detailed non-LTE excitation analysis of circumstellar H12CN and H13CN line emission around R Scl, observed with ALMA and APEX, using a radiative transfer code, ALI. The spatial extent of the molecular distribution for both isotopologues is constrained based on the spatially resolved H13CN(4-3) ALMA observations. We find fractional abundances of H12CN/H2 = (5.0 +\- 2.0) x 10^{-5} and H13CN/H2 = (1.9 +\- 0.4) x 10^{-6} in the inner wind (r < (2.0 +\- 0.25) x 10^{15} cm) of R Scl. The derived circumstellar isotopologue ratio of H12CN/H13CN = 26.3 +\- 11.9 is consistent with the photospheric ratio of 12C/13C ~ 19 \pm 6. We show that the circumstellar H12CN/H13CN ratio traces the photospheric 12C/13C ratio. These results support the previously proposed explanation that CO isotopologue selective-shielding is the main factor responsible for the observed discrepancy between 12C/13C and 12CO/13CO ratios in the inner CSE of R Scl. This indicates that UV radiation impacts on the CO isotopologue ratio.
We have mapped 12CO J=3-2 and other molecular lines from the "water-fountain" bipolar pre-planetary nebula (PPN) IRAS 16342-3814 with ~0."35 resolution using ALMA. We find (i) two very high-speed knotty, jet-like molecular outflows, (ii) a central high-density (> few x 10^6 cm^{-3}), expanding torus of diameter 1300 AU, and (iii) the circumstellar envelope of the progenitor AGB, generated by a sudden, very large increase in the mass-loss rate to >3.5 x 10^{-4} Msun/yr in the past ~455 yr. Strong continuum emission at 0.89 mm from a central source (690 mJy), if due to thermally-emitting dust, implies a substantial mass (0.017 Msun) of very large (~mm-sized) grains. The measured expansion ages of the above structural components imply that the torus (age~160 yr) and the younger high-velocity outflow (age~110 yr) were formed soon after the sharp increase in the AGB mass-loss rate. Assuming a binary model for the jets in IRAS 16342, the high momentum rate for the dominant jet-outflow in IRAS 16342 implies a high minimum accretion rate, ruling out standard Bondi-Hoyle-Lyttleton wind accretion and wind Roche lobe overflow (RLOF) models with white-dwarf or main-sequence companions. Most likely, enhanced RLOF from the primary or accretion modes operating within common envelope evolution are needed.
Links to: arXiv, form interface, find, astro-ph, recent, 1612, contact, help (Access key information)
We present a suite of 3D multi-physics MHD simulations following star formation in isolated turbulent molecular gas disks ranging from 5 to 500 parsecs in radius. These simulations are designed to survey the range of surface densities between those typical of Milky Way GMCs ($\sim 10^2 M_\odot\,pc^{-2}$) and extreme ULIRG environments ($\sim 10^4 M_\odot\,pc^{-2}$) so as to map out the scaling of star formation efficiency (SFE) between these two regimes. The simulations include prescriptions for supernova, stellar wind, and radiative feedback, which we find to be essential in determining both the instantaneous ($\epsilon_{ff}$) and integrated ($\epsilon_{int}$) star formation efficiencies. In all simulations, the gas disks form stars until a critical stellar mass has been reached and the remaining gas is blown out by stellar feedback. We find that surface density is a good predictor of $\epsilon_{int}$, as suggested by analytic force balance arguments from previous works. Furthermore, SFE eventually saturates to $\sim1$ at high surface density, with very good agreement across different spatial scales. We also find a roughly proportional relationship between $\epsilon_{ff}$ and $\epsilon_{int}$. These results have implications for star formation in galactic disks, the nature and fate of nuclear starbursts, and the formation of bound star clusters. The scaling of $\epsilon_{ff}$ also contradicts star formation models in which $\epsilon_{ff} \sim 1\%$ universally, including popular subgrid models for galaxy simulations.
We use N-body chemo-dynamic simulations to study the coupling between morphology, kinematics and metallicity of the bar/bulge region of our Galaxy. We make qualitative comparisons of our results with available observations and find very good agreement. We conclude that this region is complex, since it comprises several stellar components with different properties -- i.e. a boxy/peanut bulge, thin and thick disc components, and, to lesser extents, a disky pseudobulge, a stellar halo and a small classical bulge -- all cohabiting in dynamical equilibrium. Our models show strong links between kinematics and metallicity, or morphology and metallicity, as already suggested by a number of recent observations. We discuss and explain these links.
We aim to explore the relationship between globular cluster total number, $N_{\rm GC}$, and central black hole mass, $M_\bullet$, in spiral galaxies, and compare it with that recently reported for ellipticals. We present results for the Sbc galaxy NGC 4258, from Canada France Hawaii Telescope data. Thanks to water masers with Keplerian rotation in a circumnuclear disk, NGC 4258 has the most precisely measured extragalactic distance and supermassive black hole mass to date. The globular cluster (GC) candidate selection is based on the ($u^*\ -\ i^\prime$) vs. ($i^\prime\ -\ K_s$) diagram, which is a superb tool to distinguish GCs from foreground stars, background galaxies, and young stellar clusters, and hence can provide the best number counts of GCs from photometry alone, virtually free of contamination, even if the galaxy is not completely edge-on. The mean optical and optical-near infrared colors of the clusters are consistent with those of the Milky Way and M 31, after extinction is taken into account. We directly identify 39 GC candidates; after completeness correction, GC luminosity function extrapolation and correction for spatial coverage, we calculate a total $N_{\rm GC} = 144\pm31^{+38}_{-36}$ (random and systematic uncertainties, respectively). We have thus increased to 6 the sample of spiral galaxies with measurements of both $M_\bullet$ and $N_{\rm GC}$. NGC 4258 has a specific frequency $S_{\rm N} = 0.4\pm0.1$ (random uncertainty), and is consistent within 2$\sigma$ with the $N_{\rm GC}$ vs. $M_\bullet$ correlation followed by elliptical galaxies. The Milky Way continues to be the only spiral that deviates significantly from the relation.
Recent observational evidence for initial mass function (IMF) variations in massive quiescent galaxies at $z = 0$ challenges the long-established paradigm of a universal IMF. While a few theoretical models relate the IMF to birth cloud conditions, the physical driver underlying these putative IMF variations is still largely unclear. Here we use post-processing analysis of the Illustris cosmological hydrodynamical simulation to investigate possible physical origins of IMF variability with galactic properties. We do so by tagging stellar particles in the simulation (each representing a stellar population of $\approx10^{6}~\mathrm{M}_{\odot}$) with individual IMFs that depend on various physical conditions, such as velocity dispersion, metallicity, or SFR, at the time and place the stars are formed. We then follow the assembly of these populations throughout cosmic time, and reconstruct the overall IMF of each $z=0$ galaxy from the many distinct IMFs it is comprised of. Our main result is that applying the observed relations between IMF and galactic properties to the conditions at the star-formation sites does not result in strong enough IMF variations between $z = 0$ galaxies. Steeper physical IMF relations are required for reproducing the observed IMF trends, and some stellar populations must form with more extreme IMFs than those observed. The origin of this result is the hierarchical nature of massive galaxy assembly, and it has implications for the reliability of the strong observed trends, for the ability of cosmological simulations to capture certain physical conditions in galaxies, and for theories of star-formation aiming to explain the physical origin of a variable IMF.
Previous studies have shown that star formation depends on the driving of molecular cloud turbulence, and differences in the driving can produce an order of magnitude difference in the star formation rate. The turbulent driving is characterised by the parameter $\zeta$, with $\zeta=0$ for compressive, curl-free driving (e.g. accretion or supernova explosions), and $\zeta=1$ for solenoidal, divergence-free driving (e.g. Galactic shear). Here we develop a new method to measure $\zeta$ from observations of synchrotron emission from molecular clouds. We calculate statistics of mock synchrotron intensity images produced from magnetohydrodynamic simulations of molecular clouds, in which the driving was controlled to produce different values of $\zeta$. We find that the mean and standard deviation of the log-normalised synchrotron intensity are sensitive to $\zeta$, for values of $\zeta$ between $0$ (curl-free driving) and $0.5$ (naturally-mixed driving). We quantify the dependence of zeta on the direction of the magnetic field relative to the line of sight. We provide best-fit formulae for $\zeta$ in terms of the log-normalised mean and standard deviation of synchrotron intensity, with which $\zeta$ can be determined for molecular clouds that have similar Alfv\'enic Mach number to our simulations. These formulae are independent of the sonic Mach number. Signal-to-noise ratios larger than $5$, and angular resolutions smaller than $5\%$ of the cloud diameter, are required to apply these formulae. Although there are no firm detections of synchrotron emission from molecular clouds, by combining Green Bank Telescope and Very Large Array observations it should be possible to detect synchrotron emission from molecular clouds, thereby constraining the value of $\zeta$.
Recently reported variations in the typical physical properties of Galactic and extra-Galactic molecular clouds (MCs), and in their ability to form stars have been attributed to local variations in the magnitude of interstellar pressure. Inferences from these surveys have called into question two long-standing beliefs that the MCs : 1 are Virialised entities and (2) have approximately constant surface density i.e., the validity of the Larson's third law. In this work we invoke the framework of cloud-formation via collisions between warm gas flows. Post-collision clouds forming in these realisations cool rapidly and evolve primarily via the interplay between the Non-linear Thin Shell Instability (NTSI), and the self-gravity. Over the course of these simulations we traced the temporal evolution of the surface density of the assembled clouds, the fraction of dense gas, the distribution of gas column density (NPDF), and the Virial nature of the assembled clouds. We conclude, these physical properties of MCs not only exhibit temporal variation, but their respective peak-magnitude also increases in proportion with the magnitude of external pressure, $P_{ext}$. The velocity dispersion in assembled clouds appears to follow the power-law, $\sigma_{gas}\propto P_{ext}^{0.23}$. Also, the power-law tail at higher densities becomes shallower with increasing magnitude of external pressure, for magnitudes, $P_{ext}/k_{B}\lesssim 10^{7}$ K cm$^{-3}$, at higher magnitudes such as those typically found in the Galactic CMZ ($P_{ext}/k_{B} > 10^{7}$ K cm$^{-3}$), the power-law shows significant steepening. Thus while our results are broadly consistent with inferences from various recent observational surveys, it appears, MCs hardly exhibit a unique set of properties, but rather a wide variety, that can be reconciled with a range of magnitudes of pressure between 10$^{4}$ K cm$^{-3}$ - 10$^{8}$ K cm$^{-3}$.
We present new multi-epoch Very Long Baseline Array (VLBA) observations of a set of TeV blazars drawn from our VLBA program to monitor all TeV-detected high-frequency peaked BL Lac objects (HBLs) at parsec scales. Most of these sources are faint in the radio, so they have not been well observed with VLBI by other surveys. Our previous measurements of apparent jet speeds in TeV HBLs showed apparent jet speeds that were subluminal or barely superluminal, suggesting jets with velocity structures at the parsec-scale. Here we present apparent jet speed measurements for eight new TeV HBLs, which for the first time show a superluminal tail to the apparent speed distribution for the TeV HBLs.
We present an analysis of the diffuse ultraviolet (UV) emission near the Taurus Molecular Cloud based on observations made by the Galaxy Evolution Explorer (GALEX). We used a Monte Carlo dust scattering model to show that about half of the scattered flux originates in the molecular cloud with 25% arising in the foreground and 25% behind the cloud. The best-fit albedo of the dust grains is 0.3 but the geometry is such that we could not constrain the phase function asymmetry factor (g).
We study trends in the slope of the total mass profiles and dark matter fractions within the central half-light radius of 258 early-type galaxies, using data from the volume-limited ATLAS$^{\mathrm{3D}}$ survey. We use three distinct sets of dynamical models, which vary in their assumptions and also allow for spatial variations in the stellar mass-to-light ratio, to test the robustness of our results. We confirm that the slopes of the total mass profiles are approximately isothermal, and investigate how the total-mass slope depends on various galactic properties. The most statistically-significant correlations we find are a function of either surface density, \(\Sigma_e\), or velocity dispersion, \(\sigma_e\). However there is evidence for a break in the latter relation, with a nearly universal logarithmic slope above \(\log_{10}[\sigma_e/(\si{km~s^{-1}})]\sim 2.1\) and a steeper trend below this value. For the 142 galaxies above that critical \(\sigma_e\) value, the total mass-density logarithmic slopes have a mean value \(\left\langle\gamma^\prime\right\rangle = -2.192 \pm 0.016\) (\(1\sigma\) error) with an observed rms scatter of only \(\sigma_{\gamma^\prime}=0.167 \pm 0.016\). Considering the observational errors, we estimate an intrinsic scatter of \(\sigma_{\gamma^\prime}^\mathrm{intr} \approx 0.15\). These values are broadly consistent with those found by strong lensing studies at similar radii and agree, within the tight errors, with values recently found at much larger radii via stellar dynamics or HI rotation curves (using significantly smaller samples than this work).
Observations of the frequencies of different rotational transitions of the methanol molecule have provided the most sensitive probe to date for changes in the proton-to-electron mass ratio, over space and time. Using methanol absorption detected in the gravitational lens system PKS B1830-211, changes in the proton-to-electron ratio over the last 7.5 billion years have been constrained to a fractional change less than 1.1e-07. Molecular absorption systems at cosmological distances present the best opportunity for constraining or measuring changes in the fundamental constants of physics over time, however, we are now at the stage where potential differences in the morphology of the absorbing systems and the background source, combined with their temporal evolution, provide the major source of uncertainty in some systems. Here we present the first milliarcsecond resolution observations of the molecular absorption system towards PKS B1830-211. We have imaged the absorption from the 12.2-GHz transition of methanol (which is redshifted to 6.45 GHz) toward the southwestern component and show that it is possibly offset from the peak of the continuum emission and partially resolved on milliarcsecond scales. Future observations of other methanol transitions with similar angular resolution offer the best prospects for reducing systematic errors in investigations of possible changes in the proton-to-electron mass ratio on cosmological scales.
The SUperluminous Supernova Host galaxIES (SUSHIES) survey aims to provide strong new constraints on the progenitors of superluminous supernovae (SLSNe) by understanding the relationship to their host galaxies. Here, we present the photometric properties of 53 H-poor and 16 H-rich SLSN host galaxies out to $z\sim4$. We model the spectral energy distributions of the hosts to derive physical properties (e.g., stellar mass and star-formation-rate distribution functions), which we compare with other established galaxy populations. At low redshift, H-poor SLSNe are preferentially found in very blue, low-mass galaxies with high average specific star-formation rates. As the redshift increases, the host population follows the general evolution of star-forming galaxies towards more luminous galaxies with higher absolute star-formation rates. After accounting for secular evolution, we find evidence for differential evolution in galaxy mass, but not in the $B$-band and the far UV luminosity ($3\sigma$ confidence). Most remarkable is the scarcity of hosts with stellar masses above $10^{10}$ $M_{\odot}$ for both classes of SLSNe at all redshifts. In the case of H-poor SLSNe, we attribute this to a metallicity cut-off at $\sim0.4$ solar metallicity, above which the production efficiency is stifled. However, we argue that, in addition to a low metallicity, a short-lived stellar population is also required as a regulating factor for the SLSN production. Although the host population of H-rich SLSNe is very diverse, and therefore they show a weaker dependence on environmental properties, the lack of massive hosts suggests a stifled production efficiency above $\sim0.8$ solar metallicity. The large dispersion of the SLSN-IIn host properties is in stark contrast to those of GRB, ordinary core-collapse SN, and SLSN-I host galaxies. We propose that multiple progenitor channels give rise to SLSNe-IIn.
We present our detailed spectroscopic analysis of the chemical composition of four red giant stars in the halo globular cluster NGC 6426. We obtained high-resolution spectra using the Magellan2/MIKE spectrograph, from which we derived equivalent widths and subsequently computed abundances of 24 species of 22 chemical elements. For the purpose of measuring equivalent widths, we developed a new semi-automated tool, called EWCODE. We report a mean Fe content of [Fe/H] = -2.34$\pm$0.05 dex (stat.) in accordance with previous studies. At a mean $\alpha$-abundance of [(Mg,Si,Ca)/3 Fe] = 0.39$\pm$0.03 dex, NGC 6426 falls on the trend drawn by the Milky Way halo and other globular clusters at comparably low metallicities. The distribution of the lighter $\alpha$-elements as well as the enhanced ratio [Zn/Fe] = 0.39 dex could originate from hypernova enrichment of the pre-cluster medium. We find tentative evidence for a spread in the elements Mg, Si, and Zn, indicating an enrichment scenario, where ejecta of evolved massive stars of a slightly older population polluted a newly born younger one. The heavy element abundances in this cluster fit well into the picture of metal-poor globular clusters, which in that respect appear to be remarkably homogeneous. The pattern of the neutron-capture elements heavier than Zn point towards an enrichment history governed by the r-process with only little -if any- sign of s-process contributions. This finding is supported by the striking similarity of our program stars to the metal-poor field star HD 108317.
We consider energy budgets and radiative history of 8 fading AGN, identified from mismatch between the ionizion of very extended (>10 kpc) ionized clouds and the luminosity of the nucleus viewed directly. All show significant fading on ~50,000-year timescales. We explore the use of minimum ionizing luminosity Q derived from photoionization balance in the brightest pixels in H-alpha at each projected radius. Tests using PG QSOs, and one target with detailed photoionization modeling, suggest that we can derive useful histories of individual AGN; the minimum ionizing luminosity is always an underestimate and subject to fine structure in the ionized material. These tests suggest that the underestimation from the upper envelope of Q values is roughly constant for a given object. These AGN show rapid drops and standstills; the common feature is a rapid drop in the last 20,000 years before our view of the nucleus. E-folding timescales are mostly thousands of years, with a few episodes as short as 400. In the limit of largely obscured AGN, we find additional evidence for fading, comparing lower limits from recombination balance and the maximum luminosities derived from from infrared fluxes. We compare these long-term light curves to simulations of AGN accretion; the strongest variations on these timespans are seen in models with strong and local feedback. Gemini integral-field optical spectroscopy shows a very limited role for outflows in these structures. While rings and loops of emission are common, their kinematic structure shows some to be in regular rotation. UGC 7342 exhibits local signatures of outflows <300 km/s, largely associated with very diffuse emission. Only in the Teacup AGN do we see outflow signatures of order 1000 km/s. Clouds around these fading AGN consist largely of tidal debris being externally illuminated but not displaced by AGN outflows. (Abridged)
We present measurements of star formation rates (SFRs) for dwarf galaxies (M*~10^8.5 Msun ) at z~1 using near-infrared slitless spectroscopy from the Hubble Space Telescope (HST) by targetting and measuring the luminosity of the H{\alpha} emission line. Our sample is derived from the Redshift One LDSS3 Emission Line Survey (ROLES), which used [O II]{\lambda} 3727 as a tracer of star formation to target very low stellar masses down to very low SFRs (~0.1 Msun yr^-1) at this epoch. Dust corrections are estimated using SED-fits and we find, by comparison with other studies using Balmer decrement dust corrections, that we require a smaller ratio between the gas phase and stellar extinction than the nominal Calzetti relation, in agreement with recent findings by other studies. By stacking the WFC3 spectra at the redshifts obtained from ground-based [O II] detections, we are able to push the WFC3 spectra to much lower SFRs and obtain the most complete spectroscopic measurement of the low mass end of the SSFR--mass relation to date. We measure a flatter low mass power law slope (-0.47 +/- 0.04) than found by other (shallower) H{\alpha}-selected samples (\approx -1), although still somewhat steeper than that predicted by the EAGLE simulation (-0.14 +/- 0.05), hinting at possible missing physics not modelled by EAGLE or remaining incompleteness for our H{\alpha} data.
Luminous high-redshift quasars can be used to probe of the intergalactic medium (IGM) in the early universe because their UV light is absorbed by the neutral hydrogen along the line of sight. They help us to measure the neutral hydrogen fraction of the high-z universe, shedding light on the end of reionization epoch. In this paper, we present a discovery of a new quasar (PSO J006.1240+39.2219) at redshift z = 6.61 +- 0.02 from Panoramic Survey Telescope & Rapid Response System 1. Including this quasar, there are nine quasars above z > 6.5 up to date. The estimated continuum brightness is M1450= 25.96 +- 0.08. PSO J006.1240+39.2219 has a strong Ly alpha emission compared with typical low-redshift quasars, but the measured near-zone region size is RNZ = 3.2 +- 1.1 proper megaparsecs, which is consistent with other quasars at z~6.
We examine the origin of the mass discrepancy-radial acceleration relation (MDAR) of disc galaxies. This is a tight empirical correlation between the disc centripetal acceleration and that expected from the baryonic component. The MDAR holds at all radii probed by disc kinematic tracers, regardless of galaxy mass or surface brightness. The relation has two characteristic accelerations; $a_0$, above which all galaxies are baryon-dominated; and $a_{\rm min}$, an effective minimum acceleration probed by disc tracers. We use a simple model to show that these arise naturally in $\Lambda$CDM. This is because: (i) disc galaxies in $\Lambda$CDM form at the centre of dark matter haloes spanning a relatively narrow range of virial mass; (ii) cold dark matter halo acceleration profiles are self-similar and have a broad maximum at the centre, reaching values bracketed precisely by $a_{\rm min}$ and $a_0$ in that mass range; and (iii) halo mass and galaxy size scale relatively tightly with the baryonic mass of a galaxy in any successful $\Lambda$CDM galaxy formation model. The MDAR relation is thus just a reflection of the self-similar nature of CDM haloes and of the physical scales introduced by the galaxy formation process.
We present high spatial resolution mid-infrared images of the nebula around the late-type carbon-rich Wolf-Rayet (WC)-OB binary system WR~112 taken by the recently upgraded VLT spectrometer and imager for the mid-infrared (VISIR) with the PAH1, NeII\_2, and Q3 filters. The observations reveal a morphology resembling a series of arc-like filaments and broken shells. Dust temperatures and masses are derived for each of the identified filamentary structures, which exhibit temperatures ranging from $179_{-6}^{+8}$ K at the exterior W2 filament to $355_{-25}^{+37}$ K in the central 3". The total dust mass summed over the features is $2.6\pm0.4\times10^{-5}$ $\mathrm{M}_\odot$. A multi-epoch analysis of mid-IR photometry of WR~112 over the past $\sim20$ yr reveals no significant variability in the observed dust temperature and mass. The morphology of the mid-IR dust emission from WR~112 also exhibits no significant expansion from imaging data taken in 2001, 2007, and 2016, which disputes the current interpretation of the nebula as a high expansion velocity ($\sim1200$ km s$^{-1}$) "pinwheel"-shaped outflow driven by the central WC-OB colliding-wind binary. An upper limit of $\lesssim120$ km s$^{-1}$ is derived for the expansion velocity assuming a distance of $4.15$ kpc. The upper limit on the average total mass-loss rate from the central 3" of WR~112 is estimated to be $\lesssim8\times10^{-6}$ $\mathrm{M}_\odot$ yr$^{-1}$. We leave its true nature as an open question, but propose that the WR~112 nebula may have formed in the outflow during a previous red or yellow supergiant phase of the central Wolf-Rayet star.
Using the classical top-hat profile, we study the non-linear growth of spherically symmetric density perturbation and structure formation in $f(T)$ gravities. In particular, three concrete models, which have been tested against the observation of large-scale evolution and linear perturbation of the universe in the cosmological scenario, are investigated in this framework, covering both minimal and nonminimal coupling cases of $f(T)$ gravities. Moreover, we consider the virialization of the overdense region in the models after they detatch from the background expanding universe and turn around to collapse. We find that there are constraints in the magnitude and occurring epoch of the initial perturbation. The existence of these constraints indicates that a perturbation that is too weak or occurs too late will not be able to stop the expanding of the overdense region. The illustration of the evolution of the perturbation shows that in $f(T)$ gravities, the initial perturbation within the constraints can eventually lead to clustering and form structure. The evolution also shows that nonminimal coupling models collapse slower than the minimal coupling one.
We investigate the role of magnetic fields in the fragmentation of
self-gravitating discs using 3D global ideal magnetohydrodynamic simulations
performed with the "phantom" smoothed particle hydrodynamics code.
For initially toroidal fields, we find two regimes. In the first, where the
cooling time is greater than five times the dynamical time, magnetic fields
reduce spiral density wave amplitudes, which in turn suppresses fragmentation.
This is the case even if the magnetic pressure is only a tenth of the thermal
pressure. The second regime occurs when the cooling time is sufficiently short
that magnetic fields cannot halt fragmentation.
We find that magnetised discs produce more massive fragments, due to both the
additional pressure exerted by the magnetic field, and the additional angular
momentum transport induced by Maxwell stresses. The fragments are confined to a
narrower range of initial semimajor axes than those in unmagnetised discs. The
orbital eccentricity and inclination distributions of unmagnetised and
magnetised disc fragments are similar. Our results suggest the fragmentation
boundary could be at cooling times a factor of two lower than predicted by
purely hydrodynamical models.
We explore the idea that the coupling between matter and spacetime is more complex than the one originally envisioned by Einstein. We propose that such coupling takes the form of a new fundamental tensor in the Einstein field equations. We then show that the introduction of this tensor can account for dark phenomenology in General Relativity, maintaining a weak field limit compatible with standard Newtonian gravitation. The same paradigm can be applied any other theory of gravitation. We show, as an example, that in the context of conformal gravity a generalised coupling is able to solve compatibility issues between the matter and the gravitational sector.
In a recent paper, Erik Verlinde has further developed the interesting possibility that spacetime and gravity may emerge from the entanglement structure of an underlying microscopic theory. In this picture dark matter arises as the response to the standard model of particle physics from the delocalized low energy degrees of freedom that build up the dark energy component of the Universe. Physics is then regulated by a characteristic acceleration scale $a_0$, identified in this model with the dark energy de Sitter radius by $a_0=cH_0\approx 5.4\times 10^{-10}\,\textrm{m/s}^2$ (using {\it Planck} data). For a point particle, or outside an extended spherically symmetric massive object, Milgrom's empirical fitting formula is recovered. However, Verlinde's theory critically departs from MOND when considering the inner structure of galaxies. For illustration, we use the the eight classical dwarf spheroidal satellites of the Milky Way. These objects are perfect testbeds for the model given their spherical symmetry, measured kinematics, and large identified missing mass. As a consistency check we show that, for reasonable stellar mass-to-light ratios, Verlinde's theory can fit the velocity dispersion profile in dwarf spheroidals with no further need of an extra dark particle component. Finally we compare our results with the recent phenomenological interpolating MOND function of McGaugh {\it et al}, and find a departure of up to 50 percent in the innermost region of these galaxies.
We derive constraints on feedback by active galactic nuclei (AGN) by setting limits on their thermal Sunyaev-Zel'dovich (SZ) imprint on the cosmic microwave background (CMB). The amplitude of any SZ signature is small and degenerate with the poorly known sub-mm spectral energy distribution of the AGN host galaxy and other unresolved dusty sources along the line of sight. Here we break this degeneracy by combining microwave and sub-mm data from Planck with all-sky far-infrared maps from the AKARI satellite. We first test our measurement pipeline using the Sloan Digital Sky Survey (SDSS) redMaPPer catalogue of galaxy clusters, finding a highly significant detection ($>$$20\sigma$) of the SZ effect together with correlated dust emission. We then constrain the SZ signal associated with spectroscopically confirmed quasi-stellar objects (QSOs) from SDSS data release 7 (DR7) and the Baryon Oscillation Spectroscopic Survey (BOSS) DR12. We obtain a low-significance ($1.6\sigma$) hint of an SZ signal, pointing towards a mean thermal energy of $\simeq 5 \times 10^{60}$ erg, lower than reported in some previous studies. A comparison of our results with high-resolution hydrodynamical simulations including AGN feedback suggests QSO host masses of $M_{200c} \sim 4 \times 10^{12}~h^{-1}M_\odot$, but with a large uncertainty. Our analysis provides no conclusive evidence for an SZ signal specifically associated with AGN feedback.
Links to: arXiv, form interface, find, astro-ph, recent, 1612, contact, help (Access key information)
We investigated the effect of photoionising feedback inside turbulent star-forming clouds, comparing the resultant star formation in both idealised profiles and more realistic cloud structures drawn from a global galaxy simulation. We performed a series of numerical simulations which compared the effect of star formation alone, photoionisation and photoionisation plus supernovae feedback. In the idealised cloud, photoionisation suppresses gas fragmentation at early times, resulting in the formation of more massive stars and an increase in the star formation efficiency. At later times, the dispersal of the dense gas causes the radiative feedback effect to switch from positive to negative as the star formation efficiency drops. In the cloud extracted from the global simulation, the initial cloud is heavily fragmented prior to the stellar feedback beginning and is largely structurally unaffected by the late injection of radiation energy. The result is a suppression of the star formation. We conclude that the efficiency of feedback is heavily dependent on the gas structure, with negative feedback dominating when the density is high.
The initial mass function of the first, Population III (Pop III), stars plays a vital role in shaping galaxy formation and evolution in the early Universe. One key remaining issue is the final fate of secondary protostars formed in the accretion disc, specifically whether they merge or survive. We perform a suite of hydrodynamic simulations of the complex interplay between fragmentation, protostellar accretion, and merging inside dark matter minihaloes. Instead of the traditional sink particle method, we employ a stiff equation of state approach, so that we can more robustly ascertain the viscous transport inside the disc. The simulations show inside-out fragmentation because the gas collapses faster in the central region. Fragments migrate on the viscous timescale, over which angular momentum is lost, enabling them to move towards the disc centre, where merging with the primary protostar can occur. This process depends on the fragmentation scale, such that there is a maximum scale of $(1 - 5) \times 10^4$ au, inside which fragments can migrate to the primary protostar. Viscous transport is active until radiative feedback from the primary protostar destroys the accretion disc. The final mass spectrum and multiplicity thus crucially depends on the effect of viscosity in the disc. In difference to Population I discs, where viscosity is smaller, almost the entire disc is subjected to efficient viscous transport in the primordial case. An important aspect of this question is the survival probability of Pop III binary systems, possible gravitational wave sources to be probed with the Advanced LIGO detectors.
We present results from MMT/Hectochelle spectroscopy of red giant candidate stars along the line of sight to the recently-discovered Galactic satellite Crater 2. Modelling the joint distribution of stellar positions, velocities and metallicities as a mixture of Crater 2 and Galactic foreground populations, we identify 62 members of Crater 2, for which we resolve line-of-sight velocity dispersion 2.7 +/- 0.3 km/s about mean velocity of 87.5 +/- 0.4 km/s. We also resolve a metallicity dispersion 0.22 about a mean of [Fe/H]=-1.98 +/- 0.1 that is 0.28 +/- 0.14 poorer than is estimated from photometry. Despite Crater 2's relatively large size (projected halflight radius R(h)=1 kpc) and intermediate luminosity (M_V =-8), its velocity dispersion is the coldest that has been resolved for any dwarf galaxy. These properties make Crater 2 the most extreme low-density outlier in dynamical as well as structural scaling relations among the Milky Way's dwarf spheroidals. Even so, under assumptions of dynamical equilibrium and negligible contamination by unresolved binary stars, the observed velocity distribution implies a gravitationally dominant dark matter halo, with dynamical mass 4.4 x 10^6 Msun and mass-to-light ratio M/L=53 enclosed within a radius of 1 kpc, where the equivalent circular velocity is 4.3 km/s.
I review briefly some dynamical models of structures in the outer parts of disc galaxies, including models of polar rings, tidal tails and bridges. I then discuss the density distribution in the outer parts of discs. For this, I compare observations to results of a model in which the disc galaxy is in fact the remnant of a major merger, and find good agreement. This comparison includes radial profiles of the projected surface density and of stellar age, as well as time evolution of the break radius and of the inner and outer disc scale lengths. I also compare the radial projected surface density profiles of dynamically motivated mono-age populations and find that, compared to older populations, younger ones have flatter density profiles in the inner region and steeper in the outer one. The break radius, however, does not vary with stellar age, again in good agreement with observations.
Star forming galaxies have long been considered the dominant sources of the cosmic ultraviolet background radiation at early epochs. However, observing and characterizing the galaxy population with significant ionizing emission has proven to be challenging. In particular, the fraction of ionizing radiation that escapes the local environment to the intergalactic medium is poorly known. We investigate the relation between the escape fraction and galaxy luminosity. We combine the deep ultraviolet observations of Hubble Ultra Deep Field (UVUDF) with the deep Multi Unit Spectroscopic Explorer (MUSE) observations of the same field, collecting a sample of 165 faint star forming galaxies in the $3 < z < 4$ redshift range with deep rest-frame observations of the Lyman continuum. We identify one galaxy as a candidate source of detected LyC radiation. We bin the galaxies in various redshift and brightness intervals and stack their images. From stacked images we estimate the relative escape fraction upper limits as a function of the luminosity. Thanks to the depth of the sample we measure meaningful 1$\sigma$ upper limits of $f_{esc,rel} < 0.07, 0.2$ and 0.6 at $L \sim L_{\rm z=3}^{*}, 0.5L_{\rm z=3}^{*}$ and $0.1L_{\rm z=3}^{*}$, respectively. We use our estimates and theoretical predictions from the literature to study a possible dependence of the escape fraction on galaxy luminosity by modelling the ionizing background with different prescriptions of $f_{\rm esc} (M_{\rm UV})$. We show that the understanding of the luminosity dependence hinges on the ability to constrain the escape fraction down to $M_{\rm UV} \sim -18$ mag in the future.
We present a quick-start guide to BOND, a statistical method to derive oxygen and nitrogen abundances in H II regions. BOND compares a set of carefully selected strong and semistrong emission lines to a grid photoionization models. The first novelty, in comparison to other statistical methods, is that BOND relies on the [Ar III]/[Ne III] emission line ratio to break the oxygen abundance bimodality. In doing so, we can measure oxygen and nitrogen abundances without assuming any a priori relation between N/O and O/H. The second novelty is that BOND takes into account changes in the hardness of the ionizing radiation field, which can come about due to the ageing of H II regions or the stochastically sampling of the IMF. We use the emission line ratio He I/Hb, in addition to commonly used strong lines, to constrain the hardness of the ionizing radiation field. Finally, we also stress the pragmatic considerations behind our Bayesian inference.
The Universe is almost totally unexplored at low surface brightness levels. In spite of great progress in the construction of large telescopes and improvements in the sensitivity of detectors, the limiting surface brightness of imaging observations has remained static for about forty years. Recent technical advances have at last begun to erode the barriers preventing progress. In this Chapter we describe the technical challenges to low surface brightness imaging, describe some solutions, and highlight some relevant observations that have been undertaken recently with both large and small telescopes. Our main focus will be on discoveries made with the Dragonfly Telephoto Array (Dragonfly), which is a new telescope concept designed to probe the Universe down to hitherto unprecedented low surface brightness levels. We conclude by arguing that these discoveries are probably only scratching the surface of interesting phenomena that are observable when the Universe is explored at low surface brightness levels.
Central galaxies (CGs) in galaxy groups and clusters are believed to form and assemble a good portion of their stellar mass at early times, but also to accrete significant mass at late times via galactic cannibalism, that is merging with companion group or cluster galaxies that experience dynamical friction against the common host dark-matter halo. The effect of these mergers on the structure and kinematics of the CG depends not only on the properties of the accreted satellites, but also on the orbital parameters of the encounters. Here we present the results of numerical simulations aimed at estimating the distribution of merging orbital parameters of satellites cannibalized by the CGs in groups and clusters. As a consequence of dynamical friction, the satellites' orbits evolve losing energy and angular momentum, with no clear trend towards orbit circularization. The distributions of the orbital parameters of the central-satellite encounters are markedly different from the distributions found for halo-halo mergers in cosmological simulations. The orbits of satellites accreted by the CGs are on average less bound and less eccentric than those of cosmological halo-halo encounters. We provide fits to the distributions of the central-satellite merging orbital parameters that can be used to study the merger-driven evolution of the scaling relations of CGs.
Based on an extensive redshift survey for galaxy cluster Abell 2029 and Coma, we measure the luminosity functions (LFs), stellar mass functions (SMFs) for the entire cluster member galaxies. Most importantly, we measure the velocity dispersion functions (VDFs) for quiescent members. The MMT/Hectospec redshift survey for galaxies in A2029 identifies 982 spectroscopic members; for 838 members we derive the central velocity dispersion from the spectroscopy. Coma is the only other cluster surveyed as densely. The LFs, SMFs and VDFs for A2029 and Coma are essentially identical. The SMFs of the clusters are consistent with simulations. The A2029 and Coma VDFs for quiescent galaxies have a significantly steeper slope than those of field galaxies for velocity dispersion $\lesssim 100$ km s$^{-1}$. The cluster VDFs also exceed the field at velocity dispersion $\gtrsim 250$ km s$^{-1}$. The differences between cluster and field VDFs are potentially important tests of simulations and of the formation of structure in the universe.
We present a radio and optical study of the double-double radio galaxy J1328+2752 based on new low-frequency GMRT observations and SDSS data. The radio data were used to investigate the morphology and to perform a spectral index analysis. In this source we find that the inner double is misaligned by $\sim$30$^\circ$ from the axis of the outer diffuse structure. The SDSS spectrum shows that the central component has double-peaked line profiles with different emission strengths. The average velocity off-set of the two components is 235$\pm$10.5 km s$^{-1}$. The misaligned radio morphology along with the double-peaked emission lines indicate that this source is a potential candidate binary supermassive black hole. This study further supports mergers as a possible explanation for repeated jet activity in radio sources.
Revealing the properties of molecular gas in starburst galaxies is key to understand the process of star formation across cosmic history. In this paper, we present our recent high-resolution (\(\sim1''\) or 50 pc) ALMA observations of the nearby, barred starburst galaxy NGC 1808, and highlight some of the main findings: (1) the discovery of a double peak revealed in the distributions of CO (3-2) and CS (2-1) in the circumnuclear disk indicating the presence of a molecular gas torus with a radius of \(r\sim32\) pc; inside the torus, we found a compact source of dust continuum at 0.87 mm and CO, coincident with the location of the low-luminosity AGN revealed by X-ray observations; (2) the excitation of CO gas, traced by the line intensity ratio of CO (3-2) to CO (1-0) is elevated to \(\sim1\) in the star-forming disk; the ratio is between 0.5 and 1 in the 500 pc ring; (3) the ratio of HCN (1-0) to HCO\(^{+}\) (1-0) in the central 1 kpc exhibits a radial gradient: while the ratio decreases from \(\sim1.6\) in the center to \(\sim0.9\) in the 500 pc disk, it exhibits a peak of \(\sim2\) in a ring-like structure at a galactocentric radius of \(r\sim100\) pc. The detection of SiO (2-1) within the CND suggests the presence of shocks, likely generated by the intense nuclear star formation (supernova explosions and stellar winds from massive stars), as well as inflows and outflows of turbulent gas in the central 100 pc.
We investigate two-dimensional image decomposition of nearby, morphologically selected early-type galaxies (ETGs). We are motivated by recent observational evidence of significant size growth of quiescent galaxies and theoretical development advocating a two-phase formation scenario for ETGs. We find that a significant fraction of nearby ETGs show changes in isophotal shape that require multi-component models. The characteristic sizes of the inner and outer component are $\sim 3$ and $\sim 15$ kpc. The inner component lies on the mass-size relation of ETGs at $z \sim 0.25-0.75$, while the outer component tends to be more elliptical and hints at a stochastic buildup process. We find real physical differences between the single- and double-component ETGs, with the double-component galaxies being younger and more metal-rich. The fraction of double component ETGs increases with increasing $\sigma$ and decreases in denser environments. We hypothesize that double-component systems were able to accrete gas and small galaxies until later times, boosting their central densities, building up their outer parts, and lowering their typical central ages. In contrast, the oldest galaxies, perhaps due to residing in richer environments, have no remaining hints of their last accretion episode.
We investigate a sample of 2293 ICRF2 extragalactic radio-loud sources with accurate positions determined by VLBI, mostly active galactic nuclei (AGN) and quasars, which are cross-matched with optical sources in the first Gaia release (Gaia DR1). The distribution of offsets between the VLBI sources and their optical counterparts is strongly non-Gaussian, with powerful wings extending beyond 1 arcsecond. Limiting our analysis to only high-confidence difference detections, we find (and publish) a list of 188 objects with normalized variances above 12 and offsets below 1 arcsecond. Pan-STARRS stacked and monochromatic images resolve some of these sources indicating the presence of double sources, confusion sources, or pronounced extended structures. Some 89 high-quality objects, however, do not show any perturbations and appear to be star-like single sources, yet displaced by multiples of the expected error from the radio-loud AGN. We conclude that a fraction of luminous AGN (more than 4%) can be physically dislodged from the optical centers of their parent galaxies.
In this work we have a closer look at the gas content or fraction and the associated star formation rate in main sequence and starburst galaxies at z=0 and z~1-2 by applying an analytical model of galactic clumpy gas disks to samples of local spiral galaxies, ULIRGs, submillimeter (smm), and high-z starforming galaxies. The model gas and dust temperatures are determined by the heating and cooling equilibrium. Dense clouds are heated by turbulent mechanical and cosmic ray heating. The molecular abundances of individual gas clouds are determined by a detailed chemical network involving the cloud lifetime, density, and temperature. Molecular line emission is calculated with an escape probability formalism. The model calculates simultaneously the total gas mass, HI/H_2 mass, the gas velocity dispersion, IR luminosity, IR spectral energy distribution, CO spectral line energy distribution (SLED), HCN(1-0), and HCO+(1-0) emission of a galaxy given its size, integrated star formation rate, stellar mass radial profile, rotation curve, and Toomre Q parameter. The model reproduces the observed CO luminosities and SLEDs of all sample galaxies within the model uncertainties (~0.3 dex). Whereas the CO emission is robust against the variation of model parameters, the HCN and HCO+ emission is sensitive to the chemistry of the interstellar medium. The CO and HCN mass-to-light conversion factors including CO-dark H_2 are given and compared to the values found in the literature. Both, the HCN and HCO+ emission trace the dense molecular gas to a factor of ~2 for the local spiral galaxies, ULIRGs and smm-galaxies. About 80% of the molecular line emission of compact starburst galaxies originates in non-selfgravitating gas clouds. The integrated Kennicutt-Schmidt law has a slope of ~1 for the local spirals, ULIRGs, and smm-galaxies, whereas the slope is 1.7 for high-z starforming galaxies.
We present high resolution radio continuum observations with the Karl G. Jansky Very Large Array at 6, 8.5, 11.5 and 15 GHz of the double-peaked emission-line galaxy 2MASXJ12032061+1319316. The radio emission has a prominent S-shaped morphology with highly symmetric radio jets that extend over a distance of $\sim1.5^{\prime\prime}$ (1.74~kpc) on either side of the core of size $\sim0.1^{\prime\prime}$(116~pc). The radio jets have a helical structure resembling the precessing jets in the galaxy NGC~326 which has confirmed dual active galactic nuclei (AGN). The nuclear bulge velocity dispersion gives an upper limit of $(1.56\pm$0.26$)\times$10$^8$~M$_{\odot}$ for the total mass of nuclear black hole(s). We present a simple model of precessing jets in 2MASXJ1203 and find that the precession timescale is around 10$^5$ years: this matches the source lifetime estimate via spectral aging. We find that the expected super massive black hole (SMBH) separation corresponding to this timescale is 0.02 pc. We used the double peaked emission lines in 2MASXJ1203 to determine an orbital speed for a dual AGN system and the associated jet precession timescale, which turns out to be more than the Hubble time, making it unfeasible. We conclude that the S-shaped radio jets are due to jet precession caused either by a binary/dual SMBH system, a single SMBH with a tilted accretion disk or a dual AGN system where a close pass of the secondary SMBH in the past has given rise to jet precession.
Aims: This work is the first stage of a campaign to search for IMBHs, in low luminosity AGN (LLAGN) and dwarf galaxies. An additional, and equally important, aim of this pilot study is to investigate the consistency between the predictions of several BH scaling relations and the fundamental plane of black hole activity (FP-BH). Methods: We use X-ray and radio luminosity relations in accreting BHs, along with the latest scaling relations between the mass of the central black hole (MBH) and the properties of its host spheroid, to predict MBH in seven LLAGN, that were previously reported to be in the IMBH regime. Namely, we use the recently re-evaluated MBH - Msph (Msph : spheroid absolute magnitude at 3.6 microm) scaling relation for spiral galaxies, the MBH - nsph (nsph : major axis Sersic index of the bulge) relation, the MBH - PA (PA: pitch angle) relation and the FP-BH for weakly accreting BHs, to independently estimate MBH in all seven galaxies. Results: We find that all LLAGN in our list have low-mass central black holes - with log MBH/Msol~6.5 on average - but they are, not IMBHs. All 4 methods used predicted consistent BH masses, in the 1sigma range. Furthermore, we report that, in contrast to previous classification, galaxy NGC 4470 is bulge-less, and cast doubts on the AGN classification of NGC 3507. Conclusions: We find that our latest, state-of-the-art techniques for bulge magnitude & Sersic index computations and the most recent updates of the MBH - Lsph, MBH - nsph , MBH - PA relations and the FP-BH produce consistent results in the low mass regime. We establish a multiple method approach for predicting BH masses in the regime where their spheres-of- grav. influence cannot be spatially resolved. Our approach mitigates against outliers from any one relation and provides a more robust average prediction. We will use our new method to revisit more IMBH candidates in LLAGN
We analyse the spectra of the archival XMM-Newton data of the Seyfert 1 AGN Zw 229.015 in the energy range $0.3 - 10.0$ keV. When fitted with a simple power-law, the spectrum shows signatures of weak soft excess below 1.0 keV. We find that both thermal comptonisation and relativistically blurred reflection models provide the most acceptable spectral fits with plausible physical explanations to the origin of the soft excess than do multicolour disc blackbody and smeared wind absorption models. This motivated us to study the variability properties of the soft and the hard X-ray emissions from the source and the relationship between them to put further constraints on the above models. Our analysis reveals that the variation in the $3.0 - 10.0$ keV band lags that in the $0.3 - 1.0$ keV by $600^{+290}_{-280}$ s, while the lag between the $1.0 - 10.0$ keV and $0.3 - 1.0$ keV is $980^{+500}_{-500}$ s. This implies that the X-ray emissions are possibly emanating from different regions within the system. From these values, we estimate the X-ray emission region to be within $20R_g$ of the central supermassive black hole (where $R_g=GM/c^2$, $M$ is the mass of black hole, $G$ the Newton gravitational constant and $c$ the speed of light). Furthermore, we use XMM-Newton and Kepler photometric lightcurves of the source to search for possible nonlinear signature in the flux variability. We find evidence that the variability in the system may be dominated by stochasticity rather than deterministic chaos which has implications for the dynamics of the accretion system.
We describe the concept of a shallow-water setup for simulation of gas accretion onto a black hole in the mode of a quasi-spherical accretion. The bottom for the shallow-water container must have the funnel-shaped curvilinear concavo-convex shape. We calculate the configuration surface of the properly shaped bottom that simulates precisely the Newtonian or pseudo-Newtonian gravitational potentials. Like the spatial part of the Schwarzchild metric, the funnel's surface metric has a (removable) singularity at the finite distance from the funnel's center and places the certain funnel's depth which we call `gravitational length'. The gravitational length is analogous to the gravitational radius and defines the equivalent of the black hole's mass in the laboratory model. The mass equivalent corresponds to $\sim 0.367\cdot 10^{12}$ g for the funnel as deep as 5 cm. We define more precisely the inviscid shallow water equations for the arbitrary bottom curvature. We show that in general case the shallow water pressure obeys the non-barotropic equation of state. We suggest the schematic course for experiments for simulation of accretion in a thick accretion disk mode as well as the Bondi-Hoyle accretion.
Broad (~10,000 km/s), double-peaked emission-line profiles of Balmer lines emitted by active galactic nuclei (AGN) are thought to originate in the outer parts of an accretion disk surrounding a nuclear supermassive black hole (SMBH), at R~1000 gravitational radii and are most frequently observed in the nuclear spectra of low-luminosity AGN (LLAGN) and radio-galaxies. In the present paper we argue that broad double-peaked profiles are present also in the spectra of other Type 1 AGN, such as Seyfert 1 galaxies, suggesting that the inner part of the broad-line region (BLR) is also the outer part of the accretion disk. We use the Palomar spectral survey of nearby galaxies to show that the only difference between Seyfert 1 BLR line profiles and "bona fide" double peakers is that, in most cases, besides a disk component, we need an additional Gaussian component attributed to non-disk clouds. The recognition that the inner and most variable part of the BLR has a disk geometry suggests that the factor 'f ' in the expression to obtain the SMBH mass (M) in Type 1 AGN M(SMBH)=f [R(BLR) (Delta V)^2]/G is f = 1/[sin^2(i)] for the disk dominated sources. Our median inclination i=27 degrees implies f=4.5, very close to the most recent value of f=4.3+/-1.05, obtained from independent studies. We derive a relation between f and the FWHM of the broad profile that may help to reduce the uncertainties in the SMBH mass determinations of AGN.
The recent low value of Planck (2016) integrated optical depth to Thomson scattering suggests that the reionization occurred fairly suddenly, disfavoring extended reionization scenarios. This will have a significant impact on the 21cm power spectrum. Using a semi-numerical framework, we improve our model from Hassan et al. (2016) to include time-integrated ionisation and recombination effects, and find that this leads to more sudden reionisation. It also yields larger HII bubbles which leads to an order of magnitude more 21cm power on large scales, while suppressing the small scale ionization power. Local fluctuations in the neutral hydrogen density play the dominant role in boosting the 21cm power spectrum on large scales, while recombinations are subdominant. We use a Monte Carlo Markov Chain approach to constrain our model to observations of the star formation rate functions at z = 6,7,8 from Bouwens et al. (2015), the Planck (2016) optical depth measurements, and the Becker & Bolton (2013) ionising emissivity data at z~5. We then use this constrained model to perform 21cm forecasting for LOFAR, HERA, and SKA in order to determine how well such data can characterise the sources driving reionisation. We find that the 21cm power spectrum alone can somewhat constrain the halo mass dependence of ionising sources, the photon escape fraction and ionising amplitude, but combining the 21cm data with other current observations enables us to separately constrain all these parameters. Our framework illustrates how 21cm data can play a key role in understanding the sources and topology of reionisation as observations improve.
We present the first full catalog and science results for the NuSTAR serendipitous survey. The catalog incorporates data taken during the first 40 months of NuSTAR operation, which provide ~20Ms of effective exposure time over 331 fields, with an areal coverage of 13 sq deg, and 497 sources detected in total over the 3-24 keV energy range. There are 276 sources with spectroscopic redshifts and classifications, largely resulting from our extensive campaign of ground-based spectroscopic followup. We characterize the overall sample in terms of the X-ray, optical, and infrared source properties. The sample is primarily comprised of active galactic nuclei (AGNs), detected over a large range in redshift from z = 0.002 - 3.4 (median of <z> = 0.56), but also includes 16 spectroscopically confirmed Galactic sources. There is a large range in X-ray flux, from log( f_3-24keV / erg s^-1 cm^-2 ) ~ -14 to -11, and in rest-frame 10-40 keV luminosity, from log( L_10-40keV / erg s^-1 ) ~ 39 to 46, with a median of 44.1. Approximately 79% of the NuSTAR sources have lower energy (<10 keV) X-ray counterparts from XMM-Newton, Chandra, and Swift/XRT. The mid-infrared (MIR) analysis, using WISE all-sky survey data, shows that MIR AGN color selections miss a large fraction of the NuSTAR-selected AGN population, from ~15% at the highest luminosities (Lx > 10^44 erg s^-1) to ~80% at the lowest luminosities (Lx < 10^43 erg s^-1). Our optical spectroscopic analysis finds that the observed fraction of optically obscured AGNs (i.e., the Type 2 fraction) is F_Type2 = 53(+14-15)%, for a well-defined subset of the 8-24 keV selected sample. This is higher, albeit at a low significance level, than the Type 2 fraction measured for redshift- and luminosity-matched AGNs selected by <10 keV X-ray missions.
Most of relatively warm, unevolved metal-poor stars ($T_{\rm eff}\gtrsim 5800\,\mathrm{K}$ and $[{\rm Fe/H}]\lesssim -1.5$) exhibit almost constant lithium abundances, irrespective of metallicity or effective temperature, and thus form the so-called Spite plateau. This was originally interpreted as arising from lithium created by the Big Bang nucleosynthesis. Recent observations, however, have revealed that ultra metal-poor stars (UMP stars; $[\mathrm{Fe/H}]<-4.0$) have significantly lower lithium abundances than that of the plateau. Since most of the UMP stars are CEMP-no stars, carbon-enhanced metal-poor stars with no excess of neutron-capture elements, a connection between the carbon enhancement and lithium depletion is suspected. A straightforward approach to this question is to investigate carbon-normal UMP stars. However only one object is known in this class. As an alternative, we have determined lithium abundances for two CEMP-no main-sequence turn-off stars with metallicities $[{\rm Fe/H}]\sim -3.0$, where there are numerous carbon-normal stars with available lithium abundances that can be considered. Our 1D LTE analysis indicates that the two CEMP-no stars have lithium abundances that are consistent with values near the plateau, which suggests that carbon enhancement and lithium depletion are not directly related. Instead, our results suggest that extremely low iron abundance is a fundamental cause to depleted lithium in UMP stars.
We take advantage of the deep and wide coverage of the VST ATLAS survey to study the line-of-sight structure of the Sagittarius stellar stream in the Southern hemisphere, only ~40{\deg} away from the progenitor. We use photometrically selected Sub-Giant Branch (SGB) stars to reveal a complex debris morphology of the trailing arm and detect at least two clear peaks in the SGB distance modulus distribution. The separation between the two line-of-sight components is at least 5 kpc at the edge of the VST ATLAS footprint, but appears to change along the stream, which allows us to conclude that these detections correspond to two physically independent stellar structures, rather than a mix of co-distant stellar populations within a single stream. Our discovery of a fork in the Sgr trailing arm is verified using Blue Horizontal Branch stars and our distance measurements are calibrated using RR Lyrae stars from the Catalina Real-Time Transient Survey. Comparing with numerical simulations of the Sgr dwarf disruption, the more distant of the two components in the fork matches perfectly with the track of the trailing debris. However, no obvious counterpart exists in the simulation for the closer line-of-sight component.
The evolution of massive stars surviving the red supergiant (RSG) stage remains unexplored due to the rarity of such objects. The yellow hypergiants (YHGs) appear to be the warm counterparts of post-RSG classes located near the Humphreys-Davidson upper luminosity limit, which are characterized by atmospheric instability and high mass-loss rates. We aim to increase the number of YHGs in M33 and thus to contribute to a better understanding of the pre-supernova evolution of massive stars. Optical spectroscopy of five dust-enshrouded YSGs selected from mid-IR criteria was obtained, with the goal of detecting evidence of extensive atmospheres. We also analyzed BVI photometry for 21 of the most luminous YSGs in M33 to identify changes in the spectral type. To explore the properties of circumstellar dust, we performed SED-fitting of multi-band photometry of the 21 YSGs. We find three luminous YSGs in our sample to be YHG candidates, as they are surrounded by hot dust and are enshrouded within extended, cold dusty envelopes. Our spectroscopy of Star 2 shows emission of more than one H$\alpha$ components, as well as emission of Ca II, implying expanding structures associated with large outflow velocities. In addition, the long-term monitoring of the star reveals a dimming in the visual light curve of amplitude larger than 0.5 mag, which caused an apparent drop in the temperature that exceeded 500 K. We suggest the observed variability to be analogous to that of the Galactic YHG $\rho$ Cas. Five less luminous YSGs are suggested as post-RSG candidates showing evidence of hot or/and cool dust emission. We demonstrate that mid-IR photometry, combined with optical spectroscopy and time-series photometry, provide a robust method for identifying candidate YHGs. Future discovery of YHGs in Local Group galaxies is critical for the study of the late evolution of intermediate-mass massive stars.
Links to: arXiv, form interface, find, astro-ph, recent, 1612, contact, help (Access key information)
To elucidate the intrinsic broadband infrared (IR) emission properties of active galactic nuclei (AGN), we analyze the spectral energy distributions (SEDs) of 87 z<0.5 Palomar-Green (PG) quasars. While the Elvis AGN template with a moderate far-IR correction can reasonably match the SEDs of the AGN components in ~60% of the sample (and is superior to alternatives such as that by Assef), it fails on two quasar populations: 1) hot-dust-deficient (HDD) quasars that show very weak emission thoroughly from the near-IR to the far-IR, and 2) warm-dust-deficient (WDD) quasars that have similar hot dust emission as normal quasars but are relatively faint in the mid- and far-IR. After building composite AGN templates for these dust-deficient quasars, we successfully fit the 0.3-500 {\mu}m SEDs of the PG sample with the appropriate AGN template, an infrared template of a star-forming galaxy, and a host galaxy stellar template. 20 HDD and 12 WDD quasars are identified from the SED decomposition, including 7 ambiguous cases. Compared with normal quasars, the HDD quasars have AGN with relatively low Eddington ratios and the fraction of WDD quasars increases with AGN luminosity. Moreover, both the HDD and WDD quasar populations show relatively stronger mid-IR silicate emission. Virtually identical SED properties are also found in some quasars from z = 0.5 to 6. We propose a conceptual model to demonstrate that the observed dust deficiency of quasars can result from a change of structures of the circumnuclear tori that can occur at any cosmic epoch.
The abundance of elements synthesized by the rapid neutron-capture process (r-process elements) of extremely metal-poor (EMP) stars in the Local Group galaxies gives us clues to clarify the early evolutionary history of the Milky Way halo. The Local Group dwarf galaxies would have similarly evolved with building blocks of the Milky Way halo. However, how the chemo-dynamical evolution of the building blocks affects the abundance of r-process elements is not yet clear. In this paper, we perform a series of simulations using dwarf galaxy models with various dynamical times and total mass, which determine star-formation histories. We find that galaxies with dynamical times longer than 100 Myr have star formation rates less than $10^{-3} M_{\odot}$ yr$^{-1}$ and slowly enrich metals in their early phase. These galaxies can explain the observed large scatters of r-process abundance in EMP stars in the Milky Way halo regardless of their total mass. On the other hand, the first neutron star merger appears at a higher metallicity in galaxies with a dynamical time shorter than typical neutron star merger times. The scatters of r-process elements mainly come from inhomogeneity of the metals in the interstellar medium whereas the scatters of $\alpha$-elements are mostly due to the difference in the yield of each supernova. Our results demonstrate that the future observations of r-process elements in EMP stars will be able to constrain the early chemo-dynamical evolution of the Local Group galaxies.
The CIII]{\lambda}{\lambda}1907, 1909 rest-frame UV emission doublet has recently been detected in galaxies during the epoch of reionization (z > 6), with high equivalent width (>10 {\AA}, rest frame). Currently, it is possible to obtain much more detailed information for star-forming galaxies at significantly lower redshift. Accordingly, studies of their far-UV spectra are useful for understanding the factors modulating the strength of CIII] emission. We present the first statistical sample of CIII] emission measurements in star-forming galaxies at z ~ 1. Our sample is drawn from the DEEP2 survey and spans the redshift 0.64 < z < 1.35 (<z> = 1.08). We find that the median equivalent width (EW) of individual CIII] detections in our sample (1.30 {\AA}) is much smaller than the typical value observed thus far at z > 6. Furthermore, out of 184 galaxies with coverage of CIII], only 40 have significant detections. Galaxies with individual CIII] detections have bluer colors and lower luminosities on average than those without, implying that strong CIII] emitters are in general young and low-mass galaxies without significant dust extinction. Using stacked spectra, we further investigate how CIII] strength correlates with multiple galaxy properties (M_B, U-B, M*, star-formation rate, specific star-formation rate) and rest-frame near-UV (FeII* and MgII) and optical [OIII] and H{\beta}) emission line strengths. These results provide a detailed picture of the physical environment in star-forming galaxies at z ~ 1, and motivate future observations of strong CIII] emitters at similar redshifts.
We present the results of clustering analyses of Lyman break galaxies (LBGs) at $z \sim 3$, $4$, and $5$ using the final data release of the Canada--France--Hawaii Telescope Legacy Survey (CFHTLS). Deep- and wide-field images of the CFHTLS Deep Survey enable us to obtain sufficiently accurate two-point angular correlation functions to apply a halo occupation distribution analysis. The mean halo masses, calculated as $\langle M_{h} \rangle = 10^{11.7} - 10^{12.8} h^{-1} M_{\odot}$, increase with stellar-mass limit of LBGs. The threshold halo mass to have a central galaxy, $M_{{\rm min}}$, follows the same increasing trend with the low-$z$ results, whereas the threshold halo mass to have a satellite galaxy, $M_{1}$, shows higher values at $z = 3 - 5$ than $z = 0.5 - 1.5$ over the entire stellar mass range. Along with considering the low satellite fractions in high-$z$, these results suggest that satellite galaxies form inefficiently within dark haloes at $z=3-5$ even for less massive satellites with $M_{\star} < 10^{10} M_{\odot}$. We compute stellar-to-halo mass ratios (SHMRs) assuming a main sequence of galaxies, which is found to provide consistent SHMRs with those derived from a spectral energy distribution fitting method. The observed SHMRs are in good agreement with the model predictions based on the abundance-matching method within $1\sigma$ confidence intervals. We show, for the first time observationally, the increasing trend of $M_{{\rm h}}^{{\rm pivot}}$, which is the halo mass with a peak SHMR, with cosmic time at $z > 3$, and with keeping the constant star-formation efficiency indicates that mass growth rates of stellar components and dark haloes are comparable at $3 < z < 5$.
Numerical simulations have become a major tool for understanding galaxy formation and evolution. Over the decades the field has made significant progress. It is now possible to simulate the formation of individual galaxies and galaxy populations from well defined initial conditions with realistic abundances and global properties. An essential component of the calculation is to correctly estimate the inflow to and outflow from forming galaxies since observations indicating low formation efficiency and strong circum-glactic presence of gas are persuasive. Energetic 'feedback' from massive stars and accreting super-massive black holes - generally unresolved in cosmological simulations - plays a major role for driving galactic outflows, which have been shown to regulate many aspects of galaxy evolution. A surprisingly large variety of plausible sub-resolution models succeeds in this exercise. They capture the essential characteristics of the problem, i.e. outflows regulating galactic gas flows, but their predictive power is limited. In this review we focus on one major challenge for galaxy formation theory: to understand the underlying physical processes that regulate the structure of the interstellar medium, star formation and the driving of galactic outflows. This requires accurate physical models and numerical simulations, which can precisely describe the multi-phase structure of the interstellar medium on the currently unresolved few hundred parsecs scales of large scale cosmological simulations. Such models ultimately require the full accounting for the dominant cooling and heating processes, the radiation and winds from massive stars and accreting black holes, an accurate treatment of supernova explosions as well as the non-thermal components of the interstellar medium like magnetic fields and cosmic rays.
We present a sample of 34 spectroscopically confirmed BzK-selected ~1e11 Msun quiescent galaxies (pBzK) in the COSMOS field. The targets were initially observed with VIMOS on the VLT to facilitate the calibration of the photometric redshifts of massive galaxies at z >~ 1.5. Here we describe the reduction and analysis of the data, and the spectrophotometric properties of these pBzK galaxies. In particular, using a spatially resolved median 2D spectrum, we find that the fraction of stellar populations with ages <1 Gyr is at least 3 times higher in the outer regions of the pBzK galaxies than in their cores. This results in a mild age gradient of ~<0.4 Gyr over ~6 kpc and suggests either the occurrence of widespread rejuvenation episodes or that inside-out quenching played a role in the passivization of this galaxy population. We also report on low-level star formation rates derived from the [OII]3727A emission line, with SFR_OII = 3.7-4.5 Msun/yr. This estimate is confirmed by an independent measurement on a separate sample of similarly-selected quiescent galaxies in the COSMOS field, using stacked far-infrared data (SFR_FIR = 2-4 Msun/yr). This second, photometric sample also displays significant excess at 1.4 GHz, suggestive of the presence of radio-mode AGN activity.
The unification scheme of active galactic nuclei (AGNs) invokes an optically thick molecular torus component hiding the broad emission line region. Assuming the presence of a thick neutral component in the molecular torus characterized by a \ion{H}{I} column density > $10^{22}{\rm\ cm^{-2}}$, we propose that far UV radiation around Ly$\alpha$ can be significantly polarized through Rayleigh scattering. Adopting a Monte Carlo technique we compute polarization of Rayleigh scattered radiation near Ly$\alpha$ in a thick neutral region in the shape of a slab and a cylindrical shell. It is found that radiation near Ly$\alpha$ Rayleigh reflected from a very thick slab can be significantly polarized in a fairly large range of wavelength $\Delta\lambda\sim 50$ \AA\ exhibiting a flux profile similar to the incident one. Rayleigh transmitted radiation in a slab is characterized by the central dip with a complicated polarization behavior. The optically thick part near Ly$\alpha$ center is polarized in the direction perpendicular to the slab normal, which is in contrast to weakly polarized wing parts in the direction parallel to the slab normal. A similar polarization flip phenomenon is also found in the case of a tall cylindrical shell, in which the spatial diffusion along the vertical direction near the inner cylinder wall for core photons leads to a tendency of the electric field aligned to the direction perpendicular to the vertical axis. Observational implications are briefly discussed including spectropolarimetry of the quasar PG~1630+377 by Koratkar et al. in 1990 where Ly$\alpha$ is strongly polarized with no other emission lines polarized.
We present a combined experimental and theoretical study focussing on the quantum tunneling of atoms in the reaction between CH4 and OH. The importance of this reaction pathway is derived by investigating isotope substituted analogs. Quantitative reaction rates needed for astrochemical models at low temperature are currently unavailable both in the solid state and in the gas phase. Here, we study tunneling effects upon hydrogen abstraction in CH4 + OH by focusing on two reactions: CH4 + OD -> CH3 + HDO and CD4 + OH -> CD3 + HDO. The experimental study shows that the solid-state reaction rate R(CH4 + OD) is higher than R(CD4 + OH) at 15 K. Experimental results are accompanied by calculations of the corresponding unimolecular and bimolecular reaction rate constants using instanton theory taking into account surface effects. From the work presented here, the unimolecular reactions are particularly interesting as these provide insight into reactions following a Langmuir-Hinshelwood process. The resulting ratio of the rate constants shows that the H abstraction (k(CH4 + OD)) is approximately ten times faster than D-abstraction (k(CD4 + OH)) at 65 K. We conclude that tunneling is involved at low temperatures in the abstraction reactions studied here. The unimolecular rate constants can be used by the modeling community as a first approach to describe OH-mediated abstraction reactions in the solid phase. For this reason we provide fits of our calculated rate constants that allow the inclusion of these reactions in models in a straightforward fashion.
The time domain is the emerging forefront of astronomical research with new facilities and instruments providing unprecedented amounts of data on the temporal behavior of astrophysical populations. Dealing with the size and complexity of this requires new techniques and methodologies. Quasars are an ideal work set for developing and applying these: they vary in a detectable but not easily quantifiable manner whose physical origins are poorly understood. In this paper, we will review how quasars are identified by their variability and how these techniques can be improved, what physical insights into their variability can be gained from studying extreme examples of variability, and what approaches can be taken to increase the number of quasars known. These will demonstrate how astroinformatics is essential to discovering and understanding this important population.
We report Spitzer Space Telescope IRAC 3.6, 4.5, 5.8 and 8 um and MIPS 24 and 70 um observations of the 32 Ori Group, a recently discovered nearby stellar association situated towards northern Orion. The proximity of the group (~93 pc) has enabled a sensitive search for circumstellar dust around group members, and its age (~20 Myr) corresponds roughly to an epoch thought to be important for terrestrial planet formation in our own solar system. We quantify infrared excess emission due to circumstellar dust among group members, utilizing available optical (e.g. Hipparcos, Tycho) and near-IR (2MASS) photometry in addition to the Spitzer IR photometry. We report 4 out of the 14 objects which exhibit 24 um excess emission more than 4\sigma above the stellar photosphere (>20%) though lacking excess emission at shorter wavelengths: HD 35656 (A0Vn), HD 36338 (F4.5), RX J0520.5+0616 (K3), and HD 35499 (F4). Two objects (HD 35656 and RX J0520.0+0612) have 70 um excesses, although the latter lacks 24 um excess emission. The 24 um disk fraction of this group is 29(+14,-9%), which is similar to previous findings for groups of comparable ages and places 32 Ori as the young stellar group with the 2nd most abundant 24 um excesses among groups lacking accreting T Tauri stars (behind only the approximately coeval Beta Pic Moving Group). We also model the infrared excess emission using circumstellar dust disk models, placing constraints on disk parameters including L_IR/L_*, T_disk, characteristic grain distance, and emitting area. The L_IR/L_* values for all the stars can be reasonably explained by steady state disk evolution.
In 2015 July 29 - September 1 the satellite XMM-Newton pointed at the BL Lac object PG 1553+133 six times, collecting data for 218 hours. During one of these epochs, simultaneous observations by the Swift satellite were requested to compare the results of the X-ray and optical-UV instruments. Optical, near-infrared and radio monitoring was carried out by the Whole Earth Blazar Telescope (WEBT) collaboration for the whole observing season. We here present the results of the analysis of all these data, together with an investigation of the source photometric and polarimetric behaviour over the last three years. The 2015 EPIC spectra show slight curvature and the corresponding light curves display fast X-ray variability with a time scale of the order of 1 hour. In contrast to previous results, during the brightest X-ray states detected in 2015 the simple log-parabolic model that best-fits the XMM-Newton data also reproduces reasonably well the whole synchrotron bump, suggesting a peak in the near-UV band. We found evidence of a wide rotation of the polarization angle in 2014, when the polarization degree was variable, but the flux remained almost constant. This is difficult to interpret with deterministic jet emission models, while it can be easily reproduced by assuming some turbulence of the magnetic field.
The detection and study of molecular gas in born-again stars would be of great importance to understand their composition and chemical evolution. In addition, the molecular emission would be an invaluable tool to explore the physical conditions, kinematics and formation of asymmetric structures in the circumstellar envelopes of these evolved stars. However, until now, all attempts to detect molecular emission from the cool material around born-again stars have failed. We carried out observations using the APEX and IRAM 30m telescopes to search for molecular emission toward four well studied born-again stars, V4334 Sgr, V605 Aql, A30 and A78, that are thought to represent an evolutionary sequence. We detected for the first time emission from HCN and H$^{13}$CN molecules toward V4334 Sgr, and CO emission in V605 Aql. No molecular emission was detected above the noise level toward A30 and A78. A first estimate of the H$^{12}$CN/H$^{13}$CN abundance ratio in the circumstellar environment of V4334 Sgr is $\approx$3, which is similar to the value of the $^{12}$C/$^{13}$C ratio measured from other observations. We derived a rotational temperature of $T_{\rm rot}$=13$\pm1$ K, and a total column density of $N_{{\rm HCN}}$=1.6$\pm0.1\times$10$^{16}$ cm$^{-2}$ for V4334 Sgr. This result sets a lower limit on the amount of hydrogen that was ejected into the wind during the born-again event of this source. For V605 Aql, we obtained a lower limit for the integrated line intensities $I_{^{12}\rm C}$/$I_{^{13}\rm C}$>4.
In this work we study the dynamics of the Local Group (LG) within the context of cosmological models beyond General Relativity (GR). Using observable kinematic quantities to identify candidate pairs we build up samples of simulated LG-like objects drawing from $f(R)$, symmetron, DGP and quintessence N-body simulations together with their $\Lambda$CDM counterparts featuring the same initial random phase realisations. The variables and intervals used to define LG-like objects are referred to as Local Group model; different models are used throughout this work and adapted to study their dynamical and kinematic properties. The aim is to determine how well the observed LG-dynamics can be reproduced within cosmological theories beyond GR, We compute kinematic properties of samples drawn from alternative theories and $\Lambda$CDM and compare them to actual observations of the LG mass, velocity and position. As a consequence of the additional pull, pairwise tangential and radial velocities are enhanced in modified gravity and coupled dark energy with respect to $\Lambda$CDM, inducing significant changes to the total angular momentum and energy of the LG. For example, in models such as $f(R)$ and the symmetron this increase can be as large as $60\%$, peaking well outside of the $95\%$ confidence region allowed by the data. This shows how simple considerations about the LG dynamics can lead to clear small-scale observational signatures for alternative scenarios, without the need of expensive high-resolution simulations.
We present the results of a study of the stellar and accretion properties of the (almost) complete sample of class II and transitional YSOs in the Lupus I, II, III and IV clouds, based on spectroscopic data acquired with the VLT/X-Shooter spectrograph. Our study combines the dataset from our previous work with new observations of 55 additional objects. We have investigated 92 YSO candidates in total, 11 of which have been definitely identified with giant stars unrelated to Lupus. The stellar and accretion properties of the 81 bona fide YSOs, which represent more than 90% of the whole class~II and transition disc YSO population in the aforementioned Lupus clouds, have been homogeneously and self-consistently derived, allowing for an unbiased study of accretion and its relationship with stellar parameters. The accretion luminosity, Lacc, increases with the stellar luminosity, Lstar, with an overall slope of ~1.6, similar but with a smaller scatter than in previous studies. There is a significant lack of strong accretors below Lstar~0.1Lsun, where Lacc is always lower than 0.01Lstar. We argue that the Lacc-Lstar slope is not due to observational biases, but is a true property of the Lupus YSOs. The logMacc-logMstar correlation shows a statistically significant evidence of a break, with a steeper relation for Mstar<0.2Msun and a flatter slope for higher masses. The bimodality of the Macc-Mstar relation is confirmed with four different evolutionary models used to derive the stellar mass. The bimodal behaviour of the observed relationship supports the importance of modelling self-gravity in the early evolution of the more massive discs, but other processes, such as photo evaporation and planet formation during the YSO's lifetime, may also lead to disc dispersal on different timescales depending on the stellar mass. We also refined the empirical Lacc vs. Lline relationships.
A mysterious type of matter is supposed to exist, because the observed rotational velocity curves of particle moving around the galactic center and the expected rotational velocity curves do not match. There are also a number of proposals in the modified gravity for this discrepancy. In this contrast, in $2008$, B$\ddot{\text{o}}$hmer et al. presented an interesting idea in (Astropart Phys 29(6):386-392, 2008) where they showed that a $f(\mathcal{R})$ gravity model could actually explain dark matter to be a geometric effect only. They solved the gravitational field equations in vacuum using generic $f(\mathcal{R})$ gravity model for constant velocity regions and found that the resulting modifications in the Einstein-Hilbert Lagrangian is of the form $\mathcal{R}^{1+m}$, where $m=V_{tg}^2/c^2$; $V_{tg}$ being the tangential velocity of the test particle moving around galactic dark matter region and, $c$, the speed of light. From observations it is known that $m\approx\mathcal{O}(10^{-6})$. In this article, we perform two things (1) We show that the form of $f(\mathcal{R})$ they claimed is not correct. In doing the calculations, we found that when the radial component of the metric for constant velocity regions is a constant then the exact solutions for $f(\mathcal{R})$ obtained is of the form of $\mathcal{R}^{1-\alpha}$ which corresponds to a negative correction rather than positive, $\alpha$ is a function of $m$. (2) We also show that we can not have an analytic solution of $f(\mathcal{R})$ for all values of tangential velocity including the observed value of tangential velocity $200-300$Km/s if the radial coefficient of the metric which describes the dark matter regions is \emph{not a constant}. Thus, we have to rely on the numerical solutions to get an approximate model for dark matter in $f(\mathcal{R})$ gravity.
Several starburst galaxies have been observed in the GeV and TeV bands; in this regime, gamma rays are mainly produced through the pionic process from cosmic-ray proton interactions with the interstellar medium. Furthermore, the dense environments of starbursts have been proposed to act as proton "calorimeters" in which collisions dominate proton losses, so that a substantial fraction of cosmic-ray proton energy input is emitted in gamma rays. Here we build a one-zone, "thick-target" model that implements calorimetry and thus places a firm upper bound to the gamma-ray emission due to cosmic-ray interactions. The model assumes that cosmic rays are accelerated by supernovae, and escape is neglected. Our model has only two free parameters: the cosmic-ray proton acceleration energy per supernova, and the proton injection spectral index. We find that in our thick-target limit, the emergent gamma-ray spectral index is the same as the injection index. We calculate the pionic gamma-ray emission from 10 MeV to 10 TeV for five starburst galaxies (M82, NGC 253, NGC 4945, NGC 1068 and Circinus) and the ultraluminous infrared galaxy Arp 220. Using Fermi, H.E.S.S., and VERITAS data, we perform chi-squared fits to derive best-fit parameters for each galaxy. We find that for the best-observed galaxies M82 and NGC 253, the cosmic-ray acceleration energetics and spectral indices are consistent with Fermi observations of Galactic supernova remnants; this suggests that cosmic-ray acceleration by supernovae is similar in starbursts and in our Galaxy. For NGC 4945 and NGC 1068, the models are not well-constrained due to the lack of TeV data, but the GeV data also are consistent with near or full proton calorimetry. Consequently, we conclude that these starbursts, and perhaps most starbursts, are consistent with near or total proton calorimetry. (Abridged)
Links to: arXiv, form interface, find, astro-ph, recent, 1612, contact, help (Access key information)
Abridged - We quantify the effect of the galaxy group environment (for 12.5 < log(M_group/Msun) < 14.0) on the star formation rates of the (morphologically-selected) population of disk-dominated local Universe spiral galaxies (z < 0.13) with stellar masses log(M*/Msun) > 9.5. Within this population, we find that, while a small minority of group satellites are strongly quenched, the group centrals, and the large majority of satellites exhibit levels of SFR indistinguishable from ungrouped "field" galaxies of the same M*, albeit with a higher scatter, and for all M*. Modelling these results, we deduce that disk-dominated satellites continue to be characterized by a rapid cycling of gas into and out of their ISM at rates similar to those operating prior to infall, with the on-going fuelling likely sourced from the group intrahalo medium (IHM) on Mpc scales, rather than from the circum-galactic medium on 100kpc scales. Consequently, the color-density relation of the galaxy population as a whole would appear to be primarily due to a change in the mix of disk- and spheroid-dominated morphologies in the denser group environment compared to the field, rather than to a reduced propensity of the IHM in higher mass structures to cool and accrete onto galaxies. We also suggest that the inferred substantial accretion of IHM gas by satellite disk-dominated galaxies will lead to a progressive reduction in their specific angular momentum, thereby representing an efficient secular mechanism to transform morphology from star-forming disk-dominated types to more passive spheroid-dominated types.
The dwarf galaxies around the Milky Way are distributed in a so-called vast polar structure (VPOS) that may be in conflict with Lambda CDM simulations. Here, we seek to determine if the VPOS poses a serious challenge to the Lambda CDM paradigm on galactic scales. Specifically, we investigate if the VPOS remains coherent as a function of time. Using the measured HST proper motions and associated uncertainties, we integrate the orbits of the classical Milky Way satellites backwards in time and find that the structure disperses well before a dynamical time. We also examine in particular Leo I and Leo II using their most recent proper motion data, both of which have extreme kinematic properties, but these satellites do not appear to drive the polar fit that is seen at the present day. We have studied the effect of the uncertainties on the HST proper motions on the coherence of the VPOS as a function of time. We find that eight of the eleven classical dwarfs have reliable proper motions; for these eight, the VPOS also loses significance in less than a dynamical time, indicating that the VPOS is not a dynamically stable structure. Obtaining more accurate proper motion measurements of Ursa Minor, Sculptor, and Carina would bolster these conclusions.
We present high-resolution (1".0) Atacama Large Millimeter/submillimeter Array (ALMA) observations of CO (1-0) and CO (2- 1) rotational transitions toward the nearby IR-luminous merger NGC 1614 supplemented with ALMA archival data of CO (3-2), and CO (6-5) transitions. The CO (6-5) emission arises from the starburst ring (central 590 pc in radius), while the lower-$J$ CO lines are distributed over the outer disk ($\sim$ 3.3 kpc in radius). Radiative transfer and photon dominated region (PDR) modeling reveal that the starburst ring has a single warmer gas component with more intense far-ultraviolet radiation field ($n_{\rm{H_2}}$ $\sim$ 10$^{4.6}$ cm$^{-3}$, $T_{\rm{kin}}$ $\sim$ 42 K, and $G_{\rm{0}}$ $\sim$ 10$^{2.7}$) relative to the outer disk ($n_{\rm{H_2}}$ $\sim$ 10$^{5.1}$ cm$^{-3}$, $T_{\rm{kin}}$ $\sim$ 22 K, and $G_{\rm{0}}$ $\sim$ 10$^{0.9}$). A two-phase molecular interstellar medium with a warm and cold ($>$ 70 K and $\sim$ 19 K) component is also an applicable model for the starburst ring. A possible source for heating the warm gas component is mechanical heating due to stellar feedback rather than PDR. Furthermore, we find evidence for non-circular motions along the north-south optical bar in the lower-$J$ CO images, suggesting a cold gas inflow. We suggest that star formation in the starburst ring is sustained by the bar-driven cold gas inflow, and starburst activities radiatively and mechanically power the CO excitation. The absence of a bright active galactic nucleus can be explained by a scenario that cold gas accumulating on the starburst ring is exhausted as the fuel for star formation, or is launched as an outflow before being able to feed to the nucleus.
We investigate the angular momentum evolution of four disk galaxies residing in Milky Way-sized halos formed in cosmological zoom-in simulations with various sub-grid physics and merging histories. We decompose these galaxies kinematically and photometrically, into their disk and bulge components. The simulated galaxies and their components lie on the observed sequences in the $j_*$--$M_*$ diagram relating the specific angular momentum and mass of the stellar component. We find that galaxies in low-density environments follow the relation $j_* \propto M_*^{\alpha}$ past major mergers, with $\alpha \sim 0.6$ in the case of strong feedback, when bulge-to-disk ratios are relatively constant, and $\alpha \sim 1.4$ in the other cases, when secular processes operate on shorter timescales. We compute the retention factors (i.e. the ratio of the specific angular momenta of stars and dark matter) for both disks and bulges and show that they vary relatively slowly after averaging over numerous but brief fluctuations. For disks, the retention factors are usually close to unity, while for bulges, they are a few times smaller. Our simulations therefore indicate that galaxies and their halos grow in a quasi-homologous way.
We present a multi-wavelength photometric catalog in the COSMOS field as part of the observations by the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS). The catalog is based on Hubble Space Telescope Wide Field Camera 3 (HST/WFC3) and Advanced Camera for Surveys (ACS) observations of the COSMOS field (centered at RA: $10^h00^m28^s$, Dec:$+02^{\circ}12^{\prime}21^{\prime\prime}$). The final catalog has 38671 sources with photometric data in forty two bands from UV to the infrared ($\rm \sim 0.3-8\,\mu m$). This includes broad-band photometry from the HST, CFHT, Subaru, VISTA and Spitzer Space Telescope in the visible, near infrared and infrared bands along with intermediate and narrow-band photometry from Subaru and medium band data from Mayall NEWFIRM. Source detection was conducted in the WFC3 F160W band (at $\rm 1.6\,\mu m$) and photometry is generated using the Template FITting algorithm. We further present a catalog of the physical properties of sources as identified in the HST F160W band and measured from the multi-band photometry by fitting the observed spectral energy distributions of sources against templates.
A number of planetary nebulae show binary central stars and significant abundance discrepancies between values estimated from colisionally excited lines when compared to the same abundances estimated from recombination lines. One approach to investigate this yet unsolved problem is using spatially resolved images of emission lines in an attempt to detect a possibly distinct metal rich component in the nebula. In this work we present results of spatially resolved bundance analysis of NGC 6778 based on data gathered from VLT VIMOS-IFU. We discuss the spatial variations found as well as possible limitations of the method in answering questions about abundance variations.
Jet powers in many radio galaxies with extended radio structures appear to exceed their associated accretion luminosities. In systems with very low accretion rates, this is likely due to the very low accretion luminosities resulting from radiatively inefficient accretion flows. In systems with high accretion rates, the accretion flows are expected to be radiatively efficient, and the production of such powerful jets may require an accretion scenario which involves magnetically arrested discs (MADs). However, numerical simulations of the MAD scenario indicate that jet production efficiency is large only for geometrically thick accretion flows and scales roughly with $(H/R)^2$, where $H$ is the disc height and $R$ is the distance from the BH. Using samples of FRII radio galaxies and quasars accreting at moderate accretion rates we show that their jets are much more powerful than predicted by the MAD scenario. We discuss possible origins of this discrepancy, suggesting that it can be related to approximations adopted in MHD simulations to treat optically thick accretion flow within the MAD-zone, or may indicate that accretion disks are geometrically thicker than the standard theory predicts.
We report the broadband X-ray spectra of the ultra-luminous infrared galaxy (ULIRG) UGC 5101 in the 0.25-100 keV band observed with Swift/Burst Alert Telescope (BAT), NuSTAR, Suzaku, XMM-Newton, and Chandra. A Compton-thick AGN obscured with a hydrogen column density of $\approx 1.3\times10^{24}$ cm$^{-2}$ is detected above 10 keV. A spectral fit with a numerical torus model favors a large half opening angle of the torus, $>41$ degrees, suggesting that the covering fraction of material heavily obscuring the X-ray source is moderate. The intrinsic 2-10 keV luminosity is determined to be $\approx 1.4\times 10^{43}$ erg s$^{-1}$, which is $\approx$2.5 times larger than the previous estimate using only data below 10 keV with a simple spectral model. We find that UGC 5101 shows the ratio between the [O IV] 26 $\mu$m line and 2-10 keV luminosities similar to those of normal Seyfert galaxies, along with other ULIRGs observed with NuSTAR, indicating that a significant portion of local ULIRGs are not really "X-ray faint" with respect to the flux of forbidden lines originating from the narrow line region (NLR). We propose a possible scenario that (1) the AGN in UGC 5101 is surrounded not only by Compton-thick matter located close to the equatorial plane but also by Compton-thin ($N_\mathrm{H} \sim 10^{21}$ cm$^{-2}$) matter in the torus-hole region and (2) it is accreting at a high Eddington rate with a steep UV to X-ray spectral energy distribution. Nevertheless, we argue that AGNs in many ULIRGs do not look extraordinary (i.e., extremely X-ray faint), as suggested by recent works, compared with normal Seyferts.
The gamma-ray blazar OJ 287 was in a high activity state during December 2015 - February 2016. Coinciding with this high brightness state, we observed this source for photometry on 40 nights in R-band and for polarimetry on 9 epochs in UBVRI bands. During the period of our observations, the source brightness varied between $13.20 \pm 0.04$ to $14.98 \pm 0.04$ mag and the degree of polarization (P ) fluctuated between $6.0 \pm 0.3$% and $28.3 \pm 0.8$% in R-band. Focusing on intra-night optical variability (INOV), we find a duty cycle of about 71% using $\chi^2$-statistics, similar to that known for blazars. From INOV data, the shortest variability time scale is estimated to be $142 \pm 38$ min yielding a lower limit of the observed Doppler factor $\delta_0 = 1.17$, the magnetic field strength $B \le 3.8$ G and the size of the emitting region Rs < $2.28 \times 10^{14}$ cm. On inter-night timescales, a significant anti-correlation between R-band flux and P is found. The observed P at U-band is generally larger than that observed at longer wavelength bands suggesting a wavelength dependent polarization. Using V -band photometric and polarimetric data from Steward Observatory obtained during our monitoring period we find a varied correlation between P and V-band brightness. While an anticorrelation is seen between P and V -band mag at sometimes, no correlation is seen at other times, thereby, suggesting the presence of more than one short-lived shock components in the jet of OJ 287.
We present the summary of our last results on the spatial distribution and relative frequencies of Supernovae (SNe) in a large number of host galaxies from the Sloan Digital Sky Survey (SDSS). We use the locations of SNe in order to study the relations between radial/azimuthal distributions of SNe and properties of their hosts and environments. On the other hand, the vertical distribution of SNe allows to study the progenitors association to the thin or thick discs, and to the stellar halo. We also propose the underlying mechanisms shaping the number ratios of SNe types. It is important to note that there were no extended studies of the 3D distribution of SNe and structural parameters of hosts. Our study is intended to fill this gap and better constrain the nature of SN progenitors.
Dust modeling is crucial to infer dust properties and budget for galaxy studies. However, there are systematic disparities between dust grain models that result in corresponding systematic differences in the inferred dust properties of galaxies. Quantifying these systematics requires a consistent fitting analysis. We compare the output dust parameters and assess the differences between two dust grain models, Compiegne et al (2011), and THEMIS (Jones et al 2013, Kohler et al 2015). In this study, we use a single fitting method applied to all the models to extract a coherent and unique statistical analysis. We fit the models to the dust emission seen by Spitzer and Herschel in the Small and Large Magellanic Clouds (SMC and LMC). The observations cover the infrared (IR) spectrum from a few microns to the sub-millimeter range. For each fitted pixel, we calculate the full n-D likelihood, based on the method described in Gordon et al (2014). The free parameters are both environmental ($U$, the interstellar radiation field strength; $\alpha_\mathrm{ISRF}$, power-law coefficient for a multi-U environment; $\Omega^*$, the starlight strength) and intrinsic to the model ($Y_\mathrm{i}$: abundances of the grain species $i$; $\alpha_\mathrm{sCM20}$, coefficient in the small carbon grain size distribution). Fractional residuals of 5 different sets of parameters show that fitting THEMIS brings a more accurate reproduction of the observations than the Compiegne model. However, independent variations of the dust species show strong model-dependencies. We find that the abundance of silicates can only be constrained to an upper-limit and the silicate/carbon ratio is different than that seen in our Galaxy. In the LMC, our fits result in dust masses slightly lower than those found in literature , by a factor lower than 2. In the SMC, we find dust masses in agreement with previous studies.
Extragalactic observations of water emission can provide valuable insights into the excitation of the interstellar medium. In addition, extragalactic megamasers are powerful probes of kinematics close to active nuclei. Therefore, it is paramount to determine the true origin of the water emission, whether it is excited by processes close to an AGN or in star-forming regions. We use ALMA Band 5 science verification observations to analyse the emission of the 183 GHz water line in Arp 220 on sub-arcsecond scales, in conjunction with new ALMA Band 7 data at 325 GHz. Specifically, the nature of the process leading to the excitation of emission at these water lines is studied in this context. Supplementary 22 GHz VLA observations are used to better constrain the parameter space in the excitation modelling of the water lines. We detect 183 GHz H2O and 325 GHz water emission towards the two compact nuclei at the center of Arp 220, being brighter in Arp 220 West. The emission at these two frequencies is compared to previous single-dish data and does not show evidence of variability. The 183 and 325 GHz lines show similar spectra and kinematics, but the 22 GHz profile is significantly different in both nuclei due to a blend with an NH3 absorption line. Our findings suggest that the most likely scenario to cause the observed water emission in Arp 220 is a large number of independent masers originating from numerous star-forming regions.
Comparison analyses between the gas emission data (HI 21 cm line and CO 2.6 mm line) and the Planck/IRAS dust emission data (optical depth at 353 GHz tau353 and dust temperature Td) allow us to estimate the amount and distribution of the hydrogen gas more accurately, and our previous studies revealed the existence of a large amount of optically-thick HI gas in the solar neighborhood. Referring to this, we discuss the neutral hydrogen gas around the Perseus cloud in the present paper. By using the J-band extinction data, we found that tau353 increases as a function of the 1.3-th power of column number density of the total hydrogen (NH), and this implies dust evolution in high density regions. This calibrated tau353-NH relationship shows that the amount of the HI gas can be underestimated to be ~60% if the optically-thin HI method is used. Based on this relationship, we calculated optical depth of the 21 cm line (tauHI), and found that <tauHI> ~ 1.1 around the molecular cloud. The effect of tauHI is still significant even if we take into account the dust evolution. We also estimated a spatial distribution of the CO-to-H2 conversion factor (XCO), and we obtained <XCO> ~ 9.5x10^19 cm-2 K-1 km-1 s. Although these results are inconsistent with some previous studies, these discrepancies can be well explained by the difference of the data and analyses methods.
Active Galactic Nuclei (AGN) are perfect laboratories to check General Relativity (GR) effects by using Broad Line Region (BLR) clouds eclipses to probe the innermost regions of the accretion disk. A new relativistic X-ray spectral model for X-ray eclipses is introduced. First we present the different observables that are involved in X-ray eclipses, including the X-ray emitting regions size, the emissivity index, the cloud's column density, ionization, size and velocity, the black hole spin, and the system's inclination. Then we highlight some theoretical predictions on the observables by using XMM-Newton simulations, finding that absorption varies depending on the photons' energy range, being maximum when the approaching side of the X-ray-emitting region is covered. Finally, we fit our relativistic model to actual XMM-Newton data from a long observation of the NLS1 galaxy SWIFT J2127.4+5654, and compare our results with a previous work, in which we addressed the BLR cloud eclipse from a non-relativistic prespective.
The Milky Way galaxy is a typical spiral galaxy which consists of a black hole in its centre, a barred bulge and a disk which contains spiral arms. The complex structure of the Galaxy makes it extremely difficult and challenging to model its mass distribution, particularly for the Galactic disk which plays the most important role in the dynamics and evolution of the Galaxy. Conventionally an axisymmetric disk model with an exponential brightness distribution and a constant mass-to-light ratio is assumed for the Galactic disk. In order to generate a flat rotation curve, a dark halo has also to be included. Here, by using the recently released Gaia billion-star map, we propose a Galactic disk mass distribution model which is based on the star density distribution rather than the brightness and mass-to-light ratio. The model is characterized by two parameters, a bulge radius and a characteristic length. Using the mass distribution model and solving the Poisson equation of the Galaxy, we obtain a flat rotation curve which reproduces the key observed features with no need for a dark halo.
Outflows in active galaxies (AGNs) are common, although their launching mechanism, location, and physical impact on the host galaxy remain controversial. We conducted a multiwavelength six-month campaign to observe the nearby Seyfert galaxy NGC7469 with several observatories in order to better understand and quantify the outflow in this AGN. We report on the time-integrated line-resolved X-ray spectrum of NGC7469 obtained with the Reflection Grating Spectrometer (RGS) on board XMM-Newton. We use the RGS spectrum to discern the many AGN outflow components. A global fit is applied to obtain their physical parameters. We find that the AGN wind can be well described by three narrow velocity components at -650, -950, and -2050 kms-1. The RGS clearly resolves the -20 50 kms-1 component in C5+ Ly, while the -650 kms-1 and -950 kms-1 velocities are blended. Similar velocities are resolved in the UV. The H-equivalent column densities of these components are, respectively, NH 7 10 20, 2.2 10 21, and 10 20 cm-2, for a total of 3 10 21 cm-2. The -650 kms-1 component shows a broad ionisation distribution. We identify a photo-ionised emission component blue-shifted by -450 kms-1 which we ascribe to the same outflow that produces the absorption lines. The elemental abundance ratios of C, N, Ne, S, and Fe to O in the outflow tend to be between 1 - 2 times solar. Preliminary estimates of the absorber distance from the AGN center suggest it is at least a few pc away from the center, but more advanced methods need to be applied in order to obtain better constraints. The complex X-ray spectrum of NGC 7469 demonstrates the richness of high energy phenomena taking place in AGN cores. The subtle spectroscopic differences between the various components require deep, high-resolution observations, such as the present RGS spectrum, if one is to resolve them and perform quantitative plasma diagnostics.
The aim of this paper is to characterise the abundance patterns of five iron-peak elements (Mn, Fe, Ni, Cu, and Zn) for which the stellar origin and chemical evolution are still debated. We automatically derived iron peak (Mn, Fe, Ni, Cu, and Zn) and alpha element (Mg) chemical abundances for 4666 stars. We used the bimodal distribution of [Mg/Fe] to chemically classify sample stars into different Galactic substructures: thin disc, metal-poor and high-alpha metal rich, high-alpha and low-alpha metal-poor populations. High-alpha and low-alpha metal-poor populations are fully distinct in Mg, Cu, and Zn. Thin disc trends of [Ni/Fe] and [Cu/Fe] are very similar and show a small increase at supersolar metallicities. Thin and thick disc trends of Ni and Cu are very similar and indistinguishable. Mn looks different from Ni and Cu. [Mn/Fe] trends of thin and thick discs actually have noticeable differences: the thin disc is slightly Mn richer than the thick disc. [Zn/Fe] trends look very similar to those of [alpha/Fe] trends. The dispersion of results in both discs is low (approx 0.05 dex for [Mg, Mn, and Cu/Fe]) and is even much lower for [Ni/Fe] (approx 0.035 dex). Zn is an alpha-like element and could be used to separate thin and thick disc stars. [Mn/Mg] ratio could also be a very good tool for tagging Galactic substructures. Some models can partially reproduce the observed Mg, Zn, and, Cu behaviours. Models mostly fail to reproduce Mn and Ni in all metallicity domains, however, models adopting yields normalised from solar chemical properties reproduce Mn and Ni better, suggesting that there is still a lack of realistic theoretical yields of some iron-peak elements. Very low scatter (approx 0.05 dex) in thin and thick disc sequences could provide an observational constrain for Galactic evolutionary models that study the efficiency of stellar radial migration.
Links to: arXiv, form interface, find, astro-ph, recent, 1612, contact, help (Access key information)