SEARCH · Engineering Papers

Results for “ion trajectory simulations”

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.

SourceResearch connectionFrom yearThrough year

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

At least 163 records · Page 9

All sources NASA NTRS DOE OSTI

NASA NTRS2014

Features of Afterbody Radiative Heating for Earth Entry

Radiative heating is identified as a major contributor to afterbody heating for Earth entry capsules at velocities above 10 km/s. Because of rate-limited electron-ion recombination processes, a large fraction of the electronically-excited N and O atoms produced in the high temperature/pressure forebody remain as they expand into the afterbody region, which results in significant afterbody radiation. Large radiative heating sensitivities to electron-impact ionization rates and escape factors are identified. Ablation products from a forebody ablator are shown to increase the afterbody radiation by as much as 40%. The tangent-slab radiation transport approach is shown to over-predict the radiative flux by as much as 40% in the afterbody, therefore making the more computationally expensive ray-tracing approach necessary for accurate radiative flux predictions. For the Stardust entry, the afterbody radiation is predicted to be nearly twice as large as the convective heating during the peak heating phase of the trajectory. Comparisons between simulations and the Stardust Echelle observation measurements, which are shown to be dominated by afterbody emission, indicate agreement within 20% for various N and O lines. Similarly, calorimeter measurements from the Fire II experiment are identified as a source of validation data for afterbody radiation. For the afterbody calorimeter measurement closest to the forebody, which experiences the largest afterbody radiative heating component, the convective heating alone is shown to under-predict the measurement, even for the fullycatalytic assumption. Agreement with the measurements is improved with the addition of afterbody radiation. These comparisons with Stardust and Fire II measurements provide validation that the significant afterbody radiation values proposed in this work are legitimate.

Johnston, Christopher O.↗

DOE OSTI2021

Algorithm advances and applications of time‐dependent first‐principles simulations for ultrafast dynamics

Abstract Far from equilibrium phenomenon is a central theme of contemporary material research. Such phenomenon can exhibit itself in atomic structure and dynamics, but very often it also happens as non‐equilibrium phenomenon in the electronic structure. In ab initio material simulation, density functional theory (DFT) has played an essential role in studying electronic ground state problems. For excited states, besides many‐body perturbation theory, another powerful tool is the time dependent DFT (TDDFT) method. In particular, the real‐time TDDFT (rt‐TDDFT) method can be used to simulate many non‐equilibrium phenomena directly. Here we introduce our works on some algorithm advances based on our recently rt‐TDDFT method. This method uses the plane‐wave basis set, and significantly accelerates its efficiency by increasing the time step from 0.1–1 as in traditional methods to 0.2–0.5 fs. The noncollinear magnetic moments and spin–orbit coupling have also been included in our rt‐TDDFT method. Furthermore, a Boltzmann‐TDDFT algorithm has been developed to solve the hot carrier overheating problem in Ehrenfest dynamics, and a natural orbital branching algorithm has been developed to overcome the mean‐field approximation in Ehrenfest dynamics nuclear trajectory, thus allows stochastic multiple paths in chemical reactions. Utilizing these methods, we have studied the photoinduced ultrafast demagnetization, ultrafast phase transition, energy transfer between plasmon and hot carriers, as well as the high‐energy ion implantation and low‐energy atomic diffusion in semiconductors. We believe the tools as the ones introduced here can enable us to study a wide range of phenomena which are of great interest in modern day material research. This article is categorized under: Structure and Mechanism > Computational Materials Science Electronic Structure Theory > Ab Initio Electronic Structure Methods Electronic Structure Theory > Density Functional Theory

Liu, Wen‐Hao↗

DOE OSTI2025

Digital twin framework for PIP-II linac: AI-driven multi-scale modeling from ion source to 800 MeV

