Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “flight state”

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 55 records · Page 3

AstroPix: A pixelated HVCMOS sensor for space-based gamma-ray measurement

A next-generation medium-energy gamma-ray telescope targeting the MeV range would address open questions in astrophysics regarding how extreme conditions accelerate cosmic-ray particles, produce relativistic jet outflows, and more. One concept, AMEGO-X, relies upon the mission-enabling CMOS Monolithic Active Pixel Sensor silicon chip AstroPix. AstroPix is designed for space-based use, featuring low noise, low power consumption, and high scalability. Desired performance of the device include an energy resolution of 5 keV (or 10% FWHM) at 122 keV and a dynamic range per-pixel of 25–700 keV, enabled by the addition of a high-voltage bias to each pixel which supports a depletion depth of 500 μ m. This work reports on the status of the AstroPix development process with emphasis on the current version under test, version three (v3), and highlights of version two (v2). Version 3 achieves energy resolution of 10.4 ± 3.2% at 59.5 keV and 94 ± 6 μ m depletion in a low-resistivity test silicon substrate.

Astrophysics instrumentation↗

A U.S. scientific community review of carbon cycle science gaps and opportunities to better support earth system science and carbon management

Greenhouse gas (GHG) emissions continue to grow, while natural carbon reservoirs are becoming increasingly vulnerable to anthropogenic pressures, climate extremes, and disturbance. These changes are impacting humans, ecosystems, and natural resources worldwide. Tracking and mitigating GHG emissions require a pivot to operational monitoring of regional carbon flux and stock changes. The current GHG observing system is addressing needs at two distinct scales: 1) Local scale (< 1 km), related to anthropogenic point source emissions, and 2) global scales (> 1000 km), related to land and ocean carbon sinks. More focus on intermediate (10–1000 km) scales is needed to more effectively monitor progress in reducing carbon emissions, enhancing removals, and maintaining sinks. Representatives from carbon cycle biomass and flux communities across United States government agencies and academic institutions met in September 2024 to discuss the rationale and scientific context for more effectively implementing an operational system for GHG monitoring in support of urban and national carbon management needs. To guide development of this system, we propose a multi-tiered global spaceborne observing framework for carbon flux and stock, prioritizing: 1) frequent GHG partial columns for carbon emissions and removals; 2) continuous time series and data fusion of biomass from Lidar and Synthetic Aperture Radar (SAR) for carbon stocks, and 3) expanded coverage of tropical, high latitude, and oceanic regions to monitor carbon cycle tipping points and feedbacks. This system should be complemented by expanded surface and airborne networks for oceanic and terrestrial/aquatic ecosystems for calibration, ground truthing, and study of under-sampled regions.

54 ENVIRONMENTAL SCIENCES↗

Remote sensing of planetary boundary layer from ground and space: Structure, thermodynamics, and boundary layer clouds

The planetary boundary layer (PBL) is a critical interface between Earth's surface and atmosphere, influencing atmospheric convection, weather patterns, and air quality. Recognized by the 2017 National Academies of Sciences, Engineering, and Medicine Earth Science Decadal Survey as an Incubation Targeted Observable, high-quality and effective observation of the PBL has become a priority. This paper reviews surface-based and satellite remote sensing techniques for characterizing PBL features, including PBL height (PBLH), boundary layer thermodynamics, turbulence, and boundary layer clouds (BLCs). These elements are treated as interconnected aspects of the PBL system, while recognizing that different instruments retrieve different physical manifestations of the PBL. The review summarizes recent advances and limitations in Micro-Pulse Lidar (MPL), Doppler lidar, Raman lidar, Differential Absorption Lidar (DIAL), ceilometers, wind profilers, GNSS Radio Occultation, radar, and hyperspectral sounders. These observations help trace dynamic processes within the PBL and link PBL structure to broader weather and climate processes. The review also highlights persistent observational gaps over oceans, remote land, and polar regions, where continuous surface-based profiling is sparse. By comparing instrument capabilities, practical limitations, and interpretation issues, this review emphasizes the need for integrated remote sensing approaches and careful definition of the retrieved PBL quantities.

Boundary layer clouds↗

ELM–Wet: Inclusion of a Wet–Landunit With Sub–Grid Representation of Eco–Hydrological Patches and Hydrological Forcing Improves Methane Emission Estimations in the E3SM Land Model (ELM)

