Search

We introduce a GPU-accelerated hybrid hydro/ N -body code (Enzo-N) designed to address the challenges of concurrently simulating star clusters and their parent galaxies. This task has been exceedingly challenging, primarily due to the considerable computational time required, which stems from the substantial scale difference between galaxies (∼0.1 Mpc) and star clusters (∼parsecs). Yet, this significant scale separation means that particles within star clusters perceive those outside the star cluster in a semistationary state. By leveraging this aspect, we integrate a direct N -body code (Nbody6++GPU) into a cosmological (magneto)hydrodynamic code (Enzo) through utilization of the semistationary background acceleration approximation. We solve the dynamics of particles within star clusters using a direct N -body solver with regularization for few-body interactions, while evolving particles outside—dark matter, gas, and stars—using a particle-mesh gravity solver and hydrodynamic methods. We demonstrate that Enzo-N successfully simulates the coevolution of star clusters and their parent galaxies, capturing phenomena such as core collapse of the star cluster and tidal stripping due to galactic tides. This comprehensive framework opens up new possibilities for studying the evolution of star clusters within galaxies, offering insights that were previously inaccessible. [ABSTRACT FROM AUTHOR]

The rising part of a tidal disruption event light curve provides unique insight into early emission and the onset of accretion. Various mechanisms are proposed to explain the pre-peak emission, including shocks from debris interaction and reprocessing of disk emission. We study the pre-peak emission and its influence on the gas circularization by a series of gray radiation hydrodynamic simulations with varying black hole mass. We find that, given a super-Eddington fallback rate of 10 M ̇ Edd , the stream–stream collision can occur multiple times and drive strong outflows of up to 9 M ̇ Edd . By dispersing gas to ≳100 r s , the outflow can delay gas circularization and leads to sub-Eddington accretion rates during the first few stream–stream collisions. The stream–stream collision shock and circularization shock can sustain a luminosity of ∼1044 erg s−1 for days. The luminosity is generally sub-Eddington and shows a weak correlation with accretion rate at early times. The outflow is optically thick, yielding a reprocessing layer with a size of ∼1014 cm and photospheric temperature of ∼4 × 104 K. [ABSTRACT FROM AUTHOR]

Satellite galaxies can be used to indicate the dynamic mass of galaxy groups and clusters. In this study, we apply the axisymmetric Jeans Anisotropic Multi-Gaussian Expansion (JAM) modeling to satellite galaxies in 28 galaxy clusters selected from the TNG300-1 simulation with halo masses of log 10 M 200 / M ⊙ > 14.3 . If using true bound satellites as tracers, the best constrained total mass within the half-mass radius of satellites, M (< r half), and the virial mass, M 200, have average biases of −0.01 and 0.03 dex, with average scatters of 0.11 dex and 0.15 dex. If selecting companions in the redshift space with a line-of-sight depth of 2000 km s−1, the biases are −0.06 and 0.01 dex, while the scatters are 0.12 and 0.18 dex for M (< r half) and M 200. By comparing the best-fitting and actual density profiles, we find that ∼29% of the best-fitting density profiles show very good agreement with the truth, ∼32% display over/underestimates at most of the radial range with biased M (< r half), and 39% show under/overestimates in central regions and over/underestimates in the outskirts, with good constraints on M (< r half); yet most of the best constraints are still consistent with the true profiles within 1 σ statistical uncertainties for the three circumstances. Using a mock DESI Bright Galaxy Survey catalog with the effect of fiber incompleteness, we find DESI fiber assignments and the choice of flux limits barely modify the velocity dispersion profiles and are thus unlikely to affect the dynamical modeling outcomes. Our results show that with current and future deep spectroscopic surveys, JAM can be a powerful tool to constrain the underlying density profiles of individual massive galaxy clusters. [ABSTRACT FROM AUTHOR]

We perform the numerical simulations of axisymmetric, relativistic, optically thin jets under the influence of the radiation field of an accretion disk. We show that starting from a very low injection velocity at the base, jets can be accelerated to relativistic terminal speeds when traveling through the radiation field. The jet gains momentum through the interaction with the radiation field. We use a relativistic equation of state for multispecies plasma, which self-consistently calculates the adiabatic index for the jet material. All the jet solutions obtained are transonic in nature. In addition to the acceleration of the jet to relativistic speeds, our results show that the radiation field also acts as a collimating agent. The jets remain well collimated under the effect of radiation pressure. We also show that if the jet starts with a rotational velocity, the radiation field will reduce the angular momentum of the jet beam. [ABSTRACT FROM AUTHOR]

The interaction between supermassive black hole (SMBH) feedback and the circumgalactic medium (CGM) continues to be an open question in galaxy evolution. In our study, we use smoothed particle hydrodynamics simulations to explore the impact of SMBH feedback on galactic metal retention and the motion of metals and gas into and through the CGM of L* galaxies. We examine 140 galaxies from the 25 Mpc cosmological volume Romulus25, with stellar masses between log(M */ M ⊙) = 9.5–11.5. We measure the fraction of metals remaining in the interstellar medium (ISM) and CGM of each galaxy and calculate the expected mass of each SMBH based on the M BH– σ relation (Kormendy & Ho 2013). The deviation of each SMBH from its expected mass, Δ M BH, is compared to the potential of its host via σ. We find that SMBHs with accreted mass above M BH– σ are more effective at removing metals from the ISM than undermassive SMBHs in star-forming galaxies. Overall, overmassive SMBHs suppress the total star formation of their host galaxies and more effectively move metals from the ISM into the CGM. However, we see little to no evacuation of gas from the CGM out of their halos, in contrast with other simulations. Finally, we predict that C iv column densities in the CGM of L* galaxies are unlikely to depend on host galaxy SMBH mass. Our results show that the scatter in the low-mass end of the M BH– σ relation may indicate how effective an SMBH is in the local redistribution of mass in its host galaxy. [ABSTRACT FROM AUTHOR]

Using cosmological hydrodynamical zoom-in simulations, we explore the properties of subhalos in Milky Way analogs that contain a subcomponent of atomic dark matter (ADM). ADM differs from cold dark matter (CDM) due to the presence of self-interactions that lead to energy dissipation, analogous to standard model baryons. This model can arise in 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 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 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. [ABSTRACT FROM AUTHOR]

Circumbinary accretion occurs throughout the universe, from the formation of stars and planets to the aftermath of major galactic mergers. We present an extensive investigation of circumbinary accretion disks, studying circular binaries with mass ratios (q ≡ M 2/ M 1) from 0.01 to 1 and at each mass ratio probing the effects of disk thickness and viscosity. We study disks with aspect ratios H / r ∈ {0.1, 0.05, 0.03} and vary both the magnitude and spatial dependence of viscosity. Although thin accretion disks have previously been found to promote rapid inspirals of equal-mass binaries, we find that gravitational torques become weaker at lower mass ratios and most binaries with 0.01 ≤ q ≤ 0.04 outspiral, which may delay the coalescence of black hole binaries formed from minor mergers and cause high-mass exoplanets to migrate outward. However, in a number of cases, the disks accreting onto binaries with mass ratios ∼0.07 fail to develop eccentric modes, leading to extremely rapid inspirals. Variability in black hole accretion correlates with disk eccentricity, and we observe variability above the ∼10% level even for mass ratios of 0.01. We demonstrate that the spatial dependence of the viscosity (e.g., α vs. constant ν) significantly affects the degree of preferential accretion onto the secondary, resolving discrepancies between previous studies. Colder circumbinary disks remain eccentric even at q ∼ 0.01 and sustain deep, asymmetric cavities. [ABSTRACT FROM AUTHOR]

