Search NASA⌕ Search

SEARCH · Search NASA

Results for “distribution shift”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

I/O in Machine Learning Applications on HPC Systems: A 360-degree Survey

Growing interest in Artificial Intelligence (AI) has resulted in a surge in demand for faster methods of Machine Learning (ML) model training and inference. This demand for speed has prompted the use of high performance computing (HPC) systems that excel in managing distributed workloads. Because data is the main fuel for AI applications, the performance of the storage and I/O subsystem of HPC systems is critical. In the past, HPC applications accessed large portions of data written by simulations or experiments or ingested data for visualizations or analysis tasks. ML workloads perform small reads spread across a large number of random files. This shift of I/O access patterns poses several challenges to modern parallel storage systems. In this paper, we survey I/O in ML applications on HPC systems, and target literature within a 6-year time window from 2019 to 2024. We define the scope of the survey, provide an overview of the common phases of ML, review available profilers and benchmarks, examine the I/O patterns encountered during offline data preparation, training, and inference, and explore I/O optimizations utilized in modern ML frameworks and proposed in recent literature. Lastly, we seek to expose research gaps that could spawn further R&D.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Mechanistic Tuning of Chemical Transformations for Coupling the Geo-mimicry of Acid Gas Storage with Design Strategies to Produce Clean Energy Carriers in Multi-Phase Reaction Environments (MATTER) (Final Report)

The need to diversify approaches to produce essential energy carriers such as H 2 motivate advances in thermodynamically downhill geo-inspired pathways. Currently, more than 80% of hydrogen is produced via steam methane reforming (SMR) followed by water gas shift reaction (WGSR) pathway with the co-production of CO 2 . As an alternative to introducing CO 2 separation strategies downstream such as the use of membranes, solvents, or sorbents, geo-inspired carbon mineralization is harnessed as an alternative crystallization pathway. The thermodynamically downhill carbon mineralization pathways are hypothesized to accelerate H2 conversion while separating CO 2 . The work conducted through this project discusses the mechanisms associated with coupling the water gas shift reaction (WGSR) with carbon mineralization. In addition to harnessing Ca- and Mg-bearing oxides or hydroxides that are known to be reactive for CO 2 capture, the use of earth abundant silicate minerals such as Ca- and Mg-bearing silicates, is also investigated. Key outcomes include approaches to architect Mg-silicate with well-controlled pore size distributions, probing the enhancement in H 2 yield with inherent CO 2 suppression using Ca- and Mg-bearing oxides, hydroxides, and silicates, elucidating the changes in silicate chemistry during carbon mineralization, and exploring direct integration of carbon mineralization with WGSR and the development of separate low temperature pathway for reactive CO 2 capture and mineralization are developed.

08 HYDROGEN↗
DOE OSTI2025

Scalability and Effectiveness of Smart Charge Management

The rise in electric vehicle (EV) adoption presents growing challenges for power grids, particularly from simultaneous residential charging, which can cause voltage fluctuations and increase feeder peak loads. Baltimore Gas and Electric (BGE), with support from the U.S. Department of Energy, initiated a pilot program to evaluate managed residential EV charging through Smart Charge Management (SCM). This study analyzes real-world charging behavior data from the pilot and feeder-level base loads from BGE to simulate residential charging scenarios through 2035 across the Washington, DC–Baltimore region. Grid impacts under unmanaged charging are compared to three SCM strategies: TOU-immediate, TOU-distributed, and Load Balancing. Results show that the magnitude of peak reduction is highly feeder-dependent. Some feeders achieve reductions of more than 40% at high enrollment levels, while others show improvements closer to 10–15%. This heterogeneity reflects differences in baseline feeder load shapes, EV penetration, and plug-in behavior across customers. Results also highlight trade-offs between shifting load away from peak periods and minimizing secondary demand peaks, offering practical insights for future utility program design.

Electric vehicle↗
DOE OSTI2026

SPRUCE FT-ICR MS, Bulk Chemistry, and Mass Loss from Litter Decomposition Study in Experimental Plots, Marcell Experimental Forest, Minnesota, 2015-2017