Wetlands are the largest emitters of biogenic methane (CH 4 ) and represent the highest source of uncertainty in global CH 4 budgets. Here, we aim to improve the realism of wetland representation in the U.S. Department of Energy's Exascale Earth System Model land surface model, ELM, thereby reducing uncertainty of CH 4 flux predictions. We develop an updated version, ELM-Wet, where we activate a separate landunit for wetlands that handles multiple wetland-specific eco-hydrological patch functional types. We introduce more realistic hydrological forcing through prescribing site-level constraints on surface water elevation, which allows resolving different sustained inundation depth for different patches, and if data exists, prescribing inundation depth. We modified the calculation of aerenchyma transport diffusivity based on observed conductance per leaf area for different vegetation types. We use Bayesian Optimization to parameterize CO 2 and CH 4 fluxes in the developed wet-landunit. Site-level simulations of a coastal non-tidal freshwater wetland in Louisiana were performed with the updated model. Eddy covariance observations of CO 2 and CH 4 fluxes from 2012 to 2013 were used to train the model and data from 2021 were used for validation. Patch-specific chamber flux observations and observations of CH 4 concentration profiles in the soil porewater from 2021 were used for evaluation of the model performance. Our results show that ELM-Wet reduces the model's CH 4 emission root mean squared error by up to 33% and is able to represent inter-daily CO 2 and CH 4 flux variability across the wetland's eco-hydrological patches, including during periods of extreme dry or wet conditions.

54 ENVIRONMENTAL SCIENCES↗

On Alfvénic turbulence of solar wind streams observed by Solar Orbiter during March 2022 perihelion and their source regions

It has been recently accepted that the standard classification of the solar wind solely according to flow speed is outdated, and particular interest has been devoted to the study of the origin and evolution of so-called Alfvénic slow solar wind streams and to what extent such streams resemble or differ from fast wind. In March 2022, Solar Orbiter completed its first nominal phase perihelion passage. During this interval, it observed several Alfvénic streams, allowing for characterization of fluctuations in three slow wind intervals (AS1-AS3) and comparison with a fast wind stream (F) at almost the same heliocentric distance. This work makes use of Solar Orbiter plasma parameters from the Solar Wind Analyzer (SWA) and magnetic field measurements from the magnetometer (MAG). The magnetic connectivity to the solar sources of selected solar wind intervals was reconstructed using a ballistic extrapolation based on measured solar wind speed down to the (spherical) source surface at 2.5 R s below which a potential field extrapolation was used to map back to the Sun. The source regions were identified using SDO/AIA observations. A spectral analysis of in situ measured magnetic field and velocity fluctuations was performed to characterize correlations, Alfvénicity, normalized cross-helicity, and residual energy in the frequency domain as well as intermittency of the fluctuations and spectral energy transfer rate estimated via mixed third-order moments. A machine learning technique was used to separate proton core, proton beam, and alpha particles and to study v − b correlations for the different ion populations in order to evaluate the role played by each population in determining the Alfvénic content of solar wind fluctuations. The comparison between fast wind and Alfvénic slow wind intervals highlights the differences between the two solar wind regimes: The fast wind is characterized by larger amplitude fluctuations, and magnetic and velocity fluctuations are closer to equipartition of energy. In fact the Alfvénic slow wind streams appear to be on a spectrum of wind types, with AS1, originating from open field lines neighboring active regions and displaying similarities with the fast wind in terms of fluctuation amplitude and turbulence characteristics, but not with respect to the alpha particles and proton beams. The other two slow streams differed both in their sources as well as plasma characteristics, with AS2 coming from the expansion of a narrow coronal hole corridor and AS3 from a region straddling a pseudostreamer. The latter displayed the coldest and highest density but the slowest stream with the smallest fluctuation amplitude and greatest magnetic energy excess. It also showed the largest scatter in proton beam speeds and the greatest difference in speed between proton beam and alpha particles. This study shows how the old fast–slow solar wind dichotomy, already called into question by the observations of slower Alfvénic solar wind streams, should further be refined, as the Alfvénic slow wind, originating in different solar wind regions, show significant differences in density, temperature, and proton and alpha-particle properties in the inner heliosphere. The observations presented here provide the starting point for a better understanding of the origin and evolution of different solar wind streams as well as the evolving turbulence contained within.

magnetohydrodynamics (MHD)↗

Flux Cube Reconstruction from Slitless Spectroscopy

