Showing 1–50 of 196 results for author: Su, K

Unraveling the role of merger histories in the population of Insitu stars: linking IllustrisTNG cosmological simulation to H3 survey

Authors: Razieh Emami, Lars Hernquist, Randall Smith, James F. Steiner, Grant Tremblay, Douglas Finkbeiner, Mark Vogelsberger, Josh Grindlay, Federico Marinacci, Kung-Yi Su, Cecilia Garraffo, Yuan-Sen Ting, Phillip A. Cargile, Rebecca L. Davies, Chloë E. Benton, Yijia Li, Letizia Bugiani, Amir H. Khoram, Sownak Bose

Abstract: We undertake a comprehensive investigation into the distribution of insitu stars within Milky Way-like galaxies, leveraging TNG50 simulations and comparing their predictions with data from the H3 survey. Our analysis reveals that 28% of galaxies demonstrate reasonable agreement with H3, while only 12% exhibit excellent alignment in their profiles, regardless of the specific spatial cut employed to… ▽ More We undertake a comprehensive investigation into the distribution of insitu stars within Milky Way-like galaxies, leveraging TNG50 simulations and comparing their predictions with data from the H3 survey. Our analysis reveals that 28% of galaxies demonstrate reasonable agreement with H3, while only 12% exhibit excellent alignment in their profiles, regardless of the specific spatial cut employed to define insitu stars. To uncover the underlying factors contributing to deviations between TNG50 and H3 distributions, we scrutinize correlation coefficients among internal drivers(e.g., virial radius, star formation rate [SFR]) and merger-related parameters (such as the effective mass-ratio, mean distance, average redshift, total number of mergers, average spin-ratio and maximum spin alignment between merging galaxies). Notably, we identify significant correlations between deviations from observational data and key parameters such as the median slope of virial radius, mean SFR values, and the rate of SFR change across different redshift scans. Furthermore, positive correlations emerge between deviations from observational data and parameters related to galaxy mergers. We validate these correlations using the Random Forest Regression method. Our findings underscore the invaluable insights provided by the H3 survey in unraveling the cosmic history of galaxies akin to the Milky Way, thereby advancing our understanding of galactic evolution and shedding light on the formation and evolution of Milky Way-like galaxies in cosmological simulations. △ Less

Submitted 9 July, 2024; originally announced July 2024.

Comments: 20 pages, 13 figures

Multi-Zone Modeling of Black Hole Accretion and Feedback in 3D GRMHD: Bridging Vast Spatial and Temporal Scales

Authors: Hyerin Cho, Ben S. Prather, Kung-Yi Su, Ramesh Narayan, Priyamvada Natarajan

Abstract: Simulating accretion and feedback from the horizon scale of supermassive black holes (SMBHs) out to galactic scales is challenging because of the vast range of scales involved. We describe and test a "multi-zone" technique which is designed to tackle this difficult problem in 3D general relativistic magnetohydrodynamic (GRMHD) simulations. We simulate accretion on a non-spinning SMBH ($a_*=0$) usi… ▽ More Simulating accretion and feedback from the horizon scale of supermassive black holes (SMBHs) out to galactic scales is challenging because of the vast range of scales involved. We describe and test a "multi-zone" technique which is designed to tackle this difficult problem in 3D general relativistic magnetohydrodynamic (GRMHD) simulations. We simulate accretion on a non-spinning SMBH ($a_*=0$) using initial conditions from a large scale galaxy simulation, and achieve steady state over 8 decades in radius. The density scales with radius as $ρ\propto r^{-1}$ inside the Bondi radius $R_B$, which is located at $R_B=2\times 10^5 \,r_g$ ($\approx 60\,{\rm pc}$ for M87) where $r_g$ is the gravitational radius of the SMBH; the plasma-$β\sim$ unity, indicating an extended magnetically arrested state; the mass accretion rate $\dot{M}$ is $\approx 1\%$ of the analytical Bondi accretion rate $\dot{M}_B$; and there is continuous energy feedback out to $\approx 100R_B$ (or beyond $>\,{\rm kpc}$) at a rate $\approx 0.02 \dot{M}c^2$. Surprisingly, any ordered rotation in the external medium does not survive as the magnetized gas flows to smaller radii, and the final steady solution is very similar to when the exterior has no rotation. Using the multi-zone method, we simulate GRMHD accretion for a wide range of Bondi radii, $R_{\rm B} \sim 10^2 - 10^7\,r_{\rm g}$, and find that $\dot{M}/\dot{M}_B\approx (R_B/6\, r_g)^{-0.5}$. △ Less

Submitted 22 May, 2024; originally announced May 2024.

Comments: 19 pages, 9 figures, submitted to ApJ

The ALMA-QUARKS Survey: II. the ACA 1.3 mm continuum source catalog and the assembly of dense gas in massive star-forming clumps

Authors: Fengwei Xu, Ke Wang, Tie Liu, Lei Zhu, Guido Garay, Xunchuan Liu, Paul Goldsmith, Qizhou Zhang, Patricio Sanhueza, Shengli Qin, Jinhua He, Mika Juvela, Anandmayee Tej, Hongli Liu, Shanghuo Li, Kaho Morii, Siju Zhang, Jianwen Zhou, Amelia Stutz, Neal J. Evans, Kim Kee-Tae, Shengyuan Liu, Diego Mardones, Guangxing Li, Leonardo Bronfman , et al. (8 additional authors not shown)

Abstract: Leveraging the high resolution, high sensitivity, and wide frequency coverage of the Atacama Large Millimeter/submillimeter Array (ALMA), the QUARKS survey, standing for "Querying Underlying mechanisms of massive star formation with ALMA-Resolved gas Kinematics and Structures", is observing 139 massive star-forming clumps at ALMA Band 6 ($λ\sim$ 1.3 mm). This paper introduces the Atacama Compact A… ▽ More Leveraging the high resolution, high sensitivity, and wide frequency coverage of the Atacama Large Millimeter/submillimeter Array (ALMA), the QUARKS survey, standing for "Querying Underlying mechanisms of massive star formation with ALMA-Resolved gas Kinematics and Structures", is observing 139 massive star-forming clumps at ALMA Band 6 ($λ\sim$ 1.3 mm). This paper introduces the Atacama Compact Array (ACA) 7-m data. Combining multi-wavelength data, we provide the first edition of QUARKS atlas, offering insights into the multiscale and multiphase interstellar medium in high-mass star formation. The ACA 1.3 mm catalog includes 207 continuum sources that are called ACA sources. Their gas kinetic temperatures are estimated using three formaldehyde (H$_2$CO) transitions with a non-LTE radiation transfer model, and the mass and density are derived from a dust emission model. The ACA sources are massive (16-84 percentile values of 6-160 $M_{\odot}$), gravity-dominated ($M\propto R^{1.1}$) fragments within massive clumps, with supersonic turbulence ($\mathcal{M}>1$) and embedded star-forming protoclusters. We find a linear correlation between the masses of the fragments and the massive clumps, with a ratio of 6% between the two. When considering the fragments as representative of dense gas, the ratio indicates a dense gas fraction (DGF) of 6%, although with a wide scatter ranging from 1% to 10%. If we consider the QUARKS massive clumps to be what is observed at various scales, then the size-independent DGF indicates a self-similar fragmentation or collapsing mode in protocluster formation. With the ACA data over four orders of magnitude of luminosity-to-mass ratio ($L/M$), we find that the DGF increases significantly with $L/M$, which indicates clump evolutionary stage. We observed a limited fragmentation at the subclump scale, which can be explained by dynamic global collapse process. △ Less

Submitted 4 April, 2024; v1 submitted 2 April, 2024; originally announced April 2024.

Comments: 24 pages, 7 figures. Accepted for publication in Research in Astronomy and Astrophysics. QUARKS atlas link: https://drive.google.com/file/d/1KTqXxCDduYepvLd9kIvZVSSytK48OmfL/view?usp=sharing

Abundant sub-micron grains revealed in newly discovered extreme debris discs

Authors: A. Moór, P. Ábrahám, K. Y. L. Su, T. Henning, S. Marino, L. Chen, Á. Kóspál, N. Pawellek, J. Varga, K. Vida

Abstract: Extreme debris discs (EDDs) are bright and warm circumstellar dusty structures around main sequence stars. They may represent the outcome of giant collisions occuring in the terrestrial region between large planetesimals or planetary bodies, and thus provide a rare opportunity to peer into the aftermaths of these events. Here, we report on results of a mini-survey we conducted with the aim to incr… ▽ More Extreme debris discs (EDDs) are bright and warm circumstellar dusty structures around main sequence stars. They may represent the outcome of giant collisions occuring in the terrestrial region between large planetesimals or planetary bodies, and thus provide a rare opportunity to peer into the aftermaths of these events. Here, we report on results of a mini-survey we conducted with the aim to increase the number of known EDDs, investigate the presence of solid-state features around 10μm in eight EDDs, and classify them into the silica or silicate dominated groups. We identify four new EDDs and derive their fundamental properties. For these, and for four other previously known discs, we study the spectral energy distribution around 10μm by means of VLT/VISIR photometry in three narrow-band filters and conclude that all eight objects likely exhibit solid-state emission features from sub-micron grains. We find that four discs probably belong to the silicate dominated subgroup. Considering the age distribution of the entire EDD sample, we find that their incidence begins to decrease only after 300 Myr, suggesting that the earlier common picture that these objects are related to the formation of rocky planets may not be exclusive, and that other processes may be involved for older objects (>100 Myr). Because most of the older EDD systems have wide, eccentric companions, we suggest that binarity may play a role in triggering late giant collisions. △ Less

Submitted 23 February, 2024; originally announced February 2024.

Comments: 20 pages, 9 figures, published in MNRAS

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 528, Issue 3, March 2024, pp.4528-4546

Observational Signatures of AGN Feedback in the Morphology and the Ionization States of Milky Way-like Galaxies

Authors: Nadia Qutob, Razieh Emami, Kung-Yi Su, Randall Smith, Lars Hernquist, Dian P. Triani, Cameron Hummels, Drummond Fielding, Philip F. Hopkins, Rachel S. Somerville, David R. Ballantyne, Mark Vogelsberger, Grant Tremblay, James F. Steiner, Douglas Finkbeiner, Ramesh Narayan, Minjung Park, Josh Grindlay, Priyamvada Natarajan, Christopher C. Hayward, Dušan Kereš, Sam B. Ponnada, Sirio Belli, Rebecca Davies, Gabriel Maheson , et al. (2 additional authors not shown)

Abstract: We make an in-depth analysis of different AGN jet models' signatures, inducing quiescence in galaxies with a halo mass of $10^{12} M_\odot$. Three jet models, including cosmic ray-dominant, hot thermal, and precessing kinetic jets, are studied at two energy flux levels each, compared to a jet-free, stellar feedback-only simulation. We examine the distribution of Mg II, O VI, and O VIII ions, along… ▽ More We make an in-depth analysis of different AGN jet models' signatures, inducing quiescence in galaxies with a halo mass of $10^{12} M_\odot$. Three jet models, including cosmic ray-dominant, hot thermal, and precessing kinetic jets, are studied at two energy flux levels each, compared to a jet-free, stellar feedback-only simulation. We examine the distribution of Mg II, O VI, and O VIII ions, alongside gas temperature and density profiles. Low-energy ions, like Mg II, concentrate in the ISM, while higher energy ions, e.g., O VIII, prevail at the AGN jet cocoon's edge. High-energy flux jets display an isotropic ion distribution with lower overall density. High-energy thermal or cosmic ray jets pressurize at smaller radii, significantly suppressing core density. The cosmic ray jet provides extra pressure support, extending cool and warm gas distribution. A break in the ion-to-mass ratio slope in O VI and O VIII is demonstrated in the ISM-to-CGM transition (between 10-30 kpc), growing smoothly towards the CGM at greater distances. △ Less

Submitted 22 December, 2023; originally announced December 2023.

Comments: 17 pages, 13 figures

Dense stellar clump formation driven by strong quasar winds in the FIRE cosmological hydrodynamic simulations

Authors: Jonathan Mercedes-Feliz, Daniel Anglés-Alcázar, Boon Kiat Oh, Christopher C. Hayward, Rachel K. Cochrane, Alexander J. Richings, Claude-André Faucher-Giguère, Sarah Wellons, Bryan A. Terrazas, Jorge Moreno, Kung Yi Su, Philip F. Hopkins

Abstract: We investigate the formation of dense stellar clumps in a suite of high-resolution cosmological zoom-in simulations of a massive, star forming galaxy at $z \sim 2$ under the presence of strong quasar winds. Our simulations include multi-phase ISM physics from the Feedback In Realistic Environments (FIRE) project and a novel implementation of hyper-refined accretion disk winds. We show that powerfu… ▽ More We investigate the formation of dense stellar clumps in a suite of high-resolution cosmological zoom-in simulations of a massive, star forming galaxy at $z \sim 2$ under the presence of strong quasar winds. Our simulations include multi-phase ISM physics from the Feedback In Realistic Environments (FIRE) project and a novel implementation of hyper-refined accretion disk winds. We show that powerful quasar winds can have a global negative impact on galaxy growth while in the strongest cases triggering the formation of an off-center clump with stellar mass ${\rm M}_{\star}\sim 10^{7}\,{\rm M}_{\odot}$, effective radius ${\rm R}_{\rm 1/2\,\rm Clump}\sim 20\,{\rm pc}$, and surface density $Σ_{\star} \sim 10^{4}\,{\rm M}_{\odot}\,{\rm pc}^{-2}$. The clump progenitor gas cloud is originally not star-forming, but strong ram pressure gradients driven by the quasar winds (orders of magnitude stronger than experienced in the absence of winds) lead to rapid compression and subsequent conversion of gas into stars at densities much higher than the average density of star-forming gas. The AGN-triggered star-forming clump reaches ${\rm SFR} \sim 50\,{\rm M}_{\odot}\,{\rm yr}^{-1}$ and $Σ_{\rm SFR} \sim 10^{4}\,{\rm M}_{\odot}\,{\rm yr}^{-1}\,{\rm kpc}^{-2}$, converting most of the progenitor gas cloud into stars in $\sim$2\,Myr, significantly faster than its initial free-fall time and with stellar feedback unable to stop star formation. In contrast, the same gas cloud in the absence of quasar winds forms stars over a much longer period of time ($\sim$35\,Myr), at lower densities, and losing spatial coherency. The presence of young, ultra-dense, gravitationally bound stellar clumps in recently quenched galaxies could thus indicate local positive feedback acting alongside the strong negative impact of powerful quasar winds, providing a plausible formation scenario for globular clusters. △ Less

