Showing 1–50 of 219 results for author: Sternberg, A

GATOS: missing molecular gas in the outflow of NGC5728 revealed by JWST

Authors: R. Davies, T. Shimizu, M. Pereira-Santaella, A. Alonso-Herrero, A. Audibert, E. Bellocchi, P. Boorman, S. Campbell, Y. Cao, F. Combes, D. Delaney, T. Diaz-Santos, F. Eisenhauer, D. Esparza Arredondo, H. Feuchtgruber, N. M. Forster Schreiber, L. Fuller, P. Gandhi, I. Garcia-Bernete, S. Garcia-Burillo, B. Garcia-Lorenzo, R. Genzel, S. Gillessen, O. Gonzalez Martin, H. Haidar , et al. (27 additional authors not shown)

Abstract: The ionisation cones of NGC5728 have a deficit of molecular gas based on millimetre observations of CO(2-1) emission. Although photoionisation from the active nucleus may lead to suppression of this transition, warm molecular gas can still be present. We report the detection of eight mid-infrared rotational H$_2$ lines throughout the central kiloparsec, including the ionisation cones, using integr… ▽ More The ionisation cones of NGC5728 have a deficit of molecular gas based on millimetre observations of CO(2-1) emission. Although photoionisation from the active nucleus may lead to suppression of this transition, warm molecular gas can still be present. We report the detection of eight mid-infrared rotational H$_2$ lines throughout the central kiloparsec, including the ionisation cones, using integral field spectroscopic observations with JWST/MIRI MRS. The H$_2$ line ratios, characteristic of a power-law temperature distribution, indicate that the gas is warmest where it enters the ionisation cone through disk rotation, suggestive of shock excitation. In the nucleus, where the data can be combined with an additional seven ro-vibrational H$_2$ transitions, we find that moderate velocity (30 km s$^{-1}$) shocks in dense ($10^5$ cm$^{-3}$) gas, irradiated by an external UV field ($G_0 = 10^3$), do provide a good match to the full set. The warm molecular gas in the ionisation cone that is traced by the H$_2$ rotational lines has been heated to temperatures $>200$ K. Outside of the ionisation cone the molecular gas kinematics are undisturbed. However, within the ionisation cone, the kinematics are substantially perturbed, indicative of a radial flow, but one that is quantitatively different from the ionised lines. We argue that this outflow is in the plane of the disk, implying a short 50 pc acceleration zone up to speeds of about 400 km s$^{-1}$ followed by an extended deceleration over $\sim$700 pc where it terminates. The deceleration is due to both the radially increasing galaxy mass, and mass-loading as ambient gas in the disk is swept up. △ Less

Submitted 24 June, 2024; originally announced June 2024.

Comments: A&A accepted; 16 pages

X-ray signatures of galactic outflows into the circumgalactic medium

Authors: Ranita Jana, Kartick C. Sarkar, Jonathan Stern, Amiel Sternberg

Abstract: We present a set of controlled hydrodynamical simulations to study the effects of strong galactic outflows on the density and temperature structures, and associated X-ray signatures, of extra-planar and circumgalactic gas. We consider three initial state models, isothermal, isentropic, and rotating cooling-flow, for the hot circumgalactic medium (CGM) into which the outflows are driven. The energy… ▽ More We present a set of controlled hydrodynamical simulations to study the effects of strong galactic outflows on the density and temperature structures, and associated X-ray signatures, of extra-planar and circumgalactic gas. We consider three initial state models, isothermal, isentropic, and rotating cooling-flow, for the hot circumgalactic medium (CGM) into which the outflows are driven. The energy sources are either stellar winds and supernovae, or active galactic nuclei. We consider energy injection rates in the range $10^{40} < \dot{E}_{\rm inj} <10^{44.5}$ erg s$^{-1}$, and compute the time-dependent soft X-ray (0.5-2 keV) surface brightness. For $\dot{E}_{\rm inj} \gtrsim 10^{41} - 10^{42}$ erg s$^{-1}$, with the exact threshold depending on the initial CGM state, the X-ray response is dominated by dense hot gas in the forward shock that eventually fades into the CGM as a sound wave. The shock surrounds an inner hot bubble leading to a radial flattening of the X-ray surface brightness. For lower energy injection rates, the X-ray surface brightness of the initial CGM state is almost unaffected. We present analytic approximations for the outflow shock propagation and the associated X-ray emissions. △ Less

Submitted 9 May, 2024; originally announced May 2024.

Comments: Accepted for publication in MNRAS

Sunyaev-Zeldovich Signals from $L^*$ Galaxies: Observations, Analytics, and Simulations

Authors: Yossi Oren, Amiel Sternberg, Christopher F. McKee, Yakov Faerman, Shy Genel

Abstract: We analyze measurements of the thermal Sunyaev-Zeldovich (tSZ) effect arising in the circumgalactic medium (CGM) of $L^*$ galaxies, reported by Bregman et al. 2022 and Das et al. 2023. In our analysis we use the Faerman et al. 2017 and Faerman et al. 2020 CGM models, a new power-law model (PLM), and the TNG100 simulation. For a given $M_{\rm vir}$, our PLM has four parameters; the fraction,… ▽ More We analyze measurements of the thermal Sunyaev-Zeldovich (tSZ) effect arising in the circumgalactic medium (CGM) of $L^*$ galaxies, reported by Bregman et al. 2022 and Das et al. 2023. In our analysis we use the Faerman et al. 2017 and Faerman et al. 2020 CGM models, a new power-law model (PLM), and the TNG100 simulation. For a given $M_{\rm vir}$, our PLM has four parameters; the fraction, $f_{\rm hCGM}$, of the halo baryon mass in hot CGM gas, the ratio, $φふぁい_T$, of the actual gas temperature at the virial radius to the virial temperature, and the power-law indicies, $a_{P,{\rm th}}$ and $a_n$ for the thermal electron pressure and the hydrogen nucleon density. The B+22 Compton-$y$ profile implies steep electron pressure slopes ($a_{P,{\rm th}}\simeq 2$). For isothermal conditions the temperature is at least $1.1\times 10^6$ K, with a hot CGM gas mass of up to $3.5\times 10^{11}$ M$_\odot$ for a virial mass of $2.75\times 10^{12}$ M$_\odot$. However, if isothermal the gas must be expanding out of the halos. An isentropic equation of state is favored for which hydrostatic equilibrium is possible. The B+22 and D+23 results are consistent with each other and with recent (0.5-2 keV) CGM X-ray observations by Zhang et al. 2024 of Milky Way mass systems. For $M_{\rm vir}\simeq 3\times 10^{12}$ M$_\odot$, the scaled Compton pressure integrals, $E(z)^{-2/3}Y_{500}/M_{\rm vir,12}^{5/3}$, lie in the narrow range, $2.5\times 10^{-4}$ to $5.0\times 10^{-4}$ kpc$^2$, for all three sets of observations. TNG100 underpredicts the tSZ parameters by factors $\sim 0.5$ dex for the $L^*$ galaxies, suggesting that the feedback strengths and CGM gas losses are overestimated in the simulated halos at these mass scales. △ Less

Submitted 14 March, 2024; originally announced March 2024.

Comments: Submitted to the Astrophysical Journal. 30 pages, 18 figures

The Molecular Cloud Lifecycle II: Formation and Destruction of Molecular Clouds Diagnosed via H$_2$ Fluorescent Emission Emission

Authors: Blakesley Burkhart, Shmuel Bialy, Daniel Seifried, Stefanie Walch, Erika Hamden, Thomas J. Haworth, Keri Hoadley, Shuo Kong, Madisen Johnson, Sarah Jeffreson, Mark R. Krumholz, Min-Young Lee, Amiel Sternberg, Neal J. Turner

Abstract: Molecular hydrogen (H$_2$) formation and dissociation are key processes that drive the gas lifecycle in galaxies. Using the SImulating the LifeCycle of Molecular Clouds (SILCC) zoom-in simulation suite, we explore the utility of future observations of H$_2$ dissociation and formation for tracking the lifecycle of molecular clouds. The simulations used in this work include non-equilibrium H$_2$ for… ▽ More Molecular hydrogen (H$_2$) formation and dissociation are key processes that drive the gas lifecycle in galaxies. Using the SImulating the LifeCycle of Molecular Clouds (SILCC) zoom-in simulation suite, we explore the utility of future observations of H$_2$ dissociation and formation for tracking the lifecycle of molecular clouds. The simulations used in this work include non-equilibrium H$_2$ formation, stellar radiation, sink particles, and turbulence. We find that, at early times in the cloud evolution, H$_2$ formation rapidly outpaces dissociation and molecular clouds build their mass from the atomic reservoir in their environment. Rapid H$_2$ formation is also associated with a higher early star formation rate. For the clouds studied here, H$_2$ is strongly out of chemical equilibrium during the early stages of cloud formation but settles into a bursty chemical steady-state about 2 Myrs after the first stars form. At the latest stage of cloud evolution, dissociation outweighs formation and the clouds enter a dispersal phase. We discuss how theories for the molecular cloud lifecycle and the star formation efficiency may be distinguished with observational measurements of H$_2$ fluorescence with a space-based high-resolution FUV spectrometer, such as the proposed Hyperion and Eos NASA Explorer missions. Such missions would enable measurements of the H$_2$ dissociation and formation rates, which we demonstrate can be connected to different phases in a molecular cloud's star-forming life, including cloud building, rapidly star-forming, H$_2$ chemical equilibrium, and cloud destruction. △ Less

Submitted 5 February, 2024; v1 submitted 2 February, 2024; originally announced February 2024.

Comments: Submitted to ApJ, comments welcome

A dynamical measure of the black hole mass in a quasar 11 billion years ago

Authors: R. Abuter, F. Allouche, A. Amorim, C. Bailet, A. Berdeu, J. -P. Berger, P. Berio, A. Bigioli, O. Boebion, M. -L. Bolzer, H. Bonnet, G. Bourdarot, P. Bourget, W. Brandner, Y. Cao, R. Conzelmann, M. Comin, Y. Clénet, B. Courtney-Barrer, R. Davies, D. Defrère, A. Delboulbé, F. Delplancke-Ströbele, R. Dembet, J. Dexter , et al. (102 additional authors not shown)

Abstract: Tight relationships exist in the local universe between the central stellar properties of galaxies and the mass of their supermassive black hole. These suggest galaxies and black holes co-evolve, with the main regulation mechanism being energetic feedback from accretion onto the black hole during its quasar phase. A crucial question is how the relationship between black holes and galaxies evolves… ▽ More Tight relationships exist in the local universe between the central stellar properties of galaxies and the mass of their supermassive black hole. These suggest galaxies and black holes co-evolve, with the main regulation mechanism being energetic feedback from accretion onto the black hole during its quasar phase. A crucial question is how the relationship between black holes and galaxies evolves with time; a key epoch to probe this relationship is at the peaks of star formation and black hole growth 8-12 billion years ago (redshifts 1-3). Here we report a dynamical measurement of the mass of the black hole in a luminous quasar at a redshift of 2, with a look back time of 11 billion years, by spatially resolving the broad line region. We detect a 40 micro-arcsecond (0.31 pc) spatial offset between the red and blue photocenters of the H$αあるふぁ$ line that traces the velocity gradient of a rotating broad line region. The flux and differential phase spectra are well reproduced by a thick, moderately inclined disk of gas clouds within the sphere of influence of a central black hole with a mass of 3.2x10$^{8}$ solar masses. Molecular gas data reveal a dynamical mass for the host galaxy of 6x10$^{11}$ solar masses, which indicates an under-massive black hole accreting at a super-Eddington rate. This suggests a host galaxy that grew faster than the supermassive black hole, indicating a delay between galaxy and black hole formation for some systems. △ Less

Submitted 25 January, 2024; originally announced January 2024.

Comments: 5 pages Main text, 8 figures, 2 tables, to be published in Nature, under embargo until 29 January 2024 16:00 (London)

The size-luminosity relation of local active galactic nuclei from interferometric observations of the broad-line region

Authors: GRAVITY Collaboration, A. Amorim, G. Bourdarot, W. Brandner, Y. Cao, Y. Clénet, R. Davies, P. T. de Zeeuw, J. Dexter, A. Drescher, A. Eckart, F. Eisenhauer, M. Fabricius, H. Feuchtgruber, N. M. Förster Schreiber, P. J. V. Garcia, R. Genzel, S. Gillessen, D. Gratadour, S. Hönig, M. Kishimoto, S. Lacour, D. Lutz, F. Millour, H. Netzer , et al. (20 additional authors not shown)

