Search NASA⌕ Search

SEARCH · Search NASA

Results for “trajectory simulation”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Calculation of ion–ion mutual neutralization rate constants using Landau–Zener theory coupled with trajectory simulations for Ar + –Cl − , Br − , I −

In this computational study, we self-consistently calculate the rate constants of mutual neutralization reactions by incorporating the electron transfer probability, using Landau–Zener state transition theory with inputs derived from ab initio quantum chemistry calculations, into classical trajectory simulations. Electronic structure calculations are done using correlation consistent basis sets with multi-reference configuration interaction to map all the molecular electronic states below the ion-dissociation limit as a function of the distance between the reacting species. Our electronic structure calculations have been significantly improved from our previous work through improved selection of molecular electronic configurations maintaining a fine grid of 1a 0 over a wide range of bond lengths and accurate treatment of spin–orbit couplings. Non-adiabatic coupling matrix elements are calculated with the three-point central difference method near each avoided crossing to estimate the exact crossing point R x and coupling parameter H if , which are inputs to the multi-channel Landau–Zener theory to calculate the electron transition probability. Our approach is applied to estimate the mutual neutralization rate constants for the following ion pairs: Ar + –Cl − , Ar + –Br − , Ar + –I − at ∼133 Pa. Furthermore, our predictions are compared against the experimental data reported. It is seen that the improvement in the electronic structure calculation results in excellent agreement between the simulation results and the available experimental data to within a factor of ∼2 or ∼±50%.

Complete-active space self-consistent field↗
DOE OSTI2025

Langevin Dynamics modeling of gas-phase ion-ion recombination (Final Technical Report)

A self-consistent trajectory simulation approach to model MN reactions (Fig. 1) which incorporates the probability of electron transfer as a Monte Carlo operator (Fig. 2) was developed and published as Liu et al. J. Chem. Phys. 159, 114111 (2023). The electron transfer probability p ET estimated using the two-state Landau-Zener (LZ) theory was incorporated into classical trajectory simulations to elicit predictions of MN reaction cross-section σ (vacuum) or rate constant β (finite pressure). Electronic structure calculations with multireference configuration interaction (MRCI) and large correlation consistent basis sets were used to derive inputs to the LZ theory. The key advance of our trajectory simulation approach is the incorporation of electron transfer probability and the inclusion of the effect of ion-neutral interactions on MN using a Langevin representation of the effect of neutral gas on ions. For H + – H - and Li + – H(D) - pairs, our approach quantitatively agrees with measured speed-dependent cross-sections for up to ~10 5 m/s. For the ion pair Ne + – Cl - , our predictions of the MN rate constant at ~1 torr are a factor of ~2 – 3 higher than the experimentally measured value. Similarly, for Xe + – F - in the pressure range of ~20000 – 80000 Pa, our predictions of the MN rate constant are ~20% lower but are in excellent qualitative agreement with experimental data. The paradigm of using trajectory simulations to self-consistently model MN reactions is the basis for inclusion of additional non-classical, and static magnetic and electric field effects. Subsequent work, published as Roy et al. focused on modeling recombination rate constant for three ion pairs (rare gas Ar + cation and halide anions): Ar + – Cl - , Ar + – Br - , Ar + – I - , 2) considering spin-orbit couplings in the electronic structure calculations to obtain high-fidelity estimates of the electron transfer probability and incorporated within the classical trajectory simulations to elicit predictions. In addition to calculations of ion-ion recombination rate constants, a classical trajectory simulation technique (published as Roy et al. J. Chem. Phys. 162(9), 094104 (2023)) that uses quaternions to represent orientation of non-spherical particles (ions or aerosol particles) was developed to simulate the recombination of diatomic or more generally, polyatomic molecules. Finally, several other ion pairs such as Ne + – Cl - , Kr + – Cl - , were explored using the developed semi-classical trajectory simulations to understand various challenges in tackling electronic structure calculations. Using empirical approaches to parameterize the electron transfer radius, trajectory simulations were also used to probe the effect of ion number density on MN rate constant.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Quasiclassical sampling and Wigner sampling of initial vibrational coordinates and momenta for polyatomic molecules in Monte Carlo molecular dynamics simulations