Slitless spectroscopy enables efficient, large-area surveys without target preselection, yet it faces challenges from source blending, higher noise, and lost spatial–spectral information. We present an advanced, nonparametric, data-driven algorithm that leverages multiple dispersion angles to reconstruct three-dimensional flux distributions, providing low-resolution integral field unit capabilities from slitless data. By treating each pixel as an independent element, our method naturally handles source confusion without requiring prior assumptions regarding redshifts, templates, or model libraries. We validate the algorithm using simulated Roman Space Telescope wide-field slitless spectroscopy images that are equivalent to what is expected from the High-Latitude Time-Domain Survey. First, we demonstrate that a host-galaxy model reconstructed from multiple dispersion angles can be used to accurately subtract host light from a transient, recovering a Type Ia supernova spectrum with minimal bias. Second, we showcase a high-fidelity flux-cube reconstruction of a complex galaxy, successfully measuring the redshift and recovering continuum, emission, and absorption features. This approach highlights the potential of multi-dispersion-angle slitless data to provide spatially resolved spectral information in a nonparametric way, which is traditionally accessible only with integral field spectroscopy, opening a new window into large, unbiased, and spatially resolved studies of galaxy evolution.

Griggio, M. [Space Telescope Science Institute, Ba↗

Unprecedented Beaufort Sea ice loss in late summer 2021 and its relationship to an extended period of unusually stormy weather

Previous case studies have linked cyclone-induced atmospheric forcing and/or upper-ocean processes to notable Arctic sea ice loss events in the summers of 2012 and 2016. This study examines a more recent and noteworthy case in late summer 2021 in which substantial sea ice loss followed a period of surface meteorological extremes in the Beaufort Sea region of the Arctic. We focus on the period from mid-August to mid-September 2021 that coincided with the Office of Naval Research THINICE Pilot Field Campaign and investigate stormy and windy conditions with respect to air-sea processes impacting sea ice conditions. We find that during the stormy first half of the campaign, cyclone-induced energy fluxes into the marginal ice zone and surrounding waters preconditioned the ice pack for more rapid melt later in the campaign. The second half of the campaign, in contrast, was marked by non-cyclone wind events that enhanced turbulent (namely sensible) heat fluxes into the ice and upper ocean that increased melt. Moreover, this latter period had enhanced advection of the Beaufort Sea ice pack into above-freezing waters, increasing bottom melt to >1 cm d −1 over the remainder of the campaign. While findings are shown to vary by surface type and at relatively small (i.e. ice-floe) scales, insights are offered on the roles of late summer coupled processes on rapid ice loss events in today’s Arctic environment.

54 ENVIRONMENTAL SCIENCES↗

NASA’s Pandora SmallSat Mission: Simulating the Impact of Stellar Photospheric Heterogeneity and Its Correction

Stellar photospheric heterogeneity is a dominant astrophysical systematic impacting exoplanet transmission spectroscopy. NASA’s Pandora SmallSat Mission is designed to address this challenge through contemporaneous visible-band photometry and near-infrared spectroscopy of exoplanet host stars. Here, we present an end-to-end simulation study quantifying Pandora’s ability to infer stellar photospheric properties and correct stellar contamination using out-of-transit observations. We construct eight representative stellar activity scenarios and generate 160 simulated Pandora datasets, incorporating time-dependent stellar spectra, instrument response, and noise. Given accurate models, Bayesian retrievals of joint visible photometry (0.4–0.7 μm) and near-infrared spectroscopy (0.9–1.6 μm, R ≈ 120) recover photospheric temperatures with typical uncertainties of ≈30 K, with no significant bias. Models with two spectral components (i.e., a quiescent photosphere and spots) are strongly favored in 95% of cases; one-component models are preferred when true spot filling factors fall below a detection threshold of ≈0.3%. We propagate the true and inferred stellar parameters to compute true, inferred, and residual contamination signals under physically motivated spot geometries. For simple spot distributions, contamination signals of 10 2 –10 3 ppm are reduced to ≲10 ppm—well below Pandora’s expected transmission spectroscopy precision (30–100 ppm). For more complex spot distributions, geometric degeneracies limit deterministic corrections, leaving residual contamination at the 10 3 ppm level that must be mitigated using additional constraints, such as spot-crossing events and joint stellar–planetary retrievals of transmission spectra. These results define regimes in which stellar contamination can be corrected from stellar observations alone and show how Pandora stellar observations can identify cases where additional information is required.

Astronomy and AstroPhysics↗

Hot Hydrogen Exposure of U x Zr 1-x C y Nuclear Fuel: The Effect of Additional Carbon

Nuclear thermal propulsion (NTP) using hydrogen as propellant in a solid-state fission reactor to reach temperatures of up to 3000K can achieve significantly higher specific impulse than chemical propulsion. A promising fuel that is again considered for NTP is U x Zr 1-x C y . To reduce the overall mass loss as well as the uranium fuel loss, a range of experiments using additional carbon in hot hydrogen (2600 K, 6 SLPM flow rate) were performed for up to 5-6 h. Both a CH 4 addition (0.1 and 0.2 vol%) in the hydrogen stream and the use of sacrificial graphite pieces upstream of the samples were tested, with sample compositions from 5 at% UC to 30 at% UC. The results showed a strong reduction of mass loss, up to a factor of 3, as well as a reduction of surface uranium losses in the presence of additional carbon in the atmosphere. Even with additional carbon in the atmosphere, material with 30 at% UC was not stable, but material with 20 at% UC was stable for up to 6h. The sample surfaces showed texturing after processing, possibly by grain growth on the surface. This led to a separation into grains with high uranium content close to the (100)-orientation, and with low uranium content in grains close to the (111)-orientation.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Beam test performance of AstroPix sensor with 120 GeV protons

AstroPix is a High-Voltage CMOS Monolithic Active Pixel Sensor (HV-CMOS MAPS) developed for precision gamma-ray imaging and spectroscopy in the medium-energy regime, as well as for precise shower imaging and tracking in the Barrel Imaging Calorimeter (BIC) of the Electron Proton/Ion Collider (ePIC) detector at the future Electron–Ion Collider (EIC). We present beam test results of the AstroPix_v3 sensor using a 120 GeV proton beam at the Fermilab Test Beam Facility (FTBF), performed as part of the broader experimental campaign for the BIC prototype calorimeter. The sensor’s 500 µm pixel pitch enabled precise measurement of the beam profile, providing important information for the calorimeter performance studies. Using the measured 120 GeV proton data, we measure the energy deposit of minimum ionizing particles (MIP) and use them to extract the corresponding effective depletion depth at a single bias voltage of −150 V .

AstroPix↗

Hot Hydrogen Exposure of U x Zr 1-x C y Nuclear Fuel: The Influence of Composition and Density

Refractory carbide nuclear fuel has been one of the most promising fuel candidates for space nuclear propulsion due to its high melting point, temperature stability, and compatibility in a hot hydrogen environment. In this study, U x Zr 1-x C y fuel was produced by means of a carbothermic reduction process in different UC compositions including 5,10, 20, and 30 at.% UC in the fuel compound. The powder feedstock was consolidated via direct current sintering with densities up to 97% of the theoretical density. The samples with different U x Zr 1-x C y compositions were exposed to hot hydrogen at 2600 K for a cumulative time of 300 min. The samples were characterized by SEM, XRD, density, and measured for mass losses. The high-density samples displayed improved performance in hot hydrogen by minimizing porous sites and reducing areas of direct contact with hydrogen gas, leading to reduced mass losses. Variations in sample density proved to induce large changes in mass loss rates, increasing them up to 90%. The compositions with higher UC content reported the largest mass losses in the study. The loss of uranium occurred primarily at the surfaces exposed to the hot hydrogen where changes in the lattice constant confirmed losses exceeding 50% of the initial UC content in higher compositions, specifically to U 0.3 Zr 0.7 C y . XRD analyses revealed the presence of UH 3 in U 0.3 Zr 0.7 C y suggesting that metallic uranium formed inside the sample as a product of carbon losses. High-density U x Zr x-1 C y fuel with UC concentrations at or below 20 at.% UC exhibited stability and negligible density changes in a high temperature hydrogen environment.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

High-count-rate effects in event processing for the XRISM/Resolve X-ray microcalorimeter. II. Energy scale and resolution in orbit

The Resolve instrument on the X-ray Imaging and Spectroscopy Mission (XRISM) uses a 36 pixel microcalorimeter designed to deliver high-resolution, non-dispersive X-ray spectroscopy. Although it is optimized for extended sources with low count rates, Resolve observations of bright point sources are still able to provide unique insights into the physics of these objects, as long as high-count-rate effects are addressed in the analysis. These effects include the loss of exposure time for each pixel, changes in the energy scale, and changes in the energy resolution. To investigate these effects under realistic observational conditions, we observed the bright X-ray source, the Crab Nebula, with XRISM at several offset positions with respect to the Resolve field of view and with continuous illumination from 55 Fe sources on the filter wheel. For the spectral analysis, we excluded data where exposure-time loss was too significant to ensure reliable spectral statistics. The energy scale at 6 keV shows a slight negative shift in the high-count-rate regime. The energy resolution at 6 keV worsens as the count rate in electrically neighboring pixels increases, but can be restored by applying a nearest-neighbor coincidence cut (“cross-talk cut”). We examined how these effects influence the observation of bright point sources, using GX 13+1 as a test case, and identified an eV-scale energy offset at 6 keV between the inner (brighter) and outer (fainter) pixels. Users who seek to analyze velocity structures on the order of tens of km s–1 should account for such high-count-rate effects. These findings will aid in the interpretation of Resolve data from bright sources and provide valuable considerations for designing and planning for future microcalorimeter missions.

X-rays: general↗

Resolve instrument onboard XRISM: design, integration, and instrument test results

The Resolve instrument onboard the X-Ray Imaging and Spectroscopy Mission (XRISM) consists of an array of 6 × 6 silicon-thermistor microcalorimeters cooled down to 50 mK and a high-throughput X-ray mirror assembly (XMA) with a focal length of 5.6 m. XRISM is a recovery mission of ASTRO-H/Hitomi, and the Resolve instrument is a rebuild of the ASTRO-H Soft X-ray spectrometer (SXS) and the Soft X-ray Telescope (SXT) that achieved energy resolution of ∼ 5 eV FWHM on orbit, with several important changes based on lessons learned from ASTRO-H. The flight models of the Dewar and the electronics boxes were fabricated, and the instrument test and calibration were conducted in 2021. By tuning the cryocooler frequencies, energy resolution better than 4.9 eV FWHM at 6 keV was demonstrated for all 36 pixels and high-resolution grade events, as well as energy-scale accuracy better than 2 eV up to 30 keV. The immunity of the detectors to microvibration, electrical conduction, and radiation was evaluated. The instrument was delivered to the spacecraft system in April 2022. The XMA was tested and calibrated separately. Its angular resolution is 1.27′, and the effective area of the mirror itself is 570 cm 2 at 1 keV and 424 cm 2 at 6 keV. We report the design and the major changes from the ASTRO-H SXS, the integration, and the results of the instrument test.

X-ray↗

Linking Lattice Strain and Fractal Dimensions to Non‐monotonic Volume Changes in Irradiated Nuclear Graphite

Graphite's resilience to high temperatures and neutron damage makes it vital for nuclear reactors, yet irradiation alters its microstructure, degrading key properties. We used small- and wide-angle X-ray scattering to study neutron-irradiated fine-grain nuclear graphite (Grade G347A) across varied temperatures and fluences. Results show significant shifts in internal strain and porosity, correlating with radiation-induced volume changes. Notably, porosity volume distribution (fractal dimensions) follows non-monotonic volume changes, suggesting a link to the Weibull distribution of fracture stress.

X-ray scattering↗

Structural Design of Bismuth Telluride Nanoplates through Process Variables

Binary pnictogen chalcogen compounds, primarily bismuth tellurides and selenides, are of great interest due to their applications in emerging quantum devices, as well as thermoelectric generators. The performance of bismuth telluride in these roles depends on its structure at the nanoscale, particularly the size, shape, and crystallinity of its nanocrystalline forms. However, current methods for controlling these features are often slow, inconsistent, or difficult to scale. Here, we demonstrate that through a solvothermal synthesis and hot injection process, precise control over the morphology of bismuth telluride nanoplates is possible with independent tuning of process variables, such as temperature and reaction time. We find that the nanoplate shape and internal porosity vary systematically with synthesis temperature and that the same morphological outcomes can be rapidly achieved at a fixed temperature by adjusting reaction duration. These results reveal that both the temperature and time can independently direct bismuth telluride morphological features, allowing for rapid, tunable synthesis strategies. Our approach offers a scalable framework, not only for bismuth telluride but also for related layered chalcogenides used in energy harvesting and quantum technologies.

Ackley, Jordan [Boise State Univ., ID (United Stat↗

Restoring Historic Forest Disturbance Frequency Would Partially Mitigate Droughts in the Central Sierra Nevada Mountains

Forest thinning and prescribed fire are expected to improve the climate resilience and water security of forests in the western U.S., but few studies have directly modeled the hydrological effects of multi-decadal landscape-scale forest disturbance. By updating a distributed process-based hydrological model (DHSVM) with vegetation maps from a distributed forest ecosystem model (LANDIS-II), we simulate the water resource impacts of forest management scenarios targeting partial or full restoration of the pre-colonial disturbance return interval in the central Sierra Nevada mountains. In a fully restored disturbance regime that includes fire, thinning, and insect mortality, reservoir inflow increases by 4%–9% total and 8%–14% in dry years. At sub-watershed scales (10–100 km2), thinning dense forests can increase streamflow by >20% in dry years. In a thinner forest, increased understory transpiration compensates for decreased overstory transpiration. Consequentially, 73% of streamflow gains are attributable to decreased overstory rain and snow interception loss. Thinner forests can increase headwater peak flows, but reservoir-scale peak flows are almost exclusively influenced by climate. Uncertainty in future precipitation causes high uncertainty in future water yield, but the additional water yield attributable to forest disturbance is about five times less sensitive to annual precipitation uncertainty. This partial decoupling of the streamflow disturbance response from annual precipitation makes disturbance especially valuable for water supply during dry years. Our study can increase confidence in the water resource benefits of restoring historic forest disturbance frequencies in the central Sierra Nevada mountains, and our modeling framework is widely applicable to other forested mountain landscapes.

Boardman, Eli N. [University of Nevada, Reno, NV (↗

A Framework for Parametric and Predictive Uncertainty Quantification in the E3SM Land Model: Assessing Site and Observable Generalizability

Quantifying parametric uncertainty using observations from individual sites provides a critical foundation for Earth system modeling, serving as a necessary first step before scaling up to regional or global applications. This study introduces a novel computational framework designed to enhance model predictability by reducing parametric uncertainty and assessing site and observable generalizability using various observational constraints. The framework integrates five components: Model Simulation, Statistical Emulation, Global Sensitivity Analysis (GSA), Model Calibration, and Model Prediction. Using the E3SM land model, we simulated site-level land-atmosphere carbon and energy fluxes from 2003 to 2007 across five evergreen needleleaf FLUXNET sites, perturbing 26 vegetation-related model parameters. Gaussian process emulators were employed to expedite GSA and model calibration. Four critical parameters that strongly influence selected land-atmosphere fluxes were identified by GSA. Bayesian approaches were used to infer parameter probability distributions leveraging synthetic data and FLUXNET observations. The results reveal that posterior parameter distributions vary significantly across different sites and observables within the same plant functional type. Probabilistic predictions indicate that parameters calibrated at one site can enhance predictive accuracy at other sites, although site heterogeneity may sometimes outweigh parametric uncertainty. Additionally, the probabilistic predictions demonstrate that calibration for one variable can also improve predictability for other variables, thereby maximizing predictive capabilities with limited observations. This framework provides a powerful approach for reducing parametric uncertainty in Earth system models and deepening our understanding of carbon dynamics and energy cycles. Its adaptability makes it a valuable tool for broader applications in Earth system modeling.

54 ENVIRONMENTAL SCIENCES↗

New particle pusher with hadronic interactions for modeling multimessenger emission from compact objects

We propose novel numerical schemes based on the Boris method in curved spacetime, incorporating both hadronic and radiative interactions for the first time. Once the proton has lost significant energy due to radiative and hadronic losses, and its gyroradius has decreased below typical scales on which the electromagnetic field varies, we apply a guiding center approximation (GCA). We fundamentally simulate collision processes either with a Monte-Carlo method or, where applicable, as a continuous energy loss, contingent on the local optical depth. To test our algorithm for the first time combining the effects of electromagnetic, gravitational, and radiation fields including hadronic interactions, we simulate highly relativistic protons traveling through various electromagnetic fields and proton backgrounds. We provide unit tests in various spatially dependent electromagnetic and gravitational fields and background photon and proton distributions, comparing the trajectory against analytic results. We propose that our method can be used to analyze hadronic interactions in black hole accretion disks, jets, and coronae to study the neutrino abundance from active galactic nuclei.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