In the realm of the Hubble tension—a review of solutions
Eleonora di Valentino, O. Mena, S. Pan
et al.
The simplest ΛCDM model provides a good fit to a large span of cosmological data but harbors large areas of phenomenology and ignorance. With the improvement of the number and the accuracy of observations, discrepancies among key cosmological parameters of the model have emerged. The most statistically significant tension is the 4σ to 6σ disagreement between predictions of the Hubble constant, H 0, made by the early time probes in concert with the ‘vanilla’ ΛCDM cosmological model, and a number of late time, model-independent determinations of H 0 from local measurements of distances and redshifts. The high precision and consistency of the data at both ends present strong challenges to the possible solution space and demands a hypothesis with enough rigor to explain multiple observations—whether these invoke new physics, unexpected large-scale structures or multiple, unrelated errors. A thorough review of the problem including a discussion of recent Hubble constant estimates and a summary of the proposed theoretical solutions is presented here. We include more than 1000 references, indicating that the interest in this area has grown considerably just during the last few years. We classify the many proposals to resolve the tension in these categories: early dark energy, late dark energy, dark energy models with 6 degrees of freedom and their extensions, models with extra relativistic degrees of freedom, models with extra interactions, unified cosmologies, modified gravity, inflationary models, modified recombination history, physics of the critical phenomena, and alternative proposals. Some are formally successful, improving the fit to the data in light of their additional degrees of freedom, restoring agreement within 1–2σ between Planck 2018, using the cosmic microwave background power spectra data, baryon acoustic oscillations, Pantheon SN data, and R20, the latest SH0ES Team Riess, et al (2021 Astrophys. J. 908 L6) measurement of the Hubble constant (H 0 = 73.2 ± 1.3 km s−1 Mpc−1 at 68% confidence level). However, there are many more unsuccessful models which leave the discrepancy well above the 3σ disagreement level. In many cases, reduced tension comes not simply from a change in the value of H 0 but also due to an increase in its uncertainty due to degeneracy with additional physics, complicating the picture and pointing to the need for additional probes. While no specific proposal makes a strong case for being highly likely or far better than all others, solutions involving early or dynamical dark energy, neutrino interactions, interacting cosmologies, primordial magnetic fields, and modified gravity provide the best options until a better alternative comes along.
High energy astrophysics: The contents of the Universe – the grand design
M. Longair
A study on sferic removal methods from actual airborne transient electromagnetic data
Ping Qi, Jiangyuan Chen, Zhihong Wang
et al.
In airborne transient electromagnetic (ATEM) surveys, sferic significantly degrades data quality due to its randomness and uncertain energy amplitude. This study explores three methods for removing sferic: wavelet decomposition and reconstruction, α-trimmed mean filter, and Hampel filter. Their respective advantages and application scenarios in signal processing are compared. Wavelet decomposition and reconstruction offer multiscale and multiresolution signal analysis, allowing simultaneous consideration of time and frequency characteristics. The choice of basis functions affects global properties and phase shifts, while the number of decomposition levels influences data smoothing, and threshold settings impact noise removal effectiveness. α-trimmed mean filtering removes noise by averaging the remaining data after excluding the maximum and minimum values within a window, thus demonstrating adaptability. The Hampel filter detects outliers using median and median absolute deviation (MAD), effectively targeting noise removal, particularly suitable for scenarios requiring signal detail preservation. This study enhances the Hampel filter, reducing computational overhead while ensuring accurate outlier detection. Simulation experiments on random signals and field-collected ATEM data validate the denoising capabilities of the three methods, revealing that the Hampel filter achieves optimal sferic reduction.
Geophysics. Cosmic physics, Astrophysics
Optimizing the potential of KM3NeT in detecting core-collapse supernovae
KM3NeT Collaboration O. Adriani, A. Albert, A. Alhebsi
et al.
Core-collapse supernovae mark the end of life of massive stars. However, despite their importance in astrophysics, their underlying mechanisms remain unclear. Neutrinos that emerge from the dense core of the star offer a promising way to study supernova dynamics. A strategy is presented to improve the potential of the KM3NeT neutrino telescope to detect core-collapse supernovae in our Galaxy or the Large Magellanic Cloud by further exploiting the properties of its optical modules equipped with multiple photomultipliers. A supernova burst is expected to produce a sudden hit rate increase in the KM3NeT detectors, which could be used to detect a supernova even in the absence of triggers from other experiments. New observables have been defined for individual optical modules that exploit the geometry and time distribution of the detected hits, enabling a better discrimination between signal and background signatures. In addition, a thorough investigation of the related systematic uncertainties is presented for the first time. When implemented, this new methodology allowed KM3NeT to probe 46% more Galactic core-collapse supernova candidates than with the previous trigger strategy, reaching the dense Galactic bulge. It is now expected that, once completed, KM3NeT will achieve full Galactic sensitivity to core-collapse supernovae independently from other experiments.
Thermodynamic Topology, Photon Spheres, and Evidence for Weak Gravity Conjecture in Charged Black Holes with Perfect Fluid within Rastall Theory
S. N. Gashti, .Izzet Sakalli, B. Pourhassan
In this paper, we explore the Weak Gravity Conjecture (WGC) within the context of photon spheres in charged black holes, framed by Perfect Fluid in Rastall Theory. We aim to validate the WGC by identifying the extremality states of these black holes. We highlight the interplay between quantum dynamics and gravitational forces, opening new avenues in high-energy physics and quantum gravity. Our analysis reveals significant system changes with varying perfect fluid intensity $(\alpha)$ and Rastall parameters ( $k$ and $\lambda$ ). For dust field $(\omega=0)$, the WGC is met in extremality $(T=0)$ with $(Q / M)_{\text {ext }}>1$, indicating a black hole due to the presence of photon spheres (PS) with a total charge-1. However, further increases in $k$ and $\lambda$ or decreases $(\alpha)$ lead to $P S=0$ that determines a singularity, not a black hole. We observed that the radiation field $(\omega=1 / 3)$, quintessence field ( $\omega=-2 / 3$ ), and phantom fields field ( $\omega=-4 / 3$ ) also confirmed the WGC and maintaining a total charge of $P S=-1$ in some regions of the free parameters. Our numerical solutions identify points satisfying the WGC, establishing a bridge between quantum and cosmic realms. The results are summarized in Table (I). We also examine Duans topological current $\phi$-mapping theory by analyzing generalized Helmholtz free energy methods to study the topological classes of our black hole. We reveal that for given values of the free parameters, the total topological numbers $(W=0)$ exist for the generalized Helmholtz free energy method for $\omega=0,1 / 3,-2 / 3,-4 / 3$.
A robust numerical method based on deep learning operator for solving the 2D acoustic wave equation
Yunfan Chang, Dinghui Yang, Xijun He
The numerical solution of wave equations plays a crucial role in computational geophysics problems, which forms the foundation of inverse problems and directly impact the high-precision imaging results of earth models. However, common numerical methods often lead to signifcant computational and storage requirements. Due to the heavy reliance on forward modeling methods in inversion techniques, particularly full waveform inversion, enhancing the computational efficiency and reducing storage demands of traditional numerical methods becomes a key issue in computational geophysics. In this paper, we present the deep Lax-Wendroff correction method (DeepLWC), a deep learning-based numerical format for solving two#xD; dimensional (2D) hyperbolic wave equations. DeepLWC combines the advantages of the traditional numerical schemes with a deep neural network. We provide a detailed comparison of this method with representative traditional Lax-Wendroff correction (LWC) method. Our numerical results indicate that the DeepLWC signifcantly improves calculation speed (by more than ten times) and reduces storage space by over 10000 times compared to traditional numerical methods. In contrast to the more popular Physics Informed Neural Network (PINN) method, DeepLWC maximizes the advantages of traditional mathematical methods in solving PDEs and employs a new sampling approach, leading to improved accuracy and faster computations. It is particularly worth pointing that, DeepLWC introduces a novel research paradigm for numerical equation-solving, which can be combined with various traditional numerical methods, enabling acceleration and reduction in storage requirements of conventional approaches.
The influence of regularization on the geomagnetic core field modeling
Juyuan Xu, Yufeng Lin, Keke Zhang
Earth's magnetic field is generated through the dynamo action in the liquid outer core. Using observational data from geomagnetic satellites and observatories, we can construct geomagnetic models that describe spatial and temporal variations of the core field. Due to the limited temporal and spatial resolution of observations and incomplete parameterization of the model, the inversion problem of geomagnetic field modeling is non-unique. Therefore, we need to incorporate regularization constraints based on prior information to alleviate the non-uniqueness of the inversion problem for geomagnetic modeling. For the core field, the non-uniqueness of modeling inversion manifests in the temporal variations of the field. Studies on the influence of regularization can help us build more reliable core field models, especially on how to suppress artificial secular variations through regularizations. In this study, we construct a series of models based on different regularization parameters within the framework of CHAOS modeling to explore the impact of regularization on the core field modeling. Our results indicate that it is necessary to use regularization constraints for the core field modeling based solely on satellite data. By comparing models with different regularization intensities of the third time derivative of the core field, we show that regularization of the third time derivative with appropriate intensity can effectively suppress artificial signals caused by overfitting of the model to observational data. However, too strong regularizations will reduce the temporal resolution of small-scale secular variations.
Geophysics. Cosmic physics, Astrophysics
Research on abutment stress distribution of roof-cutting coalface: numerical simulation and field measurement
Yangyang Guo, Xiaoli Liu, Weitao Li
et al.
Abstract During the processing of deep mining, revealing the distribution of abutment pressure is significant for controlling stability of the entry. In this study, the abutment pressure distribution of roof-cutting coalface was investigated by FLAC3D and self-developed flexible detection unit (FDU). In the numerical simulation, the double-yield model was built to analyze the goaf abutment pressure under the fracturing roofs to maintain entry (FRME). Compared with the non-fracturing side, the peak value of the advanced abutment pressure on the fracturing side is reduced by 19.29% on average, the influence range (span) increases by 30.78% and the distance between the peak value and the working face increases by 66.7%. The goaf abutment pressure within 23m near the cutting side is significantly higher than other areas along the dip. The FDU was employed in the coalface to record the change of advanced abutment stress. And the field measured results are in well agreement with the numerical results.
Geophysics. Cosmic physics
Sulfide Copper‐Iron Isotopic Fractionation During Formation of the Kalatongke Magmatic Cu‐Ni Sulfide Deposit in the Central Asian Orogenic Belt
Dongmei Tang, Kezhang Qin, Noreen J. Evans
et al.
Abstract Copper and iron isotopic signatures in sulfide and silicate minerals are important genetic indicators in magmatic sulfide deposits. Kalatongke is a large‐scale magmatic Cu‐Ni sulfide deposit in the Central Asian Orogenic Belt, and one that experienced multiple stages of magmatism and contamination. It is an ideal deposit in which to study Cu‐Fe isotopic fractionation during multiple stages of magmatism and sulfide mineralization processes. The Kalatongke sulfide orebodies are hosted by three small mafic intrusions in which pyroxene and sulfides (pyrrhotite, pentlandite, and chalcopyrite) are the most common Fe‐rich minerals, and chalcopyrite is the dominant Cu‐rich mineral. Sulfide liquid and silicate melt ▵56FeSul‐Sil (0.03–0.19‰) and ▵65CuCcp‐Sil (−0.78–0.74‰) values are indicative of non‐equilibrium fractionation. Most of the Cu isotope compositions in the sulfide ores at Kalatongke can be modeled as subduction‐ metasomatized, oxidized mantle source‐derived silicate melt (initial δ57Fe = 0.15‰, δ65Cu = −0.07‰) that underwent lower crustal contamination, and then reacted with silicate melt, having an R factor of 100–1,000. Rapid silicate melt and sulfide liquid Fe isotope exchange and re‐equilibration between chalcopyrite and pyrrhotite in the massive ores is reflected in the similarity of their δ56Fe values. Sulfide in disseminated ores shows a range of Fe isotope ratios, influenced by the proportions of monosulfide solid solution (MSS) and intermediate solid solution (ISS) formed. Copper isotopes can be utilized to characterize crustal contamination and silicate melt‐sulfide liquid interaction, while the Fe isotope ratios of sulfide minerals record sulfide liquid segregation and evolution in magmatic sulfide deposits.
Geophysics. Cosmic physics, Geology
Particulate Nitrate Photolysis as a Possible Driver of Rising Tropospheric Ozone
Viral Shah, Christoph A. Keller, K. Emma Knowland
et al.
Abstract Tropospheric ozone is an air pollutant and a greenhouse gas whose anthropogenic production is limited principally by the supply of nitrogen oxides (NOx) from combustion. Tropospheric ozone in the northern hemisphere has been rising despite the flattening of NOx emissions in recent decades. Here we propose that this sustained increase could result from the photolysis of nitrate particles (pNO3−) to regenerate NOx. Including pNO3− photolysis in the GEOS‐Chem atmospheric chemistry model improves the consistency with ozone observations. Our simulations show that pNO3− concentrations have increased since the 1960s because of rising ammonia and falling SO2 emissions, augmenting the increase in ozone in the northern extratropics by about 50% to better match the observed ozone trend. pNO3− will likely continue to increase through 2050, which would drive a continued increase in ozone even as NOx emissions decrease. More work is needed to better understand the mechanism and rates of pNO3− photolysis.
Geophysics. Cosmic physics
A sustainable clean energy source for mitigating CO2 emissions: numerical simulation of Hamit granitoid, Central Anatolian Massif
Tolga Ayzit, Mrityunjay Singh, Dornadula Chandrasekharam
et al.
Abstract Türkiye relies on coal-fired power plants for approximately 18 GW of annual electricity generation, with significantly higher CO2 emissions compared to geothermal power plants. On the other hand, geothermal energy resources, such as Enhanced Geothermal Systems (EGS) and hydrothermal systems, offer low CO2 emissions and baseload power, making them attractive clean energy sources. Radiogenic granitoid, with high heat generation capacity, is a potential and cleaner energy source using EGS. The Anatolian plateau hosts numerous tectonic zones with plutonic rocks containing high concentrations of radioactive elements, such as the Central Anatolian Massif. This study evaluates the power generation capacity of the Hamit granitoid (HG) and presents a thermo-hydraulic-mechanical (THM) model for a closed-loop geothermal well for harnessing heat from this granitoid. A sensitivity analysis based on fluid injection rates and wellbore length emphasizes the importance of fluid resident time for effective heat extraction. Closed-loop systems pose fewer geomechanical risks than fractured systems and can be developed through site selection, system design, and monitoring. Geothermal wellbore casing material must withstand high temperatures, corrosive environments, and should have low thermal conductivity. The HG exhibits the highest heat generation capacity among Anatolian granitoid intrusions and offers potential for sustainable energy development through EGS, thereby reducing CO2 emissions.
Geophysics. Cosmic physics
A resolution of the monopole problem in the <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" display="inline" id="d1e24" altimg="si4.svg"><mml:mrow><mml:msub><mml:mrow><mml:mi>R</mml:mi></mml:mrow><mml:mrow><mml:mi mathvariant="normal">h</mml:mi></mml:mrow></mml:msub><mml:mo linebreak="goodb
F. Melia
Spontaneous symmetry breaking in grand unified theories is thought to have produced an exceedingly large number of magnetic monopoles in the early Universe. In the absence of suppression or annihilation, these very massive particles should be dominating the cosmic energy budget today, but none has ever been found. Inflation was invented in part to dilute their number, thereby rendering their density undetectable by current instruments. Should the inflationary paradigm not survive, however, the ensuing disagreement between theory and observation would constitute a cosmological `monopole problem' and create further tension for any extension to the standard model of particle physics. But as is also true for all horizon problems, a monopole overabundance emerges only in cosmologies with an initial period of deceleration. We show that the alternative Friedmann-Lemaitre-Robertson-Walker cosmology known as the R_h=ct universe completely eliminates all such anomalies rather trivially and naturally, without the need for an inflated expansion. We find that the monopole energy density today would be completely undetectable in R_h=ct. Evidence continues to grow that the zero active mass condition from general relativity ought to be an essential ingredient in LCDM.
Plasma heating and particle acceleration in collisionless shocks through astrophysical observations
M. Miceli
Supernova remnants (SNRs), the products of stellar explosions, are powerful astrophysical laboratories, which allow us to study the physics of collisionless shocks, thanks to their bright electromagnetic emission. Blast wave shocks generated by supernovae (SNe) provide us with an observational window to study extreme conditions, characterized by high Mach (and Alfvénic Mach) numbers, together with powerful nonthermal processes. In collisionless shocks, temperature equilibration between different species may not be reached at the shock front. In this framework, different particle species may be heated at different temperatures (depending on their mass) in the post-shock medium of SNRs. SNRs are also characterized by broadband nonthermal emission stemming from the shock front as a result of nonthermal populations of leptons and hadrons. These particles, known as cosmic rays, are accelerated up to ultrarelativistic energies via diffusive shock acceleration. If SNRs lose a significant fraction of their ram energy to accelerate cosmic rays, the shock dynamics should be altered with respect to the adiabatic case. This shock modification should result in an increase in the total shock compression ratio with respect to the Rankine–Hugoniot value of 4. Here, I show that the combination of x-ray high resolution spectroscopy (to measure ion temperatures) and moderate resolution spectroscopy (for a detailed diagnostic of the post-shock density) can be exploited to study both the heating mechanism and the particle acceleration in collisionless shocks. I report on new results on the temperatures measured for different ion species in the remnant of the SN observed in 1987 in the Large Magellanic Cloud (SN 1987A). I also discuss evidence of shock modification recently obtained in the remnant of SN 1006 a. D., where the shock compression ratio increases significantly as the angle between the shock velocity and the ambient magnetic field is reduced.
Constraining primordial black hole masses through f(R) gravity scalarons in Big Bang Nucleosynthesis
Abhijit Talukdar, Sanjeev Kalita, Nirmali Das
et al.
Big Bang Nucleosynthesis (BBN) is a strong probe for constraining new physics including gravitation. f(R) gravity theory is an interesting alternative to general relativity which introduces additional degrees of freedom known as scalarons. In this work we demonstrate the existence of black hole solutions in f(R) gravity and develop a relation between scalaron mass and black hole mass. We have used observed bound on the freezeout temperature to constrain scalaron mass range by modifying the cosmic expansion rate at the BBN epoch. The mass range of primordial black holes (PBHs) which are astrophysical dark matter candidates is deduced. The range of scalaron mass which does not spoil the BBN era is found to be 10-16–104 eV for both relativistic and non-relativistic scalarons. The window 10-16–10-14 eV of scalaron mass obtained from solar system constraint on PPN parameter is compatible with the BBN bound derived in this work. The PBH mass range is obtained as 106–10-14 M ⊙. Scalarons constrained by BBN are also eligible to accommodate axion like dark matter particles. The problem of ultra-light PBHs (M ≤ 10-24 M ⊙) not constrained by the present study of BBN is still open. Estimation of deuterium (D) fraction and relative D+3He abundance in the f(R) gravity scenario shows that the BBN history mimics that of general relativity. While the PBH mass range is eligible for non-baryonic dark matter, the BBN bounded scalarons provide with an independent strong field test of f(R) gravity. The PBH mass range obtained in the study is discussed in relation to future astronomical measurements.
Nonlinear interaction of two non-uniform vortex sheets and large vorticity amplification in Richtmyer–Meshkov instability
C. Matsuoka, K. Nishihara
Vortex dynamics is an important research subject for geophysics, astrophysics, engineering, and plasma physics. Regarding vortex interactions, only limited problems, such as point vortex interactions, have been studied. Here, the nonlinear interaction of two non-uniform vortex sheets with density stratification is investigated using the vortex sheet model. These non-uniform vortex sheets appear, for example, in the Richtmyer–Meshkov instability that occurs when a shock wave crosses a layer with two corrugated interfaces. When a strong vortex sheet approaches a weaker vortex sheet with opposite-signed vorticity, a locally peaked secondary vorticity is induced on the latter sheet. This emerging secondary vorticity results in a remarkable vorticity amplification on the stronger sheet, and a strong vortex core is formed involving the weak vortex sheet. The amplified vortices with opposite signs on the two vortex sheets form pseudo-vortex pairs, which cause an intense rolling-up of the two sheets. We also investigated the dependence of distance and initial phase difference of vorticity perturbations between two vortex sheets on the vorticity amplification and vortex sheet dynamics.
Lyman-α polarization from cosmological ionization fronts. Part II. Implications for intensity mapping
Emily Koivu, Heyang Long, Yuanyuan Yang
et al.
This is the second paper in a series whose aim is to predict the power spectrum of intensity and polarized intensity from cosmic reionization fronts. After building the analytic models for intensity and polarized intensity calculations in paper I, here we apply these models to simulations of reionization. We construct a geometric model for identifying front boundaries, calculate the intensity and polarized intensity for each front, and compute a power spectrum of these results. This method was applied to different simulation sizes and resolutions, so we ensure that our results are convergent. We find that the power spectrum of fluctuations at z = 8 in a bin of width Δz = 0.5 (λ/Δλ = 18) is Δℓ ≡ [ℓ(ℓ + 1)C ℓ/2π]1/2 is 3.2 × 10-11 erg s-1 cm-2 sr-1 for the intensity I, 7.6 × 10-13 erg s-1 cm-2 sr-1 for the E-mode polarization, and 5.8 × 10-13 erg s-1 cm-2 sr-1 for the B-mode polarization at ℓ = 1.5 × 104. After computing the power spectrum, we compare results to detectable scales and discuss implications for observing this signal based on a proposed experiment. We find that, while fundamental physics does not exclude this kind of mapping from being attainable, an experiment would need to be highly ambitious and require significant advances to make mapping Lyman-α polarization from cosmic reionization fronts a feasible goal.
Intelligent Retrieval of Radar Reflectivity Factor With Privacy Protection Under Meteorological Satellite Remote Sensing
Huichao Lin, Xiaolong Xu, Muhammad Bilal
et al.
Meteorological radar data are essential for meteorological monitoring, forecasting, and research, and it plays a crucial role in observing and warning of extreme weather risks. However, meteorological radars have some limitations, such as uneven distribution and severe topographical influence. Meteorological remote sensing satellites can partially overcome these limitations by providing larger observational scope and high spatial and temporal resolution. Using data from meteorological remote sensing satellites to train radar reflectivity factor retrieval models can effectively compensate for the missing and poor quality of radar data. However, there are still some challenges, such as extracting the features of intense convective weather with unclear coverage from complex multichannel meteorological remote sensing satellite data and removing the interference caused by nonprecipitation clouds on retrieval models. Moreover, the privacy and security of remote sensing data transmission need to be ensured. In this article, we propose a novel method that combines the advanced encryption standard method to protect the transmission of remote sensing data, a multiscale feature fusion module to extract multiscale features from multichannel meteorological remote sensing satellite data, and an attention technique to reduce the interference of nonprecipitation clouds on retrieval models. We conduct comparison experiments with multiple indicators to demonstrate that our method has certain advantages in retrieving radar reflectivity values of different sizes. Our method achieves 0.63, 0.36, 0.49, 0.55, and 0.99 on probability of detection, false alarm ratio, critical success index, Heidke skill score, and accuracy scores, respectively.
Ocean engineering, Geophysics. Cosmic physics
Probing the scale of grand unification with gravitational waves
W. Buchmuller, V. Domcke, H. Murayama
et al.
Abstract The spontaneous breaking of U ( 1 ) B − L around the scale of grand unification can simultaneously account for hybrid inflation, leptogenesis, and neutralino dark matter, thus resolving three major puzzles of particle physics and cosmology in a single predictive framework. The B − L phase transition also results in a network of cosmic strings. If strong and electroweak interactions are unified in an S O ( 10 ) gauge group, containing U ( 1 ) B − L as a subgroup, these strings are metastable. In this case, they produce a stochastic background of gravitational waves that evades current pulsar timing bounds, but features a flat spectrum with amplitude h 2 Ω GW ∼ 10 − 8 at interferometer frequencies. Ongoing and future LIGO observations will hence probe the scale of B − L breaking.
Millicharged particles from the heavens: single- and multiple-scattering signatures
Carlos A. Argüelles, K. Kelly, V. Muñoz
For nearly a century, studying cosmic-ray air showers has driven progress in our understanding of elementary particle physics. In this work, we revisit the production of millicharged particles in these atmospheric showers and provide new constraints for XENON1T and Super-Kamiokande and new sensitivity estimates of current and future detectors, such as JUNO. We discuss distinct search strategies, specifically studies of single-energy-deposition events, where one electron in the detector receives a relatively large energy transfer, as well as multiple-scattering events consisting of (at least) two relatively small energy depositions. We demonstrate that these atmospheric search strategies — especially the multiple-scattering signature — provide significant room for improvement beyond existing searches, in a way that is complementary to anthropogenic, beam-based searches for MeV-GeV millicharged particles. Finally, we also discuss the implementation of a Monte Carlo simulation for millicharged particle detection in large-volume neutrino detectors, such as IceCube.
Estimating Inner Magnetospheric Radial Diffusion Using a Hybrid-Vlasov Simulation
H. George, A. Osmane, E. K. J. Kilpua
et al.
Radial diffusion coefficients quantify non-adiabatic transport of energetic particles by electromagnetic field fluctuations in planetary radiation belts. Theoretically, radial diffusion occurs for an ensemble of particles that experience irreversible violation of their third adiabatic invariant, which is equivalent to a change in their Roederer L* parameter. Thus, the Roederer L* coordinate is the fundamental quantity from which radial diffusion coefficients can be computed. In this study, we present a methodology to calculate the Lagrangian derivative of L* from global magnetospheric simulations, and test it with an application to Vlasiator, a hybrid-Vlasov model of near-Earth space. We use a Hamiltonian formalism for particles confined to closed drift shells with conserved first and second adiabatic invariants to compute changes in the guiding center drift paths due to electric and magnetic field fluctuations. We investigate the feasibility of this methodology by computing the time derivative of L* for an equatorial ultrarelativistic electron population travelling along four guiding center drift paths in the outer radiation belt during a 5 minute portion of a Vlasiator simulation. Radial diffusion in this simulation is primarily driven by ultralow frequency waves in the Pc3 range (10–45 s period range) that are generated in the foreshock and transmitted through the magnetopause to the outer radiation belt environment. Our results show that an alternative methodology to compute detailed radial diffusion transport is now available and could form the basis for comparison studies between numerical and observational measurements of radial transport in the Earth’s radiation belts.
Astronomy, Geophysics. Cosmic physics