Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Molecular statics”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Performant automatic differentiation of local coupled cluster theories: Response properties and ab initio molecular dynamics

In this work, we introduce a differentiable implementation of the local natural orbital coupled cluster (LNO-CC) method within the automatic differentiation framework of the PySCFAD package. The implementation is comprehensively tuned for enhanced performance, which enables the calculation of first-order static response properties on medium-sized molecular systems using coupled cluster theory with single, double, and perturbative triple excitations [CCSD(T)]. We evaluate the accuracy of our method by benchmarking it against the canonical CCSD(T) reference for nuclear gradients, dipole moments, and geometry optimizations. In addition, we demonstrate the possibility of property calculations for chemically interesting systems through the computation of bond orders and Mössbauer spectroscopy parameters for a [NiFe]-hydrogenase active site model, along with the simulation of infrared spectra via ab initio LNO-CC molecular dynamics for a protonated water hexamer.

Chemistry↗
NASA NTRS2003

Photonic Interrogation and Control of Nano Processes

My research activities for the summer of 2003 consisted of two projects: One project was concerned with determining a method for predicting the static and dynamic assembly properties of nano-structures using laser tweezers. The other project was to investigate the generation of Laguerre-Gaussian modes using a spatial light modulator incorporated into an optical tweezers system. Concerning the first project, I initially pursued the approach suggested by my NASA colleague Dr. Art Decker. This approach involved mimicking the model of the structure of atomic nucleus for the assembly of 1 to 100 atoms using allowed quadruple transitions induced by orbital angular momentums of a Laguerre- Gaussian (Doughnut) laser mode. After realizing the inaptness of the nuclear model with the nanostructure model as far as the binding forces and transitions were concerned, I focused on using quantum dot modei. This model was not attuned also for the host lattice influences the electronic structure of the quantum dot. Thus one other option that I decided to pursue was the approach of molecular quantum mechanics. In this approach the nanostructure is treated as a large (10-100 nm) molecule constructed from single element or multi-elements. Subsequent to consultation with Dr. Fred Morales, a chemical engineer at NASA GRC, and Dr. David Ball, a computational chemist at Cleveland State University, I acquired a molecular-quantum computation software, Hyperchem 7.0. This software allows simulation of different molecular structures as far as their static and dynamic behaviors are concerned. The time that I spent on this project was about eight weeks. Once this suitable approach was identified, I realized the need to collaborate with a computational quantum chemist to pursue searching for stable nanostructures in the range of 10-100 nm that we can be assembled using laser tweezers. The second project was about generating laser tweezers that possess orbital angular momentum. As shown, we were able to generate laser tweezers modes of different orbital angular momentum using a spatial light modulator incorporated into a laser tweezers system. The motivation for investigating these types of modes stems from being able to spin particles at high speeds and also to orient two particles in separate traps and then join them together. Also, there has been recent intense interest on fundamental physics research on orbital angular momentum of light. The fact that circularly polarized light may have associated with it angular momentum that relates to the spin of individual photons (spin 0 for the plane polarized light, spin +1 for the right-circularly polarized light and spin -1 for the left-circularly polarized light) was first demonstrated by Beth in 1936. Orbital angular momentum is, however, distinct from spin in that the spin angular momentum of light is intrinsically linked to the behavior of the electric field in the light whereas orbital angular momentum is a consequence of inclined wavefronts. In 1992 L. Allen, et al showed that the Laguerre-Gaussian (LG) modes could possess well-defined orbital angular momentum that can exceed 1 planck's constant, i.e. l plancks constant per photon, where l is the azimuthal index of the mode.

Jassemnejad, Baha↗
DOE OSTI2024

GW with hybrid functionals for large molecular systems

A low-cost approach for stochastically sampling static exchange during time-dependent Hartree–Fock-type propagation is presented. This enables the use of an excellent hybrid density functional theory (DFT) starting point for stochastic GW quasiparticle energy calculations. Generalized Kohn–Sham molecular orbitals and energies, rather than those of a local-DFT calculation, are used for building the Green function and effective Coulomb interaction. The use of an optimally tuned hybrid diminishes the starting point dependency in one-shot stochastic GW, effectively avoiding the need for self-consistent GW iterations.

Chemistry↗
NASA NTRS1991

Laser studies of the photodissociation dynamics of cometary radicals

