Search NASA⌕ Search

SEARCH · Search NASA

Results for “Kinetic methods”

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 577 records · Page 32
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Hydrodynamic characterization of the redox chemistry of crown-encapsulated uranyl complexes

The redox properties of actinide-containing species strongly influence their reactivity, speciation, and interfacial behavior, but the experimental quantification of the electrochemical characteristics of molecular actinide complexes in nonaqueous media has not received the attention it deserves. Here, results from hydrodynamic methods and electrochemical simulations of U(VI)/U(V) redox are reported, including quantification of heterogeneous electron-transfer kinetics and estimation of chemical reversibility of U(VI)/U(V) interconversion at electrodes in acetonitrile-based electrolyte. The complexes investigated are recently reported U(VI) and U(V) complexes in which the uranyl ion (UO 2 n+ ) is encapsulated in a macrocyclic 18-crown-6-like moiety templated by a Pt(II) center. These complexes feature the most positive value U VI /U V reduction potential yet reported and are thus particularly relevant to study of facile U(V) generation from U(VI) precursors as well as uranium electroanalysis. Rotating disk electrode (RDE) studies have been used to quantify the diffusion coefficients of the U(VI) and U(V) complexes, and standard heterogeneous electron transfer rate constants ( k 0 ) for the redox have been determined using a conventional Koutecký-Levich analysis. Rotating ring-disk electrode (RRDE) studies have been used to directly interrogate the chemical reversibility of U(VI)-U(V) interconversion, confirming that reduction of the U(VI) complex at an Au disk is associated with formation of the U(V) analogue that can be readily re-oxidized at a Pt ring under hydrodynamic (rotating) conditions. Because measurements of the type reported here are generally associated with current flows that are larger than those found in corresponding quiescent (unstirred) conditions, our findings suggest that hydrodynamic methods could be advantageous for design of electroanalytical approaches to detection of actinide species and study of their redox properties.

actinides↗
DOE OSTI2024

Dynamical defects in a two-dimensional Wigner crystal: Self-doping and kinetic magnetism

We study the quantum dynamics of interstitials and vacancies in a two-dimensional Wigner crystal (WC) using a semi-classical instanton method that is asymptotically exact at low density, i.e., in the r s → ∞ limit. Here, the dynamics of these point defects mediates magnetism with much higher energy scales than the exchange energies of the pure WC. Via exact diagonalization of the derived effective Hamiltonians in the single-defect sectors, we find the dynamical corrections to the defect energies. The resulting expression for the interstitial (vacancy) energy extrapolates to 0 at r s = r mit ≈ 70 (r s ≈ 30), suggestive of a self-doping instability to a partially melted WC for some range of r s below r mit . We thus propose a “metallic electron crystal” phase of the two-dimensional electron gas at intermediate densities between a low density insulating WC and a high density Fermi fluid.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1990

Production of oxygen from lunar ilmenite

The following subjects are addressed: (1) the mechanism and kinetics of carbothermal reduction of simulated lunar ilmenite using carbon and, particularly, CO as reducing agents; (2) the determination of the rate-limiting steps; (3) the investigation of the effect of impurities, particularly magnesium; (4) the search for catalysts suitable for enhancement of the rate-limiting step; (5) the comparison of the kinetics of carbothermal reduction with those of hydrogen reduction; (6) the study of the combined use of CO and hydrogen as products of gasification of carbonaceous solids; (7) the development of reduction methods based on the use of waste carbonaceous compounds for the process; (8) the development of a carbothermal reaction path that utilizes gasification of carbonaceous solids to reducing gaseous species (hydrocarbons and carbon monoxide) to facilitate the reduction reaction kinetics and make the process more flexible in using various forms of carbonaceous feeds; (9) the development of advanced gas separation techniques, including the use of high-temperature ceramic membranes; (10) the development of an optimum process flow sheet for carbothermal reduction, and comparison of this process with the hydrogen reduction scheme, as well as a general comparison with other leading oxygen production schemes; and (11) the use of new and advanced material processing and separation techniques.

Zhao, Y.↗
NASA NTRS2005

A Kinetic and Product Study of the Cl + HO2 Reaction