Submitted 17 April, 2024; v1 submitted 30 October, 2023; originally announced October 2023.

Comments: 14 pages, 12 figures. Accepted for publication in MNRAS; doi:10.1093/mnras/stae1021

Bridging Scales in Black Hole Accretion and Feedback: Magnetized Bondi Accretion in 3D GRMHD

Authors: Hyerin Cho, Ben S. Prather, Ramesh Narayan, Priyamvada Natarajan, Kung-Yi Su, Angelo Ricarte, Koushik Chatterjee

Abstract: Fueling and feedback couple supermassive black holes (SMBHs) to their host galaxies across many orders of magnitude in spatial and temporal scales, making this problem notoriously challenging to simulate. We use a multi-zone computational method based on the general relativistic magneto-hydrodynamic (GRMHD) code KHARMA that allows us to span $7$ orders of magnitude in spatial scale, to simulate ac… ▽ More Fueling and feedback couple supermassive black holes (SMBHs) to their host galaxies across many orders of magnitude in spatial and temporal scales, making this problem notoriously challenging to simulate. We use a multi-zone computational method based on the general relativistic magneto-hydrodynamic (GRMHD) code KHARMA that allows us to span $7$ orders of magnitude in spatial scale, to simulate accretion onto a non-spinning SMBH from an external medium with Bondi radius $R_B\approx 2\times 10^5 \,GM_\bullet/c^2$, where $M_\bullet$ is the SMBH mass. For the classic idealized Bondi problem, spherical gas accretion without magnetic fields, our simulation results agree very well with the general relativistic analytic solution. Meanwhile, when the accreting gas is magnetized, the SMBH magnetosphere becomes saturated with a strong magnetic field. The density profile varies as $\sim r^{-1}$ rather than $r^{-3/2}$ and the accretion rate $\dot{M}$ is consequently suppressed by over 2 orders of magnitude below the Bondi rate $\dot{M}_B$. We find continuous energy feedback from the accretion flow to the external medium at a level of $\sim10^{-2}\dot{M}c^2 \sim 5\times 10^{-5} \dot{M}_B c^2$. Energy transport across these widely disparate scales occurs via turbulent convection triggered by magnetic field reconnection near the SMBH. Thus, strong magnetic fields that accumulate on horizon scales transform the flow dynamics far from the SMBH and naturally explain observed extremely low accretion rates compared to the Bondi rate, as well as at least part of the energy feedback. △ Less

Submitted 1 December, 2023; v1 submitted 29 October, 2023; originally announced October 2023.

Comments: 8 pages with 5 pages of Appendices, 4 figures, accepted to ApJ Letters

Unraveling Jet Quenching Criteria Across L* Galaxies and Massive Cluster Ellipticals

Authors: Kung-Yi Su, Greg L. Bryan, Christopher C. Hayward, Rachel S. Somerville, Philip F. Hopkins, Razieh Emami, Claude-André Faucher-Giguère, Eliot Quataert, Sam B. Ponnada, Drummond Fielding, Dušan Kereš

Abstract: In the absence of supplementary heat, the radiative cooling of halo gas around massive galaxies (Milky Way mass and above) leads to an excess of cold gas or stars beyond observed levels. AGN jet-induced heating is likely essential, but the specific properties of the jets remain unclear. Our previous work (Su et al. 2021) concludes from simulations of a halo with $10^{14} M_\odot$ that a successful… ▽ More In the absence of supplementary heat, the radiative cooling of halo gas around massive galaxies (Milky Way mass and above) leads to an excess of cold gas or stars beyond observed levels. AGN jet-induced heating is likely essential, but the specific properties of the jets remain unclear. Our previous work (Su et al. 2021) concludes from simulations of a halo with $10^{14} M_\odot$ that a successful jet model should have an energy flux comparable to the free-fall energy flux at the cooling radius and should inflate a sufficiently wide cocoon with a long enough cooling time. In this paper, we investigate three jet modes with constant fluxes satisfying the criteria, including high-temperature thermal jets, cosmic ray (CR)-dominant jets, and widely precessing kinetic jets in $10^{12}-10^{15}\,{\rm M}_{\odot}$ halos using high-resolution, non-cosmological MHD simulations with the FIRE-2 (Feedback In Realistic Environments) stellar feedback model, conduction, and viscosity. We find that scaling the jet energy according to the free-fall energy at the cooling radius can successfully suppress the cooling flows and quench galaxies without obviously violating observational constraints. We investigate an alternative scaling method in which we adjust the energy flux based on the total cooling rate within the cooling radius. However, we observe that the strong interstellar medium (ISM) cooling dominates the total cooling rate in this scaling approach, resulting in a jet flux that exceeds the amount needed to suppress the cooling flows. With the same energy flux, the CR-dominant jet is most effective in suppressing the cooling flow across all the surveyed halo masses due to the enhanced CR pressure support. We confirm that the criteria for a successful jet model, which we proposed in Su et al. (2021), work across a much wider range, encompassing halo masses of $10^{12}-10^{15} {\rm M_\odot}$. △ Less

Submitted 2 November, 2023; v1 submitted 26 October, 2023; originally announced October 2023.

Comments: 15 pages, 10 figures, Submitted to MNRAS; comments welcome!

Searching for Planets Orbiting Fomalhaut with JWST/NIRCam

Authors: Marie Ygouf, Charles Beichman, Jorge Llop-Sayson, Geoffrey Bryden, Jarron Leisenring, Andras Gaspar, John Krist, Marcia Rieke, George Rieke, Schuyler Wolff, Thomas Roellig, Kate Su, Kevin Hainline, Klaus Hodapp, Thomas Greene, Michael Meyer, Doug Kelly, Karl Misselt, John Stansberry, Martha Boyer, Doug Johnstone, Scott Horner, Alexandra Greenbaum

Abstract: We report observations with the JWST/NIRCam coronagraph of the Fomalhaut system. This nearby A star hosts a complex debris disk system discovered by the IRAS satellite. Observations in F444W and F356W filters using the round 430R mask achieve a contrast ratio of ~ 4 x 10-7 at 1'' and ~ 4 x 10-8 outside of 3''. These observations reach a sensitivity limit <1 MJup across most of the disk region. Con… ▽ More We report observations with the JWST/NIRCam coronagraph of the Fomalhaut system. This nearby A star hosts a complex debris disk system discovered by the IRAS satellite. Observations in F444W and F356W filters using the round 430R mask achieve a contrast ratio of ~ 4 x 10-7 at 1'' and ~ 4 x 10-8 outside of 3''. These observations reach a sensitivity limit <1 MJup across most of the disk region. Consistent with the hypothesis that Fomalhaut b is not a massive planet but is a dust cloud from a planetesimal collision, we do not detect it in either F356W or F444W (the latter band where a Jovian-sized planet should be bright). We have reliably detected 10 sources in and around Fomalhaut and its debris disk, all but one of which are coincident with Keck or HST sources seen in earlier coronagraphic imaging; we show them to be background objects, including the "Great Dust Cloud" identified in MIRI data. However, one of the objects, located at the edge of the inner dust disk seen in the MIRI images, has no obvious counterpart in imaging at earlier epochs and has a relatively red [F356W]-[F444W]>0.7 mag (Vega) color. Whether this object is a background galaxy, brown dwarf, or a Jovian mass planet in the Fomalhaut system will be determined by an approved Cycle 2 follow-up program. Finally, we set upper limits to any scattered light from the outer ring, placing a weak limit on the dust albedo at F356W and F444W. △ Less

Submitted 23 October, 2023; originally announced October 2023.

Comments: 24 pages, 17 figures

RZ Piscium Hosts a Compact and Highly Perturbed Debris Disk

Authors: Kate Y. L. Su, Grant M. Kennedy, George H. Rieke, A. Meredith Hughes, Yu-Chia Lin, Jamar Kittling, Alan P. Jackson, Ramya M. Anche, Hauyu Baobab Liu

Abstract: RZ Piscium (RZ Psc) is well-known in the variable star field because of its numerous, irregular optical dips in the past five decades, but the nature of the system is heavily debated in the literature. We present multiyear infrared monitoring data from Spitzer and WISE to track the activities of the inner debris production, revealing stochastic infrared variability as short as weekly timescales th… ▽ More RZ Piscium (RZ Psc) is well-known in the variable star field because of its numerous, irregular optical dips in the past five decades, but the nature of the system is heavily debated in the literature. We present multiyear infrared monitoring data from Spitzer and WISE to track the activities of the inner debris production, revealing stochastic infrared variability as short as weekly timescales that is consistent with destroying a 90-km-size asteroid every year. ALMA 1.3 mm data combined with spectral energy distribution modeling show that the disk is compact ($\sim$0.1--13 au radially) and lacks cold gas. The disk is found to be highly inclined and has a significant vertical scale height. These observations confirm that RZ Psc hosts a close to edge-on, highly perturbed debris disk possibly due to migration of recently formed giant planets which might be triggered by the low-mass companion RZ Psc B if the planets formed well beyond the snowlines. △ Less

Submitted 18 October, 2023; originally announced October 2023.

Comments: 16 pages, 5 figures, accepted for publication in ApJ

The first white dwarf debris disk observed by JWST

Authors: Andrew Swan, Jay Farihi, Kate Y. L. Su, Steven J. Desch

Abstract: This letter reports the first JWST spectroscopy of a white dwarf debris disk, giving a preliminary assessment of the salient features, and recommendations for future observations. The polluted and dusty star WD 0145+234 experienced a major collisional event in its circumstellar disk in 2018, accompanied by an infrared outburst, and subsequently a gradual decrease in thermal emission. Time-series N… ▽ More This letter reports the first JWST spectroscopy of a white dwarf debris disk, giving a preliminary assessment of the salient features, and recommendations for future observations. The polluted and dusty star WD 0145+234 experienced a major collisional event in its circumstellar disk in 2018, accompanied by an infrared outburst, and subsequently a gradual decrease in thermal emission. Time-series NIRSPEC observations demonstrate that the circumstellar disk is returning to a quiescent state with a T~1000 K infrared excess similar to the bulk of known dusty white dwarfs. MIRI spectroscopy reveals a 9-12 micron solid-state emission feature consistent with silicate minerals as observed in debris disks observed with Spitzer IRS. The strength and morphology of the silicate feature appear unchanged relative to the continuum in spectra taken over a year apart, consistent with steady-state collisional evolution of the circumstellar debris. A tentative emission feature around 7 microns may be due to carbonates, and if confirmed would indicate aqueous alteration in the parent body. △ Less

Submitted 19 December, 2023; v1 submitted 13 October, 2023; originally announced October 2023.

Comments: Accepted to MNRAS Letters

An Analytic Model For Magnetically-Dominated Accretion Disks

Authors: Philip F. Hopkins, Jonathan Squire, Eliot Quataert, Norman Murray, Kung-Yi Su, Ulrich P. Steinwandel, Kyle Kremer, Claude-Andre Faucher-Giguere, Sarah Wellons

Abstract: Recent numerical cosmological radiation-magnetohydrodynamic-thermochemical-star formation simulations have resolved the formation of quasar accretion disks with Eddington or super-Eddington accretion rates onto supermassive black holes (SMBHs) down to a few hundred gravitational radii. These 'flux-frozen' and hyper-magnetized disks appear to be qualitatively distinct from classical $α$ disks and m… ▽ More Recent numerical cosmological radiation-magnetohydrodynamic-thermochemical-star formation simulations have resolved the formation of quasar accretion disks with Eddington or super-Eddington accretion rates onto supermassive black holes (SMBHs) down to a few hundred gravitational radii. These 'flux-frozen' and hyper-magnetized disks appear to be qualitatively distinct from classical $α$ disks and magnetically-arrested disks: the midplane pressure is dominated by toroidal magnetic fields with plasma $β\ll 1$ powered by advection of magnetic flux from the interstellar medium (ISM), and they are super-sonically and trans-Alfvenically turbulent with cooling times short compared to dynamical times yet remain gravitationally stable owing to magnetic support. In this paper, we present a simple analytic similarity model for such disks. For reasonable assumptions, the model is entirely specified by the boundary conditions (inflow rate at the BH radius of influence [BHROI]). We show that the scalings from this model are robust to various detailed assumptions, agree remarkably well with the simulations (given their simplicity), and demonstrate the self-consistency and gravitational stability of such disks even in the outer accretion disk (approaching the BHROI) at hyper-Eddington accretion rates. △ Less

Submitted 12 March, 2024; v1 submitted 6 October, 2023; originally announced October 2023.