The PIP-II superconducting linac at Fermilab is designed to deliver multi-megawatt proton beams for neutrino physics and other high-intensity applications. To expedite commissioning and enhance operational reliability, we have developed an EPICS-based data flow framework that seamlessly integrates digital twins (DT) with physical twins (PT). These digital twins comprise high-fidelity beam dynamics models or data-driven surrogate models connected to their physical counterparts through real-time diagnostics and advanced machine-learning algorithms.Central to this framework is Linac_Gen, an accelerated simulation tool that incorporates convolutional neural networks, random forests, and genetic algorithms to provide up to a tenfold speedup in optimizing the accelerator geometry model. An EPICS translator layer ensures interoperability by efficiently mapping lattice parameters across diverse simulation platforms.Our EPICS-based framework supports multiple operational modes—monitoring, passive learning, closed-loop control, and online learning—covering the entire machine lifecycle. By leveraging HPC resources and multi-objective optimization techniques, the digital twin enables adaptive trajectory correction, real-time fault detection, and predictive modeling of beam stability. This comprehensive approach paves the way for robust, high-intensity operation and data-driven accelerator R&D at Fermilab.

Pathak, Abhishek [Fermilab]↗

DOE OSTI2021

Elastic forces drive nonequilibrium pattern formation in a model of nanocrystal ion exchange

The widely used process of nanocrystal ion exchange operates out of thermodynamic equilibrium and can require mixing components of varying sizes. Here we use theory and computer simulation to study a simple model which captures these two basic features of ion exchange reactions. We show that a strong driving force for exchange among different-sized species creates nonequilibrium patterns within model nanocrystals. We further demonstrate that such patterns can be thermodynamically stable within core/shell nanocrystals. These results help us understand the heterostructures formed in ion-exchanged nanocrystals and suggest strategies for leveraging elasticity to design patterns in nanoscale materials. Chemical transformations, such as ion exchange, are commonly employed to modify nanocrystal compositions. Yet the mechanisms of these transformations, which often operate far from equilibrium and entail mixing diverse chemical species, remain poorly understood. Here we explore an idealized model for ion exchange in which a chemical potential drives compositional defects to accumulate at a crystal’s surface. These impurities subsequently diffuse inward. We find that the nature of interactions between sites in a compositionally impure crystal strongly impacts exchange trajectories. In particular, elastic deformations which accompany lattice-mismatched species promote spatially modulated patterns in the composition. These same patterns can be produced at equilibrium in core/shell nanocrystals, whose structure mimics transient motifs observed in nonequilibrium trajectories. Moreover, the core of such nanocrystals undergoes a phase transition—from modulated to unstructured—as the thickness or stiffness of the shell is decreased. Our results help explain the varied patterns observed in heterostructured nanocrystals produced by ion exchange and suggest principles for the rational design of compositionally patterned nanomaterials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS1993

Stormtime transport of ring current and radiation belt ions

A dynamical guiding-center simulation model is used to study the stormtime ion transport which leads to the formation of the ring current and diffusion in the radiation belts. Representative ions guiding-center motion in response to model storm-associated impulses in the convection electric field is traced for a range of ion mu values. The present numerical results are compared with previously formulated limiting idealization of particle transport in order to assess the limits of validity of these approximations. For ions having drift periods that exceed the duration of the main phase of the storm, their inward transport to form the stormtime ring current is appropriately described as direct convective access. For ions having drift periods comparable to the duration of the main phase of the storm, there is a transition between direct convective access and transport that resembles radial diffusion. Lower-energy ring-current ions at L of about 3 are freshly injected there from open adiabatic trajectories, whereas the higher-energy ring-current population consists of a mixture of freshly injected and previously trapped ions.

Chen, Margaret W.↗

DOE OSTI2020

High precision magnetic field measurement and mapping of the LEReC 180 degree bending magnet using very low field NMR with hall combined probe (140 - 350 gauss)