We study the coupling of hydrodynamics and reactions in simulations of the double-detonation model for Type Ia supernovae. When assessing the convergence of simulations, the focus is usually on spatial resolution; however, the method of coupling the physics together as well as the tolerances used in integrating a reaction network also play an important role. In this paper, we explore how the choices made in both coupling and integrating the reaction portion of a simulation (operator/Strang splitting versus the simplified spectral deferred corrections method we introduced previously) influences the accuracy, efficiency, and nucleosynthesis of simulations of double detonations. We find no need to limit reaction rates or reduce the simulation time step to the reaction timescale. The entire simulation methodology used here is GPU-accelerated and made freely available as part of the Castro simulation code. [ABSTRACT FROM AUTHOR]

Long-duration gamma-ray bursts (LGRBs), thought to be produced during core-collapse supernovae, may have a prominent neutron component in the outflow material. If present, neutrons can change how photons scatter in the outflow by reducing its opacity, thereby allowing the photons to decouple sooner than if there were no neutrons present. Understanding the details of this process could therefore allow us to probe the central engine of LGRBs, which is otherwise hidden. Here, we present results of the photospheric emission from an LGRB jet, using a combination of relativistic hydrodynamic simulations and radiative transfer postprocessing using Monte Carlo radiation transfer code. We control the size of the neutron component in the jet material by varying the equilibrium electron fraction Y e , and we find that the presence of neutrons in the GRB fireball affects the Band parameters α and E 0, while the picture with the β parameter is less clear. In particular, the break energy E 0 is shifted to higher energies. Additionally, we find that increasing the size of the neutron component also increases the total radiated energy of the outflow across multiple viewing angles. Our results not only shed light on LGRBs but are also relevant to short-duration gamma-ray bursts associated with binary neutron star mergers due to the likelihood of a prominent neutron component in such systems. [ABSTRACT FROM AUTHOR]

The detailed structure of core-collapse supernova progenitors is crucial for studying supernova explosion engines and the corresponding multimessenger signals. In this paper, we investigate the influence of stellar rotation on binary systems consisting of a 30 M ⊙ donor star and a 20 M ⊙ accretor using the MESA stellar evolution code. We find that through mass transfer in binary systems, fast-rotating red- and blue-supergiant progenitors can be formed within a certain range of the initial orbital periods, although the correlation is not linear. We also find that even with the same initial mass ratio of the binary system, the resulting final masses of the collapsars, the iron core masses, the compactness parameters, and the final rotational rates can vary widely and are sensitive to the initial orbital periods. For instance, the progenitors with strong convection form a thinner Si shell and a wider O shell compared to those in single-star systems. In addition, we conduct 2D self-consistent core-collapse supernova simulations with neutrino transport for these rotating progenitors derived from binary stellar evolution. We find that the neutrino and gravitational-wave signatures of these binary progenitors could exhibit significant variations. Progenitors with larger compactness parameters produce more massive proto-neutron stars, have higher mass accretion rates, and emit brighter neutrino luminosity and louder gravitational emissions. Finally, we observe stellar-mass black hole formation in some of our failed exploding models. [ABSTRACT FROM AUTHOR]

We perform radiation hydrodynamics simulations to study the structure and evolution of a photoevaporating protoplanetary disk. Ultraviolet and X-ray radiation from the host star heats the disk surface, where H2 pumping also operates efficiently. We run a set of simulations in which we varied the number of dust grains or the dust-to-gas mass ratio, which determines the relative importance between photoelectric heating and H2 pumping. We show that H2 pumping and X-ray heating contribute more strongly to the mass loss of the disk when the dust-to-gas mass ratio is D ≤ 10 − 3 . The disk mass-loss rate decreases with a lower dust amount, but remains around 10−10−11 M ⊙yr−1. In these dust-deficient disks, H2 pumping enhances photoevaporation from the inner disk region and shapes the disk mass-loss profile. We thus argue that the late-stage disk evolution is affected by the ultraviolet H2 pumping effect. The mass-loss rates derived from our simulations can be used in the study of long-term disk evolution. [ABSTRACT FROM AUTHOR]

We present the first results from a high-resolution simulation with a focus on galactic wind driving for an isolated galaxy with a halo mass of ∼1011 M ⊙ (similar to the Large Magellanic Cloud) and a total gas mass of ∼6 × 108 M ⊙, resulting in ∼108 gas cells at ∼4 M ⊙ mass resolution. We adopt a resolved stellar feedback model with nonequilibrium cooling and heating, including photoelectric heating and photoionizing radiation, as well as supernovae, coupled to the second-order meshless finite-mass method for hydrodynamics. These features make this the largest resolved interstellar medium (ISM) galaxy model run to date. We find mean star formation rates around 0.05 M ⊙ yr−1 and evaluate typical time-averaged loading factors for mass (η M ∼ 1.0, in good agreement with recent observations) and energy (η E ∼ 0.01). The bulk of the mass of the wind is transported by the warm (T < 5 × 105 K) phase, while there is a similar amount of energy transported in the warm and the hot phases (T > 5 × 105 K). We find an average opening angle of 30° for the wind, decreasing with higher altitude above the midplane. The wind mass loading is decreasing (flat) for the warm (hot) phase as a function of the star formation surface rate density ΣSFR, while the energy loading shows inverted trends with ΣSFR, decreasing for the warm wind and increasing for the hot wind, although with very shallow slopes. These scalings are in good agreement with previous simulations of resolved wind driving in the multiphase ISM. [ABSTRACT FROM AUTHOR]

We investigate the utility of deep learning for modeling the clustering of particles that are aerodynamically coupled to turbulent fluids. Using a Lagrangian particle module within the Athena++ hydrodynamics code, we simulate the dynamics of particles in the Epstein drag regime within a periodic domain of isotropic forced hydrodynamic turbulence. This setup is an idealized model relevant to the collisional growth of micron- to millimeter-sized dust particles in early-stage planet formation. The simulation data are used to train a U-Net deep-learning model to predict gridded three-dimensional representations of the particle density and velocity fields, given as input the corresponding fluid fields. The trained model qualitatively captures the filamentary structure of clustered particles in a highly nonlinear regime. We assess model fidelity by calculating metrics of the density field (the radial distribution function) and of the velocity field (the relative velocity and the relative radial velocity between particles). Although trained only on the spatial fields, the model predicts these statistical quantities with errors that are typically <10%. Our results suggest that, given appropriately expanded training data, deep learning could complement direct numerical simulations in predicting particle clustering within turbulent flows. [ABSTRACT FROM AUTHOR]

