Search

Modern high-resolution simulations of the interstellar medium (ISM) have shown that key factors in governing star formation are the competing influences of radiative dissipation, pressure support driven by stellar feedback, and the relentless pull of gravity. Cosmological simulations of galaxy formation, such as IllustrisTNG or ASTRID, are however not able to resolve this physics in detail and therefore need to rely on approximate treatments. These have often taken the form of empirical subgrid models of the ISM expressed in terms of an effective equation of state (EOS) that relates the mean ISM pressure to the mean gas density. Here we seek to improve these heuristic models by directly fitting their key ingredients to results of the high-resolution TIGRESS simulations, which have shown that the dynamical equilibrium of the ISM can be understood in terms of a pressure-regulated, feedback modulated (PRFM) model for star formation. Here we explore a simple subgrid model that draws on the PRFM concept but uses only local quantities. It accurately reproduces PRFM for pure gas disks, while it predicts slightly less star formation than PRFM in the presence of an additional thin stellar disk. We compare the properties of this model with the older Springel and Hernquist and TNG prescriptions, and apply all three to isolated simulations of disk galaxies as well as to a set of high-resolution zoom-in simulations carried out with a novel 'multi-zoom' technique that we introduce in this study. The softer EOS implied by TIGRESS produces substantially thinner disk galaxies, which has important ramifications for disk stability and galaxy morphology. The total stellar mass of galaxies is however hardly modified at low redshift, reflecting the dominating influence of large-scale gaseous inflows and outflows to galaxies, which are not sensitive to the EOS itself, Comment: 22 pages, 21 figures, to be submitted to MNRAS. This is a Learning the Universe publication

Galactic double white dwarfs will be prominent gravitational-wave sources for the Laser Interferometer Space Antenna (LISA). While previous studies have primarily focused on formation scenarios in which binaries form and evolve in isolation, we present the first detailed study of the role of triple stellar evolution in forming the population of LISA double white dwarfs. In this work, we present the first detailed study of the role of triple stellar evolution in forming the population of LISA double white dwarfs. We use the multiple stellar evolution code (MSE) to model the stellar evolution, binary interactions, and the dynamics of triple star systems then use a Milky Way-like galaxy from the TNG50 simulations to construct a representative sample of LISA double white dwarfs. In our simulations about $7\times10^6$ Galactic double white dwarfs in the LISA frequency bandwidth originate from triple systems, whereas $\sim4\times10^6$ form from isolated binary stars. The properties of double white dwarfs formed in triples closely resemble those formed from isolated binaries, but we also find a small number of systems $\sim\mathcal{O}(10)$ that reach extreme eccentricities $(>0.9)$, a feature unique to the dynamical formation channels. Our population produces $\approx 10^{4} $ individually resolved double white dwarfs (from triple and binary channels) and an unresolved stochastic foreground below the level of the LISA instrumental noise. About $57\,\%$ of double white dwarfs from triple systems retain a bound third star when entering the LISA frequency bandwidth. However, we expect the tertiary stars to be too distant to have a detectable imprint in the gravitational-wave signal of the inner binary., Comment: 16 pages, 13 figures

The circumgalactic medium (CGM) plays a critical role in galaxy evolution, influencing gas flows, feedback processes, and galactic dynamics. Observations show a substantial cold gas reservoir in the CGM, but the mechanisms driving its formation and evolution remain unclear. Cosmic rays (CRs), as a source of non-thermal pressure, are increasingly recognized as key regulators of cold gas dynamics. This study explores how CRs affect cold clouds that condense from the hot CGM via thermal instability (TI). Using 3D CR-magnetohydrodynamic (CRMHD) simulations with AREPO, we assess the impact of various CR transport models on cold gas evolution. Under purely advective CR transport, CR pressure significantly suppresses the collapse of thermally unstable regions, altering the CGM's structure. In contrast, realistic CR transport models reveal that CRs escape collapsing regions via streaming and diffusion along magnetic fields, diminishing their influence on the thermal and dynamic structure of the cold CGM. The ratio of the CR escape time to the cloud collapse time emerges as a critical factor in determining the impact of CRs on TI. CRs remain confined within cold clouds when effective CR diffusion is slow which maximizes their pressure support and inhibits collapse. Fast effective CR diffusion, as realized in our 2-moment CRMHD model, facilitates rapid CR escape, reducing their stabilizing effect. This realistic CR transport model shows a wide dynamic range of the effective CR diffusion coefficient, ranging from $10^{29}$ to $10^{30}\,\mathrm{cm^{2}\,s^{-1}}$ for thermally- to CR-dominated atmospheres, respectively. In addition to these CR transport-related effects, we demonstrate that high numerical resolution is crucial to avoid spuriously large clouds formed in low-resolution simulations, which would result in overly long CR escape times and artificially amplified CR pressure support., Comment: 20 pages and 17 figures in main text; submitted to A&A

The stellar haloes of dwarf galaxies are becoming an object of interest in the extragalactic community due to their detection in some recent observations. Additionally, new cosmological simulations of very high resolution were performed, allowing their study. These stellar haloes could help shed light on our understanding of the assembly of dwarf galaxies and their evolution, and allow us to test the hierarchical model for the formation of structures at small scales. We aim to characterise the stellar haloes of simulated dwarf galaxies and analyse their evolution and accretion history. We use a sample of 17 simulated galaxies from the Auriga Project with a stellar mass range from 3.28x10^8 Msun to 2.08x10^10 Msun. We define the stellar halo as the stellar material located outside an ellipsoid with semi-major axes equal to 4 times the half light radius (Rh) of each galaxy. We find that the inner regions of the stellar halo (4 to 6 times the Rh) are dominated by in-situ material. For the less massive simulated dwarfs (M*<=4.54x10^8 Msun), this dominance extends to all radii. We find that this in-situ stellar halo is mostly formed in the inner regions of the galaxies and then ejected into the outskirts during interactions and merger events. In ~50% of the galaxies, the stripped gas from satellites contributed to the formation of this in-situ halo. The stellar haloes of the galaxies more massive than M*>=1x10^9 Msun are dominated by the accreted component beyond 6 Rh. We find that the more massive dwarf galaxies accrete stellar material until later times (t90~4.44 Gyr ago, being t90 the formation time) than the less massive ones (t90~8.17 Gyr ago), impacting on the formation time of the accreted stellar haloes. The galaxies have a range of 1 to 7 significant progenitors contributing to their accreted component but there is no correlation between this quantity and the galaxies' accreted mass., Comment: Submitted to A&A. 21 pages, 1 table, 15 figures and 2 appendices

The double detonation is a widely discussed explosion mechanism for Type Ia supernovae, whereby a helium shell detonation ignites a secondary detonation in the carbon/oxygen core of a white dwarf. Even for modern models that invoke relatively small He shell masses, many previous studies have found that the products of the helium shell detonation lead to discrepancies with normal Type Ia supernovae, such as strong Ti II absorption features, extremely red light curves and too large a variation with viewing direction. It has been suggested that non local thermodynamic equilibrium (non-LTE) effects may help to reduce these discrepancies with observations. Here we carry out full non-LTE radiative transfer simulations for a recent double detonation model with a relatively small helium shell mass of 0.05 M$_\odot$. We construct 1D models representative of directions in a 3D explosion model to give an indication of viewing angle dependence. The full non-LTE treatment leads to improved agreement between the models and observations. The light curves become less red, due to reduced absorption by the helium shell detonation products, since these species are more highly ionised. Additionally, the expected variation with observer direction is reduced. The full non-LTE treatment shows promising improvements, and reduces the discrepancies between the double detonation models and observations of normal Type Ia supernovae., Comment: Submitted to MNRAS

Compact binaries containing hot subdwarfs and white dwarfs have the potential to evolve into a variety of explosive transients. These systems could also explain hypervelocity runaway stars such as US 708. We use the detailed binary evolution code MESA to evolve hot subdwarf and white dwarf stars interacting in binaries. We explore their evolution towards double detonation supernovae, helium novae, or double white dwarfs. Our grid of 3120 models maps from initial conditions such as orbital period and masses of hot subdwarf and white dwarf to these outcomes. The minimum amount of helium required to ignite the helium shell that leads to a double detonation supernova in our grid is $\approx 0.05 \, \mathrm{M_{\odot}}$, likely too large to produce spectra similar to normal type Ia supernovae, but compatible with inferred helium shell masses from some observed peculiar type I supernovae. We also provide the helium shell masses for our double white dwarf systems, with a maximum He shell mass of $\approx 0.18\,\mathrm{M_{\odot}}$. In our double detonation systems, the orbital velocity of the surviving donor star ranges from $\approx 450 \, \mathrm{km\,s^{-1}}$ to $\approx 1000 \, \mathrm{km\,s^{-1}}$. Among the surviving donors, we also estimate the runaway velocities of proto-white dwarfs, which have higher runaway velocities than hot subdwarf stars of the same mass. Our grid will provide a first-order estimate of the potential outcomes for the observation of binaries containing hot subdwarfs and white dwarfs from future missions like Gaia, LSST, and LISA., Comment: 15 pages, 11 figures, comments welcome