The Relativistic Heavy Ion Collider (RHIC) at BNL are using the Low Energy RHIC Electron Cooling (LEReC) to conduct experiments that search for the quantum chromodynamic (QCD) critical point. The first ever electron cooling based on the RF acceleration of electron beams was experimentally demonstrated on April 5, 2019 at LEReC at BNL. The first critical step in obtaining successful 3D non-magnetized cooling of the Au ion bunches in the RHIC cooling section was matching the electron beam energy with a relative error less than 5*10 -4 to the ion beam energy. Part of the LEReC beamline is a dipole magnet that bends the electron beam 180 degree. One of the most outstanding measurement challenges is that the dipole field is so low (≈200 G). Most of the existing NMR probes can only measure fields >400 G. Lower signal-to-noise ratio at low fields is requires the use of larger sample volumes. Working with CAYLAR, a NMR probe has been redesigned and optimized for these low field measurements with high resolution. We report the methods, challenges, and results for extensive magnetic field mappings of the 180 dipole magnet. A combination of NMR and Hall sensors has been successfully implemented to measure uniform field regimes inside the magnet center area and non-uniform field regimes at the magnet ends. Detailed measurement and mapping have been performed at five radii and five heights along the beam trajectory. Meanwhile a finite element magnetic modeling simulation of the magnet using Opera software has been performed. The calculated and measured data are compared, and the calculated data is good reference for measured data over long length mapping from magnet edge to center. The measured magnetic measurement data is directly useful for beam instrumentation, diagnostics and operation.

Magnetic Field Measurement↗

NASA NTRS2018

Modeling of Orbital and Attitude Dynamics of Nanosatellites Controlled via Active Electrostatic Charging

The large-scale exploration of airless bodies, such as asteroids and moons, is gaining interest, however it is limited by mobility issues: the lack of atmosphere, low gravity, and unknown soil properties pose difficult challenges for many forms of traditional locomotion. The environment in proximity of these bodies is also electrically charged due to interactions with solar wind and UV radiation. The EGlider (Electrostatic Glider) concept would be able to overcome these mobility issues by leveraging the natural environment, allowing operations in close proximity of the surface, while enabling long duration missions by minimizing propellant consumption. The EGlider is an advanced concept for small satellite mobility and propulsion, which relies on the electric fields naturally present around airless bodies in order to generate forces and torques useful for maneuvering. It does so by extending electrically charged appendages, which enable it to electrostatically soar above the surface. By differentially charging its electrodes it can also produce torques to control its attitude. The charges are maintained by continuous active ion or electron emission from the spacecraft, in order to cancel out the neutralizing influx of charges from surrounding plasma. An investigation of the spacecraft-plasma interaction was carried out. This included studying the effect of electrode geometry and calculating the charge-to-mass ratios required to enable several mission scenarios. Long, thin wire electrodes were identified to be the most power-efficient and would allow power-to-weight ratios achievable with current nanosatellite technologies. High electrode potential represents the main limiting factor for the system design. In order to test the feasibility of active control by means of differential charging, a simple 2D interaction model was developed, and a feedback controller to stabilize the vehicle was tested in a simulation environment. The results confirmed that good performance could be obtained for both position and attitude control. Finally, a dedicated software was developed for future simulation and testing of control strategies for the EGlider. This software allows to study the trajectory and attitude of an arbitrarily configured spacecraft in the proximity of an arbitrarily defined main airless body. The spacecraft can be assembled from basic parts, each with specified electrical, mass and optical properties. Efficient models allow to calculate gravitational and electrical interactions with the rotating main body and the local plasma field, as well as solar radiation pressure effects. Control models can be implemented as simple plug-in functions and easily tested. The preliminary validation campaign showed good matching with the reference cases that have been analyzed.

Corpino, Sabrina↗

DOE OSTI2023

Does Carrier Envelope Phase Affect the Ionization Site in a Neutral Diatomic Molecule?

A recent work shows how to extract the ionization site of a neutral diatomic molecule by comparing Quantum Trajectory Monte Carlo (QTMC) simulations with experimental measurements of the final electron momenta distribution. This method was applied to an experiment using a 40-femtosecond infrared pulse, finding that a downfield atom is roughly twice as likely to be ionized as an upfield atom in a neutral nitrogen molecule. However, an open question remains as to whether an assumption of the zero carrier envelope phase (CEP) used in the above work is still valid for short, few-cycle pulses where the CEP can play a large role. Given experimentalists’ limited control over the CEP and its dramatic effect on electron momenta after ionization, it is desirable to see what influence the CEP may have in determining the ionization site. In this paper, we employ QTMC techniques to simulate strong-field ionization and electron propagation from neutral N2 using an intense 6-cycle laser pulse with various CEP values. Comparing simulated electron momenta to experimental data indicates that the ratio of down-to-upfield ions remains roughly 2:1 regardless of the CEP. This confirms that the ionization site of a neutral molecule is determined predominantly by the laser frequency and intensity, as well as the ground-state molecular wavefunction, and is largely independent of the CEP.

