Magnetohydrodynamic Simulations Research Articles

Magnetohydrodynamic Simulation and Particle Simulation for Vacuum Circuit Breakers

Magnetohydrodynamic Simulation and Particle Simulation for Vacuum Circuit Breakers

Dynamics and emission properties of flux ropes from two-temperature GRMHD simulations with multiple magnetic loops

Context. Flux ropes erupting in the vicinity of a black hole are thought to be a potential model for the flares observed in Sagittarius A*. Aims. In this study, we examine the radiative properties of flux ropes that emerged in the vicinity of a black hole. Methods. We performed three-dimensional two-temperature general relativistic magnetohydrodynamic (GRMHD) simulations of magnetized accretion flows with alternating multiple magnetic loops and general relativistic radiation transfer (GRRT) calculations. In the GRMHD simulations, we implemented two different sizes of initial magnetic loops. Results. In the small loop case, magnetic dissipation leads to a weaker excitement of magneto-rotational instability inside the torus, which generates a lower accretion rate compared to the large loop case. However, the small loop case generates more flux ropes due to frequent reconnection by magnetic loops with different polarities. By calculating the thermal synchrotron emission, we found that the variability of light curves and the emitting region are tightly related. At 230 GHz and higher, the emission from the flux ropes is relatively stronger compared to the background, which is responsible for the filamentary structure in the images. At lower frequencies (e.g., 43 GHz), emission comes from more extended regions, which have a less filamentary structure in the image. Conclusions. Our study shows that self-consistent electron temperature models are essential for the calculation of thermal synchrotron radiation and the morphology of the GRRT images. Flux ropes contribute considerable emission at frequencies ≳230 GHz.

The Impact of Stellar Radiative Feedback on Formation of Young Massive Clusters via Fast H i Gas Collisions

Young massive clusters (YMCs) are dense aggregates of young stars and are often speculated as potential precursors to globular clusters. However, the formation mechanism of massive and compact gas clumps that precede YMCs remains unknown. In this paper, we study the formation of such massive clumps via fast H i gas collisions (∼100 km s−1) as suggested by recent observations and their subsequent evolution into YMCs by using three-dimensional magnetohydrodynamics simulations involving self-gravity and detailed thermal/chemical processes. In particular, the impact of ionization feedback from stellar radiation is included in an approximate fashion where the temperature within the H ii regions is elevated to 10,000 K, while supernova feedback is not included. We examine whether the resulting massive clumps can survive this ionization feedback and evolve into YMCs. Our simulations reveal the emergence of gas clumps that not only possess substantial mass (∼105 M ⊙) but also sufficient compactness (∼5 pc). Notably, these clumps exhibit significantly higher escape velocities compared to the sound speed of the H ii region, indicating effective gravitational retention of gas against feedback-induced evaporation. Consequently, these conditions foster efficient star formation within the massive gas clumps, ultimately leading to their evolution into YMCs. We also perform simulations involving lower-velocity gas collisions, approximately 15 km s−1, typical shock velocities induced by galactic superbubbles. In contrast to the high-velocity collisions, we find that molecular cloud formation does not occur in the case of 1 cm−3 gas collision, while YMC formation is observed in the presence of denser gas of 10 cm−3. However, the formation of YMCs requires compression periods exceeding 10 Myr in these cases, indicating a potential preference for gas collisions driven by intergalactic interactions rather than galactic superbubbles for YMC formation.

Observations of Locally Excited Waves in the Low Solar Atmosphere Using the Daniel K. Inouye Solar Telescope

We present an interpretation of the recent Daniel K. Inouye Solar Telescope (DKIST) observations of propagating wave fronts in the lower solar atmosphere. Using MPS/University of Chicago MHD radiative magnetohydrodynamic simulations spanning the solar photosphere, the overshoot region, and the lower chromosphere, we identify three acoustic-wave source mechanisms, each occur at a different atmospheric height. We synthesize the DKIST Visible Broadband Imager G-band, blue-continuum, and Ca ii K signatures of these waves at high spatial and temporal resolution, and conclude that the wave fronts observed by DKIST likely originate from acoustic sources at the top of the solar photosphere overshoot region and in the chromosphere proper. The overall importance of these local sources to the atmospheric energy and momentum budget of the solar atmosphere is unknown, but one of the excitation mechanisms identified (upward propagating shock interaction with down-welling chromospheric plasma resulting in acoustic radiation) may be an important shock dissipation mechanism. Additionally, the observed wave fronts may prove useful for ultralocal helioseismological inversions and promise to play an important diagnostic role at multiple atmospheric heights.

