Showing 1–50 of 506 results for author: Nelson, D

Zooming in on the Circumgalactic Medium with GIBLE: Tracing the Origin and Evolution of Cold Clouds

Authors: Rahul Ramesh, Dylan Nelson, Drummond Fielding, Marcus Brüggen

Abstract: We use the GIBLE suite of cosmological zoom-in simulations of Milky Way-like galaxies with additional super-Lagrangian refinement in the circumgalactic medium (CGM) to quantify the origin and evolution of CGM cold gas clouds. The origin of $z$\,$=$\,$0$ clouds can be traced back to recent ($\lesssim$\,$2$\,Gyr) outflows from the central galaxy ($\sim$\,45\,$\%$), condensation out of the hot phase… ▽ More We use the GIBLE suite of cosmological zoom-in simulations of Milky Way-like galaxies with additional super-Lagrangian refinement in the circumgalactic medium (CGM) to quantify the origin and evolution of CGM cold gas clouds. The origin of $z$\,$=$\,$0$ clouds can be traced back to recent ($\lesssim$\,$2$\,Gyr) outflows from the central galaxy ($\sim$\,45\,$\%$), condensation out of the hot phase of the CGM in the same time frame ($\sim$\,45\,$\%$), and to a lesser degree to satellite galaxies ($\lesssim$\,5\,$\%$). We find that in-situ condensation results from rapid cooling around local over-densities primarily seeded by the dissolution of the previous generation of clouds into the hot halo. About $\lesssim$\,10\,$\%$ of the cloud population is long lived, with their progenitors having already assembled $\sim$\,$2$\,Gyr ago. Collective cloud-cloud dynamics are crucial to their evolution, with coalescence and fragmentation events occurring frequently ($\gtrsim$\,20\,Gyr$^{-1}$). These interactions are modulated by non-vanishing pressure imbalances between clouds and their interface layers. The gas content of clouds is in a constant state of flux, with clouds and their surroundings exchanging mass at a rate of \mbox{$\gtrsim$\,$10^3$\,M$_\odot$\,Myr$^{-1}$}, depending on cloud relative velocity and interface vorticity. Furthermore, we find that a net magnetic tension force acting against the density gradient is capable of inhibiting cloud-background mixing. Our results show that capturing the distinct origins of cool CGM clouds, together with their physical evolution, requires high-resolution, cosmological galaxy formation simulations with both stellar and supermassive black hole feedback-driven outflows. △ Less

Submitted 28 June, 2024; originally announced July 2024.

Comments: Submitted to A&A

LBT SHARK-VIS Observes a Major Resurfacing Event on Io

Authors: Al Conrad, Fernando Pedichini, Gianluca Li Causi, Simone Antoniucci, Imke de Pater, Ashley Gerard Davies, Katherine de Kleer, Roberto Piazzesi, Vincenzo Testa, Piero Vaccari, Martina Vicinanza, Jennifer Power, Steve Ertel, Joseph C. Shields, Sam Ragland, Fabrizio Giorgi, Stuart M. Jefferies, Douglas Hope, Jason Perry, David A. Williams, David M. Nelson

Abstract: Since volcanic activity was first discovered on Io from Voyager images in 1979, changes on Io's surface have been monitored from both spacecraft and ground-based telescopes. Here, we present the highest spatial resolution images of Io ever obtained from a ground-based telescope. These images, acquired by the SHARK-VIS instrument on the Large Binocular Telescope, show evidence of a major resurfacin… ▽ More Since volcanic activity was first discovered on Io from Voyager images in 1979, changes on Io's surface have been monitored from both spacecraft and ground-based telescopes. Here, we present the highest spatial resolution images of Io ever obtained from a ground-based telescope. These images, acquired by the SHARK-VIS instrument on the Large Binocular Telescope, show evidence of a major resurfacing event on Io's trailing hemisphere. When compared to the most recent spacecraft images, the SHARK-VIS images show that a plume deposit from a powerful eruption at Pillan Patera has covered part of the long-lived Pele plume deposit. Although this type of resurfacing event may be common on Io, few have been detected due to the rarity of spacecraft visits and the previously low spatial resolution available from Earth-based telescopes. The SHARK-VIS instrument ushers in a new era of high resolution imaging of Io's surface using adaptive optics at visible wavelengths. △ Less

Submitted 29 May, 2024; originally announced May 2024.

Comments: 15 pages, 4 figures

Defect ground states for liquid crystals on cones and hyperbolic cones

Authors: Farzan Vafa, Grace H. Zhang, David R. Nelson

Abstract: This contribution is intended for Journal of Physics A: Mathematical and Theoretical Special issue on Non-equilibrium Dynamics in Complex Systems: Celebrating the Contributions of Uwe Täuber on his 60th Birthday. Cones with orientational order in the local tangent plane provide a soft matter analog of the Aharonov-Bohm effect. We first review recent work on 2D liquid crystals with $p$-fold rotatio… ▽ More This contribution is intended for Journal of Physics A: Mathematical and Theoretical Special issue on Non-equilibrium Dynamics in Complex Systems: Celebrating the Contributions of Uwe Täuber on his 60th Birthday. Cones with orientational order in the local tangent plane provide a soft matter analog of the Aharonov-Bohm effect. We first review recent work on 2D liquid crystals with $p$-fold rotational symmetry on cones. We exploit an analogy with electrostatics to determine the ground state as a function of both the cone angle and the liquid crystal symmetry for both free and tangential boundary conditions applied at the cone base. There is an effective topological charge $-χ$ at the apex, where $2πχ$ is the deficit angle. The ground states are in general frustrated due to parallel transport along the azimuthal direction on the cone. In the case of tangential boundary conditions, the ground state changes as a function of $χ$, where the apex absorbs and emits quantized defect charges, with intricate dependence on both the deficit angle and $p$. We check our predictions numerically for a set of commensurate cone angles, whose surfaces can be polygonized as a perfect triangular or square mesh, and find excellent agreement. Cones with both free and tangential boundary conditions can also exhibit metastable states distinguished by quantized screening of the apex charge. We also present preliminary work on hyperbolic cones, where the Gaussian curvature singularity is negative at the apex. When free boundary conditions are applied at the base, the ground states are characterized by an effective topological charge at the apex, similarly to the case of conventional cones. However, when tangential boundary conditions are applied, decreasing the deficit angle (which is now negative) induces neutral defect pair nucleation at the apex followed by emission of a positive defect. △ Less

Submitted 6 May, 2024; originally announced May 2024.

Comments: 22 pages, 19 figures, Journal of Physics A: Mathematical and Theoretical Special issue on Non-equilibrium Dynamics in Complex Systems: Celebrating the Contributions of Uwe Täuber on his 60th Birthday

MAGAZ3NE: Massive, Extremely Dusty Galaxies at $z\sim2$ Lead to Photometric Overestimation of Number Densities of the Most Massive Galaxies at $3<z<4$

Authors: Ben Forrest, M. C. Cooper, Adam Muzzin, Gillian Wilson, Danilo Marchesini, Ian McConachie, Percy Gomez, Marianna Annunziatella, Z. Cemile Marsan, Joey Braspenning, Wenjun Chang, Gabriella de Lucia, Fabio Fontanot, Michaela Hirschmann, Dylan Nelson, Annalisa Pillepich, Joop Schaye, Stephanie M. Urbano Stawinski, Mauro Stefanon, Lizhi Xie

