Global Magnetohydrodynamic Simulations Research Articles

Unveiling the Connection between the Alfvénic Oval and the Open–Closed Field Line Boundary at Ganymede

Alfvén waves generated by the interactions between the Ganymede magnetosphere and the fast-corotating Jovian plasma carry a significant amount of energy. This energy is partially transferred to precipitating particles, which correspond to a feature called the Alfvénic oval. However, the spatial distribution of the Alfvénic oval in the Ganymede system is unclear, and the physical correlation between the Alfvénic oval and the open–closed field line boundary (OCFB) has not been studied quantitatively. Using a high-resolution, three-dimensional global magnetohydrodynamics simulation of Ganymede's magnetosphere, we investigate the spatial distribution of the Alfvénic oval mapped on the ionosphere surface of Ganymede based on upstream measurements from Juno's PJ34 flyby. Our results show that (1) the simulated Alfvénic oval of Ganymede closely resembles the observed aurora oval, and although without kinetic physics, the magnitude of the simulated Afvénic Poynting flux is significantly lower than the estimation based on auroral luminosity; (2) the OCFB aligns closely with the poleward (equatorial) boundary of the enhanced Alfvénic power on the downstream (upstream) side; and (3) the magnetospheric flux tube convection pattern can guide the Alfvénic wave packet to propagate toward latitudes either equatorward or poleward of the OCFB on the downstream/upstream side. These results, presented for the first time in this study, provide testable predictions for both the ongoing Juno mission and future Jovian missions such as JUICE probing the minimagnetosphere system.

What Controls the Dawn‐Dusk Asymmetry of Dayside Alfvénic Power: Interplanetary Magnetic Field or Ionospheric Conductance?

AbstractAlfvén waves are important carriers of solar wind energy into the Earth's magnetosphere and upper atmosphere. The observed distribution of Alfvénic power at low altitude exhibits significant dawn‐dusk asymmetry with different causes and impacts on the dayside and nighside. The nightside asymmetry has been shown to be controlled by the meridional gradient in the ionospheric Hall conductance. It is not clear what controls the dayside asymmetry. Both the ionospheric conductance gradient and orientation of the y‐component of the interplanetary magnetic field (IMF) are possible factors. We examined both possibilities using global magnetohydrodynamic (MHD) simulations. Our results show: (a) The dayside distribution of Alfvénic power is very sensitive to the orientation of IMF By, with enhanced power in the dawn (dusk) sector of the northern hemisphere for positive (negative) By (and opposite in the southern hemisphere); and (b) gradients in ionospheric conductance exert a weaker influence on the dayside asymmetry.

Numerical investigation of magnetic reconnection in Jupiter’s dayside magnetodisc

Dayside magnetodisc reconnection is a unique process that may occur in rapidly rotating magnetospheres with internal plasma sources, and it is not expected at Earth. Several observations suggest that dayside magnetodisc reconnection could be driven in the Kronian magnetosphere. This raises the question as to whether it can occur at Jupiter, the fastest rotating planet in the Solar System. Recent reports have suggested its occurrence, but the limited number of identified events leaves the overall global picture in the Jovian magnetosphere uncertain. The primary objective of this study is to answer whether dayside magnetodisc reconnection could exist within the Jovian magnetosphere through the use of global magnetohydrodynamics (MHD) simulations. Additionally, we aim to conduct a numerical investigation to reveal its general distribution for possible future exploration. The numerical simulations for Jovian magnetospheric dynamics are based on three-dimensional MHD calculations, which enable magnetic reconnection through numerical resistivity. In post-processing procedures, the simulated magnetodisc reconnection event is primarily identified by equatorial $B_ reversals, accompanied by an analysis of local plasma flows and global magnetic field lines. Our global MHD simulations reveal the existence of magnetodisc reconnection within the dayside Jovian magnetosphere. The simulation results indicate that dayside magnetodisc reconnection is more likely to occur during solar wind rarefaction, primarily due to the rapid expansion of the dayside magnetospheric volume. This expansion allows a rotation-driven centrifugal force to generate sufficient magnetic field line stretching for current sheet thinning and reconnection. In simulated dayside events, $B_ reversals may occur locally or non-locally through rotational transport from reconnection sites at earlier magnetic local times. The existence and the features derived from simulations are testable as Galileo/Juno has accumulated a considerable amount of dayside magnetospheric data. These results are valuable for future Jovian missions, such as JUICE and Europa Clipper, and provide new insights into interpreting measurements from other rapidly rotating systems.

