skip to main content


The NSF Public Access Repository (NSF-PAR) system and access will be unavailable from 11:00 PM ET on Thursday, May 23 until 2:00 AM ET on Friday, May 24 due to maintenance. We apologize for the inconvenience.

Title: A New Covariant Formalism for Kinetic Plasma Simulations in Curved Spacetimes

Low-density plasmas are characterized by a large-scale separation between the gyromotion of particles around local magnetic fields and the macroscopic scales of the system, often making global kinetic simulations computationally intractable. The guiding center formalism has been proposed as a powerful tool to bridge the gap between these scales. Despite its usefulness, the guiding center approach has been formulated successfully only in flat spacetimes, limiting its applicability in astrophysical settings. Here, we present a new covariant formalism that leads to kinetic equations in the guiding center limit that are valid in arbitrary spacetimes. Through a variety of experiments, we demonstrate that our equations capture all known gyrocenter drifts while overcoming one severe limitation imposed on numerical algorithms by the fast timescales of the particle gyromotion. This formalism will enable explorations of a variety of global plasma kinetic phenomena in the curved spacetimes around black holes and neutron stars.

more » « less
Author(s) / Creator(s):
; ; ; ;
Publisher / Repository:
DOI PREFIX: 10.3847
Date Published:
Journal Name:
The Astrophysical Journal Letters
Medium: X Size: Article No. L6
["Article No. L6"]
Sponsoring Org:
National Science Foundation
More Like this
  1. Over the late Holocene, a variety of hydroclimate-sensitive proxies have identified substantial, multidecadal changes in Indian summer monsoon (ISM) precipitation, the most prominent of which is the “4.2 ka event”. This interval, dated to ~4.2-3.9 ka, is associated with severe droughts across South Asia that are linked to societal change. Given the absence of the 4.2 ka event in polar records, the 4.2 ka event is generally associated with low latitude forcings, but no clear consensus on its origins has been reached. We investigated the ISM response to the 4.2 ka event through analysis of aragonite stalagmites from Siddha cave, formed in the lower Paleozoic Dhading dolomite in the Pokhara Valley of central Nepal (28.0˚N, 84.1˚E; ~850 m.a.s.l.). The climate of this region is dominated by small monthly variations in air temperature (21±5˚C) but strong precipitation seasonality associated with the ISM: ~80% of the annual 3900 mm of rainfall occurs between June and September. High uranium and low detrital thorium abundances in these stalagmites yield precise U/Th ages that all fall within stratigraphic order. These dates reveal continuous growth from 4.30-2.26 ka, interrupted only by a hiatus from 3.27-3.10 ka. Overlap with coeval aragonite stalagmites reveals generally consistent trends in carbon and oxygen isotope ratios, suggesting that these stalagmites reflect environmental variability and not secondary (e.g., kinetic) effects. Many stalagmite-based paleomonsoon reconstructions rely on oxygen isotope ratios, which track amount effects in regional rainfall. However, our on-going rainwater collection and analysis program, as well as a previous study conducted in Kathmandu, 120 km the east of Siddha cave, reveals that amount effects in precipitation are weak in this region, particularly during the monsoon season, and thus we rely instead on carbon isotope ratios, which have been demonstrated to track site-specific effective precipitation. Siddha cave stalagmite carbon isotopes, in contrast to other South Asian proxy records, indicate that ISM rainfall increased at Siddha cave from 4.13-3.91 ka. As a further test of this result, we analyzed uranium abundances in the section spanning 4.3-3.4 ka. Uranium serves as an indicator of prior aragonite precipitation and thus of hydroclimate, and like carbon isotopes, suggests increased ISM rainfall coincident with the 4.2 ka event. This precipitation anomaly is nearly identical in timing and structure but anti-phased with stalagmites from Mawmluh cave, northeastern India. We investigated the climatic origins of this precipitation dipole using observational data from the Global Precipitation Climatology Centre (GPCC) and Hadley Center Sea Ice and Sea Surface Temperature (HadISST) products. Preliminary spatial composites suggest that large precipitation differences between Mawmluh and Siddha caves are associated with SST anomalies in the equatorial Pacific. Additionally, superposed Epoch Analysis shows relatively rapid eastern Indian Ocean cooling during the summer monsoon season coeval with large precipitation differences between these sites. Our findings lend support to a tropical Indo-Pacific origin of the 4.2 ka event. 
    more » « less
  2. Abstract

    We presentGRaM-X(GeneralRelativisticacceleratedMagnetohydrodynamics on AMReX), a new GPU-accelerated dynamical-spacetime general relativistic magnetohydrodynamics (GRMHD) code which extends the GRMHD capability of Einstein Toolkit to GPU-based exascale systems.GRaM-Xsupports 3D adaptive mesh refinement (AMR) on GPUs via a new AMR driver for the Einstein Toolkit calledCarpetXwhich in turn leveragesAMReX, an AMR library developed for use by the United States DOE’s Exascale Computing Project. We use the Z4c formalism to evolve the Einstein equations and the Valencia formulation to evolve the equations of GRMHD.GRaM-Xsupports both analytic as well as tabulated equations of state. We implement TVD and WENO reconstruction methods as well as the HLLE Riemann solver. We test the accuracy of the code using a range of tests on static spacetime, e.g. 1D magnetohydrodynamics shocktubes, the 2D magnetic rotor and a cylindrical explosion, as well as on dynamical spacetimes, i.e. the oscillations of a 3D Tolman-Oppenheimer-Volkhof star. We find excellent agreement with analytic results and results of other codes reported in literature. We also perform scaling tests and find thatGRaM-Xshows a weak scaling efficiency of ∼40%–50% on 2304 nodes (13824 NVIDIA V100 GPUs) with respect to single-node performance on OLCF’s supercomputer Summit.

    more » « less
  3. Abstract

    The dissipative mechanism in weakly collisional plasma is a topic that pervades decades of studies without a consensus solution. We compare several energy dissipation estimates based on energy transfer processes in plasma turbulence and provide justification for the pressure–strain interaction as a direct estimate of the energy dissipation rate. The global and scale-by-scale energy balances are examined in 2.5D and 3D kinetic simulations. We show that the global internal energy increase and the temperature enhancement of each species are directly tracked by the pressure–strain interaction. The incompressive part of the pressure–strain interaction dominates over its compressive part in all simulations considered. The scale-by-scale energy balance is quantified by scale filtered Vlasov–Maxwell equations, a kinetic plasma approach, and the lag dependent von Kármán–Howarth equation, an approach based on fluid models. We find that the energy balance is exactly satisfied across all scales, but the lack of a well-defined inertial range influences the distribution of the energy budget among different terms in the inertial range. Therefore, the widespread use of the Yaglom relation in estimating the dissipation rate is questionable in some cases, especially when the scale separation in the system is not clearly defined. In contrast, the pressure–strain interaction balances exactly the dissipation rate at kinetic scales regardless of the scale separation.

    more » « less
  4. Abstract

    The expanding-box model of the solar wind has been adopted in the literature within the context of magnetohydrodynamics, hybrid, and full particle-in-cell simulations to investigate the dynamic evolution of the solar wind. The present paper extends such a method to the framework of self-consistent quasilinear kinetic theory. It is shown that the expanding-box quasilinear methodology is largely equivalent to the inhomogeneous steady-state quasilinear model discussed earlier in the literature, but a distinction regarding the description of wave dynamics between the two approaches is also found. The expanding-box quasilinear formalism is further extended to include the effects of a spiraling solar-wind magnetic field as well as collisional age effects. The present finding shows that the expanding-box quasilinear approach and the steady-state global-kinetic models may be employed interchangeably in order to address other more complex problems associated with the solar-wind dynamics.

    more » « less
  5. Abstract

    Knowledge silos emerge when structural properties of organizational interaction networks limit the diffusion of information. These structural barriers are known to take many forms at different scales—hubs in otherwise sparse organizations, large dense teams, or global core-periphery structure—but we lack an understanding of how these different structures interact and shape dynamics. Here we take a first theoretical step in bridging the gap between the mathematical literature on localization of spreading dynamics and the more applied literature on knowledge silos in organizational interaction networks. To do so, we introduce a new model that considers a layered structure of teams to unveil a new form of hierarchical localization (i.e. the localization of information at the top or center of an organization) and study its interplay with known phenomena of mesoscopic localization (i.e. the localization of information in large groups),k-core localization (i.e. around denser subgraphs) and hub localization (i.e. around high degree stars). We also include a complex contagion mechanism by considering a general infection kernel which can depend on hierarchical level (influence), degree (popularity), infectious neighbors (social reinforcement) or team size (importance). This very general model allows us to explore the multifaceted phenomenon of information siloing in complex organizational interaction networks and opens the door to new optimization problems to promote or hinder the emergence of different localization regimes.

    more » « less