Schimmoller, Alex↗

NASA NTRS2004

Simulations of Solar Wind Plasma Flow Around a Simple Solar Sail

In recent years, a number of solar sail missions of various designs and sizes have been proposed (e.g., Geostorm). Of importance to these missions is the interaction between the ambient solar wind plasma environment and the sail. Assuming a typical 1 AU solar wind environment of 400 km/s velocity, 3.5 cu cm density, ion temperature of approx.10 eV, electron temperature of 40 eV, and an ambient magnetic field strength of 10(exp -4) G, a first order estimate of the plasma interaction with square solar sails on the order of the sizes being considered for a Geostorm mission (50 m x 50 m and 75 m x 75 m corresponding to approx.2 and approx.3 times the Debye length in the plasma) is carried out. First, a crude current balance for the sail surface immersed in the plasma environment and in sunlight was used to estimate the surface potential of the model sails. This gave surface potentials of approx.10 V positive relative to the solar wind plasma. A 3-D, Electrostatic Particle-in-Cell (PIC) code was then used to simulate the solar wind flowing around the solar sail. It is assumed in the code that the solar wind protons can be treated as particles while the electrons follow a Boltzmann distribution. Next, the electric field and particle trajectories are solved self-consistently to give the proton flow field, the electrostatic field around the sail, and the plasma density in 3-D. The model sail was found to be surrounded by a plasma sheath within which the potential is positive compared to the ambient plasma and followed by a separate plasma wake which is negative relative to the plasma. This structure departs dramatically from a negatively charged plate such as might be found in the Earth s ionosphere on the night side where both the plate and its negative wake are contiguous. The implications of these findings are discussed as they apply to the proposed Geostorm solar sail mission.

Garrett, Henry B.↗

NASA NTRSDate not supplied

Wind-Optimal Trajectories for Multirotor eVTOL Aircraft on UAM Missions

Distributed electric propulsion powered electric vertical takeoff and landing aircraft are expected to enable urban air mobility. Low specific energy of onboard lithium-ion polymer batteries and wind conditions impose constraints on flight endurance. Therefore, from the safety and efficiency perspective, planning and flying minimum energy trajectories are important.The primary motivation for this paper is to determine if there is an operational benefit to flying wind-optimal trajectories for short flights (less than 60 miles) in the urban environment. This study employs wind-optimal trajectories for a NASA-proposed conceptual multirotor aircraft flight in the urban environment. The optimal control model presented includes a wind model for quantifying the effect of wind on the trajectory. The optimal control problem is numerically solved using the direct method. Energy consumption and flight duration flying wind-optimal trajectories between origins and destinations in a metropolitan area are compared with energy consumption and flight duration flying on the great-circle paths between the same origin and destination pairs to determine the operational benefit of wind-optimal routing for short flights.The flight duration results for different scenarios are validated using a simulation tool designed and developed by NASA for exploring advanced air traffic management concepts. In summary,this research study suggests that for short flights in an urban environment, the wind-optimal trajectories have a very slight (insignificant) operational advantage over the great-circle trajectories in terms of energy consumption and flight duration. As expected, headwinds and tailwinds along the route of flight have a significant impact on energy consumption and flight duration.

Urban Air Mobility↗

DOE OSTI2023

Stochastic and mixed density functional theory within the projector augmented wave formalism for simulation of warm dense matter