Generation of Field‐Aligned Currents in Response to Sudden Enhancement of Solar Wind Dynamic Pressure

AbstractWe investigated the generation of field‐aligned currents (FACs) in response to the sudden enhancement of the solar wind dynamic pressure by tracing backward in time a packet of the Alfvén wave in the global magnetohydrodynamic (MHD) simulation. The generation region is identified from three perspectives, including the continuity of the current, the energy conservation and the time rate of change in the FACs. The generation mechanism is found to be related with the tailward motion of the compressional wave, which is excited when the magnetopause is compressed due to the solar wind dynamic pressure pulse. The compressional wave with a high magnetic pressure center interacts with the Earth's dipole field and forms a protruding part of the wavefront near the equatorial plane. The leading edge of the merged magnetic pressure that pertains to the compressional wave and the Earth starts to generate preliminary impulse (PI) FACs off the equator, as the magnetic pressure force accelerates the plasma and magnetic field lines are bent (FAC dynamo 1). Main impulse (MI) FACs are generated behind the leading edge in the equatorial plane due to the enhanced magnetic tension force (FAC dynamo 2). Magnetic field lines would be extremely curved during the passage due to the increasing magnetic tension force. The polarity of PI FACs is decided by the parallel vorticity of plasma flow, and MI FACs are the result of enhanced perpendicular currents together with the plasma flow in the equatorial plane.

Magnetorotational dynamo can generate large-scale vertical magnetic fields in 3D GRMHD simulations of accreting black holes

ABSTRACT Jetted astrophysical phenomena with black hole engines, including binary mergers, jetted tidal disruption events, and X-ray binaries, require a large-scale vertical magnetic field for efficient jet formation. However, a dynamo mechanism that could generate these crucial large-scale magnetic fields has not been identified and characterized. We have employed three-dimensional global general relativistic magnetohydrodynamical simulations of accretion discs to quantify, for the first time, a dynamo mechanism that generates large-scale magnetic fields. This dynamo mechanism primarily arises from the non-linear evolution of the magnetorotational instability (MRI). In this mechanism, large non-axisymmetric MRI-amplified shearing wave modes, mediated by the axisymmetric azimuthal magnetic field, generate and sustain the large-scale vertical magnetic field through their non-linear interactions. We identify the advection of magnetic loops as a crucial feature, transporting the large-scale vertical magnetic field from the outer regions to the inner regions of the accretion disc. This leads to a larger characteristic size of the, now advected, magnetic field when compared to the local disc height. We characterize the complete dynamo mechanism with two time-scales: one for the local magnetic field generation, $t_{\rm gen}$, and one for the large-scale scale advection, $t_{\rm adv}$. Whereas the dynamo we describe is non-linear, we explore the potential of linear mean field models to replicate its core features. Our findings indicate that traditional $\alpha$-dynamo models, often computed in stratified shearing box simulations, are inadequate and that the effective large-scale dynamics is better described by the shear current effects or stochastic $\alpha$-dynamos.

Dynamics near the inner dead-zone edges in a proprotoplanetary disk

Abstract We perform three-dimensional global non-ideal magnetohydrodynamic simulations of a protoplanetary disk containing the inner dead-zone edge. We take into account realistic diffusion coefficients of the Ohmic resistivity and ambipolar diffusion based on detailed chemical reactions with single-size dust grains. We found that the conventional dead zone identified by the Elsässer numbers of the Ohmic resistivity and ambipolar diffusion is divided into two regions: “the transition zone” and “the coherent zone.” The coherent zone has the same properties as the conventional dead zone, and extends outside of the transition zone in the radial direction. Between the active and coherent zones, we discover the transition zone, the inner edge of which is identical to that of the conventional dead zone. The transition zone extends out over the regions where thermal ionization determines diffusion coefficients. The transition zone has completely different physical properties than the conventional dead zone, the so-called undead zone, and the zombie zone. The combination of amplification of the radial magnetic field owing to the ambipolar diffusion and a steep radial gradient of the Ohmic diffusivity causes the efficient evacuation of the net vertical magnetic flux from the transition zone within several rotations. Surface gas accretion occurs in the coherent zone but not in the transition zone. The presence of the transition zone prohibits mass and magnetic flux transport from the coherent zone to the active zone. Mass accumulation occurs at both edges of the transition zone as a result of mass supply from the active and coherent zones.