Abstract: We present rest-frame optical spectra from Keck/MOSFIRE and Keck/NIRES of 16 candidate ultramassive galaxies targeted as part of the Massive Ancient Galaxies at $z>3$ Near-Infrared (MAGAZ3NE) Survey. These candidates were selected to have photometric redshifts $3\lesssim z_{\rm phot}<4$, photometric stellar masses log($M$/M$_\odot$)$>11.7$, and well-sampled photometric spectral energy distribution… ▽ More We present rest-frame optical spectra from Keck/MOSFIRE and Keck/NIRES of 16 candidate ultramassive galaxies targeted as part of the Massive Ancient Galaxies at $z>3$ Near-Infrared (MAGAZ3NE) Survey. These candidates were selected to have photometric redshifts $3\lesssim z_{\rm phot}<4$, photometric stellar masses log($M$/M$_\odot$)$>11.7$, and well-sampled photometric spectral energy distributions (SEDs) from the UltraVISTA and VIDEO surveys. In contrast to previous spectroscopic observations of blue star-forming and post-starburst ultramassive galaxies, candidates in this sample have very red SEDs implying significant dust attenuation, old stellar ages, and/or active galactic nuclei (AGN). Of these galaxies, eight are revealed to be heavily dust-obscured $2.0<z<2.7$ galaxies with strong emission lines, some showing broad features indicative of AGN, three are Type I AGN hosts at $z>3$, one is a $z\sim1.2$ dusty galaxy, and four galaxies do not have a confirmed spectroscopic redshift. In fact, none of the sample has |$z_{\rm spec}-z_{\rm phot}$|$<0.5$, suggesting difficulties for photometric redshift programs in fitting similarly red SEDs. The prevalence of these red interloper galaxies suggests that the number densities of high-mass galaxies are overestimated at $z\gtrsim3$ in large photometric surveys, helping to resolve the `impossibly early galaxy problem' and leading to much better agreement with cosmological galaxy simulations. A more complete spectroscopic survey of ultramassive galaxies is required to pin down the uncertainties on their number densities in the early universe. △ Less

Submitted 29 April, 2024; originally announced April 2024.

Comments: Submitted to ApJ. Comments welcome!

Local integral transforms and global spectral decomposition

Authors: Valentin Blomer, Subhajit Jana, Paul D. Nelson

Abstract: We establish an explicit global spectral decomposition of shifted convolution sums and the second moment of automorphic L-functions for Maass forms with explicit integral transforms as well as explicit inversion formulae over every local field. We establish an explicit global spectral decomposition of shifted convolution sums and the second moment of automorphic L-functions for Maass forms with explicit integral transforms as well as explicit inversion formulae over every local field. △ Less

Submitted 16 April, 2024; originally announced April 2024.

Comments: 41 pages

Local analysis of the Kuznetsov formula and the density conjecture

Authors: Edgar Assing, Valentin Blomer, Paul D. Nelson

Abstract: We prove Sarnak's spherical density conjecture for the principal congruence subgroup of SL(n, Z) of arbitrary level. Applications include a complete version of Sarnak's optimal lifting conjecture for principal congruence subgroups of SL(n, Z), as well as a transfer of the density theorem to certain co-compact situations. The main ingredients are new lower bounds for Whittaker functions and strong… ▽ More We prove Sarnak's spherical density conjecture for the principal congruence subgroup of SL(n, Z) of arbitrary level. Applications include a complete version of Sarnak's optimal lifting conjecture for principal congruence subgroups of SL(n, Z), as well as a transfer of the density theorem to certain co-compact situations. The main ingredients are new lower bounds for Whittaker functions and strong estimates for the cardinality of ramified Kloosterman sets. △ Less

Submitted 8 April, 2024; originally announced April 2024.

Comments: 48 pages

Zooming in on the Circumgalactic Medium with GIBLE: the Topology and Draping of Magnetic Fields around Cold Clouds

Authors: Rahul Ramesh, Dylan Nelson, Drummond Fielding, Marcus Brüggen

Abstract: We use a cosmological zoom-in simulation of a Milky Way-like galaxy to study and quantify the topology of magnetic field lines around cold gas clouds in the circumgalactic medium (CGM). This simulation is a new addition to Project GIBLE, a suite of cosmological magnetohydrodynamical simulations of galaxy formation with preferential super-Lagrangian refinement in the CGM, reaching an unprecedented… ▽ More We use a cosmological zoom-in simulation of a Milky Way-like galaxy to study and quantify the topology of magnetic field lines around cold gas clouds in the circumgalactic medium (CGM). This simulation is a new addition to Project GIBLE, a suite of cosmological magnetohydrodynamical simulations of galaxy formation with preferential super-Lagrangian refinement in the CGM, reaching an unprecedented (CGM) gas mass resolution of $\sim$ $225$ M$_\odot$. To maximize statistics and resolution, we focus on a sample of $\sim$ $200$ clouds with masses of $\sim$ $10^6$ M$_\odot$. The topology of magnetic field lines around clouds is diverse, from threading to draping, and there is large variation in the magnetic curvature ($κ$) within cloud-background interfaces. We typically find little variation of $κ$ between upstream and downstream cloud faces, implying that strongly draped configurations are rare. In addition, $κ$ correlates strongly with multiple properties of the interface and the ambient background, including cloud overdensity and relative velocity, suggesting that cloud properties impact the topology of interface magnetic fields. △ Less

Submitted 1 April, 2024; originally announced April 2024.

Comments: Accepted for publication (A&A Letters); 5 pages, 4 figures

Journal ref: A&A 684, L16 (2024)

Atacama Large Aperture Submillimeter Telescope (AtLAST) Science: Resolving the Hot and Ionized Universe through the Sunyaev-Zeldovich effect

Authors: Luca Di Mascolo, Yvette Perrott, Tony Mroczkowski, Stefano Andreon, Stefano Ettori, Aurora Simionescu, Srinivasan Raghunathan, Joshiwa van Marrewijk, Claudia Cicone, Minju Lee, Dylan Nelson, Laura Sommovigo, Mark Booth, Pamela Klaassen, Paola Andreani, Martin A. Cordiner, Doug Johnstone, Eelco van Kampen, Daizhong Liu, Thomas J. Maccarone, Thomas W. Morris, Amélie Saintonge, Matthew Smith, Alexander E. Thelen, Sven Wedemeyer

Abstract: An omnipresent feature of the multi-phase ``cosmic web'' is that warm/hot (>$10^5$ K) ionized gas pervades it. This gas constitutes a relevant contribution to the overall universal matter budget across multiple scales, from the several tens of Mpc-scale IGM filaments, to the Mpc ICM, all the way down to the CGM surrounding individual galaxies, on scales from ~1 kpc up to their respective virial ra… ▽ More An omnipresent feature of the multi-phase ``cosmic web'' is that warm/hot (>$10^5$ K) ionized gas pervades it. This gas constitutes a relevant contribution to the overall universal matter budget across multiple scales, from the several tens of Mpc-scale IGM filaments, to the Mpc ICM, all the way down to the CGM surrounding individual galaxies, on scales from ~1 kpc up to their respective virial radii (~100 kpc). The study of the hot baryonic component of cosmic matter density represents a powerful means for constraining the intertwined evolution of galactic populations and large-scale cosmological structures, for tracing the matter assembly in the Universe and its thermal history. To this end, the SZ effect provides the ideal observational tool for measurements out to the beginnings of structure formation. The SZ effect is caused by the scattering of the photons from the cosmic microwave background off the hot electrons embedded within cosmic structures, and provides a redshift-independent perspective on the thermal and kinematic properties of the warm/hot gas. Still, current and future (sub)mm facilities have been providing only a partial view of the SZ Universe due to any combination of: limited angular resolution, spectral coverage, field of view, spatial dynamic range, sensitivity. In this paper, we motivate the development of a wide-field, broad-band, multi-chroic continuum instrument for the Atacama Large Aperture Submillimeter Telescope (AtLAST) by identifying the scientific drivers that will deepen our understanding of the complex thermal evolution of cosmic structures. On a technical side, this will necessarily require efficient multi-wavelength mapping of the SZ signal with an unprecedented spatial dynamic range (from arcsecond to degree scales) and we employ theoretical forecasts to determine the key instrumental constraints for achieving our goals. [abridged] △ Less

Submitted 1 March, 2024; originally announced March 2024.

Comments: 29 pages, 17 figures, 1 table. Submitted to Open Research Europe as part of the AtLAST Design Study collection: https://open-research-europe.ec.europa.eu/collections/atlast/about. Comments are welcome

scida: scalable analysis for scientific big data

Authors: Chris Byrohl, Dylan Nelson

Abstract: scida is a Python package for reading and analyzing large scientific data sets with support for various cosmological and galaxy formation simulations out-of-the-box. Data access is provided through a hierarchical dictionary-like data structure after a simple load() function. Using the dask library for scalable, parallel and out-of-core computation, all computation requests from a user session are… ▽ More scida is a Python package for reading and analyzing large scientific data sets with support for various cosmological and galaxy formation simulations out-of-the-box. Data access is provided through a hierarchical dictionary-like data structure after a simple load() function. Using the dask library for scalable, parallel and out-of-core computation, all computation requests from a user session are first collected in a task graph. Arbitrary custom analysis, as well as all available dask (array) operations, can be performed. The subsequent computation is executed only upon request, on a target resource (e.g. a HPC cluster). △ Less

Submitted 27 February, 2024; originally announced February 2024.

Comments: recommended for acceptance in the Journal of Open Source Software; open-source development at https://github.com/cbyrohl/scida

On the relationship between speech and hearing

Authors: Srinivasan Umesh, Leon Cohen, Douglas Nelson

Abstract: We present a framework for experimentally linking speech production and hearing. Using this approach, we describe experimental results, that lead to the concept that sounds made by different individuals and perceived to be the same can be transformed into each other by a "speech scale". The speech scale is empirically determined using only speech data. We show the similarity of the speech scale to… ▽ More We present a framework for experimentally linking speech production and hearing. Using this approach, we describe experimental results, that lead to the concept that sounds made by different individuals and perceived to be the same can be transformed into each other by a "speech scale". The speech scale is empirically determined using only speech data. We show the similarity of the speech scale to the MEL scale of Stevens and Volkmann, which was derived only from hearing experiments. We thus experimentally link speech production and hearing. △ Less

Submitted 19 February, 2024; originally announced February 2024.

Kerr optical frequency division with integrated photonics for stable microwave and mmWave generation

Authors: Shuman Sun, Mark W. Harrington, Fatemehsadat Tabatabaei, Samin Hanifi, Kaikai Liu, Jiawei Wang, Beichen Wang, Zijiao Yang, Ruxuan Liu, Jesse S. Morgan, Steven M. Bowers, Paul A. Morton, Karl D. Nelson, Andreas Beling, Daniel J. Blumenthal, Xu Yi

Abstract: Optical frequency division (OFD) has revolutionized microwave and mmWave generation and set spectral purity records owing to its unique capability to transfer high fractional stability from optical to electronic frequencies. Recently, rapid developments in integrated optical reference cavities and microresonator-based optical frequency combs (microcombs) have created a path to transform OFD techno… ▽ More Optical frequency division (OFD) has revolutionized microwave and mmWave generation and set spectral purity records owing to its unique capability to transfer high fractional stability from optical to electronic frequencies. Recently, rapid developments in integrated optical reference cavities and microresonator-based optical frequency combs (microcombs) have created a path to transform OFD technology to chip scale. Here, we demonstrate an ultra-low phase noise mmWave oscillator by leveraging integrated photonic components and Kerr optical frequency division. The oscillator derives its stability from an integrated CMOS-compatible SiN coil cavity, and the optical frequency division is achieved spontaneously through Kerr interaction between the injected reference lasers and soliton microcombs in the integrated SiN microresonator. Besides achieving record-low phase noise for integrated mmWave oscillators, our demonstration greatly simplifies the implementation of integrated OFD oscillators and could be useful in applications of Radar, spectroscopy, and astronomy. △ Less

Submitted 18 February, 2024; originally announced February 2024.

Comments: 8 pages, 5 figures

Inverse Designed WS2 Planar Chiral Metasurface with Geometric Phase

Authors: Jaegang Jo, Sangbin Lee, Munseong Bae, Damian Nelson, Kenneth B. Crozier, Nanfang Yu, Haejun Chung, Sejeong Kim

Abstract: Increasing attention is being paid to chiral metasurfaces due to their ability to selectively manipulate right-hand circularly polarized light or left-hand circularly polarized light. The thin nature of metasurfaces, however, poses a challenge in creating a device with effective phase modulation. Plasmonic chiral metasurfaces have attempted to address this issue by increasing light-matter interact… ▽ More Increasing attention is being paid to chiral metasurfaces due to their ability to selectively manipulate right-hand circularly polarized light or left-hand circularly polarized light. The thin nature of metasurfaces, however, poses a challenge in creating a device with effective phase modulation. Plasmonic chiral metasurfaces have attempted to address this issue by increasing light-matter interaction, but they suffer from metallic loss. Dielectric metasurfaces made from high index materials enable phase modulation while being thin. Very few materials, however, have high refractive index and low loss at visible wavelengths. Recently, some 2D materials have been shown to exhibit high refractive index and low loss in the visible wavelengths, positioning them as promising platform for meta-optics. This study introduces and details a planar chiral metasurface with geometric phase composed of WS2 meta-units. By employing adjoint optimization techniques, we achieved broadband circular dichroism. △ Less

Submitted 6 February, 2024; originally announced February 2024.

Comments: 9 pages, 5 figures

Radio relics in massive galaxy cluster mergers in the TNG-Cluster simulation

Authors: W. Lee, A. Pillepich, J. ZuHone, D. Nelson, M. J. Jee, D. Nagai, K. Finner

Abstract: Radio relics are diffuse synchrotron sources in the outskirts of merging galaxy clusters energized by the merger shocks. In this paper, we present an overview of the radio relics in massive cluster mergers identified in the new TNG-Cluster simulation. This is a suite of magnetohydrodynamical cosmological zoom-in simulations of 352 massive galaxy clusters with… ▽ More Radio relics are diffuse synchrotron sources in the outskirts of merging galaxy clusters energized by the merger shocks. In this paper, we present an overview of the radio relics in massive cluster mergers identified in the new TNG-Cluster simulation. This is a suite of magnetohydrodynamical cosmological zoom-in simulations of 352 massive galaxy clusters with $M_{\rm 500c}= 10^{14.0-15.3}\rm~M_{\odot}$ sampled from a 1 Gpc-size cosmological box. The simulations are performed using the moving-mesh code AREPO with the galaxy formation model and high numerical resolution consistent with the TNG300 run of the IllustrisTNG series. We post-process the shock properties obtained from the on-the-fly shock finder to estimate the diffuse radio emission generated by cosmological shockwaves for a total of $\sim300$ radio relics at redshift $z=0-1$. TNG-Cluster returns a variety of radio relics with diverse morphologies, encompassing textbook examples of double radio relics, single relics, and ``inverted" radio relics that are convex to the cluster center. Moreover, the simulated radio relics reproduce both the abundance and statistical relations of observed relics. We find that extremely large radio relics ($>$ 2 Mpc) are predominantly produced in massive cluster mergers with $M_{\rm 500c}\gtrsim8\times10^{14}~\rm~M_{\odot}$. This underscores the significance of simulating massive mergers to study giant radio relics similar to those found in observations. We release a library of radio relics from the TNG-Cluster simulation, which will serve as a crucial reference for upcoming next-generation surveys. △ Less

Submitted 11 March, 2024; v1 submitted 10 November, 2023; originally announced November 2023.

Comments: 24 pages, 16 figures, 2 appendices, Accepted for publication in A&A. See the TNG-Cluster website at www.tng-project.org/cluster/

An Atlas of Gas Motions in the TNG-Cluster Simulation: from Cluster Cores to the Outskirts

Authors: Mohammadreza Ayromlou, Dylan Nelson, Annalisa Pillepich, Eric Rohr, Nhut Truong, Yuan Li, Aurora Simionescu, Katrin Lehle, Wonki Lee

Abstract: Galaxy clusters are unique laboratories for studying astrophysical processes and their impact on gas kinematics. Despite their importance, the full complexity of gas motion within and around clusters remains poorly known. This paper is part of a series presenting first results from the new TNG-Cluster simulation, a suite of 352 massive clusters including the full cosmological context, mergers, acc… ▽ More Galaxy clusters are unique laboratories for studying astrophysical processes and their impact on gas kinematics. Despite their importance, the full complexity of gas motion within and around clusters remains poorly known. This paper is part of a series presenting first results from the new TNG-Cluster simulation, a suite of 352 massive clusters including the full cosmological context, mergers, accretion, baryonic processes, feedback, and magnetic fields. Studying the dynamics and coherence of gas flows, we find that gas motions in cluster cores and intermediate regions are largely balanced between inflows and outflows, exhibiting a Gaussian distribution centered at zero velocity. In the outskirts, even the net velocity distribution becomes asymmetric, featuring a double peak where the second peak reflects cosmic accretion. Across all cluster regions, the resulting net flow distribution reveals complex gas dynamics. These are strongly correlated with halo properties: at a given total cluster mass, unrelaxed, late-forming halos with less massive black holes and lower accretion rates exhibit a more dynamic behavior. Our analysis shows no clear relationship between line-of-sight and radial gas velocities, suggesting that line-of-sight velocity alone is insufficient to distinguish between inflowing and outflowing gas. Additional properties, such as temperature, can help break this degeneracy. A velocity structure function (VSF) analysis indicates more coherent gas motion in the outskirts and more disturbed kinematics towards halo centers. In all cluster regions, the VSF shows a slope close to the theoretical models of Kolmogorov (1/3), except within 50 kpc of the cluster cores, where the slope is significantly steeper. The outcome of TNG-Cluster broadly aligns with observations of the VSF of multiphase gas across different scales in galaxy clusters, ranging from 1 kpc to Megaparsec scales. △ Less

Submitted 10 November, 2023; originally announced November 2023.

Comments: Submitted to A&A. See the TNG-Cluster website at https://www.tng-project.org/cluster/

Introducing the TNG-Cluster Simulation: overview and physical properties of the gaseous intracluster medium

Authors: Dylan Nelson, Annalisa Pillepich, Mohammadreza Ayromlou, Wonki Lee, Katrin Lehle, Eric Rohr, Nhut Truong

Abstract: We introduce the new TNG-Cluster project, an addition to the IllustrisTNG suite of cosmological magnetohydrodynamical simulations of galaxy formation. Our objective is to significantly increase the statistical sampling of the most massive and rare objects in the Universe: galaxy clusters with log(M_200c / Msun) > 14.3 - 15.4 at z=0. To do so, we re-simulate 352 cluster regions drawn from a 1 Gpc v… ▽ More We introduce the new TNG-Cluster project, an addition to the IllustrisTNG suite of cosmological magnetohydrodynamical simulations of galaxy formation. Our objective is to significantly increase the statistical sampling of the most massive and rare objects in the Universe: galaxy clusters with log(M_200c / Msun) > 14.3 - 15.4 at z=0. To do so, we re-simulate 352 cluster regions drawn from a 1 Gpc volume, thirty-six times larger than TNG300, keeping entirely fixed the IllustrisTNG physical model as well as the numerical resolution. This new sample of hundreds of massive galaxy clusters enables studies of the assembly of high-mass ellipticals and their supermassive black holes (SMBHs), brightest cluster galaxies (BCGs), satellite galaxy evolution and environmental processes, jellyfish galaxies, intracluster medium (ICM) properties, cooling and active galactic nuclei (AGN) feedback, mergers and relaxedness, magnetic field amplification, chemical enrichment, and the galaxy-halo connection at the high-mass end, with observables from the optical to radio synchrotron and the Sunyaev-Zeldovich (SZ) effect, to X-ray emission, as well as their cosmological applications. We present an overview of the simulation, the cluster sample, selected comparisons to data, and a first look at the diversity and physical properties of our simulated clusters and their hot ICM. △ Less

Submitted 6 March, 2024; v1 submitted 10 November, 2023; originally announced November 2023.

Comments: A&A accepted. See companion papers (Ayromlou, Lee, Lehle, Rohr, Truong). Additional information and visuals are available on the TNG-Cluster website at https://www.tng-project.org/cluster/

The hot circumgalactic media of massive cluster satellites in the TNG-Cluster simulation: existence and detectability

Authors: Eric Rohr, Annalisa Pillepich, Dylan Nelson, Mohammadreza Ayromlou, Elad Zinger

Abstract: The most massive galaxy clusters in the Universe host hundreds of massive satellite galaxies~$\mstar\sim10^{10-12.5} msun$, but it is unclear if these satellites are able to retain their own gaseous halos. We analyze the evolution of $\approx90,000$ satellites of stellar mass $\sim10^{9-12.5} msun$ around 352 galaxy clusters of mass $\mvir\sim10^{14.3-15.4} msun$ at $z=0$ from the new TNG-Cluster… ▽ More The most massive galaxy clusters in the Universe host hundreds of massive satellite galaxies~$\mstar\sim10^{10-12.5} msun$, but it is unclear if these satellites are able to retain their own gaseous halos. We analyze the evolution of $\approx90,000$ satellites of stellar mass $\sim10^{9-12.5} msun$ around 352 galaxy clusters of mass $\mvir\sim10^{14.3-15.4} msun$ at $z=0$ from the new TNG-Cluster suite of cosmological magneto-hydrodynamical galaxy cluster simulations. The number of massive satellites per host increases with host mass, and the mass--richness relation broadly agrees with observations. A halo of mass $\mvirhost\sim10^{14.5} (10^{15}) msun$ hosts $\sim100 (300)$ satellites today. Only a minority of satellites retain some gas, hot or cold, and this fraction increases with stellar mass. Lower-mass satellites $\sim10^{9-10} msun$ are more likely to retain part of their cold interstellar medium, consistent with ram pressure preferentially removing hot extended gas first. At higher stellar masses $\sim10^{10.5-12.5} msun$, the fraction of gas-rich satellites increases to unity, and nearly all satellites retain a portion of their hot, spatially extended circumgalactic medium (CGM), despite the ejective activity of their supermassive black holes. According to TNG-Cluster, the CGM of these gaseous satellites can be seen in soft X-ray emission (0.5-2.0 keV) that is $\gtrsim10$~times brighter than the local background. This X-ray surface brightness excess around satellites extends to $\approx30-100$ kpc, and is strongest for galaxies with higher stellar masses and larger host-centric distances. Approximately $10$~percent of the soft X-ray emission in cluster outskirts $\approx0.75-1.5\rvir$ originates from satellites. The CGM of member galaxies reflects the dynamics of cluster-satellite interactions and contributes to the observationally inferred properties of the intracluster medium. △ Less

Submitted 16 April, 2024; v1 submitted 10 November, 2023; originally announced November 2023.

Comments: Submitted to A&A, comments welcome. See companion papers by Nelson et al., Lehle et al., Truong et al., Lee et al., and Ayromlou et al., and see the TNG-Cluster website at www.tng-project.org/cluster/. Reuploaded after acceptance for publication

Journal ref: A&A 686, A86 (2024)

X-ray-inferred kinematics of the core ICM in Perseus-like clusters: Insights from the TNG-Cluster simulation

Authors: Nhut Truong, Annalisa Pillepich, Dylan Nelson, Irina Zhuravleva, Wonki Lee, Mohammadreza Ayromlou, Katrin Lehle

Abstract: The intracluster medium (ICM) of galaxy clusters encodes the impact of the physical processes that shape these massive halos, including feedback from central supermassive black holes (SMBHs). In this study, we examine the gas thermodynamics, kinematics, and the effects of SMBH feedback on the core of Perseus-like galaxy clusters with a new simulation suite: TNG-Cluster. We first make a selection o… ▽ More The intracluster medium (ICM) of galaxy clusters encodes the impact of the physical processes that shape these massive halos, including feedback from central supermassive black holes (SMBHs). In this study, we examine the gas thermodynamics, kinematics, and the effects of SMBH feedback on the core of Perseus-like galaxy clusters with a new simulation suite: TNG-Cluster. We first make a selection of simulated clusters similar to Perseus based on the total mass and inner ICM properties, such as their cool-core nature. We identify 30 Perseus-like systems among the 352 TNG-Cluster halos at $z=0$. Many exhibit thermodynamical profiles and X-ray morphologies with disturbed features such as ripples, bubbles, and shock fronts that are qualitatively similar to X-ray observations of Perseus. To study observable gas motions, we generate XRISM mock X-ray observations and conduct a spectral analysis of the synthetic data. In agreement with existing Hitomi measurements, TNG-Cluster predicts subsonic gas turbulence in the central regions of Perseus-like clusters, with a typical line-of-sight velocity dispersion of 200 km/s. This implies that turbulent pressure contributes $< 10\%$ to the dominant thermal pressure. In TNG-Cluster, such low (inferred) values of ICM velocity dispersion coexist with high-velocity outflows and bulk motions of relatively small amounts of super-virial hot gas, moving up to thousands of km/s. However, detecting these outflows in observations may prove challenging due to their anisotropic nature and projection effects. Driven by SMBH feedback, such outflows are responsible for many morphological disturbances in the X-ray maps of cluster cores. They also increase both the inferred and intrinsic ICM velocity dispersion. This effect is somewhat stronger when velocity dispersion is measured from higher-energy lines. △ Less

Submitted 28 March, 2024; v1 submitted 10 November, 2023; originally announced November 2023.

Comments: 14 pages, 8 figures. Accepted to A&A. See companion papers (Ayromlou, Lee, Lehle, Nelson, Rohr) and additional information on the TNG-Cluster website www.tng-project.org/cluster/

Journal ref: A&A 686, A200 (2024)

The heart of galaxy clusters: demographics and physical properties of cool-core and non-cool-core halos in the TNG-Cluster simulation

Authors: Katrin Lehle, Dylan Nelson, Annalisa Pillepich, Nhut Truong, Eric Rohr

Abstract: We analyze the physical properties of the gaseous intracluster medium (ICM) at the center of massive galaxy clusters with TNG-Cluster, a new cosmological magnetohydrodynamical simulation. Our sample contains 352 simulated clusters spanning a halo mass range of $10^{14} < {\rm M}_{\rm 500c} / M_\odot < 2 \times 10^{15}$ at $z=0$. We focus on the proposed classification of clusters into cool-core (C… ▽ More We analyze the physical properties of the gaseous intracluster medium (ICM) at the center of massive galaxy clusters with TNG-Cluster, a new cosmological magnetohydrodynamical simulation. Our sample contains 352 simulated clusters spanning a halo mass range of $10^{14} < {\rm M}_{\rm 500c} / M_\odot < 2 \times 10^{15}$ at $z=0$. We focus on the proposed classification of clusters into cool-core (CC) and non-cool-core (NCC) populations, the $z=0$ distribution of cluster central ICM properties, and the redshift evolution of the CC cluster population. We analyze resolved structure and radial profiles of entropy, temperature, electron number density, and pressure. To distinguish between CC and NCC clusters, we consider several criteria: central cooling time, central entropy, central density, X-ray concentration parameter, and density profile slope. According to TNG-Cluster and with no a-priori cluster selection, the distributions of these properties are unimodal, whereby CCs and NCCs represent the two extremes. Across the entire TNG-Cluster sample at $z=0$ and based on central cooling time, the strong CC fraction is $f_{\rm SCC} = 24\%$, compared to $f_{\rm WCC} = 60\% $ and $f_{\rm NCC} = 16\%$ for weak and non-cool-cores, respectively. However, the fraction of CCs depends strongly on both halo mass and redshift, although the magnitude and even direction of the trends vary with definition. The abundant statistics of simulated high-mass clusters in TNG-Cluster enables us to match observational samples and make a comparison with data. The CC fractions from $z=0$ to $z=2$ are in broad agreement with observations, as are radial profiles of thermodynamical quantities, globally as well as split for CC versus NCC halos. TNG-Cluster can therefore be used as a laboratory to study the evolution and transformations of cluster cores due to mergers, AGN feedback, and other physical processes. △ Less

Submitted 10 November, 2023; originally announced November 2023.

Comments: Submitted to A&A, comments welcome. See the TNG-Cluster website at www.tng-project.org/cluster/ for more details

Journal ref: A&A 687, A129 (2024)

ERGO-ML: Comparing IllustrisTNG and HSC galaxy images via contrastive learning

Authors: Lukas Eisert, Connor Bottrell, Annalisa Pillepich, Rhythm Shimakawa, Vicente Rodriguez-Gomez, Dylan Nelson, Eirini Angeloudi, Marc Huertas-Company

Abstract: Modern cosmological hydrodynamical galaxy simulations provide tens of thousands of reasonably realistic synthetic galaxies across cosmic time. However, quantitatively assessing the level of realism of simulated universes in comparison to the real one is difficult. In this paper of the ERGO-ML series (Extracting Reality from Galaxy Observables with Machine Learning), we utilize contrastive learning… ▽ More Modern cosmological hydrodynamical galaxy simulations provide tens of thousands of reasonably realistic synthetic galaxies across cosmic time. However, quantitatively assessing the level of realism of simulated universes in comparison to the real one is difficult. In this paper of the ERGO-ML series (Extracting Reality from Galaxy Observables with Machine Learning), we utilize contrastive learning to directly compare a large sample of simulated and observed galaxies based on their stellar-light images. This eliminates the need to specify summary statistics and allows to exploit the whole information content of the observations. We produce survey-realistic galaxy mock datasets resembling real Hyper Suprime-Cam (HSC) observations using the cosmological simulations TNG50 and TNG100. Our focus is on galaxies with stellar masses between $10^9$ and $10^{12} M_\odot$ at $z=0.1-0.4$. This allows us to evaluate the realism of the simulated TNG galaxies in comparison to actual HSC observations. We apply the self-supervised contrastive learning method NNCLR to the images from both simulated and observed datasets (g, r, i - bands). This results in a 256-dimensional representation space, encoding all relevant observable galaxy properties. Firstly, this allows us to identify simulated galaxies that closely resemble real ones by seeking similar images in this multi-dimensional space. Even more powerful, we quantify the alignment between the representations of these two image sets, finding that the majority ($\gtrsim 70$ per cent) of the TNG galaxies align well with observed HSC images. However, a subset of simulated galaxies with larger sizes, steeper Sersic profiles, smaller Sersic ellipticities, and larger asymmetries appears unrealistic. We also demonstrate the utility of our derived image representations by inferring properties of real HSC galaxies using simulated TNG galaxies as the ground truth. △ Less

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

Comments: 26 pages, 17 figures. Comments are welcome. Highlights in Fig. 3, 7 and 13

Journal ref: MNRAS, Volume 528, Issue 4, March 2024, Pages 7411-7439

The physical origins of gas in the circumgalactic medium using observationally-motivated TNG50 mocks

Authors: Simon Weng, Celine Peroux, Rahul Ramesh, Dylan Nelson, Elaine M. Sadler, Martin Zwaan, Victoria Bollo, Benedetta Casavecchia

Abstract: Absorbers in the spectrum of background objects probe the circumgalactic medium (CGM) surrounding galaxies, but its physical properties remain unconstrained. We use the cosmological hydrodynamical simulation TNG50 to statistically trace the origins of HI Ly-$α$ absorbers around galaxies at $z = 0.5$ with stellar masses ranging from 10$^8$ to 10$^{11}$ M$_\odot$. We emulate observational CGM studie… ▽ More Absorbers in the spectrum of background objects probe the circumgalactic medium (CGM) surrounding galaxies, but its physical properties remain unconstrained. We use the cosmological hydrodynamical simulation TNG50 to statistically trace the origins of HI Ly-$α$ absorbers around galaxies at $z = 0.5$ with stellar masses ranging from 10$^8$ to 10$^{11}$ M$_\odot$. We emulate observational CGM studies by considering all gas within a line of sight velocity range of $\pm 500$ km s$^{-1}$ from the central, to quantitatively assess the impact of other galaxy haloes and overdense gas in the IGM that intersect sightlines. The impact of satellites to the total absorber fraction is most significant at impact parameters $0.5 R_{\rm vir} < b < R_{\rm vir}$ and satellites with masses below typical detection limits ($M_* < 10^8$ M$_\odot$) account for 10 (40) per cent of absorbers that intersect any satellite bound to $10^{10}$ and $10^{11}$ $(10^9)$ M$_\odot$ centrals. After confirming outflows are more dominant along the minor axis, we additionally show that at least 20 per cent of absorbers exhibit no significant radial movement, indicating that absorbers can also trace quasi-static gas. The metallicity of absorbers also depends on the azimuthal angle, but this signal is largely driven by enriched inflowing and quasi-static gas. Our work shows that determining the stellar mass of galaxies at $z_{\rm abs}$ is essential to constrain the physical origin of the gas traced in absorption, which in turn is key to characterising the kinematics and distribution of gas and metals in the CGM. △ Less

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

Comments: 23 pages, 13 figures. Accepted for publication in MNRAS

Io's polar volcanic thermal emission indicative of magma ocean and shallow tidal heating models

Authors: Ashley Gerard Davies, Jason Perry, David A. Williams, David M. Nelson

Abstract: The distribution of Io's volcanic activity likely reflects the position and magnitude of internal tidal heating. We use new observations of Io's polar regions by the Juno spacecraft Jovian Infrared Auroral Mapper (JIRAM) to complete near-infrared global coverage, revealing the global distribution and magnitude of thermal emission from Io's currently erupting volcanoes. We show that the distributio… ▽ More The distribution of Io's volcanic activity likely reflects the position and magnitude of internal tidal heating. We use new observations of Io's polar regions by the Juno spacecraft Jovian Infrared Auroral Mapper (JIRAM) to complete near-infrared global coverage, revealing the global distribution and magnitude of thermal emission from Io's currently erupting volcanoes. We show that the distribution of volcanic heat flow from 266 active hot spots is consistent with the presence of a global magma ocean, and/or shallow asthenospheric heating. We find that Io's polar volcanoes are less energetic but about the same in number per unit area than at lower latitudes. We also find that volcanic heat flow in the north polar cap is greater than that in the south. The low volcanic advection seen at Io's poles is therefore at odds with measurements of background temperature showing Io's poles are anomalously warm. We suggest that the differences in volcanic thermal emission from Io's poles compared to that at lower latitudes is indicative of lithospheric dichotomies that inhibit volcanic advection towards Io's poles, particularly in the south polar region. △ Less

Submitted 18 October, 2023; originally announced October 2023.

Comments: 17 pages, two tables, 7 figures

VERTICO and IllustrisTNG: The spatially resolved effects of environment on galactic gas

Authors: Adam R. H. Stevens, Toby Brown, Benedikt Diemer, Annalisa Pillepich, Lars Hernquist, Dylan Nelson, Yannick M. Bahé, Alessandro Boselli, Timothy A. Davis, Pascal J. Elahi, Sara L. Ellison, María J. Jiménez-Donaire, Ian D. Roberts, Kristine Spekkens, Vicente Villanueva, Adam B. Watts, Christine D. Wilson, Nikki Zabel

Abstract: It has been shown in previous publications that the TNG100 simulation quantitatively reproduces the observed reduction in each of the total atomic and total molecular hydrogen gas for galaxies within massive halos, i.e.~dense environments. In this Letter, we study how well TNG50 reproduces the resolved effects of a Virgo-like cluster environment on the gas surface densities of satellite galaxies w… ▽ More It has been shown in previous publications that the TNG100 simulation quantitatively reproduces the observed reduction in each of the total atomic and total molecular hydrogen gas for galaxies within massive halos, i.e.~dense environments. In this Letter, we study how well TNG50 reproduces the resolved effects of a Virgo-like cluster environment on the gas surface densities of satellite galaxies with $m_* > \! 10^9\,{\rm M}_\odot$ and ${\rm SFR} \! > 0.05\,{\rm M}_\odot\,{\rm yr}^{-1}$. We select galaxies in the simulation that are analogous to those in the HERACLES and VERTICO surveys, and mock-observe them to the common specifications of the data. Although TNG50 does not quantitatively match the observed gas surface densities in the centers of galaxies, the simulation does qualitatively reproduce the trends of gas truncation and central density suppression seen in VERTICO in both HI and H$_2$. This result promises that modern cosmological hydrodynamic simulations can be used to reliably model the post-infall histories of cluster satellite galaxies. △ Less

Submitted 11 October, 2023; originally announced October 2023.

Comments: 9 pages, 4 figures, accepted in ApJL

Periodic orbits, pair nucleation, and unbinding of active nematic defects on cones

Authors: Farzan Vafa, David R. Nelson, Amin Doostmohammadi

Abstract: Geometric confinement and topological constraints present promising means of controlling active materials. By combining analytical arguments derived from the Born-Oppenheimer approximation with numerical simulations, we investigate the simultaneous impact of confinement together with curvature singularity by characterizing the dynamics of an active nematic on a cone. Here, the Born-Oppenheimer app… ▽ More Geometric confinement and topological constraints present promising means of controlling active materials. By combining analytical arguments derived from the Born-Oppenheimer approximation with numerical simulations, we investigate the simultaneous impact of confinement together with curvature singularity by characterizing the dynamics of an active nematic on a cone. Here, the Born-Oppenheimer approximation means that textures can follow defect positions rapidly on the time scales of interest. Upon imposing strong anchoring boundary conditions at the base of a cone, we find a a rich phase diagram of multi-defect dynamics including exotic periodic orbits of one or two $+1/2$ flank defects, depending on activity and non-quantized geometric charge at the cone apex. By characterizing the transitions between these ordered dynamical states, we can understand (i) defect unbinding, (ii) defect absorption and (iii) defect pair nucleation at the apex. Numerical simulations confirm theoretical predictions of not only the nature of the circular orbits but also defect unbinding from the apex. △ Less

Submitted 9 October, 2023; originally announced October 2023.

Comments: 17 pages, 13 figures, 5 tables

Mapping the Intracluster Medium in the Era of High-resolution X-ray Spectroscopy

Authors: Congyao Zhang, Irina Zhuravleva, Maxim Markevitch, John ZuHone, François Mernier, Veronica Biffi, Ákos Bogdán, Priyanka Chakraborty, Eugene Churazov, Klaus Dolag, Stefano Ettori, William R. Forman, Christine Jones, Ildar Khabibullin, Caroline Kilbourne, Ralph Kraft, Erwin T. Lau, Sheng-Chieh Lin, Daisuke Nagai, Dylan Nelson, Anna Ogorzałek, Elena Rasia, Arnab Sarkar, Aurora Simionescu, Yuanyuan Su , et al. (2 additional authors not shown)

Abstract: High-resolution spectroscopy in soft X-rays will open a new window to map multiphase gas in galaxy clusters and probe physics of the intracluster medium (ICM), including chemical enrichment histories, circulation of matter and energy during large-scale structure evolution, stellar and black hole feedback, halo virialization, and gas mixing processes. An eV-level spectral resolution, large field-of… ▽ More High-resolution spectroscopy in soft X-rays will open a new window to map multiphase gas in galaxy clusters and probe physics of the intracluster medium (ICM), including chemical enrichment histories, circulation of matter and energy during large-scale structure evolution, stellar and black hole feedback, halo virialization, and gas mixing processes. An eV-level spectral resolution, large field-of-view, and effective area are essential to separate cluster emissions from the Galactic foreground and efficiently map the cluster outskirts. Several mission concepts that meet these criteria have been proposed recently, e.g., LEM, HUBS, and SuperDIOS. This theoretical study explores what information on ICM physics could be recovered with such missions and the associated challenges. We emphasize the need for a comprehensive comparison between simulations and observations to interpret the high-resolution spectroscopic observations correctly. Using Line Emission Mapper (LEM) characteristics as an example, we demonstrate that it enables the use of soft X-ray emission lines (e.g., O VII/VIII and Fe-L complex) from the cluster outskirts to measure the thermodynamic, chemical, and kinematic properties of the gas up to $r_{200}$ and beyond. By generating mock observations with full backgrounds, analysing their images/spectra with observational approaches, and comparing the recovered characteristics with true ones from simulations, we develop six key science drivers for future missions, including the exploration of multiphase gas in galaxy clusters (e.g., temperature fluctuations, phase-space distributions), metallicity, ICM gas bulk motions and turbulence power spectra, ICM-cosmic filament interactions, and advances for cluster cosmology. △ Less

Submitted 3 October, 2023; originally announced October 2023.

Comments: 24 pages, 26 figures, submitted to MNRAS. Comments are welcome

Geometric frustration of hard-disk packings on cones

Authors: Jessica H. Sun, Abigail Plummer, Grace H. Zhang, David R. Nelson, Vinothan N. Manoharan

Abstract: Conical surfaces pose an interesting challenge to crystal growth: a crystal growing on a cone can wrap around and meet itself at different radii. We use a disk-packing algorithm to investigate how this closure constraint can geometrically frustrate the growth of single crystals on cones with small opening angles. By varying the crystal seed orientation and cone angle, we find that -- except at spe… ▽ More Conical surfaces pose an interesting challenge to crystal growth: a crystal growing on a cone can wrap around and meet itself at different radii. We use a disk-packing algorithm to investigate how this closure constraint can geometrically frustrate the growth of single crystals on cones with small opening angles. By varying the crystal seed orientation and cone angle, we find that -- except at special commensurate cone angles -- crystals typically form a seam that runs along the axial direction of the cone, while near the tip, a disordered particle packing forms. We show that the onset of disorder results from a finite-size effect that depends strongly on the circumference and not on the seed orientation or cone angle. This finite-size effect occurs also on cylinders, and we present evidence that on both cylinders and cones, the defect density increases exponentially as circumference decreases. We introduce a simple model for particle attachment at the seam that explains the dependence on the circumference. Our findings suggest that the growth of single crystals can become frustrated even very far from the tip when the cone has a small opening angle. These results may provide insights into the observed geometry of conical crystals in biological and materials applications. △ Less

Submitted 25 September, 2023; originally announced September 2023.

Comments: 12 pages, 13 figures

Journal ref: Physical Review E 108, 054608 (2023)

Subconvex bounds for $U_{n+1}\times U_n$ in horizontal aspects

Authors: Yueke Hu, Paul D. Nelson

Abstract: For $L$-functions attached to automorphic representations of unitary groups $U_{n+1}\times U_n$, we establish a subconvex bound valid in certain horizontal aspects, where the set of ramified places is allowed to vary. For $L$-functions attached to automorphic representations of unitary groups $U_{n+1}\times U_n$, we establish a subconvex bound valid in certain horizontal aspects, where the set of ramified places is allowed to vary. △ Less

Submitted 15 December, 2023; v1 submitted 12 September, 2023; originally announced September 2023.

Comments: 45 pages

MSC Class: 11F67 (Primary) 11F70; 11M99 (Secondary)

Imposters among us: globular cluster kinematics and the halo mass of ultra-diffuse galaxies in clusters

Authors: Jessica E. Doppel, Laura V. Sales, José A. Benavides, Elisa Toloba, Eric W. Peng, Dylan Nelson, Julio F. Navarro

Abstract: The velocity dispersion of globular clusters (GCs) around ultra-diffuse galaxies (UDGs) in the Virgo cluster spans a wide range, including cases where GC kinematics suggest halos as massive as (or even more massive than) that of the Milky Way around these faint dwarfs. We analyze the catalogs of GCs derived in post-processing from the TNG50 cosmological simulation to study the GC system kinematics… ▽ More The velocity dispersion of globular clusters (GCs) around ultra-diffuse galaxies (UDGs) in the Virgo cluster spans a wide range, including cases where GC kinematics suggest halos as massive as (or even more massive than) that of the Milky Way around these faint dwarfs. We analyze the catalogs of GCs derived in post-processing from the TNG50 cosmological simulation to study the GC system kinematics and abundance of simulated UDGs in galaxy groups and clusters. UDGs in this simulation reside exclusively in dwarf-mass halos with $M_{200} \sim 10^{11}$ M$_{\odot}$. When considering only GCs gravitationally bound to simulated UDGs, we find GCs properties that overlap well with several observational measurements for UDGs. In particular, no bias towards overly massive halos is inferred from the study of bound GCs, confirming that GCs are good tracers of UDG halo mass. However, we find that contamination by intra-cluster GCs may, in some cases, substantially increase velocity dispersion estimates when performing projected mock observations of our sample. We caution that targets with less than $10$ GC tracers are particularly prone to severe uncertainties.Measuring the stellar kinematics of the host galaxy should help confirm the unusually massive halos suggested by GC kinematics around some UDGs △ Less

Submitted 6 September, 2023; originally announced September 2023.

Comments: 12 pages, 8 figures, submitted to MNRAS

IllustrisTNG in the HSC-SSP: image data release and the major role of mini mergers as drivers of asymmetry and star formation

Authors: Connor Bottrell, Hassen M. Yesuf, Gergö Popping, Kiyoaki Christopher Omori, Shenli Tang, Xuheng Ding, Annalisa Pillepich, Dylan Nelson, Lukas Eisert, Hua Gao, Andy D. Goulding, Boris S. Kalita, Wentao Luo, Jenny E. Greene, Jingjing Shi, John D. Silverman

Abstract: At fixed galaxy stellar mass, there is a clear observational connection between structural asymmetry and offset from the star forming main sequence, $Δ$SFMS. Herein, we use the TNG50 simulation to investigate the relative roles of major mergers (stellar mass ratios $μ\geq0.25$), minor ($0.1 \leq μ< 0.25$), and mini mergers ($0.01 \leq μ< 0.1$) in driving this connection amongst star forming galaxi… ▽ More At fixed galaxy stellar mass, there is a clear observational connection between structural asymmetry and offset from the star forming main sequence, $Δ$SFMS. Herein, we use the TNG50 simulation to investigate the relative roles of major mergers (stellar mass ratios $μ\geq0.25$), minor ($0.1 \leq μ< 0.25$), and mini mergers ($0.01 \leq μ< 0.1$) in driving this connection amongst star forming galaxies (SFGs). We use dust radiative transfer post-processing with SKIRT to make a large, public collection of synthetic Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) images of simulated TNG galaxies over $0.1\leq z \leq 0.7$ with $\log M_{\star} / \mathrm{M}_{\odot}\geq9$ ($\sim750$k images). Using their instantaneous SFRs, known merger histories/forecasts, and HSC-SSP asymmetries, we show (1) that TNG50 SFGs qualitatively reproduce the observed trend between $Δ$SFMS and asymmetry and (2) a strikingly similar trend emerges between $Δ$SFMS and the time-to-coalescence for mini mergers. Controlling for redshift, stellar mass, environment, and gas fraction, we show that individual mini merger events yield small enhancements in SFRs and asymmetries that are sustained on long timescales (at least $\sim3$ Gyr after coalescence, on average) -- in contrast to major/minor merger remnants which peak at much greater amplitudes but are consistent with controls only $\sim1$ Gyr after coalescence. Integrating the boosts in SFRs and asymmetries driven by $μ\geq0.01$ mergers since $z=0.7$ in TNG50 SFGs, we show that mini mergers are responsible for (i) $55$ per cent of all merger-driven star formation and (ii) $70$ per cent of merger-driven asymmetric structure. Due to their relative frequency and prolonged boost timescales, mini mergers dominate over their minor and major counterparts in driving star formation and asymmetry in SFGs. △ Less

Submitted 7 October, 2023; v1 submitted 28 August, 2023; originally announced August 2023.

Comments: 31 pages; 19 figures; MNRAS accepted; Image data available via the TNG website: www.tng-project.org/bottrell23

Zooming in on the circumgalactic medium: resolving small-scale gas structure with the GIBLE cosmological simulations

Authors: Rahul Ramesh, Dylan Nelson

Abstract: We introduce Project GIBLE (Gas Is Better resoLved around galaxiEs), a suite of cosmological zoom-in simulations where gas in the circumgalactic medium (CGM) is preferentially simulated at ultra-high numerical resolution. Our initial sample consists of eight galaxies, all selected as Milky Way-like galaxies at $z=0$ from the TNG50 simulation. Using the same galaxy formation model as IllustrisTNG,… ▽ More We introduce Project GIBLE (Gas Is Better resoLved around galaxiEs), a suite of cosmological zoom-in simulations where gas in the circumgalactic medium (CGM) is preferentially simulated at ultra-high numerical resolution. Our initial sample consists of eight galaxies, all selected as Milky Way-like galaxies at $z=0$ from the TNG50 simulation. Using the same galaxy formation model as IllustrisTNG, and the moving-mesh code AREPO, we re-simulate each of these eight galaxies maintaining a resolution equivalent to TNG50-2 ($m_{\rm{gas}}$ $\sim$ $8 \times 10^5 {\rm M}_{\odot}$). However, we use our super-Lagrangian refinement scheme to more finely resolve gas in the CGM around these galaxies. Our highest resolution runs achieve 512 times better mass resolution ($\sim$ $10^3 {\rm M}_{\odot}$). This corresponds to a median spatial resolution of $\sim$ $75$ pc at $0.15~R_{\rm{200,c}}$, which coarsens with increasing distance to $\sim$ $700$ pc at the virial radius. We make predictions for the covering fractions of several observational tracers of multi-phase CGM gas: HI, MgII, CIV and OVII. We then study the impact of improved resolution on small scale structure. While the abundance of the smallest cold, dense gas clouds continues to increase with improving resolution, the number of massive clouds is well converged. We conclude by quantifying small scale structure with the velocity structure function and the auto-correlation function of the density field, assessing their resolution dependence. The GIBLE cosmological hydrodynamical simulations enable us to improve resolution in a computationally efficient manner, thereby achieving numerical convergence of a subset of key CGM gas properties and observables. △ Less

Submitted 20 July, 2023; originally announced July 2023.

Comments: Submitted to MNRAS. Comments welcome

X-ray metal line emission from the hot circumgalactic medium: probing the effects of supermassive black hole feedback

Authors: Nhut Truong, Annalisa Pillepich, Dylan Nelson, Ákos Bogdán, Gerrit Schellenberger, Priyanka Chakraborty, William R. Forman, Ralph Kraft, Maxim Markevitch, Anna Ogorzalek, Benjamin D. Oppenheimer, Arnab Sarkar, Sylvain Veilleux, Mark Vogelsberger, Q. Daniel Wan, Norbert Werner, Irina Zhuravleva, John Zuhone

Abstract: We derive predictions from state-of-the-art cosmological galaxy simulations for the spatial distribution of the hot circumgalactic medium (CGM, ${\rm [0.1-1]R_{200c}}$) through its emission lines in the X-ray soft band ($[0.3-1.3]$ keV). In particular, we compare IllustrisTNG, EAGLE, and SIMBA and focus on galaxies with stellar mass $10^{10-11.6}\, \MSUN$ at $z=0$. The three simulation models retu… ▽ More We derive predictions from state-of-the-art cosmological galaxy simulations for the spatial distribution of the hot circumgalactic medium (CGM, ${\rm [0.1-1]R_{200c}}$) through its emission lines in the X-ray soft band ($[0.3-1.3]$ keV). In particular, we compare IllustrisTNG, EAGLE, and SIMBA and focus on galaxies with stellar mass $10^{10-11.6}\, \MSUN$ at $z=0$. The three simulation models return significantly different surface brightness radial profiles of prominent emission lines from ionized metals such as OVII(f), OVIII, and FeXVII as a function of galaxy mass. Likewise, the three simulations predict varying azimuthal distributions of line emission with respect to the galactic stellar planes, with IllustrisTNG predicting the strongest angular modulation of CGM physical properties at radial range ${\gtrsim0.3-0.5\,R_{200c}}$. This anisotropic signal is more prominent for higher-energy lines, where it can manifest as X-ray eROSITA-like bubbles. Despite different models of stellar and supermassive black hole (SMBH) feedback, the three simulations consistently predict a dichotomy between star-forming and quiescent galaxies at the Milky-Way and Andromeda mass range, where the former are X-ray brighter than the latter. This is a signature of SMBH-driven outflows, which are responsible for quenching star formation. Finally, we explore the prospect of testing these predictions with a microcalorimeter-based X-ray mission concept with a large field-of-view. Such a mission would probe the extended hot CGM via soft X-ray line emission, determine the physical properties of the CGM, including temperature, from the measurement of line ratios, and provide critical constraints on the efficiency and impact of SMBH feedback on the CGM. △ Less

Submitted 26 August, 2023; v1 submitted 3 July, 2023; originally announced July 2023.

Comments: 22 pages, 16 figures. Accepted for publication in MNRAS

Properties of the Line-of-Sight Velocity Field in the Hot and X-ray Emitting Circumgalactic Medium of Nearby Simulated Disk Galaxies

Authors: J. A. ZuHone, G. Schellenberger, A. Ogorzalek, B. D. Oppenheimer, J. Stern, A. Bogdan, N. Truong, M. Markevitch, A. Pillepich, D. Nelson, J. N. Burchett, I. Khabibullin, C. A. Kilbourne, R. P. Kraft, P. E. J. Nulsen, S. Veilleux, M. Vogelsberger, Q. D. Wang, I. Zhuravleva

Abstract: The hot, X-ray-emitting phase of the circumgalactic medium of massive galaxies is believed to be the reservoir of baryons from which gas flows onto the central galaxy and into which feedback from AGN and stars inject mass, momentum, energy, and metals. These effects shape the velocity fields of the hot gas, which can be observed via the Doppler shifting and broadening of emission lines by X-ray IF… ▽ More The hot, X-ray-emitting phase of the circumgalactic medium of massive galaxies is believed to be the reservoir of baryons from which gas flows onto the central galaxy and into which feedback from AGN and stars inject mass, momentum, energy, and metals. These effects shape the velocity fields of the hot gas, which can be observed via the Doppler shifting and broadening of emission lines by X-ray IFUs. In this work, we analyze the gas kinematics of the hot circumgalactic medium of Milky Way-mass disk galaxies from the TNG50 simulation with synthetic observations to determine how future instruments can probe this velocity structure. We find that the hot phase is often characterized by outflows from the disk driven by feedback processes, radial inflows near the galactic plane, and rotation, though in some systems the velocity field is more disorganized and turbulent. With a spectral resolution of $\sim$1 eV, fast and hot outflows ($\sim$200-500 km s$^{-1}$) can be measured, depending on the orientation of the galaxy on the sky. The rotation velocity of the hot phase ($\sim$100-200 km s$^{-1}$) can be measured using line shifts in edge-on galaxies, and is slower than that of colder gas phases but similar to stellar rotation velocities. By contrast, the slow inflows ($\sim$50-100 km s$^{-1}$) are difficult to measure in projection with these other components, but may be detected in multi-component spectral fits. We find that the velocity measured is sensitive to which emission lines are used. Measuring these flows will constrain theories of how the gas in these galaxies evolves. △ Less

Submitted 20 May, 2024; v1 submitted 3 July, 2023; originally announced July 2023.

Comments: 46 pages, 32 figures, accepted to ApJ

Mapping the imprints of stellar and AGN feedback in the circumgalactic medium with X-ray microcalorimeters

Authors: Gerrit Schellenberger, Ákos Bogdán, John A. ZuHone, Benjamin D. Oppenheimer, Nhut Truong, Ildar Khabibullin, Fred Jennings, Annalisa Pillepich, Joseph Burchett, Christopher Carr, Priyanka Chakraborty, Robert Crain, William Forman, Christine Jones, Caroline A. Kilbourne, Ralph P. Kraft, Maxim Markevitch, Daisuke Nagai, Dylan Nelson, Anna Ogorzalek, Scott Randall, Arnab Sarkar, Joop Schaye, Sylvain Veilleux, Mark Vogelsberger , et al. (2 additional authors not shown)

Abstract: The Astro2020 Decadal Survey has identified the mapping of the circumgalactic medium (CGM, gaseous plasma around galaxies) as a key objective. We explore the prospects for characterizing the CGM in and around nearby galaxy halos with a future, large grasp X-ray microcalorimeter. We create realistic mock observations from hydrodynamical simulations (EAGLE, IllustrisTNG, and Simba) that demonstrate… ▽ More The Astro2020 Decadal Survey has identified the mapping of the circumgalactic medium (CGM, gaseous plasma around galaxies) as a key objective. We explore the prospects for characterizing the CGM in and around nearby galaxy halos with a future, large grasp X-ray microcalorimeter. We create realistic mock observations from hydrodynamical simulations (EAGLE, IllustrisTNG, and Simba) that demonstrate a wide range of potential measurements, which will address the open questions in galaxy formation and evolution. By including all background and foreground components in our mock observations, we show why it is impossible to perform these measurements with current instruments, such as X-ray CCDs, and only microcalorimeters will allow us to distinguish the faint CGM emission from the bright Milky Way (MW) foreground emission lines. We find that individual halos of MW mass can, on average and depending on star formation rate, be traced out to large radii, around R500, and for larger galaxies even out to R200, using prominent emission lines, such as OVII, or OVIII. Furthermore, we show that emission line ratios for individual halos can reveal the radial temperature structure. Substructure measurements show that it will be possible to relate azimuthal variations to the feedback mode of the galaxy. We demonstrate the ability to construct temperature, velocity, and abundance ratio maps from spectral fitting for individual galaxy halos, which reveal rotation features, AGN outbursts, and enrichment. △ Less

Submitted 29 April, 2024; v1 submitted 3 July, 2023; originally announced July 2023.

Comments: 41 pages, 19 figures, accepted for publication in ApJ

Formation of Galactic Disks II: the Physical Drivers of Disk Spin-up

Authors: Vadim A. Semenov, Charlie Conroy, Vedant Chandra, Lars Hernquist, Dylan Nelson

Abstract: Using a representative sample of Milky Way (MW)-like galaxies from the TNG50 cosmological-volume simulation, we investigate physical processes driving the formation of galactic disks. A disk forms as a result of the interplay between inflow and outflow carrying angular momentum in and out of the galaxy. Interestingly, the inflow and outflow have remarkably similar distributions of angular momentum… ▽ More Using a representative sample of Milky Way (MW)-like galaxies from the TNG50 cosmological-volume simulation, we investigate physical processes driving the formation of galactic disks. A disk forms as a result of the interplay between inflow and outflow carrying angular momentum in and out of the galaxy. Interestingly, the inflow and outflow have remarkably similar distributions of angular momentum, suggesting an exchange of angular momentum and/or outflow recycling, leading to continuous feeding of pre-aligned material from the co-rotating circumgalactic medium. We show that disk formation in TNG50 is correlated with stellar bulge formation, in qualitative agreement with a recent theoretical model of disk formation facilitated by steep gravitational potentials. Disk formation is also correlated with the formation of a hot circumgalactic halo with a significant fraction of the inflow occurring at sub- and transonic velocities. In the context of recent theoretical works connecting disk settling and hot halo formation, our results imply that the subsonic part of the inflow may settle into a disk while the remaining supersonic inflow will perturb this disk via the chaotic cold accretion. We find that disks tend to form when the host halos become more massive than $\sim (1-2) \times 10^{11} M_\odot$, consistent with previous theoretical findings and observational estimates of the pre-disk protogalaxy remnant in the MW. Our results do not prove that either co-rotating outflow recycling, gravitational potential steepening, or hot halo formation cause disk formation but they show that all these processes occur concurrently and may play an important role in disk growth. △ Less

Submitted 22 June, 2023; originally announced June 2023.

Comments: 22 pages, 15 figures; submitted to ApJ; comments are welcome

Formation of Galactic Disks I: Why Did the Milky Way's Disk Form Unusually Early?

Authors: Vadim A. Semenov, Charlie Conroy, Vedant Chandra, Lars Hernquist, Dylan Nelson

Abstract: Recent results from spectroscopic and astrometric surveys of nearby stars suggest that the stellar disk of our Milky Way (MW) was formed quite early, within the first few billion years of its evolution. Chemokinematic signatures of disk formation in cosmological zoom-in simulations appear to be in tension with these data, implying that MW-like disk formation is delayed in simulations. We investiga… ▽ More Recent results from spectroscopic and astrometric surveys of nearby stars suggest that the stellar disk of our Milky Way (MW) was formed quite early, within the first few billion years of its evolution. Chemokinematic signatures of disk formation in cosmological zoom-in simulations appear to be in tension with these data, implying that MW-like disk formation is delayed in simulations. We investigate the formation of galactic disks using a representative sample of MW-like galaxies from the cosmological-volume simulation TNG50. We find that on average MW-mass disks indeed form later than the local data suggest. However, their formation time and metallicity exhibit a substantial scatter, such that $\sim$10% of MW-mass galaxies form disks early, similar to the MW. Thus, although the MW is unusual, it is consistent with the overall population of MW-mass disk galaxies. The direct MW analogs assemble most of their mass early, $\gtrsim 10$ Gyr ago, and are not affected by destructive mergers after that. In addition, these galaxies form their disks during the early enrichment stage when the interstellar medium metallicity increases rapidly, with only $\sim$25% of early-forming disks being as metal-poor as the MW was at the onset of disk formation, [Fe/H] $\approx -1.0$. In contrast, most MW-mass galaxies either form disks from already enriched material or experience late destructive mergers that reset the signatures of galactic disk formation to later times and higher metallicities. Finally, we also show that earlier disk formation leads to more dominant rotationally supported stellar disks at redshift zero. △ Less

Submitted 22 January, 2024; v1 submitted 15 June, 2023; originally announced June 2023.

Comments: 19 pages, 13 figures + appendix; accepted for publication in ApJ

Resonant scattering of the OVII X-ray emission line in the circumgalactic medium of TNG50 galaxies

Authors: Dylan Nelson, Chris Byrohl, Anna Ogorzalek, Maxim Markevitch, Ildar Khabibullin, Eugene Churazov, Irina Zhuravleva, Akos Bogdan, Priyanka Chakraborty, Caroline Kilbourne, Ralph Kraft, Annalisa Pillepich, Arnab Sarkar, Gerrit Schellenberger, Yuanyuan Su, Nhut Truong, Stephan Vladutescu-Zopp, Nastasha Wijers

Abstract: We study the impact of resonantly scattered X-ray line emission on the observability of the hot circumgalactic medium (CGM) of galaxies. We apply a Monte Carlo radiative transfer post-processing analysis to the high-resolution TNG50 cosmological magnetohydrodynamical galaxy formation simulation. This allows us to model the resonant scattering of OVII(r) X-ray photons within the complex, multi-phas… ▽ More We study the impact of resonantly scattered X-ray line emission on the observability of the hot circumgalactic medium (CGM) of galaxies. We apply a Monte Carlo radiative transfer post-processing analysis to the high-resolution TNG50 cosmological magnetohydrodynamical galaxy formation simulation. This allows us to model the resonant scattering of OVII(r) X-ray photons within the complex, multi-phase, multi-scale CGM. The resonant transition of the OVII He-like triplet is one of the brightest, and most promising, X-ray emission lines for detecting the hot CGM and measuring its physical properties. We focus on galaxies with stellar masses 10 < log(M*/Msun) < 11 at z ~ 0. After constructing a model for OVII(r) emission from the central galaxy as well as from CGM gas, we forward model these intrinsic photons to derive observable surface brightness maps. We find that scattering significantly boosts the observable OVII(r) surface brightness of the extended and diffuse CGM. This enhancement can be large -- an order of magnitude on average at a distance of 200 projected kpc for high-mass M* = 10^10.7 Msun galaxies. The enhancement is larger for lower mass galaxies, and can even reach a factor of 100, across the extended CGM. Galaxies with higher star formation rates, AGN luminosities, and central OVII(r) luminosities all have larger scattering enhancements, at fixed stellar mass. Our results suggest that next-generation X-ray spectroscopic missions including XRISM, LEM, ATHENA, and HUBS -- which aim to detect the hot CGM in emission -- could specifically target halos with significant enhancements due to resonant scattering. △ Less

Submitted 8 June, 2023; originally announced June 2023.

Comments: Published in MNRAS. See https://www.lem-observatory.org/ and https://www.tng-project.org/ for more details; 2023MNRAS.522.3665N

Circumgalactic Medium on the Largest Scales: Detecting X-ray Absorption Lines with Large-Area Microcalorimeters

Authors: Akos Bogdan, Ildar Khabibullin, Orsolya Kovacs, Gerrit Schellenberger, John ZuHone, Joseph Burchett, Klaus Dolag, Eugene Churazov, William Forman, Christine Jones, Caroline Kilbourne, Ralph Kraft, Erwin Lau, Maxim Markevitch, Dan McCammon, Daisuke Nagai, Dylan Nelson, Anna Ogorzalek, Benjamin Oppenheimer, Arnab Sarkar, Yuanyuan Su, Nhut Truong, Sylvain Veilleux, Stephan Vladutescu-Zopp, Irina Zhuravleva

Abstract: The circumgalactic medium (CGM) plays a crucial role in galaxy evolution as it fuels star formation, retains metals ejected from the galaxies, and hosts gas flows in and out of galaxies. For Milky Way-type and more massive galaxies, the bulk of the CGM is in hot phases best accessible at X-ray wavelengths. However, our understanding of the CGM remains largely unconstrained due to its tenuous natur… ▽ More The circumgalactic medium (CGM) plays a crucial role in galaxy evolution as it fuels star formation, retains metals ejected from the galaxies, and hosts gas flows in and out of galaxies. For Milky Way-type and more massive galaxies, the bulk of the CGM is in hot phases best accessible at X-ray wavelengths. However, our understanding of the CGM remains largely unconstrained due to its tenuous nature. A promising way to probe the CGM is via X-ray absorption studies. Traditional absorption studies utilize bright background quasars, but this method probes the CGM in a pencil beam, and, due to the rarity of bright quasars, the galaxy population available for study is limited. Large-area, high spectral resolution X-ray microcalorimeters offer a new approach to exploring the CGM in emission and absorption. Here, we demonstrate that the cumulative X-ray emission from cosmic X-ray background sources can probe the CGM in absorption. We construct column density maps of major X-ray ions from the Magneticum simulation and build realistic mock images of nine galaxies to explore the detectability of X-ray absorption lines arising from the large-scale CGM. We conclude that the OVII absorption line is detectable around individual massive galaxies at the $3σ-6σ$ confidence level. For Milky Way-type galaxies, the OVII and OVIII absorption lines are detectable at the $\sim\,6σ$ and $\sim\,3σ$ levels even beyond the virial radius when co-adding data from multiple galaxies. This approach complements emission studies, does not require additional exposures, and will allow probing of the baryon budget and the CGM at the largest scales. △ Less

Submitted 8 June, 2023; originally announced June 2023.

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

Remembrances of Michael E. Fisher

Authors: David R. Nelson

Abstract: This contribution will be published in '50 years of the renormalization group', dedicated to the memory of Michael E. Fisher, edited by Amnon Aharony, Ora Entin-Wohlman, David Huse, and Leo Radzihovsky, World Scientific. These personal remembrances come in three parts. The first contains a brief personal perspective on Michael E. Fisher's contributions to science. The second tells how I came to wo… ▽ More This contribution will be published in '50 years of the renormalization group', dedicated to the memory of Michael E. Fisher, edited by Amnon Aharony, Ora Entin-Wohlman, David Huse, and Leo Radzihovsky, World Scientific. These personal remembrances come in three parts. The first contains a brief personal perspective on Michael E. Fisher's contributions to science. The second tells how I came to work with Michael and describes events while I was under his supervision during my graduate years at Cornell University. The third part summarizes recent work, done in collaboration with Suraj Shankar, on thermalized buckling of isotopically compressed thin (perhaps atomically thin) sheets of materials such as graphene or MoS$_2$. These investigations were inspired by Michael's beautiful work on the effect of constraints at critical points, with fluctuations at all length scales, which leads to 'Fisher renormalization' of critical exponents. [1] Thin fluctuating sheets embedded in three dimensions, when they are tensionless as in a cantilever or 'diving board' geometry, are automatically at a critical point everywhere in a low temperature flat phase. However, when we consider thin sheets supported on multiple sides in various ways, Fisher's ideas lead to the inequivalence of isotensional and isometric thermodynamic ensembles, which triggers dramatic differences in some of the critical exponents associated with the two types of boundary conditions. [2] Readers not interested in my experiences while a student at Cornell University may only want to read parts I and III. There are also some concluding remarks. △ Less

Submitted 1 June, 2023; originally announced June 2023.

Comments: 31 pages, 5 figures, contribution to '50 years of the renormalization group', book honoring Michael E. Fisher

The BarYon CYCLE Project (ByCycle): Identifying and Localizing MgII Metal Absorbers with Machine Learning

Authors: Roland Szakacs, Céline Péroux, Dylan Nelson, Martin A. Zwaan, Daniel Grün, Simon Weng, Alejandra Y. Fresco, Victoria Bollo, Benedetta Casavecchia

Abstract: The upcoming ByCycle project on the VISTA/4MOST multi-object spectrograph will offer new prospects of using a massive sample of $\sim 1$ million high spectral resolution ($R$ = 20,000) background quasars to map the circumgalactic metal content of foreground galaxies (observed at $R$ = 4000 - 7000), as traced by metal absorption. Such large surveys require specialized analysis methodologies. In the… ▽ More The upcoming ByCycle project on the VISTA/4MOST multi-object spectrograph will offer new prospects of using a massive sample of $\sim 1$ million high spectral resolution ($R$ = 20,000) background quasars to map the circumgalactic metal content of foreground galaxies (observed at $R$ = 4000 - 7000), as traced by metal absorption. Such large surveys require specialized analysis methodologies. In the absence of early data, we instead produce synthetic 4MOST high-resolution fibre quasar spectra. To do so, we use the TNG50 cosmological magnetohydrodynamical simulation, combining photo-ionization post-processing and ray tracing, to capture MgII ($\lambda2796$, $\lambda2803$) absorbers. We then use this sample to train a Convolutional Neural Network (CNN) which searches for, and estimates the redshift of, MgII absorbers within these spectra. For a test sample of quasar spectra with uniformly distributed properties ($λ_{\rm{MgII,2796}}$, $\rm{EW}_{\rm{MgII,2796}}^{\rm{rest}} = 0.05 - 5.15$ Å, $\rm{SNR} = 3 - 50$), the algorithm has a robust classification accuracy of 98.6 per cent and a mean wavelength accuracy of 6.9 Å. For high signal-to-noise spectra ($\rm{SNR > 20}$), the algorithm robustly detects and localizes MgII absorbers down to equivalent widths of $\rm{EW}_{\rm{MgII,2796}}^{\rm{rest}} = 0.05$ Å. For the lowest SNR spectra ($\rm{SNR=3}$), the CNN reliably recovers and localizes EW$_{\rm{MgII,2796}}^{\rm{rest}}$ $\geq$ 0.75 Å\, absorbers. This is more than sufficient for subsequent Voigt profile fitting to characterize the detected MgII absorbers. We make the code publicly available through GitHub. Our work provides a proof-of-concept for future analyses of quasar spectra datasets numbering in the millions, soon to be delivered by the next generation of surveys. △ Less

Submitted 29 May, 2023; originally announced May 2023.

Comments: 13 pages, 9 figures, 1 table. Accepted for publication in MNRAS

A Virgo Environmental Survey Tracing Ionised Gas Emission (VESTIGE) XV. The Halpha luminosity function of the Virgo cluster

Authors: A. Boselli, M. Fossati, P. Cote, J. C. Cuillandre, L. Ferrarese, S. Gwyn, P. Amram, M. Ayromlou, M. Balogh, G. Bellusci, M. Boquien, G. Gavazzi, G. Hensler, A. Longobardi, D. Nelson, A. Pillepich, J. Roediger, R. Sanchez-Jansen, M. Sun, G. Trinchieri

Abstract: We use a complete set of deep narrow-band imaging data for 384 galaxies gathered during the VESTIGE survey to derive the first Halpha luminosity function (LF) of the Virgo cluster within R200. The data allow us to cover the whole dynamic range of the Halpha LF (10^36<LHa<10^42 erg s^-1). After they are corrected for [NII] contamination and dust attenuation, the data are used to derive the SFR func… ▽ More We use a complete set of deep narrow-band imaging data for 384 galaxies gathered during the VESTIGE survey to derive the first Halpha luminosity function (LF) of the Virgo cluster within R200. The data allow us to cover the whole dynamic range of the Halpha LF (10^36<LHa<10^42 erg s^-1). After they are corrected for [NII] contamination and dust attenuation, the data are used to derive the SFR function in the range 10^-4<SFR<10 Mo yr^-1. These LF are compared to those derived at other frequencies or using different tracers of star formation in Virgo, in other nearby and high-z clusters, in the field, and to those predicted by the IllustrisTNG cosmological hydrodynamical simulations. The Halpha LF of the Virgo cluster is fairly flat (a=-1.07) in the range 10^38.5<LHa<10^40.5 erg s^-1, and it abruptly decreases at lower luminosities. When compared to those derived for other nearby clusters and for the field, the slope and the characteristic luminosity of the Schechter function change as a function of the dynamical mass of the system, of the temperature of the X-rays gas, and of the dynamical pressure exerted on the interstellar medium of galaxies moving at high velocity within the intracluster medium. All these trends can be explained in a scenario in which the activity of SF is reduced in massive clusters due to their hydrodynamical interaction with the surrounding medium, suggesting once again that ram-pressure stripping is the dominant mechanism affecting galaxy evolution in local clusters of dynamical mass M200>10^14 Mo. The comparison with the IllustrisTNG cosmological hydrodynamical simulations shows a more pronounced decrease at the faint end of the distribution. If Virgo is representative of typical nearby clusters of similar mass, this difference suggests that the stripping process in simulated galaxies in these environments is more efficient than observed. △ Less

Submitted 25 May, 2023; originally announced May 2023.

Comments: Accepted for publication on A&A

Journal ref: A&A 675, A123 (2023)

MUSE-ALMA Halos XI: Gas flows in the circumgalactic medium

Authors: Simon Weng, Céline Péroux, Arjun Karki, Ramona Augustin, Varsha P. Kulkarni, Aleksandra Hamanowicz, Martin Zwaan, Elaine M. Sadler, Dylan Nelson, Matthew J. Hayes, Glenn G. Kacprzak, Andrew J. Fox, Victoria Bollo, Benedetta Casavecchia, Roland Szakacs

Abstract: The flow of gas into and out of galaxies leaves traces in the circumgalactic medium which can then be studied using absorption lines towards background quasars. We analyse 27 log(N_HI) > 18.0 HI absorbers at z = 0.2 to 1.4 from the MUSE-ALMA Halos survey with at least one galaxy counterpart within a line of sight velocity of +/-500 km s^{-1}. We perform 3D kinematic forward modelling of these asso… ▽ More The flow of gas into and out of galaxies leaves traces in the circumgalactic medium which can then be studied using absorption lines towards background quasars. We analyse 27 log(N_HI) > 18.0 HI absorbers at z = 0.2 to 1.4 from the MUSE-ALMA Halos survey with at least one galaxy counterpart within a line of sight velocity of +/-500 km s^{-1}. We perform 3D kinematic forward modelling of these associated galaxies to examine the flow of dense, neutral gas in the circumgalactic medium. From the VLT/MUSE, HST broadband imaging and VLT/UVES and Keck/HIRES high-resolution UV quasar spectroscopy observations, we compare the impact parameters, star-formation rates and stellar masses of the associated galaxies with the absorber properties. We find marginal evidence for a bimodal distribution in azimuthal angles for strong HI absorbers, similar to previous studies of the MgII and OVI absorption lines. There is no clear metallicity dependence on azimuthal angle and we suggest a larger sample of absorbers are required to fully test the relationship predicted by cosmological hydrodynamical simulations. A case-by-case study of the absorbers reveals that ten per cent of absorbers are consistent with gas accretion, up to 30 per cent trace outflows while the remainder trace gas in the galaxy disk, the intragroup medium and low-mass galaxies below the MUSE detection limit. Our results highlight that the baryon cycle directly affects the dense neutral gas required for star-formation and plays a critical role in galaxy evolution. △ Less

Submitted 18 May, 2023; originally announced May 2023.

Comments: 13 pages, 6 figures, 12 pages of appendix. Accepted for publication in MNRAS

Azimuthal Anisotropy of Magnetic Fields in the Circumgalactic Medium Driven by Galactic Feedback Processes

Authors: Rahul Ramesh, Dylan Nelson, Volker Heesen, Marcus Brüggen

Abstract: We use the TNG50 cosmological magnetohydrodynamical simulation of the IllustrisTNG project to show that magnetic fields in the circumgalactic medium (CGM) have significant angular structure. This azimuthal anisotropy at fixed distance is driven by galactic feedback processes that launch strong outflows into the halo, preferentially along the minor axes of galaxies. These feedback-driven outflows e… ▽ More We use the TNG50 cosmological magnetohydrodynamical simulation of the IllustrisTNG project to show that magnetic fields in the circumgalactic medium (CGM) have significant angular structure. This azimuthal anisotropy at fixed distance is driven by galactic feedback processes that launch strong outflows into the halo, preferentially along the minor axes of galaxies. These feedback-driven outflows entrain strong magnetic fields from the interstellar medium, dragging fields originally amplified by small-scale dynamos into the CGM. At the virial radius, $z=0$ galaxies with M$_\star \sim 10^{10}\,\rm{M_\odot}$ show the strongest anisotropy ($\sim 0.35$ dex). This signal weakens with decreasing impact parameter, and is also present but weaker for lower mass as well as higher mass galaxies. Creating mock Faraday rotation measure (RM) sightlines through the simulated volume, we find that the angular RM trend is qualitatively consistent with recent observational measurements. We show that rich structure is present in the circumgalactic magnetic fields of galaxies. However, TNG50 predicts small RM amplitudes in the CGM that make detection difficult as a result of other contributions along the line of sight. △ Less

Submitted 10 October, 2023; v1 submitted 18 May, 2023; originally announced May 2023.

Comments: Accepted for publication (MNRAS)

The signature of galaxy formation models in the power spectrum of the hydrogen 21cm line during reionization

Authors: Joseph S. W. Lewis, Annalisa Pillepich, Dylan Nelson, Ralf S. Klessen, Simon C. O. Glover

Abstract: Observations of the 21cm line of neutral hydrogen are poised to revolutionize our knowledge of cosmic reionization and the high-redshift population of galaxies. However, harnessing such information requires robust and comprehensive theoretical modeling. We study the non-linear effects of hydrodynamics and astrophysical feedback processes, including stellar and AGN feedback, on the 21cm signal by p… ▽ More Observations of the 21cm line of neutral hydrogen are poised to revolutionize our knowledge of cosmic reionization and the high-redshift population of galaxies. However, harnessing such information requires robust and comprehensive theoretical modeling. We study the non-linear effects of hydrodynamics and astrophysical feedback processes, including stellar and AGN feedback, on the 21cm signal by post-processing three existing cosmological hydrodynamical simulations of galaxy formation: Illustris, IllustrisTNG, and Eagle. Overall and despite their different underlying galaxy-formation models, the three simulations return similar predictions for the global 21cm rightness temperature and its power spectrum. At fixed redshift, most differences are attributable to differences in the history of reionization, in turn driven by differences in the build-up of stellar sources of radiation. However, the impact of astrophysics is imprinted in the 21cm power spectrum through several unique signatures. First, we find significant small scale ($k \geq 10\, \rm {Mpc}^{-1}$) differences between Illustris and IllustrisTNG, where higher velocity winds generated by supernova feedback soften density peaks and lead to lower 21cm power in TNG. Second, we find more 21cm power at intermediate scales ($k \approx 0.8\, \rm {Mpc}^{-1}$) in Eagle, due to differences in ionization driven by highly effective stellar feedback, leading to lower star formation, older and redder stellar populations, and thus lower ionizing luminosities. Though subtle, these features could allow future observations of the 21cm signal, in conjunction with other reionization observables, to constrain theoretical models for galactic feedback at high redshift. △ Less

Submitted 18 December, 2023; v1 submitted 16 May, 2023; originally announced May 2023.

Comments: submitted to MNRAS (12/05/23); Accepted 07/12/23 after minor revision

Vibrations and tunneling of strained nanoribbons at finite temperature

Authors: Paul Z. Hanakata, Sourav S. Bhabesh, David Yllanes, David R. Nelson, Mark J. Bowick

Abstract: Crystalline sheets (e.g., graphene and transition metal dichalcogenides) liberated from a substrate are a paradigm for materials at criticality because flexural phonons can fluctuate into the third dimension. Although studies of static critical behaviors (e.g., the scale-dependent elastic constants) are plentiful, investigations of dynamics remain limited. Here, we use molecular dynamics to study… ▽ More Crystalline sheets (e.g., graphene and transition metal dichalcogenides) liberated from a substrate are a paradigm for materials at criticality because flexural phonons can fluctuate into the third dimension. Although studies of static critical behaviors (e.g., the scale-dependent elastic constants) are plentiful, investigations of dynamics remain limited. Here, we use molecular dynamics to study the time dependence of the midpoint (the height center-of-mass) of doubly clamped nanoribbons, as prototypical graphene resonators, under a wide range of temperature and strain conditions. By treating the ribbon midpoint as a Brownian particle confined to a nonlinear potential (which assumes a double-well shape beyond the buckling transition), we formulate an effective theory describing the ribbon's tunneling rate across the two wells and its oscillations inside a given well. We find that, for nanoribbbons compressed above the Euler buckling point and thermalized above a temperature at which the non-linear effects due to thermal fluctuations become significant, the exponential term (the ratio between energy barrier and temperature) depends only on the geometry, but not the temperature, unlike the usual Arrhenius behavior. Moreover, we find that the natural oscillation time for small strain shows a non-trivial scaling $τ_{\rm o}\sim L_0^{\,z}T^{-η/4}$, with $L_0$ being the ribbon length, $z=2-η/2$ being the dynamic critical exponent, $η=0.8$ being the scaling exponent describing scale-dependent elastic constants, and $T$ being the temperature. These unusual scale- and temperature-dependent dynamics thus exhibit dynamic criticality and could be exploited in the development of graphene-based nanoactuators. △ Less

Submitted 11 May, 2023; originally announced May 2023.

Journal ref: Phys. Rev. Materials 8, 016001 (2024)

Jellyfish galaxies with the IllustrisTNG simulations -- Citizen-science results towards large distances, low-mass hosts, and high redshifts

Authors: Elad Zinger, Gandhali Joshi, Annalisa Pillepich, Eric Rohr, Dylan Nelson

Abstract: We present the ``Cosmological Jellyfish'' project - a citizen-science classification program to identify jellyfish galaxies within the IllustrisTNG cosmological simulations. Jellyfish (JF) are satellite galaxies that exhibit long trailing gas features -- `tails' -- extending from their stellar body. Their distinctive morphology arises due to ram-pressure stripping (RPS) as they move through the ba… ▽ More We present the ``Cosmological Jellyfish'' project - a citizen-science classification program to identify jellyfish galaxies within the IllustrisTNG cosmological simulations. Jellyfish (JF) are satellite galaxies that exhibit long trailing gas features -- `tails' -- extending from their stellar body. Their distinctive morphology arises due to ram-pressure stripping (RPS) as they move through the background gaseous medium. Using the TNG50 and TNG100 simulations, we construct a sample of $\sim 80,000$ satellite galaxies spanning an unprecedented range of stellar masses, $10^{8.3-12.3}\,\mathrm{M_\odot}$, and host masses of $M_\mathrm{200,c}=10^{10.4-14.6}\,\mathrm{M_\odot}$ back to $z=2$ \citep[extending the work of][]{yun_jellyfish_2019}. Based on this sample, $\sim 90,000$ galaxy images were presented to volunteers in a citizen-science project on the Zooniverse platform who were asked to determine if each galaxy image resembles a jellyfish. Based on volunteer votes, each galaxy was assigned a score determining if it is a JF or not. This paper describes the project, the inspected satellite sample, the methodology, and the classification process that resulted in a dataset of $5,307$ visually-identified jellyfish galaxies. We find that JF galaxies are common in nearly all group- and cluster-sized systems, with the JF fraction increasing with host mass and decreasing with satellite stellar mass. We highlight JF galaxies in three relatively unexplored regimes: low-mass hosts of $M_\mathrm{200,c}\sim10^{11.5-13}\,\mathrm{M_\odot}$, radial positions within hosts exceeding the virial radius $R_\mathrm{200,c}$, and at high redshift up to $z=2$. The full dataset of our jellyfish scores is publicly available and can be used to select and study JF galaxies in the IllustrisTNG simulations. △ Less