This dataset contains molecular, bulk chemical, and mass loss measurements from a litter decomposition study at the Spruce and Peatland Responses Under Changing Environments (SPRUCE) experimental site within the Marcell Experimental Forest in northern Minnesota, USA. This site is in a Sphagnum spp. ombrotrophic bog forest. Litterbags were deployed into the peat in September 2015 across three warming levels (+0, +4.5, and +9°C) under ambient and elevated carbon dioxide (CO₂ - +500 ppm) and retrieved after roughly 0.5, 1, and 2 years of field incubation (2015-09-23 to 2017-08-02). Litterbags containing six peatland litter types: black spruce needles (Picea mariana - SPL), spruce fine roots (SPR), Sphagnum angustifolium (ANG), Sphagnum magellanicum (MAG), Labrador tea leaves (Rhododendron groenlandicum - LTL), and Labrador tea roots (LTR). Molecular composition of water-soluble organic matter extracts was characterized using Fourier Transform Ion Cyclotron Resonance Mass Spectrometry (FT-ICR MS) at 9.4 Tesla, operated in negative ion mode with electrospray ionization, providing molecular formula assignments and compound-class distributions across the decomposition time series. Bulk chemical characterization included elemental analysis (percent carbon, nitrogen, and phosphorus) and Fourier Transform Infrared Spectroscopy (FTIR) to quantify functional group composition. Litter mass loss was tracked gravimetrically at each retrieval interval, expressed as percent mass remaining relative to initial dry mass for each litter type and treatment combination. These data are valuable for understanding how vegetation shifts driven by increased atmospheric CO2 and temperature in peatlands alter litter inputs and organic matter stabilization trajectories, with implications for projecting and modeling peatland carbon cycling. This dataset contains two data files in comma-separated value (.csv) format. Additional metadata are provided: two data dictionaries and a file-level metadata file in comma separate (.csv) format and a user guide in PDF (*.pdf) format.

decomposition↗
DOE OSTI2025

Mapping the structural–mechanical landscape of amorphous carbon with ReaxFF molecular dynamics

We use ReaxFF molecular dynamics (MD) to investigate the relationship between structural and mechanical properties in bulk and nanostructured amorphous carbon (a-C). The liquid-quench MD method is used to generate isotropic bulk samples with mass densities ranging from 0.96 to 3.29 g/cm3. Structural analysis identifies two types of structures with distinct short- and medium-range order: lower-density sp2-dominated a-C, which is characterized by a bimodal ring-size distribution, and higher-density sp3-dominated tetrahedral amorphous carbon (ta-C), exhibiting a unimodal ring-size distribution. Stress–strain MD simulations and analysis reveal how an atomistic structure impacts elastic properties and post-yield atomic rearrangements. All stretched structures demonstrate elastic isotropy and plasticity driven by a ring-size expansion mechanism reflected in changes in ring statistics. The plastic region is substantially larger in ta-C than in a-C due to the post-yield shift from sp3 to sp2 C dominant bonding. In both a-C and ta-C, ultimate failure occurs when a reactive crack, traversed by long sp chains, forms and propagates predominantly perpendicular to the direction of the applied strain. Oxygen infiltration into the fractured region significantly reduces stress resistance, primarily through the early rupture of long sp chains. MD simulations and analysis are extended to a-C slabs, a-C nanotubes, and partially a-C nanotubes. The latter nanostructure highlights the differences between the elastically isotropic a-C walls, which develop circumferential cracking, and the crystalline walls, which tear along crystallographic directions. These results provide a strong foundation for further computational characterization of a-C materials.

Dernov, A. (ORCID:0009000004220973)↗
DOE OSTI2024

Studies on the compositional dependent structural and electrical properties of CaTiO3-modified K0.5Na0.5NbO3 piezoelectric system