We investigate how patchy reionization affects the star formation history (SFH) and stellar metallicity of ultrafaint dwarf galaxies (UFDs). Patchy reionization refers to varying ultraviolet background strengths depending on a galaxy's environment. Recent observations highlight the significance of this effect on UFDs, as UFDs can have different SFHs depending on their relative position with respect to their host halo during the period of reionization. However, most cosmological hydrodynamic simulations do not consider environmental factors such as patchy reionization, and the effect of reionization is typically applied homogeneously. Using a novel approach to implement patchy reionization, we show how SFHs of simulated UFDs can change. Our cosmological hydrodynamic zoom-in simulations focus on UFD analogs with M vir ∼ 109 M ⊙, M * ≲ 105 M ⊙ at z = 0. We find that patchy reionization can weaken the effect of reionization by 2 orders of magnitude up to z = 3, enabling late star formation in half of the simulated UFDs, with quenching times ∼460 Myr later than those with homogeneous reionization. We also show that halo merger and mass assembly can affect the SFHs of simulated UFDs, in addition to patchy reionization. The average stellar iron-to-hydrogen ratio, [Fe/H], of the simulated UFDs with patchy reionization increases by 0.22–0.42 dex. Finally, our findings suggest that patchy reionization could be responsible for the extended SFHs of Magellanic UFDs compared to non-Magellanic UFDs. [ABSTRACT FROM AUTHOR]

Direct numerical simulations of three-dimensional compressible magnetohydrodynamic (MHD) turbulence have been performed in order to study the relation between wave modes and coherent structures and the consequent energization of test particles. Moreover, the question of which is the main mechanism of this particle energization is rigorously discussed. In particular, using the same initial conditions, we analyzed the nonlinear and linear evolution of a turbulent state along with the case of randomized phases. Then, the behaviors of the linear and nonlinear simulations were compared through the study of the time evolution of particle kinetic energy and preferential concentration. Also, spatiotemporal spectra were used to identify the presence of wave modes and quantify the fraction of energy around the MHD modes in linear and nonlinear simulations. Finally, the variation of the correlation time of the external forcing is studied in detail along with the effect on the particle energization (and clustering) and the presence of wave modes. More specifically, particle energization tends to decrease when the fraction of linear energy increases, supporting the idea that energization by structures is the dominant mechanism for particle energization instead of resonance with wave modes as suggested by Fermi energization theory. [ABSTRACT FROM AUTHOR]

Streaming instability (SI) in dust has long been thought to be a promising process in triggering planetesimal formation in the protoplanetary disks (PPDs). In this study, we present the first numerical investigation that models the SI in the vertically stratified disk together with the dust coagulation process. Our simulations reveal that, even with the initially small dust sizes, because dust coagulation promotes dust size growth, SI can eventually still be triggered. Specifically, dust coagulation, limited only by dust fragmentation, broadens the parameter boundaries obtained from previous SI studies using single dust species. We describe the various stages of dust dynamics along with their size evolution and explore the impact of different dust fragmentation velocities. Implications of these results for realistic PPDs are also discussed.

Shadows are often observed in transition disks, which can result from obscuring by materials closer to the star, such as a misaligned inner disk. While shadows leave apparent darkened emission as observational signatures, they have significant dynamical impact on the disk. We carry out 3D radiation-hydrodynamical simulations to study shadows in transition disks and find that the temperature drop due to the shadow acts as an asymmetric driving force, leading to spirals in the cavity. These spirals have zero pattern speed following the fixed shadow. The pitch angle is given by tan ^−1 ( c _s / v _ϕ ) (6° if h / r = 0.1). These spirals transport mass through the cavity efficiently, with α ∼ 10 ^−2 in our simulation. Besides spirals, the cavity edge can also form vortices and flocculent streamers. When present, these features could disturb the shadow-induced spirals. By carrying out Monte Carlo radiative transfer simulations, we show that these features resemble those observed in near-infrared scattered light images. In the vertical direction, the vertical gravity is no longer balanced by the pressure gradient alone. Instead, an azimuthal convective acceleration term balances the gravity–pressure difference, leading to azimuthally periodic upward and downward gas motion reaching 10% of the sound speed, which can be probed by Atacama Large Millimeter/submillimeter Array line observations.

Recent observational constraints on the internal structure of Jupiter and Saturn suggest that these planets have “fuzzy” cores, i.e., gradients of the concentration of heavy elements that might span a large fraction of the planet’s radius. These cores could be composed of a semiconvective staircase, i.e., multiple convective layers separated by diffusive interfaces arising from double-diffusive instabilities. However, to date, no study has demonstrated how such staircases can avoid layer mergers and persist over evolutionary timescales. In fact, previous work has found that these mergers occur rapidly, leading to only a single convective layer. Using 3D simulations, we demonstrate that rotation prolongs the lifetime of a convective staircase by increasing the timescale for both layer merger and erosion of the interface between the final two layers. We present an analytic model for the erosion phase, predicting that rotation increases the erosion time by a factor of approximately Ro ^−1/2 , where Ro is the Rossby number of the convective flows (the ratio of the rotation period to the convective turnover time). For Jovian conditions at early times after formation (when convection is vigorous enough to mix a large fraction of the planet), we find the erosion time to be roughly 10 ^9 yr in the nonrotating case and 10 ^11 yr in the rotating case. If these timescales are confirmed with a larger suite of numerical simulations, the existence of convective staircases within the deep interiors of giant planets is a strong possibility, and rotation could be an important factor in the preservation of their fuzzy cores.

The gas surrounding first-generation (Population III, hereafter Pop III) stars is expected to emit a distinct signature in the form of the He ii recombination line at 1640 Å(He ii λ 1640). Here we explore the challenges and opportunities in identifying this elusive stellar population via the He ii λ 1640 in M _⋆ > 10 ^7.5 M _⊙ galaxies during the Epoch of Reionization ( z ≃ 6–10), using JWST/NIRSpec. With this aim in mind, we combine cosmological dustyGadget simulations with analytical modeling of the intrinsic He ii emission. While tentative candidates with bright He ii emission like GN-z11 have been proposed in the literature, the prevalence of such bright systems remains unclear due to significant uncertainties involved in the prediction of the He ii luminosity. In fact, similar Pop III clumps might be almost 2 orders of magnitude fainter, primarily depending on the assumed Pop III formation efficiency and initial mass function in star-forming clouds, while the effect of stellar mass loss is responsible for a factor of order unity. Moreover, up to ∼90% of these clumps might be missed with NIRSpec multi-object spectroscopy due to the limited field of view, while this problem appears to be less severe with NIRSpec's integral field unit. We investigate the potential of deep spectroscopy targeting peripheral Pop III clumps around bright, massive galaxies to achieve a clear detection of the first stars.

