Astrophysical Systems Research Articles

Three-dimensional magnetic reconnection in complex multiple X-point configurations in an ancient solar–lunar terrestrial system

Magnetic reconnection processes in three-dimensional (3D) complex field configurations have been investigated in different magneto-plasma systems in space, laboratory, and astrophysical systems. Two-dimensional (2D) features of magnetic reconnection have been well developed and applied successfully to systems with symmetrical property, such as toroidal fusion plasmas and laboratory experiments with an axial symmetry. But in asymmetric systems, the 3D features are inevitably different from those in the 2D case. Magnetic reconnection structures in multiple celestial body systems, particularly star–planet–Moon systems, bring fresh insights to the understanding of the 3D geometry of reconnection. Thus, we take magnetic reconnection in an ancient solar–lunar terrestrial magneto-plasma system as an example by using its crucial parameters approximately estimated already and also some specific applications in pathways for energy and matter transports among Earth, ancient Moon, and the interplanetary magnetic field (IMF). Then, magnetic reconnection of the ancient lunar–terrestrial magnetospheres with the IMF is investigated numerically in this work. In a 3D simulation for the Earth–Moon–IMF system, topological features of complex magnetic reconnection configurations and dynamical characteristics of magnetic reconnection processes are studied. It is found that a coupled lunar-terrestrial magnetosphere is formed, and under various IMF orientations, multiple X-points emerge at distinct locations, showing three typical magnetic reconnection structures in such a geometry, i.e., the X-line, the triple current sheets, and the A–B null pairs. The results can conduce to further understanding of reconnection physics in 3D for plasmas in complex magnetic configurations, and also a possible mechanism for energy and matters transport in evolutions of similar astrophysical systems.

Irregular Proton Injection to High Energies at Interplanetary Shocks

How thermal particles are accelerated to suprathermal energies is an unsolved issue, crucial for many astrophysical systems. We report novel observations of irregular, dispersive enhancements of the suprathermal particle population upstream of a high-Mach-number interplanetary shock. We interpret the observed behavior as irregular “injections” of suprathermal particles resulting from shock front irregularities. Our findings, directly compared to self-consistent simulation results, provide important insights for the study of remote astrophysical systems where shock structuring is often neglected.

Ion acoustic soliton with thermal ions and non-thermal electrons in a high-relativistic electron-positron-ion plasma

In this investigation, both compressive and rarefactive solitons are shown to exist in electron-positron-ion plasma consisting of thermal positrons, non-thermal electrons, and high relativistic thermal ions. The reductive perturbation method (RPM) is used to derive the nonlinear Korteweg-de Vries (KdV) equation from the governing normalized basic set of equations. It is shown that only the fast ion-acoustic mode may produce the small amplitude rarefactive and compressive KdV solitons. It has been observed that the addition of variable ion species temperature dramatically alters the fundamental characteristics of amplitude and width. Moreover, the amplitude decreases as the ion to electron temperature ratio rises. The soliton's amplitude is also increased by the relativistic factor. High-relativistic electron-positron-ion plasmas are found in extreme astrophysical environments like pulsar magnetospheres, gamma-ray bursts, and active galactic nuclei. These plasmas’ wave and soliton activity can shed light on the dynamics of these astrophysical systems.

Assessment of a new sub-grid model for magnetohydrodynamical turbulence – II. Kelvin–Helmholtz instability

ABSTRACT The modelling of astrophysical systems such as binary neutron star mergers or the formation of magnetars from the collapse of massive stars involves the numerical evolution of magnetized fluids at extremely large Reynolds numbers. This is a major challenge for (unresolved) direct numerical simulations which may struggle to resolve highly dynamical features as, e.g. turbulence, magnetic field amplification, or the transport of angular momentum. Sub-grid models offer a means to overcome those difficulties. In a recent paper we presented MInIT, an MHD-instability-induced-turbulence mean-field, sub-grid model based on the modelling of the turbulent (Maxwell, Reynolds, and Faraday) stress tensors. While in our previous work MInIT was assessed within the framework of the magnetorotational instability, in this paper we further evaluate the model in the context of the Kelvin–Helmholtz instability (KHI). The main difference with other sub-grid models (as e.g. the alpha-viscosity model or the gradient model) is that in MInIT, we track independently the turbulent energy density at sub-grid scales, which is used, via a simple closure relation, to compute the different turbulent stresses relevant for the dynamics. The free coefficients of the model are calibrated using well-resolved box simulations of magnetic turbulence generated by the KHI. We test the model against these simulations and show that it yields order-of-magnitude accurate predictions for the evolution of the turbulent Reynolds and Maxwell stresses.