Comments: 12 pages, 1 figure. Replaced with accepted version. Companion paper to 'FORGE'D IN FIRE II' (arXiv:2310.04506, part of a series with arXiv:2309.13115 -- animations of the simulations referred to here can be viewed at http://www.tapir.caltech.edu/~phopkins/Site/animations/Movies_zoom.html)

FORGE'd in FIRE II: The Formation of Magnetically-Dominated Quasar Accretion Disks from Cosmological Initial Conditions

Authors: Philip F. Hopkins, Jonathan Squire, Kung-Yi Su, Ulrich P. Steinwandel, Kyle Kremer, Yanlong Shi, Michael Y. Grudic, Sarah Wellons, Claude-Andre Faucher-Giguere, Daniel Angles-Alcazar, Norman Murray, Eliot Quataert

Abstract: In a companion paper, we reported the self-consistent formation of quasar accretion disks with inflow rates $\sim 10\,{\rm M_{\odot}\,yr^{-1}}$ down to <300 Schwarzschild radii from cosmological radiation-magneto-thermochemical-hydrodynamical galaxy and star formation simulations. We see the formation of a well-defined, steady-state accretion disk which is stable against star formation at sub-pc s… ▽ More In a companion paper, we reported the self-consistent formation of quasar accretion disks with inflow rates $\sim 10\,{\rm M_{\odot}\,yr^{-1}}$ down to <300 Schwarzschild radii from cosmological radiation-magneto-thermochemical-hydrodynamical galaxy and star formation simulations. We see the formation of a well-defined, steady-state accretion disk which is stable against star formation at sub-pc scales. The disks are optically thick, with radiative cooling balancing accretion, but with properties that are distinct from those assumed in most previous accretion disk models. The pressure is strongly dominated by (primarily toroidal) magnetic fields, with a plasma $β\sim 10^{-4}$ even in the disk midplane. They are qualitatively distinct from magnetically elevated or arrested disks. The disks are strongly turbulent, with trans-Alfvenic and highly super-sonic turbulence, and balance this via a cooling time that is short compared to the disk dynamical time, and can sustain highly super-Eddington accretion rates. Their surface and 3D densities at $\sim 10^{3}-10^{5}$ gravitational radii are much lower than in a Shakura-Sunyaev disk, with important implications for their thermo-chemistry and stability. We show how the magnetic field strengths and geometries arise from rapid advection of flux with the inflow from much weaker galaxy-scale fields in these 'flux-frozen' disks, and how this stabilizes the disk and gives rise to efficient torques. Re-simulating without magnetic fields produces catastrophic fragmentation with a vastly smaller, lower-$\dot{M}$ Shakura-Sunyaev-like disk. △ Less

Submitted 18 January, 2024; v1 submitted 6 October, 2023; originally announced October 2023.

Comments: 46 pages, 30 figures. Updated to match accepted version from The Open Journal of Astrophysics. Part of a series with arXiv:2309.13115 . Animations of the simulations can be viewed at http://www.tapir.caltech.edu/~phopkins/Site/animations/Movies_zoom.html

Seeding the CGM: How Satellites Populate the Cold Phase of Milky Way Halos

Authors: Manami Roy, Kung-Yi Su, Stephanie Tonnesen, Drummond B. Fielding, Claude-André Faucher-Giguère

Abstract: The origin of the cold phase in the CGM is a highly debated question. We investigate the contribution of satellite galaxies to the cold gas budget in the circumgalactic medium (CGM)of a Milky Way-like host galaxy. We perform controlled experiments with three different satellite mass distributions and identify several mechanisms by which satellites can add cold gas to the CGM, including ram pressur… ▽ More The origin of the cold phase in the CGM is a highly debated question. We investigate the contribution of satellite galaxies to the cold gas budget in the circumgalactic medium (CGM)of a Milky Way-like host galaxy. We perform controlled experiments with three different satellite mass distributions and identify several mechanisms by which satellites can add cold gas to the CGM, including ram pressure stripping and induced cooling in the mixing layer of the stripped cold gas. These two mechanisms contribute a comparable amount of cold gas to the host CGM. We find that the less massive satellites ($\leq 10^9 M_\odot$) not only lose all of their cold gas in a short period ($\sim$ 0.5-1 Gyr), but their stripped cold clouds also mix with the hot CGM gas and get heated up quickly. However, stellar feedback from these less massive satellites can hugely alter the fate of their stripped gas. Feedback speeds up the destruction of the stripped cold clouds from these satellites by making them more diffuse with more surface area. On the other hand, the more massive satellites (LMC or SMC-like $\sim 10^{10} M_\odot$) can add cold gas to the total gas budget of the host CGM for several Gyrs. △ Less

Submitted 6 October, 2023; originally announced October 2023.

Comments: 17 pages, 12 figures, 2 tables Accepted for publication in MNRAS on Oct 6, 2023

Synchrotron Signatures of Cosmic Ray Transport Physics in Galaxies

Authors: Sam B. Ponnada, Iryna S. Butsky, Raphael Skalidis, Philip F. Hopkins, Georgia V. Panopoulou, Cameron Hummels, Dušan Kereš, Eliot Quataert, Claude-André Faucher-Giguère, Kung-Yi Su