We present the first measurement of the Ly α forest baryon acoustic oscillations (BAO) shift parameter from cosmological simulations. In particular, we generate a suite of 1000 accurate effective field-level bias-based Ly α forest simulations of volume $V={(1\,{h}^{-1}\,\mathrm{Gpc})}^{3}$ at z = 2, both in real and redshift space, calibrated upon two fixed-and-paired cosmological hydrodynamic simulations. To measure the BAO, we stack the 3D power spectra of the 1000 different realizations, compute the average, and use a model accounting for a proper smooth-peak component decomposition of the power spectrum, to fit it via an efficient Markov Chain Monte Carlo scheme estimating the covariance matrices directly from the simulations. We report the BAO shift parameters to be $\alpha ={0.9969}_{-0.0014}^{+0.0014}$ and $\alpha ={0.9905}_{-0.0027}^{+0.0027}$ in real and redshift space, respectively. We also measure the bias b _lya and the BAO broadening parameter Σ _nl , finding ${b}_{\mathrm{lya}}=-{0.1786}_{-0.0001}^{+0.0001}$ and ${{\rm{\Sigma }}}_{\mathrm{nl}}={3.87}_{-0.20}^{+0.20}$ in real space, and ${b}_{\mathrm{lya}}=-{0.073}_{-0.004}^{+0.005}$ and ${{\rm{\Sigma }}}_{\mathrm{nl}}={6.55}_{-0.22}^{+0.23}$ in redshift space. Moreover, we measure the linear Kaiser factor ${\beta }_{\mathrm{lya}}={1.39}_{-0.18}^{+0.24}$ from the isotropic redshift space fit. Overall, we find evidence for a negative shift of the BAO peak at the ∼2.2 σ and ∼3.5 σ levels in real and redshift space, respectively. This work sets new important theoretical constraints on the Ly α forest BAO scale and offers a potential solution to the tension emerging from previous observational analysis, in light of ongoing and upcoming Ly α forest spectroscopic surveys, such as DESI, the Prime Focus Spectrograph Survey, and WEAVE-QSO.

Type Ia supernovae (SNe Ia) are standardizable cosmological candles that led to the discovery of the accelerating Universe. However, the physics of how white dwarfs (WDs) explode and lead to SNe Ia is still poorly understood. The initiation of the detonation front that rapidly disrupts the WD is a crucial element of the puzzle, and global 3D simulations of SNe Ia cannot resolve the requisite length scales to capture detonation initiation. In this work, we elucidate a theoretical criterion for detonation initiation in the distributed burning regime. We test this criterion against local 3D driven turbulent hydrodynamical simulations within electron-degenerate WD matter consisting initially of pure helium. We demonstrate a novel pathway for detonation, in which strong turbulent dissipation rapidly heats the helium, and forms carbon nuclei sufficient to lead to a detonation through accelerated burning via α captures. Simulations of strongly driven turbulent conditions lead to detonations at a mean density of 10 ^6 g cm ^−3 and mean temperature of 1.4–1.8 × 10 ^9 K, but fail to detonate at a lower density of 10 ^5 g cm ^−3 , in excellent agreement with theoretical predictions.

A kinematic misalignment of the stellar and gas components is a phenomenon observed in a significant fraction of galaxies. However, the underlying physical mechanisms are not well understood. A commonly proposed scenario for the formation of a misaligned component requires any preexisting gas disk to be removed, via flybys or ejective feedback from an active galactic nucleus. In this Letter, we study the evolution of a Milky Way mass galaxy in the FIREbox cosmological volume that displays a thin, counterrotating gas disk with respect to its stellar component at low redshift. In contrast to scenarios involving gas ejection, we find that preexisting gas is mainly removed via the conversion into stars in a central starburst, triggered by a merging satellite galaxy. The newly accreted, counterrotating gas eventually settles into a kinematically misaligned disk. About 4% (8 out of 182) of FIREbox galaxies with stellar masses larger than 5 × 10 ^9 M _⊙ at z = 0 exhibit gas–star kinematic misalignment. In all cases, we identify central starburst-driven depletion as the main reason for the removal of the preexisting corotating gas component, with no need for feedback from, e.g., a central active black hole. However, during the starburst, the gas is funneled toward the central regions, likely enhancing black hole activity. By comparing the fraction of misaligned discs between FIREbox and other simulations and observations, we conclude that this channel might have a non-negligible role in inducing kinematic misalignment in galaxies.

The James Webb Space Telescope (JWST) is capable of probing extremely early eras of our Universe, when the supersonic relative motions between dark matter and baryonic overdensities modulate structure formation ( z ≳ 10). We study low-mass galaxy formation, including this “stream velocity,” using high-resolution AREPO hydrodynamics simulations and present theoretical predictions of the UV luminosity function (UVLF) and galaxy stellar mass function down to extremely faint and low-mass galaxies ( M _UV ≳ −15, 10 ^4 M _⊙ ≤ M _* ≤ 10 ^8 M _⊙ ). We show that, although the stream velocity suppresses early star formation overall, it induces a short period of rapid star formation in some larger dwarfs, leading to an enhancement in the faint end of the UVLF at z = 12. We demonstrate that JWST observations are close to this enhanced regime and propose that the UVLF may constitute an important probe of the stream velocity at high redshift for JWST and future observatories.

Current models of galaxy formation require strong feedback from active galactic nuclei (AGN) to explain the observed lack of star formation in massive galaxies since z ≈ 2, but direct evidence of this energy input is limited. We use the SIMBA cosmological galaxy formation simulations to assess the ability of thermal Sunyaev–Zel'dovich (tSZ) measurements to provide such evidence, by mapping the pressure structure of the circumgalactic medium around massive z ≈ 0.2–1.5 galaxies. We undertake a stacking approach to calculate the total tSZ signal and its radial profile in simulations with varying assumptions of AGN feedback, and we assess its observability with current and future telescopes. By convolving our predictions with the 2.′1 beam of the Atacama Cosmology Telescope, we show that current observations at z ≈ 1 are consistent with SIMBA's fiducial treatment of AGN feedback and inconsistent with SIMBA models without feedback. At z ≈ 0.5, observational signals lie between SIMBA run with and without AGN feedback, suggesting AGN in SIMBA may inject too much energy at late times. By convolving our data with a 9.″5 beam corresponding to the TolTEC camera on the Large Millimeter Telescope Alfonso Serrano, we predict a unique profile for AGN feedback that can be distinguished with future higher-resolution measurements. Finally, we explore a novel approach to quantify the nonspherically symmetric features surrounding our galaxies by plotting radial profiles representing the component of the stack with m-fold symmetry. [ABSTRACT FROM AUTHOR]

We examine the quenching characteristics of 328 isolated dwarf galaxies 10 8 < M star / M ⊙ < 10 10 within the Romulus25 cosmological hydrodynamic simulation. Using mock-observation methods, we identify isolated dwarf galaxies with quenched star formation and make direct comparisons to the quenched fraction in the NASA Sloan Atlas (NSA). Similar to other cosmological simulations, we find a population of quenched, isolated dwarf galaxies below M star < 109 M ⊙ not detected within the NSA. We find that the presence of massive black holes (MBHs) in Romulus25 is largely responsible for the quenched, isolated dwarfs, while isolated dwarfs without an MBH are consistent with quiescent fractions observed in the field. Quenching occurs between z = 0.5–1, during which the available supply of star-forming gas is heated or evacuated by MBH feedback. Mergers or interactions seem to play little to no role in triggering the MBH feedback. At low stellar masses, M star ≲ 109.3 M ⊙, quenching proceeds across several Gyr as the MBH slowly heats up gas in the central regions. At higher stellar masses, M star ≳ 109.3 M ⊙, quenching occurs rapidly within 1 Gyr, with the MBH evacuating gas from the central few kpc of the galaxy and driving it to the outskirts of the halo. Our results indicate the possibility of substantial star formation suppression via MBH feedback within dwarf galaxies in the field. On the other hand, the apparent overquenching of dwarf galaxies due to MBH suggests that higher-resolution and/or better modeling is required for MBHs in dwarfs, and quenched fractions offer the opportunity to constrain current models. [ABSTRACT FROM AUTHOR]