Submitted 18 April, 2023; originally announced April 2023.

Comments: submitted to MNRAS ; See additional jellyfish companion papers today on astro-ph: Rohr et al. and Goeller et al.; Jellyfish image gallery: https://www.tng-project.org/explore/gallery/zinger23/

Jellyfish galaxies with the IllustrisTNG simulations -- When, where, and for how long does ram pressure stripping of cold gas occur?

Authors: Eric Rohr, Annalisa Pillepich, Dylan Nelson, Elad Zinger, Gandhali Joshi, Mohommadreza Ayromlou

Abstract: Jellyfish galaxies are prototypical examples of satellite galaxies undergoing strong ram pressure stripping (RPS). We analyze the evolution of 512 unique, first-infalling jellyfish galaxies from the TNG50 cosmological simulation. These have been visually inspected to be undergoing RPS sometime in the past 5 billion years (since $z=0.5$), have satellite stellar masses… ▽ More Jellyfish galaxies are prototypical examples of satellite galaxies undergoing strong ram pressure stripping (RPS). We analyze the evolution of 512 unique, first-infalling jellyfish galaxies from the TNG50 cosmological simulation. These have been visually inspected to be undergoing RPS sometime in the past 5 billion years (since $z=0.5$), have satellite stellar masses $\mstarsat\sim10^{8-10.5}\,M_\odot$, and live in hosts with $\mvir\sim10^{12-14.3}\,M_\odot$ at $z=0$. We quantify the cold gas ($T\leq10^{4.5}$ K) removal using the tracer particles, confirming that for these jellyfish, RPS is the dominant driver of cold gas loss after infall. Half of these jellyfish are completely gas-less by $z=0$, and these galaxies have earlier infall times and smaller satellite-to-host mass ratios than their gaseous counterparts. RPS can act on jellyfish galaxies over long time scales of $\approx1.5-8$ Gyr. Jellyfish in more massive hosts are impacted by RPS for a shorter time span and, at a fixed host mass, jellyfish with less cold gas at infall and lower stellar masses at $z=0$ have shorter RPS time spans. While RPS may act for long periods of time, the peak RPS period -- where at least 50 pecent of the total RPS occurs -- begins within $\approx1$ Gyr of infall and lasts $\lesssim2$ Gyr. During this period, the jellyfish are at host-centric distances $\sim0.2-2\rvir$, illustrating that much of RPS occurs at large distances from the host galaxy. Interestingly, jellyfish continue forming stars until they have lost $\approx98$ percent of their cold gas. For groups and clusters in TNG50 $(\mvirhost\sim10^{13-14.3}\,M_\odot)$, jellyfish galaxies deposit more cold gas ($\sim10^{11-12}\,M_\odot$) into halos than exist in them at $z=0$, demonstrating that jellyfish, and in general satellite galaxies, are a significant source of cold gas accretion. △ Less