Abstract: Cosmic rays (CRs) may drive outflows and alter the phase structure of the circumgalactic medium, with potentially important implications on galaxy formation. However, these effects ultimately depend on the dominant mode of transport of CRs within and around galaxies, which remains highly uncertain. To explore potential observable constraints on CR transport, we investigate a set of cosmological FI… ▽ More Cosmic rays (CRs) may drive outflows and alter the phase structure of the circumgalactic medium, with potentially important implications on galaxy formation. However, these effects ultimately depend on the dominant mode of transport of CRs within and around galaxies, which remains highly uncertain. To explore potential observable constraints on CR transport, we investigate a set of cosmological FIRE-2 CR-MHD simulations of L$_{\ast}$ galaxies which evolve CRs with transport models motivated by self-confinement (SC) and extrinsic turbulence (ET) paradigms. To first order, the synchrotron properties diverge between SC and ET models due to a CR physics driven hysteresis. SC models show a higher tendency to undergo `ejective' feedback events due to a runaway buildup of CR pressure in dense gas due to the behavior of SC transport scalings at extremal CR energy densities. The corresponding CR wind-driven hysteresis results in brighter, smoother, and more extended synchrotron emission in SC runs relative to ET and constant diffusion runs. The differences in synchrotron arise from different morphology, ISM gas and \textbf{B} properties, potentially ruling out SC as the dominant mode of CR transport in typical star-forming L$_{\ast}$ galaxies, and indicating the potential for non-thermal radio continuum observations to constrain CR transport physics. △ Less

Submitted 28 February, 2024; v1 submitted 28 September, 2023; originally announced September 2023.

Comments: 5 pages, 3 figures. Accepted version

The ALMA Survey of Star Formation and Evolution in Massive Protoclusters with Blue Profiles (ASSEMBLE): Core Growth, Cluster Contraction, and Primordial Mass Segregation

Authors: Fengwei Xu, Ke Wang, Tie Liu, Mengyao Tang, Neal J. Evans II, Aina Palau, Kaho Morii, Jinhua He, Patricio Sanhueza, Hong-Li Liu, Amelia Stutz, Qizhou Zhang, Xi Chen, Pak Shing Li, Gilberto C. Gómez, Enrique Vázquez-Semadeni, Shanghuo Li, Xiaofeng Mai, Xing Lu, Meizhu Liu, Li Chen, Chuanshou Li, Hongqiong Shi, Zhiyuan Ren, Di Li , et al. (18 additional authors not shown)

Abstract: The ALMA Survey of Star Formation and Evolution in Massive Protoclusters with Blue Profiles (ASSEMBLE) aims to investigate the process of mass assembly and its connection to high-mass star formation theories in protoclusters in a dynamic view. We observed 11 massive (Mclump>1000 Msun), luminous (Lbol>10,000 Lsun), and blue-profile (infall signature) clumps by ALMA with resolution of 2200-5500 au a… ▽ More The ALMA Survey of Star Formation and Evolution in Massive Protoclusters with Blue Profiles (ASSEMBLE) aims to investigate the process of mass assembly and its connection to high-mass star formation theories in protoclusters in a dynamic view. We observed 11 massive (Mclump>1000 Msun), luminous (Lbol>10,000 Lsun), and blue-profile (infall signature) clumps by ALMA with resolution of 2200-5500 au at 350 GHz (870 um) in continuum and line emission. 248 dense cores were identified, including 106 cores showing protostellar signatures and 142 prestellar core candidates. Compared to early-stage infrared dark clouds (IRDCs) by ASHES, the core mass and surface density within the ASSEMBLE clumps exhibited significant increment, suggesting concurrent core accretion during the evolution of the clumps. The maximum mass of prestellar cores was found to be 2 times larger than that in IRDCs, indicating evolved protoclusters have the potential to harbor massive prestellar cores. The mass relation between clumps and their most massive core (MMCs) is observed in ASSEMBLE but not in IRDCs, which is suggested to be regulated by multiscale mass accretion. The mass correlation between the core clusters and their MMCs has a steeper slope compared to that observed in stellar clusters, which can be due to fragmentation of the MMC and stellar multiplicity. We observe a decrease in core separation and an increase in central concentration as protoclusters evolve. We confirm primordial mass segregation in the ASSEMBLE protoclusters, possibly resulting from gravitational concentration and/or gas accretion. △ Less

Submitted 26 September, 2023; originally announced September 2023.

Comments: 37 pages, 13 figures, 5 tables; accepted for publication in ApJS

FORGE'd in FIRE: Resolving the End of Star Formation and Structure of AGN Accretion Disks from Cosmological Initial Conditions

Authors: Philip F. Hopkins, Michael Y. Grudic, Kung-Yi Su, Sarah Wellons, Daniel Angles-Alcazar, Ulrich P. Steinwandel, David Guszejnov, Norman Murray, Claude-Andre Faucher-Giguere, Eliot Quataert, Dusan Keres

Abstract: It has recently become possible to zoom-in from cosmological to sub-pc scales in galaxy simulations to follow accretion onto supermassive black holes (SMBHs). However, at some point the approximations used on ISM scales (e.g. optically-thin cooling and stellar-population-integrated star formation [SF] and feedback [FB]) break down. We therefore present the first cosmological radiation-magnetohydro… ▽ More It has recently become possible to zoom-in from cosmological to sub-pc scales in galaxy simulations to follow accretion onto supermassive black holes (SMBHs). However, at some point the approximations used on ISM scales (e.g. optically-thin cooling and stellar-population-integrated star formation [SF] and feedback [FB]) break down. We therefore present the first cosmological radiation-magnetohydrodynamic (RMHD) simulation which self-consistently combines the FIRE physics (relevant on galactic/ISM scales where SF/FB are ensemble-averaged) and STARFORGE physics (relevant on small scales where we track individual (proto)stellar formation and evolution), together with explicit RMHD (including non-ideal MHD and multi-band M1-RHD) which self-consistently treats both optically-thick and thin regimes. This allows us to span scales from ~100 Mpc down to <100 au (~300 Schwarzschild radii) around a SMBH at a time where it accretes as a bright quasar, in a single simulation. We show that accretion rates up to $\sim 10-100\,{\rm M_{\odot}\,yr^{-1}}$ can be sustained into the accretion disk at $\ll 10^{3}\,R_{\rm schw}$, with gravitational torques between stars and gas dominating on sub-kpc scales until star formation is shut down on sub-pc scales by a combination of optical depth to cooling and strong magnetic fields. There is an intermediate-scale, flux-frozen disk which is gravitoturbulent and stabilized by magnetic pressure sustaining strong turbulence and inflow with persistent spiral modes. In this paper we focus on how gas gets into the small-scale disk, and how star formation is efficiently suppressed. △ Less

Submitted 12 March, 2024; v1 submitted 22 September, 2023; originally announced September 2023.

Comments: 38 pages, 18 figures. Replaced with accepted version

Synchrotron Emission on FIRE: Equipartition Estimators of Magnetic Fields in Simulated Galaxies with Spectrally-Resolved Cosmic Rays

Authors: Sam B. Ponnada, Georgia V. Panopoulou, Iryna S. Butsky, Philip F. Hopkins, Raphael Skalidis, Cameron Hummels, Eliot Quataert, Dušan Kereš, Claude-André Faucher-Giguère, Kung-Yi Su

Abstract: Synchrotron emission is one of few observable tracers of galactic magnetic fields (\textbf{B}) and cosmic rays (CRs). Much of our understanding of \textbf{B} in galaxies comes from utilizing synchrotron observations in conjunction with several simplifying assumptions of equipartition models, however it remains unclear how well these assumptions hold, and what \textbf{B} these estimates physically… ▽ More Synchrotron emission is one of few observable tracers of galactic magnetic fields (\textbf{B}) and cosmic rays (CRs). Much of our understanding of \textbf{B} in galaxies comes from utilizing synchrotron observations in conjunction with several simplifying assumptions of equipartition models, however it remains unclear how well these assumptions hold, and what \textbf{B} these estimates physically represent. Using FIRE simulations which self consistently evolve CR proton, electron, and positron spectra from MeV to TeV energies, we present the first synthetic synchrotron emission predictions from simulated L$_{*}$ galaxies with "live" spectrally-resolved CR-MHD. We find that synchrotron emission can be dominated by relatively cool and dense gas, resulting in equipartition estimates of \textbf{B} with fiducial assumptions underestimating the "true" \textbf{B} in the gas that contributes the most emission by factors of 2-3 due to small volume filling factors. Motivated by our results, we present an analytic framework that expands upon equipartition models for estimating \textbf{B} in a multi-phase medium. Comparing our spectrally-resolved synchrotron predictions to simpler spectral assumptions used in galaxy simulations with CRs, we find that spectral evolution can be crucial for accurate synchrotron calculations towards galactic centers, where loss terms are large. △ Less

Submitted 9 January, 2024; v1 submitted 8 September, 2023; originally announced September 2023.

Comments: 10 pages, 5 figures (with 3 additional figures in the Appendix). Accepted version (MNRAS December 21, 2023)

The science case for a far-infrared interferometer in the era of JWST and ALMA

Authors: David Leisawitz, Matteo Bonato, Duncan Farrah, T. Tupper Hyde, Aláine Lee, Joshua Bennett Lovell, Brenda Matthews, Lee G. Mundy, Conor Nixon, Petr Pokorny, Berke V. Ricketti, Giorgio Savini, Jeremy Scott, Irene Shivaei, Locke Spencer, Kate Su, C. Megan Urry, David Wilner

Abstract: A space-based far-infrared interferometer could work synergistically with the James Webb Space Telescope (JWST) and the Atacama Large Millimeter Array (ALMA) to revolutionize our understanding of the astrophysical processes leading to the formation of habitable planets and the co-evolution of galaxies and their central supermassive black holes. Key to these advances are measurements of water in it… ▽ More A space-based far-infrared interferometer could work synergistically with the James Webb Space Telescope (JWST) and the Atacama Large Millimeter Array (ALMA) to revolutionize our understanding of the astrophysical processes leading to the formation of habitable planets and the co-evolution of galaxies and their central supermassive black holes. Key to these advances are measurements of water in its frozen and gaseous states, observations of astronomical objects in the spectral range where most of their light is emitted, and access to critical diagnostic spectral lines, all of which point to the need for a far-infrared observatory in space. The objects of interest - circumstellar disks and distant galaxies - typically appear in the sky at sub-arcsecond scales, which rendered all but a few of them unresolvable with the successful and now-defunct 3.5-m Herschel Space Observatory, the largest far-infrared telescope flown to date. A far-infrared interferometer with maximum baseline length in the tens of meters would match the angular resolution of JWST at 10x longer wavelengths and observe water ice and water-vapor emission, which ALMA can barely do through the Earth's atmosphere. Such a facility was conceived and studied two decades ago. Here we revisit the science case for a space-based far-infrared interferometer in the era of JWST and ALMA and summarize the measurement capabilities that will enable the interferometer to achieve a set of compelling scientific objectives. Common to all the science themes we consider is a need for sub-arcsecond image resolution. △ Less

Submitted 29 August, 2023; originally announced August 2023.

Comments: 11 pages, 12 figures. To be published in Proc. SPIE Conf. 12686 "Instruments: Scientific Returns and Conceptual Designs"

V445 Puppis -- Dustier than a Thousand Novae

Authors: D. P. K. Banerjee, A. Evans, C. E. Woodward, S. Starrfield, K. Y. L. Su, N. M. Ashok, R. M. Wagner

Abstract: V445 Puppis, the only known Galactic helium nova, is a unique test-bed to verify supernova (SN) theories in the single degenerate channel that involve a white dwarf (WD) accreting matter from a helium-rich donor. An estimate of the mass of the helium shell on the WD is crucial to deciding whether or not it will undergo a SN detonation. In this context, this study estimates the dust and ejecta mass… ▽ More V445 Puppis, the only known Galactic helium nova, is a unique test-bed to verify supernova (SN) theories in the single degenerate channel that involve a white dwarf (WD) accreting matter from a helium-rich donor. An estimate of the mass of the helium shell on the WD is crucial to deciding whether or not it will undergo a SN detonation. In this context, this study estimates the dust and ejecta masses in the 2000 November eruption of V445 Pup. Subsequent to its outburst, the star became cocooned in a dust envelope. An analysis of the spectral energy distribution (SED) of the dust using infrared data shows that V445 Pup produced at least $10^{-3}$ solar masses of dust which is unprecedented for a classical or recurrent nova. The SED can be explained by a combination of a cold dust component at 105 +/- 10 K, mass (1.9 +/- 0.8) $\times 10^{-3}$ solar masses, and a warm dust component at 255 +/- 10 K, mass (2.2 +/- 1.2) $\times 10^{-5}$ solar masses. For a conservative choice of the gas-to-dust mass ratio in the range 10--100, the mass of the ejecta is 0.01--0.1 solar masses. Such a high mass range raises the question: why did V445 Pup not detonate as a Type 1a SN as is predicted in certain double-detonation sub-Chandrasekhar supernovae formalisms? We re-examine the nature of V445 Pup and discuss its role as a potential SN progenitor. △ Less

Submitted 18 June, 2023; originally announced June 2023.

Comments: 8 page, 3 figures, 1 table, Accepted to ApJ Letters

Accretion onto disk galaxies via hot and rotating CGM inflows

Authors: Jonathan Stern, Drummond Fielding, Zachary Hafen, Kung-Yi Su, Nadav Naor, Claude-André Faucher-Giguère, Eliot Quataert, James Bullock

Abstract: Observed accretion rates onto the Milky-Way and other local spirals fall short of that required to sustain star formation for cosmological timescales. A potential avenue for this unseen accretion is an inflow in the volume-filling hot phase ($\sim10^6$ K) of the circumgalactic medium (CGM), as suggested by some cosmological simulations. Using hydrodynamic simulations and a new analytic solution va… ▽ More Observed accretion rates onto the Milky-Way and other local spirals fall short of that required to sustain star formation for cosmological timescales. A potential avenue for this unseen accretion is an inflow in the volume-filling hot phase ($\sim10^6$ K) of the circumgalactic medium (CGM), as suggested by some cosmological simulations. Using hydrodynamic simulations and a new analytic solution valid in the slow-rotation limit, we show that a hot inflow spins up as it approaches the galaxy, while remaining hot, subsonic and quasi-spherical. At the radius of angular momentum support ($\approx15$ kpc for the Milky-Way) the hot flow flattens into a disk geometry and then cools from $\sim10^6$ K to $\sim10^4$ K at the disk-halo interface. Cooling affects all hot gas, rather than just a subset of individual gas clouds, implying that accretion via hot inflows does not rely on local thermal instability in contrast with 'precipitation' models for galaxy accretion. Prior to cooling and accretion the inflow completes $\sim t_{\rm cool}/t_{\rm ff}$ radians of rotation, where $t_{\rm cool}/t_{\rm ff}$ is the cooling time to free-fall time ratio in hot gas immediately outside the galaxy. The ratio $t_{\rm cool}/t_{\rm ff}$ may thus govern the development of turbulence and enhancement of magnetic fields in gas accreting onto low-redshift spirals. We argue that accretion via hot inflows can explain the observed truncation of nearby thin stellar disks at $\approx4$ disk radii. We also show that if rotating hot inflows are common in Milky-Way size disk galaxies, as predicted, then signatures should be observable with X-ray telescopes and FRB surveys. △ Less

Submitted 15 February, 2024; v1 submitted 31 May, 2023; originally announced June 2023.

Comments: 21 pages, 12 figures, resubmitted to MNRAS following referee report

Spatially resolved imaging of the inner Fomalhaut disk using JWST/MIRI

Authors: Andras Gaspar, Schuyler Grace Wolff, George H. Rieke, Jarron M. Leisenring, Jane Morrison, Kate Y. L. Su, Kimberly Ward-Duong, Jonathan Aguilar, Marie Ygouf, Charles Beichman, Jorge Llop-Sayson, Geoffrey Bryden

Abstract: Planetary debris disks around other stars are analogous to the Asteroid and Kuiper belts in the Solar System. Their structure reveals the configuration of small bodies and provides hints for the presence of planets. The nearby star Fomalhaut hosts one of the most prominent debris disks, resolved by HST, Spitzer, Herschel, and ALMA. Images of this system at mid-infrared wavelengths using JWST/MIRI… ▽ More Planetary debris disks around other stars are analogous to the Asteroid and Kuiper belts in the Solar System. Their structure reveals the configuration of small bodies and provides hints for the presence of planets. The nearby star Fomalhaut hosts one of the most prominent debris disks, resolved by HST, Spitzer, Herschel, and ALMA. Images of this system at mid-infrared wavelengths using JWST/MIRI not only show the narrow Kuiper-Belt-analog outer ring, but also that (1) what was thought from indirect evidence to be an asteroid-analog structure is instead broad, extending outward into the outer system; (2) there is an intermediate belt, probably shepherded by an unseen planet. The newly discovered belt is demarcated by an inner gap, located at ~ 78 au, and it is misaligned relative to the outer belt. The previously known collisionally generated dust cloud, Fomalhaut b, could have originated from this belt, suggesting increased dynamical stirring and collision rates there. We also discovered a large dust cloud within the outer ring, possible evidence of another dust-creating collision. Taken together with previous observations, Fomalhaut appears to be the site of a complex and possibly dynamically active planetary system. △ Less

Submitted 5 May, 2023; originally announced May 2023.

Comments: 17 Figures, 3 tables; includes Supplementary Information package; Published in Nature Astronomy

The impact of AGN-driven winds on physical and observable galaxy sizes

Authors: R. K. Cochrane, D. Anglés-Alcázar, J. Mercedes-Feliz, C. C. Hayward, C. -A. Faucher-Giguère, S. Wellons, B. A. Terrazas, A. Wetzel, P. F. Hopkins, J. Moreno, K. -Y. Su, R. S. Somerville

Abstract: Without AGN feedback, simulated massive, star-forming galaxies become too compact relative to observed galaxies at z<2. In this paper, we perform high-resolution re-simulations of a massive (M_star~10^11 M_sol) galaxy at z~2.3, drawn from the Feedback in Realistic Environments (FIRE) project. In the simulation without AGN feedback, the galaxy experiences a rapid starburst and shrinking of its half… ▽ More Without AGN feedback, simulated massive, star-forming galaxies become too compact relative to observed galaxies at z<2. In this paper, we perform high-resolution re-simulations of a massive (M_star~10^11 M_sol) galaxy at z~2.3, drawn from the Feedback in Realistic Environments (FIRE) project. In the simulation without AGN feedback, the galaxy experiences a rapid starburst and shrinking of its half-mass radius. We experiment with driving mechanical AGN winds, using a state-of-the-art hyper-Lagrangian refinement technique to increase particle resolution. These winds reduce the gas surface density in the inner regions of the galaxy, suppressing the compact starburst and maintaining an approximately constant half-mass radius. Using radiative transfer, we study the impact of AGN feedback on the magnitude and extent of the multi-wavelength continuum emission. When AGN winds are included, the suppression of the compact, dusty starburst results in lowered flux at FIR wavelengths (due to decreased star formation) but increased flux at optical-to-near-IR wavelengths (due to decreased dust attenuation, in spite of the lowered star formation rate), relative to the case without AGN winds. The FIR half-light radius decreases from ~1 kpc to ~0.1 kpc in <40 Myr when AGN winds are not included, but increases to ~2 kpc when they are. Interestingly, the half-light radius at optical-NIR wavelengths remains approximately constant over 35 Myr, for simulations with and without AGN winds. In the case without winds, this occurs despite the rapid compaction, and is due to heavy dust obscuration in the inner regions of the galaxy. This work highlights the importance of forward-modelling when comparing simulated and observed galaxy populations. △ Less

Submitted 17 May, 2023; v1 submitted 22 March, 2023; originally announced March 2023.

Comments: 14 pages, 10 figures. Accepted for publication in MNRAS

Arkenstone I: A Novel Method for Robustly Capturing High Specific Energy Outflows In Cosmological Simulations

Authors: Matthew C. Smith, Drummond B. Fielding, Greg L. Bryan, Chang-Goo Kim, Eve C. Ostriker, Rachel S. Somerville, Jonathan Stern, Kung-Yi Su, Rainer Weinberger, Chia-Yu Hu, John C. Forbes, Lars Hernquist, Blakesley Burkhart, Yuan Li

Abstract: Arkenstone is a new model for multiphase, stellar feedback driven galactic winds designed for inclusion in coarse resolution cosmological simulations. In this first paper of a series, we describe the features that allow Arkenstone to properly treat high specific energy wind components and demonstrate them using idealised non-cosmological simulations of a galaxy with a realistic CGM, using the Arep… ▽ More Arkenstone is a new model for multiphase, stellar feedback driven galactic winds designed for inclusion in coarse resolution cosmological simulations. In this first paper of a series, we describe the features that allow Arkenstone to properly treat high specific energy wind components and demonstrate them using idealised non-cosmological simulations of a galaxy with a realistic CGM, using the Arepo code. Hot, fast gas phases with low mass loadings are predicted to dominate the energy content of multiphase outflows. In order to treat the huge dynamic range of spatial scales involved in cosmological galaxy formation at feasible computational expense, cosmological volume simulations typically employ a Lagrangian code or else use adaptive mesh refinement with a quasi-Lagrangian refinement strategy. However, it is difficult to inject a high specific energy wind in a Lagrangian scheme without incurring artificial burstiness. Additionally, the low densities inherent to this type of flow result in poor spatial resolution. Arkenstone addresses these issues with a novel scheme for coupling energy into the ISM/CGM transition region which also provides the necessary level of refinement at the base of the wind. In the absence of our improvements, we show that poor spatial resolution near the sonic point of a hot, fast outflow leads to an underestimation of gas acceleration as the wind propagates. We explore the different mechanisms by which low and high specific energy winds can regulate the SFR of galaxies. In future work, we will demonstrate other aspects of the Arkenstone model. △ Less

Submitted 7 March, 2024; v1 submitted 17 January, 2023; originally announced January 2023.

Comments: Published MNRAS, 27 pages, 17 figures; updated for consistency with journal version

Journal ref: MNRAS, 527, 1216 (2024)

Local positive feedback in the overall negative: the impact of quasar winds on star formation in the FIRE cosmological simulations

Authors: Jonathan Mercedes-Feliz, Daniel Anglés-Alcázar, Christopher C. Hayward, Rachel K. Cochrane, Bryan A. Terrazas, Sarah Wellons, Alexander J. Richings, Claude-André Faucher-Giguère, Jorge Moreno, Kung Yi Su, Philip F. Hopkins, Eliot Quataert, Dušan Kereš

Abstract: Negative feedback from accreting supermassive black holes is regarded as a key ingredient in suppressing star formation and quenching massive galaxies. However, several models and observations suggest that black hole feedback may have a positive effect, triggering star formation by compressing interstellar medium gas to higher densities. We investigate the dual role of black hole feedback using co… ▽ More Negative feedback from accreting supermassive black holes is regarded as a key ingredient in suppressing star formation and quenching massive galaxies. However, several models and observations suggest that black hole feedback may have a positive effect, triggering star formation by compressing interstellar medium gas to higher densities. We investigate the dual role of black hole feedback using cosmological hydrodynamic simulations from the Feedback In Realistic Environments (FIRE) project, including a novel implementation of hyper-refined accretion-disc winds. Focusing on a massive, star-forming galaxy at $z \sim 2$ ($M_{\rm halo} \sim 10^{12.5} \, {\rm M}_{\odot}$), we show that strong quasar winds with kinetic power $\sim$10$^{46}$ erg/s acting for $>$20$\,$Myr drive the formation of a central gas cavity and can dramatically reduce the star formation rate surface density across the galaxy disc. The suppression of star formation is primarily driven by reducing the amount of gas that can become star-forming, compared to directly evacuating the pre-existing star-forming gas reservoir (preventive feedback dominates over ejective feedback). Despite the global negative impact of quasar winds, we identify several plausible signatures of local positive feedback, including: (1) spatial anti-correlation of wind-dominated regions and star-forming clumps, (2) higher local star formation efficiency in compressed gas near the edge of the cavity, and (3) increased local contribution of outflowing material to star formation. Stars forming under the presence of quasar winds tend to do so at larger radial distances. Our results suggest that positive and negative AGN feedback can coexist in galaxies, but local positive triggering of star formation plays a minor role in global galaxy growth. △ Less

Submitted 1 August, 2023; v1 submitted 4 January, 2023; originally announced January 2023.

Comments: 17 pages, 12 figures. Accepted for publication in MNRAS

Active galactic nucleus jet feedback in hydrostatic halos

Authors: Rainer Weinberger, Kung-Yi Su, Kristian Ehlert, Christoph Pfrommer, Lars Hernquist, Greg L. Bryan, Volker Springel, Yuan Li, Blakesley Burkhart, Ena Choi, Claude-André Faucher-Giguère

Abstract: Feedback driven by jets from active galactic nuclei is believed to be responsible for reducing cooling flows in cool-core galaxy clusters. We use simulations to model feedback from hydrodynamic jets in isolated halos. While the jet propagation converges only after the diameter of the jet is well resolved, reliable predictions about the effects these jets have on the cooling time distribution funct… ▽ More Feedback driven by jets from active galactic nuclei is believed to be responsible for reducing cooling flows in cool-core galaxy clusters. We use simulations to model feedback from hydrodynamic jets in isolated halos. While the jet propagation converges only after the diameter of the jet is well resolved, reliable predictions about the effects these jets have on the cooling time distribution function only require resolutions sufficient to keep the jet-inflated cavities stable. Comparing different model variations, as well as an independent jet model using a different hydrodynamics code, we show that the dominant uncertainties are the choices of jet properties within a given model. Independent of implementation, we find that light, thermal jets with low momentum flux tend to delay the onset of a cooling flow more efficiently on a $50$ Myr timescale than heavy, kinetic jets. The delay of the cooling flow originates from a displacement and boost in entropy of the central gas. If the jet kinetic luminosity depends on accretion rate, collimated, light, hydrodynamic jets are able to reduce cooling flows in halos, without a need for jet precession or wide opening angles. Comparing the jet feedback with a `kinetic wind' implementation shows that equal amounts of star formation rate reduction can be achieved by different interactions with the halo gas: the jet has a larger effect on the hot halo gas while leaving the denser, star forming phase in place, while the wind acts more locally on the star forming phase, which manifests itself in different time-variability properties. △ Less