Lead-free piezoelectric ceramics of (1 − x)K0.5Na0.5NbO3-xCaTiO3 were fabricated, and their crystal structure, microstructure, and electrical properties were systematically studied. Rietveld refinement of the x-ray diffraction data and Raman spectroscopic analyses revealed a composition-dependent structural phase transition: three phase transitions, namely, from a pure orthorhombic phase for x ≤ 0.02 to a mixed phase of orthorhombic and tetragonal phases (0.03 ≤ x ≤ 0.08) and finally another mixed phase of tetragonal + cubic for x = 0.10 and 0.15 at room temperature (RT). The morphological study reveals a decrease in grain size along with a more uniform distribution of grains as the concentration of CaTiO3 (CT) increases; notably, a homogeneous distribution of grains is observed for x = 0.05. The temperature-dependent dielectric properties show two phase transitions, from orthorhombic to tetragonal (TO-T) and tetragonal to cubic (TC), for unmodified K0.5Na0.5NbO3 (KNN). However, both the phase transition temperatures (TO-T and TC) decrease, and the transition peaks broaden with an increase in CT substitution, and for x > 0.06, the TO-T shifted below RT. The broadening of the transition peak at TO-T may be due to the relaxation behavior. Among the prepared samples, the 5 mol. % CT-modified KNN shows the optimum electrical properties (d33 = 114 pC/N, ɛr = 412, and 2Pr = 15.25 μC/cm2) at RT. The enhanced electrical properties for x = 0.05 are due to the coexistence of orthorhombic and tetragonal phases, facilitating easy polarization rotation and flattening of the free energy profile. A phase diagram has been constructed based on the information gathered from the temperature-dependent dielectric measurements, RT x-ray diffraction, and Raman spectroscopy data and is discussed in detail.

Physics↗
DOE OSTI2025

Excitation of whistler and slow-X waves by runaway electrons in a collisional plasma

Runaway electrons are known to provide robust ideal or collisionless kinetic drive for plasma wave instabilities in both the whistler and slow-X branches, via the anomalous Doppler-shifted cyclotron resonances. In a cold and dense post-thermal-quench plasma, collisional damping of the plasma waves can compete with the collisionless drive. Previous studies have found that, due to their higher wavelength and frequency, slow-X waves suffer stronger collisional damping than the whistlers, while the ideal growth rate of slow-X modes is higher. Here, we study runaway avalanche distributions that maintain the same eigen distribution and increase only in magnitude over time. The distributions are computed from the relativistic Fokker–Planck–Boltzmann solver, upon which a linear dispersion analysis is performed to search for the most unstable or least damped slow-X and whistler modes. Taking into account the effect of plasma density, plasma temperature, and effective charge number, we find that the slow-X modes tend to be excited before the whistlers in a runaway current ramp-up. Furthermore, even when the runaway current density is sufficiently high that both branches are excited, the most unstable slow-X mode has a much higher growth rate than the most unstable whistler mode. The qualitative and quantitative trends uncovered in the current study indicate that even though past experiments and modeling efforts have concentrated on whistler modes, there is a compelling case that slow-X modes should also be a key area of focus in the runaway self-mediation through wave instabilities.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

The Dark Energy Bedrock All-sky Supernova Program: Cross Calibration, Simulations, and Cosmology Forecasts

Type Ia supernovae (SNe Ia) have been essential for probing the nature of dark energy; however, most SN analyses rely on the same low-redshift sample, which may lead to shared systematics. In a companion paper, we introduce the Dark Energy Bedrock All-Sky Supernova (DEBASS) program, which has already collected more than 500 low-redshift SNe Ia on the Dark Energy Camera, and present an initial release of 77 SNe Ia within the Dark Energy Survey (DES) footprint observed between 2021 and 2024. Here, we examine the systematics, including photometric calibration and selection effects. We find agreement at the 10 mmag level among the tertiary standard stars of DEBASS, DES, and Pan-STARRS1. Our simulations reproduce the observed distributions of DEBASS SN light-curve properties, and we measure a bias-corrected Hubble residual scatter of 0.08 mag, which, while small, is found in 10% of our simulations. We compare the DEBASS SN distances to the Foundation sample and find consistency with a median residual offset of 0.016 ± 0.019 mag. Selection effects have negligible impacts on distances, but a different photometric calibration solution shifts the median residual −0.015 ± 0.019 mag, highlighting calibration sensitivity. Using conservative simulations, we forecast that replacing historical low-redshift samples with the full DEBASS sample will improve the statistical uncertainties on dark energy parameters w 0 and w a by 30% and 24%, respectively, enhance the dark energy Figure of Merit by up to 60%, and enable a measurement of fσ 8 at the 25% level.