Observations of Faraday rotation and synchrotron emission in galaxy clusters imply large-scale magnetic fields with $\mu\mathrm{G}$ strengths possibly extending back to $z=4$. Non-radiative cosmological simulations of galaxy clusters show a comparably slow magnetic field growth that only saturates at late times. We include galaxy formation physics and find a significantly accelerated magnetic field growth. We identify three crucial stages in the magnetic field evolution. 1) At high redshift, the central dominant galaxy serves as the prime agent that magnetizes not only its immediate vicinity but also most of the forming protocluster through a combination of a small-scale dynamo induced by gravitationally driven, compressive turbulence and stellar and active galactic nuclei feedback that distribute the magnetic field via outflows. 2) This process continues as other galaxies merge into the forming cluster in subsequent epochs, thereby transporting their previously amplified magnetic field to the intracluster medium (ICM) through ram pressure stripping and galactic winds. 3) At lower redshift, gas accretion and frequent cluster mergers trigger additional small-scale dynamo processes, thereby preventing the decay of the magnetic field and fostering the increase of the magnetic coherence scale. We show that the magnetic field observed today in the weakly collisional ICM is consistently amplified on collisional scales. Initially, this occurs in the collisional interstellar medium during protocluster assembly, and later in the ICM on the magnetic coherence scale, which always exceeds the particle mean free path, supporting the use of magneto-hydrodynamics for studying the cluster dynamo. We generate synthetic Faraday rotation measure observations of protoclusters, highlighting the potential for studying magnetic field growth during the onset of cluster formation at cosmic dawn., Comment: 23 pages, 15 figures, submitted to A&A

Hypervelocity white dwarfs (HVWDs) are stellar remnants moving at speeds exceeding the Milky Way's escape velocity. The origins of the fastest HVWDs are enigmatic, with proposed formation scenarios facing challenges explaining both their extreme velocities and observed properties. Here we report a three-dimensional hydrodynamic simulation of a merger between two hybrid helium-carbon-oxygen white dwarfs (HeCO WDs with masses of 0.68 and 0.62 M$_\odot$). We find that the merger leads to a partial disruption of the secondary WD, coupled with a double-detonation explosion of the primary WD. This launches the remnant core of the secondary WD at a speed of $~2000$ km s$^{-1}$, consistent with observed HVWDs. The low mass of the ejected remnant and its heating from the primary WD's ejecta explain the observed luminosities and temperatures of hot HVWDs, which are otherwise difficult to reconcile with previous models. This discovery establishes a new formation channel for HVWDs and points to a previously unrecognized pathway for producing peculiar Type Ia supernovae and faint explosive transients., Comment: 13 pages, 6 figures. Comments welcome!

In weakly collisional, strongly magnetised plasmas such as the intracluster medium (ICM), hot accretion flows and the solar corona, the transport of heat and momentum occurs primarily along magnetic field lines. In this paper we present a new scheme for modelling anisotropic thermal conduction which we have implemented in the moving mesh code AREPO. Our implementation uses a semi-implicit time integration scheme which works accurately and efficiently with individual timestepping, making the scheme highly suitable for use in cosmological simulations. We apply the scheme to a number of test-problems including the diffusion of a hot patch of gas in a circular magnetic field, the progression of a point explosion in the presence of thermal conduction, and the evolution and saturation of buoyancy instabilities in anisotropically conducting plasmas. We use these idealised tests to demonstrate the accuracy and stability of the solver and highlight the ways in which anisotropic conduction can fundamentally change the behaviour of the system. Finally, we demonstrate the solver's capability when applied to highly non-linear problems with deep timestep hierarchies by performing high-resolution cosmological zoom-in simulations of a galaxy cluster with conduction. We show that anisotropic thermal conduction can have a significant impact on the temperature distribution of the ICM and that whistler suppression may be relevant on cluster scales. The new scheme is, therefore, well suited for future work which will explore the role of anisotropic thermal conduction in a range of astrophysical contexts including the ICM of clusters and the circumgalactic medium of galaxies., Comment: 20 pages, 14 figures, submitted to MNRAS, comments welcome

Magnetic fields have been shown both observationally and through theoretical work to be an important factor in the formation of protostars and their accretion disks. Accurate modelling of the evolution of the magnetic field in low-ionization molecular cloud cores requires the inclusion of non-ideal magnetohydrodynamics (MHD) processes, specifically Ohmic and ambipolar diffusion and the Hall effect. These have a profound influence on the efficiency of magnetic removal of angular momentum from protostellar disks and simulations that include them can avoid the `magnetic-braking catastrophe' in which disks are not able to form. However, the impact of the Hall effect, in particular, is complex and remains poorly studied. In this work, we perform a large suite of simulations of the collapse of cloud cores to protostars with several non-ideal MHD chemistry models and initial core geometries using the moving-mesh code {\small AREPO}. We find that the efficiency of angular momentum removal is significantly reduced with respect to ideal MHD, in line with previous results. The Hall effect has a varied influence on the evolution of the disk which depends on the initial orientation of the magnetic field. This extends to the outflows seen in a subset of the models, where this effect can act to enhance or suppress them and open up new outflow channels. We conclude, in agreement with a subset of the previous literature, that the Hall effect is the dominant non-ideal MHD process in some collapse scenarios and thus should be included in simulations of protostellar disk formation., Comment: 23 pages, 20 Figures; version accepted by MNRAS

Blue straggler stars in stellar clusters are a subset of stars that are bluer and appear younger than other cluster members, seemingly straggling behind in their evolution. They offer a unique opportunity to understand the stellar dynamics and populations within their hosts. In the collisional formation scenario, a persistent challenge is the excessive angular momentum in the collision product. The consequent significant mass loss during transition to a stable state leads to a star with too low a mass to be a blue straggler, unless it spins down efficiently. While many proposed spin-down mechanisms involve boosted angular momentum loss via magnetic braking within the collision product, the existence or strength of these magnetic fields has not been confirmed. Here, we report three-dimensional magnetohydrodynamical simulations of collisions between two low-mass main-sequence stars and investigate magnetic field amplification. Magnetic field energy is amplified by a factor of $10^{8}-10^{10}$, resulting in the magnetic field strength of $10^{7}-10^{8}$G at the core of the collision product, independent of collision parameters. The surface magnetic field strengths have increased up to $10-10^{4}$ G. In addition, a distinctly flattened, rotating gas structure appears around the collision products in off-axis collisions, which may be a hint of possible disk formation. Such significant magnetic amplification and potential disk formation suggest the possibility of efficient spin-down of collision products via magnetic braking and magnetic disk locking, which can result in their appearance as blue stragglers., Comment: 15 pages, 11 figures, 1 table. Accepted for publication in ApJL. Videos: https://www.youtube.com/watch?v=lkpTEJTI1Ew&list=PLxLK3qI02cQef6wX5wi-QMSFyxBGK4tGL

Magnetic fields are ubiquitous in the universe and an important component of the interstellar medium. It is crucial to accurately model and understand their properties in different environments and across all mass ranges to interpret observables related to magnetic fields correctly. However, the assessment of the role of magnetic fields in galaxy evolution is often hampered by limited numerical resolution in cosmological simulations, in particular for satellite galaxies. To this end, we study the magnetic fields in high-resolution cosmological zoom simulations of disk galaxies (with $M_{200}\approx10^{10}$ to $10^{13}\,\mathrm{M_\odot}$) and their satellites within the Auriga galaxy formation model including cosmic rays. We find significantly higher magnetic field strengths in satellite galaxies compared to isolated dwarfs with a similar mass or star-formation rate, in particular after they had their first close encounter with their host galaxy. While this result is ubiquitous and independent of resolution in the satellites that are past their first infall, there seems to be a wide range of amplification mechanisms acting together. Our result highlights the importance of considering the environment of dwarf galaxies when interpreting their magnetic field properties as well as related observables such as their gamma-ray and radio emission, the latter being particularly relevant for future observations such as with the SKA observatory., Comment: 9 pages, 5 figures, submitted to MNRAS. Comments are welcome!

The hot phase of the circumgalactic medium (CGM) allows us to probe the inflow and outflow of gas within a galaxy, which is responsible for dictating the evolution of the galaxy. Studying the hot CGM sheds light on a better understanding of gas physics, which is crucial to inform and constrain simulation models. With the recent advances in observational measurements probing the hot CGM in X-rays and tSZ, we have a new avenue for widening our knowledge of gas physics and feedback by exploiting the information from current/future observations. In this paper, we use the TNG300 hydrodynamical simulations to build a fully self-consistent forward model for the hot CGM. We construct a lightcone and generate mock X-ray observations. We quantify the projection effects, namely the locally correlated large-scale structure in X-rays and the effect due to satellite galaxies misclassified as centrals which affects the measured hot CGM galactocentric profiles in stacking experiments. We present an analytical model that describes the intrinsic X-ray surface brightness profile across the stellar and halo mass bins. The increasing stellar mass bins result in decreasing values of $\beta$, the exponent quantifying the slope of the intrinsic galactocentric profiles. We carry forward the current state-of-the-art by also showing the impact of the locally correlated environment on the measured X-ray surface brightness profiles. We also present, for the first time, the effect of misclassified centrals in stacking experiments for three stellar mass bins: $10^{10.5-11}\ M_\odot$, $10^{11-11.2}\ M_\odot$, and $10^{11.2-11.5}\ M_\odot$. We find that the contaminating effect of the misclassified centrals on the stacked profiles increases when the stellar mass decreases., Comment: 14 pages, 10 figures, Submitted to A&A, comments welcome

We study the dependence of the physical and observable properties of the CGM on its halo mass. We analyse 22 simulations from the Auriga suite of high resolution cosmological `zoom-in' simulations at $z=0$ with halo masses $10^{10}~\text{M}_{\odot}\leq\text{M}_{\mathrm{200c}}\leq10^{12}~\text{M}_{\odot}$. We find a larger scatter in temperature and smaller scatter in metallicity in more massive haloes. The scatter of temperature and metallicity as a function of radius increases out to larger radii. The median and scatter of the volume-weighted density and mass-weighted radial velocity show no significant dependence on halo mass. Our results highlight that the CGM is more multiphase in haloes of higher mass. We additionally investigate column densities for HI and the metal ions CIV, OVI, MgII and SiII as a function of stellar mass and radius. We find the HI and metal ion column densities increase with stellar mass at sufficiently large radii ($R\gtrsim{0.2}$R$_{\mathrm{200c}}$). We find good agreement between our HI column densities and observations outside $20$% of the virial radius and overpredict within $20$%. MgII and SiII are similarly overpredicted within $20$% of the virial radius, but drop off steeply at larger radii. Our OVI column densities underpredict observations for stellar masses between $10^{9.7}~\text{M}_{\odot}\leq\text{M}_{\star}<10^{10.8}~\text{M}_{\odot}$ with reasonable agreement at $10^{10.8}~\text{M}_{\odot}$. CIV column densities agree with observational detections above a halo mass of $10^{9.7}~\text{M}_{\odot}$. We find that OVI (MgII) traces the highest (lowest) temperatures, and lowest (highest) density and metallicity. OVI and CIV are photo-ionized (collisionally ionized) at low (high) halo masses with a transition to higher temperatures at $10^{11}~\text{M}_{\odot}$. However, there is no clear trend for the radial velocity of the ions., Comment: 15 pages, 8 figures. Submitted for publication in MNRAS