Submitted 16 April, 2024; v1 submitted 18 April, 2023; originally announced April 2023.

Comments: 20 pages, 11 figures + 3 appendices with 4 figures. Submitted to MNRAS. Key figures are 2, 8, 9, 11. See additional jellyfish companion papers today on astro-ph: Zinger+ and Goeller+. All data used in this publication, including the Cosmological Jellyfish Project results, are publicaly available. Reuploaded after publication

Journal ref: MNRAS, 524, 3502-3525 (2023)

Active topological defect absorption by a curvature singularity

Authors: Farzan Vafa, David R. Nelson, Amin Doostmohammadi

Abstract: Using the Born-Oppenheimer approximation, we present a general description of topological defects dynamics in $p$-atic materials on curved surfaces, and simplify it in the case of active nematics. We find that activity induces a geometric contribution to the motility of the $+1/2$ defect. Moreover, in the case of a cone, the simplest example of a geometry with curvature singularity, we find that t… ▽ More Using the Born-Oppenheimer approximation, we present a general description of topological defects dynamics in $p$-atic materials on curved surfaces, and simplify it in the case of active nematics. We find that activity induces a geometric contribution to the motility of the $+1/2$ defect. Moreover, in the case of a cone, the simplest example of a geometry with curvature singularity, we find that the motility depends on the deficit angle of the cone and changes sign when the deficit angle is bigger than $π$, leading to the change in active behavior from contractile (extensile) to extensile (contractile) behavior. Using our analytical framework, we then identify for positively charged defects the basin of attraction to the cone apex and present closed-form predictions for defect trajectories near the apex. The analytical results are quantitatively corroborated against full numerical simulations. Provided the capture radius is small compared to the cone size, the agreement is excellent. △ Less