In a quasiclassical trajectory simulation, the vibrational modes are initialised with quantised vibrational energies, but vibrational phases are sampled by Monte Carlo. This requires an algorithm to assign coordinates and momenta to the various atoms. In this work, we present two methods for implementing this for nonrotating polyatomic molecules, namely, fixed-energy vibrational-state-selected initial conditions and thermal initial conditions. We also present a method for initiating classical trajectories with a ground-state Wigner distribution. These vibrational treatments are sufficient to initialise trajectories for unimolecular processes, and we also show how they can be applied to simulate bimolecular collision processes. The treatments of unimolecular and bimolecular collision processes are available in two Python codes called wigner_state_selected.py and bimolecular_collision.py, respectively, which will generate initial condition files that are recognisable by the SHARC and SHARC-MN computer programs for dynamics calculations. Both codes are available as standalone programs, as well as being included in SHARC-MN, and they will be included in future versions of SHARC. Here, the methods implemented in these codes are mostly also available in the ANT computer program, and those that are not available in ANT will be incorporated in future versions of ANT.

Wigner distribution↗
DOE OSTI2024

Data Science Shows that Entropy Correlates with Accelerated Zeolite Crystallization in Monte Carlo Simulations

We have performed a data science study of Monte Carlo simulation trajectories to understand factors that can accelerate formation of zeolite nanoporous crystals, a process that can take days or even weeks. In previous work, Monte Carlo simulations predicted and experiments confirmed that using a secondary organic structure-directing agent (OSDA) accelerates crystallization of all-silica LTA zeolite, with experiments finding a three-fold speedup [PCCP 24, 142-148 (2022)]. However, it remains unclear what physical factors cause the speed-up. Here, we apply data science to analyze the simulation trajectories to discover what drives accelerated zeolite crystallization in Monte Carlo going from a one-OSDA synthesis (1OSDA) to a two-OSDA version (2OSDA). We encoded simulation snapshots using the Smooth Overlap of Atomic Positions approach, which represents all 2- and 3-body correlations within a given cutoff distance. Principal component analyses failed to discriminate datasets of structures from 1OSDA and 2OSDA simulations, while the Support Vector Machine (SVM) approach succeeded at classifying such structures with an area-under-curve (AUC) score of 0.99 (where AUC = 1 is a perfect classification) with all 3-body correlations, and as high as 0.94 with only 2-body correlations. SVM decision functions reveal relatively broad / narrow histograms for 1OSDA / 2OSDA datasets, suggesting that the two simulations differ strongly in information heterogeneity. Informed by these results, we performed pair (2-body) entropy calculations during crystallization, resulting in entropy differences that semi-quantitatively account for the speedup observed in the previous Monte Carlo simulations. We conclude that altering synthesis conditions in ways that substantially changes the entropy of labile silica networks may accelerate zeolite crystallization, and we discuss possible approaches for achieving such acceleration.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2024

On the instrument-dependent appearance of ion dissociation events in atom probe tomography mass spectra

The successful application of atom probe tomography (APT) relies on the accurate interpretation of the mass spectrum (i.e. m/z histogram) from a sample. Some materials yield mass spectra that are amenable to a straightforward peak assignment/ranging, however, there are many materials that produce mass spectra with features that defy simple interpretation. One such example is Ga 2 O 3 which yields mass spectra containing several broad and difficult to interpret features. Herein, we study the GaO 2+ → O 1+ + Ga 1+ dissociation and we explain how this dissociation process gives rise to broad and previously unassigned features in the mass spectrum. Trajectory simulations are performed for the dissociation reaction utilizing realistic electrostatic models and compared to experiments using commercially available straight flight and reflectron based local electrode (LE) APT instruments. It is shown that the appearance of these features is strongly dependent on the specific design of the time-of-flight (ToF) mass analyzer. Additionally, we explore how various experimental parameters can affect the appearance of the dissociation process in the one-dimensional (1D) mass spectrum and in the two-dimensional (2D) correlation histogram. While the focus of this work is on a particular dissociation process related to Ga 2 O 3 , the understanding gained in the course of these simulations and experiments should be applicable to the interpretation of dissociation processes in other materials.