We investigate the redshift evolution of the concentration-mass relationship of dark matter haloes in state-of-the-art cosmological hydrodynamic simulations and their dark-matter-only counterparts. By combining the IllustrisTNG suite and the novel MillenniumTNG simulation, our analysis encompasses a wide range of box size ($50 - 740 \: \rm cMpc$) and mass resolution ($8.5 \times 10^4 - 3.1 \times 10^7 \: \rm M_{\odot}$ per baryonic mass element). This enables us to study the impact of baryons on the concentration-mass relationship in the redshift interval $0

We discuss the cold dark matter plus massive neutrinos simulations of the MillenniumTNG (MTNG) project, which aim to improve understanding of how well ongoing and future large-scale galaxy surveys will measure neutrino masses. Our largest simulations, $3000\,{\rm Mpc}$ on a side, use $10240^3$ particles of mass $m_{p} = 6.66\times 10^{8}\,h^{-1}{\rm M}_\odot$ to represent cold dark matter, and $2560^3$ to represent a population of neutrinos with summed mass $M_\nu = 100\,{\rm meV}$. Smaller volume runs with $\sim 630\,{\rm Mpc}$ also include cases with $M_\nu = 0\,\textrm{and}\, 300\,{\rm meV}$. All simulations are carried out twice using the paired-and-fixed technique for cosmic variance reduction. We evolve the neutrino component using the particle-based $\delta f$ importance sampling method, which greatly reduces shot noise in the neutrino density field. In addition, we modify the GADGET-4 code to account both for the influence of relativistic and mildly relativistic components on the expansion rate and for non-Newtonian effects on the largest represented simulation scales. This allows us to quantify accurately the impact of neutrinos on basic statistical measures of nonlinear structure formation, such as the matter power spectrum and the halo mass function. We use semi-analytic models of galaxy formation to predict the galaxy population and its clustering properties as a function of summed neutrino mass, finding significant ($\sim 10\%$) impacts on the cosmic star formation rate history, the galaxy mass function, and the clustering strength. This offers the prospect of identifying combinations of summary statistics that are optimally sensitive to the neutrino mass., Comment: 19 pages, 15 figures

The nature of progenitors of Type Ia supernovae has long been debated, primarily due to the elusiveness of the progenitor systems to traditional electromagnetic observation methods. We argue that gravitational wave observations with the upcoming Laser Interferometer Space Antenna (LISA) offer the most promising way to test one of the leading progenitor scenarios - the double-degenerate scenario, which involves a binary system of two white dwarf stars. In this study we review published results, supplementing them with additional calculations for the context of Type Ia supernovae. We discuss the fact that LISA will be able to provide a complete sample of double white dwarf Type Ia supernova progenitors with orbital periods shorter than 16-11 minutes (gravitational wave frequencies above 2-3 millihertz). Such a sample will enable a statistical validation of the double-degenerate scenario by simply counting whether LISA detects enough double white dwarf binaries to account for the measured Type Ia merger rate in Milky Way-like galaxies. Additionally, we illustrate how LISA's capability to measure the chirp mass will set lower bounds on the primary mass, revealing whether detected double white dwarf binaries will eventually end up as a Type Ia supernova. We estimate that the expected LISA constraints on the Type Ia merger rate for the Milky Way will be 4-9%. We also discuss the potential gravitational wave signal from a Type Ia supernova assuming a double-detonation mechanism and explore how multi-messenger observations could significantly advance our understanding of these transient phenomena., Comment: accepted version

Recent observations have found a growing number of hypervelocity stars with speeds of $\approx 1500-2500\,$km\,s$^{-1}$ which could have only been produced through thermonuclear supernovae in white dwarf binaries. Most of the observed hypervelocity runaways in this class display a surprising inflated structure: their current radii are roughly an order of magnitude greater than they would have been as white dwarfs filling their Roche lobe. While many simulations exist studying the dynamical phase leading to supernova detonation in these systems, no detailed calculations of the long-term structure of the runaways have yet been performed. We use an existing \textsc{Arepo} hydrodynamical simulation of a supernova in a white dwarf binary as a starting point for the evolution of these stars with the 1 dimensional stellar evolution code MESA. We show that the supernova shock is not enough to inflate the white dwarf over timescales longer than a few thousand years, significantly shorter than the $10^{5-6}$ year lifetimes inferred for observed hypervelocity runaways. Despite experiencing a shock from a supernova less than $\approx 0.02\,R_\odot$ away, our models do not experience significant interior heating, and all contract back to radii around $0.01\,R_\odot$ within about $10^4$\,years. Explaining the observed inflated states requires either an additional source of significant heating or some other physics that is not yet accounted for in the subsequent evolution., Comment: Accepted for publication by A\&A. 15 pages, 17 figures