Submitted 13 April, 2023; originally announced April 2023.

Comments: 14 + 4 pages, 7 figures

Disk flaring with TNG50: diversity across Milky Way and M31 analogs

Authors: Diego Sotillo-Ramos, Martina Donnari, Annalisa Pillepich, Neige Frankel, Dylan Nelson, Volker Springel, Lars Hernquist

Abstract: We use the sample of 198 Milky Way (MW) and Andromeda (M31) analogs from TNG50 to quantify the level of disk flaring predicted by a modern, high-resolution cosmological hydrodynamical simulation. Disk flaring refers to the increase of vertical stellar disk height with galactocentric distance. The TNG50 galaxies are selected to have stellar disky morphology, a stellar mass in the range of… ▽ More We use the sample of 198 Milky Way (MW) and Andromeda (M31) analogs from TNG50 to quantify the level of disk flaring predicted by a modern, high-resolution cosmological hydrodynamical simulation. Disk flaring refers to the increase of vertical stellar disk height with galactocentric distance. The TNG50 galaxies are selected to have stellar disky morphology, a stellar mass in the range of $M_* = 10^{10.5 - 11.2}~\rm{M_{\odot}}$, and a MW-like Mpc-scale environment at $z=0$. The stellar disks of such TNG50 MW/M31 analogs exhibit a wide diversity of structural properties, including a number of galaxies with disk scalelength and thin and thick disk scaleheights that are comparable to those measured or inferred for the Galaxy and Andromeda. With one set of physical ingredients, TNG50 returns a large variety of flaring flavours and amounts, also for mono-age stellar populations. With this paper, we hence propose a non-parametric characterization of flaring. The typical MW/M31 analogs exhibit disk scaleheights that are $1.5-2$ times larger in the outer than in the inner regions of the disk for both old and young stellar populations, but with a large galaxy-to-galaxy variation. Which stellar population flares more, and by how much, also varies from galaxy to galaxy. TNG50 de facto brackets existing observational constraints for the Galaxy and all previous numerical findings. A link between the amount of flaring and the $z=0$ global galaxy structural properties or merger history is complex. However, a connection between the scaleheights and the local stellar vertical kinematics and gravitational potential is clearly in place. △ Less

