Xianjin Shen, Hong-Li Liu, Zhiyuan Ren, Anandmayee Tej, Di Li, Hauyu Baobab Liu, Gary A. Fuller, Jinjin Xie, Sihan Jiao, Aiyuan Yang, Patrick M. Koch, Fengwei Xu, Patricio Sanhueza, Pham N. Diep, Nicolas Peretto, Ram K. Yadav, Busaba H. Kramer, Koichiro Sugiyama, Mark Rawlings, Chang Won Lee, et al (15) Filaments are believed to play a key role in high-mass star formation. We present a systematic study of the filaments and their hosting clumps in the G35 molecular complex using JCMT SCUBA-2 850 $\micron$ continuum data. We identified five clouds in the complex and 91 filaments within them, some of which form 10 hub-filament systems (HFSs), each with at least 3 hub-composing filaments. We also compiled a catalogue of 350 dense clumps, 183 of which are associated with the filaments. We investigated the physical properties of the filaments and clumps, such as mass, density, and size, and their relation to star formation. We find that the global mass-length trend of the filaments is consistent with a turbulent origin, while the hub-composing filaments of high line masses ($m_{\rm l}\,>$\u2009230\,$\mathrm{M_{\odot}~pc^{-1}}$) in HFSs deviate from this relation, possibly due to feedback from massive star formation. We also find that the most massive and densest clumps (R\,$>$\u20090.2\u2009pc, M\,$>35\,\mathrm{M_{\odot}}$, $\mathrm{\Sigma}>\,0.05\,\mathrm{g~cm^{-2}}$) are located in the filaments and in the hubs of HFS with the latter bearing a higher probability of occurrence of high-mass star-forming signatures, highlighting the preferential sites of HFSs for high-mass star formation. We do not find significant variation in the clump mass surface density across different evolutionary environments of the clouds, which may reflect the balance between mass accretion and stellar feedback.
Tri Nguyen, Francisco Villaescusa-Navarro, Siddharth Mishra-Sharma, Carolina Cuesta-Lazaro, Paul Torrey, Arya Farahi, Alex M. Garcia, Jonah C. Rose, Stephanie O'Neil, Mark Vogelsberger, Xuejian Shen, Cian Roche, Daniel Anglés-Alcázar, Nitya Kallivayalil, Julian B. Muñoz, Francis-Yan Cyr-Racine, Sandip Roy, Lina Necib, Kassidy E. Kollmann The connection between galaxies and their host dark matter (DM) halos is critical to our understanding of cosmology, galaxy formation, and DM physics. To maximize the return of upcoming cosmological surveys, we need an accurate way to model this complex relationship. Many techniques have been developed to model this connection, from Halo Occupation Distribution (HOD) to empirical and semi-analytic models to hydrodynamic. Hydrodynamic simulations can incorporate more detailed astrophysical processes but are computationally expensive; HODs, on the other hand, are computationally cheap but have limited accuracy. In this work, we present NeHOD, a generative framework based on variational diffusion model and Transformer, for painting galaxies/subhalos on top of DM with an accuracy of hydrodynamic simulations but at a computational cost similar to HOD. By modeling galaxies/subhalos as point clouds, instead of binning or voxelization, we can resolve small spatial scales down to the resolution of the simulations. For each halo, NeHOD predicts the positions, velocities, masses, and concentrations of its central and satellite galaxies. We train NeHOD on the TNG-Warm DM suite of the DREAMS project, which consists of 1024 high-resolution zoom-in hydrodynamic simulations of Milky Way-mass halos with varying warm DM mass and astrophysical parameters. We show that our model captures the complex relationships between subhalo properties as a function of the simulation parameters, including the mass functions, stellar-halo mass relations, concentration-mass relations, and spatial clustering. Our method can be used for a large variety of downstream applications, from galaxy clustering to strong lensing studies.
We present the first suite of cosmological hydrodynamical zoom-in simulations of isolated dwarf galaxies for a dark sector that consists of Cold Dark Matter and a strongly-dissipative sub-component. The simulations are implemented in GIZMO and include standard baryons following the FIRE-2 galaxy formation physics model. The dissipative dark matter is modeled as Atomic Dark Matter (aDM), which forms a dark hydrogen gas that cools in direct analogy to the Standard Model. Our suite includes seven different simulations of $\sim 10^{10} M_{\odot}$ systems that vary over the aDM microphysics and the dwarf's evolutionary history. We identify a region of aDM parameter space where the cooling rate is aggressive and the resulting halo density profile is universal. In this regime, the aDM gas cools rapidly at high redshifts and only a small fraction survives in the form of a central dark gas disk; the majority collapses centrally into collisionless dark "clumps", which are clusters of sub-resolution dark compact objects. These dark clumps rapidly equilibrate in the inner galaxy, resulting in an approximately isothermal distribution that can be modeled with a simple fitting function. Even when only a small fraction ($\sim 5\%$) of the total dark matter is strongly dissipative, the central densities of classical dwarf galaxies can be enhanced by over an order of magnitude, providing a sharp prediction for observations.
Desika Narayanan, Daniel P. Stark, Steven L. Finkelstein, Paul Torrey, Qi Li, Fergus Cullen, Micheal W. Topping, Federico Marinacci, Laura V. Sales, Xuejian Shen, Mark Vogelsberger JWST has enabled the detection of the UV continuum of galaxies at z>10, evidencing a population of extremely blue, potentially dust-free galaxies. Interpreting the UV spectra of galaxies as they redden is complicated by the well-known degeneracy between stellar ages, dust, and nebular continuum. The main goal of this paper is to develop a theoretical model for the relationship between galaxy UV slopes, bursty star formation histories, dust evolution, and the contribution from nebular regions. We accomplish this via cosmological zoom-in simulations, and in specific, build a layered model where we simulate the UV slopes of galaxies with increasingly complex physics. Our main results follow. (i) Unattenuated stellar populations with no nebular emission exhibit a diverse range of intrinsic UV slopes, with values ranging from beta ~ -3 --> -2.2 due to long delays between bursts. This is manifested by an inverse correlation between the intrinsic UV slope and sSFR for early galaxies such that higher sSFR corresponds to bluer UV slopes. (ii) When including dust, our model galaxies demonstrate a rapid rise in dust obscuration between z ~ 8-10. This increase in dust mass is due to high grain-grain shattering rates, and enhanced growth per unit dust mass in very small grains, resulting in UV-detected galaxies at z ~ 12 descending into ALMA-detectable galaxies by z ~ 6. The rapid rise in dust content at z ~ 8-10 leads to a systematic reddening of the UV slopes during this redshift range. (iii) The inclusion of nebular continuum reddens the UV slope by a median factor Delta beta ~ 0.2-0.4. However, when including nebular continuum, our highest redshift galaxies (z~12) are insufficiently blue compared to observations; this may imply an evolving escape fraction from HII regions with redshift.
Robert Feldmann, Michael Boylan-Kolchin, James S. Bullock, Onur Çatmabacak, Claude-André Faucher-Giguère, Christopher C. Hayward, Dušan Kereš, Alexandres Lazar, Lichen Liang, Jorge Moreno, Pascal A. Oesch, Eliot Quataert, Xuejian Shen, Guochao Sun Recent observations with the James Webb Space Telescope (JWST) have uncovered unexpectedly high cosmic star formation activity in the early Universe, mere hundreds of millions of years after the Big Bang. These observations are often understood to reflect an evolutionary shift in star formation efficiency (SFE) caused by changing galactic conditions during these early epochs. We present FIREbox-HR, a high-resolution, cosmological hydrodynamical simulation from the Feedback in Realistic Environments project, which offers insights into the SFE of galaxies during the first billion years of cosmic time. FIREbox-HR re-simulates the cosmic volume (L = 22.1 cMpc) of the original FIREbox run with eight times higher mass resolution (m_b ~ 7800 M_sun), but with identical physics, down to z ~ 6. FIREbox-HR predicts ultraviolet (UV) luminosity functions in good agreement with available observational data. The simulation also successfully reproduces the observed cosmic UV luminosity density at z ~ 6 - 14, demonstrating that relatively high star formation activity in the early Universe is a natural outcome of the baryonic processes encoded in the FIRE-2 model. According to FIREbox-HR, the SFE - halo mass relation for intermediate mass halos (M_halo ~ 10^9 - 10^11 M_sun) does not significantly evolve with redshift and is only weakly mass-dependent. These properties of the SFE - halo mass relation lead to a larger contribution from lower mass halos at higher z, driving the gradual evolution of the observed cosmic UV luminosity density. A theoretical model based on the SFE - halo mass relation inferred from FIREbox-HR allows us to explore implications for galaxy evolution. Future observations of UV faint galaxies at z > 12 will provide an opportunity to further test these predictions and deepen our understanding of star formation during Cosmic Dawn.
JWST has revealed a large population of ultra-violet (UV)-bright galaxies at $z\gtrsim 10$ and possibly overly massive galaxies at $z\gtrsim 7$, challenging standard galaxy formation models in the $\Lambda$CDM cosmology. We use an empirical galaxy formation model to explore the potential of alleviating these tensions through an Early Dark Energy (EDE) model, originally proposed to solve the Hubble tension. Our benchmark model demonstrates excellent agreement with the UV luminosity functions (UVLFs) at $4\lesssim z \lesssim10$ in both $\Lambda$CDM and EDE cosmologies. In the EDE cosmology, the UVLF measurements at $z\simeq 12$ based on spectroscopically confirmed galaxies exhibit no tension with the benchmark model. Photometric constraints at $12 \lesssim z\lesssim 16$ can be fully explained within EDE via either moderately increased star formation efficiencies ($\epsilon_{\ast}\sim 3-10\%$ at $M_{\rm halo}\sim 10^{10.5}\,{\rm M}_\odot$) or enhanced UV variabilities ($\sigma_{\rm UV}\sim 0.8-1.3$ mag at $M_{\rm halo}\sim 10^{10.5}\,{\rm M}_\odot$) that are within the scatter of hydrodynamical simulation predictions. A similar agreement is difficult to achieve in $\Lambda$CDM, especially at $z\gtrsim 14$, where the required $\sigma_{\rm UV}$ exceeds the maximum value seen in simulations. Furthermore, the implausibly large cosmic stellar mass densities inferred from some JWST observations are no longer in tension with cosmology when the EDE is considered. Our findings highlight EDE as an intriguing unified solution to a fundamental problem in cosmology and the recent tensions raised by JWST observations. Data at the highest redshifts reached by JWST ($z \sim 14-16$) will be crucial for differentiating modified galaxy formation physics from new cosmological physics.
We investigate a self-interacting dark matter (SIDM) model featuring a velocity-dependent cross section with an order-of-magnitude resonant enhancement of the cross section at $\sim 16\,{\rm km}\,{\rm s}^{-1}$. To understand the implications for the structure of dark matter halos, we perform N-body simulations of isolated dark matter halos of mass $\sim 10^8\,{\rm M}_\odot$, a halo mass selected to have a maximum response to the resonance. We track the core formation and the gravothermal collapse phases of the dark matter halo in this model and compare the halo evolving with the resonant cross section with halos evolving with velocity-independent cross sections. We show that dark matter halo evolution with the resonant cross section exhibits a deviation from universality that characterizes halo evolution with velocity-independent cross sections. The halo evolving under the influence of the resonance reaches a lower minimum central density during core formation. It subsequently takes about $20\%$ longer to reach its initial central density during the collapse phase. These results motivate a more detailed exploration of halo evolution in models with pronounced resonances.
Isabel S. Sands, Philip F. Hopkins, Xuejian Shen, Michael Boylan-Kolchin, James Bullock, Claude-Andre Faucher-Giguere, Francisco J. Mercado, Jorge Moreno, Lina Necib, Xiaowei Ou, Sarah Wellons, Andrew Wetzel While galaxy rotation curves provide one of the most powerful methods for measuring dark matter profiles in the inner regions of rotation-supported galaxies, at the dwarf scale there are factors that can complicate this analysis. Given the expectation of a universal profile in dark matter-only simulations, the diversity of observed rotation curves has become an often-discussed issue in Lambda Cold Dark Matter cosmology on galactic scales. We analyze a suite of Feedback in Realistic Environments (FIRE) simulations of $10^{10}-10^{12}$ $M_\odot$ halos with standard cold dark matter, and compare the true circular velocity to rotation curve reconstructions. We find that, for galaxies with well-ordered gaseous disks, the measured rotation curve may deviate from true circular velocity by at most 10% within the radius of the disk. However, non-equilibrium behavior, non-circular motions, and non-thermal and non-kinetic stresses may cause much larger discrepancies of 50% or more. Most rotation curve reconstructions underestimate the true circular velocity, while some reconstructions transiently over-estimate it in the central few kiloparsecs due to dynamical phenomena. We further demonstrate that the features that contribute to these failures are not always visibly obvious in HI observations. If such dwarf galaxies are included in galaxy catalogs, they may give rise to the appearance of "artificial" rotation curve diversity that does not reflect the true variation in underlying dark matter profiles.
Xuejian Shen, Mark Vogelsberger, Josh Borrow, Yongao Hu, Evan Erickson, Rahul Kannan, Aaron Smith, Enrico Garaldi, Lars Hernquist, Takahiro Morishita, Sandro Tacchella, Oliver Zier, Guochao Sun, Anna-Christina Eilers, Hui Wang We investigate galaxy sizes at redshift $z\gtrsim 6$ with the cosmological radiation-magneto-hydrodynamic simulation suite THESAN(-HR). These simulations simultaneously capture the reionization of the large-scale intergalactic medium and resolved galaxy properties. The intrinsic size ($r^{\ast}_{1/2}$) of simulated galaxies increases moderately with stellar mass at $M_{\ast} \lesssim 10^{8}\,{\rm M}_{\odot}$ and decreases fast at larger masses, resulting in a hump feature at $M_{\ast}\sim 10^{8}\,{\rm M}_{\odot}$ that is insensitive to redshift. Low-mass galaxies are in the initial phase of size growth and are better described by a spherical shell model with feedback-driven gas outflows competing with the cold inflows. In contrast, massive galaxies fit better with the disk formation model. They generally experience a phase of rapid compaction and gas depletion, likely driven by internal disk instability rather than external processes. We identify four compact quenched galaxies in the $(95.5\,{\rm cMpc})^{3}$ volume of THESAN-1 at $z\simeq 6$, and their quenching follows reaching a characteristic stellar surface density akin to the massive compact galaxies at cosmic noon. Compared to observations, we find that the median UV effective radius ($R^{\rm UV}_{\rm eff}$) of simulated galaxies is at least three times larger than the observed ones at $M_{\ast}\lesssim 10^{9}\,{\rm M}_{\odot}$ or $M_{\rm UV}\gtrsim -20$ at $6 \lesssim z \lesssim 10$. This inconsistency, related to the hump feature of the intrinsic size--mass relation, persists across many other cosmological simulations with different galaxy formation models and demonstrates the potential of using galaxy morphology to constrain the physics of galaxy formation at high redshifts.
The distribution of offsets between the brightest cluster galaxies of galaxy clusters and the centroid of their dark matter distributions is a promising probe of the underlying dark matter physics. In particular, since this distribution is sensitive to the shape of the potential in galaxy cluster cores, it constitutes a test of dark matter self-interaction on the largest mass scales in the universe. We examine these offsets in three suites of modern cosmological simulations; IllustrisTNG, MillenniumTNG and BAHAMAS. For clusters above $10^{14}\rm{M_\odot}$, we examine the dependence of the offset distribution on gravitational softening length, the method used to identify centroids, redshift, mass, baryonic physics, and establish the stability of our results with respect to various nuisance parameter choices. We find that offsets are overwhelmingly measured to be smaller than the minimum converged length scale in each simulation, with a median offset of $\sim1\rm{kpc}$ in the highest resolution simulation considered, TNG300-1, which uses a gravitational softening length of $1.48\rm{kpc}$. We also find that centroids identified via source extraction on smoothed dark matter and stellar particle data are consistent with the potential minimum, but that observationally relevant methods sensitive to cluster strong gravitational lensing scales, or those using gas as a tracer for the potential can overestimate offsets by factors of $\sim10$ and $\sim30$, respectively. This has the potential to reduce tensions with existing offset measurements which have served as evidence for a nonzero dark matter self-interaction cross section.
Self-interacting dark matter (SIDM) is increasingly studied as a potential solution to small-scale discrepancies between simulations of cold dark matter (CDM) and observations. We examine a physically motivated two-state SIDM model with both elastic and inelastic scatterings. In particular, endothermic, exothermic, and elastic scattering occur with equal probability at high relative velocities ($v_{\rm rel}\gtrsim400~{\rm km/s})$. In a suite of cosmological zoom-in simulation of Milky Way-size haloes, we vary the primordial state fractions to understand the impact of inelastic dark matter self-interactions on halo structure and evolution. In particular, we test how the initial conditions impact the present-day properties of dark matter haloes. Depending on the primordial state fraction, scattering reactions will be dominated by either exothermic or endothermic effects for high and low initial excited state fractions respectively. We find that increasing the initial excited fraction reduces the mass of the main halo, as well as the number of subhaloes on all mass scales. The main haloes are cored, with lower inner densities and higher outer densities compared with CDM. Additionally, we find that the shape of the main halo becomes more spherical the higher the initial excited state fraction is. Finally, we show that the number of satellites steadily decreases with initial excited state fraction across all satellite masses.
Using cosmological hydrodynamical zoom-in simulations, we explore the properties of subhalos in Milky Way analogs that contain a sub-component of Atomic Dark Matter (ADM). ADM differs from Cold Dark Matter (CDM) due to the presence of self interactions that lead to energy dissipation and bound-state formation, analogous to Standard Model baryons. This model can arise in complex dark sectors that are natural and theoretically-motivated extensions to the Standard Model. The simulations used in this work were carried out using GIZMO and utilize the FIRE-2 galaxy formation physics in the Standard Model baryonic sector. For the parameter points we consider, the ADM gas cools efficiently, allowing it to collapse to the center of subhalos. This increases a subhalo's central density and affects its orbit, with more subhalos surviving small pericentric passages. The subset of subhalos that host visible satellite galaxies have cuspier density profiles and smaller stellar-half-mass radii relative to CDM. The entire population of dwarf galaxies produced in the ADM simulations is much more compact than those seen in CDM simulations, unable to reproduce the entire diversity of observed dwarf galaxy structures. Additionally, we also identify a population of highly compact subhalos that consist nearly entirely of ADM and form in the central region of the host, where they can leave distinctive imprints in the baryonic disk. This work presents the first detailed exploration of subhalo properties in a strongly dissipative dark matter scenario, providing intuition for how other regions of ADM parameter space, as well as other dark sector models, would impact galactic-scale observables.
Francesco Palma, Matteo Martucci, Coralie Neubüser, Alessandro Sotgiu, Francesco Maria Follega, Pietro Ubertini, Angela Bazzano, James Rodi, Roberto Ammendola, Davide Badoni, Simona Bartocci, Roberto Battiston, Stefania Beolè, Igor Bertello, William Jerome Burger, Donatella Campana, Antonio Cicone, Piero Cipollone, Silvia Coli, Livio Conti, et al (51) In this paper we report the detection of five strong Gamma-Ray Bursts (GRBs) by the High-Energy Particle Detector (HEPD-01) mounted on board the China Seismo-Electromagnetic Satellite (CSES-01), operational since 2018 on a Sun-synchronous polar orbit at a $\sim$ 507 km altitude and 97$^\circ$ inclination. HEPD-01 was designed to detect high-energy electrons in the energy range 3 - 100 MeV, protons in the range 30 - 300 MeV, and light nuclei in the range 30 - 300 MeV/n. Nonetheless, Monte Carlo simulations have shown HEPD-01 is sensitive to gamma-ray photons in the energy range 300 keV - 50 MeV, even if with a moderate effective area above $\sim$ 5 MeV. A dedicated time correlation analysis between GRBs reported in literature and signals from a set of HEPD-01 trigger configuration masks has confirmed the anticipated detector sensitivity to high-energy photons. A comparison between the simultaneous time profiles of HEPD-01 electron fluxes and photons from GRB190114C, GRB190305A, GRB190928A, GRB200826B and GRB211211A has shown a remarkable similarity, in spite of the different energy ranges. The high-energy response, with peak sensitivity at about 2 MeV, and moderate effective area of the detector in the actual flight configuration explain why these five GRBs, characterised by a fluence above $\sim$ 3 $\times$ 10$^{-5}$ erg cm$^{-2}$ in the energy interval 300 keV - 50 MeV, have been detected.
Enrico Garaldi, Rahul Kannan, Aaron Smith, Josh Borrow, Mark Vogelsberger, Rüdiger Pakmor, Volker Springel, Lars Hernquist, Daniela Galárraga-Espinosa, Jessica Y.-C. Yeh, Xuejian Shen, Clara Xu, Meredith Neyer, Benedetta Spina, Mouza Almualla, Yu Zhao Cosmological simulations serve as invaluable tools for understanding the Universe. However, the technical complexity and substantial computational resources required to generate such simulations often limit their accessibility within the broader research community. Notable exceptions exist, but most are not suited for simultaneously studying the physics of galaxy formation and cosmic reionization during the first billion years of cosmic history. This is especially relevant now that a fleet of advanced observatories (e.g. James Webb Space Telescope, Nancy Grace Roman Space Telescope, SPHEREx, ELT, SKA) will soon provide an holistic picture of this defining epoch. To bridge this gap, we publicly release all simulation outputs and post-processing products generated within the THESAN simulation project at https://thesan-project.com. This project focuses on the $z \geq 5.5$ Universe, combining a radiation-hydrodynamics solver (AREPO-RT), a well-tested galaxy formation model (IllustrisTNG) and cosmic dust physics to provide a comprehensive view of the Epoch of Reionization. The THESAN suite includes 16 distinct simulations, each varying in volume, resolution, and underlying physical models. This paper outlines the unique features of these new simulations, the production and detailed format of the wide range of derived data products, and the process for data retrieval. Finally, as a case study, we compare our simulation data with a number of recent observations from the James Webb Space Telescope, affirming the accuracy and applicability of THESAN. The examples also serve as prototypes for how to utilise the released dataset to perform comparisons between predictions and observations.
We present a comprehensive search and analysis of high-redshift galaxies in a suite of nine public JWST extragalactic fields taken in Cycle 1, covering a total effective search area of $\sim358{\rm arcmin^2}$. Through conservative ($8\sigma$) photometric selection, we identify 341 galaxies at $5<z<14$, with 109 having spectroscopic redshift measurements from the literature, including recent JWST NIRSpec observations. Our regression analysis reveals that the rest-frame UV size-stellar mass relation follows $R_{\rm eff}\propto M_*^{0.19\pm0.03}$, similar to that of star-forming galaxies at $z\sim3$, but scaled down in size by $\sim0.7$dex. We find a much slower rate for the average size evolution over the redshift range, $R_{\rm eff}\propto(1+z)^{-0.4\pm0.2}$, than that derived in the literature. A fraction ($\sim13\,\%$) of our sample are marginally resolved even in the NIRCam imaging ($<100$pc), located at $>1.5\,\sigma$ below the derived size-mass slope. These compact sources exhibit a high star formation surface density $\Sigma_{\rm SFR}>10\,M_\odot\,{\rm yr^{-1}\,kpc^{-2}}$, a range in which only $<0.01\,\%$ of the local star-forming galaxy sample is found. For those with available NIRSpec data, no evidence of ongoing supermassive black hole accretion is observed. A potential explanation for the observed high [OIII]-to-Hbeta ratios could be high shock velocities, likely originating within intense star-forming regions characterized by high $\Sigma_{\rm SFR}$. Lastly, we find that the rest-frame UV and optical sizes of our sample are comparable. Our results are consistent with these early galaxies building up their structures inside-out and yet to exhibit the strong color gradient seen at lower redshift.
Recent discoveries of a significant population of bright galaxies at cosmic dawn $\left(z \gtrsim 10\right)$ have enabled critical tests of cosmological galaxy formation models. In particular, the bright end of the galaxy UV luminosity function (UVLF) appears higher than predicted by many models. Using approximately 25,000 galaxy snapshots at $8 \leq z \leq 12$ in a suite of FIRE-2 cosmological "zoom-in'' simulations from the Feedback in Realistic Environments (FIRE) project, we show that the observed abundance of UV-bright galaxies at cosmic dawn is reproduced in these simulations with a multi-channel implementation of standard stellar feedback processes, without any fine-tuning. Notably, we find no need to invoke previously suggested modifications such as a non-standard cosmology, a top-heavy stellar initial mass function, or a strongly enhanced star formation efficiency. We contrast the UVLFs predicted by bursty star formation in these original simulations to those derived from star formation histories (SFHs) smoothed over prescribed timescales (e.g., 100 Myr). The comparison demonstrates that the strongly time-variable SFHs predicted by the FIRE simulations play a key role in correctly reproducing the observed, bright-end UVLFs at cosmic dawn: the bursty SFHs induce order-or-magnitude changes in the abundance of UV-bright ($M_\mathrm{UV} \lesssim -20$) galaxies at $z \gtrsim 10$. The predicted bright-end UVLFs are consistent with both the spectroscopically confirmed population and the photometrically selected candidates. We also find good agreement between the predicted and observationally inferred integrated UV luminosity densities, which evolve more weakly with redshift in FIRE than suggested by some other models.
Transient electron dynamics near the interface of counterstreaming plasmas at the onset of a relativistic collisionless shock (RCS) is investigated using particle-in-cell simulations. We identify a slingshot-like injection process induced by the drifting electric field sustained by the flowing focus of backwards-moving electrons, which is distinct from the well-known stochastic acceleration. The flowing focus signifies the plasma kinetic transition from a preturbulent laminar motion to a chaotic turbulence. We find a characteristic correlation between the electron dynamics in the slingshot acceleration and the photon emission features. In particular, the integrated radiation from the RCS exhibits a counterintuitive non-monotonic dependence of the photon polarization degree on the photon energy, which originates from a polarization degradation of relatively high-energy photons emitted by the slingshot-injected electrons. Our results demonstrate the potential of photon polarization as an essential information source in exploring intricate transient dynamics in RCSs with relevance for earth-based plasma and astrophysical scenarios.
Dongting Yang, Hong-Li Liu, Anandmayee Tej, Tie Liu, Patricio Sanhueza, Sheng-Li Qin, Xing Lu, Ke Wang, Sirong Pan, Feng-Wei Xu, Enrique Vazquez-Semadeni, Shanghuo Li, Gilberto C. Gomez, Aina Palau, Guido Garay, Paul F. Goldsmith, Mika Juvela, Anindya Saha, Leonardo Bronfman, Chang Won Lee, et al (11) There is growing evidence that high-mass star formation and hub-filament systems (HFS) are intricately linked. The gas kinematics along the filaments and the forming high-mass star(s) in the central hub are in excellent agreement with the new generation of global hierarchical high-mass star formation models. In this paper, we present an observational investigation of a typical HFS cloud, G310.142+0.758 (G310 hereafter) which reveals unambiguous evidence of mass inflow from the cloud scale via the filaments onto the forming protostar(s) at the hub conforming with the model predictions. Continuum and molecular line data from the ATOMS and MALT90 surveys are used that cover different spatial scales. Three filaments (with total mass $5.7\pm1.1\times 10^3~M_{\odot}$) are identified converging toward the central hub region where several signposts of high-mass star formation have been observed. The hub region contains a massive clump ($1280\pm260~M_{\odot}$) harbouring a central massive core. Additionally, five outflow lobes are associated with the central massive core implying a forming cluster. The observed large-scale, smooth and coherent velocity gradients from the cloud down to the core scale, and the signatures of infall motion seen in the central massive clump and core, clearly unveil a nearly-continuous, multi-scale mass accretion/transfer process at a similar mass infall rate of $\sim 10^{-3}~M_{\odot}~yr^{-1}$ over all scales, feeding the central forming high-mass protostar(s) in the G310 HFS cloud.
JWST observations have revealed a population of galaxies bright enough that potentially challenge standard galaxy formation models in the $\Lambda$CDM cosmology. Using a minimal empirical framework, we investigate the influence of variability on the rest-frame ultra-violet (UV) luminosity function (UVLF) of galaxies at $z\geq 9$. Our study differentiates between the $\textit{median UV radiation yield}$ and the $\textit{variability of UV luminosities}$ of galaxies at a fixed dark matter halo mass. We primarily focus on the latter effect, which depends on halo assembly and galaxy formation processes and can significantly increase the abundance of UV-bright galaxies due to the upscatter of galaxies in lower-mass haloes. We find that a relatively low level of variability, $\sigma_{\rm UV} \approx 0.75$ mag, matches the observational constraints at $z\approx 9$. However, increasingly larger $\sigma_{\rm UV}$ is necessary when moving to higher redshifts, reaching $\sigma_{\rm UV} \approx 2.0\,(2.5)\,{\rm mag}$ at $z\approx 12$ ($16$). This implied variability is consistent with expectations of physical processes in high-redshift galaxies such as bursty star formation and cycles of dust clearance. Photometric constraints from JWST at $z\gtrsim 9$ therefore can be reconciled with a standard $\Lambda$CDM-based galaxy formation model calibrated at lower redshifts without the need for adjustments to the median UV radiation yield.
Both observations and simulations have shown strong evidence for highly time-variable star formation in low-mass and/or high-redshift galaxies, which has important observational implications because high-redshift galaxy samples are rest-UV selected and therefore particularly sensitive to the recent star formation. Using a suite of cosmological "zoom-in" simulations at $z>5$ from the Feedback in Realistic Environments (FIRE) project, we examine the implications of bursty star formation histories for observations of high-redshift galaxies with JWST. We characterize how the galaxy observability depends on the star formation history. We also investigate selection effects due to bursty star formation on the physical properties measured, such as the gas fraction, specific star formation rate, and metallicity. We find the observability to be highly time-dependent for galaxies near the survey's limiting flux due to the SFR variability: as the star formation rate fluctuates, the same galaxy oscillates in and out of the observable sample. The observable fraction $f_\mathrm{obs} = 50\%$ at $z \sim 7$ and $M_{\star} \sim 10^{8.5}$ to $10^{9}\,M_{\odot}$ for a JWST/NIRCam survey reaching a limiting magnitude of $m^\mathrm{lim}_\mathrm{AB} \sim 29$-$30$, representative of surveys such as JADES and CEERS. JWST-detectable galaxies near the survey limit tend to have properties characteristic of galaxies in the bursty phase: on average, they show approximately 2.5 times higher cold, dense gas fractions and 20 times higher specific star formation rates at a given stellar mass than galaxies below the rest-UV detection threshold. Our study represents a first step in quantifying selection effects and the associated biases due to bursty star formation in studying high-redshift galaxy properties.
Dark sector theories naturally lead to multi-component scenarios for dark matter where a sub-component can dissipate energy through self-interactions, allowing it to efficiently cool inside galaxies. We present the first cosmological hydrodynamical simulations of Milky Way analogues where the majority of dark matter is collisionless Cold Dark Matter (CDM), but a sub-component (6%) is strongly dissipative minimal Atomic Dark Matter (ADM). The simulations, implemented in GIZMO and utilizing FIRE-2 galaxy formation physics to model the standard baryonic sector, demonstrate that the addition of even a small fraction of dissipative dark matter can significantly impact galactic evolution despite being consistent with current cosmological constraints. We show that ADM gas with roughly Standard-Model-like masses and couplings can cool to form a rotating "dark disk" with angular momentum closely aligned with the visible stellar disk. The morphology of the disk depends sensitively on the parameters of the ADM model, which affect the cooling rates in the dark sector. The majority of the ADM gas gravitationally collapses into dark "clumps" (regions of black hole or mirror star formation), which form a prominent bulge and a rotating thick disk in the central galaxy. These clumps form early and quickly sink to the inner ~kpc of the galaxy, affecting the galaxy's star-formation history and present-day baryonic and CDM distributions.
Xuejian Shen, Josh Borrow, Mark Vogelsberger, Enrico Garaldi, Aaron Smith, Rahul Kannan, Sandro Tacchella, Jesús Zavala, Lars Hernquist, Jessica Y.-C. Yeh, Chunyuan Zheng Using high-resolution cosmological radiation-hydrodynamic (RHD) simulations (THESAN-HR), we explore the impact of alternative dark matter (altDM) models on galaxies during the Epoch of Reionization. The simulations adopt the IllustrisTNG galaxy formation model. We focus on altDM models that exhibit small-scale suppression of the matter power spectrum, namely warm dark matter (WDM), fuzzy dark matter (FDM), and interacting dark matter (IDM) with strong dark acoustic oscillations (sDAO). In altDM scenarios, both the halo mass functions and the UV luminosity functions at $z\gtrsim 6$ are suppressed at the low-mass/faint end, leading to delayed global star formation and reionization histories. However, strong non-linear effects enable altDM models to "catch up" with cold dark matter (CDM) in terms of star formation and reionization. The specific star formation rates are enhanced in halos below the half-power mass in altDM models. This enhancement coincides with increased gas abundance, reduced gas depletion times, more compact galaxy sizes, and steeper metallicity gradients at the outskirts of the galaxies. These changes in galaxy properties can help disentangle altDM signatures from a range of astrophysical uncertainties. Meanwhile, it is the first time that altDM models have been studied in RHD simulations of galaxy formation. We uncover significant systematic uncertainties in reionization assumptions on the faint-end luminosity function. This underscores the necessity of accurately modeling the small-scale morphology of reionization in making predictions for the low-mass galaxy population. Upcoming James Webb Space Telescope (JWST) imaging surveys of deep, lensed fields hold potential for uncovering the faint, low-mass galaxy population, which could provide constraints on altDM models.
Philip F. Hopkins, Alexander B. Gurvich, Xuejian Shen, Zachary Hafen, Michael Y. Grudic, Shalini Kurinchi-Vendhan, Christopher C. Hayward, Fangzhou Jiang, Matthew E. Orr, Andrew Wetzel, Dusan Keres, Jonathan Stern, Claude-Andre Faucher-Giguere, James Bullock, Coral Wheeler, Kareem El-Badry, Sarah R. Loebman, Jorge Moreno, Michael Boylan-Kolchin, Eliot Quataert As they grow, galaxies can transition from irregular/spheroidal with 'bursty' star formation histories (SFHs), to disky with smooth SFHs. But even in simulations, the direct physical cause of such transitions remains unclear. We therefore explore this in a large suite of numerical experiments re-running portions of cosmological simulations with widely varied physics, further validated with existing FIRE simulations. We show that gas supply, cooling/thermodynamics, star formation model, Toomre scale, galaxy dynamical times, and feedback properties do not have a direct causal effect on these transitions. Rather, both the formation of disks and cessation of bursty star formation are driven by the gravitational potential, but in different ways. Disk formation is promoted when the mass profile becomes sufficiently centrally-concentrated in shape (relative to circularization radii): we show that this provides a well-defined dynamical center, ceases to support the global 'breathing modes' which can persist indefinitely in less-concentrated profiles and efficiently destroy disks, promotes orbit mixing to form a coherent angular momentum, and stabilizes the disk. Smooth SF is promoted by the potential or escape velocity (not circular velocity) becoming sufficiently large at the radii of star formation that cool, mass-loaded (momentum-conserving) outflows are trapped/confined near the galaxy, as opposed to escaping after bursts. We discuss the detailed physics, how these conditions arise in cosmological contexts, their relation to other correlated phenomena (e.g. inner halo virialization, vertical disk 'settling'), and observations.
Hong-Li Liu, Anandmayee Tej, Tie Liu, Patricio Sanhueza, Shengli Qin, Jinhua He, Paul F. Goldsmith, Guido Garay, Sirong Pan, Kaho Morii, Shanghuo Li, Amelia Stutz, Keníchi Tatematsu, Feng-Wei Xu, Leonardo Bronfman, Anindya Saha, Namitha Issac, Tapas Baug, L. Viktor Toth, Lokesh Dewangan, et al (13) We present a statistical study of a sample of 17 hub-filament-system (HFS) clouds of high-mass star formation using high-angular resolution ($\sim$1-2 arcsecond) ALMA 1.3mm and 3mm continuum data. The sample includes 8 infrared (IR)-dark and 9 IR-bright types, which correspond to an evolutionary sequence from the IR-dark to IR-bright stage. The central massive clumps and their associated most massive cores are observed to follow a trend of increasing mass ($M$) and mass surface density ($\Sigma$) with evolution from IR-dark to IR-bright stage. In addition, a mass-segregated cluster of young stellar objects (YSOs) are revealed in both IR-dark and IR-bright HFSs with massive YSOs located in the hub and the population of low-mass YSOs distributed over larger areas. Moreover, outflow feedback in all HFSs are found to escape preferentially through the inter-filamentary diffuse cavities, suggesting that outflows would render a limited effect on the disruption of the HFSs and ongoing high-mass star formation therein. From the above observations, we suggest that high-mass star formation in the HFSs can be described by a multi-scale mass accretion/transfer scenario, from hub-composing filaments through clumps down to cores, that can naturally lead to a mass-segregated cluster of stars.
Modeling self-gravity of collisionless fluids (e.g. ensembles of dark matter, stars, black holes, dust, planetary bodies) in simulations is challenging and requires some force softening. It is often desirable to allow softenings to evolve adaptively, in any high-dynamic range simulation, but this poses unique challenges of consistency, conservation, and accuracy, especially in multi-physics simulations where species with different softening laws may interact. We therefore derive a generalized form of the energy-and-momentum conserving gravitational equations of motion, applicable to arbitrary rules used to determine the force softening, together with consistent associated timestep criteria, interaction terms between species with different softening laws, and arbitrary maximum/minimum softenings. We also derive new methods to maintain better accuracy and conservation when symmetrizing forces between particles. We review and extend previously-discussed adaptive softening schemes based on the local neighbor particle density, and present several new schemes for scaling the softening with properties of the gravitational field, i.e. the potential or acceleration or tidal tensor. We show that the tidal softening scheme not only represents a physically-motivated, translation and Galilean invariant and equivalence-principle respecting (and therefore conservative) method, but imposes negligible timestep or other computational penalties, ensures that pairwise two-body scattering is small compared to smooth background forces, and can resolve outstanding challenges in properly capturing tidal disruption of substructures (minimizing artificial destruction) while also avoiding excessive N-body heating. We make all of this public in the GIZMO code.
Thomas K. Waters, Colton Peterson, Razieh Emami, Xuejian Shen, Lars Hernquist, Randall Smith, Mark Vogelsberger, Charles Alcock, Grant Tremblay, Matthew Liska, John C. Forbes, Jorge Moreno We present an in-depth analysis of gas morphologies for a sample of 25 Milky Way-like galaxies from the IllustrisTNG TNG50 simulation. We constrain the morphology of cold, warm, hot gas, and gas particles as a whole using a Local Shell Iterative Method (LSIM) and explore its observational implications by computing the hard-to-soft X-ray ratio, which ranges between $10^{-3}$-$10^{-2}$ in the inner $\sim 50 \rm kpc$ of the distribution and $10^{-5}$-$10^{-4}$ at the outer portion of the hot gas distribution. We group galaxies into three main categories: simple, stretched, and twisted. These categories are based on the radial reorientation of the principal axes of the reduced inertia tensor. We find that a vast majority ($77\%$) of the galaxies in our sample exhibit twisting patterns in their radial profiles. Additionally, we present detailed comparisons between 1) the gaseous distributions belonging to individual temperature regimes, 2) the cold gas distributions and stellar distributions, and 3) the gaseous distributions and dark matter (DM) halos. We find a strong correlation between the morphological properties of the cold gas and stellar distributions. Furthermore, we find a correlation between gaseous distributions with DM halo that increases with gas temperature, implying that we may use the warm-hot gaseous morphology as a tracer to probe the DM morphology. Finally, we show gaseous distributions exhibit significantly more prolate morphologies than the stellar distributions and DM halos, which we hypothesize is due to stellar and AGN feedback.
Many theories of dark matter beyond the Weakly Interacting Massive Particles (WIMP) paradigm feature an enhanced matter power spectrum on sub-parsec scales, leading to the formation of dense dark matter minihaloes. Future local observations are promising to search for and constrain such substructures. The survival probability of these dense minihaloes in the Milky Way environment is crucial for interpreting local observations. In this work, we investigate two environmental effects: stellar disruption and (smooth) tidal disruption. These two mechanisms are studied using semi-analytic models and idealized N-body simulations. For stellar disruption, we perform a series of N-body simulations of isolated minihalo-star encounters to test and calibrate analytic models of stellar encounters before applying the model to the realistic Milky Way disk environment. For tidal disruption, we perform N-body simulations to confirm the effectiveness of the analytic treatment. Finally, we propose a framework to combine the hierarchical assembly and infall of minihaloes to the Milky Way with the late-time disruption mechanisms. We make predictions for the mass functions of minihaloes in the Milky Way. The survival fraction of dense dark matter minihaloes, e.g. for axion miniclusters and minihaloes from Early Matter Domination, is $\sim 60\%$ with the relatively low-mass, compact population surviving. The survival fraction is insensitive to the detailed model parameters. We discuss various implications of the framework and future direct detection prospects.
Fangzhou Jiang, Andrew Benson, Philip F. Hopkins, Oren Slone, Mariangela Lisanti, Manoj Kaplinghat, Annika H. G. Peter, Zhichao Carton Zeng, Xiaolong Du, Shengqi Yang, Xuejian Shen We combine the isothermal Jeans model and the model of adiabatic halo contraction into a simple semi-analytic procedure for computing the density profile of self-interacting dark-matter (SIDM) haloes with the gravitational influence from the inhabitant galaxies. We show that the model agrees well with cosmological SIDM simulations over the entire core-forming stage and up to the onset of gravothermal core-collapse. Using this model, we show that the halo response to baryons is more diverse in SIDM than in CDM and depends sensitively on galaxy size, a desirable link in the context of the structural diversity of bright dwarf galaxies. The fast speed of the method facilitates analyses that would be challenging for numerical simulations -- notably, 1) we quantify the SIDM halo response as functions of the baryonic properties, on a fine mesh grid spanned by the baryon-to-total-mass ratio, $M_{\rm b}/M_{\rm vir}$, and galaxy compactness, $r_{1/2}/R_{\rm vir}$; 2) we show with high statistical precision that for typical Milky-Way-like systems, the SIDM profiles are similar to their CDM counterparts; and 3) we delineate the regime of gravothermal core-collapse in the $M_{\rm b}/M_{\rm vir}$-$r_{1/2}/R_{\rm vir}$ space, for a given cross section and a given halo concentration. Finally, we compare the isothermal Jeans model with the more sophisticated gravothermal fluid model, and show that the former yields faster core formation and agrees better with cosmological simulations. We attribute the difference to whether the target CDM halo is used as a boundary condition or as the initial condition for the gravothermal evolution, and thus comment on possible future improvements of the fluid model. We have made our programs for the model publicly available at https://github.com/JiangFangzhou/SIDM.
We analyze the first set of cosmological baryonic zoom-in simulations of galaxies in dissipative self-interacting dark matter (dSIDM). The simulations utilize the FIRE-2 galaxy formation physics with the inclusion of dissipative dark matter self-interactions modelled as a constant fractional energy dissipation ($f_{\rm diss}=0.5$). In this paper, we examine the properties of dwarf galaxies with $M_{\ast} \sim 10^{5}\operatorname{-}10^{9}\,{\rm M}_{\odot}$ in both isolation and within Milky Way-mass hosts. For isolated dwarfs, we find more compact galaxy sizes and promotion of stellar/neutral gas disk formation in dSIDM with $(\sigma/m)\leq 1\,{\rm cm^2\,g^{-1}}$ but they are still consistent with observed galaxy sizes and masses. In addition, as a result of the steeper central density profiles developed in dSIDM, the sub-kpc circular velocities of isolated dwarfs in models with $(\sigma/m)\geq 0.1\,{\rm cm^2\,g^{-1}}$ are enhanced by about a factor of two, which are still consistent with the measured stellar velocity dispersions of Local Group dwarfs but in tension with the HI rotation curves of more massive field dwarfs. Meanwhile, for satellites of the simulated Milky Way-mass hosts, the median circular velocity profiles are marginally affected by dSIDM physics, but dSIDM may help address the missing compact dwarf satellites in CDM. The number of satellites is slightly enhanced in dSIDM, but the differences are small compared with the large host-to-host variations. In conclusion, the dSIDM models with constant cross-section $(\sigma/m) \gtrsim 0.1\,{\rm cm^2\,g^{-1}}$ (assuming $f_{\rm diss}=0.5$) are effectively ruled out in bright dwarfs ($M_{\rm halo}\sim 10^{11}\,{\rm M}_{\odot}$) by circular velocity constraints. However, models with lower effective cross-sections (at this halo mass/velocity scale) are still viable and can give rise to non-trivial observable signatures.
Razieh Emami, Lars Hernquist, Mark Vogelsberger, Xuejian Shen, Joshua S. Speagle, Jorge Moreno, Charles Alcock, Shy Genel, John C. Forbes, Federico Marinacci, Paul Torrey We analyze the velocity anisotropy of stars in real and energy space for a sample of Milky Way-like galaxies in the TNG50 simulation. We employ different selection criteria, including spatial, kinematic and metallicity cuts, and make three halo classes ($\mathcal{A}$-$\mathcal{C}$) which show mild-to-strong sensitivity to different selections. The above classes cover 48%, 16% and 36% of halos, respectively. We analyze the $\beta$ radial profiles and divide them into either monotonically increasing radial profiles or ones with peaks and troughs. We demonstrate that halos with monotonically increasing $\beta$ profiles are mostly from class $\mathcal{A}$, whilst those with peaks/troughs are part of classes $\mathcal{B}$-$\mathcal{C}$. This means that care must be taken as the observationally reported peaks/troughs might be a consequence of different selection criteria. We infer the anisotropy parameter $\beta$ energy space and compare that against the $\beta$ radial profile. It is seen that 65% of halos with very mild sensitivity to different selections in real space, are those for which the $\beta$ radial and energy profiles are closely related. Consequently, we propose that comparing the $\beta$ radial and energy profiles might be a novel way to examine the sensitivity to different selection criteria and thus examining the robustness of the anisotropy parameter in tracing stellar kinematics. We compare simulated $\beta$ radial profiles against various observations and demonstrate that, in most cases, the model diversity is comparable with the error bars from different observations, meaning that the TNG50 models are in good overall agreement with observations.
We perform cosmological zoom-in simulations of $19$ relaxed cluster-mass haloes with the inclusion of adiabatic gas in the cold dark matter (CDM) and self-interacting dark matter (SIDM) models. These clusters are selected as dynamically relaxed clusters from a parent simulation with $M_{\rm 200} \simeq 1\operatorname{-}3\times 10^{15}\,{\rm M}_{\odot}$. Both the dark matter and the intracluster gas distributions in SIDM appear more spherical than their CDM counterparts. Mock X-ray images are generated based on the simulations and are compared to the real X-ray images of $84$ relaxed clusters selected from the Chandra and ROSAT archives. We perform ellipse fitting for the isophotes of mock and real X-ray images and obtain the ellipticities at cluster-centric radii of $r\simeq 0.1\operatorname{-}0.2\,R_{\rm 200}$. The X-ray isophotes in SIDM models with increasing cross-sections are rounder than their CDM counterparts, which manifests as a systematic shift in the distribution function of ellipticities. Unexpectedly, the X-ray morphology of the observed non-cool-core clusters agrees better with SIDM models with cross-section $(\sigma/m)= 0.5\operatorname{-}1~{\rm cm}^2/{\rm g}$ than CDM and SIDM with $(\sigma/m)=0.1\,{\rm cm}^2/{\rm g}$. Our statistical analysis indicates that the latter two models are disfavored at the $68\%$ confidence level (as conservative estimates). This conclusion is not altered by shifting the radial range of measurements or applying temperature selection criterion. However, the primary uncertainty originates from the lack of baryonic physics in the adiabatic model, such as cooling, star formation and feedback effects, which still have the potential to reconcile CDM simulations with observations.
In the time domain, the radio sky in particular along the Galactic plane direction may vary significantly because of various energetic activities associated with stars, stellar and supermassive black holes. Using multi-epoch Very Large Array surveys of the Galactic plane at 5.0 GHz, Becker et al. (2010) presented a catalogue of 39 variable radio sources in the flux density range 1-70 mJy. To probe their radio structures and spectra, we observed 17 sources with the very-long-baseline interferometric (VLBI) imaging technique and collected additional multi-frequency data from the literature. We detected all of the sources at 5 GHz with the Westerbork Synthesis Radio Telescope, but only G23.6644-0.0372 with the European VLBI Network (EVN). Together with its decadal variability and multi-frequency radio spectrum, we interpret it as an extragalactic peaked-spectrum source with a size of <~10 pc. The remaining sources were resolved out by the long baselines of the EVN because of either strong scatter broadening at the Galactic latitude <1 deg or intrinsically very extended structures on centi-arcsec scales. According to their spectral and structural properties, we find that the sample has a diverse nature. We notice two young H II regions and spot a radio star and a candidate planetary nebula. The rest of the sources are very likely associated with radio active galactic nuclei (AGN). Two of them also displays arcsec-scale faint jet activity. The sample study indicates that AGN are commonplace even among variable radio sources in the Galactic plane.
Line intensity mapping (LIM) is rapidly emerging as a powerful technique to study galaxy formation and cosmology in the high-redshift Universe. We present LIM estimates of select spectral lines originating from the interstellar medium (ISM) of galaxies and 21 cm emission from neutral hydrogen gas in the Universe using the large volume, high resolution THESAN reionization simulations. A combination of sub-resolution photo-ionization modelling for HII regions and Monte Carlo radiative transfer calculations is employed to estimate the dust-attenuated spectral energy distributions (SEDs) of high-redshift galaxies ($z\gtrsim5.5$). We show that the derived photometric properties such as the ultraviolet (UV) luminosity function and the UV continuum slopes match observationally inferred values, demonstrating the accuracy of the SED modelling. We provide fits to the luminosity--star formation rate relation (L-SFR) for the brightest emission lines and find that important differences exist between the derived scaling relations and the widely used low-$z$ ones because the interstellar medium of reionization era galaxies is generally less metal-enriched than in their low redshift counterparts. We use these relations to construct line intensity maps of nebular emission lines and cross correlate with the 21 cm emission. Interestingly, the wavenumber at which the correlation switches sign ($k_\mathrm{transition}$) depends heavily on the reionization model and to a lesser extent on the targeted emission line, which is consistent with the picture that $k_\mathrm{transition}$ probes the typical sizes of ionized regions. The derived scaling relations and intensity maps represent a timely state-of-the-art framework for forecasting and interpreting results from current and upcoming LIM experiments.
We post-process galaxies in the IllustrisTNG simulations with SKIRT radiative transfer calculations to make predictions for the rest-frame near-infrared (NIR) and far-infrared (FIR) properties of galaxies at $z\geq 4$. The rest-frame $K$- and $z$-band galaxy luminosity functions from TNG are overall consistent with observations, despite a $\sim 0.5\,\mathrm{dex}$ underprediction at $z=4$ for $M_{\rm K}\lesssim -25$ and $M_{\rm z}\lesssim -24$. Predictions for the JWST MIRI observed galaxy luminosity functions and number counts are given. Based on theoretical estimations, we show that the next-generation survey conducted by JWST can detect 500 (30) galaxies in F1000W in a survey area of $500\,{\rm arcmin}^{2}$ at $z=6$ ($z=8$). As opposed to the consistency in the UV, optical and NIR, we find that TNG, combined with our dust modelling choices, significantly underpredicts the abundance of most dust-obscured and thus most luminous FIR galaxies. As a result, the obscured cosmic star formation rate density (SFRD) and the SFRD contributed by optical/NIR dark objects are underpredicted. The discrepancies discovered here could provide new constraints on the sub-grid feedback models, or the dust contents, of simulations. Meanwhile, although the TNG predicted dust temperature and its relations with IR luminosity and redshift are qualitatively consistent with observations, the peak dust temperature of $z\geq 6$ galaxies are overestimated by about $20\,{\rm K}$. This could be related to the limited mass resolution of our simulations to fully resolve the porosity of the interstellar medium (or specifically its dust content) at these redshifts.
The existence of supermassive black holes (SMBHs) with masses greater than $\sim 10^{9}M_{\odot}$ at high redshift ($z\gtrsim 7$) is difficult to accommodate in standard astrophysical scenarios. We study the possibility that (nearly) totally dissipative self-interacting dark matter (tdSIDM)--in rare, high density dark matter fluctuations in the early Universe--produces SMBH seeds through catastrophic collapse. We use a semi-analytic model, tested and calibrated by a series of N-body simulations of isolated dark matter halos, to compute the collapse criteria and timescale of tdSIDM halos, where dark matter loses nearly all of its kinetic energy in a single collision in the center-of-momentum frame. Applying this model to halo merger trees, we empirically assign SMBH seeds to halos and trace the formation and evolution history of SMBHs. We make predictions for the quasar luminosity function, the $M_{\rm BH}$-$\sigma_{\rm v}^{\ast}$ relation, and cosmic SMBH mass density at high redshift and compare them to observations. We find that a dissipative dark matter interaction cross-section of $\sigma/m \sim 0.05~\rm cm^2/g$ is sufficient to produce the SMBHs observed in the early Universe while remaining consistent with ordinary SMBHs in the late Universe.
We present the first set of cosmological baryonic zoom-in simulations of galaxies including dissipative self-interacting dark matter (dSIDM). These simulations utilize the Feedback In Realistic Environments (FIRE-2) galaxy formation physics, but allow the dark matter to have dissipative self-interactions analogous to Standard Model forces, parameterized by the self-interaction cross-section per unit mass, $(\sigma/m)$, and the dimensionless degree of dissipation, $0<f_{\rm diss}<1$. We survey this parameter space, including constant and velocity-dependent cross-sections, and focus on structural and kinematic properties of dwarf galaxies with $M_{\rm halo} \simeq 10^{10-11} {\rm M}_{\odot}$. Central density profiles of simulated dwarfs become cuspy when $(\sigma/m)_{\rm eff} \gtrsim 0.1\,{\rm cm^{2}\,g^{-1}}$ (and $f_{\rm diss}=0.5$ as fiducial). The power-law slopes asymptote to $\alpha \approx -1.5$ in low-mass dwarfs independent of cross-section, which arises from a dark matter "cooling flow". Through comparisons with dark matter only simulations, we find the profile in this regime is insensitive to the inclusion of baryons. However, when $(\sigma/m)_{\rm eff} \ll 0.1\,{\rm cm^{2}\,g^{-1}}$, baryonic effects can produce cored density profiles comparable to non-dissipative cold dark matter (CDM) runs but at smaller radii. Simulated galaxies with $(\sigma/m) \gtrsim 10\,{\rm cm^{2}\,g^{-1}}$ develop significant coherent rotation of dark matter, accompanied by halo deformation, but this is unlike the well-defined thin "dark disks" often attributed to baryon-like dSIDM. The density profiles in this high cross-section model exhibit lower normalizations given the onset of halo deformation. For our surveyed dSIDM parameters, halo masses and galaxy stellar masses do not show appreciable difference from CDM, but dark matter kinematics and halo concentrations/shapes can differ.
Amongst the available technologies for earthquake research, remote sensing has been commonly used due to its unique features such as fast imaging and wide image-acquisition range. Nevertheless, early studies on pre-earthquake and remote-sensing anomalies are mostly oriented towards anomaly identification and analysis of a single physical parameter. Many analyses are based on singular events, which provide a lack of understanding of this complex natural phenomenon because usually, the earthquake signals are hidden in the environmental noise. The universality of such analysis still is not being demonstrated on a worldwide scale. In this paper, we investigate physical and dynamic changes of seismic data and thereby develop a novel machine learning method, namely Inverse Boosting Pruning Trees (IBPT), to issue short-term forecast based on the satellite data of 1,371 earthquakes of magnitude six or above due to their impact on the environment. We have analyzed and compared our proposed framework against several states of the art machine learning methods using ten different infrared and hyperspectral measurements collected between 2006 and 2013. Our proposed method outperforms all the six selected baselines and shows a strong capability in improving the likelihood of earthquake forecasting across different earthquake databases.
Razieh Emami, Lars Hernquist, Charles Alcock, Shy Genel, Sownak Bose, Rainer Weinberger, Mark Vogelsberger, Xuejian Shen, Joshua S. Speagle, Federico Marinacci, John C. Forbes, Paul Torrey We investigate the morphology of the stellar distribution in a sample of Milky Way (MW) like galaxies in the TNG50 simulation. Using a local in shell iterative method (LSIM) as the main approach, we explicitly show evidence of twisting (in about 52% of halos) and stretching (in 48% of them) in the real space. This is matched with the re-orientation observed in the eigenvectors of the inertia tensor and gives us a clear picture of having a re-oriented stellar distribution. We make a comparison between the shape profile of dark matter (DM) halo and stellar distribution and quite remarkably see that their radial profiles are fairly close, especially at small galactocentric radii where the stellar disk is located. This implies that the DM halo is somewhat aligned with stars in response to the baryonic potential. The level of alignment mostly decreases away from the center. We study the impact of substructures in the orbital circularity parameter. It is demonstrated that in some cases, far away substructures are counter-rotating compared with the central stars and may flip the sign of total angular momentum and thus the orbital circularity parameter. Truncating them above 150 kpc, however, retains the disky structure of the galaxy as per initial selection. Including the impact of substructures in the shape of stars, we explicitly show that their contribution is subdominant. Overlaying our theoretical results to the observational constraints from previous literature, we establish fair agreement.
We present predictions for high redshift ($z=2-10$) galaxy populations based on the IllustrisTNG simulation suite and a full Monte Carlo dust radiative transfer post-processing. Specifically, we discuss the ${\rm H}_{\alpha}$ and ${\rm H}_{\beta}$ + $[\rm O \,III]$ luminosity functions up to $z=8$. The predicted ${\rm H}_{\beta}$ + $[\rm O \,III]$ luminosity functions are consistent with present observations at $z\lesssim 3$ with $\lesssim 0.1\,{\rm dex}$ differences in luminosities. However, the predicted ${\rm H}_{\alpha}$ luminosity function is $\sim 0.3\,{\rm dex}$ dimmer than the observed one at $z\simeq 2$. Furthermore, we explore continuum spectral indices, the Balmer break at $4000$Å (D4000) and the UV continuum slope $\beta$. The median D4000 versus sSFR relation predicted at $z=2$ is in agreement with the local calibration despite a different distribution pattern of galaxies in this plane. In addition, we reproduce the observed $A_{\rm UV}$ versus $\beta$ relation and explore its dependence on galaxy stellar mass, providing an explanation for the observed complexity of this relation. We also find a deficiency in heavily attenuated, UV red galaxies in the simulations. Finally, we provide predictions for the dust attenuation curves of galaxies at $z=2-6$ and investigate their dependence on galaxy colors and stellar masses. The attenuation curves are steeper in galaxies at higher redshifts, with bluer colors, or with lower stellar masses. We attribute these predicted trends to dust geometry. Overall, our results are consistent with present observations of high redshift galaxies. Future JWST observations will further test these predictions.
In this paper, we provide updated constraints on the bolometric quasar luminosity function (QLF) from $z=0$ to $z=7$. The constraints are based on an observational compilation that includes observations in the rest-frame IR, B band, UV, soft and hard X-ray in past decades. Our method follows Hopkins et al. 2007 with an updated quasar SED model and bolometric and extinction corrections. The new best-fit bolometric quasar luminosity function behaves qualitatively different from the Hopkins et al. 2007 model at high redshift. Compared with the old model, the number density normalization decreases towards higher redshift and the bright-end slope is steeper at $z\gtrsim 2$. Due to the paucity of measurements at the faint end, the faint end slope at $z\gtrsim 5$ is quite uncertain. We present two models, one featuring a progressively steeper faint-end slope at higher redshift and the other featuring a shallow faint-end slope at $z\gtrsim 5$. Further multi-band observations of the faint-end QLF are needed to distinguish between these models. The evolutionary pattern of the bolometric QLF can be interpreted as an early phase likely dominated by the hierarchical assembly of structures and a late phase likely dominated by the quenching of galaxies. We explore the implications of this model on the ionizing photon production by quasars, the CXB spectrum, the SMBH mass density and mass functions. The predicted hydrogen photoionization rate contributed by quasars is subdominant during the epoch of reionization and only becomes important at $z\lesssim 3$. The predicted CXB spectrum, cosmic SMBH mass density and SMBH mass function are generally consistent with existing observations.
Philip Mocz, Anastasia Fialkov, Mark Vogelsberger, Fernando Becerra, Xuejian Shen, Victor H. Robles, Mustafa A. Amin, Jesús Zavala, Michael Boylan-Kolchin, Sownak Bose, Federico Marinacci, Pierre-Henri Chavanis, Lachlan Lancaster, Lars Hernquist Bose-Einstein Condensate Dark Matter (BECDM; also known as Fuzzy Dark Matter) is motivated by fundamental physics and has recently received significant attention as a serious alternative to the established Cold Dark Matter (CDM) model. We perform cosmological simulations of BECDM gravitationally coupled to baryons and investigate structure formation at high redshifts ($z \gtrsim 5$) for a boson mass $m=2.5\cdot 10^{-22}~{\rm eV}$, exploring the dynamical effects of its wavelike nature on the cosmic web and the formation of first galaxies. Our BECDM simulations are directly compared to CDM as well as to simulations where the dynamical quantum potential is ignored and only the initial suppression of the power spectrum is considered -- a Warm Dark Matter-like ("WDM") model often used as a proxy for BECDM. Our simulations confirm that "WDM" is a good approximation to BECDM on large cosmological scales even in the presence of the baryonic feedback. Similarities also exist on small scales, with primordial star formation happening both in isolated haloes and continuously along cosmic filaments; the latter effect is not present in CDM. Global star formation and metal enrichment in these first galaxies are delayed in BECDM/"WDM" compared to the CDM case: in BECDM/"WDM" first stars form at $z\sim 13$/$13.5$ while in CDM star formation starts at $z\sim 35$. The signature of BECDM interference, not present in "WDM", is seen in the evolved dark matter power spectrum: although the small scale structure is initially suppressed, power on kpc scales is added at lower redshifts. Our simulations lay the groundwork for realistic simulations of galaxy formation in BECDM.
Yunchong Wang, Mark Vogelsberger, Dandan Xu, Xuejian Shen, Shude Mao, David Barnes, Hui Li, Federico Marinacci, Paul Torrey, Volker Springel, Lars Hernquist We study the evolutionary trend of the total density profile of early-type galaxies (ETGs) in IllustrisTNG. To this end, we trace ETGs from $z=0$ to $z=4$ and measure the power-law slope $\gamma^{\prime}$ of the total density profile for their main progenitors. We find that their $\gamma^{\prime}$ steepen on average during $z\sim4-2$, then becoming shallower until $z=1$, after which they remain almost constant, aside from a residual trend of becoming shallower towards $z=0$. We also compare to a statistical sample of ETGs at different redshifts, selected based on their luminosity profiles and stellar masses. Due to different selection effects, the average slopes of the statistical samples follow a modified evolutionary trend. They monotonically decrease since $z=3$, and after $z\approx 1$, they remain nearly invariant with a mild increase towards $z=0$. These evolutionary trends are mass-dependent for both samples, with low-mass galaxies having in general steeper slopes than their more massive counterparts. Galaxies that transitioned to ETGs more recently have steeper mean slopes as they tend to be smaller and more compact at any given redshift. By analyzing the impact of mergers and AGN feedback on the progenitors' evolution, we conjecture a multi-phase path leading to isothermality in ETGs: dissipation associated with rapid wet mergers tends to steepen $\gamma^{\prime}$ from $z=4$ to $z=2$, whereas subsequent AGN feedback (especially in the kinetic mode) makes $\gamma^{\prime}$ shallower again from $z=2$ to $z=1$. Afterwards, passive evolution from $z=1$ to $z=0$, mainly through gas-poor mergers, mildly decreases $\gamma^{\prime}$ and maintains the overall mass distribution close to isothermal.
Mark Vogelsberger, Dylan Nelson, Annalisa Pillepich, Xuejian Shen, Federico Marinacci, Volker Springel, Ruediger Pakmor, Sandro Tacchella, Rainer Weinberger, Paul Torrey, Lars Hernquist The James Webb Space Telescop (JWST) promises to revolutionise our understanding of the early Universe, and contrasting its upcoming observations with predictions of the $\Lambda$CDM model requires detailed theoretical forecasts. Here, we exploit the large dynamic range of the IllustrisTNG simulation suite, TNG50, TNG100, and TNG300, to derive multi-band galaxy luminosity functions from $z=2$ to $z=10$. We put particular emphasis on the exploration of different dust attenuation models to determine galaxy luminosity functions for the rest-frame ultraviolet (UV), and apparent wide NIRCam bands. Our most detailed dust model is based on continuum Monte Carlo radiative transfer calculations employing observationally calibrated dust properties. This calibration results in constraints on the redshift evolution of the dust attenuation normalisation and dust-to-metal ratios yielding a stronger redshift evolution of the attenuation normalisation compared to most previous theoretical studies. Overall we find good agreement between the rest-frame UV luminosity functions and observational data for all redshifts, also beyond the regimes used for the dust-model calibrations. Furthermore, we also recover the observed high redshift ($z=4-6$) UV luminosity versus stellar mass relation, the H$\alpha$ versus star formation rate relation, and the H$\alpha$ luminosity function at $z=2$. The bright end ($M_{\rm UV}>-19.5$) cumulative galaxy number densities are consistent with observational data. For the F200W NIRCam band, we predict that JWST will detect $\sim 80$ ($\sim 200$) galaxies with a signal-to-noise ratio of $10$ ($\sim 5$) within the NIRCam field of view, $2.2\times2.2 \,{\rm arcmin}^{2}$, for a total exposure time of $10^5{\rm s}$ in the redshift range $z=8 \pm 0.5$. These numbers drop to $\sim 10$ ($\sim 40$) for an exposure time of $10^4{\rm s}$.
We contrast predictions for the high-redshift galaxy population and reionization history between cold dark matter (CDM) and an alternative self-interacting dark matter model based on the recently developed ETHOS framework that alleviates the small-scale CDM challenges within the Local Group. We perform the highest resolution hydrodynamical cosmological simulations (a 36~Mpc$^3$ volume with gas cell mass of $\sim10^5\mathrm{M}_{\odot}$ and minimum gas softening of $\sim180$~pc) within ETHOS to date -- plus a CDM counterpart -- to quantify the abundance of galaxies at high redshift and their impact on reionization. We find that ETHOS predicts galaxies with higher ultraviolet (UV) luminosities than their CDM counterparts and a faster build-up of the faint end of the UV luminosity function. These effects, however, make the optical depth to reionization less sensitive to the power spectrum cut-off: the ETHOS model differs from the CDM $\tau$ value by only 10 per cent and is consistent with Planck limits if the effective escape fraction of UV photons is 0.1-0.5. We conclude that current observations of high-redshift luminosity functions cannot differentiate between ETHOS and CDM models, but deep JWST surveys of strongly-lensed, inherently faint galaxies have the potential to test non-CDM models that offer attractive solutions to CDM's Local Group problems.
Radiative emission lines from nitrogen and its ions are often observed in nebulae spectra, where the N$^{2+}$ abundance can be inferred from lines of the 2p4f configuration. In addition, intensity ratios between lines of the 2p3p -- 2p3s and 2p4f -- 2p3d transition arrays can serve as temperature diagnostics. To aid abundance determinations and plasma diagnostics, wavelengths and oscillator strengths were calculated with high-precision for electric-dipole (E1) transitions from levels in the 2p4f configuration of N$^{+}$. Electron correlation and relativistic effects, including the Breit interaction, were systematically taken into account within the framework of the multiconfiguration Dirac-Hartree-Fock (MCDHF) method. Except for the 2p4f - 2p4d transitions with quite large wavelengths and the two-electron-one-photon 2p4f -2s2p$^3$ transitions, the uncertainties of the present calculations were controlled to within 3% and 5% for wavelengths and oscillator strengths, respectively. We also compared our results with other theoretical and experimental values when available. Discrepancies were found between our calculations and previous calculations due to the neglect of relativistic effects in the latter.
This work describes a study of elemental abundances for 30 metal-poor stars whose chemical abundances provide excellent information for setting constraints on models of neutron-capture processes. Based on the abundances of main r- and weak r-process stars, the abundance patterns of main r-process and weak r-process are obtained. The two r-process component coefficients are defined to determine the relative contributions from individual neutron-capture process to abundances of metal-poor stars. Based on the component coefficients, we find that metal-poor stars BD+4 2621 and HD 4306 are also weak r-process stars, which means that the abundance pattern produced by weak r-process is stable. All metal-poor star abundances contain the contributions of both main r-process and weak r-process. The elements produced by weak r-process have increased along with Fe over the polluted history. Most of the metal-poor star abundances do not follow the pattern observed in solar system, but there is a small fraction that do. For the low-[Sr/Fe] star BD-18 5550 ([Sr/Fe]$\lesssim-1$), neutron-capture element abundances can be explained by the mixture of two r-process components. Since lighter elements in this star cannot be fitted by the two components, the abundance pattern of P-component is estimated from those abundances.
In this paper, we propose a new class of parametrization of the equation of state of dark energy. In contrast with the famous CPL parametrization, these new parametrization of the equation of state does not divergent during the evolution of the Universe even in the future. Also, we perform a observational constraint on two simplest dark energy models belonging to this new class of parametrization, by using the Markov Chain Monte Carlo (MCMC) method and the combined latest observational data from the type Ia supernova compilations including Union2(557), cosmic microwave background, and baryon acoustic oscillation.
Recently, the vacuum energy of the QCD ghost in a time-dependent background is proposed as a kind of dark energy candidate to explain the acceleration of the universe. In this model, the energy density of the dark energy is proportional to the Hubble parameter $H$, which is the Hawking temperature on the Hubble horizon of the Friedmann-Robertson-Walker (FRW) universe. In this paper, we perform a constraint on the ghost dark energy model with and without bulk viscosity, by using the Markov Chain Monte Carlo (MCMC) method and the combined latest observational data from the type Ia supernova compilations including Union2.1(580) and Union2(557), cosmic microwave background, baryon acoustic oscillation, and the observational Hubble parameter data.
Recently, the vacuum energy of the QCD ghost in a time-dependent background is proposed as a kind of dark energy candidate to explain the acceleration of the Universe. In this model, the energy density of the dark energy is proportional to the Hubble parameter $H$, which is the Hawking temperature on the Hubble horizon of the Friedmann-Robertson-Walker (FRW) Universe. In this paper, we generalized this model and choice the Hawking temperature on the so-called trapping horizon, which will coincides with the Hubble temperature in the context of flat FRW Universe dominated by the dark energy component. We study the thermodynamics of Universe with this kind of dark energy and find that the entropy-area relation is modified, namely, there is an another new term besides the area term.
Correlated imaging through atmospheric turbulence is studied, and the analytical expressions describing turbulence effects on image resolution are derived. Compared with direct imaging, correlated imaging can reduce the influence of turbulence to a certain extent and reconstruct high-resolution images. The result is backed up by numerical simulations, in which turbulence-induced phase perturbations are simulated by random phase screens inserting propagation paths.