We investigate the impact of binary orbital motions on the dynamical modeling of dwarf galaxies with intrinsic line-of-sight velocity dispersions ( σ v r ) of 1–9 km s−1. Using dwarf galaxies from the auriga level-2 and level-3 simulations, we apply the Jeans Anisotropic Multi-Gaussian Expansion modeling to tracer stars before and after including binaries to recover the dynamical masses. The recovered total masses within the half-mass radius of tracers, M (< r half), are always inflated due to binary motions, with greater inflations occurring for smaller σ v r . However, many dwarf galaxies experience central density deflated due to binary motions, with little dependence on σ v r . This is due to the negative radial gradients in the velocity dispersion profiles, with the fractional inflation in σ v r due to binaries more significant in outskirts. An extreme binary fraction of 70% can lead to central density deflation of up to 10%–20% at 3 km s−1 < σ v r < 8 km s−1, with M (< r half) inflated by 4% at 9 km s−1 and up to 15% at 3 km s−1. A lower binary fraction of 36% leads to similar deflations, with the inflations decreasing to approximately 10% at 3 km s−1 and becoming statistically insignificant. The choice of binary orbit distribution models does not result in significant differences, and observational errors tend to slightly weaken the deflations in the recovered central density. Two observations separated by 1 yr to exclude binaries lead to almost zero inflations/deflations for a binary fraction of 36% over 3 km s−1 < σ v r < 9 km s−1. For σ v r ∼ 1 km s−1 to 3 km s−1, a binary fraction of 70% (36%) still results in 60% (30%) to 10% (1%) of inflations in M (< r half), even with two-epoch observation. [ABSTRACT FROM AUTHOR]

The Milky Way (MW) stellar disk has both a thin and a thick component. The thin disk is composed mostly of younger stars (≲8 Gyr) with a lower abundance of α -elements, while the thick disk contains predominantly older stars (≳8–12 Gyr) with a higher α abundance, giving rise to an α -bimodality most prominent at intermediate metallicities. A proposed explanation for the bimodality is an episode of clumpy star formation, where high- α stars form in massive clumps that appear in the first few billion years of the MW's evolution, while low- α stars form throughout the disk and over a longer time span. To better understand the evolution of clumps, we track them and their constituent stars in two clumpy MW simulations that reproduce the α -abundance bimodality, one with 10% and the other with 20% supernova feedback efficiency. We investigate the paths that these clumps take in the chemical space ([O/Fe]–[Fe/H]) as well as their mass, star formation rate (SFR), formation location, lifetime, and merger history. The clumps in the simulation with lower feedback last longer on average, with several lasting hundreds of millions of years. Some of the clumps do not reach high- α, but the ones that do on average have a higher SFR, longer lifetime, greater mass, and form closer to the Galactic center than the ones that do not. Most clumps that reach high- α merge with others and eventually spiral into the Galactic center, but shed stars along the way to form most of the thick-disk component. [ABSTRACT FROM AUTHOR]

We train graph neural networks to perform field-level likelihood-free inference using galaxy catalogs from state-of-the-art hydrodynamic simulations of the CAMELS project. Our models are rotational, translational, and permutation invariant and do not impose any cut on scale. From galaxy catalogs that only contain 3D positions and radial velocities of ∼1000 galaxies in tiny (25 h − 1 Mpc) 3 volumes our models can infer the value of Ωm with approximately 12% precision. More importantly, by testing the models on galaxy catalogs from thousands of hydrodynamic simulations, each having a different efficiency of supernova and active galactic nucleus feedback, run with five different codes and subgrid models—IllustrisTNG, SIMBA, Astrid, Magneticum, SWIFT-EAGLE—we find that our models are robust to changes in astrophysics, subgrid physics, and subhalo/galaxy finder. Furthermore, we test our models on 1024 simulations that cover a vast region in parameter space—variations in five cosmological and 23 astrophysical parameters—finding that the model extrapolates really well. Our results indicate that the key to building a robust model is the use of both galaxy positions and velocities, suggesting that the network has likely learned an underlying physical relation that does not depend on galaxy formation and is valid on scales larger than ∼10 h −1 kpc. [ABSTRACT FROM AUTHOR]

We study the relation between the metallicities of ionized and atomic gas in star-forming galaxies at z = 0–3 using the Evolution and Assembly of GaLaxies and their Environments (EAGLE) cosmological, hydrodynamical simulations. This is done by constructing a dense grid of sight lines through the simulated galaxies and obtaining the star formation rate- and H i column density-weighted metallicities, Z SFR and Z H I, for each sightline as proxies for the metallicities of ionized and atomic gas, respectively. We find Z SFR ≳ Z H I for almost all sight lines, with their difference generally increasing with decreasing metallicity. The stellar masses of galaxies do not have a significant effect on this trend, but the positions of the sight lines with respect to the galaxy centers play an important role: the difference between the two metallicities decreases when moving toward the galaxy centers, and saturates to a minimum value in the central regions of galaxies, irrespective of redshift and stellar mass. This implies that the mixing of the two gas phases is most efficient in the central regions of galaxies where sight lines generally have high column densities of H i. However, a high H i column density alone does not guarantee a small difference between the two metallicities. In galaxy outskirts, the inefficiency of the mixing of star-forming gas with H i seems to dominate over the dilution of heavy elements in H i through mixing with the pristine gas. We find good agreement between the available observational data and the Z SFR– Z H I relation predicted by the EAGLE simulations. Though, observed regions with a nuclear starburst mode of star formation appear not to follow the same relation. [ABSTRACT FROM AUTHOR]

Studying compact-star binaries and their mergers is integral to determining progenitors for observable transients. Today, compact-star mergers are typically studied via state-of-the-art computational fluid dynamics codes. One such numerical technique, smoothed particle hydrodynamics (SPH), is frequently chosen for its excellent mass, energy, and momentum conservation. The natural treatment of vacuum and the ability to represent highly irregular morphologies make SPH an excellent tool for the study of compact-star binaries and mergers. For many scenarios, including binary systems, the outcome of simulations is only as accurate as the initial conditions. For SPH, it is essential to ensure that the particles are distributed regularly, representing the initial density profile but without long-range correlations. Particle noise in the form of high-frequency local motion and low-frequency global dynamics must be damped out. Damping the latter can be as computationally intensive as the actual simulation. We discuss a new and straightforward relaxation method, halted-pendulum relaxation (HPR), to remove global oscillation modes of SPH particle configurations. In combination with effective external potentials representing gravitational and orbital forces, we show that HPR has an excellent performance in efficiently relaxing SPH particles to the desired density distribution and removing global oscillation modes. We compare the method to frequently used relaxation approaches and test it on a white dwarf binary model at its Roche-lobe overflow limit. We highlight the importance of our method in achieving accurate initial conditions and its effect on achieving circular orbits and realistic accretion rates when compared with other general relaxation methods. [ABSTRACT FROM AUTHOR]