Submitted 28 March, 2023; originally announced March 2023.

Comments: Submitted to MNRAS. Main figures: 5, 13. See presentation and data release of TNG50 MW/M31 analogs by Pillepich et al. and see also Ramesh et al. on astro-ph today

Milky Way and Andromeda analogs from the TNG50 simulation

Authors: Annalisa Pillepich, Diego Sotillo-Ramos, Rahul Ramesh, Dylan Nelson, Christoph Engler, Vicente Rodriguez-Gomez, Martin Fournier, Martina Donnari, Volker Springel, Lars Hernquist

Abstract: We present the properties of Milky Way- and Andromeda-like (MW/M31-like) galaxies simulated within TNG50, the highest-resolution run of the IllustrisTNG suite of $Λ$CDM magneto-hydrodynamical simulations. We introduce our fiducial selection for MW/M31 analogs, which we propose for direct usage as well as for reference in future analyses. TNG50 contains 198 MW/M31 analogs, i.e. galaxies with stella… ▽ More We present the properties of Milky Way- and Andromeda-like (MW/M31-like) galaxies simulated within TNG50, the highest-resolution run of the IllustrisTNG suite of $Λ$CDM magneto-hydrodynamical simulations. We introduce our fiducial selection for MW/M31 analogs, which we propose for direct usage as well as for reference in future analyses. TNG50 contains 198 MW/M31 analogs, i.e. galaxies with stellar disky morphology, with a stellar mass in the range of $M_* = 10^{10.5 - 11.2}$ Msun, and within a MW-like Mpc-scale environment at z=0. These are resolved with baryonic (dark matter) mass resolution of $8.5\times10^4$ Msun ($4.5\times10^5$ Msun) and $\sim150$ pc of average spatial resolution in the star-forming regions: we therefore expand by many factors (2 orders of magnitude) the sample size of cosmologically-simulated analogs with similar ($\times 10$ better) numerical resolution. The majority of TNG50 MW/M31 analogs at $z=0$ exhibit a bar, 60 per cent are star-forming, the sample includes 3 Local Group (LG)-like systems, and a number of galaxies host one or more satellites as massive as e.g. the Magellanic Clouds. Even within such a relatively narrow selection, TNG50 reveals a great diversity in galaxy and halo properties, as well as in past histories. Within the TNG50 sample, it is possible to identify several simulated galaxies whose integral and structural properties are consistent, one or more at a time, with those measured for the Galaxy and Andromeda. With this paper, we document and release a series of broadly applicable data products that build upon the IllustrisTNG public release and aim to facilitate easy access and analysis by public users. These include datacubes across snapshots ($0 \le z \le 7$) for each TNG50 MW/M31-like galaxy, and a series of value-added catalogs that will be continually expanded to provide a convenient and up to date community resource. △ Less