47 OTHER INSTRUMENTATION↗
DOE OSTI2025

Active Learning Surrogates for Integrating Electron Microscopy and Computational Insights from Simulations in Autonomous Experiments

Artificial Intelligence (AI) combined with simulations and experiments has great potential to accelerate scientific discovery across technology and pharmaceuticals. However, the gap between simulations and experiments is challenging due to disparities in time and scale, making it difficult to estimate properties like energy and electronic states from experiments, and to provide feedback based on theoretical insights.Our research addresses the challenge by developing unique deep kernel based surrogate models that learns from microscopic images, mapping structural features to energy differences from defect formation. We start with full-training using simulated images to determine optimal settings, establishing a baseline for active learning. Using these settings from the baseline, active learning is trained, and predicts structures along simulation trajectories based on uncertainty and energetic stability, thus reducing data requirements, simulation time and computational costs. The results demonstrate that the model achieves a low average error margin of approximately 0.03 meV, indicating good performance. To enhance feature extraction and reconstruction capabilities, we developed an autoencoder-decoder as additional surrogate to create latent space to capture essential features, enabling precise comparisons between simulations and experiments. The results from this model achieved a reconstruction loss of around 0.2 and accurately reconstructed molecular structures.Overall, this work advances the steering of experiments through computational simulations by employing a surrogate models that actively predicts the trajectories of structural evolution, achieving time-to-solution comparable to experimental measurements.

Saranathan, Gayathri [Hewlett-Packard]↗
DOE OSTI2025

Selective capture and recovery of uranium oxide colloids from aqueous soil suspensions using high gradient magnetic filtration

High Gradient Magnetic Filtration (HGMF) is a promising method for the selective capture and recovery of uranium oxide from surface soils. To date, however, magnetic filtration of uranium oxide has only been demonstrated at a proof-of-principle scale using relatively small filters (<5 cm 3 ) at low flowrates (<60 mL/min). Here, to explore the efficacy of magnetic filtration of uranium oxide at a larger scale, a newly designed HGMF apparatus that is more than an order of magnitude larger than our earlier filters (106 cm 3 ) was designed, fabricated, and tested at relatively high flowrates. Filtration experiments were performed using aqueous uranium oxide particle suspensions with Arizona Road Dust (ARD) as a soil simulant. At a flowrate of 125 mL/min, the apparatus’ uranium capture rate was exceptionally high (96 %), but selectivity was poor due to the high rate of capture for diamagnetic soil constituents (e.g., 77 % for silicon). All particles were captured at a lower rate when the flowrate was increased to 250 mL/min, but uranium selectivity was significantly increased due to the more substantial reduction in diamagnetic particle capture (i.e., capture rate of 77 % and 15 % for uranium and silicon, respectively). When backwashing the apparatus at the same flowrates used during filtration experiments, the rate of uranium recovery tended to be fairly low. Nevertheless, higher flowrates (1 L/min) and sonication were both shown to be highly effective methods of increasing uranium recovery. Magnetic field simulations were also performed to investigate potential optimizations to the design of the apparatus. These simulations showed that the intensity of the applied magnetic field could be increased by increasing the thickness of the steel magnetic housing. Additionally, stochastic trajectory simulations were performed to investigate the potential mechanisms of particle capture.

HGMF↗
DOE OSTI2025

SynopFrame: Multiscale time-dependent visual abstraction framework for analyzing DNA nanotechnology simulations

