SEARCH · Engineering Papers
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.
Results of quantum statistical mechanical calculations and thermodynamic evaluation of the structure of H2CO3 and its stability against dissociation are reported. Under temperature and pressure conditions near the surface of Venus, carbonic acid would predominatly dissociate into H2O and CO2 and, hence, could not contribute to any significant absorption there.
A density-gradient (DG) model is used to calculate quantum-mechanical corrections to classical carrier transport in MOS (Metal Oxide Semiconductor) inversion/accumulation layers. The model is compared to measured data and to a fully self-consistent coupled Schrodinger and Poisson equation (SCSP) solver. Good agreement is demonstrated for MOS capacitors with gate oxide as thin as 21 A. It is then applied to study carrier distribution in ultra short MOSFETs (Metal Oxide Semiconductor Field Effect Transistor) with surface roughness. This work represents the first implementation of the DG formulation on multidimensional unstructured meshes. It was enabled by a powerful scripting approach which provides an easy-to-use and flexible framework for solving the fourth-order PDEs (Partial Differential Equation) of the DG model.
A multidisciplinary, collaborative simulation has been performed on a Grid of geographically distributed PC clusters. The multiscale simulation approach seamlessly combines i) atomistic simulation backed on the molecular dynamics (MD) method and ii) quantum mechanical (QM) calculation based on the density functional theory (DFT), so that accurate but less scalable computations are performed only where they are needed. The multiscale MD/QM simulation code has been Grid-enabled using i) a modular, additive hybridization scheme, ii) multiple QM clustering, and iii) computation/communication overlapping. The Gridified MD/QM simulation code has been used to study environmental effects of water molecules on fracture in silicon. A preliminary run of the code has achieved a parallel efficiency of 94% on 25 PCs distributed over 3 PC clusters in the US and Japan, and a larger test involving 154 processors on 5 distributed PC clusters is in progress.
Collisional quenching of molecular species is an important process in a variety of astrophysical environments including interstellar clouds, photodissociation regions, and cool stellar/planetary atmospheres. In this work, quantum mechanical scattering calculations are presented for the rotational and vibrational relaxation of rotationally-excited CO due to collisions with H, He and H2 for collision energies between 10(exp -6) and approx.15000/cm. The calculations were performed using the close-coupling approach and the l-labeled form of the coupled-states approximation. Cross sections and rate coefficients for the quenching of the v=0-2, j=0-6 levels of CO are presented and comparisons with previous calculations and measurements, where available, are provided.
The fracture of carbon materials is a complex process, the understanding of which is critical to the development of next generation high performance materials. While quantum mechanical (QM) calculations are the most accurate way to model fracture, the fracture behavior of many carbon-based composite engineering materials, such as carbon nanotube (CNT) composites, is a multi-scale process that occurs on time and length scales beyond the practical limitations of QM methods. The Reax Force Field (ReaxFF) is capable of predicting mechanical properties involving strong deformation, bond breaking and bond formation in the classical molecular dynamics framework. This has been achieved by adding to the potential energy function a bond-order term that varies continuously with distance. The use of an empirical bond order potential, such as ReaxFF, enables the simulation of failure in molecular systems that are several orders of magnitude larger than would be possible in QM techniques. In this work, the fracture behavior of an amorphous carbon (AC) matrix reinforced with CNTs was modeled using molecular dynamics with the ReaxFF reactive forcefield. Care was taken to select the appropriate simulation parameters, which can be different from those required when using traditional fixed-bond force fields. The effect of CNT arrangement was investigated with three systems: a single-wall nanotube (SWNT) array, a multi-wall nanotube (MWNT) array, and a SWNT bundle system. For each arrangement, covalent bonds are added between the CNTs and AC, with crosslink fractions ranging from 0-25% of the interfacial CNT atoms. The SWNT and MWNT array systems represent ideal cases with evenly spaced CNTs; the SWNT bundle system represents a more realistic case because, in practice, van der Waals interactions lead to the agglomeration of CNTs into bundles. The simulation results will serve as guidance in setting experimental processing conditions to optimize the mechanical properties of CNT composites.
Here, this work applies a reduced basis method to study the continuum physics of a finite quantum system—either few or many-body. Specifically, I develop reduced-order models, or emulators, for the underlying inhomogeneous Schrödinger equation and train the emulators against the equation's bound-state-like solutions at complex energies. The emulators rapidly and accurately interpolate and extrapolate the matrix elements of the Hamiltonian resolvent operator (Green's function) across a parameter space that includes both complex energy and other real-valued physical inputs in the Schrödinger equation. The spectra, discretized and compressed as the result of emulation, and the associated resolvent matrix elements (or amplitudes), have the defining characteristics of non-Hermitian quantum mechanics calculations, featuring complex eigenenergies with negative imaginary parts and branch cuts moved below the real axis in the complex energy plane. Therefore, one now has a method that extracts continuum physics from bound-state-like calculations and emulates those extractions in the input parameter space. Building on a prior Letter [Zhang, Phys. Rev. Lett. 135, 242501 (2025)], this article provides the full theoretical details, a comprehensive analysis of the method's performance, and a brief discussion of how it can be coupled with existing continuum approaches to perform emulations in their input parameter spaces.
Exact numerical calculations are made for scattering of quantum mechanical particles hitting a square two-dimensional potential barrier (an exact analog of the Goos-Haenchen optical experiments). Quantum mechanical streamlines are plotted and found to be smooth and continuous, to have continuous first derivatives even through the classical forbidden region, and to form quantized vortices around each of the nodal points. A comparison is made between the present numerical calculations and the stationary wave approximation, and good agreement is found between both the Goos-Haenchen shifts and the reflection coefficients. The time-independent Schroedinger equation for real wavefunctions is reduced to solving a nonlinear first-order partial differential equation, leading to a generalization of the Prager-Hirschfelder perturbation scheme. Implications of the hydrodynamical formulation of quantum mechanics are discussed, and cases are cited where quantum and classical mechanical motions are identical.
Quantum mechanical perturbation theory calculation of upper and lower bounds of energy eigenvalues using partitioning methods
The frequency distributions of the first six Lyman lines of hydrogenlike carbon, oxygen, neon, magnesium, aluminum, and silicon ions broadened by the local fields of both ions and electrons are calculated for dense plasmas. The electron collisions are treated by an impact theory allowing (approximately) for level splittings caused by the ion fields, finite duration of the collisions, and screening of the electron fields. These calculations are fully quantum mechanical and include the full Coulomb interaction. Ion effects are calculated in the quasistatic, linear-Stark-effect approximation, using distribution functions of Hooper and Tighe which include correlation and shielding effects. Theoretical uncertainties from the various approximations are estimated, and the scaling of the profiles with density, temperature, and nuclear charge is discussed. A correction for the effects caused by low-frequency field fluctuations is suggested.
Wave functions with embedded screening constants for quantum mechanical ground state calculations on two-electron atom with nuclear charge
State-to-state reaction probabilities are found to be highly final-state specific at state-selected threshold energies for the reactions O + H2 yield OH + H and H + H2 yield H2 + H. The study includes initial rotational states with quantum numbers 0-15, and the specificity is especially dramatic for the more highly rotationally excited reactants. The analysis is based on accurate quantum mechanical reactive scattering calculations. Final-state specificity is shown in general to increase with the rotational quantum number of the reactant diatom, and the trends are confirmed for both zero and nonzero values of the total angular momentum.
A time-independent quantum mechanical approach to the calculation of photodissociation cross sections is developed. The method is based on the use of a discrete variable representation and the application of absorbing boundary conditions. Both total photodissociation cross sections and partial cross sections (hence product state distributions) are obtained using the same basic technique, but the calculation is particularly efficient when only the former quantity is required. The method is applied to the photodissociation of HCl(+) for which accurate potential energy curves and dipole moment functions are available, to the photodissociation of ClCN which is a direct process, involving a single excited electrosonic state, and to the photodissociation of ICN, which involves several strongly coupled excited electronic states.
The interaction of the thermosphere and ionosphere is largely governed by collisions between ions and neutral particles. On Venus and the Earth, O(+) is a dominant ion, and atomic O dominates throughout much of the thermosphere; therefore an accurate O(+)-O cross section is an important prerequisite for understanding the dynamics of planetary upper atmospheres. The cross section and momentum transfer collision frequency are calculated with a quantum mechanical code which includes resonance charge exchange, polarization, and charge-quadrupole effects. Our results agree well with earlier calculations of Stubbe (1968) and Stallcop et al. (1991).
This report provides a description and performance evaluation of the optimization techniques for high performance computing (HPC) implementation of the open source Computational Materials mini-application Aladyn (https://github.com/nasa/aladyn). Aladyn is a basic molecular dynamics code written in FORTRAN 2003, which is designed to demonstrate the use of adaptive neural networks (ANNs) in atomistic simulations. The role of ANNs is to efficiently reproduce the very complex energy landscape resulting from the atomic interactions in materials with the accuracy of the more expensive quantum mechanics-based calculations. The ANN is trained on a large set of atomic structures calculated using the density functional theory (DFT) method. While achieving orders of magnitude faster computational performance than DFT, the ANN-based approach was still very computationally demanding compared to the conventional approach of using empirically fitted energy functions. After its initial development, Aladyn was evaluated and optimized by experts at the NASA Advanced Supercomputing (NAS) division to exploit modern supercomputer architectures. The code has been optimized for execution on multicore central processing units (CPUs), including Intel® Skylake microarchitecture, and on graphic accelerators, such as Nvidia® V100 graphic processing units (GPUs), using Open Multi-Processing (OpenMP) and Open Accelerators (OpenACC) programming interfaces. The optimization achieved a speedup of 4.7 times the baseline version on CPU performance and an additional 2.4 times on CPU+GPU performance. Atomistic computer simulations are a fundamental tool in materials research to model material properties form physics-based first principles. Atomic interaction, governed by Quantum Mechanics (QM) require sophisticated and highly computationally demanding mathematical models to calculate [1]. Classical methods use approximate functional forms, empirically fitted through a set of variable parameters to emulate atomic energies as direct functions of atomic coordinates [2]. While empirical potentials are computationally much simpler, allowing simulations of large-scale systems of up to a trillion (1012) atoms [3], they are substantially less accurate compared to quantum calculations and applicable only to very specific atomic configurations or predefined crystallographic phases. A recently suggested approach is to use heuristic machine learning methods [4], such as those based on Adaptive Neural Networks (ANNs) to predict atomic energies, after being trained on a sufficiently large database of QM-calculated structures [5,6]. This approach reduces significantly the computational complexity, allowing for simulations of orders of magnitude larger systems compared to QM-based methods without compromising accuracy. Still, compared to classical methods using empirical energy functions, ANN methods remain two- to three orders of magnitude more computationally demanding. Hence, the computational cost of simulations, together with the need for extensive training of ANNs, still makes the practical implementation of ANN-based methods quite challenging. The purpose of the Aladyn mini-application software [7], available as open source at https://github.com/nasa/aladyn, is to be a testbed for exploring possible optimization strategies to develop highly scalable parallel algorithms for ANN-based atomistic simulations. Aladyn is aimed at utilizing the architecture of the high-end modern highperformance computing (HPC) hardware based on multicore central processing units (CPUs) equipped with graphic processing unit (GPU) accelerators. Specifically, the goal is to optimize the performance on a single HPC compute node, before implementing scaling to multi-node parallelization using message passing interface (MPI). At the same time, the open source code of Aladyn can serve as a training model for students and professors in academia.
Histidine can be protonated on either or both of the two N atoms of the imidazole moiety. Each of the three possible forms occurs as a result of the stereochemical environment of the histidine side chain. In an atomic model, comparing the possible protonation states in situ, looking at possible hydrogen bonding and metal coordination, it is possible to predict which is most likely to be correct. A more direct method is described that uses quantum-mechanical methods to calculate, also in situ, the minimum geometry and energy for comparison, and therefore to more accurately identify the most likely protonation state.
Abstract A quantum magnetic state due to magnetic charges is never observed, even though they are treated as quantum mechanical variables in theoretical calculations. Here, the occurrence of a novel quantum disordered state of magnetic charges in a nanoengineered magnetic honeycomb lattice of ultra‐small connecting elements is demonstrated. The experimental research, performed using spin resolved neutron scattering, reveals a massively degenerate ground state, comprised of low integer and energetically forbidden high integer magnetic charges, that manifests cooperative paramagnetism at low temperature. The system tends to preserve the degenerate configuration even under large magnetic field application. It exemplifies the robustness of disordered correlation of magnetic charges in a 2D honeycomb lattice. The realization of quantum disordered ground state elucidates the dominance of exchange energy, which is enabled due to the nanoscopic magnetic element size in nanoengineered honeycomb. Consequently, an archetypal platform is envisaged to study quantum mechanical phenomena due to emergent magnetic charges.
A major effort was expended in developing procedures for calculating by quantum-mechanical methods the cross sections for the photodissociation of the hydroxyl radical OH. The molecule OH is the primary dissociation product of water and observations of OH when combined with a determination of the lifetime of OH in the solar radiation field can be used to calculate the rate of release of water from comets. There are many possible photodissociation channels whose efficiency depends upon the details of the spectrum of the radiation. We have calculated the potential energy curves of all the states of OH that can be reached by electric dipole absorption from the ground X(exp 2)(Pi) state by photons with energies up to about 12 eV.
In an earlier paper, Borucki et al (1987) calculated the electrical conductivity and electrical charge on aerosols in Titan's atmosphere due to the ionization by galactic cosmic rays and electron precipitation from Saturn's magnetosphere. The lower atmosphere was predicted to be substantially more conducting than the atmospheres of Earth and Venus because of the high concentration of free electrons. The prediction of a high conductivity is based on the lack of electrophillic species which form negative ions with low mobility and which reduce the number of free electrons. At that time, no molecular species capable of forming negative ions in concentrations sufficient to perturb the atmospheric conductivity were identified. Recently, E. Bakes and her colleagues have been investigating the formation of nitrogenous macromolecules using quantum mechanical methods. Their calculations indicate that the molecules will be highly electrophillic and are likely to be present in the atmosphere at mixing ratios of order 10(exp -7). This mixing ratio is sufficiently large that a substantial reduction in the conductivity is expected at altitudes below 100 km. Revision of the atmospheric model to accommodate the presence of negative ions and to increase the fidelity of the modeling will be described.