Submitted 28 March, 2023; originally announced March 2023.

Comments: Submitted to MNRAS. Feedback from the community encouraged. This is also a data-release paper: see visuals and documentation at https://www.tng-project.org/data/milkyway+andromeda/. See other TNG50 MW/M31-based papers also on astro-ph today: Ramesh+ and Sotillo-Ramos+

The Circumgalactic Medium of Milky Way-like Galaxies in the TNG50 Simulation -- II: Cold, Dense Gas Clouds and High-Velocity Cloud Analogs

Authors: Rahul Ramesh, Dylan Nelson, Annalisa Pillepich

Abstract: We use the TNG50 simulation of the IllustrisTNG project to study cold, dense clouds of gas in the circumgalactic media (CGM) of Milky Way-like galaxies. We find that their CGM is typically filled with of order one hundred (thousand) reasonably (marginally) resolved clouds, possible analogs of high-velocity clouds (HVCs). There is a large variation in cloud abundance from galaxy to galaxy, and the… ▽ More We use the TNG50 simulation of the IllustrisTNG project to study cold, dense clouds of gas in the circumgalactic media (CGM) of Milky Way-like galaxies. We find that their CGM is typically filled with of order one hundred (thousand) reasonably (marginally) resolved clouds, possible analogs of high-velocity clouds (HVCs). There is a large variation in cloud abundance from galaxy to galaxy, and the physical properties of clouds that we explore -- mass, size, metallicity, pressure, and kinematics -- are also diverse. We quantify the distributions of cloud properties and cloud-background contrasts, providing cosmological inputs for idealized simulations. Clouds characteristically have sub-solar metallicities, diverse shapes, small overdensities ($χ= n_{\rm cold} / n_{\rm hot} \lesssim 10$), are mostly inflowing, and have sub-virial rotation. At TNG50 resolution, resolved clouds have median masses of $\sim 10^6\,\rm{M_\odot}$ and sizes of $\sim 10$ kpc. Larger clouds are well converged numerically, while the abundance of the smallest clouds increases with resolution, as expected. In TNG50 MW-like haloes, clouds are slightly (severely) under-pressurised relative to their surroundings with respect to total (thermal) pressure, implying that magnetic fields may be important. Clouds are not distributed uniformly throughout the CGM, but are clustered around other clouds, often near baryon-rich satellite galaxies. This suggests that at least some clouds originate from satellites, via direct ram-pressure stripping or otherwise. Finally, we compare with observations of intermediate and high velocity clouds from the real Milky Way halo. TNG50 shows a similar cloud velocity distribution as observations, and predicts a significant population of currently difficult-to-detect low velocity clouds. △ Less

