SEARCH · Engineering Papers

Results for “Interaction force”

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 865 records · Page 48

All sources NASA NTRS DOE OSTI

DOE OSTI2026

Shadow molecular dynamics for flexible multipole models

Shadow molecular dynamics provide an efficient and stable atomistic simulation framework for flexible charge models with long-range electrostatic interactions. Shadow molecular dynamics simulations are driven by approximate “shadow” Born–Oppenheimer potentials for which the exact charges and forces are directly accessible without relying on costly (and approximate) iterative solvers. While previous implementations have been limited to atomic monopole charge distributions, we extend this approach to flexible multipole models. We derive detailed expressions for the shadow energy functions, potentials, and force terms, explicitly incorporating monopole–monopole, dipole–monopole, and dipole–dipole interactions. In our formulation, both atomic monopoles and atomic dipoles are treated as extended dynamical variables alongside the propagation of the nuclear degrees of freedom. We demonstrate that introducing the additional dipole degrees of freedom preserves the stability and accuracy previously seen in monopole-only shadow molecular dynamics simulations. In addition, we present a shadow molecular dynamics scheme where the monopole charges are held fixed while the dipoles remain flexible. Our extended shadow dynamics provide a framework for stable, computationally efficient, and versatile molecular dynamics simulations involving long-range interactions between flexible multipoles. This is of particular current interest in combination with machine-learned interatomic potentials, including long-range electrostatic interactions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS1978

The driving force for magnetospheric convection

Viscously driven magnetospheric models, as well as a model involving interconnection between the geomagnetic field and the magnetic field in the solar wind, have been proposed to describe the driving force for magnetospheric convection. Lack of a satisfactory theory for the interconnection in the latter model and, in the case of the viscous interaction models, inadequacies in predicting the quantity of the driving force, make these two classes of models less than successful. Accordingly, a mechanically driven magnetospheric model is proposed: solar wind plasma enters the magnetosphere around the neutral points, covers the inner surface of the magnetopause and subsequently expands, driving convection as it escapes from the open tail.

Johnson, F. S.↗

NASA NTRS2011

Overview of Experimental Investigations for Ares I Launch Vehicle Development

Another concern for the vehicle during its design trajectory was the separation of the first stage solid rocket booster from the upper stage component after it had depleted its solid fuel propellant. There has been some concern about the interstage of the first stage from clearing the nozzle of the J2-X engine. A detailed separation aerodynamic wind tunnel investigation was conducted in the AEDC VKF Tunnel A to help to investigate the interaction aerodynamic effects5. A comparison of the separation plane details between the Ares I architecture and the Ares I-X demonstration flight architecture is shown in figure 12. The Ares I design requires a more complex separation sequence and requires better control in order to avoid contact with the nozzle of the upper stage engine. The interstage, which houses the J2-X engine for the Ares I vehicle, must be able to separate cleanly to avoid contact of the J2-X engine. There is only about approximately 18 inches of buffer inside the interstage on each size of the nozzle so this is a challenging controlled separation event. This complex experimental investigation required two separate Ares I models (upper stage and first stage with interstage attached) with independent strain gauge balances installed in each model. It also required the Captive Trajectory System (CTS) that was needed to precisely locate the components in space relative to each other to fill out the planned test matrix. The model setup in the AEDC VKF Tunnel A is shown in figure 13. The CTS remotely positioned the first stage at the required x, y, and z positions and was able to provide interactions within 0.2" of the upper stage. A sample of the axial force on the first stage booster is shown in figure 14. These results, as a function of separation distance between the two stages, are compared to pre-test CFD results. Since this is a very challenging, highly unsteady flow field for CFD to correctly model, the experimental results have been utilized by GN&C discipline to more accurately represent the interaction aerodynamics. In addition to the integrated forces and moments obtained from the test, flow visualization data was obtained from this test in the form of Schlieren photographs, as shown in figure 15, which show the shock structure and interaction effects after the two stages separate during flight. This separation test was crucial in the successful flight test of the Ares I-X vehicle and provided the GN&C discipline with the unpowered proximity aerodynamic effect for a separation of the Ares I vehicle.