Stellar mergers are one important path to highly magnetised stars. Mergers of two low-mass white dwarfs may create up to every third hot subdwarf star. The merging process is usually assumed to dramatically amplify magnetic fields. However, so far only four highly magnetised hot subdwarf stars have been found, suggesting a fraction of less than $1\%$. We present two high-resolution magnetohydrodynamical (MHD) simulations of the merger of two helium white dwarfs in a binary system with the same total mass of $0.6\,M_\odot$. We analysed an equal-mass merger with two $0.3\,M_\odot$ white dwarfs, and an unequal-mass merger with white dwarfs of $0.25\,M_\odot$ and $0.35\,M_\odot$. We simulated the inspiral, merger, and further evolution of the merger remnant for about $50$ rotations. We found efficient magnetic field amplification in both mergers via a small-scale dynamo, reproducing previous results of stellar merger simulations. The magnetic field saturates at a similar strength for both simulations. We then identified a second phase of magnetic field amplification in both merger remnants that happens on a timescale of several tens of rotational periods of the merger remnant. This phase generates a large-scale ordered azimuthal field via a large-scale dynamo driven by the magneto-rotational instability. Finally, we speculate that in the unequal-mass merger remnant, helium burning will initially start in a shell around a cold core, rather than in the centre. This forms a convection zone that coincides with the region that contains most of the magnetic energy, and likely destroys the strong, ordered field. Ohmic resistivity might then quickly erase the remaining small-scale field. Therefore, the mass ratio of the initial merger could be the selecting factor that decides if a merger remnant will stay highly magnetised long after the merger., Comment: 13 pages, 10 figures, accepted by A&A, comments welcome

Simulations of galaxy formation are mostly unable to resolve the energy-conserving phase of individual supernova events, having to resort to subgrid models to distribute the energy and momentum resulting from stellar feedback. However, the properties of these simulated galaxies, including the morphology, stellar mass formed and the burstiness of the star formation history, are highly sensitive to numerical choices adopted in these subgrid models. Using the {\small SMUGGLE} stellar feedback model, we compute idealized simulations of a $M_{\rm vir} \sim 10^{10} \, \msun$ dwarf galaxy, a regime where most simulation codes predict significant burstiness in star formation, resulting in strong gas flows that lead to the formation of dark matter cores. We find that by varying only the directional distribution of momentum imparted from supernovae to the surrounding gas, while holding the total momentum per supernova constant, bursty star formation may be amplified or completely suppressed, and the total stellar mass formed can vary by as much as a factor of $\sim 3$. In particular, when momentum is primarily directed perpendicular to the gas disk, less bursty and lower overall star formation rates result, yielding less gas turbulence, more disky morphologies and a retention of cuspy dark matter density profiles. An improved understanding of the non-linear coupling of stellar feedback into inhomogeneous gaseous media is thus needed to make robust predictions for stellar morphologies and dark matter core formation in dwarfs independent of uncertain numerical choices in the baryonic treatment., Comment: Published ApJ; 18 pages, 12 figures; comments welcome