Submitted 28 March, 2023; originally announced March 2023.

Comments: Accepted for publication (MNRAS). Part of a set of papers based on TNG50 MW/M31-like galaxies. Additional visuals and data products at www.tng-project.org/ramesh23b

The atomic-to-molecular hydrogen transition in the TNG50 simulation: Using realistic UV fields to create spatially resolved HI maps

Authors: Andrea Gebek, Maarten Baes, Benedikt Diemer, W. J. G. de Blok, Dylan Nelson, Anand Utsav Kapoor, Peter Camps, Omphile Rabyang, Lerothodi Leeuw

Abstract: Cold gas in galaxies provides a crucial test to evaluate the realism of cosmological hydrodynamical simulations. To extract the atomic and molecular hydrogen properties of the simulated galaxy population, postprocessing methods taking the local UV field into account are required. We improve upon previous studies by calculating realistic UV fields with the dust radiative transfer code SKIRT to mode… ▽ More Cold gas in galaxies provides a crucial test to evaluate the realism of cosmological hydrodynamical simulations. To extract the atomic and molecular hydrogen properties of the simulated galaxy population, postprocessing methods taking the local UV field into account are required. We improve upon previous studies by calculating realistic UV fields with the dust radiative transfer code SKIRT to model the atomic-to-molecular transition in TNG50, the highest-resolution run of the IllustrisTNG suite. Comparing integrated quantities such as the HI mass function, we study to what detail the UV field needs to be modelled in order to calculate realistic cold gas properties. We then evaluate new, spatially resolved comparisons for cold gas in galaxies by exploring synthetic maps of atomic hydrogen at redshift zero and compare them to 21-cm observations of local galaxies from the WHISP survey. In terms of non-parametric morphologies, we find that TNG50 HI maps are less concentrated than their WHISP counterparts (median $ΔC\approx0.3$), due in part to central HI deficits related to the ejective character of supermassive black hole feedback in TNG. In terms of the HI column density distribution function, we find discrepancies between WHISP and IllustrisTNG that depend on the total HI abundance in these datasets as well as the postprocessing method. To fully exploit the synergy between cosmological simulations and upcoming deep HI/H2 data, we advocate the use of accurate methods to estimate the UV radiation field and to generate mock maps. △ Less

Submitted 13 March, 2023; originally announced March 2023.

Comments: 21 pages, 10 figures (main text). Accepted to MNRAS