Tomek, William G.↗

DOE OSTI2023

Coexistence of Merons with Skyrmions in the Centrosymmetric Van Der Waals Ferromagnet Fe 5– x GeTe 2

Fe 5-x GeTe 2 is a centrosymmetric, layered van der Waals (vdW) ferromagnet that displays Curie temperatures T c (270–330 K) that are within the useful range for spintronic applications. However, little is known about the interplay between its topological spin textures (e.g., merons, skyrmions) with technologically relevant transport properties such as the topological Hall effect (THE) or topological thermal transport. Here, via high-resolution Lorentz transmission electron microscopy, it is shown that merons and anti-meron pairs coexist with Néel skyrmions in Fe 5-x GeTe 2 over a wide range of temperatures and probe their effects on thermal and electrical transport. A THE is detected, even at room T, that senses merons at higher T's, as well as their coexistence with skyrmions as T is lowered, indicating an on-demand thermally driven formation of either type of spin texture. Remarkably, an unconventional THE is also observed in absence of Lorentz force, and it is attributed to the interaction between charge carriers and magnetic field-induced chiral spin textures. Here these results expose Fe 5-x GeTe 2 as a promising candidate for the development of applications in skyrmionics/meronics due to the interplay between distinct but coexisting topological magnetic textures and unconventional transport of charge/heat carriers.

two-dimensional magnets↗

DOE OSTI2025

Layers Can Be Deceiving: A Hopping Model for Small Molecule Diffusion in TATB Crystal

Sorption of small molecule gases in materials can play a significant role in their long‐term stability and compatibility within multi‐material assemblies. While many material properties of the insensitive high explosive TATB (1,3,5‐triamino‐2,4,6‐trinitrobenzene) are well understood, very little is known regarding its permeability to gases. TATB crystal exhibits a graphitic‐like layered packing structure that evokes a mental schema in which the layers form nanoscopic channels, but it is unclear whether this structure promotes gas transport. Here, we use molecular dynamics (MD) simulations to predict transport of small molecules through TATB single crystal. An approach to fit classical force fields is developed to model TATB interactions with H 2 O, He, Ne, and Ar, which is then combined with steered MD to probe gas transport along selected directions in the crystal. We find that small molecule transport occurs via a hopping mechanism that exhibits distinct jumps between interstitial sites and is substantially faster normal to the layers as compared to through them. This result stems from the finding that intralayer junctions between adjacent TATB molecules are the most stable interstitial sites and that energetic barriers are lower for hopping between adjacent layers. An empirical model for diffusion rate based on the MD data shows that the rate decays exponentially with increasing molecular radius and is negligibly small for all molecules larger than He, including common atmospheric gases. These findings have implications for the interpretation of experiments that measure surface area, material response to extreme conditions, and are expected to help constrain models for material aging.

36 MATERIALS SCIENCE↗

DOE OSTI2019

Editorial: The World Is Nonlocal

Nonlocal modeling has come a long way. Researchers in the continuum mechanics and computational mechanics communities increasingly recognize that nonlocality is critical in realistic mathematical models of many aspects of the physical world. Physical interaction over a finite distance is fundamental at the atomic and nanoscale level, in which atoms and molecules interact through multibody potentials. Long-range forces partially determine the mechanics of surfaces and the behavior of dissolved molecules and suspended particles in a fluid. Nonlocality is therefore a vital feature of any continuum model that represents these physical systems at small length scales.

99 GENERAL AND MISCELLANEOUS↗

DOE OSTI2025

Atomistic Insights into the Reactive Diffusion of CO 2 in Guanidine-Based Facilitated Transport Membranes