Absolute rate data and product branching ratios for the reactions Cl + HO2 to HCl + O2 (k1a) and Cl + HO2 to OH + ClO (k1b) have been measured from 226 to 336 K at a total pressure of 1 Torr of helium using the discharge flow resonance fluorescence technique coupled with infrared diode laser spectroscopy. For kinetic measurements, pseudo-first-order conditions were used with both reagents in excess in separate experiments. HO2 was produced by two methods: through the termolecular reaction of H atoms with O2 and also by the reaction of F atoms with H2O2. Cl atoms were produced by a microwave discharge of Cl2 in He. HO2 radicals were converted to OH radicals prior to detection by resonance fluorescence at 308 nm. Cl atoms were detected directly at 138 nm also by resonance fluorescence. Measurement of the consumption of HO2 in excess Cl yielded k1a and measurement of the consumption of Cl in excess HO2 yielded the total rate coefficient, k1. Values of k1a and k1 derived from kinetic experiments expressed in Arrhenius form are (1.6 +/- 0.2) x 10-11 exp[(249 +/- 34)/T] and (2.8 +/- 0.1) x 10-11 exp[(123 +/- 15)/T] cm3 molecule-1 s-1, respectively. As the expression for k1 is only weakly temperature dependent, we report a temperature-independent value of k1 = (4.5 +/- 0.4) x 10-11 cm3 molecule-1 s-1. Additionally, an Arrhenius expression for k1b can also be derived: k1b = (7.7 +/- 0.8) x 10-11 exp[-(708 +/- 29)/T] cm3 molecule-1 s-1. These expressions for k1a and k1b are valid for 226 K T 336 and 256 K T 296 K, respectively. The cited errors are at the level of a single standard deviation. For the product measurements, an excess of Cl was added to known concentrations of HO2 and the reaction was allowed to reach completion. HCl product concentrations were determined by IR absorption yielding the ratio k1a/k1 over the temperature range 236 K T 296 K. OH product concentrations were determined by resonance fluorescence giving rise to the ratio k1b/k1 over the temperature range 226 K T 336 K. Both of these ratios were subsequently converted to absolute numbers. Values of k1a and k1b from the product experiments expressed in Arrhenius form are (1.5 +/- 0.1) x 10-11 exp[(222 +/- 17)/T] and (10.6 +/- 1.5) x 10-11 exp[-(733 +/- 41)/T] cm3 molecule-1 s-1, respectively. These expressions for k1a and k1b are valid for 256 K T 296 and 226 K T 336 K, respectively. A combination of the kinetic and product data results in the following Arrhenius expressions for k1a and k1b of (1.4 +/- 0.3) x 10-11 exp[(269 +/- 58)/T] and (12.7 +/- 4.1) x 10-11 exp[-(801 +/- 94)/T] cm3 molecule-1 s-1, respectively. Numerical simulations were used to check for interferences from secondary chemistry in both the kinetic and product experiments and also to quantify the losses incurred during the conversion process HO2 to OH for detection purposes.

rate constants↗
NASA NTRS1996

Investigation of ion-beam machining methods for replicated x-ray optics

The final figuring step in the fabrication of an optical component involves imparting a specified contour onto the surface. This can be expensive and time consuming step. The recent development of ion beam figuring provides a method for performing the figuring process with advantages over standard mechanical methods. Ion figuring has proven effective in figuring large optical components. The process of ion beam figuring removes material by transferring kinetic energy from impinging neutral particles. The process utilizes a Kaufman type ion source, where a plasma is generated in a discharge chamber by controlled electric potentials. Charged grids extract and accelerate ions from the chamber. The accelerated ions form a directional beam. A neutralizer outside the accelerator grids supplies electrons to the positive ion beam. It is necessary to neutralize the beam to prevent charging workpieces and to avoid bending the beam with extraneous electro-magnetic fields. When the directed beam strikes the workpiece, material sputters in a predicable manner. The amount and distribution of material sputtered is a function of the energy of the beam, material of the component, distance from the workpiece, and angle of incidence of the beam. The figuring method described here assumes a constant beam removal, so that the process can be represented by a convolution operation. A fixed beam energy maintains a constant sputtering rate. This temporally and spatially stable beam is held perpendicular to the workpiece at a fixed distance. For non-constant removal, corrections would be required to model the process as a convolution operation. Specific figures (contours) are achieved by rastering the beam over the workpiece at varying velocities. A unique deconvolution is performed, using series-derivative solution developed for the system, to determine these velocities.