Stochastic density functional theory (DFT) and mixed stochastic-deterministic DFT are burgeoning approaches for the calculation of the equation of state and transport properties in materials under extreme conditions. In the intermediate warm dense matter regime, a state between correlated condensed matter and kinetic plasma, electrons can range from being highly localized around nuclei to delocalized over the whole simulation cell. The plane-wave basis pseudopotential approach is thus the typical tool of choice for modeling such systems at the DFT level. Unfortunately, stochastic DFT methods scale as the square of the maximum plane-wave energy in this basis. To reduce the effect of this scaling and improve the overall description of the electrons within the pseudopotential approximation, we present stochastic and mixed DFT approaches developed and implemented within the projector augmented wave formalism. In conclusion, we compare results between the different DFT approaches for both single-point and molecular dynamics trajectories and present calculations of self-diffusion coefficients of solid density carbon from 1 to 50 eV.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

NASA NTRS2020

Let Dependent Low Dose and Synergistic Inhibition of Human Angiogenesis by Charged Particles: Validation of Micrornas That Drive Inhibition

The microvasculature permeates all tissues at the microscopic level making the whole body a target for Space Radiation. Depending on the trajectory a single heavy ion particle would traverse many micro-vessels as it passes through the body causing a widespread response. Studies on the effects of charged particles on the development of human 3D micro-vessel models show that angiogenesis is inhibited by all ions detectible at a dose of 1.25c Gy. However, there are two distinct mechanisms of inhibition depending on the linear energy transfer (LET). We show here that the inhibition of angiogenesis is demarcated into two groups. The blockage of the early motile stage of micro-vessel development occurs after exposure to light ions (LET 3 KeV/n or lower) whereas the prevention of the later stages of development when the microvessels widen into tubes occurs after exposure to heavy ions (LET 8 KeV/n or higher). Strikingly, the combined effect of each ion has a synergistic effect detectible as low as 0.6 cGy. The low fluence of these doses indicates a bystander effect where the response is transmitted to other cells by the secretion of molecules by targeted cells. Prime candidates are microRNAs (miRNAs) and we have previously identified three angiogenesis-associated miRNA regulators from serum, cardiac tissue, and micro-vessels exposed to a simplified Galactic Cosmic Radiation (GCR) simulation – a mixture of heavy and light ions found in Space. Synthetic molecules (Antagomirs) designed to inhibit the predicted miRNAs were successfully used to reverse the inhibition of angiogenesis by Space Radiation. With the simple application of the antagomirs as a countermeasure we have for the first time shown an amelioration of the complex response(s) of angiogenesis to Space Radiation. This opens a new avenue of possible mitigation strategies based on the regulation and understanding of master regulators. These results not only have impact with Space Radiation research, but also provide clues of possible therapeutics for terrestrial diseases in the clinical setting.

Yen-Ruh Wuu↗

NASA NTRS1999

Computer Simulation Studies of Ion Channel Gating: Characteristics of the M2 Channel of Influenza-A Virus in a Phospholipid Bilayer

The 25 amino acids long, transmembrane fragment of the Influenza virus M2 protein forms a homotetrameric channel that transports protons across lipid bilayers. It has been postulated that high efficiency and selectivity of this process is due to gating by four histidine residues that occlude the channel lumen in the closed state. Two mechanisms of gating have been postulated. In one mechanism, the proton is "shuttled" through the gate by attaching to the delta nitrogen atom on the extracellular side of the imidazole ring, followed by the release of the proton attached to the epsilon nitrogen atom on the opposite side. In the second mechanism, the four histidines move away from each other due to electrostatic repulsion upon protonation, thus opening the gate sufficiently that a wire of water molecules can penetrate the gate. Then, protons are transported by "hopping" along the wire. In this paper, both mechanisms are evaluated in a series of molecular dynamics simulations by investigating stability of different protonation states of the channel that are involved in these mechanisms. For the shuttle mechanism, these are states with all epsilon protonated histidines, one biprotonated residue or one histidine protonated in the delta position. For the gate opening mechanism, this is the state in which all four histidines are biprotonated. In addition, a state with two biprotonated histidines is considered. For each system, composed of the protein channel embedded in phospholipid bilayer located between two water lamellae, a molecular dynamics trajectory of approximately 1.3 ns (after equilibration) was obtained. It is found that the states involved in the shuttle mechanism are stable during the simulations. Furthermore, the orientations and dynamics of water molecules near the gate are conducive to proton transfers involved in the shuttle. In contract, the fully biprotonated state, implicated in the gate opening mechanism, is not stable and the channel looses its structural integrity. If only two histidines are biprotonated the channel deforms but remains intact with the gate mostly closed. In summary, the results of this study lend support to the shuttle mechanism but not to the gate opening mechanism of proton gating in M2.