Submitted 5 May, 2023; v1 submitted 21 November, 2022; originally announced November 2022.

Comments: 22 pages, 20 figures, accepted for publication in MNRAS, comments welcome

The Geometry of the G29-38 White Dwarf Dust Disk from Radiative Transfer Modeling

Authors: Nicholas P. Ballering, Colette I. Levens, Kate Y. L. Su, L. Ilsedore Cleeves

Abstract: Many white dwarfs host disks of dust produced by disintegrating planetesimals and revealed by infrared excesses. The disk around G29-38 was the first to be discovered and is now well-observed, yet we lack a cohesive picture of its geometry and dust properties. Here we model the G29-38 disk for the first time using radiative transfer calculations that account for radial and vertical temperature and… ▽ More Many white dwarfs host disks of dust produced by disintegrating planetesimals and revealed by infrared excesses. The disk around G29-38 was the first to be discovered and is now well-observed, yet we lack a cohesive picture of its geometry and dust properties. Here we model the G29-38 disk for the first time using radiative transfer calculations that account for radial and vertical temperature and optical depth gradients. We arrive at a set of models that can match the available infrared measurements well, although they overpredict the width of the 10 $μm$ silicate feature. The resulting set of models has a disk inner edge located at 92-100 $R_\text{WD}$ (where $R_\text{WD}$ is the white dwarf radius). This is farther from the star than inferred by previous modeling efforts due to the presence of a directly illuminated front edge to the disk. The radial width of the disk is narrow ($\leq$10 $R_\text{WD}$); such a feature could be explained by inefficient spreading or the proximity of the tidal disruption radius to the sublimation radius. The models have a half-opening angle of $\geq$1.4$^\circ$. Such structure would be in strong contradiction with the commonly employed flat-disk model analogous to the rings of Saturn, and in line with the vertical structure of main-sequence debris disks. Our results are consistent with the idea that disks are collisionally active and continuously fed with new material, rather than evolving passively after the disintegration of a single planetesimal. △ Less

Submitted 31 October, 2022; originally announced November 2022.

Comments: 12 pages, 5 figures, accepted for publication in ApJ

Mid-infrared time-domain study of recent dust production events in the extreme debris disc of TYC 4209-1322-1

Authors: A. Moór, P. Ábrahám, Á. Kóspál, K. Y. L. Su, G. H. Rieke, G. Cataldi, A. Bódi, Zs. Bognár, B. Cseh, G. Csörnyei, N. Egei, A. Farkas, O. Hanyecz, B. Ignácz, Cs. Kalup, R. Könyves-Tóth, L. Kriskovics, L. Mészáros, A. Pál, A. Ordasi, K. Sárneczky, B. Seli, Á. Sódor, R. Szakáts, J. Vinkó , et al. (1 additional authors not shown)

Abstract: Extreme debris discs are characterized by unusually strong mid-infrared excess emission, which often proves to be variable. The warm dust in these discs is of transient nature and is likely related to a recent giant collision occurring close to the star in the terrestrial region. Here we present the results of a 877 days long, gap-free photometric monitoring performed by the Spitzer Space Telescop… ▽ More Extreme debris discs are characterized by unusually strong mid-infrared excess emission, which often proves to be variable. The warm dust in these discs is of transient nature and is likely related to a recent giant collision occurring close to the star in the terrestrial region. Here we present the results of a 877 days long, gap-free photometric monitoring performed by the Spitzer Space Telescope of the recently discovered extreme debris disc around TYC 4209-1322-1. By combining these observations with other time-domain optical and mid-infrared data, we explore the disc variability of the last four decades with particular emphasis on the last 12 yr. During the latter interval the disc showed substantial changes, the most significant was the brightening and subsequent fading between 2014 and 2018 as outlined in WISE data. The Spitzer light curves outline the fading phase and a subsequent new brightening of the disc after 2018, revealing an additional flux modulation with a period of ~39 days on top of the long-term trend. We found that all these variations can be interpreted as the outcome of a giant collision that happened at an orbital radius of ~0.3 au sometime in 2014. Our analysis implies that a collision on a similar scale could have taken place around 2010, too. The fact that the disc was already peculiarly dust rich 40 yr ago, as implied by IRAS data, suggests that these dust production events belong to a chain of large impacts triggered by an earlier even more catastrophic collision. △ Less

Submitted 21 October, 2022; originally announced October 2022.

Comments: 18 pages, 12 figures, published in MNRAS

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 516, Issue 4, November 2022, Pages 5684-5701

Self-regulation of black hole accretion via jets in early protogalaxies

Authors: Kung-Yi Su, Greg L. Bryan, Zoltán Haiman, Rachel S. Somerville, Christopher C. Hayward, Claude-André Faucher-Giguère

Abstract: The early growth of black holes (BHs) in high-redshift galaxies is likely regulated by their feedback on the surrounding gas. While radiative feedback has been extensively studied, the role of mechanical feedback has received comparatively less scrutiny to date. Here we use high-resolution parsec-scale hydrodynamical simulations to study jet propagation and its effect on BH accretion onto 100… ▽ More The early growth of black holes (BHs) in high-redshift galaxies is likely regulated by their feedback on the surrounding gas. While radiative feedback has been extensively studied, the role of mechanical feedback has received comparatively less scrutiny to date. Here we use high-resolution parsec-scale hydrodynamical simulations to study jet propagation and its effect on BH accretion onto 100 ${\rm M_\odot}$ BHs in the dense, low-metallicity gas expected in early protogalaxies. As the jet propagates, it shocks the surrounding gas and forms a jet cocoon. The cocoon consists of a rapidly-cooling cold phase at the interface with the background gas and an over-pressured subsonic phase of reverse shock-heated gas filling the cocoon interior. We systematically vary the background gas density and temperature, BH feedback efficiency, and the jet model. We found that the jet cocoon width roughly follows a scaling derived by assuming momentum conservation in the jet propagation direction, and energy conservation in the lateral directions. Depending on the assumed gas and jet properties, the cocoon either stays elongated out to a large radius or isotropizes before reaching the Bondi radius, forming a nearly spherical bubble. Lower jet velocities and higher background gas densities result in self-regulation to higher momentum fluxes and elongated cocoons. In all cases, the outward momentum flux of the cocoon balances the inward momentum flux of the inflowing gas near the Bondi radius, which ultimately regulates BH accretion. The larger the distance the jet cocoon reaches, the longer the variability timescale of the BH accretion rate. Overall, the average accretion rate always remains below the Bondi rate, and exceeds the Eddington rate only if the ambient medium is dense and cold, and/or the jet is weak. We derive the combination of jet and ambient gas parameters yielding super-Eddington growth. △ Less

Submitted 22 July, 2022; originally announced July 2022.

Comments: 16 pages, 11 figures

The James Webb Space Telescope Absolute Flux Calibration. I. Program Design and Calibrator Stars

Authors: Karl D. Gordon, Ralph Bohlin, G. C. Sloan, George Rieke, Kevin Volk, Martha Boyer, James Muzerolle, Everett Schlawin, Susana E. Deustua, Dean C. Hines, Kathleen E. Kraemer, Susan E. Mullally, Kate Y. L. Su

Abstract: It is critical for James Webb Space Telescope (JWST) science that instrumental units are converted to physical units. We detail the design of the JWST absolute flux calibration program that has the core goal of ensuring a robust flux calibration internal to and between all the science instruments for both point and extended source science. This program will observe a sample of calibration stars th… ▽ More It is critical for James Webb Space Telescope (JWST) science that instrumental units are converted to physical units. We detail the design of the JWST absolute flux calibration program that has the core goal of ensuring a robust flux calibration internal to and between all the science instruments for both point and extended source science. This program will observe a sample of calibration stars that have been extensively vetted based mainly on Hubble Space Telescope, Spitzer Space Telescope, and Transiting Exoplanet Survey Satellite observations. The program uses multiple stars of three different, well understood types (hot stars, A dwarfs, and solar analogs) to allow for the statistical (within a type) and systematic (between types) uncertainties to be quantified. The program explicitly includes observations to calibrate every instrument mode, further vet the set of calibration stars, measure the instrumental repeatability, measure the relative calibration between subarrays and full frame, and check the relative calibration between faint and bright stars. For photometry, we have set up our calibration to directly support both the convention based on the band average flux density and the convention based on the flux density at a fixed wavelength. △ Less

Submitted 13 April, 2022; originally announced April 2022.

Comments: 13 pages, 3 figures, AJ, in press

Exploring supermassive black hole physics and galaxy quenching across halo mass in FIRE cosmological zoom simulations

Authors: Sarah Wellons, Claude-André Faucher-Giguère, Philip F. Hopkins, Eliot Quataert, Daniel Anglés-Alcázar, Robert Feldmann, Christopher C. Hayward, Dušan Kereš, Kung-Yi Su, Andrew Wetzel