First Principles Description of Plasma Expansion Using the Expanding Box Model

Multi-scale modeling of expanding plasmas is crucial for understanding the dynamics and evolution of various astrophysical plasma systems such as the solar and stellar winds. In this context, the Expanding Box Model (EBM) provides a valuable framework to mimic plasma expansion in a non-inertial reference frame, co-moving with the expansion but in a box with a fixed volume, which is especially useful for numerical simulations. Here, fundamentally based on the Vlasov equation for magnetized plasmas and the EBM formalism for coordinates transformations, for the first time, we develop a first principles description of radially expanding plasmas in the EB frame. From this approach, we aim to fill the gap between simulations and theory at microscopic scales to model plasma expansion at the kinetic level. Our results show that expansion introduces non-trivial changes in the Vlasov equation (in the EB frame), especially affecting its conservative form through non-inertial forces purely related to the expansion. In order to test the consistency of the equations, we also provide integral moments of the modified Vlasov equation, obtaining the related expanding moments (i.e., continuity, momentum, and energy equations). Comparing our results with the literature, we obtain the same fluids equations (ideal-MHD), but starting from a first principles approach. We also obtained the tensorial form of the energy/pressure equation in the EB frame. These results show the consistency between the kinetic and MHD descriptions. Thus, the expanding Vlasov kinetic theory provides a novel framework to explore plasma physics at both micro and macroscopic scales in complex astrophysical scenarios.

Simulations of weakly magnetized turbulent mixing layers

ABSTRACT Radiative turbulent mixing layers (TMLs) are expected to form pervasively at the phase boundaries in multiphase astrophysical systems. This inherently small-scale structure is dynamically crucial because it directly regulates the mass, momentum, and energy exchanges between adjacent phases. Previous studies on hydrodynamic TMLs have revealed the interactions between cold and hot phases in the context of the circumgalactic medium, offering important insight into the fate of cold clouds traveling through hot galactic winds. However, the role of magnetic field has only been sparsely investigated. We perform a series of 3D magnetohydrodynamics simulations of such mixing layers in the presence of weak to modest background magnetic field. We find that due to field amplification, even relatively weak background magnetic fields can significantly reduce the surface brightness and inflow velocity of the hot gas in the mixing layer. This reduction is attributed to a combination of magnetic pressure support and direct suppression of turbulent mixing, both of which alter the phase structures. Our results are largely independent of thermal conduction and converged with resolution, offering insights on the survival of cold gas in multiphase systems.

The Origin of Power-law Spectra in Relativistic Magnetic Reconnection

Magnetic reconnection is often invoked as a source of high-energy particles, and in relativistic astrophysical systems it is regarded as a prime candidate for powering fast and bright flares. We present a novel analytical model—supported and benchmarked with large-scale three-dimensional kinetic particle-in-cell simulations in electron–positron plasmas—that elucidates the physics governing the generation of power-law energy spectra in relativistic reconnection. Particles with Lorentz factor γ ≳ 3σ (here, σ is the magnetization) gain most of their energy in the inflow region, while meandering between the two sides of the reconnection layer. Their acceleration time is , where η rec ≃ 0.06 is the inflow speed in units of the speed of light and ω c = eB 0/mc is the gyrofrequency in the upstream magnetic field. They leave the region of active energization after t esc, when they get captured by one of the outflowing flux ropes of reconnected plasma. We directly measure t esc in our simulations and find that t esc ∼ t acc for σ ≳ few. This leads to a universal (i.e., σ-independent) power-law spectrum for the particles undergoing active acceleration, and for the overall particle population. Our results help to shed light on the ubiquitous presence of power-law particle and photon spectra in astrophysical nonthermal sources.

Hydrodynamical simulations of circumbinary accretion: balance between heating and cooling

ABSTRACT Hydrodynamical interaction in circumbinary discs (CBDs) plays a crucial role in various astrophysical systems, ranging from young stellar binaries to supermassive black hole binaries in galactic centres. Most previous simulations of binary-disc systems have adopted locally isothermal equation of state. In this study, we use the grid-based code Athena++ to conduct a suite of two-dimensional viscous hydrodynamical simulations of circumbinary accretion on a Cartesian grid, resolving the central cavity of the binary. The gas thermodynamics is treated by thermal relaxation towards an equilibrium temperature (based on the constant − β cooling ansatz, where β is the cooling time in units of the local Keplerian time). Focusing on equal mass, circular binaries in CBDs with (equilibrium) disc aspect ratio H/R = 0.1, we find that the cooling of the disc gas significantly influences the binary orbital evolution, accretion variability, and CBD morphology, and the effect depends sensitively on the disc viscosity prescriptions. When adopting a constant kinematic viscosity, a finite cooling time (β ≳ 0.1) leads to a binary inspiral as opposed to an outspiral and the CBD cavity becomes more symmetric. When adopting a dynamically varying α-viscosity, binary inspiral only occurs within a narrow range of cooling time (corresponding to β around 0.5).