Precipitation is potentially a mechanism through which the circumgalactic medium (CGM) can regulate a galaxy's star formation. Here, we present idealized simulations of isolated Milky Way–like galaxies intended to examine the ability of galaxies to self-regulate their star formation, in particular via precipitation. We also examine the impact of rotation in the CGM. Using six simulations, we explore variations in the initial CGM t cool/ t ff ratio and rotation profile. Those variations affect the amount of gas accretion and star formation within the galactic disk. To encourage this accretion and better study its dependence on CGM structure, we gradually increase the efficiency of stellar feedback during the first half of our simulations. Yet despite this gradual increase, the resulting outflows quickly evacuate large, hot cavities within the CGM and even beyond r 200. Some of the CGM gas avoids interacting with the cavities and is able to feed the disk along its midplane, but the cooling of feedback-heated gas far from the midplane is too slow to supply the disk with additional gas. Our simulations illustrate the importance of physical mechanisms in the outer CGM and IGM for star formation regulation in Milky Way–scale halos. [ABSTRACT FROM AUTHOR]

One of the leading explanations for the origin of Fermi Bubbles is past jet activity in the Galactic center supermassive black hole Sgr A*. The claimed jets are often assumed to be perpendicular to the Galactic plane. Motivated by the orientation of pc-scale nuclear stellar disk and gas streams, as well as a low inclination of the accretion disk around Sgr A* inferred by the Event Horizon Telescope, we perform hydrodynamical simulations of nuclear jets significantly tilted relative to the Galactic rotation axis. The observed axisymmetry and hemisymmetry (north–south symmetry) of Fermi/eROSITA bubbles (FEBs) due to quasi-steady jets in Sgr A* could be produced if the jet had a super-Eddington power (≳5 × 1044 erg s−1) for a short time (jet active period ≲6 kyr) for a reasonable jet opening angle (≲10°). Such powerful explosions are, however, incompatible with the observed O viii /O vii line ratio toward the bubbles, even after considering electron–proton temperature nonequilibrium. We argue that the only remaining options for producing FEBs are (i) a low-luminosity (≈1040.5–41 erg s−1) magnetically dominated jet or accretion wind from the Sgr A*, or (ii) a supernovae or tidal disruption event driven wind of a similar luminosity from the Galactic center. [ABSTRACT FROM AUTHOR]

Observed breakBRD ("break bulges in red disks") galaxies are a nearby sample of face-on disk galaxies with particularly centrally concentrated star formation: they have red disks but recent star formation in their centers as measured by the D n 4000 spectral index. In Kopenhafer et al., a comparable population of breakBRD analogs was identified in the TNG simulation, in which the central concentration of star formation was found to reflect a central concentration of dense, star-forming gas caused by a lack of dense gas in the galaxy outskirts. In this paper, we examine the circumgalactic medium of the central breakBRD analogs to determine if the extended halo gas also shows differences from that around comparison galaxies with comparable stellar mass. We examine the circumgalactic medium gas mass, specific angular momentum, and metallicity in these galaxy populations. We find less gas in the circumgalactic medium of breakBRD galaxies, and that the breakBRD circumgalactic medium is slightly more concentrated than that of comparable M * galaxies. In addition, we find that the angular momentum in the circumgalactic medium of breakBRD galaxies tends to be low for their stellar mass, and shows more misalignment to the angular momentum vector of the stellar disk. Finally, we find that the circumgalactic medium metallicity of breakBRD galaxies tends to be high for their stellar mass. Together with their low star formation rate, we argue that these circumgalactic medium properties indicate a small amount of disk feeding concentrated in the central regions and a lack of low-metallicity gas accretion from the intergalactic medium. [ABSTRACT FROM AUTHOR]

We present the Cosmic Sands suite of cosmological zoom-in simulations based on the simba galaxy formation model in order to study the buildup of the first massive and dusty galaxies in the early universe. Residing in the most massive halos, we find that the compact proto-massive galaxies undergo nearly continuous mergers with smaller subhalos, boosting star formation rates (SFRs) and the buildup of stellar mass. The galaxies are already appreciably chemically evolved by z = 7, with modeled dust masses comparable to those inferred from observations in the same epoch, except for the most extreme systems. We track gas accretion onto the galaxies to understand how extreme SFRs can be sustained by these early systems. We find that smooth gas accretion can maintain SFRs above 250 M ⊙ yr−1, but to achieve SFRs that boost galaxies well above the main sequence, a larger perturbation like a gas-rich major merger is necessary to trigger a starburst episode. Post-processing the Cosmic Sands simulations with dust RT, we find that, while the infrared luminosities of the most-dust-rich galaxies are comparable to local ULIRGs, they are substantially dimmer than classical z = 2 submillimeter galaxies. We end with a discussion on the possible reasons for this discrepancy at the highest masses and the future work we intend to carry out to study the chemical enrichment of the earliest dusty galaxies. [ABSTRACT FROM AUTHOR]

We investigated dusty and dust-free gas dynamics for a radiation-driven sub-parsec-scale outflow in an active galactic nucleus (AGN) associated with a supermassive black hole 107 M ⊙ and bolometric luminosity 1044 erg s−1 based on the two-dimensional radiation-hydrodynamic simulations. A radiation-driven "lotus-like" multi-shell outflow is launched from the inner part (r ≲ 0.04 pc) of the geometrically thin disk, and it repeatedly and steadily produces shocks as mass accretion continues through the disk to the center. The shape of the dust sublimation radius is not spherical and depends on the angle (θ) from the disk plane, reflecting the nonspherical radiation field and nonuniform dust-free gas. Moreover, we found that the sublimation radius of θ ∼ 20°–60° varies on a timescale of several years. The "inflow-induced outflow" contributes to the obscuration of the nucleus in the sub-parsec region. The column density of the dust-free gas is N H ≳ 1022 cm−2 for r ≲ 0.04 pc. Gases near the disk plane (θ ≲ 30°) can be the origin of the Compton-thick component, which was suggested by the recent X-ray observations of AGNs. The dusty outflow from the sub-parsec region can be also a source of material for the radiation-driven fountain for a larger scale. [ABSTRACT FROM AUTHOR]

Low-ionization-state Mg ii gas has been extensively studied in quasar sightline observations to understand the cool, ∼104 K gas in the circumgalactic medium. Motivated by recent observations showing that the Mg ii gas around low-redshift galaxies has significant angular momentum, we use the high-resolution EAGLE cosmological simulation to analyze the morphological and rotation structures of the z 0.3 circumgalactic Mg ii gas and examine how they change with the host galaxy properties. Around star-forming galaxies, we find that the Mg ii gas has an axisymmetric instead of a spherical distribution, and the axis of symmetry aligns with that of the Mg ii gas rotation. A similar rotating structure is less commonly found in the small sample of simulated quiescent galaxies. We also examine how often Mg ii gas around galaxies selected using a line-of-sight velocity cut includes gas physically outside of the virial radius (r vir). For example, we show that at an impact parameter of 100 pkpc, a ±500 km s-1 velocity cut around galaxies with stellar masses of 109-109.5 M o (1010-1010.5 M o) selects Mg ii gas beyond the virial radius 80% (6%) of the time. Because observers typically select Mg ii gas around target galaxies using such a velocity cut, we discuss how this issue affects the study of circumgalactic Mg ii gas properties, including the detection of corotation. While the corotating Mg ii gas generally extends beyond 0.5r vir, the Mg ii gas outside of the virial radius contaminates the corotation signal and makes observers less likely to conclude that gas at large impact parameters (e.g., 0.25r vir) is corotating.