Abstract: By using the GRAVITY instrument with the near-infrared (NIR) Very Large Telescope Interferometer (VLTI), the structure of the broad (emission-)line region (BLR) in active galactic nuclei (AGNs) can be spatially resolved, allowing the central black hole (BH) mass to be determined. This work reports new NIR VLTI/GRAVITY interferometric spectra for four type 1 AGNs (Mrk 509, PDS 456, Mrk 1239, and IC… ▽ More By using the GRAVITY instrument with the near-infrared (NIR) Very Large Telescope Interferometer (VLTI), the structure of the broad (emission-)line region (BLR) in active galactic nuclei (AGNs) can be spatially resolved, allowing the central black hole (BH) mass to be determined. This work reports new NIR VLTI/GRAVITY interferometric spectra for four type 1 AGNs (Mrk 509, PDS 456, Mrk 1239, and IC 4329A) with resolved broad-line emission. Dynamical modelling of interferometric data constrains the BLR radius and central BH mass measurements for our targets and reveals outflow-dominated BLRs for Mrk 509 and PDS 456. We present an updated radius-luminosity (R-L) relation independent of that derived with reverberation mapping (RM) measurements using all the GRAVITY-observed AGNs. We find our R-L relation to be largely consistent with that derived from RM measurements except at high luminosity, where BLR radii seem to be smaller than predicted. This is consistent with RM-based claims that high Eddington ratio AGNs show consistently smaller BLR sizes. The BH masses of our targets are also consistent with the standard $M_\mathrm{BH}$-$σしぐま_*$ relation. Model-independent photocentre fitting shows spatial offsets between the hot dust continuum and the BLR photocentres (ranging from $\sim$17 $μみゅー$as to 140 $μみゅー$as) that are generally perpendicular to the alignment of the red- and blueshifted BLR photocentres. These offsets are found to be related to the AGN luminosity and could be caused by asymmetric K-band emission of the hot dust, shifting the dust photocentre. We discuss various possible scenarios that can explain this phenomenon. △ Less

Submitted 15 January, 2024; originally announced January 2024.

Comments: Accepted for publication at A&A

HI in Molecular Clouds: Irradiation by FUV plus Cosmic Rays

Authors: Amiel Sternberg, Shmuel Bialy, Alon Gurman

Abstract: We extend the analytic theory presented by Sternberg et al. (2014) and Bialy & Sternberg (2016) for the production of atomic hydrogen (HI) via FUV photodissociation at the boundaries of dense interstellar molecular (H$_2$) clouds, to also include the effects of penetrating (low-energy) cosmic-rays for the growth of the total HI column densities. We compute the steady-state abundances of the HI and… ▽ More We extend the analytic theory presented by Sternberg et al. (2014) and Bialy & Sternberg (2016) for the production of atomic hydrogen (HI) via FUV photodissociation at the boundaries of dense interstellar molecular (H$_2$) clouds, to also include the effects of penetrating (low-energy) cosmic-rays for the growth of the total HI column densities. We compute the steady-state abundances of the HI and H$_2$ in one-dimensional gas slabs in which the FUV photodissociation rates are reduced by depth-dependent H$_2$ self-shielding and dust absorption, and in which the cosmic-ray ionization rates are either constant or reduced by transport effects. The solutions for the HI and H$_2$ density profiles and the integrated HI columns, depend primarily on the ratios $I_{\rm UV}/Rn$ and $ζぜーた/Rn$, where $I_{\rm UV}$ is the intensity of the photodissociating FUV field, $ζぜーた$ is the H$_2$ cosmic-ray ionization rate, $n$ is the hydrogen gas density, and $R$ is the dust-surface H$_2$ formation rate coefficient. We present computations for a wide range of FUV field strengths, cosmic-ray ionization rates, and dust-to-gas ratios. We develop analytic expressions for the growth of the HI column densities. For Galactic giant molecular clouds (GMCs) with multiphased (warm/cold) HI envelopes, the interior cosmic-ray zones will dominate the production of the HI only if $ζぜーた\gtrsim 4.5\times 10^{-16} \times (M_{\rm GMC}/10^6 \ M_{\odot})^{-1/2}$~s$^{-1}$, where $M_{\rm GMC}$ is the GMC mass, and including attenuation of the cosmic-ray fluxes. For most Galactic GMCs and conditions, FUV photodissociation dominates over cosmic-ray ionization for the production of the HI column densities. Furthermore, the cosmic-rays do not affect the HI-to-H$_2$ transition points. △ Less

Submitted 16 October, 2023; v1 submitted 26 August, 2023; originally announced August 2023.

Comments: Accepted for publication in ApJ

[CII] Emission in a Self-Regulated Interstellar Medium

Authors: Alon Gurman, Chia-Yu Hu, Amiel Sternberg, Ewine F. van Dishoeck

Abstract: The [CII] 157.74 $μみゅー$m fine structure transition is one of the brightest and most well-studied emission lines in the far-infrared (FIR), produced in the interstellar medium (ISM) of galaxies. We study its properties in sub-pc resolution hydrodynamical simulations for an ISM patch with gas surface density of $Σしぐま_{\rm{g}}=10\;M_{\odot}\;\rm{pc}^{-2}$, coupled with time-dependent chemistry, far-ultravi… ▽ More The [CII] 157.74 $μみゅー$m fine structure transition is one of the brightest and most well-studied emission lines in the far-infrared (FIR), produced in the interstellar medium (ISM) of galaxies. We study its properties in sub-pc resolution hydrodynamical simulations for an ISM patch with gas surface density of $Σしぐま_{\rm{g}}=10\;M_{\odot}\;\rm{pc}^{-2}$, coupled with time-dependent chemistry, far-ultraviolet (FUV) dust and gas shielding, star formation, photoionization and supernova (SN) feedback, and full line-radiative transfer. We find a [CII]-to-H$_2$ conversion factor that scales weakly with metallicity $X_{\rm{[CII]}}=6.31\times 10^{19} \;Z^{\prime\;0.17}\; \rm{cm}^{-2}\;(\rm{K}\;\rm{km}\;\rm{s}^{-1})^{-1}$, where $Z^{\prime}$ is the normalized metallicity relative to solar. {The majority of [CII] originates from atomic gas with hydrogen number density $n\sim 10~{\rm cm^{-3}}$.} The [CII] line intensity positively correlates with the star formation rate (SFR), with a normalization factor that scales linearly with metallicity. We find that this is broadly consistent with $z\sim0$ observations. As such, [CII] is a good SFR tracer even in metal-poor environments where molecular lines might be undetectable. Resolving the clumpy structure of the dense ($n=10-10^3\;\rm{cm}^{-3}$) interstellar medium (ISM) is important as it dominates [CII] 157.74 $μみゅー$m emission. We compare our full radiative transfer computation with the optically-thin limit and find that the [CII] line becomes marginally optically thick only at super-solar metallicity for our assumed gas surface density. △ Less

Submitted 3 March, 2024; v1 submitted 11 August, 2023; originally announced August 2023.

Comments: Accepted for publication in ApJ

Evidence for Large-scale, Rapid Gas Inflows in z~2 Star-forming Disks

Authors: R. Genzel, J. -B. Jolly, D. Liu, S. H. Price, L. L. Lee, N. M. Förster Schreiber, L. J. Tacconi, R. Herrera-Camus, C. Barfety, A. Burkert, Y. Cao, R. I. Davies, A. Dekel, M. M. Lee, D. Lutz, T. Naab, R. Neri, A. Nestor Shachar, S. Pastras, C. Pulsoni, A. Renzini, K. Schuster, T. T. Shimizu, F. Stanley, A. Sternberg , et al. (1 additional authors not shown)

Abstract: We report high-quality H$αあるふぁ$/CO, imaging spectroscopy of nine massive (log median stellar mass = 10.65 $M_{\odot}$), disk galaxies on the star-forming, main sequence (henceforth `SFGs'), near the peak of cosmic galaxy evolution ($z\sim$1.1-2.5), taken with the ESO-Very Large Telescope, IRAM-NOEMA and Atacama Large Millimeter/submillimeter Array. We fit the major axis position-velocity cuts with bea… ▽ More We report high-quality H$αあるふぁ$/CO, imaging spectroscopy of nine massive (log median stellar mass = 10.65 $M_{\odot}$), disk galaxies on the star-forming, main sequence (henceforth `SFGs'), near the peak of cosmic galaxy evolution ($z\sim$1.1-2.5), taken with the ESO-Very Large Telescope, IRAM-NOEMA and Atacama Large Millimeter/submillimeter Array. We fit the major axis position-velocity cuts with beam-convolved, forward models with a bulge, a turbulent rotating disk, and a dark matter (DM) halo. We include priors for stellar and molecular gas masses, optical light effective radii and inclinations, and DM masses from our previous rotation curve analyses of these galaxies. We then subtract the inferred 2D model-galaxy velocity and velocity dispersion maps from those of the observed galaxies. We investigate whether the residual velocity and velocity dispersion maps show indications for radial flows. We also carry out kinemetry, a model-independent tool for detecting radial flows. We find that all nine galaxies exhibit significant non-tangential flows. In six SFGs, the inflow velocities ($v_r\sim$30-90 km s$^{-1}$, 10%-30% of the rotational component) are along the minor axis of these galaxies. In two cases the inflow appears to be off the minor axis. The magnitudes of the radial motions are in broad agreement with the expectations from analytic models of gravitationally unstable, gas-rich disks. Gravitational torques due to clump and bar formation, or spiral arms, drive gas rapidly inward and result in the formation of central disks and large bulges. If this interpretation is correct, our observations imply that gas is transported into the central regions on ~10 dynamical time scales. △ Less

Submitted 27 September, 2023; v1 submitted 4 May, 2023; originally announced May 2023.

Comments: 31 pages, 27 figures, 3 tables. The Astrophysical Journal in press

Co-evolution of Dust and Chemistry in Galaxy Simulations with a Resolved Interstellar Medium

Authors: Chia-Yu Hu, Amiel Sternberg, Ewine F. van Dishoeck

Abstract: Nearby dwarf irregular galaxies are ideal laboratories for studying the interstellar medium (ISM) at low metallicity, which is expected to be common for galaxies at very high redshift being observed by the James Webb Space Telescope. We present the first high-resolution ($\sim 0.2$~pc) hydrodynamical simulations of an isolated low-metallicity ($0.1~Z_\odot$) dwarf galaxy coupled with a time-depend… ▽ More Nearby dwarf irregular galaxies are ideal laboratories for studying the interstellar medium (ISM) at low metallicity, which is expected to be common for galaxies at very high redshift being observed by the James Webb Space Telescope. We present the first high-resolution ($\sim 0.2$~pc) hydrodynamical simulations of an isolated low-metallicity ($0.1~Z_\odot$) dwarf galaxy coupled with a time-dependent chemistry network and a dust evolution model where dust is locally produced and destroyed by various processes. To accurately model carbon monoxide (CO), we post-process the simulations with a detailed chemistry network including the time-dependent effect of molecular hydrogen (H$_2$). Our model successfully reproduces the observed star formation rate and CO(1-0) luminosity ($L_{\rm CO}$). We find that dust growth in dense gas is required to reproduce the observed $L_{\rm CO}$ as otherwise CO would be completely photodissociated. In contrast, the H$_2$ abundance is extremely small and is insensitive to dust growth, leading to a CO-to-H$_2$ conversion factor that is only slightly higher than the Milky Way value despite the low metallicity. An observationally inferred dust-to-gas ratio is thus underestimated if adopting the metallicity-dependent CO-to-H$_2$ conversion factor. The newly-produced dust in dense gas mixes with the ISM through supernova feedback without being completely destroyed by sputtering, which leads to galactic outflows 20% - 50% dustier than the ISM, providing a possible source for intergalactic dust. △ Less

Submitted 10 July, 2023; v1 submitted 12 January, 2023; originally announced January 2023.

Comments: ApJ accepted. 19 pages, 10 figures

Hyperion: The origin of the stars A far-UV space telescope for high-resolution spectroscopy over wide fields

Authors: Erika Hamden, David Schiminovich, Shouleh Nikzad, Neal J. Turner, Blakesley Burkhart, Thomas J. Haworth, Keri Hoadley, Jinyoung Serena Kim, Shmuel Bialyh, Geoff Bryden, Haeun Chung, Nia Imara, Rob Kennicutt, Jorge Pineda, Shuo Konga, Yasuhiro Hasegawa, Ilaria Pascucci, Benjamin Godard, Mark Krumholz, Min-Young Lee, Daniel Seifried, Amiel Sternberg, Stefanie Walch, Miles Smith, Stephen C. Unwin , et al. (8 additional authors not shown)

Abstract: We present Hyperion, a mission concept recently proposed to the December 2021 NASA Medium Explorer announcement of opportunity. Hyperion explores the formation and destruction of molecular clouds and planet-forming disks in nearby star-forming regions of the Milky Way. It does this using long-slit, high-resolution spectroscopy of emission from fluorescing molecular hydrogen, which is a powerful fa… ▽ More We present Hyperion, a mission concept recently proposed to the December 2021 NASA Medium Explorer announcement of opportunity. Hyperion explores the formation and destruction of molecular clouds and planet-forming disks in nearby star-forming regions of the Milky Way. It does this using long-slit, high-resolution spectroscopy of emission from fluorescing molecular hydrogen, which is a powerful far-ultraviolet (FUV) diagnostic. Molecular hydrogen (H2) is the most abundant molecule in the universe and a key ingredient for star and planet formation, but is typically not observed directly because its symmetric atomic structure and lack of a dipole moment mean there are no spectral lines at visible wavelengths and few in the infrared. Hyperion uses molecular hydrogen's wealth of FUV emission lines to achieve three science objectives: (1) determining how star formation is related to molecular hydrogen formation and destruction at the boundaries of molecular clouds; (2) determining how quickly and by what process massive star feedback disperses molecular clouds; and (3) determining the mechanism driving the evolution of planet-forming disks around young solar-analog stars. Hyperion conducts this science using a straightforward, highly-efficient, single-channel instrument design. Hyperion's instrument consists of a 48 cm primary mirror, with an f/5 focal ratio. The spectrometer has two modes, both covering 138.5-161.5 nm bandpasses. A low resolution mode has a spectral resolution of R>10,000 with a slit length of 65 arcmin, while the high resolution mode has a spectral resolution of R>50,000 over a slit length of 5 armin. Hyperion occupies a 2 week long, high-earth, Lunar resonance TESS-like orbit, and conducts 2 weeks of planned observations per orbit, with time for downlinks and calibrations. Hyperion was reviewed as Category I, which is the highest rating possible, but was not selected. △ Less

Submitted 13 December, 2022; originally announced December 2022.

Comments: Accepted to JATIS, 9 Figures

Towards measuring supermassive black hole masses with interferometric observations of the dust continuum

Authors: GRAVITY Collaboration, A. Amorim, G. Bourdarot, W. Brandner, Y. Cao, Y. Clénet, R. Davies, P. T. de Zeeuw, J. Dexter, A. Drescher, A. Eckart, F. Eisenhauer, M. Fabricius, N. M. Förster Schreiber, P. J. V. Garcia, R. Genzel, S. Gillessen, D. Gratadour, S. Hönig, M. Kishimoto, S. Lacour, D. Lutz, F. Millour, H. Netzer, T. Ott , et al. (18 additional authors not shown)

Abstract: This work focuses on active galactic nuclei (AGNs), and the relation between the sizes of the hot dust continuum and the broad-line region (BLR). We find that the continuum size measured using optical/near-infrared interferometry (OI) is roughly twice that measured by reverberation mapping (RM). Both OI and RM continuum sizes show a tight relation with the H$βべーた$ BLR size with only an intrinsic scat… ▽ More This work focuses on active galactic nuclei (AGNs), and the relation between the sizes of the hot dust continuum and the broad-line region (BLR). We find that the continuum size measured using optical/near-infrared interferometry (OI) is roughly twice that measured by reverberation mapping (RM). Both OI and RM continuum sizes show a tight relation with the H$βべーた$ BLR size with only an intrinsic scatter of 0.25 dex. The masses of supermassive black holes (BHs) can hence be simply derived from a dust size in combination with a broad line width and virial factor. Since the primary uncertainty of these BH masses comes from the virial factor, the accuracy of the continuum-based BH masses is close to those based on the RM measurement of the broad emission line. Moreover, the necessary continuum measurements can be obtained on a much shorter timescale than those required monitoring for RM, and are also more time efficient than those needed to resolve the BLR with OI. The primary goal of this work is to demonstrate measuring the BH mass based on the dust continuum size with our first calibration of the $R_\mathrm{BLR}$-$R_\mathrm{d}$ relation. The current limitation and caveats are discussed in detail. Future GRAVITY observations are expected to improve the continuum-based method and have the potential to measure BH masses for a large sample of AGNs in the low-redshift Universe. △ Less

Submitted 28 September, 2022; originally announced September 2022.

Comments: Accepted for publication in A&A; 11 pages, 5 figures, 3 tables

Journal ref: A&A 669, A14 (2023)

RC100: Rotation Curves of 100 Massive Star-Forming Galaxies at z=0.6-2.5 Reveal Little Dark Matter on Galactic Scales

Authors: A. Nestor Shachar, S. H. Price, N. M. Förster Schreiber, R. Genzel, T. T. Shimizu, L. J. Tacconi, H. Übler, A. Burkert, R. I. Davies, A. Deke, R. Herrera-Camus, L. L. Lee, D. Liu, D. Lutz, T. Naab, R. Neri, A. Renzini, R. Saglia, K. Schuster, A. Sternberg, E. Wisnioski, S. Wuyts

Abstract: We analyze Ha or CO rotation curves (RCs) extending out to several galaxy effective radii for 100 massive, large, star-forming disk galaxies (SFGs) across the peak of cosmic galaxy star formation (z~0.6-2.5), more than doubling the previous sample presented by Genzel et al. (2020) and Price et al. (2021). The observations were taken with SINFONI and KMOS integral-field spectrographs at ESO-VLT, LU… ▽ More We analyze Ha or CO rotation curves (RCs) extending out to several galaxy effective radii for 100 massive, large, star-forming disk galaxies (SFGs) across the peak of cosmic galaxy star formation (z~0.6-2.5), more than doubling the previous sample presented by Genzel et al. (2020) and Price et al. (2021). The observations were taken with SINFONI and KMOS integral-field spectrographs at ESO-VLT, LUCI at LBT, NOEMA at IRAM, and ALMA. We fit the major axis kinematics with beam-convolved, forward models of turbulent rotating disks with bulges embedded in dark matter (DM) halos, including the effects of pressure support. The fraction of dark to total matter within the disk effective radius ($R_e ~ 5 kpc$), $f_DM (R_e)=V_{DM}^2 (R_e)/V_{circ}^2 (R_e)$, decreases with redshift: At z~1 (z~2) the median DM fraction is $0.38\pm 0.23$ ($0.27\pm 0.18$), and a third (half) of all galaxies are "maximal" disks with $f_{DM} (R_e)<0.28$. Dark matter fractions correlate inversely with the baryonic surface density, and the low DM fractions require a flattened, or cored, inner DM density distribution. At z~2 there is ~40% less dark matter mass on average within $R_e$ compared to expected values based on cosmological stellar-mass halo-mass relations. The DM deficit is more evident at high star formation rate (SFR) surface densities ($Σしぐま_{SFR}>2.5 M_{\odot} yr^{-1} kpc^{-2}$) and galaxies with massive bulges ($M_{bulge}>10^{10} M_{\odot}$). A combination of stellar or active galactic nucleus (AGN) feedback, and/or heating due to dynamical friction, either from satellite accretion or clump migration, may drive the DM from cuspy into cored mass distributions. The observations plausibly indicate an efficient build-up of massive bulges and central black holes at z~2 SFGs. △ Less

Submitted 25 September, 2022; originally announced September 2022.

Comments: Submitted to ApJ (34 pages, 12 figures, 4 tables)

Self-ionizing galactic winds

Authors: Kartick C. Sarkar, Amiel Sternberg, Orly Gnat

Abstract: We present hydrodynamical simulations of galactic winds from star-forming galaxies including non-equilibrium ionization and frequency-dependent radiative transfer, processes that have remained largely unaccounted for in galactic wind studies. We consider radiation from massive stars, the metagalactic UV/X-ray background, and the self-radiation of the supernovae heated gas. We compare our results t… ▽ More We present hydrodynamical simulations of galactic winds from star-forming galaxies including non-equilibrium ionization and frequency-dependent radiative transfer, processes that have remained largely unaccounted for in galactic wind studies. We consider radiation from massive stars, the metagalactic UV/X-ray background, and the self-radiation of the supernovae heated gas. We compare our results to classical galactic wind solutions and show the importance of our newly included physical processes toward observations of ions such as O III, O VI, O VII and O VIII plus the observable soft X-ray spectra. Non-equilibrium ionization is reflected in over-ionized gas compared to equilibrium solutions, leading to much enhanced column densities of highly ionized species. The wind produces excess soft X-ray ($E\gtrsim 100$ eV) radiation that is several orders of magnitude higher compared to the metagalactic background. This radiation ionizes the higher ions (such as O VII) somewhat, but affects the lower ions (such as O III) significantly. We predict that the observable X-ray spectra should contain the signatures of such non-equilibrium effects, especially in X-ray lines such as O VII and O VIII. Simple estimates suggest that both the temperature and density of the winds may be overestimated by factors of a few to almost 2 orders of magnitude using simple equilibrium models. We conclude that both the non-equilibrium ionization and the radiation from the wind itself need to be considered for proper modeling of the optical/UV/X-ray emitting plasma in galactic winds. △ Less

Submitted 11 October, 2022; v1 submitted 29 March, 2022; originally announced March 2022.

Comments: 19 pages, 11 figures. Accepted for publication in ApJ. Added new estimates and figures for ISM absorption of the wind radiation

The Life and Times of Giant Molecular Clouds

Authors: Mélanie Chevance, Mark R. Krumholz, Anna F. McLeod, Eve C. Ostriker, Erik W. Rosolowsky, Amiel Sternberg

Abstract: Giant molecular clouds (GMCs) are the sites of star formation and stellar feedback in galaxies. Their properties set the initial conditions for star formation and their lifecycles determine how feedback regulates galaxy evolution. In recent years, the advent of high-resolution telescopes has enabled systematic GMC-scale studies of the molecular interstellar medium in nearby galaxies, now covering… ▽ More Giant molecular clouds (GMCs) are the sites of star formation and stellar feedback in galaxies. Their properties set the initial conditions for star formation and their lifecycles determine how feedback regulates galaxy evolution. In recent years, the advent of high-resolution telescopes has enabled systematic GMC-scale studies of the molecular interstellar medium in nearby galaxies, now covering a wide range of physical conditions and allowing for the first studies of how GMC properties depend on galactic environment. These observational developments have been accompanied by numerical simulations of improving resolution that are increasingly accurately accounting for the effects of the galactic-scale environment on GMCs, while simultaneously improving the treatment of the small-scale processes of star-formation and stellar feedback within them. The combination of these recent developments has greatly improved our understanding of the formation, evolution, and destruction of GMCs. We review the current state of the field, highlight current open questions, and discuss promising avenues for future studies. △ Less

Submitted 17 March, 2022; originally announced March 2022.

Comments: 37 pages, 5 figures; To appear in Protostars and Planets VII; Editors: Shu-ichiro Inutsuka, Yuri Aikawa, Takayuki Muto, Kengo Tomida, and Motohide Tamura

Kiloparsec view of a typical star-forming galaxy when the Universe was $\sim$1 Gyr old II. Regular rotating disk and evidence for baryon dominance on galactic scales

Authors: R. Herrera-Camus, N. M. Förster Schreiber, S. H. Price, H. Übler, A. D. Bolatto, R. L. Davies, D. Fisher, R. Genzel, D. Lutz, T. Naab, A. Nestor, T. Shimizu, A. Sternberg, L. Tacconi, K. Tadaki

Abstract: We present a kinematic analysis of the main-sequence galaxy HZ4 at $z=5.5$. Our study is based on deep, spatially resolved observations of the [CII] 158 $μみゅー$m transition obtained with the Atacama Large Millimeter/Submillimeter Array (ALMA). From the combined analysis of the disk morphology, the two-dimensional velocity structure, and forward-modeling of the one-dimensional velocity and velocity dis… ▽ More We present a kinematic analysis of the main-sequence galaxy HZ4 at $z=5.5$. Our study is based on deep, spatially resolved observations of the [CII] 158 $μみゅー$m transition obtained with the Atacama Large Millimeter/Submillimeter Array (ALMA). From the combined analysis of the disk morphology, the two-dimensional velocity structure, and forward-modeling of the one-dimensional velocity and velocity dispersion profiles, we conclude that HZ4 has a regular rotating disk in place. The intrinsic velocity dispersion in HZ4 is high ($σしぐま_{0}=65.8^{+2.9}_{-3.3}$ km s$^{-1}$), and the ratio between the rotational velocity and the intrinsic velocity dispersion is $V_{\rm rot}/σしぐま_{0}=2.2$. These values are consistent with the expectations from the trends of increasing $σしぐま_{0}$ and decreasing $V_{\rm rot}/σしぐま_{0}$ as a function of redshift observed in main-sequence galaxies up to $z\approx4$. Galaxy evolution models suggest that the high level of turbulence observed in HZ4 can only be achieved if, in addition to stellar feedback, there is radial transport of gas within the disk. Finally, we find that HZ4 is baryon dominated on galactic scales ($\lesssim2\times R_{\rm e}$), with a dark matter fraction at one effective radius of $f_{\rm DM}(R_{\rm e})=0.41^{+0.25}_{-0.22}$. This value is comparable to the dark matter fractions found in lower redshift galaxies that could be the descendants of HZ4: massive ($M_{\star}\approx10^{11}~M_{\odot}$), star-forming galaxies at $z\sim2$, and passive, early type galaxies at $z\approx0$. △ Less

Submitted 1 March, 2022; originally announced March 2022.

Comments: Accepted for publication in A&A Letters, 10 pages, 6 figures

Journal ref: A&A 665, L8 (2022)

Dependence of $X_{\rm CO}$ on metallicity, intensity, and spatial scale in a self-regulated interstellar medium

Authors: Chia-Yu Hu, Andreas Schruba, Amiel Sternberg, Ewine F. van Dishoeck

Abstract: We study the CO(1-0)-to-H$_2$ conversion factor ($X_{\rm CO}$) and the line ratio of CO(2-1)-to-CO(1-0) ($R_{21}$) across a wide range of metallicity ($0.1 \leq Z/Z_\odot \leq 3$) in high-resolution (~0.2 pc) hydrodynamical simulations of a self-regulated multiphase interstellar medium. We construct synthetic CO emission maps via radiative transfer and systematically vary the "observational" beam… ▽ More We study the CO(1-0)-to-H$_2$ conversion factor ($X_{\rm CO}$) and the line ratio of CO(2-1)-to-CO(1-0) ($R_{21}$) across a wide range of metallicity ($0.1 \leq Z/Z_\odot \leq 3$) in high-resolution (~0.2 pc) hydrodynamical simulations of a self-regulated multiphase interstellar medium. We construct synthetic CO emission maps via radiative transfer and systematically vary the "observational" beam size to quantify the scale dependence. We find that the kpc-scale $X_{\rm CO}$ can be over-estimated at low $Z$ if assuming steady-state chemistry or assuming that the star-forming gas is H$_2$-dominated. On parsec scales, $X_{\rm CO}$ varies by orders of magnitude from place to place, primarily driven by the transition from atomic carbon to CO. The pc-scale $X_{\rm CO}$ drops to the Milky Way value of $2\times 10^{20}\ {\rm cm^{-2}~(K~km~s^{-1})^{-1}}$ once dust shielding becomes effective, independent of $Z$. The CO lines become increasingly optically thin at lower $Z$, leading to a higher $R_{21}$. Most cloud area is filled by diffuse gas with high $X_{\rm CO}$ and low $R_{21}$, while most CO emission originates from dense gas with low $X_{\rm CO}$ and high $R_{21}$. Adopting a constant $X_{\rm CO}$ strongly over- (under-)estimates H$_2$ in dense (diffuse) gas. The line intensity negatively (positively) correlates with $X_{\rm CO}$ ($R_{21}$) as it is a proxy of column density (volume density). On large scales, $X_{\rm CO}$ and $R_{21}$ are dictated by beam averaging, and they are naturally biased towards values in dense gas. Our predicted $X_{\rm CO}$ is a multivariate function of $Z$, line intensity, and beam size, which can be used to more accurately infer the H$_2$ mass. △ Less

Submitted 23 May, 2022; v1 submitted 11 January, 2022; originally announced January 2022.

Comments: Published in ApJ (typos in the abstract fixed). Code for interpolating Lagrangian (particle) data onto an adaptive mesh (and auxiliary scripts for RADMC-3D) available at https://github.com/huchiayu/ParticleGridMapper.jl

Rotation Curves in z~1-2 Star-Forming Disks: Comparison of Dark Matter Fractions and Disk Properties for Different Fitting Methods

Authors: S. H. Price, T. T. Shimizu, R. Genzel, H. Übler, N. M. Förster Schreiber, L. J. Tacconi, R. I. Davies, R. T. Coogan, D. Lutz, S. Wuyts, E. Wisnioski, A. Nestor, A. Sternberg, A. Burkert, R. Bender, A. Contursi, R. L. Davies, R. Herrera-Camus, M. -J. Lee, T. Naab, R. Neri, A. Renzini, R. Saglia, A. Schruba, K. Schuster

Abstract: We present a follow-up analysis examining the dynamics and structures of 41 massive, large star-forming galaxies at z~0.67-2.45 using both ionized and molecular gas kinematics. We fit the galaxy dynamics with models consisting of a bulge, a thick, turbulent disk, and a NFW dark matter halo, using code that fully forward models the kinematics, including all observational and instrumental effects. W… ▽ More We present a follow-up analysis examining the dynamics and structures of 41 massive, large star-forming galaxies at z~0.67-2.45 using both ionized and molecular gas kinematics. We fit the galaxy dynamics with models consisting of a bulge, a thick, turbulent disk, and a NFW dark matter halo, using code that fully forward models the kinematics, including all observational and instrumental effects. We explore the parameter space using Markov Chain Monte Carlo (MCMC) sampling, including priors based on stellar and gas masses and disk sizes. We fit the full sample using extracted 1D kinematic profiles. For a subset of 14 well-resolved galaxies, we also fit the 2D kinematics. The MCMC approach robustly confirms the results from least-squares fitting presented in Paper I (Genzel et al. 2020): the sample galaxies tend to be baryon-rich on galactic scales (within one effective radius). The 1D and 2D MCMC results are also in good agreement for the subset, demonstrating that much of the galaxy dynamical information is captured along the major axis. The 2D kinematics are more affected by the presence of non-circular motions, which we illustrate by constructing a toy model with constant inflow for one galaxy that exhibits residual signatures consistent with radial motions. This analysis, together with results from Paper I and other studies, strengthens the finding that massive, star-forming galaxies at z~1-2 are baryon-dominated on galactic scales, with lower dark matter fractions towards higher baryonic surface densities. Finally, we present details of the kinematic fitting code used in this analysis. △ Less

Submitted 6 September, 2021; originally announced September 2021.

Comments: Accepted for publication in ApJ (23 pages, 8 figures, 5 tables)

A geometric distance to the supermassive black Hole of NGC 3783

Authors: GRAVITY Collaboration, A. Amorim, M. Bauböck, M. C. Bentz, W. Brandner, M. Bolzer, Y. Clénet, R. Davies, P. T. de Zeeuw, J. Dexter, A. Drescher, A. Eckart, F. Eisenhauer, N. M. Förster Schreiber, P. J. V. Garcia, R. Genzel, S. Gillessen, D. Gratadour, S. Hönig, D. Kaltenbrunner, M. Kishimoto, S. Lacour, D. Lutz, F. Millour, H. Netzer , et al. (23 additional authors not shown)

Abstract: The angular size of the broad line region (BLR) of the nearby active galactic nucleus (AGN) NGC 3783 has been spatially resolved by recent observations with VLTI/GRAVITY. A reverberation mapping (RM) campaign has also recently obtained high quality light curves and measured the linear size of the BLR in a way that is complementary to the GRAVITY measurement. The size and kinematics of the BLR can… ▽ More The angular size of the broad line region (BLR) of the nearby active galactic nucleus (AGN) NGC 3783 has been spatially resolved by recent observations with VLTI/GRAVITY. A reverberation mapping (RM) campaign has also recently obtained high quality light curves and measured the linear size of the BLR in a way that is complementary to the GRAVITY measurement. The size and kinematics of the BLR can be better constrained by a joint analysis that combines both GRAVITY and RM data. This, in turn, allows us to obtain the mass of the supermassive black hole in NGC3783 with an accuracy that is about a factor of two better than that inferred from GRAVITY data alone. We derive $M_\mathrm{BH}=2.54_{-0.72}^{+0.90}\times 10^7\,M_\odot$. Finally, and perhaps most notably, we are able to measure a geometric distance to NGC 3783 of $39.9^{+14.5}_{-11.9}$ Mpc. We are able to test the robustness of the BLR-based geometric distance with measurements based on the Tully-Fisher relation and other indirect methods. We find the geometric distance is consistent with other methods within their scatter. We explore the potential of BLR-based geometric distances to directly constrain the Hubble constant, $H_0$, and identify differential phase uncertainties as the current dominant limitation to the $H_0$ measurement precision for individual sources. △ Less

Submitted 29 July, 2021; originally announced July 2021.

Comments: 9 pages and 5 figures in main text, Accepted for publication in A&A

Journal ref: A&A 654, A85 (2021)

Exploring the Milky Way Circumgalactic Medium in a Cosmological Context with a Semi-Analytic Model

Authors: Yakov Faerman, Viraj Pandya, Rachel S. Somerville, Amiel Sternberg

Abstract: We combine the Santa-Cruz Semi-Analytic Model (SAM) for galaxy formation and evolution with the circumgalactic medium (CGM) model presented in Faerman et al. (2020) to explore the CGM properties of $L^{*}$ galaxies. We use the SAM to generate a sample of galaxies with halo masses similar to the Milky Way (MW) halo, $M_{\rm vir} \approx 10^{12}~{\rm M_{sun}}$, and find that the CGM mass and mean me… ▽ More We combine the Santa-Cruz Semi-Analytic Model (SAM) for galaxy formation and evolution with the circumgalactic medium (CGM) model presented in Faerman et al. (2020) to explore the CGM properties of $L^{*}$ galaxies. We use the SAM to generate a sample of galaxies with halo masses similar to the Milky Way (MW) halo, $M_{\rm vir} \approx 10^{12}~{\rm M_{sun}}$, and find that the CGM mass and mean metallicity in the sample are correlated. We use the CGM masses and metallicities of the SAM galaxies as inputs for the FSM20 model, and vary the amount of non-thermal support. The density profiles in our models can be approximated by power-law functions with slopes in the range of $0.75 < a_n < 1.25$, with higher non-thermal pressure resulting in flatter distributions. We explore how the gas pressure, dispersion measure, OVI-OVIII column densities, and cooling rates behave with the gas distribution and total mass. We show that for CGM masses below $\sim 3 \times 10^{10}~{\rm M_{sun}}$, photoionization has a significant effect on the column densities of OVI and OVIII. The combination of different MW CGM observations favors models with similar fractions in thermal pressure, magnetic fields/cosmic rays, and turbulent support, and with $M_{\rm gas} \sim 3-10 \times 10^{10}~{\rm M_{sun}}$. The MW OVI column requires $t_{\rm cool}/t_{\rm dyn} \sim 4$, independent of the gas distribution. The AGN jet-driven heating rates in the SAM are enough to offset the CGM cooling, although exact balance is not required in star-forming galaxies. We provide predictions for the columns densities of additional metal ions - NV, NeVIII, and MgX. △ Less

Submitted 5 July, 2021; originally announced July 2021.

Comments: Submitted to ApJ. 19 pages, 15 figures (+ appendix)

Neutral CGM as damped Lyαあるふぁ absorbers at high redshift

Authors: Jonathan Stern, Amiel Sternberg, Claude-André Faucher-Giguère, Zachary Hafen, Drummond Fielding, Eliot Quataert, Andrew Wetzel, Daniel Anglés-Alcázar, Kareem El-Badry, Dušan Kereš, Philip F. Hopkins

Abstract: Recent searches for the hosts of high-redshift ($z \sim 4$) damped Ly$αあるふぁ$ absorbers (DLAs) have detected bright galaxies at distances of tens of kpc from the DLA. Using the FIRE-2 cosmological zoom simulations, we argue that these relatively large distances are due to a predominantly cool and neutral inner circumgalactic medium (CGM) surrounding high-redshift galaxies. The inner CGM is cool because… ▽ More Recent searches for the hosts of high-redshift ($z \sim 4$) damped Ly$αあるふぁ$ absorbers (DLAs) have detected bright galaxies at distances of tens of kpc from the DLA. Using the FIRE-2 cosmological zoom simulations, we argue that these relatively large distances are due to a predominantly cool and neutral inner circumgalactic medium (CGM) surrounding high-redshift galaxies. The inner CGM is cool because of the short cooling time of hot gas in $\lesssim10^{12}$ Msun halos, which implies that accretion and feedback energy are radiated quickly, while it is neutral due to the high volume densities and column densities at high redshift which shield cool gas from photoionization. Our analysis predicts large DLA covering factors ($\gtrsim50\%$) out to impact parameters $\sim0.3((1 + z)/5)^{3/2}\ R_{\rm vir}$ from the central galaxies at $z > 1$, equivalent to a physical distance of $\sim 21 M_{12}^{1/3} ((1 + z)/5)^{1/2}$ kpc ($R_{\rm vir}$ and $M_{12}$ are the halo virial radius and mass in units of $10^{12}$ Msun, respectively). This implies that DLA covering factors at $z \sim 4$ may be comparable to unity out to a distance $\sim 10$ times larger than stellar half-mass radii. A predominantly neutral inner CGM in the early universe suggests that its mass and metallicity can be directly constrained by CGM absorption surveys, without resorting to large ionization corrections as required for ionized CGM. △ Less

Submitted 30 July, 2021; v1 submitted 13 May, 2021; originally announced May 2021.

Comments: 16 pages, 13 figures, accepted to MNRAS

HI-to-H$_2$ Transitions in Dust-Free Interstellar Gas

Authors: Amiel Sternberg, Alon Gurman, Shmuel Bialy

Abstract: We present numerical computations and analysis of atomic to molecular (HI-to-H$_2$) transitions in cool ($\sim$100 K) low-metallicity dust-free (primordial) gas, in which molecule formation occurs via cosmic-ray driven negative ion chemistry, and removal is by a combination of far-UV photodissociation and cosmic-ray ionization and dissociation. For any gas temperature, the behavior depends on the… ▽ More We present numerical computations and analysis of atomic to molecular (HI-to-H$_2$) transitions in cool ($\sim$100 K) low-metallicity dust-free (primordial) gas, in which molecule formation occurs via cosmic-ray driven negative ion chemistry, and removal is by a combination of far-UV photodissociation and cosmic-ray ionization and dissociation. For any gas temperature, the behavior depends on the ratio of the Lyman-Werner (LW) band FUV intensity to gas density, $I_{\rm LW}/n$, and the ratio of the cosmic-ray ionization rate to the gas density, $ζぜーた/n$. We present sets of HI-to-H$_2$ abundance profiles for a wide range of $ζぜーた/n$ and $I_{\rm LW}/n$, for dust-free gas. We determine the conditions for which H$_2$ absorption line self-shielding in optically thick clouds enables a transition from atomic to molecular form for ionization-driven chemistry. We also examine the effects of cosmic-ray energy losses on the atomic and molecular density profiles and transition points. For a unit Galactic interstellar FUV field intensity ($I_{\rm LW}=1$) with LW flux $2.07\times 10^7$ photons cm$^{-2}$ s$^{-1}$, and a uniform cosmic-ray ionization rate $ζぜーた=10^{-16}$ s$^{-1}$, an HI-to-H$_2$ transition occurs at a total hydrogen gas column density of $4\times 10^{21}$ cm$^{-2}$, within $3\times 10^7$ yr, for a gas volume density of $n=10^6$ cm$^{-3}$ at 100 K. For these parameters, the dust-free limit obtains for a dust-to-gas ratio Z$^\prime_d \lesssim 10^{-5}$, which may be reached for overall metallicities $Z^\prime\lesssim 0.01$ relative to Galactic solar values. △ Less

Submitted 22 July, 2021; v1 submitted 4 May, 2021; originally announced May 2021.

Comments: Accepted for publication in ApJ

Metallicity dependence of the H/H$_2$ and C$^+$/C/CO distributions in a resolved self-regulating interstellar medium

Authors: Chia-Yu Hu, Amiel Sternberg, Ewine F. van Dishoeck

Abstract: We study the metallicity dependence of the H/H$_2$ and C$^+$/C/CO distributions in a self-regulated interstellar medium (ISM) across a broad range of metallicities ($0.1 < Z/Z_\odot < 3$). To this end, we conduct high-resolution (particle mass of $1\ {\rm M_\odot}$) hydrodynamical simulations coupled with a time-dependent H$_2$ chemistry network. The results are then post-processed with an accurat… ▽ More We study the metallicity dependence of the H/H$_2$ and C$^+$/C/CO distributions in a self-regulated interstellar medium (ISM) across a broad range of metallicities ($0.1 < Z/Z_\odot < 3$). To this end, we conduct high-resolution (particle mass of $1\ {\rm M_\odot}$) hydrodynamical simulations coupled with a time-dependent H$_2$ chemistry network. The results are then post-processed with an accurate chemistry network to model the associated C$^+$/C/CO abundances, based on the time-dependent non-steady-state (``non-equilibrium'') H$_2$ abundances. We find that the time-averaged star formation rate and the ISM structure are insensitive to metallicity. The column densities relevant for molecular shielding appear correlated with the volume densities in gravitationally unstable gas. As metallicity decreases, H$_2$ progressively deviates from steady state (``equilibrium'') and shows shallow abundance profiles until they sharply truncate at the photodissociation fronts. In contrast, the CO profile is sharp and controlled by photodissociation as CO quickly reaches steady state. We construct effective one-dimensional cloud models that successfully capture the time-averaged chemical distributions in simulations. At low metallicities, the steady-state model significantly overestimates the abundance of H$_2$ in the diffuse medium. The overestimated H$_2$, however, has little impact on CO. Consequently, the mass fraction of CO-dark H$_2$ gas is significantly lower than what a fully steady-state model predicts. The mass ratios of H$_2$/C$^+$ and H$_2$/C both show a weaker dependence on $Z^{\prime}$ than H$_2$/CO, which potentially indicates that C$^+$ and C could be alternative tracers for H$_2$ at low $Z^{\prime}$ in terms of mass budget. Our chemistry code for post-processing is publicly available. △ Less

Submitted 22 June, 2021; v1 submitted 5 March, 2021; originally announced March 2021.

Comments: ApJ accepted, 32 pages (main text 24 pages), chemistry code available at https://github.com/huchiayu/AstroChemistry.jl

The central parsec of NGC 3783: a rotating broad emission line region, asymmetric hot dust structure, and compact coronal line region

Authors: GRAVITY Collaboration, A. Amorim, M. Bauböck, W. Brandner, M. Bolzer, Y. Clénet, R. Davies, P. T. de Zeeuw, J. Dexter, A. Drescher, A. Eckart, F. Eisenhauer, N. M. Förster Schreiber, F. Gao, P. J. V. Garcia, R. Genzel, S. Gillessen, D. Gratadour, S. Hönig, D. Kaltenbrunner, M. Kishimoto, S. Lacour, D. Lutz, F. Millour, H. Netzer , et al. (25 additional authors not shown)

Abstract: Using VLTI/GRAVITY and SINFONI data, we investigate the sub-pc gas and dust structure around the nearby type 1 AGN hosted by NGC 3783. The K-band coverage of GRAVITY uniquely allows a simultaneous analysis of the size and kinematics of the broad line region (BLR), the size and structure of the near-IR continuum emitting hot dust, and the size of the coronal line region (CLR). We find the BLR probe… ▽ More Using VLTI/GRAVITY and SINFONI data, we investigate the sub-pc gas and dust structure around the nearby type 1 AGN hosted by NGC 3783. The K-band coverage of GRAVITY uniquely allows a simultaneous analysis of the size and kinematics of the broad line region (BLR), the size and structure of the near-IR continuum emitting hot dust, and the size of the coronal line region (CLR). We find the BLR probed through broad Br$γがんま$ emission is well described by a rotating, thick disk with a radial distribution of clouds peaking in the inner region. In our BLR model the physical mean radius of 16 light days is nearly twice the 10 day time lag that would be measured, which matches very well the 10 day time lag that has been measured by reverberation mapping. We measure a hot dust FWHM size of 0.74 mas (0.14 pc) and further reconstruct an image of the hot dust which reveals a faint (5% of the total flux) offset cloud which we interpret as an accreting cloud heated by the central AGN. Finally, we directly measure the FWHM size of the nuclear CLR as traced by the [CaVIII] and narrow Br$γがんま$ line. We find a FWHM size of 2.2 mas (0.4 pc), fully in line with the expectation of the CLR located between the BLR and narrow line region. Combining all of these measurements together with larger scale near-IR integral field unit and mid-IR interferometry data, we are able to comprehensively map the structure and dynamics of gas and dust from 0.01--100 pc. △ Less

Submitted 29 January, 2021; originally announced February 2021.

Comments: 20 pages and 18 figures in main text, Accepted for publication in A&A

Journal ref: A&A 648, A117 (2021)

The KMOS$^{\rm 3D}$ Survey: Investigating the Origin of the Elevated Electron Densities in Star-Forming Galaxies at $1\lesssim{z}\lesssim{3}$

Authors: Rebecca L. Davies, N. M. Förster Schreiber, R. Genzel, T. T. Shimizu, R. I. Davies, A. Schruba, L. J. Tacconi, H. Übler, E. Wisnioski, S. Wuyts, M. Fossati, R. Herrera-Camus, D. Lutz, J. T. Mendel, T. Naab, S. H. Price, A. Renzini, D. Wilman, A. Beifiori, S. Belli, A. Burkert, J. Chan, A. Contursi, M. Fabricius, M. M. Lee , et al. (2 additional authors not shown)

Abstract: We investigate what drives the redshift evolution of the typical electron density ($n_e$) in star-forming galaxies, using a sample of 140 galaxies drawn primarily from KMOS$^{\rm 3D}$ ($0.6\lesssim{z}\lesssim{2.6}$) and 471 galaxies from SAMI ($z<0.113$). We select galaxies that do not show evidence of AGN activity or outflows, to constrain the average conditions within H II regions. Measurements… ▽ More We investigate what drives the redshift evolution of the typical electron density ($n_e$) in star-forming galaxies, using a sample of 140 galaxies drawn primarily from KMOS$^{\rm 3D}$ ($0.6\lesssim{z}\lesssim{2.6}$) and 471 galaxies from SAMI ($z<0.113$). We select galaxies that do not show evidence of AGN activity or outflows, to constrain the average conditions within H II regions. Measurements of the [SII]$λらむだ$6716/[SII]$λらむだ$6731 ratio in four redshift bins indicate that the local $n_e$ in the line-emitting material decreases from 187$^{+140}_{-132}$ cm$^{-3}$ at $z\sim$ 2.2 to 32$^{+4}_{-9}$ cm$^{-3}$ at $z\sim$ 0; consistent with previous results. We use the H$αあるふぁ$ luminosity to estimate the root-mean-square (rms) $n_e$ averaged over the volumes of star-forming disks at each redshift. The local and volume-averaged $n_e$ evolve at similar rates, hinting that the volume filling factor of the line-emitting gas may be approximately constant across $0\lesssim{z}\lesssim{2.6}$. The KMOS$^{\rm 3D}$ and SAMI galaxies follow a roughly monotonic trend between $n_e$ and star formation rate, but the KMOS$^{\rm 3D}$ galaxies have systematically higher $n_e$ than the SAMI galaxies at fixed offset from the star-forming main sequence, suggesting a link between the $n_e$ evolution and the evolving main sequence normalization. We quantitatively test potential drivers of the density evolution and find that $n_e$(rms) $\simeq{n_{H_2}}$, suggesting that the elevated $n_e$ in high-$z$ H II regions could plausibly be the direct result of higher densities in the parent molecular clouds. There is also tentative evidence that $n_e$ could be influenced by the balance between stellar feedback, which drives the expansion of H II regions, and the ambient pressure, which resists their expansion. △ Less

Submitted 18 December, 2020; originally announced December 2020.

Comments: Main text 24 pages, 8 figures. Accepted for publication in ApJ

Non-equilibrium ionisation and radiative transport in an evolving supernova remnant

Authors: Kartick Chandra Sarkar, Orly Gnat, Amiel Sternberg

Abstract: We present numerical simulations of the evolution of a supernova (SN) remnant expanding into a uniform background medium with density $n_H = 1.0$ cm$^{-3}$ and temperature of $10^4$ K. We include a dynamically evolving non-equilibrium ionisation (NEI) network (consisting of all the ions of H, He, C, N, O, Ne, Mg, Si, S, Fe), frequency dependent radiation transfer (RT), thermal conduction, and a si… ▽ More We present numerical simulations of the evolution of a supernova (SN) remnant expanding into a uniform background medium with density $n_H = 1.0$ cm$^{-3}$ and temperature of $10^4$ K. We include a dynamically evolving non-equilibrium ionisation (NEI) network (consisting of all the ions of H, He, C, N, O, Ne, Mg, Si, S, Fe), frequency dependent radiation transfer (RT), thermal conduction, and a simple dust evolution model, all intra-coupled to each other and to the hydrodynamics. We assume spherical symmetry. Photo-ionisation, radiation losses, photo-heating, charge-exchange heating/cooling and radiation pressure are calculated on-the-fly depending on the local radiation field and ion fractions. We find that the dynamics and energetics (but not the emission spectra) of the SN remnants can be well modelled by collisional equilibrium cooling curves even in the absence of non-equilibrium cooling and radiative transport. We find that the effect of precursor ionising radiation at different stages of SN remnant are dominated by rapid cooling of the shock and differ from steady state shocks. The predicted column densities of different ions such as N II, C IV and N V, can be higher by up to several orders of magnitude compared to steady state shocks. We also present some high esolution emission spectra that can be compared with the observed remnants to obtain important information about the physical and chemical states of the remnant, as well as constrain the background ISM. △ Less

Submitted 1 March, 2021; v1 submitted 1 October, 2020; originally announced October 2020.

Comments: 18 pages, accepted by MNRAS. Added implications towards 1D oscillatory instability in SN remnants

A new ionisation network and radiation transport module in PLUTO

Authors: Kartick Chandra Sarkar, Amiel Sternberg, Orly Gnat

Abstract: We introduce a new general-purpose time-dependent ionisation network (IN) and a radiation transport (RT) module in the magneto-hydrodynamic (MHD) code PLUTO. Our ionisation network is reliable for temperatures ranging from 5e3 to 3e8 K, and includes all ionisation states of H, He, C, N, O, Ne, Mg, Si, S and Fe making it suitable for studying a variety of astrophysical scenarios. Radiation loss for… ▽ More We introduce a new general-purpose time-dependent ionisation network (IN) and a radiation transport (RT) module in the magneto-hydrodynamic (MHD) code PLUTO. Our ionisation network is reliable for temperatures ranging from 5e3 to 3e8 K, and includes all ionisation states of H, He, C, N, O, Ne, Mg, Si, S and Fe making it suitable for studying a variety of astrophysical scenarios. Radiation loss for each ion-electron pair is calculated using CLOUDY-17 data on-the-fly. Photo-ionisation and charge exchange are the chemical heating mechanisms. The IN is fully coupled to the radiation transport module over a very large range of opacities at different frequencies. The RT module employs a method of short characteristics assuming spherical symmetry. The radiation module requires the assumption of spherical symmetry, while the IN is compatible with full 3D. We also include a simple prescription for dust opacity, grain destruction, and the dust contribution to radiation pressure. We present numerical tests to show the reliability and limitations of the new modules. We also present a post-processing tool to calculate projected column densities and emission spectra. △ Less

Submitted 1 March, 2021; v1 submitted 1 October, 2020; originally announced October 2020.

Comments: 19 pages, accepted by MNRAS. Added error estimates for radiative transfer, and runtime plots for different components of the code. Added late time (>1Myr) dynamics of the HII region

The spatially resolved broad line region of IRAS 09149-6206

Authors: GRAVITY Collaboration, A. Amorim, W. Brandner, Y. Clénet, R. Davies, P. T. de Zeeuw, J. Dexter, A. Eckart, F. Eisenhauer, N. M. Förster Schreiber, F. Gao, P. J. V. Garcia, R. Genzel, S. Gillessen, D. Gratadour, S. Hönig, M. Kishimoto, S. Lacour, D. Lutz, F. Millour, H. Netzer, T. Ott, T. Paumard, K. Perraut, G. Perrin , et al. (19 additional authors not shown)

Abstract: We present new near-infrared VLTI/GRAVITY interferometric spectra that spatially resolve the broad Br$γがんま$ emission line in the nucleus of the active galaxy IRAS 09149-6206. We use these data to measure the size of the broad line region (BLR) and estimate the mass of the central black hole. Using an improved phase calibration method that reduces the differential phase uncertainty to 0.05 degree per… ▽ More We present new near-infrared VLTI/GRAVITY interferometric spectra that spatially resolve the broad Br$γがんま$ emission line in the nucleus of the active galaxy IRAS 09149-6206. We use these data to measure the size of the broad line region (BLR) and estimate the mass of the central black hole. Using an improved phase calibration method that reduces the differential phase uncertainty to 0.05 degree per baseline across the spectrum, we detect a differential phase signal that reaches a maximum of ~0.5 degree between the line and continuum. This represents an offset of ~120 $μみゅー$as (0.14 pc) between the BLR and the centroid of the hot dust distribution traced by the 2.3 $μみゅー$m continuum. The offset is well within the dust sublimation region, which matches the measured ~0.6 mas (0.7 pc) diameter of the continuum. A clear velocity gradient, almost perpendicular to the offset, is traced by the reconstructed photocentres of the spectral channels of the Br$γがんま$ line. We infer the radius of the BLR to be ~65 $μみゅー$as (0.075 pc), which is consistent with the radius-luminosity relation of nearby active galactic nuclei derived based on the time lag of the H$βべーた$ line from reverberation mapping campaigns. Our dynamical modelling indicates the black hole mass is $\sim 1\times10^8\,M_\odot$, which is a little below, but consistent with, the standard $M_{\rm BH}$-$σしぐま_*$ relation. △ Less

Submitted 17 September, 2020; originally announced September 2020.

Comments: Accepted for publication in A&A. 14 pages, 11 figures in main text

Journal ref: A&A 643, A154 (2020)

The Kinematics and Dark Matter Fractions of TNG50 Galaxies at z=2 from an Observational Perspective

Authors: Hannah Übler, Shy Genel, Amiel Sternberg, Reinhard Genzel, Sedona H. Price, Natascha M. Förster Schreiber, Taro T. Shimizu, Annalisa Pillepich, Dylan Nelson, Andreas Burkert, Ric Davies, Lars Hernquist, Philipp Lang, Dieter Lutz, Rüdiger Pakmor, Linda J. Tacconi

Abstract: We contrast the gas kinematics and dark matter contents of $z=2$ star-forming galaxies (SFGs) from state-of-the-art cosmological simulations within the $Λらむだ$CDM framework to observations. To this end, we create realistic mock observations of massive SFGs ($M_*>4\times10^{10} M_{\odot}$, SFR $>50~M_{\odot}$ yr$^{-1}$) from the TNG50 simulation of the IllustrisTNG suite, resembling near-infrared, adap… ▽ More We contrast the gas kinematics and dark matter contents of $z=2$ star-forming galaxies (SFGs) from state-of-the-art cosmological simulations within the $Λらむだ$CDM framework to observations. To this end, we create realistic mock observations of massive SFGs ($M_*>4\times10^{10} M_{\odot}$, SFR $>50~M_{\odot}$ yr$^{-1}$) from the TNG50 simulation of the IllustrisTNG suite, resembling near-infrared, adaptive-optics assisted integral-field observations from the ground. Using observational line fitting and modeling techniques, we analyse in detail the kinematics of seven TNG50 galaxies from five different projections per galaxy, and compare them to observations of twelve massive SFGs by Genzel et al. (2020). The simulated galaxies show clear signs of disc rotation but mostly exhibit more asymmetric rotation curves, partly due to large intrinsic radial and vertical velocity components. At identical inclination angle, their one-dimensional velocity profiles can vary along different lines of sight by up to $Δでるたv=200$ km s$^{-1}$. From dynamical modelling we infer rotation speeds and velocity dispersions that are broadly consistent with observational results. We find low central dark matter fractions compatible with observations ($f_{\rm DM}^v(<R_e)=v_{\rm DM}^2(R_e)/v_{\rm circ}^2(R_e)\sim0.32\pm0.10$), however for disc effective radii $R_e$ that are mostly too small: at fixed $R_e$ the TNG50 dark matter fractions are too high by a factor of $\sim2$. We speculate that the differences in gas kinematics and dark matter content compared to the observations may be due to physical processes that are not resolved in sufficient detail with the numerical resolution available in current cosmological simulations. △ Less

Submitted 4 November, 2020; v1 submitted 12 August, 2020; originally announced August 2020.

Comments: 25 pages, 16 figures, accepted for publication in MNRAS

Rotation Curves in z~1-2 Star-Forming Disks: Evidence for Cored Dark Matter Distributions

Authors: R. Genzel, S. H. Price, H. Übler, N. M. Förster Schreiber, T. T. Shimizu, L. J. Tacconi, R. Bender, A. Burkert, A. Contursi, R. Coogan, R. L. Davies, R. I. Davies, A. Dekel, R. Herrera-Camus, M. Lee, D. Lutz, T. Naab, R. Neri, A. Nestor, A. Renzini, R. Saglia, K. Schuster, A. Sternberg, E. Wisnioski, S. Wuyts

Abstract: We report high quality, Halpha or CO rotation curves (RCs) to several Re for 41 large, massive, star-forming disk galaxies (SFGs), across the peak of cosmic galaxy evolution (z~0.67-2.45), taken with the ESO-VLT, the LBT and IRAM-NOEMA. Most RC41 SFGs have reflection symmetric RCs plausibly described by equilibrium dynamics. We fit the major axis position-velocity cuts with beam-convolved, forward… ▽ More We report high quality, Halpha or CO rotation curves (RCs) to several Re for 41 large, massive, star-forming disk galaxies (SFGs), across the peak of cosmic galaxy evolution (z~0.67-2.45), taken with the ESO-VLT, the LBT and IRAM-NOEMA. Most RC41 SFGs have reflection symmetric RCs plausibly described by equilibrium dynamics. We fit the major axis position-velocity cuts with beam-convolved, forward modeling with a bulge, a turbulent rotating disk, and a dark matter (DM) halo. We include priors for stellar and molecular gas masses, optical light effective radii and inclinations, and DM masses from abundance matching scaling relations. Two-thirds or more of the z>1.2 SFGs are baryon dominated within a few Re of typically 5.5 kpc, and have DM fractions less than maximal disks (<fDM (Re)>=0.12). At lower redshift (z<1.2) that fraction is less than one-third. DM fractions correlate inversely with the baryonic angular momentum parameter, baryonic surface density and bulge mass. Inferred low DM fractions cannot apply to the entire disk & halo but more plausibly reflect a flattened, or cored, inner DM density distribution. The typical central 'DM deficit' in these cores relative to NFW distributions is ~30% of the bulge mass. The observations are consistent with rapid radial transport of baryons in the first generation massive gas rich halos forming globally gravitationally unstable disks, and leading to efficient build-up of massive bulges and central black holes. A combination of heating due to dynamical friction and AGN feedback may drive DM out of the initial cusps. △ Less

Submitted 19 August, 2020; v1 submitted 4 June, 2020; originally announced June 2020.

Comments: 15 figures, 4 tables, accepted for publication in ApJ

Detecting the Cosmic Web: Lyαあるふぁ Emission from Simulated Filaments at z=3

Authors: Lydia M. Elias, Shy Genel, Amiel Sternberg, Julien Devriendt, Adrianne Slyz, Eli Visbal, Nicolas Bouché

Abstract: The standard cosmological model ($Λらむだ$CDM) predicts the existence of the cosmic web: a distribution of matter into sheets and filaments connecting massive halos. However, observational evidence has been elusive due to the low surface brightness of the filaments. Recent deep MUSE/VLT data and upcoming observations offer a promising avenue for Ly$αあるふぁ$ detection, motivating the development of modern theo… ▽ More The standard cosmological model ($Λらむだ$CDM) predicts the existence of the cosmic web: a distribution of matter into sheets and filaments connecting massive halos. However, observational evidence has been elusive due to the low surface brightness of the filaments. Recent deep MUSE/VLT data and upcoming observations offer a promising avenue for Ly$αあるふぁ$ detection, motivating the development of modern theoretical predictions. We use hydrodynamical cosmological simulations run with the AREPO code to investigate the potential detectability of large-scale filaments, excluding contributions from the halos embedded in them. We focus on filaments connecting massive ($M_{200c}\sim(1-3)\times10^{12} M_\odot$) halos at z=3, and compare different simulation resolutions, feedback levels, and mock-image pixel sizes. We find increasing simulation resolution does not substantially improve detectability notwithstanding the intrinsic enhancement of internal filament structure. By contrast, for a MUSE integration of 31 hours, including feedback increases the detectable area by a factor of $\simeq$5.5 on average compared with simulations without feedback, implying that even the non-bound components of the filaments have substantial sensitivity to feedback. Degrading the image resolution from the native MUSE scale of (0.2")$^2$ per pixel to (5.3")$^2$ apertures has the strongest effect, increasing the detectable area by a median factor of $\simeq$200 and is most effective when the size of the pixel roughly matches the width of the filament. Finally, we find the majority of Ly$αあるふぁ$ emission is due to electron impact collisional excitations, as opposed to radiative recombination. △ Less

Submitted 20 April, 2020; originally announced April 2020.

From Nuclear to Circumgalactic: Zooming in on AGN-Driven Outflows at z~2.2 with SINFONI

Authors: Rebecca L. Davies, N. M. Förster Schreiber, D. Lutz, R. Genzel, S. Belli, T. T. Shimizu, A. Contursi, R. I. Davies, R. Herrera-Camus, M. M. Lee, T. Naab, S. H. Price, A. Renzini, A. Schruba, A. Sternberg, L. J. Tacconi, H. Übler, E. Wisnioski, S. Wuyts

Abstract: We use deep adaptive optics assisted integral field spectroscopy from SINFONI on the VLT to study the spatially resolved properties of ionized gas outflows driven by active galactic nuclei (AGN) in three galaxies at z~2.2 -- K20-ID5, COS4-11337 and J0901+1814. These systems probe AGN feedback from nuclear to circumgalactic scales, and provide unique insights into the different mechanisms by which… ▽ More We use deep adaptive optics assisted integral field spectroscopy from SINFONI on the VLT to study the spatially resolved properties of ionized gas outflows driven by active galactic nuclei (AGN) in three galaxies at z~2.2 -- K20-ID5, COS4-11337 and J0901+1814. These systems probe AGN feedback from nuclear to circumgalactic scales, and provide unique insights into the different mechanisms by which AGN-driven outflows interact with their host galaxies. K20-ID5 and COS4-11337 are compact star forming galaxies with powerful $\sim$1500 km s$^{-1}$ AGN-driven outflows that dominate their nuclear H$αあるふぁ$ emission. The outflows do not appear to have any impact on the instantaneous star formation activity of the host galaxies, but they carry a significant amount of kinetic energy which could heat the halo gas and potentially lead to a reduction in the rate of cold gas accretion onto the galaxies. The outflow from COS4-11337 is propagating directly towards its companion galaxy COS4-11363, at a projected separation of 5.4 kpc. COS4-11363 shows signs of shock excitation and recent truncation of star formation activity, which could plausibly have been induced by the outflow from COS4-11337. J0901+1814 is gravitationally lensed, giving us a unique view of a compact (R = 470 $\pm$ 70 pc), relatively low velocity ($\sim$650 km s$^{-1}$) AGN-driven outflow. J0901+1814 has a similar AGN luminosity to COS4-11337, suggesting that the difference in outflow properties is not related to the current AGN luminosity, and may instead reflect a difference in the evolutionary stage of the outflow and/or the coupling efficiency between the AGN ionizing radiation field and the gas in the nuclear regions. △ Less

Submitted 6 April, 2020; originally announced April 2020.

Comments: Accepted for publication in ApJ. Main text 23 pages, 15 figures and 4 tables, plus Appendix (3 pages, 3 figures, 1 table)

The Evolution of the Star-forming Interstellar Medium across Cosmic Time

Authors: Linda J. Tacconi, Reinhard Genzel, Amiel Sternberg

Abstract: Over the past decade increasingly robust estimates of the dense molecular gas content in galaxy populations between redshift 0 and the peak of cosmic galaxy/star formation from redshift 1-3 have become available. This rapid progress has been possible due to the advent of powerful ground-based, and space telescopes for combined study of several millimeter to far-IR, line or continuum tracers of the… ▽ More Over the past decade increasingly robust estimates of the dense molecular gas content in galaxy populations between redshift 0 and the peak of cosmic galaxy/star formation from redshift 1-3 have become available. This rapid progress has been possible due to the advent of powerful ground-based, and space telescopes for combined study of several millimeter to far-IR, line or continuum tracers of the molecular gas and dust components. The main conclusions of this review are: 1. Star forming galaxies contained much more molecular gas at earlier cosmic epochs than at the present time. 2. The galaxy integrated depletion time scale for converting the gas into stars depends primarily on z or Hubble time, and at a given z, on the vertical location of a galaxy along the star-formation rate versus stellar mass "main-sequence" (MS) correlation. 3. Global rates of galaxy gas accretion primarily control the evolution of the cold molecular gas content and star formation rates of the dominant MS galaxy population, which in turn vary with the cosmological expansion. A second key driver may be global disk fragmentation in high-z, gas rich galaxies, which ties local free-fall time scales to galactic orbital times, and leads to rapid radial matter transport and bulge growth. Third, the low star formation efficiency inside molecular clouds is plausibly set by super-sonic streaming motions, and internal turbulence, which in turn may be driven by conversion of gravitational energy at high-z, and/or by local feedback from massive stars at low-z. 4. A simple 'gas regulator' model is remarkably successful in predicting the combined evolution of molecular gas fractions, star formation rates, galactic winds, and gas phase metallicities. △ Less

Submitted 13 March, 2020; originally announced March 2020.

Comments: To be published in Annual Reviews of Astronomy and Astrophysics

Ionized outflows in local luminous AGN: what are the real densities and outflow rates?

Authors: R. Davies, D. Baron, T. Shimizu, H. Netzer, L. Burtscher, P. T. de Zeeuw, R. Genzel, E. K. S. Hicks, M. Koss, M. -Y. Lin, D. Lutz, W. Maciejewski, F. Müller-Sánchez, G. Orban de Xivry, C. Ricci, R. Riffel, R. A. Riffel, D. Rosario, M. Schartmann, A. Schnorr-Müller, J. Shangguan, A. Sternberg, E. Sturm, T. Storchi-Bergmann, L. Tacconi , et al. (1 additional authors not shown)

Abstract: We report on the determination of electron densities, and their impact on the outflow masses and rates, measured in the central few hundred parsecs of 11 local luminous active galaxies. We show that the peak of the integrated line emission in the AGN is significantly offset from the systemic velocity as traced by the stellar absorption features, indicating that the profiles are dominated by outflo… ▽ More We report on the determination of electron densities, and their impact on the outflow masses and rates, measured in the central few hundred parsecs of 11 local luminous active galaxies. We show that the peak of the integrated line emission in the AGN is significantly offset from the systemic velocity as traced by the stellar absorption features, indicating that the profiles are dominated by outflow. In contrast, matched inactive galaxies are characterised by a systemic peak and weaker outflow wing. We present three independent estimates of the electron density in these AGN, discussing the merits of the different methods. The electron density derived from the [SII] doublet is significantly lower than than that found with a method developed in the last decade using auroral and transauroral lines, as well as a recently introduced method based on the ionization parameter. The reason is that, for gas photoionized by an AGN, much of the [SII] emission arises in an extended partially ionized zone where the implicit assumption that the electron density traces the hydrogen density is invalid. We propose ways to deal with this situation and we derive the associated outflow rates for ionized gas, which are in the range 0.001--0.5 M$_{\odot}$yr$^{-1}$ for our AGN sample. We compare these outflow rates to the relation between $\dot{M}_{out}$ and $L_{AGN}$ in the literature, and argue that it may need to be modified and rescaled towards lower mass outflow rates. △ Less

Submitted 11 August, 2020; v1 submitted 13 March, 2020; originally announced March 2020.

Comments: accepted by MNRAS; 30 pages (the last 8 of which are extra figures in an appendix). This is an update of the original version posted

Modeling the orbital motion of Sgr A*'s near-infrared flares

Authors: The GRAVITY Collaboration, M. Bauböck, J. Dexter, R. Abuter, A. Amorim, J. P. Berger, H. Bonnet, W. Brandner, Y. Clénet, V. Coudé du Foresto, P. T. de Zeeuw, G. Duvert, A. Eckart, F. Eisenhauer, N. M. Förster Schreiber, F. Gao, P. Garcia, E. Gendron, R. Genzel, O. Gerhard, S. Gillessen, M. Habibi, X. Haubois, T. Henning, S. Hippler , et al. (31 additional authors not shown)

Abstract: Infrared observations of Sgr A* probe the region close to the event horizon of the black hole at the Galactic center. These observations can constrain the properties of low-luminosity accretion as well as that of the black hole itself. The GRAVITY instrument at the ESO VLTI has recently detected continuous circular relativistic motion during infrared flares which has been interpreted as orbital mo… ▽ More Infrared observations of Sgr A* probe the region close to the event horizon of the black hole at the Galactic center. These observations can constrain the properties of low-luminosity accretion as well as that of the black hole itself. The GRAVITY instrument at the ESO VLTI has recently detected continuous circular relativistic motion during infrared flares which has been interpreted as orbital motion near the event horizon. Here we analyze the astrometric data from these flares, taking into account the effects of out-of-plane motion and orbital shear of material near the event horizon of the black hole. We have developed a new code to predict astrometric motion and flux variability from compact emission regions following particle orbits. Our code combines semi-analytic calculations of timelike geodesics that allow for out-of-plane or elliptical motions with ray tracing of photon trajectories to compute time-dependent images and light curves. We apply our code to the three flares observed with GRAVITY in 2018. We show that all flares are consistent with a hotspot orbiting at R$\sim$9 gravitational radii with an inclination of $i\sim140^\circ$. The emitting region must be compact and less than $\sim5$ gravitational radii in diameter. We place a further limit on the out-of-plane motion during the flare. △ Less

Submitted 19 February, 2020; originally announced February 2020.

Comments: Accepted in A&A

Molecular Gas Inflows and Outflows in Ultraluminous Infrared Galaxies at $z\sim0.2$ and one QSO at $z=6.1$

Authors: R. Herrera-Camus, E. Sturm, J. Graciá-Carpio, S. Veilleux, T. Shimizu, D. Lutz, M. Stone, E. González-Alfonso, R. Davies, J. Fischer, R. Genzel, R. Maiolino, A. Sternberg, L. Tacconi, A. Verma

Abstract: Aims. We aim to search and characterize inflows and outflows of molecular gas in four ultraluminous infrared galaxies (ULIRGs) at $z\sim0.2-0.3$ and one distant QSO at $z=6.13$. Methods. We use Herschel PACS and ALMA Band 7 observations of the hydroxyl molecule (OH) line at rest-frame wavelength 119 $μみゅー$m which in absorption can provide unambiguous evidence for inflows or outflows of molecular ga… ▽ More Aims. We aim to search and characterize inflows and outflows of molecular gas in four ultraluminous infrared galaxies (ULIRGs) at $z\sim0.2-0.3$ and one distant QSO at $z=6.13$. Methods. We use Herschel PACS and ALMA Band 7 observations of the hydroxyl molecule (OH) line at rest-frame wavelength 119 $μみゅー$m which in absorption can provide unambiguous evidence for inflows or outflows of molecular gas in nuclear regions of galaxies. Our study contributes to double the number of OH observations of luminous systems at $z\sim0.2-0.3$, and push the search for molecular outflows based on the OH transition to $z\sim6$. Results. We detect OH high-velocity absorption wings in three of the four ULIRGs. In two cases, IRAS F20036-1547 and IRAS F13352+6402, the blueshifted absorption profiles indicate the presence of powerful and fast molecular gas outflows. Consistent with an inside-out quenching scenario, these outflows are depleting the central reservoir of molecular gas at a similar rate than the intense star formation activity. In the case of the starburst-dominated system IRAS 10091+4704, we detect an inverted P-Cygni profile that is unique among ULIRGs and indicates the presence of a fast ($\sim400$ km s$^{-1}$) inflow of molecular gas at a rate of $\sim100~M_{\odot}~{\rm yr}^{-1}$ towards the central region. Finally, we tentatively detect ($\sim3σしぐま$) the OH doublet in absorption in the $z=6.13$ QSO ULAS J131911+095051. The OH feature is blueshifted with a median velocity that suggests the presence of a molecular outflow, although characterized by a modest molecular mass loss rate of $\sim200~M_{\odot}~{\rm yr}^{-1}$. This value is comparable to the small mass outflow rates found in the stacking of the [CII] spectra of other $z\sim6$ QSOs and suggests that ejective feedback in this phase of the evolution of ULAS J131911+095051 has subsided. △ Less

Submitted 11 December, 2019; originally announced December 2019.

Comments: 9 pages, 6 figures, accepted for publication in A&A as a Letter to the Editor

An image of the dust sublimation region in the nucleus of NGC 1068

Authors: GRAVITY Collaboration, O. Pfuhl, R. Davies, J. Dexter, H. Netzer, S. Hoenig, D. Lutz, M. Schartmann, E. Sturm, A. Amorim, W. Brandner, Y. Clenet, P. T. de Zeeuw, A. Eckart, F. Eisenhauer, N. M. Foerster Schreiber, F. Gao, P. J. V. Garcia, R. Genzel, S. Gillessen, D. Gratadour, M. Kishimoto, S. Lacour, F. Millour, T. Ott , et al. (18 additional authors not shown)

Abstract: We present near-infrared interferometric data on the Seyfert 2 galaxy NGC 1068, obtained with the GRAVITY instrument on the European Southern Observatory Very Large Telescope Interferometer. The extensive baseline coverage from 5 to 60 Mλらむだallowed us to reconstruct a continuum image of the nucleus with an unrivaled 0.2 pc resolution in the K-band. We find a thin ring-like structure of emission with… ▽ More We present near-infrared interferometric data on the Seyfert 2 galaxy NGC 1068, obtained with the GRAVITY instrument on the European Southern Observatory Very Large Telescope Interferometer. The extensive baseline coverage from 5 to 60 Mλらむだallowed us to reconstruct a continuum image of the nucleus with an unrivaled 0.2 pc resolution in the K-band. We find a thin ring-like structure of emission with a radius r = 0.24+/-0.03 pc, inclination i = 70+/-5 deg, position angle PA = -50+/-4 deg, and h/r < 0.14, which we associate with the dust sublimation region. The observed morphology is inconsistent with the expected signatures of a geometrically and optically thick torus. Instead, the infrared emission shows a striking resemblance to the 22 GHz maser disc, which suggests they share a common region of origin. The near-infrared spectral energy distribution indicates a bolometric luminosity of (0.4-4.7) x 10^45 erg/s, behind a large A_K ~ 5.5 (A_V ~ 90) screen of extinction that also appears to contribute significantly to obscuring the broad line region. △ Less

Submitted 3 December, 2019; originally announced December 2019.

Comments: Accepted for publication in Astronomy and Astrophysics. 13 pages with 7 figures

Journal ref: A&A 634, A1 (2020)

Molecular outflows in local galaxies: Method comparison and a role of intermittent AGN driving

Authors: D. Lutz, E. Sturm, A. Janssen, S. Veilleux, S. Aalto, C. Cicone, A. Contursi, R. I. Davies, C. Feruglio, J. Fischer, S. Garcia-Burillo, R. Genzel, E. González-Alfonso, J. Gracía-Carpio R. Herrera-Camus, R. Maiolino, A. Schruba, T. Shimizu, A. Sternberg, L. J. Tacconi, A. Weiß

Abstract: We report new detections and limits from a NOEMA and ALMA CO(1-0) search for molecular outflows in 13 local galaxies with high FIR surface brightness, and combine with results from the literature. CO line ratios and outflow structure provide some constraints on the conversion from observables to quantities such as molecular mass outflow rates. Ratios between outflow emission in higher J CO transit… ▽ More We report new detections and limits from a NOEMA and ALMA CO(1-0) search for molecular outflows in 13 local galaxies with high FIR surface brightness, and combine with results from the literature. CO line ratios and outflow structure provide some constraints on the conversion from observables to quantities such as molecular mass outflow rates. Ratios between outflow emission in higher J CO transitions and in CO(1-0) typically are consistent with excitation Ri1<~1. For IRAS 13120-5453, however, R31=2.10 indicates optically thin CO in the outflow. Like much of the outflow literature, we use alpha(CO) = 0.8, and we present arguments for using C=1 in deriving molecular mass outflow rates Mdot = C*M*v/R. We compare the two main methods for molecular outflow detection: CO mm interferometry and Herschel OH spectroscopy. For 26 sources studied with both methods, we find 80% agreement in detecting vout>~150km/s outflows, and non-matches can be plausibly ascribed to outflow geometry and SNR. For 12 bright ULIRGs with detailed OH-based outflow modeling, CO outflows are detected in all but one. Outflow masses, velocities, and sizes for these 11 sources agree well between the two methods, and modest remaining differences may relate to the different but overlapping regions sampled by CO emission and OH absorption. Outflow properties correlate better with AGN luminosity and with bolometric luminosity than with FIR surface brightness. The most massive outflows are found for systems with current AGN activity, but significant outflows in non-AGN systems must relate to star formation or to AGN activity in the recent past. We report scaling relations for the increase of outflow mass, rate, momentum rate, and kinetic power with bolometric luminosity. Short ~10^6yr flow times and some sources with resolved multiple outflow episodes support a role of intermittent driving, likely by AGN. (abridged) △ Less

Submitted 13 November, 2019; originally announced November 2019.

Comments: 37 pages. Accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 633, A134 (2020)

The resolved size and structure of hot dust in the immediate vicinity of AGN

Authors: GRAVITY Collaboration, J. Dexter, J. Shangguan, S. Hönig, M. Kishimoto, D. Lutz, H. Netzer, R. Davies, E. Sturm, O. Pfuhl, A. Amorim, M. Bauböck, W. Brandner, Y. Clénet, P. T. de Zeeuw, A. Eckart, F. Eisenhauer, N. M. Förster Schreiber, F. Gao, P. J. V. Garcia, R. Genzel, S. Gillessen, D. Gratadour, A. Jiménez-Rosales, S. Lacour , et al. (20 additional authors not shown)

Abstract: We use VLTI/GRAVITY near-infrared interferometry measurements of 8 bright, Type 1 AGN to study the size and structure of hot dust heated by the central engine. We partially resolve each source, and report Gaussian FWHM sizes in the range 0.3-0.8 milliarcseconds. In all but one object, we find no evidence for significant elongation or asymmetry (closure phases < 1 deg). The effective physical radiu… ▽ More We use VLTI/GRAVITY near-infrared interferometry measurements of 8 bright, Type 1 AGN to study the size and structure of hot dust heated by the central engine. We partially resolve each source, and report Gaussian FWHM sizes in the range 0.3-0.8 milliarcseconds. In all but one object, we find no evidence for significant elongation or asymmetry (closure phases < 1 deg). The effective physical radius increases with bolometric luminosity as found from past reverberation and interferometry measurements. The measured sizes for Seyfert galaxies are systematically larger than for the two quasars in our sample when measured relative to the previously reported R ~ L^1/2 relationship explained by emission at the sublimation radius. This could be evidence of evolving near-infrared emission region structure as a function of central luminosity. △ Less

Submitted 1 October, 2019; originally announced October 2019.

Comments: 13 pages, 9 figures, submitted to A&A

Journal ref: A&A 635, A92 (2020)

Chemical Abundances in a Turbulent Medium -- H$_2$, OH$^+$, H$_2$O$^+$, ArH$^+$

Authors: Shmuel Bialy, David Neufeld, Mark Wolfire, Amiel Sternberg, Blakesley Burkhart

Abstract: Supersonic turbulence results in strong density fluctuations in the interstellar medium (ISM), which have a profound effect on the chemical structure. Particularly useful probes of the diffuse ISM are the ArH$^+$, OH$^+$, H$_2$O$^+$ molecular ions, which are highly sensitive to fluctuations in the density and the H$_2$ abundance. We use isothermal magnetohydrodynamic (MHD) simulations of various s… ▽ More Supersonic turbulence results in strong density fluctuations in the interstellar medium (ISM), which have a profound effect on the chemical structure. Particularly useful probes of the diffuse ISM are the ArH$^+$, OH$^+$, H$_2$O$^+$ molecular ions, which are highly sensitive to fluctuations in the density and the H$_2$ abundance. We use isothermal magnetohydrodynamic (MHD) simulations of various sonic Mach numbers, $\mathcal{M}_s$, and density decorrelation scales, $y_{\rm dec}$, to model the turbulent density field. We post-process the simulations with chemical models and obtain the probability density functions (PDFs) for the H$_2$, ArH$^+$, OH$^+$ and H$_2$O$^+$ abundances. We find that the PDF dispersions increases with increasing $\mathcal{M}_s$ and $y_{\rm dec}$, as the magnitude of the density fluctuations increases, and as they become more coherent. Turbulence also affects the median abundances: when $\mathcal{M}_s$ and $y_{\rm dec}$ are high, low-density regions with low H$_2$ abundance become prevalent, resulting in an enhancement of ArH$^+$ compared to OH$^+$ and H$_2$O$^+$. We compare our models with Herschel observations. The large scatter in the observed abundances, as well as the high observed ArH$^+$/OH$^+$ and ArH$^+$/H$_2$O$^+$ ratios are naturally reproduced by our supersonic $(\mathcal{M}_s=4.5)$, large decorrelation scale $(y_{\rm dec}=0.8)$ model, supporting a scenario of a large-scale turbulence driving. The abundances also depend on the UV intensity, CR ionization rate, and the cloud column density, and the observed scatter may be influenced by fluctuations in these parameters. △ Less

Submitted 9 October, 2019; v1 submitted 26 September, 2019; originally announced September 2019.

Comments: Accepted for publication in the Astrophysical Journal

Journal ref: The Astrophysical Journal, Volume 885, Number 2, Pages 109, 2019

Massive Warm/Hot Galaxy Coronae: II. Isentropic Model

Authors: Yakov Faerman, Amiel Sternberg, Christopher F. McKee

Abstract: We construct a new analytic phenomenological model for the extended circumgalactic material (CGM) of $L^*$ galaxies. Our model reproduces the OVII/OVIII absorption observations of the Milky Way (MW) and the OVI measurements reported by the COS-Halos and eCGM surveys. The warm/hot gas is in hydrostatic equilibrium in a MW gravitational potential, and we adopt a barotropic equation of state, resulti… ▽ More We construct a new analytic phenomenological model for the extended circumgalactic material (CGM) of $L^*$ galaxies. Our model reproduces the OVII/OVIII absorption observations of the Milky Way (MW) and the OVI measurements reported by the COS-Halos and eCGM surveys. The warm/hot gas is in hydrostatic equilibrium in a MW gravitational potential, and we adopt a barotropic equation of state, resulting in a temperature variation as a function of radius. A pressure component with an adiabatic index of $γがんま=4/3$ is included to approximate the effects of a magnetic field and cosmic rays. We introduce a metallicity gradient motivated by the enrichment of the inner CGM by the Galaxy. We then present our fiducial model for the corona, tuned to reproduce the observed OVI-OVIII column densities, and with a total mass of $M_{\rm gas} \approx 5.5 \times 10^{10}~{\rm M_{\odot}}$ inside $r_{\rm cgm} \approx 280$ kpc. The gas densities in the CGM are low ($n_{\rm H} = 10^{-5} - 3 \times 10^{-4}~{\rm cm^{-3}}$) and its collisional ionization state is modified by the metagalactic radiation field (MGRF). We show that for OVI-bearing warm/hot gas with typical observed column densities $N_{\rm OVI} \sim 3 \times 10^{14}~{\rm cm^{-2}}$ at large ($\gtrsim 100$ kpc) impact parameters from the central galaxies, the ratio of the cooling to dynamical times, $t_{\rm cool}/t_{\rm dyn}$, has a model-independent upper limit of $\lesssim 4$. In our model, $t_{\rm cool}/t_{\rm dyn}$ at large radii is $\sim 2-3$. We present predictions for a wide range of future observations of the warm/hot CGM, from UV/X-ray absorption and emission spectroscopy, to dispersion measure (DM) and Sunyaev-Zeldovich CMB measurements. We provide the model outputs in machine-readable data files, for easy comparison and analysis. △ Less

Submitted 21 March, 2020; v1 submitted 19 September, 2019; originally announced September 2019.

Comments: 24 pages, 15 figures. Accepted to ApJ

The multiphase gas structure and kinematics in the circumnuclear region of NGC 5728

Authors: T. Taro Shimizu, R. I. Davies, D. Lutz, L. Burtscher, M. Lin, D. Baron, R. L. Davies, R. Genzel, E. K. S. Hicks, M. Koss, W. Maciejewski, F. Mu ller-Sanchez, G. Orban de Xivry, S. H. Price, C. Ricci, R. Riffel, R. A. Riffel, D. Rosario, M. Schartmann, A. Schnorr-Muller, A. Sternberg, E. Sturm, T. Storchi-Bergmann, L. Tacconi, S. Veilleux

Abstract: We report on our combined analysis of HST, VLT/MUSE, VLT/SINFONI, and ALMA observations of the local Seyfert 2 galaxy, NGC 5728 to investigate in detail the feeding and feedback of the AGN. The datasets simultaneously probe the morphology, excitation, and kinematics of the stars, ionized gas, and molecular gas over a large range of spatial scales (10 pc--10 kpc). NGC 5728 contains a large stellar… ▽ More We report on our combined analysis of HST, VLT/MUSE, VLT/SINFONI, and ALMA observations of the local Seyfert 2 galaxy, NGC 5728 to investigate in detail the feeding and feedback of the AGN. The datasets simultaneously probe the morphology, excitation, and kinematics of the stars, ionized gas, and molecular gas over a large range of spatial scales (10 pc--10 kpc). NGC 5728 contains a large stellar bar which is driving gas along prominent dust lanes to the inner 1 kpc where the gas settles into a circumnuclear ring. The ring is strongly star forming and contains a substantial population of young stars as indicated by the lowered stellar velocity dispersion and gas excitation consistent with HII regions. We model the kinematics of the ring using the velocity field of the CO (2--1) emission and stars and find it is consistent with a rotating disk. The outer regions of the disk, where the dust lanes meet the ring, show signatures of inflow at a rate of 1 M$_{\sun}$ yr$^{-1}$. Inside the ring, we observe three molecular gas components corresponding to the circular rotation of the outer ring, a warped disk, and the nuclear stellar bar. The AGN is driving an ionized gas outflow that reaches a radius of 250 pc with a mass outflow rate of 0.08 M$_{\sun}$ yr$^{-1}$ consistent with its luminosity and scaling relations from previous studies. While we observe distinct holes in CO emission which could be signs of molecular gas removal, we find that largely the AGN is not disrupting the structure of the circumnuclear region. △ Less

Submitted 2 October, 2019; v1 submitted 8 July, 2019; originally announced July 2019.

Comments: accepted for publication in MNRAS; 22 pages + Appendix, 16 figures

A geometric distance measurement to the Galactic Center black hole with 0.3% uncertainty

Authors: R. Abuter, A. Amorim, M. Bauboeck, J. P. Berger, H. Bonnet, W. Brandner, Y. Clenet, V. Coude du Foresto, P. T. de Zeeuw, J. Dexter, G. Duvert, A. Eckart, F. Eisenhauer, N. M. Foerster Schreiber, P. Garcia, F. Gao, E. Gendron, R. Genzel, O. Gerhard, S. Gillessen, M. Habibi, X. Haubois, T. Henning, S. Hippler, M. Horrobin , et al. (29 additional authors not shown)

Abstract: We present a 0.16% precise and 0.27% accurate determination of R0, the distance to the Galactic Center. Our measurement uses the star S2 on its 16-year orbit around the massive black hole Sgr A* that we followed astrometrically and spectroscopically for 27 years. Since 2017, we added near-infrared interferometry with the VLTI beam combiner GRAVITY, yielding a direct measurement of the separation v… ▽ More We present a 0.16% precise and 0.27% accurate determination of R0, the distance to the Galactic Center. Our measurement uses the star S2 on its 16-year orbit around the massive black hole Sgr A* that we followed astrometrically and spectroscopically for 27 years. Since 2017, we added near-infrared interferometry with the VLTI beam combiner GRAVITY, yielding a direct measurement of the separation vector between S2 and Sgr A* with an accuracy as good as 20 micro-arcsec in the best cases. S2 passed the pericenter of its highly eccentric orbit in May 2018, and we followed the passage with dense sampling throughout the year. Together with our spectroscopy, in the best cases with an error of 7 km/s, this yields a geometric distance estimate: R0 = 8178 +- 13(stat.) +- 22(sys.) pc. This work updates our previous publication in which we reported the first detection of the gravitational redshift in the S2 data. The redshift term is now detected with a significance level of 20 sigma with f_redshift = 1.04 +- 0.05. △ Less

Submitted 11 April, 2019; originally announced April 2019.

Comments: 9 pages, 7 figures, submitted to A&A

Journal ref: 2019, A&A, 625, L10

Thermal Phases of the Neutral Atomic Interstellar Medium -- from Solar Metallicity to Primordial Gas

Authors: Shmuel Bialy, Amiel Sternberg

Abstract: We study the thermal structure of the neutral atomic (H {\small I}) interstellar medium across a wide range of metallicities, from supersolar down to vanishing metallicity, and for varying UV intensities and cosmic-ray ionization rates. We calculate self-consistently the gas temperature and species abundances (with a special focus on the residual H$_2$), assuming thermal and chemical steady-state.… ▽ More We study the thermal structure of the neutral atomic (H {\small I}) interstellar medium across a wide range of metallicities, from supersolar down to vanishing metallicity, and for varying UV intensities and cosmic-ray ionization rates. We calculate self-consistently the gas temperature and species abundances (with a special focus on the residual H$_2$), assuming thermal and chemical steady-state. For solar metallicity, $Z' \equiv 1$, we recover the known result that there exists a pressure range over which the gas is multiphased, with the warm ($\sim 10^4$ K, WNM) and cold ($\sim 100$ K, CNM) phases coexisting at the same pressure. At a metallicity $Z' \approx 0.1$, the CNM is colder (compared to $Z'=1$) due to the reduced efficiency of photoelectric heating. For $Z' \lesssim 0.1$, cosmic-ray ionization becomes the dominant heating mechanism and the WNM-to-CNM transition shifts to ever increasing pressure/density as the metallicity is reduced. For metallicities $Z' \lesssim 0.01$, H$_2$ cooling becomes important, lowering the temperature of the WNM (down to $\approx 600$ K), and smoothing out the multiphase phenomenon. At vanishing metallicities, H$_2$ heating becomes effective and the multiphase phenomenon disappears entirely. We derive analytic expressions for the critical densities for the warm-to-cold phase transition in the different regimes, and the critical metallicities for H$_2$ cooling and heating. We discuss potential implications on the star-formation rates of galaxies and self-regulation theories. △ Less

Submitted 20 June, 2019; v1 submitted 18 February, 2019; originally announced February 2019.

Comments: Accepted for publication in ApJ

Journal ref: The Astrophysical Journal, Volume 881, Number 2, 2019

Test of Einstein equivalence principle near the Galactic center supermassive black hole

Authors: A. Amorim, M. Bauböck, J. P. Berger, W. Brandner, Y. Clénet, V. Coudé du Foresto, P. T. de Zeeuw, J. Dexter, G. Duvert, M. Ebert, A. Eckart, F. Eisenhauer, N. M. Förster Schreiber, P. Garcia, F. Gao, E. Gendron, R. Genzel, S. Gillessen, M. Habibi, X. Haubois, Th. Henning, S. Hippler, M. Horrobin, Z. Hubert, A. Jiménez Rosales , et al. (27 additional authors not shown)

Abstract: During its orbit around the four million solar mass black hole Sagittarius A* the star S2 experiences significant changes in gravitational potential. We use this change of potential to test one part of the Einstein equivalence principle: the local position invariance (LPI). We study the dependency of different atomic transitions on the gravitational potential to give an upper limit on violations o… ▽ More During its orbit around the four million solar mass black hole Sagittarius A* the star S2 experiences significant changes in gravitational potential. We use this change of potential to test one part of the Einstein equivalence principle: the local position invariance (LPI). We study the dependency of different atomic transitions on the gravitational potential to give an upper limit on violations of the LPI. This is done by separately measuring the redshift from hydrogen and helium absorption lines in the stellar spectrum during its closest approach to the black hole. For this measurement we use radial velocity data from 2015 to 2018 and combine it with the gravitational potential at the position of S2, which is calculated from the precisely known orbit of S2 around the black hole. This results in a limit on a violation of the LPI of $|βべーた_{He}-βべーた_{H}| = (2.4 \pm 5.1) \cdot 10^{-2}$. The variation in potential that we probe with this measurement is six magnitudes larger than possible for measurements on Earth, and a factor ten larger than in experiments using white dwarfs. We are therefore testing the LPI in a regime where it has not been tested before. △ Less

Submitted 25 March, 2019; v1 submitted 11 February, 2019; originally announced February 2019.

Comments: Accepted for publication in Physical Review Letters

Journal ref: Phys. Rev. Lett. 122, 101102 (2019)

PHIBSS2: survey design and z=0.5-0.8 results. Molecular gas reservoirs during the winding-down of star formation

Authors: J. Freundlich, F. Combes, L. J. Tacconi, R. Genzel, S. Garcia-Burillo, R. Neri, T. Contini, A. Bolatto, S. Lilly, P. Salomé, I. C. Bicalho, J. Boissier, F. Boone, N. Bouché, F. Bournaud, A. Burkert, M. Carollo, M. C. Cooper, P. Cox, C. Feruglio, N. M. Förster Schreiber, S. Juneau, M. Lippa, D. Lutz, T. Naab , et al. (7 additional authors not shown)

Abstract: Following the success of the Plateau de Bure high-z Blue Sequence Survey (PHIBSS), we present the PHIBSS2 legacy program, a survey of the molecular gas properties of star-forming galaxies on and around the star formation main sequence (MS) at different redshifts using NOEMA. This survey significantly extends the existing sample of star-forming galaxies with CO molecular gas measurements, probing t… ▽ More Following the success of the Plateau de Bure high-z Blue Sequence Survey (PHIBSS), we present the PHIBSS2 legacy program, a survey of the molecular gas properties of star-forming galaxies on and around the star formation main sequence (MS) at different redshifts using NOEMA. This survey significantly extends the existing sample of star-forming galaxies with CO molecular gas measurements, probing the peak epoch of star formation (z=1-1.6) as well as its building-up (z=2-3) and winding-down (z=0.5-0.8) phases. The targets are drawn from the GOODS, COSMOS, and AEGIS deep fields and uniformly sample the MS in the stellar mass (M*) - star formation rate (SFR) plane with log(M*/Msun) = 10-11.8. We describe the survey strategy and sample selection before focusing on the results obtained at z=0.5-0.8, where we report 60 CO(2-1) detections out of 61 targets. We determine their molecular gas masses and separately obtain disc sizes and bulge-to-total (B/T) luminosity ratios from HST I-band images. The median molecular gas-to-stellar mass ratio, gas fraction, and depletion time as well as their dependence with M* and offset from the MS follow published scaling relations for a much larger sample of galaxies spanning a wider range of redshifts. The galaxy-averaged Kennicutt-Schmidt (KS) relation between molecular gas and SFR surface densities is strikingly linear, pointing towards similar star formation timescales within galaxies at any given epoch. In terms of morphology, the molecular gas content, the SFR, the disc stellar mass, and the disc molecular gas fraction do not seem to correlate with B/T and the stellar surface density, which suggests an ongoing supply of fresh molecular gas to compensate for the build-up of the bulge. Our measurements do not yield any significant variation of the depletion time with B/T and hence no strong evidence for morphological quenching within the scatter of the MS. △ Less

Submitted 17 January, 2019; v1 submitted 19 December, 2018; originally announced December 2018.

Comments: Accepted for publication in A&A (28 pages, 9 figures, 4 tables)

Journal ref: A&A 622, A105 (2019)

Detection of orbital motions near the last stable circular orbit of the massive black hole SgrA*

Authors: GRAVITY Collaboration, R. Abuter, A. Amorim, M. Bauböck, J. P. Berger, H. Bonnet, W. Brandner, Y. Clénet, V. Coudé du Foresto, P. T. de Zeeuw, C. Deen, J. Dexter, G. Duvert, A. Eckart, F. Eisenhauer, N. M. Förster Schreiber, P. Garcia, F. Gao, E. Gendron, R. Genzel, S. Gillessen, P. Guajardo, M. Habibi, X. Haubois, Th. Henning , et al. (35 additional authors not shown)

Abstract: We report the detection of continuous positional and polarization changes of the compact source SgrA* in high states ('flares') of its variable near- infrared emission with the near-infrared GRAVITY-Very Large Telescope Interferometer (VLTI) beam-combining instrument. In three prominent bright flares, the position centroids exhibit clockwise looped motion on the sky, on scales of typically 150 mic… ▽ More We report the detection of continuous positional and polarization changes of the compact source SgrA* in high states ('flares') of its variable near- infrared emission with the near-infrared GRAVITY-Very Large Telescope Interferometer (VLTI) beam-combining instrument. In three prominent bright flares, the position centroids exhibit clockwise looped motion on the sky, on scales of typically 150 micro-arcseconds over a few tens of minutes, corresponding to about 30% the speed of light. At the same time, the flares exhibit continuous rotation of the polarization angle, with about the same 45(+/-15)-minute period as that of the centroid motions. Modelling with relativistic ray tracing shows that these findings are all consistent with a near face-on, circular orbit of a compact polarized 'hot spot' of infrared synchrotron emission at approximately six to ten times the gravitational radius of a black hole of 4 million solar masses. This corresponds to the region just outside the innermost, stable, prograde circular orbit (ISCO) of a Schwarzschild-Kerr black hole, or near the retrograde ISCO of a highly spun-up Kerr hole. The polarization signature is consistent with orbital motion in a strong poloidal magnetic field. △ Less

Submitted 30 October, 2018; originally announced October 2018.

Comments: accepted by A&A; 16 pages

Kiloparsec Scale Properties of Star-Formation Driven Outflows at z~2.3 in the SINS/zC-SINF AO Survey

Authors: Rebecca L. Davies, Natascha M. Förster Schreiber, Hannah Übler, Reinhard Genzel, Dieter Lutz, Alvio Renzini, Sandro Tacchella, Linda J. Tacconi, Sirio Belli, Andreas Burkert, C. Marcella Carollo, Richard I. Davies, Rodrigo Herrera-Camus, Simon J. Lilly, Chiara Mancini, Thorsten Naab, Erica J. Nelson, Sedona H. Price, Thomas Taro Shimizu, Amiel Sternberg, Emily Wisnioski, Stijn Wuyts

Abstract: We investigate the relationship between star formation activity and outflow properties on kiloparsec scales in a sample of 28 star forming galaxies at $z\sim$ 2-2.6, using adaptive optics assisted integral field observations from SINFONI on the VLT. The narrow and broad components of the H$αあるふぁ$ emission are used to simultaneously determine the local star formation rate surface density (… ▽ More We investigate the relationship between star formation activity and outflow properties on kiloparsec scales in a sample of 28 star forming galaxies at $z\sim$ 2-2.6, using adaptive optics assisted integral field observations from SINFONI on the VLT. The narrow and broad components of the H$αあるふぁ$ emission are used to simultaneously determine the local star formation rate surface density ($Σしぐま_{\rm SFR}$), and the outflow velocity $v_{\rm out}$ and mass outflow rate $\dot{M}_{\rm out}$, respectively. We find clear evidence for faster outflows with larger mass loading factors at higher $Σしぐま_{\rm SFR}$. The outflow velocities scale as $v_{\rm out}$ $\propto$ $Σしぐま_{\rm SFR}^{0.34 \pm 0.10}$, which suggests that the outflows may be driven by a combination of mechanical energy released by supernova explosions and stellar winds, as well as radiation pressure acting on dust grains. The majority of the outflowing material does not have sufficient velocity to escape from the galaxy halos, but will likely be re-accreted and contribute to the chemical enrichment of the galaxies. In the highest $Σしぐま_{\rm SFR}$ regions the outflow component contains an average of $\sim$45% of the H$αあるふぁ$ flux, while in the lower $Σしぐま_{\rm SFR}$ regions only $\sim$10% of the H$αあるふぁ$ flux is associated with outflows. The mass loading factor, $ηいーた$ = $\dot{M}_{\rm out}$/SFR, is positively correlated with $Σしぐま_{\rm SFR}$ but is relatively low even at the highest $Σしぐま_{\rm SFR}$: $ηいーた\lesssim$ 0.5 $\times$ (380 cm$^{-3}$/n$_e$). This may be in tension with the $ηいーた$ $\gtrsim$ 1 required by cosmological simulations, unless a significant fraction of the outflowing mass is in other gas phases and has sufficient velocity to escape the galaxy halos. △ Less

Submitted 14 February, 2019; v1 submitted 31 August, 2018; originally announced August 2018.

Comments: 28 pages, 12 figures, 2 tables. Accepted for publication in ApJ

The KMOS^3D Survey: Demographics and Properties of Galactic Outflows at z = 0.6 - 2.7

Authors: N. M. Förster Schreiber, H. Übler, R. L. Davies, R. Genzel, E. Wisnioski, S. Belli, T. Shimizu, D. Lutz, M. Fossati, R. Herrera-Camus, J. T. Mendel, L. J. Tacconi, D. Wilman, A. Beifiori, G. Brammer, A. Burkert, C. M. Carollo, R. I. Davies, F. Eisenhauer, M. Fabricius, S. J. Lilly, I. Momcheva, T. Naab, E. J. Nelson, S. Price , et al. (5 additional authors not shown)

Abstract: We present a census of ionized gas outflows in 599 normal galaxies at redshift 0.6<z<2.7, mostly based on integral field spectroscopy of Ha, [NII], and [SII] line emission. The sample fairly homogeneously covers the main sequence of star-forming galaxies with masses 9.0<log(M*/Msun)<11.7, and probes into the regimes of quiescent galaxies and starburst outliers. About 1/3 exhibits the high-velocity… ▽ More We present a census of ionized gas outflows in 599 normal galaxies at redshift 0.6<z<2.7, mostly based on integral field spectroscopy of Ha, [NII], and [SII] line emission. The sample fairly homogeneously covers the main sequence of star-forming galaxies with masses 9.0<log(M*/Msun)<11.7, and probes into the regimes of quiescent galaxies and starburst outliers. About 1/3 exhibits the high-velocity component indicative of outflows, roughly equally split into winds driven by star formation (SF) and active galactic nuclei (AGN). The incidence of SF-driven winds correlates mainly with star formation properties. These outflows have typical velocities of ~450 km/s, local electron densities of n_e~380 cm^-3, modest mass loading factors of ~0.1-0.2 at all galaxy masses, and energetics compatible with momentum driving by young stellar populations. The SF-driven winds may escape from log(M*/Msun)<10.3 galaxies but substantial mass, momentum, and energy in hotter and colder outflow phases seem required to account for low galaxy formation efficiencies in the low-mass regime. Faster AGN-driven outflows (~1000-2000 km/s) are commonly detected above log(M*/Msun)~10.7, in up to ~75% of log(M*/Msun)>11.2 galaxies. The incidence, strength, and velocity of AGN-driven winds strongly correlates with stellar mass and central concentration. Their outflowing ionized gas appears denser (n_e~1000 cm^-3), and possibly compressed and shock-excited. These winds have comparable mass loading factors as the SF-driven winds but carry ~10 (~50) times more momentum (energy). The results confirm our previous findings of high duty cycle, energy-driven outflows powered by AGN above the Schechter mass, which may contribute to star formation quenching. △ Less

Submitted 26 February, 2019; v1 submitted 12 July, 2018; originally announced July 2018.

Comments: Accepted for publication by the The Astrophysical Journal. 36 pages, 13 Figures

Modeling the atomic-to-molecular transition in cosmological simulations of galaxy formation

Authors: Benedikt Diemer, Adam R. H. Stevens, John C. Forbes, Federico Marinacci, Lars Hernquist, Claudia del P. Lagos, Amiel Sternberg, Annalisa Pillepich, Dylan Nelson, Gergö Popping, Francisco Villaescusa-Navarro, Paul Torrey, Mark Vogelsberger

Abstract: Large-scale cosmological simulations of galaxy formation currently do not resolve the densities at which molecular hydrogen forms, implying that the atomic-to-molecular transition must be modeled either on the fly or in postprocessing. We present an improved postprocessing framework to estimate the abundance of atomic and molecular hydrogen and apply it to the IllustrisTNG simulations. We compare… ▽ More Large-scale cosmological simulations of galaxy formation currently do not resolve the densities at which molecular hydrogen forms, implying that the atomic-to-molecular transition must be modeled either on the fly or in postprocessing. We present an improved postprocessing framework to estimate the abundance of atomic and molecular hydrogen and apply it to the IllustrisTNG simulations. We compare five different models for the atomic-to-molecular transition, including empirical, simulation-based, and theoretical prescriptions. Most of these models rely on the surface density of neutral hydrogen and the ultraviolet (UV) flux in the Lyman-Werner band as input parameters. Computing these quantities on the kiloparsec scales resolved by the simulations emerges as the main challenge. We show that the commonly used Jeans length approximation to the column density of a system can be biased and exhibits large cell-to-cell scatter. Instead, we propose to compute all surface quantities in face-on projections and perform the modeling in two dimensions. In general, the two methods agree on average, but their predictions diverge for individual galaxies and for models based on the observed midplane pressure of galaxies. We model the UV radiation from young stars by assuming a constant escape fraction and optically thin propagation throughout the galaxy. With these improvements, we find that the five models for the atomic-to-molecular transition roughly agree on average but that the details of the modeling matter for individual galaxies and the spatial distribution of molecular hydrogen. We emphasize that the estimated molecular fractions are approximate due to the significant systematic uncertainties. △ Less

Submitted 23 October, 2018; v1 submitted 6 June, 2018; originally announced June 2018.

Comments: 22 pages, 13 figures

Journal ref: 2018 ApJS 238, 33