Statistics of pressure fluctuations in turbulent kinetic plasmas

ABSTRACTIn this study, we explore the statistics of pressure fluctuations in kinetic collisionless turbulence. A 2.5D kinetic particle-in-cell simulation of decaying turbulence is used to investigate pressure balance via the evolution of thermal and magnetic pressure in a plasma with β of order unity. We also discuss the behaviour of thermal, magnetic, and total pressure structure functions and their corresponding wavenumber spectra. The total pressure spectrum exhibits a slope of −7/3 extending for about a decade in the ion-inertial range. In contrast, shallower −5/3 spectra are characteristic of the magnetic pressure and thermal pressure. The steeper total pressure spectrum is a consequence of cancellation caused by density-magnetic field magnitude anti-correlation. Further, we evaluate higher order total pressure structure functions in an effort to discuss intermittency and compare the power exponents with higher order structure functions of velocity and magnetic fluctuations. Finally, applications to astrophysical systems are also discussed.

Magnetohydrodynamic Shocks Revisited: Magnetically Constraining the Upstream Solar Wind Condition

Magnetohydrodynamic (MHD) shocks are commonly found in various space and astrophysical systems. One often encounters the problem that only the magnetic field data are available when studying interplanetary shocks and planetary bow shocks in the solar system, and not the plasma data due to instrumental and/or telemetry limitations. We aim to reveal the mathematical structure of the MHD Rankine–Hugoniot relation under the condition that magnetic field changes across the shock are known, which can be devised into a diagnosis tool to constrain the shock parameters. The perturbative solution is already known for the Rankine–Hugoniot equation. The solution is combined with the condition of the de Hoffmann–Teller frame for the vanishing convective electric field, and we obtain a two-dimensional matrix equation for the density jump and the Alfvén Mach number. The lesson is that the shock parameters (density jump and Alfvén Mach number) are analytically obtained by inverting the MHD Rankine–Hugoniot equation under the condition that the plasma beta is either a priori set or known. Further numerical studies indicate that the presented shock jump estimation works particularly well for quasi-perpendicular shocks. Moreover, a test case with the Earth bow shock crossing data supports our theoretical development. The presented method opens the door to determining or constraining the shock parameters even if the data are limited to the magnetic field only.

Alpha clustering in nuclear astrophysics and topology

When we think of clustering in nuclear physics, the astrophysical importance within light nuclei and structural manifestations with classical analogs immediately come to mind. 4He, also known as the alpha particle, is the most abundant nucleus in the Universe, being quite tightly bound for its mass, with a first excited state of over 20 MeV. The nature of the alpha particle places it in a unique position within nuclear astrophysics and structure (including geometry). The plurality of energy release from stellar hydrogen fusion—whether quiescent or explosive—comes from the conversion of hydrogen to helium. Within more complex nuclei, the alpha particles are continuously arranged, leading to fascinating phenomena such as excited rotational bands, Borromean ring ground states, and linear structures. Nuclei with an equal and even number of protons and neutrons are colloquially referred to as “alpha conjugate nuclei,” where such special properties are the most pronounced and easiest to spot. However, when a single nucleon or a pair of nucleons is added to the system, alpha clustering not only remains evident but it may also be enhanced. Excited states with large alpha partial widths are a signature of clustering behavior, and these states can have a profound effect on the reaction rates in astrophysical systems when the excitation energy aligns with the so-called Gamow energy—the preferential thermal energy to statistically overcome the Coulomb barrier. In this article, we will consider in detail the specific ramifications of alpha clustering in selected scenarios for both nuclear astrophysics and topology. In particular, we discussed the astrophysical reactions of 7Li (α, γ), 7Be+α, 11C (α, p), and 30S (α, p), where α-clusters may increase the reaction rates from 10% to an order of magnitude; large α resonances make the astrophysical rate of 18F (p, α) quite uncertain. We also focused on the α rotational bands of both positive and negative parities of 11B and 11C, and finally on the strongest evidence for the linear-chain cluster state observed in 14C.

Simple Convective Accretion Flows (SCAFs): Explaining the ≈ −1 Density Scaling of Hot Accretion Flows around Compact Accretors

