Search NASA⌕ Search

SEARCH · Search NASA

Results for “strong field processes”

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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Magnetism in curved geometries

Curvature impacts physical properties across multiple length scales, ranging from the macroscopic scale, where the shape and size vary drastically with the curvature, to the nanoscale at interfaces and inhomogeneities in materials with structural, chemical, electronic, and magnetic short-range order. In quantum materials, where correlations, entanglement, and topology dominate, the curvature opens the path to novel characteristics and phenomena that have recently emerged and could have a dramatic impact on future fundamental and applied studies of materials. Particularly, magnetic systems hosting non-collinear and topological states and 3D magnetic nanostructures strongly benefit from treating curvature as a new design parameter to explore prospective applications in the magnetic field and stress sensing, microrobotics, and information processing and storage. This Perspective gives an overview of recent progress in synthesis, theory, and characterization studies and discusses future directions, challenges, and application potential of the harnessing curvature for 3D nanomagnetism.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Magnetic field induced antiferromagnetic cone structure in multiferroic BiFeO 3

Neutron diffraction measurements were performed under high magnetic fields up to 17 T in a multiferroic BiFeO 3 single crystal, in which an intermediate magnetic (IM) phase has been found between the cycloid and canted antiferromagnetic phases [S. Kawachi et al., Phys. Rev. Mater. 1, 024408 (2017)]. We clearly found that the incommensurate magnetic peaks, which split perpendicular to the magnetic field in the cycloid phase, rotate by 90 deg to align parallel to the field in the IM phase. The magnetic structure in the IM phase can be best described by an antiferromagnetic cone (AF cone) structure. The transition from the cycloid to AF cone is of first order and the direction of the magnetic wave vector and the easy plane of the cycloidal component are rotated by 90 deg without changing the cycloidal modulation period, whereas the transition from the AF cone to canted antiferromagnetic phase is gradual and the cone angle becomes smaller gradually without changing the modulation period. Interestingly, the cycloidal component as well as the cone angle in the IM phase shows a large hysteresis between the field increasing and decreasing processes. Finally, this result, combined with the magnetostriction with a large hysteresis previously reported in the IM phase, suggests a strong magnetoelastic coupling.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Anisotropic Heating and Parallel Heat Flux in Electron-only Magnetic Reconnection with Intense Guide Fields

Electron-only reconnection (E-REC) is a process recently observed in the Earth’s magnetosheath, where magnetic reconnection occurs at electron kinetic scales, and ions do not couple to the reconnection process. Electron-only reconnection is likely to have a significant impact on the energy conversion and dissipation of turbulence cascades at kinetic scales in some settings. This paper investigates E-REC under different intensities of strong guide fields (the ratio between the guide field and the in-plane asymptotic field strength is 5, 10 and 20, respectively) via two-dimensional fully kinetic particle-in-cell simulations, focusing on electron heating. The simulations are initialized with a force-free current sheet equilibrium under various intensities of strong guide fields. Similarly to previous experimental studies, electron temperature anisotropy along separatrices is observed, which is found to be mainly caused by the variations of parallel temperature. Both regions of anisotropy and parallel temperature increase/decrease along separatrices become thinner with increasing guide fields. Besides, we find a transition from a quadrupolar to a hexapolar (six-polar) to an octopolar (eight-polar) structure in temperature anisotropy and parallel temperature as the guide field intensifies. Non-Maxwellian electron velocity distribution functions (EVDFs) at different locations in the three simulations are observed. Our results show that parallel electron velocity varies notably with different guide field intensities and finite parallel electron heat flux density is observed. The three simulations exhibit features of the Chew–Goldberger–Low theory, with the level of consistency increasing as the guide field strength increases. This explains the electron parallel temperature variations and the shape of the EVDFs observed along the separatrices. This work may provide insights into the understanding of electron heating and parallel heat flux density in E-REC observed in the turbulent magnetosheath.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2015

Effects of Atomic-Scale Structure on the Fracture Properties of Amorphous Carbon - Carbon Nanotube Composites