Abstract: Feedback from accreting supermassive black holes (SMBHs) is thought to be a primary driver of quenching in massive galaxies, but the best way to implement SMBH physics into galaxy formation simulations remains ambiguous. As part of the Feedback in Realistic Environments (FIRE) project, we explore the effects of different modeling choices for SMBH accretion and feedback in a suite of $\sim500$ cosm… ▽ More Feedback from accreting supermassive black holes (SMBHs) is thought to be a primary driver of quenching in massive galaxies, but the best way to implement SMBH physics into galaxy formation simulations remains ambiguous. As part of the Feedback in Realistic Environments (FIRE) project, we explore the effects of different modeling choices for SMBH accretion and feedback in a suite of $\sim500$ cosmological zoom-in simulations across a wide range of halo mass (10^10-10^13 Msun). Within the suite, we vary the numerical schemes for BH accretion and feedback, the accretion efficiency, and the strength of mechanical, radiative, and cosmic ray feedback independently. We then compare the outcomes to observed galaxy scaling relations. We find several models that satisfy the observational constraints, and for which the energetics in different feedback channels are physically plausible. Interestingly, cosmic rays accelerated by SMBHs play an important role in many successful models. However, it is non-trivial to reproduce scaling relations across halo mass, and many model variations produce qualitatively incorrect results regardless of parameter choices. The growth of stellar and BH mass are closely related: for example, over-massive BHs tend to over-quench galaxies. BH mass is most strongly affected by the choice of accretion efficiency in high-mass halos, but by feedback efficiency in low-mass halos. The amount of star formation suppression by SMBH feedback in low-mass halos is determined primarily by the time-integrated feedback energy. For massive galaxies, the "responsiveness" of a model (i.e. how quickly and powerfully the BH responds to gas available for accretion) is an additional important factor for quenching. △ Less

Submitted 11 March, 2022; originally announced March 2022.

Comments: 18 pages, 9 figures. Submitted to MNRAS, comments welcome

A Star-sized Impact-Produced Dust Clump in the Terrestrial Zone of the HD 166191 System

Authors: Kate Y. L. Su, Grant M. Kennedy, Everett Schlawin, Alan P. Jackson, George H. Rieke

Abstract: We report on five years of 3-5 micron photometry measurements obtained by warm Spitzer to track the dust debris emission in the terrestrial zone of HD 166191 in combination with simultaneous optical data. We show that the debris production in this young (~10 Myr) system increased significantly in early 2018 and reached a record high level (almost double by mid 2019) by the end of the Spitzer missi… ▽ More We report on five years of 3-5 micron photometry measurements obtained by warm Spitzer to track the dust debris emission in the terrestrial zone of HD 166191 in combination with simultaneous optical data. We show that the debris production in this young (~10 Myr) system increased significantly in early 2018 and reached a record high level (almost double by mid 2019) by the end of the Spitzer mission (early 2020), suggesting intense collisional activity in its terrestrial zone likely due to either initial assembling of terrestrial planets through giant impacts or dynamical shake-up from unseen planet-mass objects or recent planet migration. This intense activity is further highlighted by detecting a star-size dust clump, passing in front of the star, in the midst of its infrared brightening. We constrain the minimum size and mass of the clump using multiwavelength transit profiles and conclude that the dust clump is most likely created by a large impact involving objects of several hundred km in size with an apparent period of 142 days (i.e., 0.62 au assuming a circular orbit). The system's evolutionary state (right after the dispersal of its gas-rich disk) makes it extremely valuable to learn about the process of terrestrial planet formation and planetary architecture through future observations. △ Less

Submitted 4 March, 2022; originally announced March 2022.

Comments: 17 pages, 5 figures, accepted by ApJ

FIRE-3: Updated Stellar Evolution Models, Yields, & Microphysics and Fitting Functions for Applications in Galaxy Simulations

Authors: Philip F. Hopkins, Andrew Wetzel, Coral Wheeler, Robyn Sanderson, Michael Y. Grudic, Omid Sameie, Michael Boylan-Kolchin, Matthew Orr, Xiangcheng Ma, Claude-Andre Faucher-Giguere, Dusan Keres, Eliot Quataert, Kung-Yi Su, Jorge Moreno, Robert Feldmann, James S. Bullock, Sarah R. Loebman, Daniel Angles-Alcazar, Jonathan Stern, Lina Necib, Christopher C. Hayward

Abstract: Increasingly, uncertainties in predictions from galaxy formation simulations (at sub-Milky Way masses) are dominated by uncertainties in stellar evolution inputs. In this paper, we present the full set of updates from the FIRE-2 version of the Feedback In Realistic Environments (FIRE) project code, to the next version, FIRE-3. While the transition from FIRE-1 to FIRE-2 focused on improving numeric… ▽ More Increasingly, uncertainties in predictions from galaxy formation simulations (at sub-Milky Way masses) are dominated by uncertainties in stellar evolution inputs. In this paper, we present the full set of updates from the FIRE-2 version of the Feedback In Realistic Environments (FIRE) project code, to the next version, FIRE-3. While the transition from FIRE-1 to FIRE-2 focused on improving numerical methods, here we update the stellar evolution tracks used to determine stellar feedback inputs, e.g. stellar mass-loss (O/B and AGB), spectra (luminosities and ionization rates), and supernova rates (core-collapse and Ia), as well as detailed mass-dependent yields. We also update the low-temperature cooling and chemistry, to enable improved accuracy at $T \lesssim 10^{4}\,$K and densities $n\gg 1\,{\rm cm^{-3}}$, and the meta-galactic ionizing background. All of these synthesize newer empirical constraints on these quantities and updated stellar evolution and yield models from a number of groups, addressing different aspects of stellar evolution. To make the updated models as accessible as possible, we provide fitting functions for all of the relevant updated tracks, yields, etc, in a form specifically designed so they can be directly 'plugged in' to existing galaxy formation simulations. We also summarize the default FIRE-3 implementations of 'optional' physics, including spectrally-resolved cosmic rays and supermassive black hole growth and feedback. △ Less

Submitted 28 February, 2022; originally announced March 2022.

Comments: 26 pages, 12 figures. Submitted to MNRAS. Comments welcome

Infrared Absolute Calibration I: Comparison of Sirius with Fainter Calibration Stars

Authors: G. H. Rieke, K. Y. L. Su, G. C. Sloan, E. Schlawin

Abstract: A challenge in absolute calibration is to relate very bright stars with physical flux measurements to faint ones within range of modern instruments, e.g. those on large groundbased telescopes or on the James Webb Space Telescope (JWST). We propose Sirius as the fiducial color standard: it is an A0V star that is slowly rotating and does not have infrared excesses either due to hot dust or a planeta… ▽ More A challenge in absolute calibration is to relate very bright stars with physical flux measurements to faint ones within range of modern instruments, e.g. those on large groundbased telescopes or on the James Webb Space Telescope (JWST). We propose Sirius as the fiducial color standard: it is an A0V star that is slowly rotating and does not have infrared excesses either due to hot dust or a planetary debris disk; it also has a number of accurate (~ 1 - 2%) absolute flux measurements. We transfer the near infrared flux from Sirius accurately to BD +60 1753, an unobscured early A-type star (A1V, V ~ 9.6, E(B-V) ~ 0.009) that is faint enough to serve as a primary absolute flux calibrator for JWST. Its near-infrared spectral energy distribution and that of Sirius should be virtually identical. We have determined its output relative to that of Sirius in a number of different ways, all of which give consistent results within ~ 1%. We also transfer the calibration to GSPC P330-E, a well calibrated close solar analog (G2V). We have emphasized the 2MASS K_S band since it represents a large number and long history of measurements, but theoretical spectra (i.e., from CALSPEC) of these stars can be used to extend this result throughout the near- and mid-infrared. △ Less

Submitted 19 November, 2021; originally announced November 2021.

Comments: Astronomical Journal, accepted

Accurate Photometry of Saturated Stars Using the Point-Spread-Function Wing Technique with Spitzer

Authors: Kate Y. L. Su, G. H. Rieke, M. Marengo, Everett Schlawin

Abstract: We report Spitzer 3.6 and 4.5 $μ$m photometry of 11 bright stars relative to Sirius, exploiting the unique optical stability of the Spitzer Space Telescope point spread function (PSF). Spitzer's extremely stable beryllium optics in its isothermal environment enables precise comparisons in the wings of the PSF from heavily saturated stars. These bright stars stand as the primary sample to improve s… ▽ More We report Spitzer 3.6 and 4.5 $μ$m photometry of 11 bright stars relative to Sirius, exploiting the unique optical stability of the Spitzer Space Telescope point spread function (PSF). Spitzer's extremely stable beryllium optics in its isothermal environment enables precise comparisons in the wings of the PSF from heavily saturated stars. These bright stars stand as the primary sample to improve stellar models, and to transfer the absolute flux calibration of bright standard stars to a sample of fainter standards useful for missions like JWST and for large groundbased telescopes. We demonstrate that better than 1% relative photometry can be achieved using the PSF wing technique in the radial range of 20--100\arcsec\ for stars that are fainter than Sirius by 8 mag (from outside the saturated core to a large radius where a high signal-to-noise profile can still be obtained). We test our results by (1) comparing the [3.6]$-$[4.5] color with that expected between the WISE W1 and W2 bands, (2) comparing with stars where there is accurate $K_{\text{S}}$ photometry, and (3) also comparing with relative fluxes obtained with the DIRBE instrument on COBE. These tests confirm that relative photometry is achieved to better than 1%. △ Less

Submitted 19 November, 2021; originally announced November 2021.

Comments: 13 pages, 5 figures, 4 tables, accepted for publication in the AAS Journal

Carbon monoxide gas produced by a giant impact in the inner region of a young system

Authors: Tajana Schneiderman, Luca Matrà, Alan P. Jackson, Grant M. Kennedy, Quentin Kral, Sebastián Marino, Karin I. Öberg, Kate Y. L. Su, David J. Wilner, Mark C. Wyatt

Abstract: Models of terrestrial planet formation predict that the final stages of planetary assembly, lasting tens of millions of years beyond the dispersal of young protoplanetary disks, are dominated by planetary collisions. It is through these giant impacts that planets like the young Earth grow to their final mass and achieve long-term stable orbital configurations. A key prediction is that these impact… ▽ More Models of terrestrial planet formation predict that the final stages of planetary assembly, lasting tens of millions of years beyond the dispersal of young protoplanetary disks, are dominated by planetary collisions. It is through these giant impacts that planets like the young Earth grow to their final mass and achieve long-term stable orbital configurations. A key prediction is that these impacts produce debris. To date, the most compelling observational evidence for post-impact debris comes from the planetary system around the nearby 23 Myr-old A star HD 172555. This system shows large amounts of fine dust with an unusually steep size distribution and atypical dust composition, previously attributed to either a hypervelocity impact or a massive asteroid belt. Here, we report the spectrally resolved detection of a CO gas ring co-orbiting with dusty debris between ~6-9 au - a region analogous to the outer terrestrial planet region of our Solar System. Taken together, the dust and CO detections favor a giant impact between large, volatile-rich bodies. This suggests that planetary-scale collisions, analogous to the Moon-forming impact, can release large amounts of gas as well as debris, and that this gas is observable, providing a window into the composition of young planets. △ Less

Submitted 28 October, 2021; originally announced October 2021.

Comments: 21 pages (double spaced, including references and figures), 3 figures, 1 table, published in Nature 20/10/2021, published version available online through SharedIt initiative at https://rdcu.be/cAlEy

Journal ref: Nature 598, 425-428 (2021)

Design of the vacuum high contrast imaging testbed for CDEEP, the Coronagraphic Debris and Exoplanet Exploring Pioneer

Authors: Erin R. Maier, Ewan S. Douglas, Daewook Kim, Kate Su, Jaren N. Ashcraft, James B. Breckinridge, Supriya Chakrabarti, Heejoo Choi, Elodie Choquet, Thomas E. Connors, Olivier Durney, John Debes, Kerry L. Gonzales, Charlotte E. Guthery, Christian A. Haughwout, James C. Heath, Justin Hyatt, Jennifer Lumbres, Jared R. Males, Elisabeth C. Matthews, Kian Milani, Oscar M. Montoya, Mamadou N'Diaye, Jamison Noenickx, Leonid Pogorelyuk , et al. (4 additional authors not shown)

Abstract: The Coronagraphic Debris Exoplanet Exploring Payload (CDEEP) is a Small-Sat mission concept for high contrast imaging of circumstellar disks. CDEEP is designed to observe disks in scattered light at visible wavelengths at a raw contrast level of 10^-7 per resolution element (10^-8 with post processing). This exceptional sensitivity will allow the imaging of transport dominated debris disks, quanti… ▽ More The Coronagraphic Debris Exoplanet Exploring Payload (CDEEP) is a Small-Sat mission concept for high contrast imaging of circumstellar disks. CDEEP is designed to observe disks in scattered light at visible wavelengths at a raw contrast level of 10^-7 per resolution element (10^-8 with post processing). This exceptional sensitivity will allow the imaging of transport dominated debris disks, quantifying the albedo, composition, and morphology of these low-surface brightness disks. CDEEP combines an off-axis telescope, microelectromechanical systems (MEMS) deformable mirror, and a vector vortex coronagraph (VVC). This system will require rigorous testing and characterization in a space environment. We report on the CDEEP mission concept, and the status of the vacuum-compatible CDEEP prototype testbed currently under development at the University of Arizona, including design development and the results of simulations to estimate performance. △ Less

Submitted 26 September, 2021; originally announced September 2021.

Comments: 17 pages, 12 figures, Published in proceedings of SPIE Astronomical Telescopes + Instrumentation 2020

Extreme Variability of the V488 Persei Debris Disk