Acevedo, Maria [Duke Univ., Durham, NC (United Sta↗
DOE OSTI2024

High-intensity polarized electron gun featuring distributed Bragg reflector GaAs photocathode

The polarized electron source is a critical component in accelerator facilities such as the electron–ion collider, which requires a polarized electron gun with higher voltage and higher bunch charge than existing sources. One challenge we faced was the surface charge limit of the distributed Bragg reflector GaAs/GaAsP superlattice (DBR-SL-GaAs) photocathode. We suppressed this effect by optimizing the surface doping and heat cleaning procedures. We achieved up to 11.6 nC bunch charge of polarized electron beam. In this report, we discuss the performance of tests of a DBR-SL-GaAs photocathode in the high voltage direct current gun. Possible reasons for the observed peak quantum efficiency wavelength shift are analyzed, and we addressed it by using a wavelength tunable laser. In addition, the impact of the DBR layer and laser on the lifetime is investigated in this paper. The optimal DBR-SL-GaAs operating zone has been proposed, which gave us a long lifetime and high polarization at 30 μA operation. Finally, the success of this polarized gun will be key to the future of the nuclear sciences.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Proton Beam Dynamics in Bare IOTA with Intense Space-Charge

We are commissioning a 2.5~MeV proton beam for the Integrable Optics Test Accelerator at Fermilab, allowing experiments in the strong space-charge regime with incoherent betatron tune shifts nearing 0.5. Accurate modelling of space-charge dynamics is vital for understanding planned experiments. We compare anticipated emittance growth and beam loss in the bare IOTA configuration using transverse space-charge models in Xsuite, PyORBIT, and MAD-X simulation codes. Our findings reveal agreement within 30\% in core phase-space density predictions up to 100~synchrotron periods at moderate beam currents, while tail distributions and beam loss show significant differences.

43 PARTICLE ACCELERATORS↗
DOE OSTI2025

Quantifications of growth rate and particle size of magnesium silicate hydrate cement using small‐angle X‐ray scattering

The focus of recent research on cement binders has shifted to alternatives such as magnesium silicate hydrate (M-S-H) to mitigate the environmental impact of conventional concrete production. However, a detailed understanding of the morphology and growth rate of M-S-H, which informs its applicability as a cementing agent, is lacking. These properties are important as they influence cement hydration and consequently the development of macroscopic properties. Here, in this study, the particle size distribution and growth rate of M-S-H synthesized by mixing aqueous stock solutions were analyzed using synchrotron-based small-angle X-ray scattering. The reaction solutions contained dissolved magnesium and silicon with concentrations ranging from 2 to 500 mM, corresponding to saturation indices with respect to M-S-H endmembers, M 0.75 -S-H and M 1.50 -S-H, between 5.8 and 10.3. The precipitation rates of M-S-H align with affinity-based kinetic models and can be described by the following equations: For M 0.75 -S-H: 𝑟 M 0.75 −S−H =10 −9.54±0.70 ⁢(Ω$^{0.59±0.08}_{M_{0.75}⁢\textrm{SH}}$−1), and for M 1.50 -S-H: 𝑟 M 1.50 −S−H =10 −7.56±0.37 ⁢(Ω$^{0.34±0.04}_{M_{1.50}⁢\textrm{SH}}$−1), where 𝑟 is the rate in mol/L/s, and Ω represents the saturation ratio with respect to M-S-H. Moreover, we observed that changing the [Mg]/[Si] ratio affects particle size, a trend that cannot be fully explained by changes in the saturation index alone. These insights are critical for improving the predictability and control of cementation in concrete that utilizes M-S-H-based binders.

M-S-H↗
DOE OSTI2024

Proton beam dynamics in bare IOTA with intense space-charge

We are commissioning a 2.5-MeV proton beam for the Integrable Optics Test Accelerator at Fermilab, allowing experiments in the strong space-charge regime with incoherent betatron tune shifts nearing 0.5. Accurate modelling of space-charge dynamics is vital for understanding planned experiments. We compare anticipated emittance growth and beam loss in the bare IOTA configuration using transverse space-charge models in Xsuite, PyORBIT, and MADX simulation codes. Our findings reveal agreement within a factor of 2 in core phase-space density predictions up to 100 synchrotron periods at moderate beam currents, while tail distributions and beam loss show significant differences.

43 PARTICLE ACCELERATORS↗
DOE OSTI2026

Effect of glutathione-coated Mn-doped ZnS quantum dots on nutrient delivery in basil ( Ocimum basilicum ) plants

Ensuring efficient nutrient delivery while minimizing environmental impacts remains a significant challenge for modern agriculture. Nanotechnology-based fertilizers offer promising strategies to improve nutrient uptake and bioavailability in plants. This research aims to evaluate the use of Glutathione-coated Manganese-doped Zinc Sulfide quantum dots (GSH-ZnS-Mn QDs) as a potential nano fertilizer for basil (Ocimum basilicum). QDs' physicochemical properties were characterized using UV–Vis spectroscopy, photoluminescence, FTIR, and energy-dispersive X-ray spectroscopy, confirming successful Mn doping and glutathione surface functionalization. Basil plants were exposed to different concentrations of GSH-ZnS-Mn QDs under soil and hydroponic conditions. Plant growth parameters, oxidative stress responses, photosynthetic pigments, and macro- and micronutrient uptake were assessed using biochemical assays and inductively coupled plasma optical emission spectrometry (ICP-OES). Elemental uptake, spatial distribution, and zinc speciation were further investigated using synchrotron-based micro-X-ray fluorescence (μ-XRF) imaging and X-ray absorption near-edge structure (XANES) spectroscopy. Results show that exposure to GSH-ZnS-Mn QDs resulted in a concentration-dependent increase in leaf and stem biomass, accompanied by enhanced Zn accumulation in plant tissues. Catalase activity decreased across all tested concentrations, suggesting a shift toward glutathione-dependent antioxidant pathways rather than oxidative damage. Chlorophyll levels exhibited moderate reductions at higher concentrations. The higher increase in macronutrient (K, Ca, and Mg) uptake was reported in plants exposed to 200 ppm of QDs. μ-XRF imaging indicated a selective accumulation of Zn in roots and stems, with partial translocation to leaves. XANES analyses revealed that Zn from QDs was mainly converted into organic Zn species, such as Zn-phytate, Zn-acetate, and Zn-cysteine, indicating transformation and complexation within the plant. The findings demonstrate that a glutathione coating on GSH-ZnS-Mn QDs improves biocompatibility and nutrient delivery efficiency. These results highlight the relevance of surface functionalization in regulating nanoparticle fate, transformation, and nutrient bioavailability, supporting the potential application of GSH–ZnS–Mn QDs as modern nano fertilizers.

Basil↗
DOE OSTI2026

Group-Additivity–Embedded Multiscale Modeling for Electric Field-Enhanced Nanocatalysis

Elucidating structure-performance relationships remains a central challenge in field-enhanced catalysis, where nanoparticles exhibit nonuniform surface sites with site-dependent responses to electric fields. Low-coordination sites (edges, corners, and tips) are particularly electric field-sensitive (EF), leading to nonuniform charge distribution, adsorption energies, and catalytic activity. Here, using ammonia decomposition on a ruthenium cluster as a model system, we develop a transferable multiscale framework integrating density functional theory, group additivity (GA), Brønsted-Evans-Polanyi scaling, and microkinetic modeling to predict EF-dependent activity across nonuniform cluster sites. Across sites and fields, the nitrogen adsorption energy (E N ) emerges as the governing descriptor, yielding robust volcano relationships whose optimum shifts systematically with field: negative fields strengthen N binding via electron accumulation, while positive fields weaken N binding via charge depletion, moving the optimal E N toward weaker binding. Microkinetic analysis shows that N≡N bond formation remains the key kinetic bottleneck over most conditions; positive fields lower the effective barrier and, critically, increase the fraction of near-optimal active sites, leading to a net enhancement in overall activity relative to zero-field and negative-field cases. By capturing EF- and site-dependent energetics with high accuracy and low computational cost, this GA-embedded multi-scale simulation workflow provides a physically interpretable route to predict and design field-enhanced nanocatalysis.

ammonia decomposition↗
DOE OSTI2024

Finding a needle in a haystack: quantitative HERFD-XRF imaging and HERFD-XANES characterization of trace platinum in gold solidi from the Late Roman and Byzantine Empires

High-Energy Resolution Fluorescence Detection X-Ray Fluorescence (HERFD-XRF) imaging and HERFD X-ray Absorption Near Edge Structure (XANES) spectroscopy are used to quantify and characterize trace platinum (Pt) in gold solidi from the Late Roman and Byzantine Empires. Historically, the elemental analysis of coins has been pivotal in distinguishing authentic artifacts from forgeries, elucidating minting practices, and understanding economic shifts. Notably, a new gold source with high platinum content appeared in the fourth century CE, transforming the Roman economy. Traditional methods struggled to detect platinum due to the overwhelming gold matrix. Here, this study demonstrates the effectiveness of HERFD techniques in resolving this challenge. Three gold solidi, minted between 654 and 659 CE, were analyzed alongside reference gold materials with known Pt concentrations. The HERFD-XRF imaging revealed spatial distributions of platinum, highlighting non-uniformities within the coins. Additionally, HERFD-XANES spectroscopy identified the oxidation states and chemical speciation of platinum. Results demonstrate that platinum in the solidi primarily exists as metallic Pt, with some surface oxidation. The findings align with previous measurements but reveal higher Pt concentrations and significant inhomogeneities. This research confirms the reliability of HERFD methods for quantifying trace elements and provides new insights into the raw material sources and minting techniques of ancient gold coins. The non-destructive nature of this approach allows for extensive analyses, offering valuable data for historical, economic, and archaeological studies. This innovative application of HERFD-XRF imaging and XANES in cultural heritage research underscores the potential for detailed material characterization and conservation, enhancing our understanding of ancient economies and trade patterns.

Van Loon, Lisa L.↗
DOE OSTI2025

Enhancing DESI DR1 full-shape analyses using HOD-informed priors

We present an analysis of DESI Data Release 1 (DR1) that incorporates Halo Occupation Distribution (HOD)-informed priors into Full-Shape (FS) modeling of the power spectrum based on cosmological perturbation theory (PT). By leveraging physical insights from the galaxy-halo connection, these HOD-informed priors on nuisance parameters substantially mitigate projection effects in extended cosmological models that allow for dynamical dark energy. The resulting credible intervals now encompass the posterior maximum from the baseline analysis using gaussian priors, eliminating a significant posterior shift observed in baseline studies. In the ΛCDM framework, a combined DESI DR1 FS information and constraints from the DESI DR1 baryon acoustic oscillations (BAO) — including Big Bang Nucleosynthesis (BBN) constraints and a weak prior on the scalar spectral index — yields Ω m = 0.2994 ± 0.0090 and σ 8 = 0.836$^{+0.024}_{-0.027}$, representing improvements of approximately 4% and 23% over the baseline analysis, respectively. For the w 0 w a CDM model, our results from various data combinations are highly consistent, with all configurations converging to a region with w 0 > -1 and w a < 0. This convergence not only suggests intriguing hints of dynamical dark energy but also underscores the robustness of our HOD-informed prior approach in delivering reliable cosmological constraints.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Gating Single Molecules with Counterions

We report atomic-scale gating and visualization of local charge distribution within individual rare-earth-based molecular complexes on a metallic surface. The complexes are formed by a positively charged lanthanum ion coordinated to a (pcam)3 molecule and a negatively charged counterion trapped underneath via electrostatic interactions on a Au(111) surface. Local gating is performed by adding an additional negatively charged counterion to one side of the complex, which results in the redistribution of charges within the complex and a positive shift of the frontier orbitals. This is caused by the internal Stark effect induced by the added counterion. This effect is directly captured using tunneling spectroscopy and spectroscopic mapping at 5 K substrate temperature. The polarizability of the complex is corroborated by density functional theory and analytical calculations based on experimental findings. Furthermore, the influence of charge polarization on nearby complexes is investigated in a cluster purposely assembled using three complexes, which reveals maintaining the charge states as in single complexes. These findings will enable the design of robust charged rare-earth complexes to be tailored for potential solid-state applications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2026

Modeling the Effect of the Heliospheric Magnetic Field on Cosmic Ray Muon Shadows

Shadows cast in the cosmic ray (CR) muon sky by the Sun were located using muon data from the MINOS far detector in Northern Minnesota. The shadows were observed independently across three time periods; near solar minimum, near solar maximum, and over the entire 13 year span of the data. A distribution of muon positions for each shadow was then sampled to simulate CR motions near the Sun using the Parker spiral model of the Heliospheric Magnetic Field (HMF) and a dipole model of the Geomagnetic Field (GMF). The resulting particle distributions were then compared to their position with respect to the Sun. Results show that the Parker spiral model is most consistent with the solar minimum shadow and least consistent with the solar maximum shadow, as expected. The simple Parker spiral is more consistent with the data for a harder CR spectrum than is actually present, indicating the need for a more detailed HMF model. Plausible modifications to the Parker spiral model which would affect the overall shift of the Sun's CR shadow are discussed.

Adamson, P. [Fermilab]↗