Recent simulations find that hot gas accretion onto compact accretors is often highly turbulent and diskless, and shows a power-law density profile with slope α ρ ≈ −1. These results are consistent with observational constraints, but do not match existing self-similar solutions of radiatively inefficient accretion flows. We develop a theory for this new class of accretion flows, which we dub simple convective accretion flows (SCAFs). We use a set of hydrodynamic simulations to provide a minimalistic example of SCAFs, and develop an analytic theory to explain and predict key flow properties. We demonstrate that the turbulence in the flow is driven locally by convection, and argue that radial momentum balance, together with an approximate up–down symmetry of convective turbulence, yields α ρ = −1 ± few 0.1. Empirically, for an adiabatic hydrodynamic flow with γ ≈ 5/3, we get α ρ ≈ −0.8; the resulting accretion rate (relative to the Bondi accretion rate), , agrees very well with the observed accretion rates in Sgr A*, M87*, and a number of wind-fed supergiant X-ray binaries. We also argue that the properties of SCAFs are relatively insensitive to additional physical ingredients such as cooling and magnetic field; this explains their common appearance across simulations of different astrophysical systems.

A Model for Nonthermal Particle Acceleration in Relativistic Magnetic Reconnection

The past decade has seen an outstanding development of nonthermal particle acceleration in magnetic reconnection in magnetically dominated systems, with clear signatures of power-law energy distributions as a common outcome of first-principles kinetic simulations. Here we propose a semianalytical model for systematically investigating nonthermal particle acceleration in reconnection. We show particle energy distributions are well determined by particle injection, acceleration, and escape processes. Using a series of kinetic simulations, we accurately evaluate the energy- and time-dependent model coefficients. The resulting spectral characteristics, including the spectral index and lower and upper bounds of the power-law distribution, agree well with the simulation results. Finally, we apply the model to predict the power-law indices and break energies in astrophysical reconnection systems.

Expansion free spherical anisotropic solutions

We investigate the behavior of expansion free collapsing fluids, as studied by L. Herrera, A. Di Prisco and J. Ospino [Symmetry 15 (2023) 754], in the framework of [Formula: see text] gravity, which represents a modification of Einstein’s general relativity by establishing a function of the Ricci scalar [Formula: see text] in the gravitational action. We explore dynamical equations from Bianchi identities that demonstrate the motion and evolution of physical systems under the influence of gravitational fields. We match the inner and outer geometries of spacetime on the hypersurface to develop junction conditions by using the Misner–Sharp formalism. This allows us to identify the connection between mass functions for the inner and outer space as well as the relationship for heat flux [Formula: see text] and radial pressure [Formula: see text]. We also investigate analytical solutions of dissipative fluid distribution that fulfill the vanishing expansion condition together with the vanishing complexity factor constraint. For this, we introduce new constraints that permit the integration of the complex system in [Formula: see text] gravity. Next, we extract a set of differential equations that explain the dynamical structure of the dissipative spheres both in geodesic and non-geodesic fluids. Furthermore, we explore the physical characteristics of the obtained solutions, such as heat flux, energy density, shear stress, fluid’s temperature along with tangential and radial pressure, to assess their viability in describing real astrophysical systems.

Controlling secondary flows in Taylor–Couette flow using axially spaced superhydrophobic surfaces

Turbulent shear flows are abundant in geophysical and astrophysical systems and in engineering-technology applications. They are often riddled with large-scale secondary flows that drastically modify the characteristics of the primary stream, preventing or enhancing mixing, mass and heat transfer. Using experiments and numerical simulations, we study the possibility of modifying these secondary flows by using superhydrophobic surface treatments that reduce the local shear. We focus on the canonical problem of Taylor–Couette flow, the flow between two coaxial and independently rotating cylinders, which has robust secondary structures called Taylor rolls that persist even at significant levels of turbulence. We generate these structures by rotating only the inner cylinder of the system, and show that an axially spaced superhydrophobic treatment can weaken the rolls through a mismatching surface heterogeneity, as long as the roll size can be fixed. The minimum hydrophobicity of the treatment required for this flow control is rationalized, and its effectiveness beyond the Reynolds numbers studied here is also discussed.

Three-dimensional modelling of the shock–turbulence interaction