The pressing need to address climate change has led to significant advancements in carbon dioxide (CO 2 ) capture technologies. Notably, facilitated transport membranes (FTMs) are distinguished by their exceptional selectivity and permeance, attributed to their reversible chemical reactions with CO 2 . This study, for the first time, sheds light on the reactive diffusion mechanism of CO 2 in FTMs, utilizing 1,1,3,3-tetramethylguanidine (TMG) as a mobile carrier. Specifically, state-ofthe- art molecular dynamics (MD) simulations, augmented by a reparameterized reactive force field (ReaxFF) capable of describing atomistic interactions and reaction pathways, are conducted to investigate the transport of CO 2 in TMG. The analysis of mean squared displacement (MSD) and diffusion coefficients reveals a clear hierarchy in the mobility of reaction components. Our findings highlight a unique hopping diffusion mechanism between bicarbonate ions and TMG molecules, increasing the diffusivity of reacted CO 2 by 1.4 times. The hopping events observed not only enhance our understanding of molecular mobility but also offer a means to boost the performance of FTMs in CO 2 capture applications. Overall, this research lays the groundwork for the future design of FTMs with optimal carrier properties.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2018

How Useful Are Common Simulants of Chemical Warfare Agents at Predicting Adsorption Behavior?

We report that nanoporous materials such as metal–organic frameworks (MOFs) have attractive properties for selective capture of chemical warfare agents (CWAs). For obvious reasons, most research on adsorption of CWAs is performed with simulant molecules rather than real agents. This paper examines how effectively common CWA simulants mimic the adsorption properties of sarin and soman. To this end, we perform molecular simulations in the dilute adsorption limit for four simulants [dimethyl methylphosphonate (DMMP), diethyl chlorophosphate (DCP), diisopropyl fluorophosphate, and dimethyl p-nitrophenyl phosphate (DMNP)] and sarin and soman in a set of 2969 MOFs with experimentally known crystal structures. To establish the robustness of the conclusions with respect to the force field used in these simulations, each system was examined with two independent force fields, a “generic” force field and a density functional theory (DFT)-derived force field we established based on extensive dispersion-corrected DFT calculations of adsorption in the well-known MOF UiO-66. Our results show that when judging the performance of adsorbents using the heat of adsorption, DCP and DMMP are the best simulants for the adsorption of sarin, while DMNP is the best simulant for soman. The adsorption properties of DCP or DMMP show a strong correlation with sarin over a range of MOFs, but the correlation between DMNP and soman is considerably weaker. Finally, comparisons of results with both force fields indicate that our main conclusions are robust with respect to the force field used to define adsorbate–MOF interactions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Uncertainty-aware molecular dynamics from Bayesian active learning for phase transformations and thermal transport in SiC

Abstract Machine learning interatomic force fields are promising for combining high computational efficiency and accuracy in modeling quantum interactions and simulating atomistic dynamics. Active learning methods have been recently developed to train force fields efficiently and automatically. Among them, Bayesian active learning utilizes principled uncertainty quantification to make data acquisition decisions. In this work, we present a general Bayesian active learning workflow, where the force field is constructed from a sparse Gaussian process regression model based on atomic cluster expansion descriptors. To circumvent the high computational cost of the sparse Gaussian process uncertainty calculation, we formulate a high-performance approximate mapping of the uncertainty and demonstrate a speedup of several orders of magnitude. We demonstrate the autonomous active learning workflow by training a Bayesian force field model for silicon carbide (SiC) polymorphs in only a few days of computer time and show that pressure-induced phase transformations are accurately captured. The resulting model exhibits close agreement with both ab initio calculations and experimental measurements, and outperforms existing empirical models on vibrational and thermal properties. The active learning workflow readily generalizes to a wide range of material systems and accelerates their computational understanding.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Coulomb-free 1S0 p − p scattering length from the quasi-free p + d → p + p + n reaction and its relation to universality