Diagnostics of Magnetohydrodynamic Modes in the Interstellar Medium through Synchrotron Polarization Statistics

One of the biggest challenges in understanding magnetohydrodynamic (MHD) turbulence is identifying the plasma mode components from observational data. Previous studies on synchrotron polarization from the interstellar medium (ISM) suggest that the dominant MHD modes can be identified via statistics of Stokes parameters, which would be crucial for studying various ISM processes such as the scattering and acceleration of cosmic rays, star formation, and dynamo. In this paper, we present a numerical study of the synchrotron polarization analysis (SPA) method through systematic investigation of the statistical properties of the Stokes parameters. We derive the theoretical basis for our method from the fundamental statistics of MHD turbulence, recognizing that the projection of the MHD modes allows us to identify the modes dominating the energy fraction from synchrotron observations. Based on the discovery, we revise the SPA method using synthetic synchrotron polarization observations obtained from 3D ideal MHD simulations with a wide range of plasma parameters and driving mechanisms, and present a modified recipe for mode identification. We propose a classification criterion based on a new SPA+ fitting procedure, which allows us to distinguish between Alfvén mode and compressible/slow mode dominated turbulence. We further propose a new method to identify fast modes by analyzing the asymmetry of the SPA+ signature and establish a new asymmetry parameter to detect the presence of fast mode turbulence. Additionally, we confirm through numerical tests that the identification of the compressible and fast modes is not affected by Faraday rotation in both the emitting plasma and the foreground.

Formation Mechanism of Fingers That Protrude Eastward From the Io Plasma Disk During the Interchange Instability

AbstractThe solar wind‐magnetosphere‐ionosphere interaction at Jupiter is reproduced numerically adopting the nine‐component magnetohydrodynamic simulation. Calculations take into account the magnetosphere‐ionosphere coupling, Jovian rotation, and Io plasma source. High‐speed rotating plasma inside restricted magnetospheric space causes expansion and contraction of magnetic field, forming super‐rotation at radial distance 20∼30 Rj and co‐rotation breakdown further outside. Field‐perpendicular current that restores co‐rotational delay beyond 30 Rj is connected via field‐aligned current to the main oval in the ionosphere. Inside 20 Rj, there is almost co‐rotation region (deviation from co‐rotation less than 20 km/s). Particularly within 10 Rj, the deviation from co‐rotation is less than 2 km/s. In the nearly co‐rotating region, the Io plasma forms a disk structure through field‐aligned redistribution. The interchange instability occurs near the outer wall of the Io plasma disk, and instability flow develops to vortex. Through this instability, a part of the centrifugal drift current supporting the Io plasma disk is connected to low‐latitude field‐aligned current that generates beads‐like spots on the lower latitude side of the main oval. Resulting interchange instability comes to satisfy the structure of convection and enables further development of vortex. The Coriolis force acting on eastward flow inside the developing vortex makes this flow protrude further outward, forming eastward bending fingers. Inside 10 Rj, Io plasma transport by the interchange instability becomes slower, despite the center of the disk. Io plasma escapes from the inner magnetosphere with a time constant of 20 days if this slow transport is taken into account.

A Revisited Equilibrium Solution of the Fishbone and Moncrief Torus for Extended General Relativistic Magnetohydrodynamic Simulations