We investigate the structure of laterally propagating flames through the highly stratified burning layer in an X-ray burst. Two-dimensional hydrodynamics simulations of flame propagation are performed through a rotating plane-parallel atmosphere, exploring the structure of the flame. We discuss the approximations needed to capture the length and timescales at play in an X-ray burst and describe the flame acceleration observed. Our studies complement other multidimensional studies of burning in X-ray bursts.

We study close encounters of a 1 M ⊙ middle-age main-sequence star (modeled using MESA) with massive black holes through hydrodynamic simulations, and explore in particular the dependence of the outcomes on the black hole mass. We consider here black holes in the intermediate-mass range, M BH = 100–104 M ⊙. Possible outcomes vary from a small tidal perturbation for weak encounters all the way to partial or full disruption for stronger encounters. We find that stronger encounters lead to increased mass loss at the first pericenter passage, in many cases ejecting the partially disrupted star on an unbound orbit. For encounters that initially produce a bound system, with only partial stripping of the star, the fraction of mass stripped from the star increases with each subsequent pericenter passage and a stellar remnant of finite mass is ultimately ejected in all cases. The critical penetration depth that separates bound and unbound remnants has a dependence on the black hole mass when M BH ≲ 103 M ⊙. We also find that the number of successive close passages before ejection decreases as we go from the stellar-mass black hole to the intermediate-mass black hole regime. For instance, after an initial encounter right at the classical tidal disruption limit, a 1 M ⊙ star undergoes 16 (5) pericenter passages before ejection from a 10 M ⊙ (100 M ⊙) black hole. Observations of periodic flares from these repeated close passages could in principle indicate signatures of a partial tidal disruption event. [ABSTRACT FROM AUTHOR]

Rotating massive stars with initial progenitor masses M prog ∼ 25–140 M ⊙ can leave rapidly rotating black holes to become collapsars. The black holes and the surrounding accretion disks may develop powerful jets by magnetohydrodynamics instabilities. The propagation of the jet in the stellar envelope provides the necessary shock heating for triggering nucleosynthesis unseen in canonical core-collapse supernovae. However, the energy budget of the jet and its effects on the final chemical abundance pattern are unclear. In this exploratory work, we present a survey on the parameter dependence of collapsar nucleosynthesis on jet energetics. We use the zero-metallicity star with M prog ∼ 40 M ⊙ as the progenitor. The parameters include the jet duration, its energy deposition rate, deposited energy, and the opening angle. We examine the correlations of the following observables: (1) the ejecta and remnant masses; (2) the energy deposition efficiency; (3) the 56Ni production and its correlation with the ejecta velocity, deposited energy, and the ejected mass; (4) the Sc–Ti–V correlation as observed in metal-poor stars; and (5) the [Zn/Fe] ratio as observed in some metal-poor stars. We also provide the chemical abundance table of these explosion models for the use of the galactic chemical evolution and stellar archeology. [ABSTRACT FROM AUTHOR]

Galaxy evolution can be dramatically affected by the environment, especially by the dense environment of a galaxy cluster. Recent observational studies show that massive galaxies undergoing strong ram-pressure stripping (RPS) have an enhanced frequency of nuclear activity. We investigate this topic using a suite of wind-tunnel hydrodynamical simulations of a massive M star = 1011 M ⊙ disk galaxy with 39 pc resolution and including star formation and stellar feedback. We find that RPS increases the inflow of gas to the galaxy center regardless of the wind impact angle. This increase is driven by the mixing of interstellar and nonrotating intracluster media at all wind angles, and by increased torque on the inner disk gas, mainly from local pressure gradients when the intracluster medium (ICM) wind has an edge-on component. In turn, the strong pressure torques are driven by rising ram pressure. We estimate the black hole (BH) accretion using Bondi–Hoyle and torque models, and compare it with the mass flux in the central 140 pc region. We find that the torque model predicts much less accretion onto the BH of a RPS galaxy than the Bondi–Hoyle estimator. We argue that both models are incomplete: the commonly used torque model does not account for torques caused by the gas distribution or local pressure gradients, while the Bondi–Hoyle estimator depends on the sound speed of the hot gas, which includes the ICM in stripped galaxies. An estimator that accounts for this missing physics is required to capture BH accretion in a RPS galaxy. [ABSTRACT FROM AUTHOR]

Observations of protoplanetary disks have shown that forming exoplanets leave characteristic imprints on the gas and dust of the disk. In the gas, these forming exoplanets cause deviations from Keplerian motion, which can be detected through molecular line observations. Our previous work has shown that machine learning can correctly determine if a planet is present in these disks. Using our machine-learning models, we identify strong, localized non-Keplerian motion within the disk HD 142666. Subsequent hydrodynamics simulations of a system with a 5 M J planet at 75 au recreate the kinematic structure. By currently established standards in the field, we conclude that HD 142666 hosts a planet. This work represents a first step toward using machine learning to identify previously overlooked non-Keplerian features in protoplanetary disks. [ABSTRACT FROM AUTHOR]

Based on independent shear measurements using the Dark Energy Camera Legacy Survey/DR8 imaging data, we measure the weak lensing signals around isolated central galaxies (ICGs) from Sloan Digital Sky Survey/DR7 at z ∼ 0.1. The projected stellar mass density profiles of satellite galaxies are further deduced, using photometric sources from the Hyper Suprime-cam survey. The signals of ICGs + their extended stellar halos are taken from Wang et al. All measurements are compared with predictions by the IllustrisTNG300-1 simulation. We find, overall, a good agreement between observation and TNG300. In particular, a correction to the stellar mass of massive observed ICGs is applied based on the calibration of He et al., which brings a much better agreement with TNG300 predicted lensing signals at log 10 M * / M ⊙ > 11.1. In real observation, red ICGs are hosted by more massive dark matter halos and have more satellites and more extended stellar halos than blue ICGs at fixed stellar mass. However, in TNG300 there are more satellites around blue ICGs at fixed stellar mass, and the outer stellar halos of red and blue ICGs are similar. The stellar halos of TNG galaxies are more extended compared with real observed galaxies, especially for blue ICGs with log 10 M * / M ⊙ > 10.8. We find the same trend for TNG100 galaxies and for true halo central galaxies. The tensions between TNG and real galaxies indicate that satellite disruptions are stronger in TNG. In both TNG300 and observation, satellites approximately trace the underlying dark matter distribution beyond 0.1 R 200, but the fraction of total stellar mass in TNG300 does not show the same radial distribution as real galaxies. [ABSTRACT FROM AUTHOR]