Drueding, Thomas W.↗
DOE OSTI2025

SPARTA: A flux adjustment methodology to interpret complex experiments

For the accurate determination of reactivity from a detector count rate, correction of spatial effects is of prime importance. This spatial correction is often provided using simulation methodologies, but this may introduce a bias if the result of the experiment is also used as input data for the simulation. Here, this work presents a flux adjustment methodology able to infer experimental reactivity and correction of spatial effects without the need for a simulation. It can process the signal from a complex experiment such as a heat balance measurement in the TREAT reactor, where control rods are continuously adjusted to maintain a constant power. In the present work, this methodology successfully computed the reactivity and the local spatial variation of the flux of a generated signal. It also proved to be robust against noise and errors on kinetic parameters and provides a credible interpretation of a heat balance experiment in TREAT. Efficiency of flux adjustment methods for complex experiment enable a better experiment interpretation less reliant on nuclear data evaluation.

73 - NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Spinel high-entropy oxides (FeNiCrMnZnX) 3 O 4 (X = Al, mg) as anode materials for high-performance lithium-ion batteries

To address the high cost, cobalt dependency, and resource constraints typical of conventional high-entropy oxide (HEO) anodes, this study reports the successful synthesis of two Co-free, six-component spinel-type HEOs(FeNiCrMnZnAl) 3 O 4 (HEO-Al) and (FeNiCrMnZnMg) 3 O 4 (HEO-Mg), via a sol-gel method. The distinct effects of Al 3+ and Mg 2+ incorporation on the electrochemical performance and lithium storage kinetics were systematically investigated. XRD, Raman, and TEM characterizations confirm that both materials possess a pure spinel phase, uniform particle size, and homogeneous elemental distribution. Notably, electrochemical evaluations reveal that HEO-Al delivers a superior reversible capacity of 480.7 mAh g −1 after 100 cycles at 0.1 A g −1 , and maintains 354.5 mAh g −1 after 1000 long-term cycles at 1 A g −1 , significantly outperforming HEO-Mg. Kinetic analysis indicates that HEO-Al exhibits lower charge transfer resistance, a higher Li + diffusion coefficient, and a pseudocapacitive contribution of up to 82%. Furthermore, these findings demonstrate that Al substitution effectively optimizes the structural stability and lithium storage kinetics of Co-free HEOs, providing a viable strategy for designing low-cost, highly stable HEO anode systems.

Anode materials↗
DOE OSTI2026

Preparation of a uranium monocarbide anode and electrochemical characterization in molten LiCl-KCl-UCl 3

Porous uranium carbide (UC) pellets possessing moderate electrical conductivity were synthesized by reaction of UO 2 with graphite at temperatures up to 1550°C under rough vacuum. Conversions as high as 98% were achieved at soak times of 2-4 hours. The electrochemistry of the UC pellets in molten LiCl-KCl-6.5 wt% UCl 3 was explored using a variety of techniques including DC polarization methods, cyclic voltammetry, chronopotentiometry and bulk electrolysis. Here, the electrode reaction for anodic dissolution was found to be kinetically controlled by dissociation of UC to a transition state complex that was hypothesized to consist of a uranium atom partially complexed by chloride ions. Precise measurements of current efficiencies using chronopotentiometry indicated upper limits of 90.9 ± 3.4% and 98.3 +1.7/-3.7% for anode and cathode, respectively, when operating at anodic overpotentials near +300 mV. Bulk electrolysis of a UC pellet performed by passing 98% of the theoretical charge resulted in nearly complete recovery of its uranium content as highly pure metal at the cathode.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Redox Couples Control Band Bending, Photovoltage, and Quasi-Fermi Levels in Tungsten Oxide (WO 3 ) Photoanodes