The fracture of carbon materials is a complex process, the understanding of which is critical to the development of next generation high performance materials. While quantum mechanical (QM) calculations are the most accurate way to model fracture, the fracture behavior of many carbon-based composite engineering materials, such as carbon nanotube (CNT) composites, is a multi-scale process that occurs on time and length scales beyond the practical limitations of QM methods. The Reax Force Field (ReaxFF) is capable of predicting mechanical properties involving strong deformation, bond breaking and bond formation in the classical molecular dynamics framework. This has been achieved by adding to the potential energy function a bond-order term that varies continuously with distance. The use of an empirical bond order potential, such as ReaxFF, enables the simulation of failure in molecular systems that are several orders of magnitude larger than would be possible in QM techniques. In this work, the fracture behavior of an amorphous carbon (AC) matrix reinforced with CNTs was modeled using molecular dynamics with the ReaxFF reactive forcefield. Care was taken to select the appropriate simulation parameters, which can be different from those required when using traditional fixed-bond force fields. The effect of CNT arrangement was investigated with three systems: a single-wall nanotube (SWNT) array, a multi-wall nanotube (MWNT) array, and a SWNT bundle system. For each arrangement, covalent bonds are added between the CNTs and AC, with crosslink fractions ranging from 0-25% of the interfacial CNT atoms. The SWNT and MWNT array systems represent ideal cases with evenly spaced CNTs; the SWNT bundle system represents a more realistic case because, in practice, van der Waals interactions lead to the agglomeration of CNTs into bundles. The simulation results will serve as guidance in setting experimental processing conditions to optimize the mechanical properties of CNT composites.

Jensen, Benjamin D.↗
DOE OSTI2021

Strong Purcell enhancement at telecom wavelengths afforded by spinel Fe 3 O 4 nanocrystals with size-tunable plasmonic properties

Developments in the field of nanoplasmonics have the potential to advance applications from information processing and telecommunications to light-based sensing. Traditionally, nanoscale noble metals such as gold and silver have been used to achieve the targeted enhancements in light-matter interactions that result from the presence of localized surface plasmons (LSPs). However, interest has recently shifted to intrinsically doped semiconductor nanocrystals (NCs) for their ability to display LSP resonances (LSPRs) over a much broader spectral range, including the infrared (IR). Among semiconducting plasmonic NCs, spinel metal oxides (sp-MOs) are an emerging class of materials with distinct advantages in accessing the telecommunications bands in the IR and affording useful environmental stability. In this work, we report the plasmonic properties of Fe 3 O 4 sp-MO NCs, known previously only for their magnetic functionality, and demonstrate their ability to modify the light-emission properties of telecom-emitting quantum dots (QDs). We establish the synthetic conditions for tuning sp-MO NC size, composition and doping characteristics, resulting in unprecedented tunability of electronic behavior and plasmonic response over 450 nm. In particular, with diameter-dependent variations in free-electron concentration across the Fe 3 O 4 NC series, we introduce a strong NC size dependency onto the optical response. In addition, our observation of plasmonics-enhanced decay rates from telecom-emitting QDs reveals Purcell enhancement factors for simple plasmonic-spacer-emitter sandwich structures up to 51-fold, which are comparable to values achieved previously only for emitters in the visible range coupled with conventional noble metal NCs.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Visualizing the strong field–induced molecular breakup of C 60 via x-ray diffraction

Laser-driven dynamics in polyatomic molecules poses a complex many-body problem. Understanding intense light-matter interaction is crucial for steering intramolecular quantum dynamical processes. Here, we record time-resolved x-ray diffraction images of C 60 molecules during and after their interaction with intense near-infrared fields, giving direct access to structural changes of the molecules and their fragmentation in real time. Tuning the intensity of the excitation pulses, we uncover a transition from a weak-field regime of excited but stable molecules to a high-field regime dominated by Coulomb explosion. In the transition region, the molecules expand by up to 50% of their initial size within just 140 fs, with major fragmentation only setting in afterward. This work demonstrates that x-ray diffractive imaging is capable of retrieving time-resolved structural information of large molecules reshaped by intense laser fields. Laser-driven fragmentation is a first step toward observing molecular processes modified by laser fields of increasing intensity.