Numerical Calculations of Adiabatic Invariants From MHD‐Driven Magnetic Fields

AbstractThe adiabatic invariants (M, J, Φ) and the related invariants (M, K, L*) have been established as effective coordinate systems for describing radiation belt dynamics at a theoretical level, and through numerical techniques, can be paired with in situ observations to order phase‐space density. To date, methods for numerical techniques to calculate adiabatic invariants have focused on empirical models such the Tsyganenko models TS05, T96, and T89. In this work, we develop methods based on numerical integration and variable step size iteration for the calculation of adiabatic invariants, applying the method to the Lyon‐Fedder‐Mobarry (LFM) global magnetohydrodynamics (MHD) simulation code, with optional coupling to the Rice Convection Model (RCM). By opening the door to adiabatic invariant modeling with MHD magnetic fields, the opportunity for exploratory modeling work of radiation belt dynamics is enabled. Calculations performed using LFM are cross‐referenced with the same code applied to the T96 and TS05 Tsyganenko models evaluated on the LFM grid. Important aspects of the adiabatic invariant calculation are reviewed and discussed, including (a) sensitivity to magnetic field model used, (b) differences in the problem between quiet and disturbed geomagnetic states, and (c) the selection of key parameters, such as the magnetic local time step size for drift shell determination. The rigorous development and documentation of this algorithm additionally acts as preliminary step for future thorough reassessment of in situ phase‐space density results using alternative magnetic field models.

Radiation MHD Simulations of Soft X-Ray Emitting Regions in Changing Look AGN

Strong soft X-ray emission called soft X-ray excess is often observed in luminous active galactic nuclei (AGN). It has been suggested that the soft X-rays are emitted from a warm (T = 106 ∼ 107 K) region that is optically thick for the Thomson scattering (warm Comptonization region). Motivated by the recent observations that soft X-ray excess appears in changing look AGN (CLAGN) during the state transition from a dim state without broad emission lines to a bright state with broad emission lines, we performed global three-dimensional radiation magnetohydrodynamic simulations, assuming that the mass accretion rate increases and becomes around 10% of the Eddington accretion rate. The simulation successfully reproduces a warm, Thomson-thick region outside the hot radiatively inefficient accretion flow near the black hole. The warm region is formed by efficient radiative cooling due to inverse Compton scattering. The calculated luminosity 0.01−0.08 L Edd is consistent with the luminosity of CLAGN. We also found that the warm Comptonization region is well described by the steady model of magnetized disks supported by azimuthal magnetic fields. When the antiparallel azimuthal magnetic fields supporting the radiatively cooled region reconnect around the equatorial plane of the disk, the temperature of the region becomes higher by releasing the magnetic energy transported to the region.

Role of magnetic fields in disc galaxies: spiral arm instability