Tungsten oxide (WO 3 ) is a well-known photoanode and photocatalyst for photoelectrochemical (PEC) water oxidation. Because the compound has a deep valence band, it can facilitate the oxygen evolution reaction without added cocatalysts, and it can drive the oxidation of species with much higher electrochemical potentials, including the conversion of water to hydrogen peroxide, sulfate to persulfate, and iodate to meta-periodate. Here, we use the liquid vibrating Kelvin probe surface photovoltage (liquid VK-SPV) technique in combination with open circuit potential (OCP) and photoelectrochemical (PEC) scans to assess the possibility of reaching such oxidizing potentials in aqueous electrolytes and at open circuit. Here, this is done by mapping the quasi-Fermi levels of electrons and holes at the interfaces as a function of the light intensity. Nanostructured WO 3 photoelectrodes for this purpose were fabricated by thermal annealing of a tungstic acid solution on fluorine-doped tin oxide. Electrochemical measurements are conducted at open circuit and 400 nm LED light illumination in electrolytes containing fast (O 2 /H 2 O 2 ), slow (O 2 /H 2 O), and very oxidizing (NaIO 4 /NaIO 3 ) redox couples. Photovoltage values scale with the light intensity and with the built-in potential for each redox couple and reach values up to 0.61 V under 20 mW cm –2 illumination for the NaIO 4 electrolyte. This shows that the photoelectrodes behave like Schottky-type diodes whose maximum possible energy output is determined mainly by the built-in voltage of each junction. For slow redox couples, the quasi-Fermi level of the holes increases with light intensity due to hole accumulation at the WO 3 –liquid interface. For example, for the O 2 /H 2 O electrolyte, interfacial hole accumulation and removal occur on the 90–300 s time scale. For the fast hole acceptor H 2 O 2 , on the other hand, the quasi-Fermi level of the photoholes is pinned to the electrochemical potential of the O 2 /H 2 O 2 couple. This limits the energy conversion efficiency of the electrode. Overall, these results reveal the influence of charge transfer thermodynamics and kinetics on the photovoltage of WO 3 . Furthermore, the work further establishes VK-SPV as a contactless method to observe the photovoltage, carrier dynamics, and quasi-Fermi levels of semiconductor-liquid junctions.

Electrodes↗
DOE OSTI2026

Transient Pulse-Response Time-of-Flight Mass Spectrometry for Complex, Deactivating Heterogeneous Catalytic Systems: Application to Ethane Dehydroaromatization

The study of complex, multistep bond-forming and -breaking reactions in heterogeneous catalytic systems often encounters challenges associated with the involvement of large numbers of intermediates among branching pathways. Kinetic information obtained from traditional steady-state measurements can be complemented with that from time-resolved methods to uncover details of the underlying chemistry. Herein, we describe an approach for tracking the complete time-resolved chemical composition (ca. 4–200 u) of a reactor effluent in response to a reactant pulse. We use a six-port rotary valve with a metered sampling loop to pulse reactants at ambient pressure into a flow reactor packed with a catalyst bed within the isothermal region of a heated furnace. The temporal evolution of effluent species is tracked using time-resolved molecular-beam time-of-flight mass spectrometry. We highlight the possibilities that this method has to offer by studying the complex bifunctional mechanism of ethane dehydroaromatization over an HZSM-5-supported platinum catalyst. We demonstrate that energy-tunable ionization sources, which facilitate isomer resolution, enable the measurement of the full mass spectral time-dependent system response. This includes the evolution of major products and mechanistically relevant reactive intermediates such as 1,3-butadiene and cyclopentadiene; these species have not previously been observed from this reaction. In additional studies, we also assess the role of platinum in the catalyst by examining temporal responses to ethane and ethylene feeds. Results show that two temporally distinct formation pathways exist for methane and benzene, and that their importance depends on both catalyst composition and reactant identity. Additionally, characteristics of catalyst deactivation are uniquely observable in the time-resolved mass spectral response, including the selective deactivation of a benzene formation pathway. The combination of time-of-flight mass spectrometry with tunable ionization enables the simultaneous observation of all effluents in a complex mixture of intermediates with isomer/isobar differentiation capabilities that can be applied to any complex reaction system.

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

Institute for Catalysis in Energy Processes (ICEP) (Final Report)