Schnorr, Kirsten [Paul Scherrer Inst. (PSI), Villi↗
DOE OSTI2025

Magnetically Driven Neutron-rich Ejecta Unleashed: Global 3D Neutrino–General Relativistic Magnetohydrodynamic Simulations of Collapsars Probe the Conditions for r -process Nucleosynthesis

Abstract Collapsars—rapidly rotating stellar cores that form black holes—can power gamma-ray bursts and are proposed to be key contributors to the production of heavy elements in the Universe via the rapid neutron capture process ( r -process). Previous neutrino-transport collapsar simulations have been unable to unbind neutron-rich material from the disk. However, these simulations have not included sufficiently strong magnetic fields and the black hole (BH), both of which are essential for launching mass outflows. We present ν h-amr , a novel neutrino-transport general relativistic magnetohydrodynamic ( ν GRMHD) code, which we use to perform the first 3D global ν GRMHD collapsar simulations. We find a self-consistent formation of a weakly magnetized dense accretion disk, which has sufficient time to neutronize. Eventually, substantial magnetic flux accumulates near the BH, becomes dynamically important, leads to a magnetically arrested disk (MAD), and unbinds some of the neutron-rich material. However, the strong flux also hinders accretion, lowers density, and increases neutrino-cooling timescale, which prevents further disk neutronization. Typical collapsar progenitors with mass accretion rates, M ̇ ∼ 0.1 − 1 M ⊙ s - 1 , do not produce significant neutron-rich ( Y e < 0.25) ejecta. However, we find that MADs at higher mass accretion rates, M ̇ ≳ few M ⊙ s - 1 (e.g., for more centrally concentrated progenitors), can unbind M ej ≲ M ⊙ of neutron-rich ejecta. The outflows inflate a shocked cocoon that mixes with the infalling neutron-poor stellar gas and raises the final outflow Y e ; however, the final r -process yield may be determined earlier at the point of neutron capture freeze-out. Future work will explore under what conditions more typical collapsar engines become r -process factories.

Issa, Danat (ORCID:0009000524787631)↗
DOE OSTI2022

Beyond-mean-field approaches for nuclear neutrinoless double beta decay in the standard mechanism

Nuclear weak decays provide important probes to fundamental symmetries in nature. A precise description of these processes in atomic nuclei requires comprehensive knowledge on both the strong and weak interactions in the nuclear medium and on the dynamics of quantum many-body systems. In particular, an observation of the hypothetical double beta decay without emission of neutrinos (0νββ) would unambiguously demonstrate the Majorana nature of neutrinos and the existence of the lepton-number-violation process. It would also provide unique information on the ordering and absolute scale of neutrino masses. The next-generation tonne-scale experiments with sensitivity up to 10 28 years after a few years of running will probably provide a definite answer to these fundamental questions based on our current knowledge on the nuclear matrix element (NME), the precise determination of which is a challenge to nuclear theory. Beyond-mean-field approaches have been frequently adapted for the study of nuclear structure and decay throughout the nuclear chart for several decades. In this review, we summarize the status of beyond-mean-field calculations of the NMEs of 0νββ decay assuming the standard mechanism of an exchange of light Majorana neutrinos. In conclusion, the challenges and prospects in the extension and application of beyond-mean-field approaches for 0νββ decay are discussed.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1999

Coherent Pump-Probe Interactions and Terahertz Intersubband Gain in Semiconductor Quantum Wells

In recent years there has been considerable interest in intersubband-transition-based infrared semiconductor quantum well (QW) lasers because of their potential applications. In the mid-infrared range, both electrically-injected quantum cascade lasers [1] and optically-pumped multiple QW lasers [2] have been experimentally realized. In these studies, optical gain is due to population inversion between the lasing subbands. It was also proposed that stimulated Raman scattering in QW systems can produce net infrared optical gain [3j. In such a nonlinear optical scheme, the appearance of optical gain that may lead to intersubband Raman lasers does not rely on the population inversion. Since, in tile resonant Raman process (Raman gain is the largest in this case), the pump field induces population redistribution among subbands in the QW s ystem, it seems that a realistic estimate of the optical gain has to include this effect. Perturbative calculations used in the previous work [3] may overestimate the Raman gain. In this paper we present a nonperturbative calculation of terahertz gain of optically-pumped semiconductor step quantum wells. Limiting optical transitions within the conduction band of QW, we solve the pump-field-induced nonequilibrium distribution function for each subband of the QW system from a set of coupled rate equations. Both intrasubband and intersubband relaxation processes in the quantum well system are included. Taking into account the coherent interactions between pump and THz (signal) waves, we we derive the susceptibility of the QW system for the THz field. For a GaAs/AlGaAs step QW, we calculate the Thz gain spectrum for different pump frequencies and intensities. Under moderately strong pumping (approximately 0.3 MW/sq cm), a significant THz gain (approximately 300/m) is predicted. It is also shown that the coherent wave interactions (resonant stimulated Raman processes) contribute significantly to the THz gain.

Liu, Ansheng↗
DOE OSTI2024

Understanding aerosol–cloud interactions using a single-column model for a cold-air outbreak case during the ACTIVATE campaign

Marine boundary layer clouds play a critical role in Earth's energy balance. Their microphysical and radiative properties are highly impacted by ambient aerosols and dynamic forcings. In this study, we evaluate the representation of these clouds and related aerosol–cloud interaction processes in the single-column version of the E3SM climate model (SCM) against field measurements collected during the NASA ACTIVATE campaign over the western North Atlantic, as well as intercompare results with high-resolution process level models. We show that E3SM SCM reproduces the macrophysical properties of post-frontal boundary layer clouds in a cold-air outbreak (CAO) case well. However, it generates fewer but larger cloud droplets compared to aircraft measurements. Further sensitivity tests show that the underestimation of both aerosol number concentration and vertical velocity variance contributes to this bias. Aerosol–cloud interactions are examined by perturbing prescribed aerosol properties in E3SM SCM with fixed dynamics. Higher aerosol number concentration or hygroscopicity leads to more numerous but smaller cloud droplets, resulting in a stronger cooling via shortwave cloud forcing. This apparent Twomey effect is consistent with prior climate model studies. The cloud liquid water path shows a weakly positive relation with cloud droplet number concentration due to precipitation suppression. This weak aerosol effect on cloud macrophysics may be attributed to the dominant impact of strong dynamical forcing associated with the CAO. Our findings indicate that the SCM framework is a key tool to bridge the gap between climate models, process level models, and field observations to facilitate process level understanding.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Dilepton Program with Time-of-Flight Detector at the STAR Experiment

Pairs of lepton and antilepton (dilepton) in a continuous mass range are one of the most experimentally challenging and golden probes of the quark-gluon plasma (QGP) produced in heavy ion collisions because they do not strongly interact with the hot and dense medium, and reflect the properties of the medium at the time the dilepton is generated. The measurements of dileptons require lepton identification with high purity and high efficiency at large detector acceptance. STAR is one of two large experiments at the relativistic heavy ion collider with a primary goal of searching for the QGP and studying its properties. The STAR experiment launched a comprehensive dielectron (e + e - ) program enabled by the time-of-flight (TOF) detector that had been fully installed in 2010. In this article, we review the decade-long R&D, the construction and performance of the STAR TOF detector, and dielectron measurements, including thermal dielectron production and dielectron production from the Breit–Wheeler process. Future perspectives are also discussed.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Strong-field Driven Sub-cycle Band Structure Modulation and Dephasing Control

In this work, we present measurements of ultrafast electric-field observables in magnesium oxide using a non-resonant nonlinear optical interaction. Using field observables, we show that strong laser fields modulate the band structure on sub-cycle timescales, thereby altering the material’s nonlinear optical response. We perform time-dependent perturbation theory calculations using a field-dependent dispersion relation and semiconductor Bloch equation calculations, both of which agree with experimental observations. Furthermore, we extract pulse decay times from the real-time signal electric field envelope and show sub-cycle control of dephasing times. Our work offers a new perspective on strong-field-driven electron dynamics in solids through electric-field observables. The demonstrated attosecond modulation of the nonlinear response could have important implications for quantum light generation and quantum spectroscopy using nonlinear optical processes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Growth and Properties of Vapor Diffused Nb3Sn Coating for Superconducting Radiofrequency Accelerator Cavity Applications

The desire to reduce construction and operating costs of future SRF accelerators motivates the search for higher-performing alternative materials. Nb3Sn (Tc ~ 18.3 K and Hsh ~ 425 mT) is the front runner. The tin vapor diffusion process is currently the technique of choice to produce promising Nb3Sn-coated cavities. Understanding Nb3Sn nucleation and growth in this process is essential to progress. Samples representing different stages of Nb3Sn formation have been produced and studied to elucidate the effects of nucleation, growth, process conditions, and impurities. Nb3Sn films with thickness from a few hundred nm up to ~15 µm were grown and characterized. The microscopic examination of samples suggests the mechanisms of thin film nucleation and growth. Broadly, nucleation deposits tin as a thin surface phase and, under some conditions, a few hundred nanometer sized particles as well. Conditions that impair nucleation promote the formation of defects, such as ?patches?, in subsequent coating growth. Analysis of coated samples is consistent with the model of Nb3Sn grown in which tin diffuses via grain boundaries to Nb3Sn-Nb interface, where the growth Nb3Sn into the niobium bulk takes place. Similar scaling laws are found for grain growth and layer thickness. Non-parabolic layer growth is consistent with significant grain growth, which reduces the number of Sn-transport channels. Examination of patchy region in Nb3Sn coating revealed it to be large single crystalline grains, pointing to impeded Nb3Sn layer growth due to low grain boundary density, resulting in a significantly thin coating in those areas. Examination of RF loss regions from a coated cavity, identified with a thermometry mapping system showed patchy regions and carbonus defects were associated with strong local field-dependent surface resistance. RF measurements of coated cavities were combined with material characterization of witness samples and coated-cavity cutouts to improve the coating process. Understanding obtained and applied to cavity coatings, resulted in single-cell Nb3Sn cavities with a quality factor of ~2 ×1010 up to 15 MV/m accelerating gradient at 4 K, without "Wuppertal" Q-slope. We have also produced Nb3Sn-coated CEBAF 5-cell cavities with accelerating gradients useful for accelerator cryomodules. This dissertation will discuss the genesis of the Nb3Sn coating in a typical tin vapor diffusion process, effects of different process parameters, and its consequences to the coating of single-cell and multi-cell SRF cavities.

Pudasaini, Uttar↗
DOE OSTI2020

Advanced Simulation and Experiments of Strongly Coupled Geomechanics and Flow for Gas Hydrate Deposits: Validation and Field Application. Final Report

The objectives of the proposed research are (1) to investigate geomechanical responses induced by depressurization experimentally and numerically; (2) to enhance the current numerical simulation technology in order to simulate complex physically coupled processes by depressurization and (3) to perform in-depth numerical analyses of two selected potential production test sites: one based on the deposits observed at the Ulleung basin UBGH2-6 site; and the other based on well-characterized accumulations from the westend Prudhoe Bay. To these ends, the recipient will have the following specific objectives: 1). Information obtained from multi-scale experiments previously conducted at the recipient’s research partner (the Korean Institute of Geoscience and Mineral Resources (KIGAM)) that were designed to represent the most promising known Ulleung Basin gas hydrate deposit as drilled at site UBGH2-6 will be evaluated (Task 2). These findings will be further tested by new experimental studies at Lawrence Berkeley National Laboratory (LBNL) and Texas A&M (TAMU) (Task 3) that are designed capture complex coupled physical processes between flow and geomechanics, such as sand production, capillarity, and formation of secondary hydrates. The findings of Tasks 2 and 3 will be used to further improve numerical codes. 2) Develop (in Tasks 4 through 6) an advanced coupled geomechanics and non-isothermal flow simulator (T+MAM) to account for large deformation and strong capillarity. This new code will be validated using data from the literature, from previous work by the project team, and with the results of the proposed experimental studies. The developed simulator will be applied to both Ulleung Basin and Prudhoe Bay sites, effectively addressing complex geomechanical and petrophysical changes induced by depressurization (e.g., frost-heave, strong capillarity, cryo-suction, induced fracturing, and dynamic permeability).