Regularly spaced star-forming regions along the spiral arms of nearby galaxies provide insight into the early stages and initial conditions of star formation. The regular separation of these star-forming regions suggests spiral arm instability as their origin. We explore the effects of magnetic fields on the spiral arm instability. We use 3D global magnetohydrodynamical simulations of isolated spiral galaxies, comparing three different initial plasma beta values (ratios of the thermal to magnetic pressure) of $ $50$, and $10$. We perform a Fourier analysis to calculate the separation of the over-dense regions that formed as a result of the spiral instability. We then compare the separations with observations. We find that the spiral arms in the hydro case ($ are unstable. The fragments are initially connected by gas streams that are reminiscent of the Kelvin-Helmholtz instability. For $ = 50$, the spiral arms also fragment, but the fragments separate earlier and tend to be slightly elongated in the direction perpendicular to the spiral arms. However, in the $ = 10$ run, the arms are stabilised against fragmentation by magnetic pressure. Despite the difference in the initial magnetic field strengths of the $ = 50$ and $10$ runs, the magnetic field is amplified to $ arm 1$ inside the spiral arms for both runs. The spiral arms in the unstable cases (hydro and $ fragment into regularly spaced over-dense regions. We determine their separation to be $ 0.5$ kpc in the hydro and $ 0.65$ kpc in the $ = 50$ case. These two values agree with the observed values found in nearby galaxies. We find a smaller median characteristic wavelength of the over-densities along the spiral arms of $ kpc in the hydro case compared to $0.98^ $ kpc in the $ = 50$ case. Moreover, we find a higher growth rate of the over-densities in the $ = 50$ run compared to the hydro run. We observe magnetic hills and valleys along the fragmented arms in the $ = 50$ run, which is characteristic of the Parker instability.

Investigating the effects of numerical algorithms on global magnetohydrodynamics simulations of solar wind in the inner heliosphere

This paper explores the effects of numerical algorithms on global magnetohydrodynamics simulations of solar wind (SW) in the inner heliosphere. To do so, we use sunRunner3D, a 3-D magnetohydrodynamics model that employs the boundary conditions generated by CORHEL and the PLUTO code to compute the plasma properties of the SW with the ideal magnetohydrodynamics approximation up to 1.1 AU in the inner heliosphere. Mainly, we define three different combinations of numerical algorithms based on their diffusion level. This diffusion level is related to the way of solving the magnetohydrodynamics equations using the finite volume formulation, and, therefore, we set in terms of the divergence-free condition methods, Riemann solvers, variable reconstruction schemes, limiters, and time-steeping algorithms. According to the simulation results, we demonstrate that sunRunner3D reproduces global features of Corotating Interaction Regions observed by Earth-based spacecraft for a set of Carrington rotations that cover a period that lays in the late declining phase of solar cycle 24, independently of the numerical algorithms. Moreover, statistical analyses between models and in-situ measurements show reasonable agreement with the observations, and remarkably, the high diffusive method matches better with in-situ data than low diffusive methods.

Magnetochronology of solar-type star dynamos

In this study, we analyse the magnetic field properties of a set of 15 global magnetohydrodynamic (MHD) simulations of solar-type star dynamos conducted using the ASH code. Our objective is to enhance our understanding of these properties by comparing theoretical results to current observations, and to finally provide fresh insights into the field. We analysed the rotational and magnetic properties as a function of various stellar parameters (mass, age, and rotation rate) in a ‘Sun in time’ approach in our extended set of 3D MHD simulations. To facilitate direct comparisons with stellar magnetism observations using various Zeeman-effect techniques, we decomposed the numerical data into vectorial spherical harmonics. A comparison of the trends we find in our simulations set reveals a promising overall agreement with the observational context of stellar magnetism, enabling us to suggest a plausible scenario for the magneto-rotational evolution of solar-type stars. In particular, we find that the magnetic field may reach a minimum amplitude at a transition value of the Rossby number near unity. This may have important consequences on the long-term evolution of solar-type stars, by impacting the relation between stellar age, rotation, and magnetism. This supports the need for future observational campaigns, especially for stars in the high Rossby number regime.

Intensification of the Electron Zebra Stripes in the Earth's Inner Magnetosphere During Geomagnetic Storms

AbstractWe examined rapid variations in the electron zebra stripe patterns, specifically at L = 1.5, over a three‐month duration, using twin Van Allen Probes within Earth's inner magnetosphere. During geomagnetically quiet intervals, these stripes exhibit a peak‐to‐valley ratio (Δj) ∼1.25 in detrended electron fluxes. However, during geomagnetic storms, they became highly prominent, with Δj > 2.5. The correlation between Δj and net field‐aligned currents (FACs) is observed to be high (0.70). Global magnetohydrodynamic (MHD) simulation results indicate that the westward electric field at midnight at low latitudes in the deep inner magnetosphere correlates well with net FACs. An increase in net FACs could amplify the dawn‐to‐dusk electric field in the deep inner magnetosphere, thereby causing the inward transport of electrons. Given that FACs are linked to the interaction between solar wind and the magnetosphere, our findings emphasize the importance of solar wind‐magnetosphere coupling in the deeper regions of the inner magnetosphere.

A 3D Parametric Venusian Bow Shock Model with the Effects of Mach Number and Interplanetary Magnetic Field

Using global magnetohydrodynamics simulations, we have developed a three-dimensional parametric model for the Venusian bow shock based on a generalized conic section function defined by six parameters, with the effects of the solar wind magnetosonic Mach number (M MS) and the interplanetary magnetic field (IMF) involved. The parametric model’s results reveal the following findings: (1) The size of the Venusian bow shock is primarily determined by M MS. An increase in M MS results in the bow shock moving closer to Venus and a reduction in its flaring angle. (2) Both the subsolar standoff distance and the bow shock’s flaring angle increase with the strength of the IMF components that are perpendicular to the solar wind flow direction (B Y and B Z in the Venus-centered solar orbital coordinate system), whereas the parallel IMF component (B X ) has a limited impact on the subsolar standoff distance but affects the flaring angle. (3) The cross section of the bow shock is elongated in the direction perpendicular to the IMF on the Y–Z plane, and the elongation degree is enhanced with increasing intensities of B Y and B Z . (4) The quasi-parallel bow shock locates closer to the planet as compared to the quasi-perpendicular bow shock. These findings are in alignment with prior empirical and theoretical models. The influences of M MS and IMF on the bow shock’s position and geometry are attributed to the propagation of fast magnetosonic waves, showing the nature of the formation of a collisionless bow shock under the interaction of magnetized flow with an atmospheric object.

Magnetopause properties at the dusk magnetospheric flank from global magnetohydrodynamic simulations, the kinetic Vlasov equilibrium, and in situ observations − Potential implications for SMILE

Magnetopause properties at the dusk magnetospheric flank from global magnetohydrodynamic simulations, the kinetic Vlasov equilibrium, and in situ observations − Potential implications for SMILE

Machine Learning‐Based Emulator for the Physics‐Based Simulation of Auroral Current System

AbstractUsing a machine learning technique called echo state network (ESN), we have developed an emulator to model the physics‐based global magnetohydrodynamic simulation results of REPPU (REProduce Plasma Universe) code. The inputs are the solar wind time series with date and time, and the outputs are the time series of the ionospheric auroral current system in the form of two‐dimensional (2D) patterns of field‐aligned current, potential, and conductivity. We mediated a principal component analysis for a dimensionality reduction of the 2D map time series. In this study, we report the latest upgraded Surrogate Model for REPPU Auroral Ionosphere version 2 (SMRAI2) with significantly improved resolutions in time and space (5 min in time, ∼1° in latitude, and 4.5° in longitude), where the dipole tilt angle is also newly added as one of the input parameters to reproduce the seasonal dependence. The fundamental dependencies of the steady‐state potential and field‐aligned current patterns on the interplanetary magnetic field directions are consistent with those obtained from empirical models. Further, we show that the ESN‐based emulator can output the AE index so that we can evaluate the performance of the dynamically changing results, comparing with the observed AE index. Since the ESN‐based emulator runs a million times faster than the REPPU simulation, it is promising that we can utilize the emulator for the real‐time space weather forecast of the auroral current system as well as to obtain large‐number ensembles to achieve future data assimilation‐based forecast.

North–South Plasma Asymmetry Across Mercury's Near‐Tail Current Sheet

AbstractAmong nearly 300 near‐Mercury tail current sheet crossings performed by the MESSENGER spacecraft, we identified 37 traversals of an asymmetric current sheet, wherein the lobe densities on opposite sides differ by a factor of three or more. These asymmetric current sheet crossings primarily occur on the dawnside. A global magnetohydrodynamic (MHD) simulation was found to be in excellent agreement with the observations. The results suggest that the north–south density asymmetry is caused by solar wind entering via an upstream‐connected window in one hemisphere. Furthermore, the Parker spiral interplanetary magnetic field (IMF) controls the near‐tail density asymmetries, whereas Mercury's offset dipole magnetic field controls those in mid‐ or distant‐tail regions. We propose that hemispheric asymmetries in Mercury's magnetospheric convection occur under strong IMF conditions.

Synchrotron Polarization Signatures of Surface Waves in Supermassive Black Hole Jets

Supermassive black holes in active galactic nuclei are known to launch relativistic jets, which are observed across the entire electromagnetic spectrum and thought to be efficient particle accelerators. Their primary radiation mechanism for radio emission is polarized synchrotron emission produced by a population of nonthermal electrons. In this Letter, we present a global general relativistic magnetohydrodynamical (GRMHD) simulation of a magnetically arrested disk (MAD). After the simulation reaches the MAD state, we show that waves are continuously launched from the vicinity of the black hole and propagate along the interface between the jet and the wind. At this interface, a steep gradient in velocity is present between the mildly relativistic wind and the highly relativistic jet. The interface is, therefore, a shear layer, and due to the shear, the waves generate roll-ups that alter the magnetic field configuration and the shear layer geometry. We then perform polarized radiation transfer calculations of our GRMHD simulation and find signatures of the waves in both total intensity and linear polarization, effectively lowering the fully resolved polarization fraction. The telltale polarization signatures of the waves could be observable by future very long baseline interferometric observations, e.g., the next-generation Event Horizon Telescope.

Identifying and Visualizing Terrestrial Magnetospheric Topology using Geodesic Level Set Method

AbstractThis study introduces a novel numerical method for identifying and visualizing the terrestrial magnetic field topology in a large‐scale three‐dimensional global MHD (Magneto‐Hydro‐Dynamic) simulation. The (un)stable two‐dimensional manifolds are generated from critical points (CPs) located north and south of the magnetosphere using an improved geodesic level set method. A boundary value problem is solved numerically using a shooting method to forward a new geodesic level set from the previous set. These sets are generated starting from a small circle whose centre is a CP. The level sets are the sets of mesh points that form the magnetic manifold, which determines the magnetic field topology. In this study, a consistent method is proposed to determine the magnetospheric topology. Using this scheme, we successfully visualize a terrestrial magnetospheric field topology and identify its two neutral lines using the global MHD simulation. Our results present a terrestrial topology that agrees well with the recent magnetospheric physics and can help us understand various nonlinear magnetospheric dynamics and phenomena. Our visualization enables us to fill the gaps between current magnetospheric physics that can be observed via satellites and nonlinear dynamics, particularly, the bifurcation theory, in the future.

Nighttime Geomagnetic Response to Jumps of Solar Wind Dynamic Pressure: A Possible Cause of Québec Blackout in March 1989

AbstractBy performing a global magnetohydrodynamic (MHD) simulation, we investigated magnetic disturbances on the ground at high‐latitudes in response to jumps in the solar wind dynamic pressure, namely a sudden commencement (SC). After the arrival of the jump, a pair of field‐aligned currents (FACs), related to the preliminary impulse, develop and travel in the anti‐sunward direction. Soon after another pair related to the main impulse (MI) appears and travels in the anti‐sunward direction. The horizontal ionospheric current associated with the MI remains strong when propagating to the nightside. On the dawnside the MI current flows sunward (anti‐sunward) resulting in northward (southward) ground magnetic disturbance at higher (lower) latitude in the post‐midnight sector. These features are similar to those observed in Canada in the high‐latitude post‐midnight sector when the Québec blackout took place on 13 March 1989. The nighttime geomagnetic perturbations associated with the MI occur regardless of the magnitude of the solar wind dynamic pressure and IMF orientation. The amplitude of the geoelectric field, which is closely related to the geomagnetically induced currents (GICs), reaches the maximum value just before and around the maximum of the southward magnetic disturbance. This is consistent with the moment at which the blackout occurred during the southward magnetic perturbation. We suggest that the blackout in Québec could be caused by the MI‐associated Hall current passing over the Hydro‐Québec power system on the nightside. The nighttime polar region is shown to be sensitive to hazardous GICs for large‐amplitude jumps in the solar wind dynamic pressure.

The complexity of the day-side X-line during southward interplanetary magnetic field

High-resolution global magnetohydrodynamics (MHD) simulations include both meso- and global-scale processes occurring at the magnetopause, which interact to determine the time-dependent orientation of the day-side x-line (DXL). This study demonstrates that the global orientation of the DXL in GAMERA global MHD simulations varies on a time scale of minutes during steady southward interplanetary magnetic field conditions. This behavior manifests in observational data when reconnection outflows indicate that the direction to the x-line is opposite to the prediction from a steady-state model of the reconnection location. Because steady-state models of the DXL do not capture dynamics that are independent of solar wind variations, particularly surface waves and flux transfer events, they represent a time-averaged state of the system.