Giant planets embedded in protoplanetary disks (PPDs) can create annulus density gaps around their orbits in the type-II regime, potentially responsible for the ubiquity of annular substructures observed in PPDs. Although a substantial amount of works studying type-II planetary migration and gap properties have been published, they have almost exclusively all been conducted under the viscous accretion disk framework. However, recent studies have established magnetized disk winds as the primary mechanism driving disk accretion and evolution, which can coexist with turbulence from the magnetorotational instability (MRI) in the outer PPDs. We conduct a series of 3D global nonideal magnetohydrodynamic (MHD) simulations of type-II planet–disk interactions applicable to the outer PPDs. Our simulations properly resolve the MRI turbulence and accommodate the MHD disk wind. We found that the planet triggers the poloidal magnetic flux concentration around its orbit. The concentrated magnetic flux strongly enhances angular momentum removal in the gap, which is along the inclined poloidal field through a strong outflow emanating from the disk surface outward to the planet gap. The resulting planet-induced gap shape is more similar to an inviscid disk, while being much deeper, which can be understood from a simple inhomogeneous wind torque prescription. The corotation region is characterized by a fast trans-sonic accretion flow that is asymmetric in azimuth about the planet and lacking the horseshoe turns, and the meridional flow is weakened. The torque acting on the planet generally drives inward migration, though the migration rate can be affected by the presence of neighboring gaps through stochastic, planet-free magnetic flux concentration. [ABSTRACT FROM AUTHOR]

We present 3D hydrodynamical simulations of core convection with a stably stratified envelope of a 25 M _⊙ star in the early phase of the main sequence. We use the explicit gas-dynamics code PPMstar , which tracks two fluids and includes radiation pressure and radiative diffusion. Multiple series of simulations with different luminosities and radiative thermal conductivities are presented. The entrainment rate at the convective boundary, internal gravity waves in and above the boundary region, and the approach to dynamical equilibrium shortly after a few convective turnovers are investigated. We perform very long simulations on 896 ^3 grids accelerated by luminosity boost factors of 1000, 3162 and 10,000. In these simulations, the growing penetrative convection reduces the initially unrealistically large entrainment. This reduction is enabled by a spatial separation that develops between the entropy gradient and the composition gradient. The convective boundary moves outward much more slowly at the end of these simulations. Finally, we present a 1D method to predict the extent and character of penetrative convection beyond the Schwarzschild boundary. The 1D model is based on a spherically averaged reduced entropy equation that takes the turbulent dissipation as input from the 3D hydrodynamic simulation and takes buoyancy and all other energy sources and sinks into account. This 1D method is intended to be ultimately deployed in 1D stellar evolution calculations and is based on the properties of penetrative convection in our simulations carried forward through the local thermal timescale.

Observed high multiplicity planetary systems are often tightly packed. Numerical studies indicate that such systems are susceptible to dynamical instabilities. Dynamical instabilities in close-in tightly packed systems, similar to those found in abundance by Kepler, often lead to planet–planet collisions. For sub-Neptunes, the dominant type of observed exoplanets, the planetary mass is concentrated in a rocky core, but the volume is dominated by a low-density gaseous envelope. For these, using the traditional perfect merger assumption (also known as the “sticky-sphere” approximation) to resolve collisions is questionable. Using both N -body integration and smoothed-particle hydrodynamics, we have simulated sub-Neptune collisions for a wide range in realistic kinematic properties such as impact parameters ( $b^{\prime} $ ) and impact velocities ( ${v}_{\mathrm{im}}$ ) to study the possible outcomes in detail. We find that the majority of the collisions with kinematic properties similar to what is expected from dynamical instabilities in multiplanet systems may not lead to mergers of sub-Neptunes. Instead, both sub-Neptunes survive the encounter, often with significant atmosphere loss. When mergers do occur, they can involve significant mass loss and can sometimes lead to complete disruption of one or both planets. Sub-Neptunes merge or disrupt if $b^{\prime} \lt {b}_{{\rm{c}}}^{{\prime} }$ , a critical value dependent on ${v}_{\mathrm{im}}$ / v _esc , where v _esc is the escape velocity from the surface of the hypothetical merged planet assuming perfect merger. For ${v}_{\mathrm{im}}/{v}_{\mathrm{esc}}\lesssim 2.5$ , ${b}_{{\rm{c}}}^{{\prime} }\propto {({v}_{\mathrm{im}}/{v}_{\mathrm{esc}})}^{-2}$ , and collisions with $b^{\prime} \lt {b}_{{\rm{c}}}^{{\prime} }$ typically leads to mergers. On the other hand, for ${v}_{\mathrm{im}}/{v}_{\mathrm{esc}}\gtrsim 2.5$ , ${b}_{{\rm{c}}}^{{\prime} }$ ∝ ${v}_{\mathrm{im}}$ / v _esc , and the collisions with $b^{\prime} \lt {b}_{{\rm{c}}}^{{\prime} }$ can result in complete destruction of one or both sub-Neptunes.

Stellar bars in disk galaxies grow as stars in near-circular orbits lose angular momentum to their environments, including their dark matter (DM) halo, and transform into elongated bar orbits. This angular momentum exchange during galaxy evolution hints at a connection between bar properties and the DM halo spin λ , the dimensionless form of DM angular momentum. We investigate the connection between halo spin λ and galaxy properties in the presence/absence of stellar bars, using the cosmological magnetohydrodynamic TNG50 simulations at multiple redshifts (0 < z _r < 1). We determine the bar strength (or bar amplitude, A _2 / A _0 ), using Fourier decomposition of the face-on stellar density distribution. We determine the halo spin for barred and unbarred galaxies (0 < A _2 / A _0 < 0.7) in the center of the DM halo, close to the galaxy’s stellar disk. At z _r = 0, there is an anticorrelation between halo spin and bar strength. Strongly barred galaxies ( A _2 / A _0 > 0.4) reside in DM halos with low spin and low specific angular momentum at their centers. In contrast, unbarred/weakly barred galaxies ( A _2 / A _0 < 0.2) exist in halos with higher central spin and higher specific angular momentum. The anticorrelation is due to the barred galaxies’ higher DM mass and lower angular momentum than the unbarred galaxies at z _r = 0, as a result of galaxy evolution. At high redshifts ( z _r = 1), all galaxies have higher halo spin compared to those at lower redshifts ( z _r = 0), with a weak anticorrelation for galaxies having A _2 / A _0 > 0.2. The formation of DM bars in strongly barred systems highlights how angular momentum transfer to the halo can influence its central spin.

Recent JWST observations with superb angular resolution have revealed the existence of clumpy galaxies at high redshift through the detection of rest-frame optical emission lines. We use the FirstLight simulation to study the properties of (sub)galactic clumps that are bright in the [O III ] λ 5007 line with flux greater than ∼10 ^−18 erg s ^−1 cm ^−2 , to be detected by JWST. For 62 simulated galaxies that have stellar masses of (0.5–6) × 10 ^10 M _⊙ at z = 5, we find clumps in 1828 snapshots in the redshift range z = 9.5–5.5. The clumps are identified by the surface density of the star formation rate (SFR). About one-tenth of the snapshots show the existence of clumpy systems with two or more components. Most of the clumps are formed by mergers and can be characterized by their ages: central clumps dominated by stellar populations older than 50 Myr, and off-centered clumps dominated by younger stellar populations with specific SFRs of ∼50 Gyr ^−1 . The latter type of young clumps is formed from gas debris in the tidal tails of major mergers with baryonic mass ratios of 1 ≤ q < 4. The merger-induced clumps are short-lived and merge within a dynamical time of several tens of million years. The number density of the clumpy systems is estimated to be ∼10 ^−5 cMpc ^−3 , which is large enough to be detected in recent JWST surveys.