Galactic bars drive the internal evolution of spiral galaxies, while their formation is tightly coupled to the properties of their host galaxy and dark matter halo. To explore what drives bar formation in the cosmological context and how these structures evolve throughout cosmic history, we use the Auriga suite of magneto-hydrodynamical cosmological zoom-in simulations. We find that bars are robust and long-lived structures, and we recover a decreasing bar fraction with increasing redshift which plateaus around $\sim20\%$ at $z\sim3$. We find that bars which form at low and intermediate redshifts grow longer with time, while bars that form at high redshifts are born `saturated' in length, likely due to their merger-induced formation pathway. This leads to a larger bar-to-disc size ratio at high redshifts as compared to the local Universe. We subsequently examine the multi-dimensional parameter space thought to drive bar formation. We find that barred galaxies tend to have lower Toomre $Q$ values at the time of their formation, while we do not find a difference in the gas fraction of barred and unbarred populations when controlling for stellar mass. Barred galaxies tend to be more baryon-dominated at all redshifts, assembling their stellar mass earlier, while galaxies that are baryon-dominated but that do not host a bar, have a higher ex-situ bulge fraction. We explore the implications of the baryon-dominance of barred galaxies on the Tully-Fisher relation, finding an offset from the unbarred relation; confirming this in observations would serve as additional evidence for dark matter, as this behaviour is not readily explained in modified gravity scenarios., Comment: 21 pages including Appendices, 17 figures, submitted to MNRAS

Galactic outflows influence the evolution of galaxies not only by expelling gas from their disks but also by injecting energy into the circumgalactic medium (CGM). This alters or even prevents the inflow of fresh gas onto the disk and thus reduces the star formation rate. Supernovae (SNe) are the engines of galactic winds as they release thermal and kinetic energy into the interstellar medium (ISM). Cosmic rays (CRs) are accelerated at the shocks of SN remnants and only constitute a small fraction of the overall SN energy budget. However, their long live-times allow them to act far away from the original injection site and thereby to participate in the galactic wind launching process. Using high-resolution simulations of an isolated Milky Way-type galaxy with the moving-mesh code Arepo and the new multi-phase ISM model Crisp (Cosmic Rays and InterStellar Physics), we investigate how SNe and CRs launch galactic outflows and how the inclusion of CR-mediated feedback boosts the energy and mass entrained in the galactic wind. We find that the majority of thermal SN energy and momentum is used for stirring turbulence either directly or indirectly by causing fountain flows, thereby self-regulating the ISM and not for efficiently driving outflows to large heights. A simulation without CRs only launches a weak galactic outflow at uniformly high temperatures and low densities by means of the thermal pressure gradient. By contrast, most of the CR energy accelerated at SN remnants ($\sim80\%$) escapes the ISM and moves into the CGM. In the inner CGM, CRs dominate the overall pressure and are able to accelerate a large mass fraction in a galactic wind. This wind is turbulent and multi-phase with cold cloudlets embedded in dilute gas at intermediate temperatures ($\sim10^5$ K) and the CGM shows enhanced OVI and CIV absorption in comparison to a simulation without CRs., Comment: 16 pages, 9 figures, submitted to A&A, comments welcome

We investigate the growth of stellar discs in Milky Way-mass galaxies using the magnetohydrodynamical simulations of the Auriga Project in a full cosmological context. We focus on the gas accretion process along the discs, calculating the net, infall and outflow rates as a function of galactocentric distance, and investigate the relation between them and the star formation activity. The stellar distributions of around 70% of the simulated galaxies exhibit an ``inside-out'' pattern, with older (younger) stellar populations preferentially located in the inner (outer) disc regions. In all cases, we find a very tight correlation between the infall, outflow and net accretion rates, as well as between these three quantities and the star formation rate. This is because the amount of gas which is ultimately available for star formation in each radial ring depends not only on the infall rates, but also on the amount of gas leaving the disc in outflows, which directly relates to the local star formation level. Therefore, any of these rates can be used to identify galaxies with inside-out growth. For these galaxies, the correlation between the dominant times of accretion/star formation and disc radius is well fitted by a linear function. We also find that, when averaged over galaxies with formation histories similar to the Milky Way, the simulated accretion rates show a similar evolution (both temporally- and radially-integrated) to the usual accretion prescriptions used in chemical evolution models, although some major differences arise at early times and in the inner disc regions., Comment: 19 pages, 13 figures

Dwarf satellites in galaxy groups are distributed in an anisotropic and asymmetric manner, which is called the ``lopsided satellite distribution''. This lopsided signal has been observed not only in galaxy pairs but also in isolated systems. However, the physical origin of the lopsided signal in isolated systems is still unknown. In this work, we investigate this in the state-of-the-art hydrodynamical simulation of the MillenniumTNG Project by tracing each system back to high redshift. We find that the lopsided signal is dominated by satellites located in the outer regions of the halo and is also dominated by recently accreted satellites. The lopsided signal originates from the anisotropic accretion of galaxies from the surrounding large-scale structure and that, after accretion, the nonlinear evolution of satellites inside the dark-matter halo weakens the lopsidedness. The signal decreases as cosmic time passes because of a competition between anisotropic accretion and internal evolution within dark matter halos. Our findings provide a useful perspective for the study of galaxy evolution, especially for the origin of the spatial satellite galaxy distributions., Comment: 10 pages, 8 figures, 1 table. Submitted to MNRAS

We study the cosmological evolution of the two-point correlation functions of galactic gas-phase metal distributions using the 28 simulated galaxies from the Auriga Project. Using mock observations of the $z = 0$ snapshots to mimic our past work, we show that the correlation functions of the simulated mock observations are well matched to the correlation functions measured from local galaxy surveys. This comparison suggests that the simulations capture the processes important for determining metal correlation lengths, the key parameter in metallicity correlation functions. We investigate the evolution of metallicity correlations over cosmic time using the true simulation data, showing that individual galaxies undergo no significant systematic evolution in their metal correlation functions from $z\sim 3$ to today. In addition, the fluctuations in metal correlation length are correlated with but lag ahead fluctuations in star formation rate. This suggests that re-arrangement of metals within galaxies occurs at a higher cadence than star formation activity, and is more sensitive to the changes of environment, such as galaxy mergers, gas inflows / outflows, and fly-bys., Comment: 12 pages, 10 figures, 1 table, accepted for publication in MNRAS

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., Comment: 16 pages, 6 figures; v3: OJA published

We present an extended suite of the Auriga cosmological gravo-magnetohydrodynamical "zoom-in" simulations of 40 Milky Way-mass halos and 26 dwarf galaxy-mass halos run with the moving-mesh code Arepo. Auriga adopts the $\Lambda$ Cold Dark Matter ($\Lambda$CDM) cosmogony and includes a comprehensive galaxy formation physics model following the coupled cosmic evolution of dark matter, gas, stars, and supermassive black holes which has been shown to produce numerically well-converged galaxy properties for Milky Way-mass systems. We describe the first public data release of this augmented suite of Auriga simulations, which includes raw snapshots, group catalogues, merger trees, initial conditions, and supplementary data, as well as public analysis tools with worked examples of how to use the data. To demonstrate the value and robustness of the simulation predictions, we analyse a series of low-redshift global properties that compare well with many observed scaling relations, such as the Tully-Fisher relation, the star-forming main sequence, and HI gas fraction/disc thickness. Finally, we show that star-forming gas discs appear to build rotation and velocity dispersion rapidly for $z\gtrsim 3$ before they "settle" into ever-increasing rotation-dispersion ratios ($V/\sigma$). This evolution appears to be in rough agreement with some kinematic measurements from H$\alpha$ observations, and demonstrates an application of how to utilise the released data., Comment: Matches final MNRAS version

We investigate the numerical performance of a Discontinuous Galerkin (DG) hydrodynamics implementation when applied to the problem of driven, isothermal supersonic turbulence. While the high-order element-based spectral approach of DG is known to efficiently produce accurate results for smooth problems (exponential convergence with expansion order), physical discontinuities in solutions, like shocks, prove challenging and may significantly diminish DG's applicability to practical astrophysical applications. We consider whether DG is able to retain its accuracy and stability for highly supersonic turbulence, characterized by a network of shocks. We find that our new implementation, which regularizes shocks at sub-cell resolution with artificial viscosity, still performs well compared to standard second-order schemes for moderately high Mach number turbulence, provided we also employ an additional projection of the primitive variables onto the polynomial basis to regularize the extrapolated values at cell interfaces. However, the accuracy advantage of DG diminishes significantly in the highly supersonic regime. Nevertheless, in turbulence simulations with a wide dynamic range that start with supersonic Mach numbers and can resolve the sonic point, the low numerical dissipation of DG schemes still proves advantageous in the subsonic regime. Our results thus support the practical applicability of DG schemes for demanding astrophysical problems that involve strong shocks and turbulence, such as star formation in the interstellar medium. We also discuss the substantial computational cost of DG when going to high order, which needs to be weighted against the resulting accuracy gain. For problems containing shocks, this favours the use of comparatively low DG order., Comment: 22 pages, 14 figures, accepted for publication in MNRAS

The origin of thick disks and their evolutionary connection with thin disks are still a matter of debate. We provide new insights into this topic by connecting the stellar populations of thick disks at redshift $z=0$ with their past formation and growth, in 24 Milky Way-mass galaxies from the AURIGA zoom-in cosmological simulations. We projected each galaxy edge on, and decomposed it morphologically into two disk components, in order to define geometrically the thin and the thick disks as usually done in observations. We produced age, metallicity and [Mg/Fe] edge-on maps. We quantified the impact of satellite mergers by mapping the distribution of ex-situ stars. Thick disks are on average $\sim 3$~Gyr older, $\sim 0.25$~dex more metal poor and $\sim 0.06$~dex more [Mg/Fe]-enhanced than thin disks. Their average ages range from $\sim 6$ to $\sim 9$~Gyr, metallicities from $\sim -0.15$ to $\sim 0.1$~dex, and [Mg/Fe] from $\sim 0.12$ to $\sim 0.16$~dex. These properties are the result of an early initial in-situ formation, followed by a later growth driven by the combination of direct accretion of stars, some in-situ star formation fueled by mergers, and dynamical heating of stars. The balance between these processes varies from galaxy to galaxy. Mergers play a key role in the mass assembly of thick disks, contributing an average accreted mass fraction of $\sim 22$\% in the analyzed thick-disk dominated regions. In two galaxies, about half of the geometric thick-disk mass was directly accreted. While primordial thick disks form at high redshift in all galaxies, young metal-rich thin disks, with much lower [Mg/Fe] abundances, start to form later but at different times (higher or lower redshift) depending on the galaxy. We conclude that thick disks result from the interplay of external processes with the internal evolution of the galaxy., Comment: 57 pages, 43 figures, 3 tables; Published in Astronomy & Astrophysics

We present SPICE, a new suite of RHD cosmological simulations targeting the epoch of reionisation. The goal of these simulations is to systematically probe a variety of stellar feedback models, including "bursty" and "smooth" forms of supernova energy injection, as well as poorly-explored scenarios such as hypernova explosions and radiation pressure. Subtle differences in the behaviour of supernova feedback drive profound differences in reionisation histories, with burstier forms of feedback causing earlier reionisation. We also find that some global galaxy properties, such as the dust-attenuated luminosity functions and star formation main sequence, remain degenerate between models. Stellar feedback and its strength determine the morphological mix of galaxies emerging by z = 5 and that the reionisation history is inextricably connected to intrinsic properties such as galaxy kinematics and morphology. While star-forming, massive disks are prevalent if supernova feedback is "smooth", "bursty" feedback preferentially generates dispersion-dominated systems. Different modes of feedback produce different strengths of outflows, altering the ISM/CGM in different ways, and in turn strongly affecting the escape of LyC photons. We establish a correlation between galaxy morphology and LyC escape fraction, revealing that dispersion-dominated systems have escape fractions 10-50 times higher than their rotation-dominated counterparts at all redshifts. Dispersion-dominated systems should thus preferentially generate large HII regions as compared to their rotation-dominated counterparts. Since dispersion-dominated systems are more prevalent if stellar feedback is more explosive, reionisation occurs earlier in our simulation with burstier feedback. Statistical samples of post-reionisation galaxy morphologies probed with JWST, ALMA and MUSE can constrain stellar feedback and models of cosmic reionisation., Comment: arXiv admin note: text overlap with arXiv:1801.07259 by other authors

In the centers of dense star clusters, close encounters between stars and compact objects are likely to occur. We study tidal disruption events of main-sequence (MS) stars by stellar-mass black holes (termed $\mu$TDEs), which can shed light on the processes occurring in these clusters, including being an avenue in the mass growth of stellar-mass BHs. Using the moving-mesh hydrodynamics code \texttt{AREPO}, we perform a suite of hydrodynamics simulations of partial $\mu$TDEs of realistic, \texttt{MESA}-generated MS stars by varying the initial mass of the star ($0.5\,{\rm M}_{\rm \odot}$ and $1\,{\rm M}_{\rm \odot}$), the age of the star (zero-age, middle-age and terminal-age), the mass of the black hole ($10\,{\rm M}_{\rm \odot}$ and $40\,{\rm M}_{\rm \odot}$) and the impact parameter (yielding almost no mass loss to full disruption). We then examine the dependence of the masses, spins, and orbital parameters of the partially disrupted remnant on the initial encounter parameters. We find that the mass lost from a star decreases exponentially with increasing distance of approach and that a $1\,{\rm M}_{\rm \odot}$ star loses lesser mass than a $0.5\,{\rm M}_{\rm \odot}$. Moreover, a more evolved star is less susceptible to mass loss. Tidal torques at the closest approach spin up the remnant by factors of $10^2$--$10^4$ depending on the impact parameter. The remnant star can be bound (eccentric) or unbound (hyperbolic) to the black hole: hyperbolic orbits occur when the star's central density concentration is relatively low and the black hole-star mass ratio is high, which is the case for the disruption of a $0.5\,{\rm M}_{\rm \odot}$ star. Finally, we provide best-fit analytical formulae for a range of parameters that can be incorporated into cluster codes to model star-black hole interaction more accurately., Comment: Accepted for publication by A&A (with some revisions)

An important characteristic of cosmic hydrogen reionization is the growth of ionized gas bubbles surrounding early luminous objects. Ionized bubble sizes are beginning to be probed using Lyman-$\alpha$ emission from high-redshift galaxies, and will also be probed by upcoming 21-cm maps. We present results from a study of bubble sizes using the state-of-the-art THESAN radiation-hydrodynamics simulation suite, which self-consistently models radiation transport and realistic galaxy formation. We employ the mean-free path method, and track the evolution of the effective ionized bubble size at each point ($R_{\rm eff}$) throughout the Epoch of Reionization. We show there is a slow growth period for regions ionized early, but a rapid "flash ionization" process for regions ionized later as they immediately enter a large, pre-existing bubble. We also find that bright sources are preferentially in larger bubbles, and find consistency with recent observational constraints at $z \gtrsim 9$, but tension with idealized Lyman-$\alpha$ damping-wing models at $z \approx 7$. We find that high overdensity regions have larger characteristic bubble sizes, but the correlation decreases as reionization progresses, likely due to runaway formation of large percolated bubbles. Finally, we compare the redshift at which a region transitions from neutral to ionized ($z_{\rm reion}$) with the time it takes to reach a given bubble size and conclude that $z_{\rm reion}$ is a reasonable local probe of small-scale bubble size statistics ($R_\text{eff} \lesssim 1\,\rm{cMpc}$). However, for larger bubbles, the correspondence between $z_{\rm reion}$ and size statistics weakens due to the time delay between the onset of reionization and the expansion of large bubbles, particularly at high redshifts., Comment: 15 pages, 15 figures. Comments welcome. Please visit https://www.thesan-project.com for more details. Accepted for publication in MNRAS

Magnetic fields are ubiquitous in the Universe. Recently, cosmological simulations of galaxies have successfully begun to incorporate magnetic fields and their evolution in galaxies and their haloes. However, so far they have mostly focused on Milky Way-like galaxies. Here we analyse a sample of high resolution cosmological zoom simulations of disc galaxies in haloes with mass $M_\mathrm{200c}$ from $10^{10}\,\mathrm{M}_\odot$ to $10^{13}\,\mathrm{M}_\odot$, simulated with the Auriga galaxy formation model. We show that with sufficient numerical resolution the magnetic field amplification and saturation is converged. The magnetic field strength reaches equipartition with turbulent energy density for galaxies in haloes with $M_\mathrm{200c}\gtrsim 10^{11.5}\,\mathrm{M_\odot}$. For galaxies in less massive haloes, the magnetic field strength saturates at a fraction of equipartition that decreases with decreasing halo mass. For our lowest mass haloes, the magnetic field saturates significantly below $10\%$ of equipartition. We quantify the resolution we need to obtain converged magnetic field strengths and discuss our resolution requirements also in the context of the IllustrisTNG cosmological box simulations. We show that, at $z=0$, rotation-dominated galaxies in our sample exhibit for the most part an ordered large scale magnetic field, with fewer field reversals in more massive galaxies. Finally, we compare the magnetic fields in our cosmological galaxies at $z=0$ with simulations of isolated galaxies in a collapsing halo setup. Our results pave the way for detailed studies of cosmic rays and other physical processes in similar cosmological galaxy simulations that crucially depend on the strength and structure of magnetic fields., Comment: 19 pages, 12 figures, accepted for publication by MNRAS, comments welcome

We present a study of the evolution of cosmic filaments across redshift with an emphasis on some important properties: filament lengths, growth rates, and radial profiles of galaxy densities. Following an observation-driven approach, we build cosmic filament catalogues at z=0,1,2,3, and 4 from the galaxy distributions of the large hydro-dynamical run of the MilleniumTNG project. We employ the extensively used DisPerSE cosmic web finder code, for which we provide a user-friendly guide, including the details of a physics-driven calibration procedure, with the hope of helping future users. We perform the first statistical measurements of the evolution of connectivity in a large-scale simulation, finding that the connectivity of cosmic nodes (defined as the number of filaments attached) globally decreases from early to late times. The study of cosmic filaments in proper coordinates reveals that filaments grow in length and radial extent, as expected from large-scale structures in an expanding Universe. But the most interesting results arise once the Hubble flow is factored out. We find remarkably stable comoving filament length functions and over-density profiles, showing only little evolution of the total population of filaments in the past ~12.25 Gyrs. However, by tracking the spatial evolution of individual structures, we demonstrate that filaments of different lengths actually follow different evolutionary paths. While short filaments preferentially contract, long filaments expand along their longitudinal direction with growth rates that are the highest in the early, matter-dominated Universe. Filament diversity at fixed redshift is also shown by the different (~$5 \sigma$) density values between the shortest and longest filaments. Our results hint that cosmic filaments can be used as additional probes for dark energy, but further theoretical work is still needed., Comment: 17 pages, accepted in Astronomy & Astrophysics

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., Comment: Data and documentation at https://www.thesan-project.com, comments and requests welcome, paper accepted for publication in MNRAS

The physical collisions of two white dwarfs (WDs) (i.e. not slow mergers) have been shown to produce type-Ia-like supernovae (SNe) explosions. Most studies of WD-collisions have focused on zero impact-parameter (direct) collisions, which can also be studied in 2D. However, the vast majority of WD collisions arising from any evolutionary channels suggested to date, are expected to be indirect, i.e. have a non-negligible impact parameter upon collision. Here we use the highest resolution 3D simulations to date (making use of the AREPO code), in order to explore both direct and indirect collisions and the conditions in which they give rise to a detonation and the production of a luminous SNe. Using our simulations, we find a detonation criterion that can provide the critical impact parameter for an explosion to occur, depending on the density profile of the colliding WDs, their composition, and their collision velocities. We find that the initial velocity has a significant impact on the amount of 56Ni production from the explosion. Furthermore, the 56Ni production is also strongly dependent on the numerical resolution. While our results from the head-on collision with large initial velocities produce more 56Ni than previous simulations, those with small and comparable velocities produce significantly less 56Ni than previous works., Comment: Comments are welcome

Nebular-phase observations of peculiar Type Ia supernovae (SNe Ia) provide important constraints on progenitor scenarios and explosion dynamics for both these rare SNe and the more common, cosmologically useful SNe Ia. We present observations from an extensive ground-based and space-based follow-up campaign to characterize SN 2022pul, a "super-Chandrasekhar" mass SN Ia (alternatively "03fg-like" SN), from before peak brightness to well into the nebular phase across optical to mid-infrared (MIR) wavelengths. The early rise of the light curve is atypical, exhibiting two distinct components, consistent with SN Ia ejecta interacting with dense carbon-oxygen rich circumstellar material (CSM). In the optical, SN 2022pul is most similar to SN 2012dn, having a low estimated peak luminosity ($M_{B}=-18.9$ mag) and high photospheric velocity relative to other 03fg-like SNe. In the nebular phase, SN 2022pul adds to the increasing diversity of the 03fg-like subclass. From 168 to 336 days after peak $B$-band brightness, SN 2022pul exhibits asymmetric and narrow emission from [O I] $\lambda\lambda 6300,\ 6364$ (${\rm FWHM} \approx 2{,}000$ km s$^{-1}$), strong, broad emission from [Ca II] $\lambda\lambda 7291,\ 7323$ (${\rm FWHM} \approx 7{,}300$ km s$^{-1}$), and a rapid Fe III to Fe II ionization change. Finally, we present the first-ever optical-to-mid-infrared (MIR) nebular spectrum of an 03fg-like SN Ia using data from JWST. In the MIR, strong lines of neon and argon, weak emission from stable nickel, and strong thermal dust emission (with $T \approx 500$ K), combined with prominent [O I] in the optical, suggest that SN 2022pul was produced by a white dwarf merger within carbon/oxygen-rich CSM., Comment: 23 pages, 11 figures, submitted to ApJ

We present an analysis of ground-based and JWST observations of SN~2022pul, a peculiar "03fg-like" (or "super-Chandrasekhar") Type Ia supernova (SN Ia), in the nebular phase at 338d post explosion. Our combined spectrum continuously covers 0.4--14 $\mu$m and includes the first mid-infrared spectrum of an 03fg-like SN Ia. Compared to normal SN Ia 2021aefx, SN 2022pul exhibits a lower mean ionization state, asymmetric emission-line profiles, stronger emission from the intermediate-mass elements (IMEs) argon and calcium, weaker emission from iron-group elements (IGEs), and the first unambiguous detection of neon in a SN Ia. Strong, broad, centrally peaked [Ne II] line at 12.81 $\mu$m was previously predicted as a hallmark of "violent merger'' SN Ia models, where dynamical interaction between two sub-$M_{ch}$ white dwarfs (WDs) causes disruption of the lower mass WD and detonation of the other. The violent merger scenario was already a leading hypothesis for 03fg-like SNe Ia; in SN 2022pul it can explain the large-scale ejecta asymmetries seen between the IMEs and IGEs and the central location of narrow oxygen and broad neon. We modify extant models to add clumping of the ejecta to better reproduce the optical iron emission, and add mass in the innermost region ($< 2000$ km s$^{-1}$) to account for the observed narrow [O I]~$\lambda\lambda6300$, 6364 emission. A violent WD-WD merger explains many of the observations of SN 2022pul, and our results favor this model interpretation for the subclass of 03fg-like SN Ia., Comment: 20 pages, 10 figures, published in ApJ

Our current knowledge of star-forming metallicity relies primarily on gas-phase oxygen abundance measurements. This may not allow one to accurately describe differences in stellar evolution and feedback driven by variations in iron abundance. $\alpha$-elements (such as oxygen) and iron are produced by sources that operate on different timescales and the link between them is not straightforward. We explore the origin of the [O/Fe] - specific SFR (sSFR) relation, linking chemical abundances to galaxy formation timescales. This relation is followed by star-forming galaxies across redshifts according to cosmological simulations and basic theoretical expectations. Its apparent universality makes it suitable for trading the readily available oxygen for iron abundance. The relation is determined by the relative iron production efficiency of core-collapse and type Ia supernovae and the delay time distribution of the latter -- uncertain factors that could be constrained empirically with the [O/Fe]-sSFR relation. We compile and homogenise a literature sample of star-forming galaxies with observational iron abundance determinations to place first constraints on the [O/Fe]-sSFR relation over a wide range of sSFR. The relation shows a clear evolution towards lower [O/Fe] with decreasing sSFR and a flattening above log(sSFR/yr)>-9. The result is broadly consistent with expectations, but better constraints are needed to inform the models. We independently derive the relation from old Milky Way stars and find a remarkable agreement between the two, as long as the recombination-line absolute oxygen abundance scale is used in conjunction with stellar metallicity measurements., Comment: submitted to A&A

Understanding the evolution of galaxies provides crucial insights into a broad range of aspects in astrophysics, including structure formation and growth, the nature of dark energy and dark matter, baryonic physics, and more. It is, however, infeasible to track the evolutionary processes of individual galaxies in real time given their long timescales. As a result, galaxy evolution analyses have been mostly based on ensembles of galaxies that are supposed to be from the same population according to usually basic and crude observational criteria. We propose a new strategy of evaluating the evolution of an individual galaxy by identifying its descendant galaxies as guided by cosmological simulations. As a proof of concept, we examined the evolution of the total mass distribution of a target strong lensing galaxy at $z=0.884$ using the proposed strategy. We selected 158 galaxies from the IllustrisTNG300 simulation that we identified as analogs of the target galaxy. We followed their descendants and found 11 observed strong lensing galaxies that match in stellar mass and size with the descendants at their redshifts. The observed and simulated results are discussed, although no conclusive assessment is made given the low statistical significance due to the small sample size. Nevertheless, the test confirms that our proposed strategy is already feasible with existing data and simulations. We expect it to play an even more important role in studying galaxy evolution as more strong lens systems and larger simulations become available with the advent of next-generation survey programs and cosmological simulations., Comment: 11 pages

The double detonation is a widely discussed mechanism to explain Type Ia supernovae from explosions of sub-Chandrasekhar mass white dwarfs. In this scenario, a helium detonation is ignited in a surface helium shell on a carbon/oxygen white dwarf, which leads to a secondary carbon detonation. Explosion simulations predict high abundances of unburnt helium in the ejecta, however, radiative transfer simulations have not been able to fully address whether helium spectral features would form. This is because helium can not be sufficiently excited to form spectral features by thermal processes, but can be excited by collisions with non-thermal electrons, which most studies have neglected. We carry out a full non-local thermodynamic equilibrium (non-LTE) radiative transfer simulation for an instance of a double detonation explosion model, and include a non-thermal treatment of fast electrons. We find a clear He I {\lambda} 10830 feature which is strongest in the first few days after explosion and becomes weaker with time. Initially this feature is blended with the Mg II {\lambda} 10927 feature but over time separates to form a secondary feature to the blue wing of the Mg II {\lambda} 10927 feature. We compare our simulation to observations of iPTF13ebh, which showed a similar feature to the blue wing of the Mg II {\lambda} 10927 feature, previously identified as C I. Our simulation shows a good match to the evolution of this feature and we identify it as high velocity He I {\lambda} 10830. This suggests that He I {\lambda} 10830 could be a signature of the double detonation scenario., Comment: 7 pages, accepted by MNRAS

Magnetic fields are of critical importance for our understanding of the origin and long-term evolution of the Milky Way. This is due to their decisive role in the dynamical evolution of the interstellar medium (ISM) and their influence on the star-formation process. Faraday rotation measures (RM) along many different sightlines across the Galaxy are a primary means to infer the magnetic field topology and strength from observations. However, the interpretation of the data has been hampered by the failure of previous attempts to explain the observations in theoretical models and to synthesize a realistic multi-scale all-sky RM map. We here utilize a cosmological magnetohydrodynamic (MHD) simulation of the formation of the Milky Way, augment it with a novel star cluster population synthesis model for a more realistic structure of the local interstellar medium, and perform detailed polarized radiative transfer calculations on the resulting model. This yields a faithful first principles prediction of the Faraday sky as observed on Earth. The results reproduce the observations of the Galaxy not only on global scales, but also on local scales of individual star-forming clouds. They also imply that the Local Bubble containing our Sun dominates the RM signal over large regions of the sky. Modern cosmological MHD simulations of the Milky Way's formation, combined with a simple and plausible model for the fraction of free electrons in the ISM, explain the RM observations remarkably well, thus indicating the emergence of a firm theoretical understanding of the genesis of magnetic fields in our Universe across cosmic time., Comment: 14 pages, 3 figures

Faraday rotation measure (RM) is arguably the most practical observational tracer of magnetic fields in the diffuse circumgalactic medium (CGM). We sample synthetic Faraday rotation skies of Milky Way-like galaxies in TNG50 of the IllustrisTNG project by placing an observer inside the galaxies at a solar circle-like position. Our synthetic RM grids emulate specifications of current and upcoming surveys; the NRAO VLA Sky Survey (NVSS), the Polarisation Sky Survey of the Universe's Magnetism (POSSUM), and a future Square Kilometre Array (SKA1-mid) polarisation survey. It has been suggested that magnetic fields regulate the survival of high-velocity clouds. However, there is only a small number of observational detections of magnetised clouds thus far. In the first part of the paper, we test conditions for the detection of magnetised circumgalactic clouds. Based on the synthetic RM samplings of clouds in the simulations, we predict upcoming polarimetric surveys will open opportunities for the detection of even low-mass and distant clouds. In the second part of the paper, we investigate the imprint of the CGM in the all-sky RM distribution. We test whether the RM variation produced by the CGM is correlated with global galaxy properties, such as distance to a satellite, specific star formation rate, neutral hydrogen covering fraction, and accretion rate to the supermassive black hole. We argue that the observed fluctuation in the RM measurements on scales less than 1 degree, which has been considered an indication of intergalactic magnetic fields, might in fact incorporate a significant contribution of the Milky Way CGM., Comment: 18 pages, 11 figures, Accepted to MNRAS

Multi-body dynamical interactions of binaries with other objects are one of the main driving mechanisms for the evolution of star clusters. It is thus important to bring our understanding of three-body interactions beyond the commonly employed point-particle approximation. To this end we here investigate the hydrodynamics of three-body encounters between star-black hole (BH) binaries and single stars, focusing on the identification of final outcomes and their long-term evolution and observational properties, using the moving-mesh hydrodynamics code AREPO. This type of encounters produces five types of outcomes: stellar disruption, stellar collision, weak perturbation of the original binary, binary member exchange, and triple formation. The two decisive parameters are the binary phase angle, which determines which two objects meet at the first closest approach, and the impact parameter, which sets the boundary between violent and non-violent interactions. When the impact parameter is smaller than the semimajor axis of the binary, tidal disruptions and star-BH collisions frequently occur when the BH and the incoming star first meet, while the two stars mostly merge when the two stars meet first instead. In both cases, the BHs accrete from an accretion disk at super-Eddington rates, possibly generating flares luminous enough to be observed. The stellar collision products either form a binary with the BH or remain unbound to the BH. Upon collision, the merged stars are hotter and larger than main sequence stars of the same mass at similar age. Even after recovering their thermal equilibrium state, stellar collision products, if isolated, would remain hotter and brighter than main sequence stars until becoming giants., Comment: 16 pages, 10 figures, 2 tables. Submitted to MNRAS. Movie here: https://studio.youtube.com/playlist/PLxLK3qI02cQf_x96DSKQWVS251Jog7IOE/videos

We present a series of high-resolution simulations generated with the moving-mesh code AREPO to model the merger of a $1.1 \, \mathrm{M_\odot}$ carbon-oxygen primary white dwarf with an outer helium layer and a $0.35\,\mathrm{M_\odot}$ secondary helium white dwarf. Our simulations lead to detonations that are consistent with the edge-lit scenario, where a helium detonation is ignited at the base of the helium layer of the primary WD, which triggers an off-centre carbon detonation. This produces an asymmetric ejecta pattern and differences in line-of-sight observables (e.g. mean atomic weight). The ejecta that are flung into space are dominated by $^{56}\mathrm{Ni}$, $^{4}\mathrm{He}$, $^{28}\mathrm{Si}$, and $^{32}\mathrm{S}$. Our simulations result in a surviving degenerate companion of mass $0.22-0.25$ $\mathrm{M_\odot}$ moving at $>1\,700$ $\mathrm{km}\,\mathrm{s}^{-1}$, consistent with the observational findings of hypervelocity WDs. The secondary's surface layers are enriched by heavy metals, with $^{56}\mathrm{Ni}$ making up approximately $0.8 \%$ of the remaining mass. We also analyse the sensitivity of the outcome on simulation parameters, including the "inspiral time", which defines a period of accelerated angular momentum loss. We find that the choice of "inspiral time" qualitatively influences the simulation result, including the survival of the secondary. We argue that the shorter inspiral cases result in qualitatively and quantitatively similar outcomes. We also investigate the sensitivity of our results on the primary's chemical profile by comparing simulations using isothermal, constant composition models with the same mass and central composition and characterised by either a bare carbon-oxygen core (no helium) or a carbon-oxygen core enveloped by a thick helium layer., Comment: The submission has been accepted for publication in the Monthly Notices of the Royal Astronomical Society (MNRAS). 20 pages, 12 figures

In the coming years, Sunyaev-Zel'dovich (SZ) measurements can dramatically improve our understanding of the Intergalactic Medium (IGM) and the role of feedback processes on galaxy formation, allowing us to calibrate important astrophysical systematics in cosmological constraints from weak lensing galaxy clustering surveys. However, the signal is only measured in a two-dimensional projection, and its correct interpretation relies on understanding the connection between observable quantities and the underlying intrinsic properties of the gas, in addition to the relation between the gas and the underlying matter distribution. One way to address these challenges is through the use of hydrodynamical simulations such as the high-resolution, large-volume MillenniumTNG suite. We find that measurements of the optical depth, $\tau$, and the Compton-y parameter, $Y$, receive large line-of-sight contributions which can be removed effectively by applying a Compensated Aperture Photometry (CAP) filter. In contrast with other $\tau$ probes (e.g., X-rays and Fast Radio Bursts), the kSZ-inferred $\tau$ receives most of its signal from a confined cylindrical region around the halo due to the velocity decorrelation along the line-of-sight. Additionally, we perform fits to the $Y-M$ and $\tau-M$ scaling relations and report best-fit parameters adopting the smoothly broken power law (SBPL) formalism. We note that subgrid physics modeling can broaden the error bar on these by 30\% for intermediate-mass halos ($\sim$$10^{13} \, {\rm M}_{\odot}$). The scatter of the scaling relations can be captured by an intrinsic dependence on concentration, and an extrinsic dependence on tidal shear. Finally, we comment on the effect of using galaxies rather than halos in real observations, which can bias the inferred SZ profiles by $\sim$20\% for $L_\ast$-galaxies., Comment: 14 pages, 6 figures

The intrinsic alignment (IA) of observed galaxy shapes with the underlying cosmic web is a source of contamination in weak lensing surveys. Sensitive methods to identify the IA signal will therefore need to be included in the upcoming weak lensing analysis pipelines. Hydrodynamical cosmological simulations allow us to directly measure the intrinsic ellipticities of galaxies and thus provide a powerful approach to predict and understand the IA signal. Here we employ the novel, large-volume hydrodynamical simulation MTNG740, a product of the MillenniumTNG (MTNG) project, to study the IA of galaxies. We measure the projected correlation functions between the intrinsic shape/shear of galaxies and various tracers of large-scale structure, $w_{+g},\ w_{+m},\ w_{++}$ over the radial range $r_{\rm p} \in [0.02 , 200]\,h^{-1}{\rm Mpc}$ and at redshifts $z=0.0$, $0.5$ and $1.0$. We detect significant signal-to-noise IA signals with the density field for both elliptical and spiral galaxies. We also find significant intrinsic shear-shear correlations for ellipticals. We further examine correlations of the intrinsic shape of galaxies with the local tidal field. Here we find a significant IA signal for elliptical galaxies assuming a linear model. We also detect a weak IA signal for spiral galaxies under a quadratic tidal torquing model. Lastly, we measure the alignment between central galaxies and their host dark-matter halos, finding small to moderate misalignments between their principal axes that decline with halo mass., Comment: 16 pages, 14 figures

We study weak gravitational lensing convergence maps produced from the MillenniumTNG (MTNG) simulations by direct projection of the mass distribution on the past backwards lightcone of a fiducial observer. We explore the lensing maps over a large dynamic range in simulation mass and angular resolution, allowing us to establish a clear assessment of numerical convergence. By comparing full physics hydrodynamical simulations with corresponding dark-matter-only runs we quantify the impact of baryonic physics on the most important weak lensing statistics. Likewise, we predict the impact of massive neutrinos reliably far into the non-linear regime. We also demonstrate that the "fixed & paired" variance suppression technique increases the statistical robustness of the simulation predictions on large scales not only for time slices but also for continuously output lightcone data. We find that both baryonic and neutrino effects substantially impact weak lensing shear measurements, with the latter dominating over the former on large angular scales. Thus, both effects must explicitly be included to obtain sufficiently accurate predictions for stage IV lensing surveys. Reassuringly, our results agree accurately with other simulation results where available, supporting the promise of simulation modelling for precision cosmology far into the non-linear regime., Comment: 15 pages, 10 figures, comments welcome

We are witnessing the dawn of gravitational wave (GW) astronomy. With currently available detectors, observations are restricted to GW frequencies in the range between ${\sim} 10\,\mathrm{Hz}$ and $10\,\mathrm{kHz}$, which covers the signals from mergers of compact objects. The launch of the space observatory LISA will open up a new frequency band for the detection of stellar interactions at lower frequencies. In this work, we predict the shape and strength of the GW signals associated with common-envelope interaction and merger events in binary stars, and we discuss their detectability. Previous studies estimated these characteristics based on semi-analytical models. In contrast, we used detailed three-dimensional magnetohydrodynamic simulations to compute the GW signals. We show that for the studied models, the dynamical phase of common-envelope events and mergers between main-sequence stars lies outside of the detectability band of the LISA mission. We find, however, that the final stages of common-envelope interactions leading to mergers of the stellar cores fall into the frequency band in which the sensitivity of LISA peaks, making them promising candidates for detection. These detections can constrain the enigmatic common-envelope dynamics. Furthermore, future decihertz observatories such as DECIGO or BBO would also be able to observe this final stage and the post-merger signal, through which we might be able to detect the formation of Thorne-\.Zytkow objects., Comment: Accepted for publication in A&A, 12 pages, 8 figures

Cosmic rays (CRs) are dynamically important for the formation and evolution of galaxies by regulating star formation and by powering galactic outflows. However, to what extent CRs regulate galaxy formation depends on the coupling strength of CRs with the ambient plasma and the effective CR transport speed along the magnetic field. Moreover, both properties sensitively depend on the CR momentum, which is largely unexplored in three-dimensional hydrodynamical simulations. We perform magneto-hydrodynamical simulations of entire galaxies with masses ranging from $10^{10}$ to $10^{12}\,\mathrm{M}_\odot$ and compare dynamically coupled CRs in the grey approximation with a spectrally resolved model that includes CR momenta from $0.1\,\mathrm{GeV}~c^{-1}$ to $100\,\mathrm{TeV}~c^{-1}$. We find that hadronic cooling of CRs dominates over Alfv\'{e}n cooling, with the latter emulating CR losses as a result of streaming of CRs down their pressure gradient. While star formation rates and galaxy morphologies are only mildly affected by the spectral CR modelling, mass loading factors of galactic outflows can differ by up to a factor of four in dwarf galaxies. All simulated low-mass halos ($M=10^{10}$, $10^{11}$, and $3\times10^{11}\,\mathrm{M}_\odot$) drive strong outflows, where CR transport is temporally dominated by advection. In contrast, the Milky Way-mass galaxy with $M=10^{12}\,\mathrm{M}_\odot$ does not drive sustained outflows, so that CR transport is entirely dominated by diffusion. The effective energy weighted diffusion coefficients vary by two orders of magnitude from the canonical energy-weighted values of $\langle{D}\rangle_{e_\mathrm{cr}}\sim10^{28}\,\mathrm{cm^2\,s^{-1}}$ in the disc up to $3\times10^{29}\,\mathrm{cm^2\,s^{-1}}$ in the circumgalactic medium, where we observe substantial temperature and CR pressure differences between our grey and spectral CR models., Comment: 28 pages, MNRAS, accepted