In the past year, it was shown that in the 193 nm photolysis of C2H, the C2 radical is produced in a variety of electronic, vibrational, and rotational states. The relative population of the vibrational and rotational states of C2(A 1 Pi u), C2(B 1 Sigma g +), and C2(A 3 Pi u) were determined in a static gas cell and in a pulsed molecular beam. It seems as though the original angular momentum of the C2H molecule appears as part of the angular momentum of the C2 radical. A attempt is being made to discover the mathematical relationship that governs this mapping. New information about the bond dissociation energy of the C2 radical was produces. C2(b 3 Sigma g -) and C2( 1 Delta g) were detected in the photolysis of C2H via time resolved infrared emission spectroscopy. In the former case, vibrational excitation up to v'' = 4 is observed. All of the results suggest that the C2 models in comets need to consider the presence of vibrationally excited C2 radicals in comets. The laser induced fluorescence spectra of the C3 was observed as a product of the 193 nm photolysis of allene and propyne. The populations of the rotational levels are identical in both cases. This result has led us to conclude that an isomerization reaction occurs in the photolysis of propyne which leads to the same C3H2 intermediate that is formed in the photolysis of C3H4. Since the former molecule is one of the most abundant in the interstellar medium it is also likely that its precursor is also present in comets. This would explain why C3 is observed in comets.

Jackson, William M.↗
DOE OSTI2025

SULI Intern Final Report: Computationally Investigating Hydrogen Thermo-Diffusion in Yttrium Hydride Using Multiscale Methods

The renaissance of nuclear energy has arrived, heralding an age of abundant inexpensive clean energy, and renewed space exploration. In nuclear-powered spacecraft and microreactors, safety and size are of utmost importance. Yttrium Hydride (YHx) is being researched for its utility as a neutron moderator in nuclear reactors; the hydrogen in YHx slows down neutrons, enabling a continuous nuclear reaction in the reactor. This has the benefit of allowing reactors to be more safe, compact, and efficient. The goal of this effort is to computationally predict the coefficient of temperature-dependent hydrogen diffusion within YHx, the Soret coefficient. This parameter is essential for determining the safe operating modes of YHx moderators. Zirconium Hydride (ZrHx) is used in the Training, Research, Isotopes, General Atomics (TRIGA) reactor, is the reference material for these calculations. In this work, nanoscale atomic modeling in the Vienna Ab initio Simulation Package (VASP) is combined with the mesoscale finite element phase-field module in the Multiphysics Object-Oriented Simulation Environment (MOOSE); this culminates in a new multiscale computational method to simulate Soret diffusion of hydrogen in YHx. This data is useful for predicting experimental outcomes. This workflow involves convergence testing followed by static, Nudged Elastic Band (NEB), Quasi-Harmonic Approximation (QHA), and Molecular Dynamics (MD) calculations - linked with phase field simulation. NEB simulates hydrogen migration, while QHA and MD predict temperature-dependent properties. The static calculations align with literature, and preliminary NEB and QHA calculations yield accurate results. Once the atomic calculations are complete, we will incorporate Electron Backscatter Diffraction (EBSD) images and VASP-generated parameters into the phase field module to simulate intra- and intergranular transport of hydrogen in ZrHx and YHx. Future research will extend our approach to fuel-moderator materials systems such as Uranium-Yttrium Hydride (U-YHx). This work contributes to the development of advanced nuclear energy solutions for space travel.

36 - MATERIALS SCIENCE↗
NASA NTRS1975

Phenomenological treatments of cross-sections for proton and hydrogen impact below 1 keV on molecular nitrogen

An analytic independent-particle model is used to construct static potentials to describe the interaction of hydrogen-like ions with atoms and molecules. Parameters of the ion-atom potential are determined from ab initio total energy minimization procedure. The elastic scattering of He(+) from Ne and Ar is investigated as a test case and comparison is made with experiment. The model is then used in conjunction with low energy H(+)-N2 experimental data to construct differential and total cross-sections for the scattering of protons and hydrogen in the energy range of 10 eV to 1 keV from molecular nitrogen. Analytic forms are used to parametrize these cross-sections to facilitate their use in the calculation of energy deposition by protons and hydrogen atoms in atmospheric gases.

Porter, H. S.↗
NASA NTRS1975

Rarefied gas dynamic effects on mass spectrometric studies of upper planetary atmospheres