The Institute for Catalysis in Energy Processes (ICEP) was a multi-PI program located at the Northwestern University Center for Catalysis and Surface Science from 2006-2024. Over several renewal cycles and several organizing themes, ICEP addressed fundamental questions in catalysis science. In turn, the scientific questions addressed were directly relevant to the efficient use of the nation’s resources to produce fuels and commodity chemicals, harness alternate energy sources for chemical reactions such as light or electricity, reduce emissions and waste, and minimize the impact of our use of plastics. Selective oxidation, (oxidative) dehydrogenation, deNOx, CO 2 reduction, hydrogen release, polymer decomposition / depolymerization, and many other reactions of intense interest to the US Department of Energy were studied in the center. Major ICEP strengths were in catalyst synthesis, in measurements that elucidated catalyst properties, in reaction mechanisms and kinetics, and in predictive theory and modeling. ICEP researchers were often the inventors or developers of materials and methods that were brought to bear on the Institute’s catalytic systems. Key characterization tools advanced through this project included resonance Raman and related techniques, sum frequency generation, and X-ray standing wave spectroscopy. Center members relied extensively on computational tools such as density functional theory and microkinetic modeling to fully investigate catalyst structures and catalytic mechanisms. ICEP researchers were early developers of atomic layer deposition (ALD) for catalyst synthesis and were early pioneers of metal organic frameworks for catalysis. ICEP innovations are found in areas of intense research such as single-atom catalysts, catalytic deconstruction of plastics, and catalytic MOFs, to name a few. Research areas initiated by ICEP indirectly led to DOE EFRCs, startup companies, and other achievements. Over its 18 years of existence, ICEP involved two different PIs (Stair and Notestein) and 24 senior investigators across several academic disciplines at Northwestern University. There were tight collaborations with Argonne National Laboratory, including specialized equipment and experiments permanently located there. The center supported ~400 person-years of effort by graduate students and postdocs, either directly through DOE support or indirectly by leveraging independent support provided to the trainees, e.g. through the NSF GRFP or internal fellowships. The project resulted in 328 manuscripts and numerous conference presentations.

36 MATERIALS SCIENCE↗
NASA NTRS1974

Isothermal transitions of a thermosetting system

A study of the curing reactions of a cycloaliphatic epoxy resin/anhydride system by torsional braid analysis showed the existence of two critical isothermal temperatures - namely, the maximum glass transition temperature of the thermoset system and the glass transition temperature of the material at its gel point. Two rheologically active kinetic transitions occur during isothermal cure which correspond to gelation and vitrification. Three types of isothermal behavior occur. Methods for determining the time to gel and the time to vitrify, and also the two above-mentioned critical isothermal temperatures, have been developed. The time to gel obeyed the Arrhenius relationship, whereas the time to vitrify passed through a minimum. Application of these results to thermosetting systems in general is discussed in terms of the influence of molecular structure on the values of the critical isothermal temperatures.

Gillham, J. K.↗
NASA NTRS1982

A study of the photocatalytic effects of aqueous suspensions of platinized semiconductor materials on the reaction rates of candidate redox reactions

The effectiveness of powdered semiconductor materials in photocatalyzing candidate redox reactions was investigated. The rate of the photocatalyzed oxidation of cyanide at platinized TiO2 was studied. The extent of the cyanide reaction was followed directly using an electroanalytical method (i.e. differential pulse polarography). Experiments were performed in natural or artificial light. A comparison was made of kinetic data obtained for photocatalysis at platinized powders with rate data for nonplatinized powders.

Miles, A. M.↗
NASA NTRS1988

A computationally efficient modelling of laminar separation bubbles

The goal of this research is to accurately predict the characteristics of the laminar separation bubble and its effects on airfoil performance. To this end, a model of the bubble is under development and will be incorporated in the analysis section of the Eppler and Somers program. As a first step in this direction, an existing bubble model was inserted into the program. It was decided to address the problem of the short bubble before attempting the prediction of the long bubble. In the second place, an integral boundary-layer method is believed more desirable than a finite difference approach. While these two methods achieve similar prediction accuracy, finite-difference methods tend to involve significantly longer computer run times than the integral methods. Finally, as the boundary-layer analysis in the Eppler and Somers program employs the momentum and kinetic energy integral equations, a short-bubble model compatible with these equations is most preferable.