We present an open-source framework, SynopFrame, that allows DNA nanotechnology (DNA-nano) experts to analyze and understand molecular dynamics simulation trajectories of their designs. We use a multiscale multi-dimensional abstraction space, connect the representations to a projected conformational space plot of the structure’s temporal sequence, and thus enable experts to analyze the dynamics of their structural designs and, specifically, failure cases of the assembly. In addition, our time-dependent abstraction representation allows the biologists, for the first time in a smooth and structurally clear way, to identify and observe temporal transitions of a DNA-nano design from one configuration to another, and to highlight important periods of the simulation for further analysis. We realize SynopFrame as a dashboard of the different synchronized 3D spatial and 2D schematic visual representations, with a color overlay to show essential properties such as the status of hydrogen bonds. The linking of the spatial, schematic, and abstract views ensures that users can effectively analyze the high-frequency motion. We also categorize the status of the hydrogen bonds into a new format to allow us to color-encode it and overlay it on the representations. To demonstrate the utility of SynopFrame, we describe example usage scenarios and report user feedback.

Abstraction space↗
DOE OSTI2025

Impact of Newly Measured 𝛽-Delayed Neutron Emitters around 78 Ni on Light Element Nucleosynthesis in the Neutrino Wind Following a Neutron Star Merger

Neutron emission probabilities and half-lives of 37 𝛽-delayed neutron emitters from 75 Ni to 92 Br were measured at the RIKEN Nishina Center in Japan, including 11 one-neutron and 13 two-neutron emission probabilities and six half-lives for the first time that supersede theoretical estimates. These nuclei lie in the path of the weak 𝑟 process occurring in neutrino-driven winds from the accretion disk formed after the merger of two neutron stars synthesizing elements in the 𝐴∼80 abundance peak. The presence of such elements dominates the accompanying kilonova emission over the first few days and have been identified in the AT2017gfo event, associated to the gravitational wave detection GW170817. Abundance calculations based on over 17,000 simulated trajectories describing the evolution of matter properties in the merger outflows show that the new data lead to an increase of 50%–70% in the abundance of Y, Zr, Nb, and Mo. This enhancement is large compared to the scatter of relative abundances observed in old very metal poor stars and thus is significant in the comparison with other possible astrophysical processes contributing to the light-element production. These results underline the importance of including experimental decay data for very neutron-rich 𝛽 -delayed neutron emitters into 𝑟 -process models.

59 ≤ A ≤ 89↗
DOE OSTI2025

Data-driven equation-free dynamics applied to many-protein complexes: The microtubule tip relaxation

Microtubules (MTs) constitute the largest components of the eukaryotic cytoskeleton and play crucial roles in various cellular processes, including mitosis and intracellular transport. The property allowing MTs to cater to such diverse roles is attributed to dynamic instability, which is coupled to the hydrolysis of GTP (guanosine-5'-triphosphate) to GDP (guanosine-5'-diphosphate) within the β-tubulin monomers. Understanding the equilibrium dynamics and the structural features of both GDP- and GTP-complexed MT tips, especially at an all-atom level, remains challenging for both experimental and computational methods because of their dynamic nature and the prohibitive computational demands of simulating large, many-protein systems. This study employs the “equation-free” multiscale computational method to accelerate the relaxation of all-atom simulations of MT tips toward their putative equilibrium conformation. Using large MT lattice systems (14 protofilaments × 8 heterodimers) comprising ~21-38 million atoms, we applied this multiscale approach to leapfrog through time and nearly double the computational efficiency in realizing relaxed all-atom conformations of GDP- and GTP-complexed MT tips. Commencing from an initial 4 μs unbiased all-atom simulation, we interleave coarse projective “equation-free” jumps with short bursts of all-atom molecular dynamics simulation to realize an additional effective simulation time of 1.875 μs. Our 5.875 μs of effective simulation trajectories for each system expose the subtle yet essential differences in the structures of MT tips as a function of whether β-tubulin monomer is complexed with GDP or GTP, as well as the lateral interactions within the MT tip, offering a refined understanding of features underlying MT dynamic instability. Furthermore, the approach presents a robust and generalizable framework for future explorations of large biomolecular systems at atomic resolution.