Results are presented of measurements obtained with an open-source mass spectrometer both in earth orbit and in a laboratory molecular beam facility. The mass spectrometer/ion source combination was developed in the laboratory to be operable in either of two modes by altering ion extraction potential on ground command, i.e., as a stagnation cavity or to respond only to the incoming molecules of the unperturbed atmosphere. Results indicate that the use of this open-source configuration and dual-mode capability, allied with both standard static calibration and dynamic calibration using high-speed molecular beam techniques in the laboratory, allows collection of useful data on aeronomy of the upper atmosphere of planets and comet tails.

French, J. B.↗
NASA NTRS1981

A molecular model for ice nucleation and growth, attachment 1

The quantum mechanical technique is used to study ionic, configurational, and impurity defects in the ice surface. In addition to static calculations of the energetics of the water monomer-ice surface interactions, molecular dynamics studies were initiated. The calculations of the monomer-ice surface interaction, molecular dynamics studies were initiated. The calculations of monomer-ice surface interactions indicate that many adsorption sites exist on the ice surfaces and that the barriers between bonding sites are relatively low. Bonding on the prism face of ice is preferentially above lattice sites.

Plummer, P. L. M.↗
NASA NTRS1997

Aerodynamics of Stardust Sample Return Capsule

Successful return of interstellar dust and cometary material by the Stardust Sample Return Capsule requires an accurate description of the Earth entry vehicle's aerodynamics. This description must span the hypersonic-rarefied, hypersonic-continuum, supersonic, transonic, and subsonic flow regimes. Data from numerous sources are compiled to accomplish this objective. These include Direct Simulation Monte Carlo analyses, thermochemical nonequilibrium computational fluid dynamics, transonic computational fluid dynamics, existing wind tunnel data, and new wind tunnel data. Four observations are highlighted: 1) a static instability is revealed in the free-molecular and early transitional-flow regime due to aft location of the vehicle s center-of-gravity, 2) the aerodynamics across the hypersonic regime are compared with the Newtonian flow approximation and a correlation between the accuracy of the Newtonian flow assumption and the sonic line position is noted, 3) the primary effect of shape change due to ablation is shown to be a reduction in drag, and 4) a subsonic dynamic instability is revealed which will necessitate either a change in the vehicle s center-of-gravity location or the use of a stabilizing drogue parachute.

Mitcheltree, R. A.↗
DOE OSTI2026

Transition from Vehicular to Structural Ionic Transport in Electrified Alkali Aqueous Solutions

A molecular understanding of the solvation and dynamics of ions under static electric fields is crucial for modelling a wide range of natural and technological processes. Yet, traditional simulation methods suffer from a trade-off that has to be made between accuracy and statistical convergence. To bridge this gap, herein we extend our recently introduced Perturbed Neural Network Potential Molecular Dynamics (PNNP MD) approach to investigate the solvation structures and ionic transport mechanisms of electrified alkali cationic solutions. We obtain ionic conductivities for Li+, Na+ and Cs+ from the field dependence of the ionic current density in good agreement with experiment. Surprisingly, the migration mechanism is found to be strikingly different for the three ions despite their similar ionic conductivities. While Li+ conducts predominantly through vehicular migration of a stable 4-fold coordinated ion at all field strengths, Cs+ conducts strictly through a structural diffusion mechanism, where 9-12 transient first shell water coordination bonds are broken and renewed. Notably, aqueous Na+ emerges as a “Goldilocks” ion: its ion-water interactions are strong enough to maintain distinct 5-6 fold coordination shells at zero field (unlike Cs+) yet labile enough to be strongly perturbed by electric fields (unlike Li+). As a consequence, we observe an electric field-induced transition from vehicular to structural ionic transport for Na+ that is accompanied by a marked increase in the ionic current density. Our results imply that the conductance mechanism of ions with moderate ion-solvent interactions can be effectively tuned by external electric fields.

Joll, Kit↗
DOE OSTI2024

Direct observation of ultrafast cluster dynamics in supercritical carbon dioxide using X-ray Photon Correlation Spectroscopy

Supercritical fluids exhibit distinct thermodynamic and transport properties, making them of particular interest for a wide range of scientific and engineering applications. These anomalous properties emerge from structural heterogeneities due to the formation of molecular clusters at conditions above the critical point. While the static behavior of these clusters and their effects on the thermodynamic response functions have been recognized, the relation between the ultrafast cluster dynamics and transport properties remains elusive. By measuring the intermediate scattering function in carbon dioxide at conditions near the critical point with X-ray photon correlation spectroscopy, we directly capture the cross-over dynamics between 4 and 13 picoseconds, revealing the transition between ballistic and diffusive motion. Complementary analysis using large-scale molecular dynamics simulations reveals that this behavior arises from collisions between unbound molecules and clusters. This study provides direct evidence of the ultrafast momentum exchange between clusters, which has significant impact on transport properties, solvation processes, and reaction kinetics in supercritical fluids.