02 PETROLEUM↗
DOE OSTI2020

Origin of the multiphoton-regime harmonic-generation plateau structure

A physical interpretation is provided for the formation of the multiphoton-regime plateau feature in the spectra of nonlinear and correlated process of high-order harmonic generation (HHG) of two-active-electron atoms interacting with an intense linearly polarized laser field. While in the strong-field tunneling regime the plateau feature is well known to be due to rescattering effects of the freed electron, its counterpart in the multiphoton regime is due to atomic resonance effects involving both singly excited states and doubly excited states. Here, we propose a strategy to uncover which kind of these states are responsible for this phenomenon. By studying the sensitivity of HHG to the initial state in helium, in this work we show that this multiphoton-regime plateau structure mostly originates from singly excited states.

74 ATOMIC AND MOLECULAR PHYSICS↗
NASA NTRS1985

Study of non-thermal photon production under different scenarios in solar flares. 2: The Compton inverse and Bremsstrahlung models and fittings

Energy spectra of photons emitted from Bremsstrahlung (BR) of energetic electrons with matter, is obtained from the deconvolution of the electron energy spectra. It can be inferred that the scenario for the production of X-rays and gamma rays in solar flares may vary from event to event. However, it is possible in many cases to associated low energy events to impulsive acceleration, and the high energy phase of some events to stochastic acceleration. In both cases, flare particles seem to be strongly modulated by local energy losses. Electric field acceleration, associated to neutral current sheets is a suitable candidate for impulsive acceleration. Finally, that the predominant radiation process of this radiation is the inverse Compton effect due to the local flare photon field.