ABSTRACTThe complex interaction between shocks and plasma turbulence is extremely important to address crucial features of energy conversion in a broad range of astrophysical systems. We study the interaction between a supercritical, perpendicular shock and pre-existing, fully developed plasma turbulence, employing a novel combination of magnetohydrodynamic and small-scale, hybrid-kinetic simulations where a shock is propagating through a turbulent medium. The variability of the shock front in the unperturbed case and for two levels of upstream fluctuations is addressed. We find that the behaviour of shock ripples, i.e. shock surface fluctuations with short (a few ion skin depths, di) wavelengths, is modified by the presence of pre-existing turbulence, which also induces strong corrugations of the shock front at larger scales. We link this complex behaviour of the shock front and the shock downstream structuring with the proton temperature anisotropies produced in the shock–turbulence system. Finally, we put our modelling effort in the context of spacecraft observations, elucidating the role of novel cross-scale, multispacecraft measurements in resolving shock front irregularities at different scales. These results are relevant for a broad range of astrophysical systems characterized by the presence of shock waves interacting with plasma turbulence.

Spectrum of kinetic-Alfvén-wave turbulence: intermittency or tearing mediation?

ABSTRACT We investigate the spectral properties of the electromagnetic fluctuations of sub-ion scale turbulence in weakly collisional, low-beta plasmas using a two-field isothermal gyrofluid model. The numerical results strongly support a description of the turbulence as a critically balanced Kolmogorov-like cascade of kinetic Alfvén wave fluctuations, as amended by previous studies to include intermittency effects. The measured universal index of the energy spectra from systems with different flux-unfreezing mechanisms excludes the role of tearing mediation in determining the spectra. The fluctuations remain isotropic in the plane perpendicular to the strong background magnetic fields as they cascade to smaller scales, which explains the absence of tearing mediation. The calculation of high-order, multipoint structure functions of magnetic fluctuations suggests that the intermittent structures have a quasi-2D, sheet-type morphology. These results are useful for explaining recent observations of the spectrum and structure of magnetic and density fluctuations in the solar wind at sub-proton scales, and are relevant for modelling the energy dissipation in a broad range of astrophysical systems.

A Generalized Double Chaplygin Model for Anisotropic Matter: The Newtonian Case

In this work, we investigate astrophysical systems in a Newtonian regime using anisotropic matter. For this purpose, we considered that both radial and tangential pressures satisfy a generalized Chaplygin-type equation of state. Using this model, we found the Lane–Emden equation for this system and solved it numerically for several sets of parameters. Finally, we explored the mass supported by this physical system and compared it with the Chandrasekhar mass.

Acceleration of cosmic rays in presence of magnetohydrodynamic fluctuations at small scales

ABSTRACT This work investigates the evolution of the distribution of charged particles (cosmic rays) due to the mechanism of stochastic turbulent acceleration (STA) in presence of small-scale turbulence with a mean magnetic field. STA is usually modelled as a biased random walk process in the momentum space of the non-thermal particles. This results in an advection-diffusion type transport equation for the non-thermal particle distribution function. Under quasi-linear approximation, and by assuming turbulent spectra with power being available only in the sub-gyroscale range, we find that the Fokker–Planck diffusion coefficients Dγγ and Dμμ scale with the Lorentz factor γ as Dγγ ∝ γ−2/3 and Dμμ ∝ γ−8/3. We consider Alfvèn and fast waves in our calculations, and find a universal trend for the momentum diffusion coefficient irrespective of the properties of the small-scale turbulence. Such universality has already been reported regarding the spatial diffusion of the cosmic rays, and, here too, we observe a universality in the momentum diffusion coefficient. Furthermore, with the calculated transport coefficients, we numerically solve the advection-diffusion-type transport equation for the non-thermal particles. We demonstrate the interplay of various mircophysical processes such as STA, synchrotron loss, and particle escape on the particle distribution by systematically varying the parameters of the problem. We observe that the effect of the small-scale turbulence is more impactful for the high-energy protons as compared to the electrons and such turbulence is capable of sustaining the energy of the protons from catastrophic radiative loss processes. Such a finding is novel and helps us to enhance our understanding about the hadronic emission processes that are typically considered as a competitor for the leptonic emission for certain astrophysical systems.

Inertial Alfvén waves with polarization force in a dusty plasma

The properties of dust inertial Alfvén waves (DIAWs) in a dusty plasma with the effects of polarization force are investigated. First, the dispersion equation is obtained from the fluid equations by considering the polarization force. Then, the energy balance equation is deduced by applying the Sagdeev pseudopotential method. It is shown that the polarization force affects the characteristics of DIAWs as well as the pseudopotential significantly. In addition, the existence of the polarization force changes the width of DIAWs. The present results may be helpful in further understanding the effect of polarization forces on the formation of linear and nonlinear structures of Alfvén waves in the laboratory plasma, astrophysical dusty plasma systems, planetary ring systems, cometary tails, etc.