carbon capture and storage↗
NASA NTRS1983

Calculation of the polarization potential for e-N2 collisions

A polarization potential V(pol) for e-N2 collisions is calculated by the generalization of the static part of the method of polarized orbitals to molecular targets. Partial differential equations (PDE) are derived for polarized orbitals, which are functions of the distance (r) from the molecular center and angle (theta) from the molecular axis. The equations are solved with the use of the noniterative PDE technique. From the polarized orbitals a polarization potential can be constructed whose r and theta dependence is found to be significantly different from the well-known phenomenological one. The two potentials are further compared by carrying out limited hybrid-theory scattering calculations. Only those scattering results based on the calculated V(pol) are in satisfactory accord with experiment.

Onda, K.↗
NASA NTRS1997

Low-Density Aerodynamics of the Stardust Sample Return Capsule

The aerodynamics of the Stardust Sample Return Capsule are analyzed in the low- density, transitional flow regime using free-molecular, Direct Simulation Monte Carlo, Navier-Stokes, and Newtonian methods to provide inputs for constructing a transitional flow bridging relation. The accuracy of this bridging relation in reconstructing the aero- dynamic coefficients given by the more exact methods is presented for a range of flight conditions and vehicle attitudes. There is good agreement between the various prediction methods, and a simple sine-squared bridging relation is shown to provide a reasonably good description of the axial force, normal force, and pitching moment over a range of Knudsen numbers from 0.001 to 10. The predictions show a static instability of the Star- dust capsule in the free-molecular regime that persists well into the transitional flow. The addition of a thin disk to the base of the capsule is shown to remove this static instability. However, the extremely high entry velocity of 12.6 km/s for the proposed trajectory introduces difficult design issues for incorporating this disk caused by the high aerothermal loads that occur even under relatively rarefied conditions.

Wilmoth, Richard G.↗
DOE OSTI2024

Comparing ab initio and quantum-kinetic approaches to electron transport in warm dense matter

Accurate knowledge of the electronic transport properties of warm dense matter is one of the main concerns of research in high-energy-density physics. Three modern approaches with vastly different levels of fidelity are reviewed and compared: the Kubo–Greenwood (KG) approach based on density-functional-theory molecular dynamics simulations (QMD), quantum kinetic theory based on average-atom models, and time-dependent density functional theory. Throughout, emphasis is placed on the connection between static properties of the electrons (e.g., density of states) and transport properties. Overall, it is found that whenever the conduction electrons can be modeled as being nearly free, fair to excellent agreement is found between QMD and kinetic theory approaches. Such a circumstance is required for modeling warm dense matter as a plasma of ions and free electrons, which is assumed in most kinetic theory approaches. The sensitivity of transport properties to the electronic structure is further highlighted by comparing different exchange–correlation approximations in QMD and KG calculations. It is found that the inclusion of exact exchange via thermal hybrid functionals can make a pronounced impact on electrical and thermal conduction in warm dense matter. We also investigate dynamic screening physics via kinetic theory and time-dependent density functional theory calculations of the mean free path of an electron in a hot dense plasma. In sum, we identify three axes along which to make progress in predicting electron transport in warm dense matter.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2008

Aerodynamics for the Mars Phoenix Entry Capsule

Pre-flight aerodynamics data for the Mars Phoenix entry capsule are presented. The aerodynamic coefficients were generated as a function of total angle-of-attack and either Knudsen number, velocity, or Mach number, depending on the flight regime. The database was constructed using continuum flowfield computations and data from the Mars Exploration Rover and Viking programs. Hypersonic and supersonic static coefficients were derived from Navier-Stokes solutions on a pre-flight design trajectory. High-altitude data (free-molecular and transitional regimes) and dynamic pitch damping characteristics were taken from Mars Exploration Rover analysis and testing. Transonic static coefficients from Viking wind tunnel tests were used for capsule aerodynamics under the parachute. Static instabilities were predicted at two points along the reference trajectory and were verified by reconstructed flight data. During the hypersonic instability, the capsule was predicted to trim at angles as high as 2.5 deg with an on-axis center-of-gravity. Trim angles were predicted for off-nominal pitching moment (4.2 deg peak) and a 5 mm off-axis center-ofgravity (4.8 deg peak). Finally, hypersonic static coefficient sensitivities to atmospheric density were predicted to be within uncertainty bounds.