Authors: G. H. Rieke, K. Y. L. Su, Carl Melis, Andras Gaspar

Abstract: V488 Persei is the most extreme debris disk known in terms of the fraction of the stellar luminosity it intercepts and reradiates. The infrared output of its disk is extremely variable, similar in this respect to the most variable disk known previously, that around ID8 in NGC 2547. We show that the variations are likely to be due to collisions of large planetesimals (> 100 km in diameter) in a bel… ▽ More V488 Persei is the most extreme debris disk known in terms of the fraction of the stellar luminosity it intercepts and reradiates. The infrared output of its disk is extremely variable, similar in this respect to the most variable disk known previously, that around ID8 in NGC 2547. We show that the variations are likely to be due to collisions of large planetesimals (> 100 km in diameter) in a belt being stirred gravitationally by a planetary or low-mass-brown-dwarf member of a planetary system around the star. The dust being produced by the resulting collisions is falling into the star due to drag by the stellar wind. The indicated planetesimal destruction rate is so high that it is unlikely that the current level of activity can persist for much longer than ~ 1000 - 10,000 years, and it may signal a major realignment of the configuration of the planetary system. △ Less

Submitted 5 August, 2021; originally announced August 2021.

Comments: ApJ, in press

LBT Reveals Large Dust Particles and a High Mass Loss Rate for K2-22 b

Authors: Everett Schlawin, Kate Y. L. Su, Terry Herter, Andrew Ridden-Harper, Daniel Apai

Abstract: The disintegrating planet candidate K2-22 b shows periodic and stochastic transits best explained by an escaping debris cloud. However, the mechanism that creates the debris cloud is unknown. The grain size of the debris as well as its sublimation rate can be helpful in understanding the environment that disintegrates the planet. Here, we present simultaneous photometry with the g band at 0.48 mic… ▽ More The disintegrating planet candidate K2-22 b shows periodic and stochastic transits best explained by an escaping debris cloud. However, the mechanism that creates the debris cloud is unknown. The grain size of the debris as well as its sublimation rate can be helpful in understanding the environment that disintegrates the planet. Here, we present simultaneous photometry with the g band at 0.48 microns and KS band at 2.1 microns using the Large Binocular Telescope. During an event with very low dust activity, we put a new upper limit on the size of the planet of 0.71 earth radii or 4500 km. We also detected a medium-depth transit which can be used to constrain the dust particle sizes. We find that the median particle size must be larger than about 0.5 to 1.0 microns, depending on the composition of the debris. This leads to a high mass loss rate of about 3e8 kg/s that is consistent with hydrodynamic escape models. If they are produced by some alternate mechanism such as explosive volcanism, it would require extraordinary geological activity. Combining our upper limits on the planet size with the high mass loss rate, we find a lifetime of the planet of less than 370 Myr. This drops to just 21 Myr when adopting the 0.02 earth masses predicted from hydrodynamical models. △ Less

Submitted 14 June, 2021; originally announced June 2021.

Comments: AJ, accepted, 13 pages

A near-infrared interferometric survey of debris-disk stars. VII. The hot/warm dust connection

Authors: O. Absil, L. Marion, S. Ertel, D. Defrère, G. M. Kennedy, A. Romagnolo, J. -B. Le Bouquin, V. Christiaens, J. Milli, A. Bonsor, J. Olofsson, K. Y. L. Su, J. -C. Augereau

Abstract: (abridged) Context. The origin of hot exozodiacal dust and its connection with outer dust reservoirs remains unclear. Aims. We aim to explore the possible connection between hot exozodiacal dust and warm dust reservoirs (> 100 K) in asteroid belts. Methods. We use precision near-infrared interferometry with VLTI/PIONIER to search for resolved emission at H band around a selected sample of nearby s… ▽ More (abridged) Context. The origin of hot exozodiacal dust and its connection with outer dust reservoirs remains unclear. Aims. We aim to explore the possible connection between hot exozodiacal dust and warm dust reservoirs (> 100 K) in asteroid belts. Methods. We use precision near-infrared interferometry with VLTI/PIONIER to search for resolved emission at H band around a selected sample of nearby stars. Results. Our observations reveal the presence of resolved near-infrared emission around 17 out of 52 stars, four of which are shown to be due to a previously unknown stellar companion. The 13 other H-band excesses are thought to originate from the thermal emission of hot dust grains. Taking into account earlier PIONIER observations, and after reevaluating the warm dust content of all our PIONIER targets through spectral energy distribution modeling, we find a detection rate of 17.1(+8.1)(-4.6)% for H-band excess around main sequence stars hosting warm dust belts, which is statistically compatible with the occurrence rate of 14.6(+4.3)(-2.8)% found around stars showing no signs of warm dust. After correcting for the sensitivity loss due to partly unresolved hot disks, under the assumption that they are arranged in a thin ring around their sublimation radius, we however find tentative evidence at the 3σ level that H-band excesses around stars with outer dust reservoirs (warm or cold) could be statistically larger than H-band excesses around stars with no detectable outer dust. Conclusions. Our observations do not suggest a direct connection between warm and hot dust populations, at the sensitivity level of the considered instruments, although they bring to light a possible correlation between the level of H-band excesses and the presence of outer dust reservoirs in general. △ Less

Submitted 29 April, 2021; originally announced April 2021.

Comments: Accepted for publication in A&A

Journal ref: A&A 651, A45 (2021)

The HOSTS survey: evidence for an extended dust disk and constraints on the presence of giant planets in the Habitable Zone of $β$ Leo

Authors: D. Defrère, P. M. Hinz, G. M. Kennedy, J. Stone, J. Rigley, S. Ertel, A. Gaspar, V. P. Bailey, W. F. Hoffmann, B. Mennesson, R. Millan-Gabet, W. C. Danchi, O. Absil, P. Arbo, C. Beichman, M. Bonavita, G. Brusa, G. Bryden, E. C. Downey, S. Esposito, P. Grenz, C. Haniff, J. M. Hill, J. M. Leisenring, J. R. Males , et al. (16 additional authors not shown)

Abstract: The young (50-400 Myr) A3V star $β$ Leo is a primary target to study the formation history and evolution of extrasolar planetary systems as one of the few stars with known hot ($\sim$1600$^\circ$K), warm ($\sim$600$^\circ$K), and cold ($\sim$120$^\circ$K) dust belt components. In this paper, we present deep mid-infrared measurements of the warm dust brightness obtained with the Large Binocular Tel… ▽ More The young (50-400 Myr) A3V star $β$ Leo is a primary target to study the formation history and evolution of extrasolar planetary systems as one of the few stars with known hot ($\sim$1600$^\circ$K), warm ($\sim$600$^\circ$K), and cold ($\sim$120$^\circ$K) dust belt components. In this paper, we present deep mid-infrared measurements of the warm dust brightness obtained with the Large Binocular Telescope Interferometer (LBTI) as part of its exozodiacal dust survey (HOSTS). The measured excess is 0.47\%$\pm$0.050\% within the central 1.5 au, rising to 0.81\%$\pm$0.026\% within 4.5 au, outside the habitable zone of $β$~Leo. This dust level is 50 $\pm$ 10 times greater than in the solar system's zodiacal cloud. Poynting-Robertson drag on the cold dust detected by Spitzer and Herschel under-predicts the dust present in the habitable zone of $β$~Leo, suggesting an additional delivery mechanism (e.g.,~comets) or an additional belt at $\sim$5.5 au. A model of these dust components is provided which implies the absence of planets more than a few Saturn masses between $\sim$5 au and the outer belt at $\sim$40 au. We also observationally constrain giant planets with the LBTI imaging channel at 3.8~$μ$m wavelength. Assuming an age of 50 Myr, any planet in the system between approximately 5 au to 50 au must be less than a few Jupiter masses, consistent with our dust model. Taken together, these observations showcase the deep contrasts and detection capabilities attainable by the LBTI for both warm exozodiacal dust and giant exoplanets in or near the habitable zone of nearby stars. △ Less

Submitted 4 March, 2021; originally announced March 2021.

Comments: 11 pages, 9 figures, accepted for publication in Astronomical Journal

Which AGN Jets Quench Star Formation in Massive Galaxies?

Authors: Kung-Yi Su, Philip F. Hopkins, Greg L. Bryan, Rachel S. Somerville, Christopher C. Hayward, Daniel Anglés-Alcázar, Claude-André Faucher-Giguère, Sarah Wellons, Jonathan Stern, Bryan A. Terrazas, T. K. Chan, Matthew E. Orr, Cameron Hummels, Robert Feldmann, Dušan Kereš

Abstract: Without additional heating, radiative cooling of gas in the halos of massive galaxies (Milky Way and above) produces cold gas or stars in excess of that observed. Previous work suggested that AGN jets are likely required, but the form of jet energy required to quench remains unclear. This is particularly challenging for galaxy simulations, in which the resolution is orders of magnitude coarser tha… ▽ More Without additional heating, radiative cooling of gas in the halos of massive galaxies (Milky Way and above) produces cold gas or stars in excess of that observed. Previous work suggested that AGN jets are likely required, but the form of jet energy required to quench remains unclear. This is particularly challenging for galaxy simulations, in which the resolution is orders of magnitude coarser than necessary to form and evolve the jet. On such scales, the uncertain parameters include: jet energy form (kinetic, thermal, and cosmic ray (CR) energy), energy, momentum, and mass flux, magnetic field strength and geometry, jet precession angle and period, opening-angle, and duty cycle. We investigate all of these parameters in a $10^{14}\,{\rm M}_{\odot}$ halo using high-resolution non-cosmological MHD simulations with the FIRE-2 (Feedback In Realistic Environments) stellar feedback model, conduction, and viscosity. We explore which scenarios match observational constraints and show that CR-dominated jets can most efficiently quench the central galaxy through a combination of CR pressure support and a modification of the thermal instability. Jets with most energy in mildly relativistic ($\sim$ MeV or $\sim10^{10}$ K) thermal plasma work, but require a factor $\sim 10$ larger energy input. For a fixed energy flux, jets with higher specific energy (longer cooling times) quench more effectively. For this halo size, kinetic jets are less efficient in quenching unless they have wide opening or precession angles. Magnetic fields play a minor role except when the magnetic flux reaches $\gtrsim 10^{44}$ erg s$^{-1}$ in a kinetic jet model, which causes the jet cocoon to significantly widen, and the quenching to become explosive. We conclude that the criteria for a successful jet model are an optimal energy flux and a sufficiently wide jet cocoon with long enough cooling time at the cooling radius. △ Less

Submitted 3 February, 2021; originally announced February 2021.

Comments: 29 pages, 20 figures

Resolving Structure in the Debris Disk around HD 206893 with ALMA

Authors: Ava Nederlander, A. Meredith Hughes, Anna J. Fehr, Kevin M. Flaherty, Kate Y. L. Su, Attila Moor, Eugene Chiang, Sean M. Andrews, David J. Wilner, Sebastian Marino

Abstract: Debris disks are tenuous, dusty belts surrounding main sequence stars generated by collisions between planetesimals. HD 206893 is one of only two stars known to host a directly imaged brown dwarf orbiting interior to its debris ring, in this case at a projected separation of 10.4 au. Here we resolve structure in the debris disk around HD 206893 at an angular resolution of 0.6" (24 au) and waveleng… ▽ More Debris disks are tenuous, dusty belts surrounding main sequence stars generated by collisions between planetesimals. HD 206893 is one of only two stars known to host a directly imaged brown dwarf orbiting interior to its debris ring, in this case at a projected separation of 10.4 au. Here we resolve structure in the debris disk around HD 206893 at an angular resolution of 0.6" (24 au) and wavelength of 1.3 mm with the Atacama Large Millimeter/submillimeter Array (ALMA). We observe a broad disk extending from a radius of <51 au to 194^{+13}_{-2} au. We model the disk with a continuous, gapped, and double power-law model of the surface density profile, and find strong evidence for a local minimum in the surface density distribution near a radius of 70 au, consistent with a gap in the disk with an inner radius of 63^{+8}_{-16} au and width 31^{+11}_{-7} au. Gapped structure has been observed in four other debris disks -- essentially every other radially resolved debris disk observed with sufficient angular resolution and sensitivity with ALMA -- and could be suggestive of the presence of an additional planetary-mass companion. △ Less

Submitted 21 January, 2021; originally announced January 2021.

Comments: 24 pages, 12 figures, Accepted for publication in ApJ

Planetary Embryo Collisions and the Wiggly Nature of Extreme Debris Disks

Authors: Lewis Watt, Zoë Leinhardt, Kate Su

Abstract: In this paper, we present results from a multi-stage numerical campaign to begin to explain and determine why extreme debris disk detections are rare, what types of impacts will result in extreme debris disks and what we can learn about the parameters of the collision from the extreme debris disks. We begin by simulating many giant impacts using a smoothed particle hydrodynamical code with tabulat… ▽ More In this paper, we present results from a multi-stage numerical campaign to begin to explain and determine why extreme debris disk detections are rare, what types of impacts will result in extreme debris disks and what we can learn about the parameters of the collision from the extreme debris disks. We begin by simulating many giant impacts using a smoothed particle hydrodynamical code with tabulated equations of state and track the escaping vapour from the collision. Using an $N$-body code, we simulate the spatial evolution of the vapour generated dust post-impact. We show that impacts release vapour anisotropically not isotropically as has been assumed previously and that the distribution of the resulting generated dust is dependent on the mass ratio and impact angle of the collision. In addition, we show that the anisotropic distribution of post-collision dust can cause the formation or lack of formation of the short-term variation in flux depending on the orientation of the collision with respect to the orbit around the central star. Finally, our results suggest that there is a narrow region of semi-major axis where a vapour generated disk would be observable for any significant amount of time implying that giant impacts where most of the escaping mass is in vapour would not be observed often but this does not mean that the collisions are not occurring. △ Less