Abstract The Coulomb-free 1 S 0 proton-proton ( p - p ) scattering length relies heavily on numerous and distinct theoretical techniques to remove the Coulomb contribution. Here, it has been determined from the half-off-the-energy-shell p - p scattering cross section measured at center-of-mass energies below 1 MeV using the quasi-free p + d → p + p + n reaction. A Bayesian data-fitting approach using the expression of the s-wave nucleon-nucleon scattering cross section returned a p - p scattering length $${a}_{pp}=-18.1{7}_{-0.58}^{+0.52}{| }_{stat}\pm 0.0{1}_{syst}$$ a p p = − 18.1 7 − 0.58 + 0.52 ∣ s t a t ± 0.0 1 s y s t fm and effective range r 0 = 2.80 ± 0.05 s t a t ± 0.001 s y s t fm. A model based on universality concepts has been developed to interpret this result. It accounts for the short-range interaction as a whole, nuclear and residual electromagnetic, according to what the s-wave phase-shift δ does in the description of low-energy nucleon-nucleon scattering data. We conclude that our parameters are representative of the short-range physics and propose to assess the charge symmetry breaking of the short-range interaction instead of just the nuclear interaction. This is consistent with the current understanding that the charge dependence of nuclear forces is due to different masses of up-down quarks and their electromagnetic interactions. This achievement suggests that these properties have a lesser than expected impact in the context of the charge symmetry breaking.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2019

Molecular mechanism and binding free energy of doxorubicin intercalation in DNA

The intercalation process of binding doxorubicin (DOX) in DNA is studied by extensive MD simulations. Many molecular factors that control the binding affinity of DOX to DNA to form a stable complex are inspected and quantified by employing continuum solvation models for estimating the binding free energy. The modified MM-PB(GB)SA methodology provides a complete energetic profile of Δ G ele , Δ G vDW , Δ G polar , Δ G non-polar , T Δ S total , Δ G deform , Δ G con , and Δ G ion . To identify the sequence specificity of DOX, two different DNA sequences, d(CGATCG) or DNA1 and d(CGTACG) or DNA2, with one molecule (1 : 1 complex) or two molecule (2 : 1 complex) configurations of DOX were selected in this study. Our results show that the DNA deformation energy (Δ G deform ), the energy cost from translational and rotational entropic contributions ( T Δ S tran+rot ), the total electrostatic interactions (Δ G polar-PB/GB + Δ G ele ) of incorporation, the intramolecular electrostatic interactions (Δ G ele ) and electrostatic polar solvation interactions (Δ G polar-PB/GB ) are all unfavorable to the binding of DOX to DNA. However, they are overcome by at least five favorable interactions: the van der Waals interactions (Δ G vDW ), the non-polar solvation interaction (Δ G non-polar ), the vibrational entropic contribution ( T Δ S vib ), and the standard concentration dependent free energies of DOX (Δ G con ) and the ionic solution (Δ G ion ). Specifically, the van der Waals interaction appears to be the major driving force to form a stable DOX–DNA complex. We also predict that DOX has stronger binding to DNA1 than DNA2. The DNA deformation penalty and entropy cost in the 2 : 1 complex are less than those in the 1 : 1 complex, thus they indicate that the 2 : 1 complex is more stable than the 1 : 1 complex. Here, we have calculated the total binding free energy (BFE) (Δ G t-sim ) using both MM-PBSA and MM-GBSA methods, which suggests a more stable DOX–DNA complex at lower ionic concentration. The calculated BFE from the modified MM-GBSA method for DOX–DNA1 and DOX–DNA2 in the 1 : 1 complex is –9.1 and –5.1 kcal mol –1 respectively. The same quantities from the modified MM-PBSA method are –12.74 and –8.35 kcal mol –1 respectively. The value of the total BFE Δ G t-sim in the 1 : 1 complex is in reasonable agreement with the experimental value of –7.7 ± 0.3 kcal mol –1 .

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Structural evolution of two-dimensional silicates using a “bond-switching” algorithm