Wu, Jiangbo [University of Chicago, IL (United Sta↗
DOE OSTI2025

Physicochemical and Molecular Insights into the Boundary Layer and Free Troposphere Aerosol Interactions over the Southern Great Plains

Ambient aerosols’ vertical profiles are critical for evaluating the role of aerosols in atmospheric chemistry and radiative transfer, but limited data on these profiles hinders our ability to fully assess their impact on the Earth's radiative balance. Here, in this study, we investigated the size-, time-, and altitude resolved composition of individual particles and bulk molecular composition of particle samples collected by an uncrewed aerial system–ArcticShark over the Southern Great Plains. Single particle microanalysis shows that, the free tropospheric (FT) samples are dominated (56-66%) by carbonaceous sulfate particles, while boundary layer (BL) samples are dominated (57-74%) by carbonaceous particles. Back trajectory simulations suggest that FT particles are likely influenced by long-range transport and have undergone aqueous-phase processing. Conversely, in-situ size distribution data shows evidence of particle growth in the upper BL and just below the FT. This observation may indicate vertical transport of particles from an elevated aerosol layer in the FT, possibly linked to a new particle formation event. This observation is further supported by high resolution molecular composition data, which reveals particle volatility increasing with increasing size, which aligns with the growth event. This study aids in fundamental understanding of the compositional and molecular specificity of vertically resolved organic aerosols to provide insights into particle size evolution for future atmospheric models.

ArcticShark↗
DOE OSTI2025

Identifying High Ionic Conductivity Compositions of Ionic Liquid Electrolytes Using Features of the Solvation Environment

Binary mixtures of ionic liquids with molecular solvents are gaining interest in electrochemical applications due to the improvement in their performance over neat ionic liquids. Dilution with suitable molecular solvents can reduce the viscosity and facilitate faster diffusion of ions, thereby yielding substantially higher ionic conductivity than that for a pure ionic liquid. Although viscosity and diffusion coefficients typically behave as monotonic functions of concentration, ionic conductivity often passes through a peak value at an optimum molar ratio of the molecular solvent to the ionic liquid. The ionic conductivity maximum is generally explained in terms of a balance between the ease of charge transport and the concentration of the charge carriers. In this work, fluctuation in the local environment surrounding an ion is invoked as a plausible explanation for the ionic conductivity mechanism with a binary mixture of 1-ethyl-3-methylimidazolium tetrafluoroborate and ethylene glycol as an example. The magnitude of the dynamism in the local environment is captured by measuring the spatial and temporal features of the solvation environment. Standard deviation in the number of ions in the solvation environment serves as a spatial feature, while the cage correlation lifetimes for oppositely charged ions within the first solvation shell serve as a temporal feature. Large standard deviations in the cluster ion population and short cage correlation lifetimes are indicators of highly dynamic ionic environment at the molecular level and consequently yield high ionic conductivity. Such compositions were found to be in good agreement with the optimum ionic liquid mole fractions obtained through experimental measurement. Short cage correlation lifetimes enable the identification of optimum mixture compositions using simulation trajectories significantly shorter than those required to implement the Nernst–Einstein or Einstein formalisms for calculating ionic conductivity. We validated the applicability of this approach across force fields and in six ionic liquid-molecular solvent electrolytes formed with combination of cations, anions, and solvents. We offer a computationally efficient approach of screening ionic liquid-molecular solvent binary mixture electrolytes to identify molar ratios that yield high ionic conductivity.

25 ENERGY STORAGE↗
DOE OSTI2025

The UV Photoinduced Ring-Closing Reaction of Cyclopentadiene Probed with Ultrafast Electron Diffraction

Conjugated cyclic organic molecules are common across many fields such as pharmaceuticals, are naturally occurring in biological systems, and are used in synthetic materials. One particular area of interest from a photochemical point of view is the formation of highly strained cyclic organics. We investigate the photoinduced reaction of cyclopentadiene, a five-membered organic ring molecule which can form strained three and four carbon rings after photoexcitation with UV light, with the gas-phase ultrafast electron diffraction instrument at the SLAC MeV-UED facility. Electron diffraction offers a direct probe sensitive to the nuclear geometry during the reaction, allowing for the determination of the distribution of products formed following photoexcitation. We observe the simultaneous formation of the highly strained ring- closed bicyclo[2.1.0]pentene and vibrationally hot cyclopentadiene within the temporal resolution of the experiment, and determine the relative yield of all reaction products. Furthermore, the experimental results are in good agreement with the predictions of trajectory simulations.

Computational Chemistry↗
DOE OSTI2026

Ultrafast Proton-Coupled Electron Transfer Reactions at Unit Quantum Yield

Proton-coupled electron transfer (PCET) reactions mediate essential biological energy-conversion processes, yet dynamical experiments that distinguish concerted electron/proton motion from sequential electron-transfer/proton-transfer pathways are unprecedented on ultrafast time scales. Here we report a tightly coupled donor–spacer–acceptor assembly, NDI-Sp-PhOH/BI, that enables direct interrogation of light-triggered and ultrafast PCET dynamics. Photoexcitation of the naphthalene diimide chromophore drives charge separation forming NDI radical anion and benzimidazolium phenoxy radical (PhO●/BIH+) products with a 440 fs time constant, followed by thermal charge recombination on a 5 ps time scale. Transient visible/near-infrared and mid-infrared spectroscopies independently track electron and proton motion and reveal coincident formation and decay of the charge-separated PCET products, excluding detectable phenol radical cation or phenoxide intermediates expected for sequential electron transfer-proton transfer (ET-PT) pathways. Surface-hopping quasiclassical trajectory simulations demonstrate rapid proton relocation coupled to evolution from locally excited to charge-transfer states, congruent with a concerted mechanism. These combined spectroscopic and computational results establish NDI-Sp-PhOH/BI as an archetype for ultrafast, concerted PCET occurring with unit quantum yield; further, because this platform provides distinct electronic and vibrational spectroscopic handles in its ground, electronically excited, and charge-separated states, it provides new opportunities to characterize mechanistic crossover in biomimetic PCET systems.

10 SYNTHETIC FUELS↗
DOE OSTI2024

Ion Mobility Separations Using Cocentric Architecture

Ion mobility separations are usually performed in linear channels, which, when extended, can have a large footprint. In this work, we explored the performance of an ion mobility device with a curved architecture which can have a more compact form. The Co-centric Ion Mobility Spectrometer (CIMS) works by manipulating ions between two co-centric surfaces, each containing a serpentine track. The mobility separation inside CIMS is achieved using traveling waveforms (TWs). We initially evaluated the device using ion trajectory simulations using SIMION, which indicated that when ions traveled circularly inside CIMS, they resulted in similar resolving powers and transmitted m/z range as traveling in a straight path in structures for lossless ion manipulations (SLIM). We then performed experimental validation of CIMS in conjunction with a TOF MS. The CIMS was made of 2 flexible printed circuit board materials folded into concentric cylinders separated by a gap of 2.8 mm. The device was about 50 mm diameter × 152 mm long and provided 1.846 m of serpentine path length. Three sets of mixtures (Agilent tune mixture, tetraalkylammonium salts, and 8 peptide mixture) and four traveling waveform profiles (square, sine, triangle, and sawtooth) were used. The sawtooth TW profile produced a slightly higher resolving power for the Agilent tuning mixture and tetraalkylammonium ions. The average resolving power for Agilent tune mixture ions ranged from 37 (using sawtooth TW) to 27 (using square TW). For tetraalkylammonium ions, the average resolving powers ranged from 45 (sawtooth TW) to 31 (square TW). For the peptide mixture ions, the resolving power was similar among the four TW profiles and ranged from 51 to 56. The average percent error in TW CCS for the peptide mixture ions ranged was about 0.4%. In conclusion, the new device showed promising results for a device made of a flexible printed circuit board material, but improvements are needed to further increase the resolving power.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Relationship between microporous structure and light gas transport through glassy polymeric membranes revealed by molecular simulations

Microporous glassy polymers are attractive materials for gas separation membranes, due to their high permeability and tailorable selectivity, provided physical aging can be delayed. The archetypal microporous glassy polymer, PTMSP, can be blended with a hyper-crosslinked isatin–triptycene porous polymer network (PPN) to delay physical aging. However, while PPN is effective at reducing physical aging, it also affects the permeability of light gases through PTMSP. Molecular dynamics simulations were used here to shed fundamental light on the mechanisms responsible for these effects. Atomistic models are developed that satisfactorily reproduce experimental observations such as matrix density and cavity size distributions for neat PTMSP as well as for PTMSP–PPN blends. Analysis of the simulation results suggests that physical aging is delayed because the PPN inclusions slow down PTMSP relaxation while reducing the connectivity between free volume pockets. To understand how PPN inclusions affect light gas permeability, the atomistic models developed are used to probe CO 2 and CH 4 diffusion and sorption. These simulations are conducted for PTMSP matrices exhibiting varying density, towards reproducing experimental permeability data. Interrogating the simulation trajectories, it is discovered that while CH 4 travels preferentially through the free volume cavities, CO 2 preferentially interacts with the available surfaces, especially in the presence of PPN. These differences help interpret experimental observations. The transport of C 2 H 6 and H 2 S through PTMSP matrices was also investigated. These gases, and in particular H 2 S, were found to absorb within PTMSP, yielding very low diffusion coefficients. The ability to predict differences in diffusion pathways for various natural gas components unveils the possibility of engineering membranes to control permeability and selectivity towards large-scale applications.

Bao Le, Tran Thi [Univ. of Oklahoma, Norman, OK (U↗
DOE OSTI2024

Heavy ion beam probe for Wendelstein 7-X measurement capabilities as projected through its design

A heavy ion beam probe (HIBP) diagnostic is being developed for studies of plasma equilibrium and turbulence in the optimized Wendelstein 7-X (W7-X) stellarator. Operation of W7-X has experimentally demonstrated that its optimized magnetic field results in improved neoclassical particle confinement and, as a result, turbulence is the predominant cause of energy transport. The HIBP will have the unique ability to provide experimental data needed to complement models of both neoclassical and turbulent transport. It will acquire direct measurements in the W7-X plasma interior of the electric potential (needed for understanding ambipolar particle flux) and fluctuations of electron density and potential (needed for understanding turbulence). The HIBP for W7-X will inject singly charged ion beams with energies of up to 2 MeV and is designed to access the upper cross section of the W7-X plasma. We use trajectory simulations to illustrate the plasma coverage that the diagnostic can achieve in the reference magnetic configurations of W7-X. We calculate beam signal levels, discuss anticipated measurement sensitivity of broadband fluctuations of electron density and plasma potential, and show how they depend on plasma density. We also discuss the diagnostic sensitivity to equilibrium plasma potential.

Instruments & Instrumentation↗
DOE OSTI2025

Exploring the Utility-Privacy Trade-Off: Impacts of Semantic and Visit Types Ambiguities on Human Mobility Simulation

Humans are in perpetual movement, constantly traversing buildings, cities, waters, oceans, and countries. Mobility stands out as a major driving force shaping our modern societies. Capturing and explaining human behavior in a world of eight billion distinct mobility agendas is a complex challenge. With the rise of interconnected devices and platforms, such as smartphones, wearables, and point-of-interest data, largescale behavioral data has become more accessible, enabling rich insights into mobility patterns. However, the widespread availability of such data introduces significant ethical challenges. Detailed mobility data can inadvertently reveal sensitive personal information, including individuals' locations, habits, social interactions, and even political or religious affiliations. Beyond privacy breaches, the ethical implications of uncovering and potentially manipulating underlying behavioral patterns demand attention. Striking a balance between the utility of mobility models and the protection of individual privacy is therefore paramount. This paper explores the utility-privacy trade-offs in human mobility modeling, focusing on the impacts of introducing semantic and visit type ambiguities. By systematically examining how these ambiguities affect the fidelity of simulated trajectories and privacy risks, we provide a framework for evaluating ethical and privacy-conscious modeling practices. Our findings emphasize the need for methods that safeguard privacy without undermining the usefulness of mobility models, contributing to the responsible advancement of mobility science in alignment with ethical standards and societal expectations.

Amichi, Licia [ORNL] (ORCID:0000000177631394)↗