Accretion physics has become more important recently due to the detection of the first horizon-scale images of the supermassive black holes of M 87* and Sgr A* by the Event Horizon Telescope. General relativistic magnetohydrodynamic (GRMHD) simulations of magnetized accretion flows onto a Kerr black hole have been used to interpret them. However, further testing the theory of gravity by using horizon-scale images requires performing consistent GRMHD simulations in non-Kerr spacetime. In this paper, we revisited the hydrodynamical equilibrium solution of the Fishbone and Moncrief (FM) torus that can be used to study any stationary, axisymmetric, vacuum, or nonvacuum spacetime. Further, we check the stability of the FM torus in non-Kerr spacetime by general relativistic hydrodynamic simulations. We find that FM torus in non-Kerr spacetime is indeed stable under long-term evolution. We conclude that the generalized FM torus solution would be very useful for creating new GRMHD libraries in extended Kerr black holes.

Joint action of phase mixing and nonlinear effects in MHD waves propagating in coronal loops

Context. The evolution of Alfvén waves in cylindrical magnetic flux tubes, which represent a basic model for loops observed in the solar corona, can be affected by phase mixing and turbulent cascade. Phase mixing results from transverse inhomogeneities in the Alfvén speed, causing different shells of the flux tube to oscillate at different frequencies, thus forming increasingly smaller spatial scales in the direction perpendicular to the guide field. Turbulent cascade also contributes to the dissipation of the bulk energy of the waves through the generation of smaller spatial scales. Both processes present characteristic timescales. Different regimes can be envisaged according to how those timescales are related and to the typical timescale at which dissipation is at work. Aims. We investigate the interplay of phase mixing and the nonlinear turbulent cascade in the evolution and dissipation of Alfvén waves using compressible magnetohydrodynamic numerical simulations. We consider perturbations in the form of torsional waves, both propagating and standing, or turbulent fluctuations, or a combination of the two. The main purpose is to study how phase mixing and nonlinear couplings jointly work to produce small scales in different regimes. Methods. We conducted a numerical campaign to explore the typical parameters, such as the loop length, the amplitude and spatial profile of the perturbations, and the dissipative coefficients. A pseudo-spectral code was employed to solve the three-dimensional compressible magnetohydrodynamic equations, modeling the evolution of perturbations propagating in a flux tube corresponding to an equilibrium configuration with cylindrical symmetry. Results. We find that phase mixing takes place for moderate amplitudes of the turbulent component even in a distorted, nonaxisymmetric configuration, building small scales that are locally transverse to the density gradient. The dissipative time decreases with increasing the percentage of the turbulent component. This behavior is verified both for propagating and standing waves. Even in the fully turbulent case, a mechanism qualitatively similar to phase mixing occurs: it actively generates small scales together with the nonlinear cascade, thus providing the shortest dissipative time. General considerations are given to identify this regime in the parameter space. The turbulent perturbation also distorts the background density, locally increasing the Alfvén velocity gradient and further contributing to accelerating the formation of small scales. Conclusions. Our campaign of simulations is relevant for the coronal plasma where Reynolds and Lundquist numbers are extremely high. For sufficiently low perturbation amplitudes, phase mixing and turbulence work synergically, speeding up the dissipation of the perturbation energy: phase mixing dominates at early times and nonlinear effects at later times. We find that the dissipative time is shorter than those of phase mixing and the nonlinear cascade when individually considered.

Resolution analysis of magnetically arrested disk simulations

Abstract Polarisation measurements by the Event Horizon Telescope from M87* and Sgr A*suggest that there is a dynamically strong, ordered magnetic field, typical of what is expected of a magnetically arrested accretion disk (MAD). In such disks the strong poloidal magnetic field can suppress the accretion flow and cause episodic flux eruptions. Recent work shows that General Relativistic Magnetohydrodynamic (GRMHD) MAD simulations feature dynamics of turbulence and mixing instabilities that are becoming resolved at higher resolutions. We perform a convergence study of MAD states exceeding the status quo by an order of magnitude in resolution. We use existing 3D simulations performed with the H-AMR code, up to resolution of 5376 × 2304 × 2304 in a logarithmic spherical-polar grid. We find consistent time-averaged disk properties across all resolutions. However, higher resolutions reveal signs of inward angular momentum transport attributed to turbulent convection, particularly evident when mixing instabilities occur at the surfaces of flux tubes during flux eruptions. Additionally, we see wave-like features in the jet sheath, which become more prominent at higher resolutions, that may induce mixing between jet and disk. At higher resolutions, we observe the sheath to be thinner, resulting in increased temperature, reduced magnetisation, and greater variability. Those differences could affect the dissipation of energy, that would eventually result in distinct observable radiative emission from high-resolution simulations. With higher resolutions, we can delve into crucial questions about horizon-scale physics and its impact on the dynamics and emission properties of larger-scale jets.