Perez-Peraza, J.↗
NASA NTRS1994

Cosmic bullets as particle accelerators and radio sources

We have simulated in two dimensions the dynamical evolution of dense gas clouds(`cosmic bullets') moving supersonically through a uniform low-density medium. The diffusive shock acceleration of relativistic protons (cosmic rays) and their dynamical feedback on the background flow are included by the two-fluid model for this process. The acceleration of relativistic electrons is approximated by a test-particle model, and a passive magnetic field is followed by a simple advection scheme. Strong bow shocks, with Mach numbers similar to that of a bullet's motion, are the most important particle accelerators in the flow, while tail shocks and shocks inside the bullets do not play generally significant roles in this regard. For our simulation parameters, approximately greater than 10% of the initial bullet kinetic energy is converted to a combination of internal energy of gas and cosmic-ray protons by the time the bullets begin to be disrupted. Characteristically, the cosmic rays gain several percent of the available kinetic energy. Bullet destruction on timescales only a little larger than the ram pressure bullet crushing time begins in response to Kelvin-Helmholtz and especially to Rayleigh-Taylor instabilities along the forward bullet surface. For dense bullets this happens before the bullet is stopped by ram pressure. According to our simple model for synchrotron emission from relativistic electrons accelerated and transported within the flows, that emission increases rapidly as the bullet begins to fragment, when it is strongly dominated by field enhancement in sheared flows. Synchrotron emission from the acceleration region within the bow shock is, by contrast, much weaker.

Jones, T. W.↗
NASA NTRS2012

Doubly-Shadowed Regions in Lunar Polar Craters: Hydrogen Accumulation in the Presence of Recursive Plasma Wakes

Permanently shadowed regions (PSRs) of the Moon have been identified as unique environments of extreme cold and comprise a natural cold trap for sequestering volatiles [Paige et al. 2010]. The diverse chemical composition of the LCROSS impact plume provided evidence for a volatile-rich and chemically-complex PSR environment [Cola prete et al. 2010, Schultz et al. 2010]. Additionally, the polar electrostatic environment is highly complex, with the possibility of strong, localized electric fields that divert solar wind ions directly into polar cold traps [Farrell et al. 2010, Zimmerman et al. 2011]. Thus, regional plasma physics processes couple directly with volatile sequestration. In the present work, kinetic simulations show that recursive plasma wake structure arises in the presence of step-like topographic features (Le. doubly-shadowed craters). Combining the plasma code with a numerical sputtering model demonstrates that solar wind protons can be either a hydrogen source via implantation or a volatile loss mechanism via sputtering, depending on properties of the regolith and solar wind. The present model provides a novel theoretical pathway toward understanding the lunar surface/solar wind physical and chemical interactions for complex topography near the poles.

Zimmerman, Michael I.↗