Submitted 18 January, 2021; v1 submitted 13 January, 2021; originally announced January 2021.

Comments: 21 pages, 13 figures, 2 tables, accepted for publication in MNRAS

First results from SMAUG: Uncovering the Origin of the Multiphase Circumgalactic Medium with a Comparative Analysis of Idealized and Cosmological Simulations

Authors: Drummond B. Fielding, Stephanie Tonnesen, Daniel DeFelippis, Miao Li, Kung-Yi Su, Greg L. Bryan, Chang-Goo Kim, John C. Forbes, Rachel S. Somerville, Nicholas Battaglia, Evan E. Schneider, Yuan Li, Ena Choi, Christopher C. Hayward, Lars Hernquist

Abstract: We examine the properties of the circumgalactic medium (CGM) at low redshift in a range of simulated Milky Way mass halos. The sample is comprised of seven idealized simulations, an adaptive mesh refinement cosmological zoom-in simulation, and two groups of 50 halos with star forming or quiescent galaxies taken from the IllustrisTNG100 simulation. The simulations have very different setups, resolu… ▽ More We examine the properties of the circumgalactic medium (CGM) at low redshift in a range of simulated Milky Way mass halos. The sample is comprised of seven idealized simulations, an adaptive mesh refinement cosmological zoom-in simulation, and two groups of 50 halos with star forming or quiescent galaxies taken from the IllustrisTNG100 simulation. The simulations have very different setups, resolution, and feedback models, but are analyzed in a uniform manner. By comparing median radial profiles and mass distributions of CGM properties, we isolate key similarities and differences. In doing so, we advance the efforts of the SMAUG (Simulating Multiscale Astrophysics to Understand Galaxies) project that aims to understand the inherently multiscale galaxy formation process. In the cosmological simulations, the CGM exhibits nearly flat temperature distributions, and broad pressure and radial velocity distributions. In the idealized simulations, similar distributions are found in the inner CGM ($\lesssim 0.5 \, r_{\rm 200c}$) when strong galactic feedback models are employed, but the outer CGM ($\gtrsim 0.5 \, r_{\rm 200c}$) has a much less prominent cold phase, and narrower pressure and velocity distributions even in models with strong feedback. This comparative analysis demonstrates the dominant role feedback plays in shaping the inner CGM and the increased importance of cosmological effects, such as nonspherical accretion and satellite galaxies, in the outer CGM. Furthermore, our findings highlight that while cosmological simulations are required to capture the multiphase structure of the CGM at large radii, idealized simulations provide a robust framework to study how galactic feedback interacts with the inner CGM and thereby provide a reliable avenue to constrain feedback prescriptions. △ Less

Submitted 29 June, 2020; originally announced June 2020.

Comments: More information on the SMAUG project and the other first result papers can be found here: https://www.simonsfoundation.org/flatiron/center-for-computational-astrophysics/galaxy-formation/smaug/papersplash1

Dust Populations in the Iconic Vega Planetary System Resolved by ALMA

Authors: Luca Matrà, William R. F. Dent, David J. Wilner, Sebastián Marino, Mark C. Wyatt, Jonathan P. Marshall, Kate Y. L. Su, Miguel Chavez, Antonio Hales, A. Meredith Hughes, Jane S. Greaves, Stuartt A. Corder

Abstract: The Vega planetary system hosts the archetype of extrasolar Kuiper belts, and is rich in dust from the sub-au region out to 100's of au, suggesting intense dynamical activity. We present ALMA mm observations that detect and resolve the outer dust belt from the star for the first time. The interferometric visibilities show that the belt can be fit by a Gaussian model or by power-law models with a s… ▽ More The Vega planetary system hosts the archetype of extrasolar Kuiper belts, and is rich in dust from the sub-au region out to 100's of au, suggesting intense dynamical activity. We present ALMA mm observations that detect and resolve the outer dust belt from the star for the first time. The interferometric visibilities show that the belt can be fit by a Gaussian model or by power-law models with a steep inner edge (at 60-80 au). The belt is very broad, extending out to at least 150-200 au. We strongly detect the star and set a stringent upper limit to warm dust emission previously detected in the infrared. We discuss three scenarios that could explain the architecture of Vega's planetary system, including the new {ALMA} constraints: no outer planets, a chain of low-mass planets, and a single giant planet. The planet-less scenario is only feasible if the outer belt was born with the observed sharp inner edge. If instead the inner edge is currently being truncated by a planet, then the planet must be $\gtrsim$6 M$_{\oplus}$ and at $\lesssim71$ au to have cleared its chaotic zone within the system age. In the planet chain scenario, outward planet migration and inward scattering of planetesimals could produce the hot and warm dust observed in the inner regions of the system. In the single giant planet scenario, an asteroid belt could be responsible for the warm dust, and mean motion resonances with the planet could put asteroids on star-grazing orbits, producing the hot dust. △ Less

Submitted 29 June, 2020; originally announced June 2020.

Comments: 18 pages, 3 figures, Accepted for publication in ApJ

Mid-infrared Studies of HD 113766 and HD 172555: Assessing Variability in the Terrestrial Zone of Young Exoplanetary Systems

Authors: Kate Y. L. Su, George H. Rieke, Carl Melis, Alan P. Jackson, Paul S. Smith, Huan Y. A. Meng, Andras Gaspar

Abstract: We present multi-epoch infrared photometry and spectroscopy obtained with warm Spitzer, Subaru and SOFIA to assess variability for the young ($\sim$20 Myr) and dusty debris systems around HD 172555 and HD 113766A. No variations (within 0.5%) were found for the former at either 3.6 or 4.5 $μ$m, while significant non-periodic variations (peak-to-peak of $\sim$10-15% relative to the primary star) wer… ▽ More We present multi-epoch infrared photometry and spectroscopy obtained with warm Spitzer, Subaru and SOFIA to assess variability for the young ($\sim$20 Myr) and dusty debris systems around HD 172555 and HD 113766A. No variations (within 0.5%) were found for the former at either 3.6 or 4.5 $μ$m, while significant non-periodic variations (peak-to-peak of $\sim$10-15% relative to the primary star) were detected for the latter. Relative to the Spitzer IRS spectra taken in 2004, multi-epoch mid-infrared spectra reveal no change in either the shape of the prominent 10 $μ$m solid-state features or the overall flux levels (no more than 20%) for both systems, corroborating that the population of sub-$μ$m-sized grains that produce the pronounced solid-state features is stable over a decadal timescale. We suggest that these sub-$μ$m-sized grains were initially generated in an optically thick clump of debris of mm-sized vapor condensates resulting from a recent violent impact between large asteroidal or planetary bodies. Because of the shielding from the stellar photons provided by this clump, intense collisions led to an over-production of fine grains that would otherwise be ejected from the system by radiation pressure. As the clump is sheared by its orbital motion and becomes optically thin, a population of very fine grains could remain in stable orbits until Poynting-Robertson drag slowly spirals them into the star. We further suggest that the 3-5 $μ$m disk variation around HD 113766A is consistent with a clump/arc of such fine grains on a modestly eccentric orbit in its terrestrial zone. △ Less

Submitted 11 June, 2020; originally announced June 2020.

Comments: to be published in ApJ

No significant correlation between radial velocity planet presence and debris disc properties

Authors: Ben Yelverton, Grant M. Kennedy, Kate Y. L. Su

Abstract: We investigate whether the tentative correlation between planets and debris discs which has been previously identified can be confirmed at high significance. We compile a sample of 201 stars with known planets and existing far infrared observations. The sample is larger than those studied previously since we include targets from an unpublished Herschel survey of planet hosts. We use spectral energ… ▽ More We investigate whether the tentative correlation between planets and debris discs which has been previously identified can be confirmed at high significance. We compile a sample of 201 stars with known planets and existing far infrared observations. The sample is larger than those studied previously since we include targets from an unpublished Herschel survey of planet hosts. We use spectral energy distribution modelling to characterise Kuiper belt analogue debris discs within the sample, then compare the properties of the discs against a control sample of 294 stars without known planets. Survival analysis suggests that there is a significant ($p\sim 0.002$) difference between the disc fractional luminosity distributions of the two samples. However, this is largely a result of the fact that the control sample contains a higher proportion of close binaries and of later-type stars; both of these factors are known to reduce disc detection rates. Considering only Sun-like stars without close binary companions in each sample greatly reduces the significance of the difference ($p\sim 0.3$). We also find no evidence for a difference in the disc fractional luminosities of stars hosting planets more or less massive than Saturn ($p\sim 0.9$). Finally, we find that the planet hosts have cooler discs than the control stars, but this is likely a detection bias, since the warmest discs in the control sample are also the faintest, and would thus be undetectable around the more distant planet hosts. Considering only discs in each sample that could have been detected around a typical planet host, we find $p\sim 0.07$ for the temperatures. △ Less

Submitted 7 May, 2020; originally announced May 2020.

Comments: MNRAS, in press

Near-infrared variability in dusty white dwarfs: tracing the accretion of planetary material

Authors: Laura K. Rogers, Siyi Xu, Amy Bonsor, Simon Hodgkin, Kate Y. L. Su, Ted von Hippel, Michael Jura

Abstract: The inwards scattering of planetesimals towards white dwarfs is expected to be a stochastic process with variability on human time-scales. The planetesimals tidally disrupt at the Roche radius, producing dusty debris detectable as excess infrared emission. When sufficiently close to the white dwarf, this debris sublimates and accretes on to the white dwarf and pollutes its atmosphere. Studying thi… ▽ More The inwards scattering of planetesimals towards white dwarfs is expected to be a stochastic process with variability on human time-scales. The planetesimals tidally disrupt at the Roche radius, producing dusty debris detectable as excess infrared emission. When sufficiently close to the white dwarf, this debris sublimates and accretes on to the white dwarf and pollutes its atmosphere. Studying this infrared emission around polluted white dwarfs can reveal how this planetary material arrives in their atmospheres. We report a near-infrared monitoring campaign of 34 white dwarfs with infrared excesses with the aim to search for variability in the dust emission. Time series photometry of these white dwarfs from the United Kingdom Infrared Telescope (Wide Field Camera) in the J, H and K bands were obtained over baselines of up to three years. We find no statistically significant variation in the dust emission in all three near-infrared bands. Specifically, we can rule out variability at ~1.3% for the 13 white dwarfs brighter than 16th mag in K band, and at ~10% for the 32 white dwarfs brighter than 18th mag over time-scales of three years. Although to date two white dwarfs, SDSS J095904.69-020047.6 and WD 1226+110, have shown K band variability, in our sample we see no evidence of new K band variability at these levels. One interpretation is that the tidal disruption events which lead to large variabilities are rare, occur on short time-scales, and after a few years the white dwarfs return to being stable in the near-infrared. △ Less

Submitted 30 March, 2020; originally announced March 2020.

Comments: 14 pages + supplementary figures, 9 figures, accepted for publication in MNRAS

The HOSTS survey for exozodiacal dust: Observational results from the complete survey

Authors: Steve Ertel, Denis Defrère, Philip M. Hinz, Bertrand Mennesson, Grant M. Kennedy, William C. Danchi, Christopher Gelino, John M. Hill, William F. Hoffmann, Johan Mazoyer, George Rieke, Andrew Shannon, Karl Stapelfeldt, Eckhart Spalding, Jordan M. Stone, Amali Vaz, Alycia J. Weinberger, Phil Willems, Olivier Absil, Paul Arbo, Vanessa P. Bailey, Charles Beichman, Geoffrey Bryden, Elwood C. Downey, Olivier Durney , et al. (21 additional authors not shown)

Abstract: The Large Binocular Telescope Interferometer (LBTI) enables nulling interferometric observations across the N band (8 to 13 um) to suppress a star's bright light and probe for faint circumstellar emission. We present and statistically analyze the results from the LBTI/HOSTS (Hunt for Observable Signatures of Terrestrial Systems) survey for exozodiacal dust. By comparing our measurements to model p… ▽ More The Large Binocular Telescope Interferometer (LBTI) enables nulling interferometric observations across the N band (8 to 13 um) to suppress a star's bright light and probe for faint circumstellar emission. We present and statistically analyze the results from the LBTI/HOSTS (Hunt for Observable Signatures of Terrestrial Systems) survey for exozodiacal dust. By comparing our measurements to model predictions based on the Solar zodiacal dust in the N band, we estimate a 1 sigma median sensitivity of 23 zodis for early type stars and 48 zodis for Sun-like stars, where 1 zodi is the surface density of habitable zone (HZ) dust in the Solar system. Of the 38 stars observed, 10 show significant excess. A clear correlation of our detections with the presence of cold dust in the systems was found, but none with the stellar spectral type or age. The majority of Sun-like stars have relatively low HZ dust levels (best-fit median: 3 zodis, 1 sigma upper limit: 9 zodis, 95% confidence: 27 zodis based on our N band measurements), while ~20% are significantly more dusty. The Solar system's HZ dust content is consistent with being typical. Our median HZ dust level would not be a major limitation to the direct imaging search for Earth-like exoplanets, but more precise constraints are still required, in particular to evaluate the impact of exozodiacal dust for the spectroscopic characterization of imaged exo-Earth candidates. △ Less

Submitted 6 March, 2020; originally announced March 2020.

Comments: accepted for publication in AJ