Maughmer, Mark D.↗
NASA NTRS1989

Simulation of collisional transport processes and the stability of planetary rings

The utility of the phase-space fluid method for the study of planetary ring dynamics is presently demonstrated through the numerical solution of a model kinetic equation for a flattened Keplerian disk. Attention is given to ringlets composed of single-sized particles, as well as to ringlets composed of two different-sized particles; in the latter case, the ringlets evolve in such a way that the lighter particles are confined by the heavier ones. The results obtained indicate that some natural process may sharpen the optical depth profile of edges even without an external forcing mechanism, and that intermediate optical depths are dynamically preferred in some cases.

Brophy, Thomas G.↗
NASA NTRS1992

Middle and upper atmosphere pressure-temperature profiles and the abundances of CO2 and CO in the upper atmosphere from ATMOS/Spacelab 3 observations

An improved method for retrieving pressure-temperature profiles is described and is used to retrieve profiles of the kinetic-temperature and atmospheric-pressure profiles between 20 and 116 km altitudes and the CO2 and CO volume-mixing ratios between 70 and 116 km, using the IR occultation spectra recorded by the Spacelab 3 atmospheric trace molecular spectroscopy (ATMOS) Fourier transform spectrometer between April 29 and May 6, 1985. Profiles are derived for six ATMOS occultations. The CO2 and CO volume-mixing profiles are compared with previous observations and model predictions. Evidence is found for vibrational non-LTE by analyzing the lines of the (nu-2 + nu-3 - nu-2) (C-12)(O-16) band. Results are used for deriving (C-12)(O-16) (010) vibrational temperatures, which are compared with the retrieved kinetic temperatures and the predictions of non-LTE effects by recent models.

Rinsland, C. P.↗
NASA NTRS1995

Bounce-Averaged Hamiltonian for Charged Particles in an Axisymmetric but Nondipolar Model Magnetosphere

In order to facilitate bounce-averaged guiding center simulations of geomagnetically trapped particles, we express the kinetic energy of a particle with magnetic coordinates (L,phi) as an analytic function of the first two adiabatic invariants (M, J) and the L value of the field line. The magnetic field model is axisymmetric, consisting of a dipolar B field plus a uniform southward magnetic field parallel to the dipole moment mu(sub E). This model magnetosphere is surrounded by a circular equatorial neutral line whose radius b is an adjustable parameter. The L value of a field line is (by definition) inversely proportional to the flux enclosed by the corresponding magnetic shell of equatorial radius r(sub 0), and the L value at the neutral line (r(sub 0) = b) is denoted L*. The azimuthal coordinate phi measures magnetic local time. The best functional representation found for the normalized difference (L(exp 3)a(exp 3)/mu(sub E))(B(sub m) - B(sub 0)) between mirror-point field B(sub m) and equatorial field B(sub 0) along any field line is a 5-term expansion in powers (2/3 through 6/3) of the quantity X equivalent to (La/mu(sub E))(exp 1/2)K, where K equivalent to (J(exp 2)/8m(sub 0)M)(exp 1/2) is an adiabatically conserved quantity independent of particle energy, m(sub 0) is the rest mass of the particle, and a is the radius of the Earth. This functional form is motivated by results for limiting cases in which particles mirror very near and very far from the magnetic equator. Expansion coefficients corresponding to various powers of X are obtained from least squares fits to numerically computed results for X as a function of L and B(sub m). These are accurately expressible as fourth-order polynomials in (r(sub 0)/b)(exp 3), hence indirectly as functions of L/L* = 3La/2b. This representation, which leads (except for a manageably small region of parameter space) to better than 1% accuracy in the specification of B(sub m) as a function of K and L, allows bounce-averaged guiding center simulations to be performed without actually tracing the bounce motions of individual particles. Bounce-averaged drifts L' (meridional) and phi' (azimuthal) are proportional to derivatives of the Hamiltonian H (sum of kinetic and potential energies) with respect to phi and L, respectively. Our formulation thus provides a computationally efficient method for tracing the bounce-averaged adiabatic motion (conserving all three invariants) and nonadiabatic transport (violating the third invariant while conserving the first two invariants) of geomagnetically trapped particles in the model magnetosphere.

Schulz, Michael↗