Refueling of field-reversed configuration core via axial plasmoids injection

This study successfully developed a refueling technique for a field-reversed configuration (FRC) via axial plasmoid injection and demonstrated it on the FAT-CM device at Nihon University. The target FRC is generated using the collisional-merging formation technique combined with conical theta-pinch formation. Plasmoids with an FRC-like configuration are coaxially injected from both ends of the FAT-CM device toward the preexisting target FRC. Postinjection, the system achieves equilibrium, resulting in increases by factors of 1.8 and 2.4 in the total inventory and plasma energy, respectively, compared to cases without injection. This method effectively accomplishes FRC refueling while preserving the intrinsic characteristics of a simply connected, axisymmetric configuration and a high beta value approaching unity. Therefore, this approach offers potential for repetitive refueling in the reactor stage having a FRC plasma core. Experimental outcomes are compared with magnetohydrodynamic simulation results. In the collisional merging process, the characteristics of the pre-collision plasmoids, such as the strong toroidal rotation and coherent FRC-like magnetic field structures of the FRC, are not preserved. Experimental environments have been constructed to investigate such unique properties of the resulting FRCs.

The role of active galactic nucleus feedback on the evolution of dwarf galaxies from cosmological simulations. Supermassive black holes suppress star formation in low-mass galaxies