Silicates are the most abundant materials in the earth's crust. In recent years, two-dimensional (2D) versions of them grown on metal supports (known as bilayer silicates) have allowed their study in detail down to the atomic scale. These structures are self-containing. Furthermore, they are not covalently bound to the metal support but interact with it through van der Waals forces. Like their three-dimensional counterparts, the 2D-silicates can form both crystalline and vitreous structures. Furthermore, the interconversion between vitreous to crystalline structures has been experimentally observed at the nanoscale. While theoretical work has been carried out to try to understand these transformations, a limitation for ab initio methods, and even molecular dynamics methods, is the computational cost of studying large systems and long timescales. In this work, we present a simple and computationally inexpensive approach, that can be used to represent the evolution of bilayer silicates using a bond-switching algorithm. This approach allows reaching equilibrium ring size distributions as a function of a parameter that can be related to the ratio between temperature and the energy required for the bond-switching event. The ring size distributions are compared to experimental data available in the literature.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

DOE OSTI2025

Polarization and dynamic phases of aligning active matter in periodic obstacle arrays

Here, we numerically examine a system of monodisperse self-propelled particles interacting with each other via simple steric forces and aligning torques moving through a periodic array of obstacles. Without obstacles, this system shows a transition to a polarized or aligned state for critical alignment parameters. In the presence of obstacles, there is still a polarization transition, but for dense enough arrays, the polarization is locked to the symmetry directions of the substrate. When the obstacle array is made anisotropic, at low densities the particles can form a quasi-isotropic state where the system can be polarized in any of the dominant symmetry directions. For intermediate anisotropy, the particles self-organize into a coherent lane state with one-dimensional polarization. In this phase, a small number of highly packed lanes are adjacent to less dense lanes that have the same polarization, but lanes further away can have the opposite polarization, so that global polarization is lost. For the highest anisotropy, hopping between lanes is suppressed, and the system forms uniformly dense uncoupled but polarized lanes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2021

Design of nematic liquid crystals to control microscale dynamics

The dynamics of small particles, both living such as swimming bacteria and inanimate, such as colloidal spheres, has fascinated scientists for centuries. If one could learn how to control and streamline their chaotic motion, that would open technological opportunities in the transformation of stored or environmental energy into systematic motion, with applications in micro-robotics, transport of matter, guided morphogenesis. This review presents an approach to command microscale dynamics by replacing an isotropic medium with a liquid crystal. Orientational order and associated properties, such as elasticity, surface anchoring, and bulk anisotropy, enable new dynamic effects, ranging from the appearance and propagation of particle-like solitary waves to self-locomotion of an active droplet. By using photoalignment, the liquid crystal can be patterned into predesigned structures. In the presence of the electric field, these patterns enable the transport of solid and fluid particles through nonlinear electrokinetics rooted in anisotropy of conductivity and permittivity. Director patterns command the dynamics of swimming bacteria, guiding their trajectories, polarity of swimming, and distribution in space. This guidance is of a higher level of complexity than a simple following of the director by rod-like microorganisms. Namely, the director gradients mediate hydrodynamic interactions of bacteria to produce an active force and collective polar modes of swimming. The patterned director could also be engraved in a liquid crystal elastomer. When an elastomer coating is activated by heat or light, these patterns produce a deterministic surface topography. The director gradients define an activation force that shapes the elastomer in a manner similar to the active stresses triggering flows in active nematics. Finally, the patterned elastomer substrates could be used to define the orientation of cells in living tissues. The liquid-crystal guidance holds a major promise in achieving the goal of commanding microscale active flows.

36 MATERIALS SCIENCE↗

DOE OSTI2018

Manipulation of individual superconducting vortices and stick-slip motion in periodic pinning arrays