Edquist, Karl T.↗
NASA NTRS2005

Time-Averaged Velocity, Temperature and Density Surveys of Supersonic Free Jets

A spectrally resolved molecular Rayleigh scattering technique was used to simultaneously measure axial component of velocity U, static temperature T, and density p in unheated free jets at Mach numbers M = 0.6,0.95, 1.4 and 1.8. The latter two conditions were achieved using contoured convergent-divergent nozzles. A narrow line-width continuous wave laser was passed through the jet plumes and molecular scattered light from a small region on the beam was collected and analyzed using a Fabry-Perot interferometer. The optical spectrum analysis air density at the probe volume was determined by monitoring the intensity variation of the scattered light using photo-multiplier tubes. The Fabry-Perot interferometer was operated in the imaging mode, whereby the fringe formed at the image plane was captured by a cooled CCD camera. Special attention was given to remove dust particles from the plume and to provide adequate vibration isolation to the optical components. The velocity profiles from various operating conditions were compared with that measured by a Pitot tube. An excellent comparison within 5m's demonstrated the maturity of the technique. Temperature was measured least accurately, within 10K, while density was measured within 1% uncertainty. The survey data consisted of centerline variations and radial profiles of time-averaged U, T and p. The static temperature and density values were used to determine static pressure variations inside the jet. The data provided a comparative study of jet growth rates with increasing Mach number. The current work is part of a data-base development project for Computational Fluid Dynamics and Aeroacoustics codes that endeavor to predict noise characteristics of high speed jets. A limited amount of far field noise spectra from the same jets are also presented. Finally, a direct experimental validation was obtained for the Crocco-Busemann equation which is commonly used to predict temperature and density profiles from known velocity profiles. Data presented in this paper are available in ASCII format upon request.

Panda, Jayanta↗
DOE OSTI2025

Contextual subspace variational quantum eigensolver calculation of the dissociation curve of molecular nitrogen on a superconducting quantum computer

Abstract We present an experimental demonstration of the Contextual Subspace Variational Quantum Eigensolver on superconducting hardware. Calculating the potential energy curve of molecular nitrogen proves challenging for many conventional quantum chemistry techniques, since static correlation dominates in the dissociation limit. Our quantum simulations retain good agreement with the Full Configuration Interaction energy, outperforming all benchmarked single-reference wavefunction techniques in capturing the bond-breaking appropriately. Moreover, our methodology is competitive with multiconfigurational approaches but at a saving of quantum resource, meaning larger active spaces can be treated for a fixed qubit allowance. To achieve this result, we deploy an error mitigation/suppression strategy comprised of Dynamical Decoupling, Measurement-Error Mitigation and Zero-Noise Extrapolation. Circuit parallelization also provides passive noise-averaging and improves the effective shot yield to reduce the measurement overhead. Furthermore, we introduce a modified adaptive ansatz construction algorithm that incorporates hardware awareness into our variational circuits, minimizing the transpilation cost for the target qubit topology.

Physics↗
DOE OSTI2024

Quantitative x-ray scattering of free molecules

Advances in x-ray free electron lasers have made ultrafast scattering a powerful method for investigating molecular reaction kinetics and dynamics. Accurate measurement of the ground-state, static scattering signals of the reacting molecules is pivotal for these pump-probe x-ray scattering experiments as they are the cornerstone for interpreting the observed structural dynamics. Here, this article presents a data calibration procedure, designed for gas-phase x-ray scattering experiments conducted at the Linac Coherent Light Source x-ray Free-Electron Laser at SLAC National Accelerator Laboratory, that makes it possible to derive a quantitative dependence of the scattering signal on the scattering vector. A self-calibration algorithm that optimizes the detector position without reference to a computed pattern is introduced. Angle-of-scattering corrections that account for several small experimental non-idealities are reported. Their implementation leads to near quantitative agreement with theoretical scattering patterns calculated with ab-initio methods as illustrated for two x-ray photon energies and several molecular test systems.

74 ATOMIC AND MOLECULAR PHYSICS↗