Schweighofer, Karl J.↗

DOE OSTI2023

Time-dependent SOLPS-ITER simulations of the tokamak plasma boundary for model predictive control using SINDy *

Abstract Time-dependent SOLPS-ITER simulations have been used to identify reduced models with the sparse identification of nonlinear dynamics (SINDy) method and develop model-predictive control of the boundary plasma state using main ion gas puff actuation. A series of gas actuation sequences are input into SOLPS-ITER to produce a dynamic response in upstream and divertor plasma quantities. The SINDy method is applied to identify reduced linear and nonlinear models for the electron density at the outboard midplane n e , s e p O M P and the electron temperature at the outer divertor T e , s e p d i v . Note that T e , s e p d i v is not necessarily the peak value of T e along the divertor. The identified reduced models are interpretable by construction (i.e. not black box), and have the form of coupled ordinary differential equations. Despite significant noise in T e , s e p d i v , the reduced models can be used to predict the response over a range of actuation levels to a maximum deviation of 0.5% in n e , s e p O M P and 5%–10% in T e , s e p d i v for the cases considered. Model retraining using time history data triggered by a preset error threshold is also demonstrated. A model predictive control strategy for nonlinear models is developed and used to perform feedback control of a SOLPS-ITER simulation to produce a setpoint trajectory in n e , s e p O M P using the integrated plasma simulator framework. The developed techniques are general and can be applied to time-dependent data from other boundary simulations or experimental data. Ongoing work is extending the approach to model identification and control for divertor detachment, which will present transient nonlinear behavior from impurity seeding, including realistic latency and synthetic diagnostic signals derived from the full SOLPS-ITER output.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Overview of the TCV tokamak program: Scientific progress and facility upgrades

The TCV tokamak is augmenting its unique historical capabilities (strong shaping, strong electron heating) with ion heating, additional electron heating compatible with high densities, and variable divertor geometry, in a multifaceted upgrade program designed to broaden its operational range without sacrificing its fundamental flexibility. The TCV program is rooted in a three-pronged approach aimed at ITER support, explorations towards DEMO, and fundamental research. A 1 MW, tangential neutral beam injector (NBI) was recently installed and promptly extended the TCV parameter range, with record ion temperatures and toroidal rotation velocities and measurable neutral-beam current drive. ITER-relevant scenario development has received particular attention, with strategies aimed at maximizing performance through optimized discharge trajectories to avoid MHD instabilities, such as peeling-ballooning and neoclassical tearing modes. Experiments on exhaust physics have focused particularly on detachment, a necessary step to a DEMO reactor, in a comprehensive set of conventional and advanced divertor concepts. The specific theoretical prediction of an enhanced radiation region between the two X-points in the low-field-side snowflake-minus configuration was experimentally confirmed. Fundamental investigations of the power decay length in the scrape-off layer (SOL) are progressing rapidly, again in widely varying configurations and in both D and He plasmas; in particular, the double decay length in L-mode limited plasmas was found to be replaced by a single length at high SOL resistivity. Experiments on disruption mitigation by massive gas injection and electron-cyclotron resonance heating (ECRH) have begun in earnest, in parallel with studies of runaway electron generation and control, in both stable and disruptive conditions; a quiescent runaway beam carrying the entire electrical current appears to develop in some cases. Developments in plasma control have benefited from progress in individual controller design and have evolved steadily towards controller integration, mostly within an environment supervised by a tokamak profile control simulator. TCV has demonstrated effective wall conditioning with ECRH in He in support of the preparations for JT-60SA operation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

Bilayer Ion Trap Design for 2D Arrays