Recent observational studies suggest that feedback from active galactic nuclei (AGNs) may play an important role in the formation and evolution of dwarf galaxies, an issue that has received little attention from a theoretical perspective. We investigated this using two sets of 12 cosmological magnetohydrodynamic simulations of the formation of dwarf galaxies (10$^ M_ M$_ M_ one set using a version of the AURIGA galaxy formation physics model including AGN feedback and a parallel set with AGN feedback turned off. We show that the full-physics AGN runs satisfactorily reproduce several scaling relations, including the black-hole-to-stellar mass (M$_ BH $-M$_ star $), the black-hole-to-sigma (M$_ BH and the baryonic Tully-Fisher relation. We find that the global star formation (SF) of galaxies run with an AGN is reduced compared to the one in which the AGN has been turned off, suggesting that AGN feedback is a viable way of suppressing SF in dwarf galaxies, even though none of our galaxies is completely quenched by $z$=$0$. Furthermore, we find a tight correlation between the median SF rates and the BH $/M$_ star $ ratio in our simulated dwarfs. Star formation is suppressed due to gas heating in the vicinity of the AGN: less HI gas is available in AGN runs, though the total amount of gas is preserved across the two settings within each galaxy. This indicates that the main effect of AGN feedback in our dwarfs is to heat up and push the gas away from the galaxy's centre rather than expelling it completely. Finally, we show that the two galaxies harbouring the largest supermassive black holes have suffered a considerable (up to sim 65$\,<!PCT!>$) reduction in their central dark matter density, pinpointing the role of AGNs in determining the final dark matter mass distribution within dwarf galaxies. This pilot paper highlights the importance of modelling AGN feedback at the lowest mass scales and the impact this can have on dwarf galaxy evolution.

The gravitational-wave emission from the explosion of a 15 solar mass star with rotation and magnetic fields

Abstract Gravitational waveform predictions from 3D simulations of explosions of non-rotating massive stars with no magnetic fields have been extensively studied. However, the impact of magnetic fields and rotation on the core-collapse supernova gravitational-wave signal is not well understood beyond the core-bounce phase. Therefore, we perform four magnetohydrodynamical simulations of the explosion of a 15 M⊙ star with the SFHx and SFHo equations of state. All of the models start with a weak magnetic field strength of 108 G, and two of the models are rapidly rotating. We discuss the impact of the rotation and magnetic fields on the gravitational-wave signals. We find that the weak pre-collapse fields do not have a significant impact on the gravitational-wave signal amplitude. With rapid rotation, the f/g-mode trajectory can change in shape, and the dominant emission band becomes broader. We include the low-frequency memory component of the gravitational-wave signal from both matter motions and neutrino emission anisotropy. We show that including the gravitational waves from anisotropic neutrino emission increases the supernova detection distances for the Einstein Telescope, and would also be detectable out to Mpc distances by a moon-based gravitational-wave detector.

A Structure-Preserving Semi-implicit IMEX Finite Volume Scheme for Ideal Magnetohydrodynamics at all Mach and Alfvén Numbers

We present a divergence-free semi-implicit finite volume scheme for the simulation of the ideal magnetohydrodynamics (MHD) equations which is stable for large time steps controlled by the local transport speed at all Mach and Alfvén numbers. An operator splitting technique allows to treat the convective terms explicitly while the hydrodynamic pressure and the magnetic field contributions are integrated implicitly, yielding two decoupled linear implicit systems. The linearity of the implicit part is achieved by means of a semi-implicit time linearization. This structure is favorable as second-order accuracy in time can be achieved relying on the class of semi-implicit IMplicit–EXplicit Runge–Kutta (IMEX-RK) methods. In space, implicit cell-centered finite difference operators are designed to discretely preserve the divergence-free property of the magnetic field on three-dimensional Cartesian meshes. The new scheme is also particularly well suited for low Mach number flows and for the incompressible limit of the MHD equations, since no explicit numerical dissipation is added to the implicit contribution and the time step is scale independent. Likewise, highly magnetized flows can benefit from the implicit treatment of the magnetic fluxes, hence improving the computational efficiency of the novel method. The convective terms undergo a shock-capturing second order finite volume discretization to guarantee the effectiveness of the proposed method even for high Mach number flows. The new scheme is benchmarked against a series of test cases for the ideal MHD equations addressing different acoustic and Alfvén Mach number regimes where the performance and the stability of the new scheme is assessed.

Hybrid GRMHD and force-free simulations of black hole accretion

ABSTRACT We present a new approach for stably evolving general relativistic magnetohydrodynamic (GRMHD) simulations in regions where the magnetization $\sigma =b^2/\rho c^2$ becomes large. GRMHD codes typically struggle to evolve plasma above $\sigma \approx 100$ in simulations of black hole accretion. To ensure stability, GRMHD codes will inject mass density artificially to the simulation as necessary to keep the magnetization below a ceiling value $\sigma _{\rm max}$. We propose an alternative approach where the simulation transitions to solving the equations of general relativistic force-free electrodynamics (GRFFE) above a magnetization $\sigma _{\rm trans}$. We augment the GRFFE equations in the highly magnetized region with approximate equations to evolve the decoupled field-parallel velocity and plasma energy density. Our hybrid scheme is explicit and easily added to the framework of standard-volume GRMHD codes. We present a variety of tests of our method, implemented in the GRMHD code koral, and we show results from a 3D hybrid GRMHD + GRFFE simulation of a magnetically arrested disc (MAD) around a spinning black hole. Our hybrid MAD simulation closely matches the average properties of a standard GRMHD MAD simulation with the same initial conditions in low magnetization regions, but it achieves a magnetization $\sigma \approx 10^6$ in the evacuated jet funnel. We present simulated horizon-scale images of both simulations at 230 GHz with the black hole mass and accretion rate matched to M87*. Images from the hybrid simulation are less affected by the choice of magnetization cut-off $\sigma _{\rm cut}$ imposed in radiative transfer than images from the standard GRMHD simulation.

Propagating kink waves in an open coronal magnetic flux tube with gravitational stratification: Magnetohydrodynamic simulation and forward modelling

In coronal open-field regions, such as coronal holes, there are many transverse waves propagating along magnetic flux tubes, which are generally interpreted as kink waves. Previous studies have highlighted their potential role in coronal heating, solar wind acceleration, and seismological diagnostics of various physical parameters. This study aims to investigate propagating kink waves, considering both vertical and horizontal density inhomogeneity, using 3D magnetohydrodynamic (MHD) simulations. We established a 3D MHD model of a gravitationally stratified open flux tube, incorporating a velocity driver at the lower boundary to excite propagating kink waves. Forward modelling was conducted to synthesise observational signatures of the Fe ix 17.1 nm line. Resonant absorption and density stratification both affect the wave amplitude. When diagnosing the relative density profile with velocity amplitude, resonant damping needs to be properly considered to avoid a possible underestimation. In addition, unlike standing modes, propagating waves are believed to be Kelvin-Helmholtz stable. In the presence of vertical stratification, however, the phase mixing of transverse motions around the tube boundary can still induce small-scale structures, partially dissipating wave energy and leading to a temperature increase, especially at higher altitudes. Moreover, we conducted forward modeling to synthesise observational signatures, which revealed the promising potential of future coronal imaging spectrometers such as MUSE in resolving these wave-induced signatures. Also, the synthesised intensity signals exhibit apparent periodic variations, offering a potential method for indirectly observing propagating kink waves with current extreme ultraviolet imagers.

Wide-binary eccentricity distribution in young star clusters: dependence on the binary separation and mass

Abstract We study the wide-binary eccentricity (e) distribution in young star clusters and the role of turbulence in setting the form of the e distribution using magnetohydrodynamical (MHD) simulations of star cluster formation. The simulations incorporate gravity, turbulence, magnetic fields, protostellar heating, and jets/outflows. We find that (1) simulations that employ purely compressive turbulence driving produce binaries with a superthermal e distribution (α > 1 in p(e)∝eα), while simulations with purely solenoidal driving or natural mixture of driving modes produce subthermal/thermal distributions (α ≤ 1), (2) the e distribution over the full range of binary separations in our simulations is set at the early stages of the star cluster formation process, (3) while binaries (separation of rpair ≤ 1000 AU) have subthermal to thermal e distributions (α ∼ 0.8), wide binaries (rpair > 1000 AU) have a superthermal distribution (α ∼ 1.8), and (4) low-mass binary systems (system masses of Msys ≤ 0.8 M⊙) have a highly superthermal distribution (α ∼ 2.4), whereas high-mass systems (Msys > 0.8 M⊙) exhibit a subthermal/thermal distribution (α ∼ 0.8). The binary eccentricity distribution is often modelled as a thermal distribution. However, our results suggest that the e distribution depends on the range of separation of the sampled binaries, which agrees with the findings from recent Gaia observations. We conclude that the dependence of the e distribution on the binary separation and mass is linked to the binary formation mechanism governed by the turbulent properties of the parent cloud.

The three-dimensional structure of black hole accretion flows within the plunging region

Abstract We analyse, using new analytical models and numerical general relativistic magnetohydrodynamic simulations, the three-dimensional properties of accretion flows inside the plunging region of black hole spacetimes (i.e., at radii smaller than the innermost stable circular orbit). These simulations are of thick discs, with aspect ratios of order unity h/r ∼ 1, and with a magnetic field geometry given by the standard low-magnetization ‘SANE’ configuration. This work represents the first step in a wider analysis of this highly relativistic region. We show that analytical expressions derived in the ‘thin disc’ limit describe the numerical results remarkably well, despite the large aspect ratio of the flow. We further demonstrate that accretion within this region is typically mediated via spiral arms, and that the geometric properties of these spiral structures can be understood with a simple analytical model. These results highlight how accretion within the plunging region is fundamentally two dimensional in character, which may have a number of observational implications. We derive a modified theoretical description of the pressure within the plunging region which accounts for turbulent heating and may be of use to black hole image modelling.

Inner dusty regions of protoplanetary discs - III. The role of non-radial radiation pressure in dust dynamics

Abstract We explore dynamical behaviour of dust particles that populate the surface of inner optically thick protoplanetary discs. This is a disc region with the hottest dust and of a great importance for planet formation and dust evolution, but we still struggle to understand all the forces that shape this environment. In our approach we combine results from two separate numerical studies - one is the wind velocity and density distributions obtained from magnetohydrodynamical simulations of accretion discs, and the other is a high-resolution multigrain dust radiation transfer. In our previous paper in the series, we described the methodology for utilising these results as an environmental input for the integration of dust trajectories driven by gravity, gas drag, and radiation pressure. Now we have two improvements - we incorporate time changes in the wind density and velocity, and we implement the non-radial radiation pressure force. We applied our analysis on Herbig Ae and T Tau stars. We confirm that the radiation pressure force can lead to dust outflow, especially in the case of more luminous stars. Additionally, it opposes dust accretion at the inner disc edge and reduces dust settling. These effects are enhanced by the disc wind, especially in the zone where the stellar and the disc magnetic fields meet. Our results suggest that dust grains can stay in the hottest disc region for an extended period and then can end up ejected into the outer disc regions.

Unravelling jet quenching criteria across L* galaxies and massive cluster ellipticals

ABSTRACT In the absence of supplementary heat, the radiative cooling of halo gas around massive galaxies (Milky Way mass and above) leads to an excess of cold gas or stars beyond observed levels. Active galactic nucleus jet-induced heating is likely essential, but the specific properties of the jets remain unclear. Our previous work concludes from simulations of a halo with $10^{14} \,\mathrm{ M}_\odot$ that a successful jet model should have an energy flux comparable to the free-fall energy flux at the cooling radius and should inflate a sufficiently wide cocoon with a long enough cooling time. In this paper, we investigate three jet modes with constant fluxes satisfying the criteria, including high-temperature thermal jets, cosmic ray (CR)-dominant jets, and widely precessing kinetic jets in $10^{12}-10^{15}\, {\rm M}_{\odot }$ haloes using high-resolution, non-cosmological magnetohydrodynamic simulations with the FIRE-2 (Feedback In Realistic Environments) stellar feedback model, conduction, and viscosity. We find that scaling the jet energy according to the free-fall energy at the cooling radius can successfully suppress the cooling flows and quench galaxies without violating observational constraints. On the contrary, if we scale the energy flux based on the total cooling rate within the cooling radius, strong interstellar medium cooling dominates this scaling, resulting in a jet flux exceeding what is needed. Among the three jet types, the CR-dominant jet is most effective in suppressing cooling flows across all surveyed halo masses due to enhanced CR pressure support. We confirm that the criteria for a successful jet model work across a wider range, encompassing halo masses of $10^{12}-10^{15} {\rm M_\odot }$.

The Dynamic Response of Jovian Magnetotail Reconnection to Enhanced Solar Wind Ram Pressure

In this study, we employ the Grid Agnostic Magnetohydrodynamic (MHD) for Extended Research Applications (GAMERA), a high-resolving-power, three-dimensional global MHD model, to simulate magnetotail reconnection in Jupiter's magnetosphere. While previous satellite observations have provided initial statistics on magnetotail reconnection properties at Jupiter, they have been limited in spacetime coverage, leaving the dynamic process of Jovian magnetotail reconnection and its response to the solar wind (SW) poorly understood. Using MHD simulations, we quantitatively analyze the temporal evolution and spatial dependence of nightside reconnection in Jupiter's magnetotail under ideal quiet and enhanced SW conditions. Our results demonstrate that magnetotail reconnection tends to occur in the midnight and postmidnight sectors, with a low occurrence in the premidnight sector, consistent with both Galileo and Juno observations and predictions by Delamere & Bagenal. The magnetic local time (MLT)–radial distribution of magnetotail reconnection is broad, indicating that Jovian magnetotail reconnection is always dynamic rather than steady state. Enhanced SW ram pressure can decrease the MLT coverage of magnetotail reconnection by compressing Jupiter's magnetosphere. However, the occurrence of magnetotail reconnection near the midnight and postmidnight sectors is enhanced by SW compression beyond 60 R J, but is not significantly impacted by SW compression within 60 R J. Conversely, SW compression suppresses reconnection in the premidnight sector, leading to a stronger dawn–dusk asymmetry in the occurrence and location of magnetotail reconnection. This study validates the applicability of the GAMERA code for simulating Jupiter’s magnetosphere and provides complementary insights into the dynamic structure and the SW response of Jupiter’s magnetosphere.