We numerically examine the manipulation of superconducting vortices interacting with a moving trap representing a magnetic force tip translating across a superconducting sample containing a periodic array of pinning sites. As a function of the tip velocity and coupling strength, we find five distinct dynamic phases, including a decoupled regime where the vortices are dragged a short distance within a pinning site, an intermediate coupling regime where vortices in neighboring pinning sites exchange places, an intermediate trapping regime where individual vortices are dragged longer distances and exchange modes of vortices occur in the surrounding pins, an intermittent multiple trapping regime where the trap switches between capturing one or two vortices, and a strong coupling regime in which the trap permanently captures and drags two vortices. In some regimes we observe the counterintuitive behavior that slow moving traps couple less strongly to vortices than faster moving traps; however, the fastest moving traps are generally decoupled. The different phases can be characterized by the distances the vortices are displaced and the force fluctuations exerted on the trap. We find different types of stick-slip motion depending on whether vortices are moving into and out of pinning sites, undergoing exchange, or performing correlated motion induced by vortices outside of the trap. Finally, our results are general to the manipulation of other types of particle-based systems interacting with periodic trap arrays, such as colloidal particles or certain types of frictional systems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2024

Compiled properties of nucleonic matter and nuclear and neutron star models from nonrelativistic and relativistic interactions

Here, this paper compiles the model parameters and zero-temperature properties of an extensive collection of published theoretical nuclear interactions, including 255 nonrelativistic (Skyrme-like) forces, 270 relativistic mean field (RMF) and point-coupling (RMF-PC) forces, and 13 Gogny-like forces. This forms the most exhaustive tabulation of model parameters to date. The properties of uniform symmetric matter and pure neutron matter at the saturation density are determined. Symmetry properties found from the second-order term of a Taylor expansion in neutron excess are compared with the energy difference of pure neutron and symmetric nuclear matter at the saturation density. Selected liquid-droplet model parameters, including the surface tension and surface symmetry energy, are determined for semi-infinite surfaces. Theoretical liquid droplet model neutron skin thicknesses of the neutron-rich closed-shell nuclei 48 Ca and 208 Pb are compared with published theoretical Hartree-Fock and experimental results. A similar comparison is made between theoretical liquid-droplet and experimental values of dipole polarizabilities. In addition, radii, binding energies, moments of inertia and tidal deformabilities of 1.2⁢M ⊙ , 1.4⁢M ⊙ and 1.6⁢M ⊙ neutron stars are computed. An extensive correlation analysis of bulk matter, nuclear structure, and low-mass neutron star properties is performed and compared with nuclear experiments and astrophysical observations.

nuclear astrophysics↗

DOE OSTI2022

Models of ultraheavy dark matter visible to macroscopic mechanical sensing arrays

In recent years, the sensitivity of optomechanical force sensors has improved leading to increased interest in using these devices as particle detectors. In this study we consider scenarios where dark matter with mass close to the Planck scale may be probed by a large array of optomechanical accelerometers. We motivate a macroscopic mechanical search for ultraheavy dark matter, exemplified by the efforts of the Windchime Collaboration, by providing a survey of the model space that would be visible to such a search. We consider two classes of models, one that invokes a new long-range unscreened force and another that is only gravitationally interacting. We identify significant regions of well-motivated, potentially visible parameter space for versatile models such as Q-balls, composite dark matter, relics of gravitational singularities, and gravitationally produced ultraheavy particles.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2024

Dark matter relic density in strongly interacting dark sectors with light vector mesons

Stable dark matter particles may arise as pseudo-Goldstone bosons from the confinement of dark quarks interacting via a non-Abelian gauge force. Their relic abundance is determined not by annihilations into visible particles but by dark pion number-changing processes within the dark sector, such as 3 π D → 2 π D . However, if the dark vector mesons ρ D are light enough for 3 π D → π D ρ D annihilations to be kinematically allowed, this process dominates and significantly delays freeze-out. As a result, the preferred dark matter mass scale increases and bounds from the Bullet Cluster can be evaded. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