Junctions are fundamental elements that support qubit locomotion in two-dimensional ion trap arrays and enhance connectivity in emerging trapped-ion quantum computers. In surface ion traps they have typically been implemented by shaping radio frequency (RF) electrodes in a single plane to minimize the disturbance to the pseudopotential. However, this method introduces issues related to RF lead routing that can increase power dissipation and the likelihood of voltage breakdown. Here, we propose and simulate a novel two-layer junction design incorporating two perpendicularly rotoreflected (rotated, then reflected) linear ion traps. The traps are vertically separated, and create a trapping potential between their respective planes. The orthogonal orientation of the RF electrodes of each trap relative to the other provides perpendicular axes of confinement that can be used to realize transport in two dimensions. While this design introduces manufacturing and operating challenges, as now two separate structures have to be precisely positioned relative to each other in the vertical direction and optical access from the top is obscured, it obviates the need to route RF leads below the top surface of the trap and eliminates the pseudopotential bumps that occur in typical junctions. Here in this paper the stability of idealized ion transfer in the new configuration is demonstrated, both by solving the Mathieu equation analytically to identify the stable regions and by numerically modeling ion dynamics. Our novel junction layout has the potential to enhance the flexibility of microfabricated ion trap control to enable large-scale trapped-ion quantum computing.

42 ENGINEERING↗

DOE OSTI2025

Molecular-Level Insights into the NMR Relaxivity of Gadobutrol Using Quantum and Classical Molecular Simulations

MRI is an indispensable diagnostic tool in modern medicine; however, understanding the molecular-level processes governing NMR relaxation of water in the presence of MRI contrast agents remains a challenge, hindering the molecular-guided development of more effective contrast agents. By using quantum-based polarizable force fields, the first-of-its-kind molecular dynamics (MD) simulations of Gadobutrol are reported where the 1 H NMR longitudinal relaxivity r 1 of the aqueous phase is determined without any adjustable parameters. The MD simulations of r 1 dispersion (i.e., frequency dependence) show good agreement with measurements at frequencies of interest in clinical MRI. Importantly, the simulations reveal key insights into the molecular level processes leading to r 1 dispersion by decomposing the NMR dipole–dipole autocorrelation function G(t) into a discrete set of molecular modes, analogous to the eigenmodes of a quantum harmonic oscillator. The molecular modes reveal important aspects of the underlying mechanisms governing r 1 , such as its multiexponential nature and the importance of the second eigenmodal decay. By simply analyzing the MD trajectories on a parameter-free approach, the Gadobutrol simulations show that the outer-shell water contributes ∼50% of the total relaxivity r 1 compared to the inner-shell water, in contrast to simulations of (nonchelated) gadolinium-aqua where the outer shell contributes only ∼15% of r 1 . The deviation between simulations and measurements of r 1 below clinical MRI frequencies is used to determine the low-frequency electron-spin relaxation time for Gadobutrol, in good agreement with independent studies.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS1994

Molecular Dynamics of a Water-Lipid Bilayer Interface

We present results of molecular dynamics simulations of a glycerol 1-monooleate bilayer in water. The total length of analyzed trajectories is 5ns. The calculated width of the bilayer agrees well with the experimentally measured value. The interior of the membrane is in a highly disordered fluid state. Atomic density profile, orientational and conformational distribution functions, and order parameters indicate that disorder increases toward the center of the bilayer. Analysis of out-of-plane thermal fluctuations of the bilayer surfaces occurring at the time scale of the present calculations reveals that the distribution of modes agrees with predictions of the capillary wave model. Fluctuations of both bilayer surfaces are uncorrelated, yielding Gaussian distribution of instantaneous widths of the membrane. Fluctuations of the width produce transient thinning defects in the bilayer which occasionally span almost half of the membrane. The leading mechanism of these fluctuations is the orientational and conformational motion of head groups rather than vertical motion of the whole molecules. Water considerably penetrates the head group region of the bilayer but not its hydrocarbon core. The total net excess dipole moment of the interfacial water points toward the aqueous phase, but the water polarization profile is non-monotonic. Both water and head groups significantly contribute to the surface potential across the interface. The calculated sign of the surface potential is in agreement with that from experimental measurements, but the value is markedly overestimated. The structural and electrical properties of the water-bilayer system are discussed in relation to membrane functions, in particular transport of ions and nonelectrolytes across membranes.

Wilson, Michael A.↗