Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “quantum mechanical calculations”

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 703 records · Page 39
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Optical measure of disorder: Why Urbach analysis works for amorphous silicon but fails for amorphous carbon

We investigate the basis for the Paradox where the Urbach-slope disorder parameter disagrees with the Raman width in amorphous diamond-like carbon (DLC) materials (lower Urbach-Eo yet greater Raman width). We examined the bandgap and Urbach-slope measurement issues. While significant errors are identified, ultimately these cannot resolve the Paradox. This resolution involves large changes in shape of DLC's absorption band resulting from large energy shifts. To solidify the understanding of band energy shifts, we examined their properties using a-Si:H, well known for its agreement between Raman and Urbach-slope disorder parameters. Using Tauc-Lorentz (T-L) dispersion functions, the Urbach method works for a-Si:H because small changes in shape occur, as measured by (i) peak position, (ii) peak width, (iii) peak amplitude, and (iv) and gap energies. Moreover, for a-Si:H, a 5th condition (v) there is no other feature below the principal band to interfere with its fall. All 5 of these conditions fail for DLCs. Because T-L dispersions accommodate asymmetry, we traced the Paradox to a steepening of the band shape as it pushes toward lower energy. This steepening is unrelated to disorder, but to a boundary value problem where energy cannot be negative. We conclude that Urbach analysis is invalid for DLC materials, but the T-L width is a good measure instead. We also examined another disorder measure, where π bonds should be considered as disorder relative to the σ bonds in an idealized “amorphous diamond” sp 3 network: superconvergence can be used to quantify these π transitions. Finally, further molecular orbital calculations are needed to properly interpret the changes in T-L widths.

36 MATERIALS SCIENCE↗
DOE OSTI2022

A method for generating moving, orthogonal, area preserving polygonal meshes

A new method for generating locally orthogonal polygonal meshes from a set of generator points is presented in which polygon areas are a constraint. The area constraint property is particularly useful for particle methods where moving polygons track a discrete portion of material. Because Voronoi polygon meshes have some very attractive mathematical and numerical properties for numerical computation, a generalization of Voronoi polygon meshes was formulated that enforces a polygon area constraint. Area constrained moving polygonal meshes allow one to develop hybrid particle-mesh numerical methods that display some of the most attractive features of each approach. It is shown that this mesh construction method can continuously reconnect a moving, unstructured polygonal mesh in a pseudo-Lagrangian fashion without change in cell area/volume, and the method's ability to simulate various physical scenarios is shown. Overall, the advantages are identified for incompressible fluid flow calculations, with demonstration cases that include material discontinuities of all three phases of matter and large density jumps.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

On the Finite Nuclear Effect and Gaussian Basis Sets for Four-Component Dirac Hartree-Fock Calculations

In this perspective we describe an overview of Dirac Hartree-Fock theory, which is the foundation of relativistic quantum chemistry methods vital in the molecular modeling of heavy element chemistry and physics. A special emphasis was placed upon the physical connection between the finite nuclear model and the kinetic-balance condition between the large and small components. Additionally, we review the consequences of basis set choices for Dirac Hartree-Fock calculations in the restricted kinetic balance formalism, and show that the basis set limit does not have a hard lower bound due to the nature of the Hamiltonian, referred to in the literature as prolapse. Standard Gaussian basis set contraction schemes can break down in Dirac Hartree-Fock calculations, and special care must be taken for the description of core orbitals near the nucleus. The phenomena described in this perspective are vital for theoreticians and practitioners to understand when selecting a basis set for heavy element quantum chemical calculations using Dirac Hartree-Fock theory.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Modified Fermi’s golden rule rate expressions

Fermi’s golden rule (FGR) serves as the basis for many expressions of spectroscopic observables and quantum transition rates. The utility of FGR has been demonstrated through decades of experimental confirmation. However, there still remain important cases where the evaluation of a FGR rate is ambiguous or ill-defined. Examples are cases where the rate has divergent terms due to the sparsity in the density of final states or time dependent fluctuations of system Hamiltonians. Strictly speaking, assumptions of FGR are no longer valid for such cases. However, it is still possible to define modified FGR rate expressions that are useful as effective rates. The resulting modified FGR rate expressions resolve a long standing ambiguity often encountered in using FGR and offer more reliable ways to model general rate processes. Finally, simple model calculations illustrate the utility and implications of new rate expressions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Stastically determined experimental uncertainties in photon Doppler velocimetry (PDV) measurements

Shock-ramp experiments place high demands on the precision, accuracy, and temporal resolution of transformed velocimetry signals. To minimize uncertainty after total propagation of a shockramp experiment, ideally multiple PDV traces are used. However, each trace has some uncertainty that is a function of the PDV system and analysis of the velocimetry data. In order to determine a reasonable experimental estimate of minimum uncertainty, we calculate observed uncertainties in as-fielded PDV systems using multiple channels to sample the velocity on each shot. We also calculate the uncertainty’s sensitivity to choice of analysis methods.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1961

THEORETICAL STUDY OF ATOMIC AND MOLECULAR GASES AND THEIR REACTIONS IN THE UPPER ATMOSPHERE. SUMMARY OF WORK DONE IN THE PERIOD NOV. 1, 1960 TO APR. 30, 1961

In the previous report it was stated that the major mathematical problems of computations were nearly solved and an attempt is being made to start application of "Tables of Integrals" to the computation of the different states of atoms and molecules. In this respect I like to report that the computation of a number of states of atoms and molecules has been carried out. The results are so far confined only to the equilibrium distance of molecules. We are extending the calculation to other distances so as to compute potential energy curves of different states of molecules. We are also developing Thomas Fermi Statistical model of the atom to verify our results of more accurate calculation. The results so far have been very encouraging and I hope to report some of these results in the next paper.

Gas↗
NASA NTRS1998

Theoretical Studies of Important Processes in Planetary and Comet Atmospheres

Using theoretical quantum chemical calculations, I have successfully described the dissociative recombination (DR) of O2(+) leading to the excited S-1 state of atomic oxygen, the upper state of the well known green line emission. The process is described by O2(+) + e(-) yields O(S-1)) + O(D-1) (1) where e(-) is an electron and the product oxygen atoms are both excited. This process is important in the atmospheres of Venus, Mars and Earth. I have shown in prior calculations that only one repulsive potential curve of O2, f(sup l)Sigma(sub u, sup +) can generate O(S-1) from DR of the lowest vibrational levels of O2(+). However, in the prior results, the calculated quantum yield (i.e. the number of O(S-1) atoms produced for every two product atoms) from the O2(+) v = 0 level was smaller than the laboratory and atmospheric measurements by more than an order of magnitude. Including only direct recombination, the calculated quantum yield for O(S-1) is only 0.0016. In a calculation that accounts for both direct and indirect recombination, the quantum yield is 0.0012. The range of experimentally determined quantum yields is between 0.01 and 0.23. Because of this large difference between the theoretical and experimental quantum yields, it was thought for some time that the ionospheric and laboratory O2(+) must be vibrationally excited since for excited levels, theory gave quantum yields that are similar to experimental yields. It was also suggested that some other process was generating O(S-1) but none could be identified. Under current NASA support, I have found that reaction (1) proceeds via an unusual mechanism. The f(sup 1)Sigma(sub u, sup,+) state does not cross the ion between the turning points of the v = 0 level of the O2(+) ground state. The lack of a favorable crossing leads to a very small calculated DR rate coefficient. However, this mechanism assumes that initial electron capture must occur into the repulsive state that leads to )(S-1). Instead, I have found that initial electron capture occurs mostly into the B(sup 3)Sigma(sup u, sub -) state which crosses the ion between the turning points of the v = 0 ion level and has a large DR rate coefficient. The B state dissociates to O(D-1) and )(P-3). After capture, some of the flux is transferred to the f(sup 1)Sigma(sup u, sub +) state via symmetry mixed intermediate Rydberg states. The neutral Rydberg states are a mixture of Sigma(sub u, sup +)-1 and Sigma(sub u, sup -)-3 symmetry.

Guberman, Steven L.↗
DOE OSTI2024

Interplay between Topological States and Rashba States as Manifested on Surface Steps at Room Temperature

The unique spin texture of quantum states in topological materials underpins many proposed spintronic applications. However, realizations of such great potential are stymied by perturbations, such as temperature and local fields imposed by impurities and defects, that can render a promising quantum state uncontrollable. Here, in this study, we report room-temperature scanning tunneling microscopy/spectroscopy observation of interaction between Rashba states and topological surface states, which manifests local electronic structure along step edges controllable by the layer thickness of thin films. The first-principles theoretical calculation elucidates the robust Rashba states coexisting with topological surface states along the surface steps with characteristic spin textures in momentum space. Furthermore, the Rashba edge states can be switched off by reducing the thickness of a topological insulator Bi 2 Se 3 to bolster their interaction with the hybridized topological surface states. The study unveils a manipulating mechanism of the spin textures at room temperature, reinforcing the necessity of thin film technology in controlling the quantum states.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Interplays of electron and nuclear motions along CO dissociation trajectory in myoglobin revealed by ultrafast X-rays and quantum dynamics calculations

Significance In this multidisciplinary study, we investigated local and global structural responses upon carbon monoxide (CO) dissociation from the heme iron(II) in myoglobin in solution using ultrafast X-ray absorption and scattering. CO dissociation was optically triggered and pulsed X-rays probed structural and electronic changes as a function of time, revealing an iron(II) spin state transition and nuclear motions propagating from the active site throughout the entire protein. Quantum mechanical computational results were used to assign the iron K-edge absorption features, and to map out electronic and nuclear trajectories during Fe–CO bond elongation and heme macrocycle doming. The study has important implications for heme–protein function and the interaction of small molecule substrates with the metal active centers in enzymes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Ab-initio simulations of coherent phonon-induced pumping of carriers in zirconium pentatelluride

Abstract Laser-driven coherent phonons can act as modulated strain fields and modify the adiabatic ground state topology of quantum materials. Here we use time-dependent first-principles and effective model calculations to simulate the effect of the coherent phonon induced by strong terahertz electric field on electronic carriers in the topological insulator ZrTe 5 . We show that a coherent A 1g Raman mode modulation can effectively pump carriers across the band gap, even though the phonon energy is about an order of magnitude smaller than the equilibrium band gap. We reveal the microscopic mechanism of this effect which occurs via Landau–Zener-Stückelberg tunneling of Bloch electrons in a narrow region in the Brillouin zone center where the transient energy gap closes when the system switches from strong to weak topological insulator. The quantum dynamics simulation results are in excellent agreement with recent pump-probe experiments in ZrTe 5 at low temperature.

36 MATERIALS SCIENCE↗
DOE OSTI2021

OrbNet Denali: A machine learning potential for biological and organic chemistry with semi-empirical cost and DFT accuracy

We present OrbNet Denali, a machine learning model for an electronic structure that is designed as a drop-in replacement for groundstate density functional theory (DFT) energy calculations. The model is a message-passing graph neural network that uses symmetryadapted atomic orbital features from a low-cost quantum calculation to predict the energy of a molecule. OrbNet Denali is trained on a vast dataset of 2.3 × 10 6 DFT calculations on molecules and geometries. This dataset covers the most common elements in biochemistry and organic chemistry (H, Li, B, C, N, O, F, Na, Mg, Si, P, S, Cl, K, Ca, Br, and I) and charged molecules. OrbNet Denali is demonstrated on several well-established benchmark datasets, and we find that it provides accuracy that is on par with modern DFT methods while offering a speedup of up to three orders of magnitude. For the GMTKN55 benchmark set, OrbNet Denali achieves WTMAD-1 and WTMAD-2 scores of 7.19 and 9.84, on par with modern DFT functionals. For several GMTKN55 subsets, which contain chemical problems that are not present in the training set, OrbNet Denali produces a mean absolute error comparable to those of DFT methods. For the Hutchison conformer benchmark set, OrbNet Denali has a median correlation coefficient of R 2 = 0.90 compared to the reference DLPNO-CCSD(T) calculation and R 2 = 0.97 compared to the method used to generate the training data (ωB97X-D3/def2- TZVP), exceeding the performance of any other method with a similar cost. Similarly, the model reaches chemical accuracy for noncovalent interactions in the S66x10 dataset. For torsional profiles, OrbNet Denali reproduces the torsion profiles of ωB97X-D3/def2-TZVP with an average mean absolute error of 0.12 kcal/mol for the potential energy surfaces of the diverse fragments in the TorsionNet500 dataset.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Role of generic scale invariance in a Mott transition from a U(1) spin-liquid insulator to a Landau Fermi-liquid metal

Highlights: • We look at the role played by generic scale invariance (GSI) in the phase transition between a U(1) spin liquid and a Fermi liquid. • We identify dangerously irrelevant operators in the renormalization group transformation, involved with the GSI. • We highlight novel mean-field dynamics of this Mott quantum criticality, resulting from the GSI. We investigate the role of generic scale invariance in a Mott transition from a U(1) spin-liquid insulator to a Landau Fermi-liquid metal, where there exist massless degrees of freedom in addition to quantum critical fluctuations. Here, the Mott quantum criticality is described by critical charge fluctuations, and additional gapless excitations are U(1) gauge-field fluctuations coupled to a spinon Fermi surface in the spin-liquid state, which turn out to play a central role in the Mott transition. An interesting feature of this problem is that the scaling dimension of effective leading local interactions between critical charge fluctuations differs from that of the coupling constant between U(1) gauge fields and matter-field fluctuations in the presence of a Fermi surface. As a result, there appear dangerously irrelevant operators, which can cause conceptual difficulty in the implementation of renormalization group (RG) transformations. Indeed, we find that the curvature term along the angular direction of the spinon Fermi surface is dangerously irrelevant at this spin-liquid Mott quantum criticality, responsible for divergence of the self-energy correction term in U(1) gauge-field fluctuations. Performing the RG analysis in the one-loop level based on the dimensional regularization method, we reveal that such extremely overdamped dynamics of U(1) gauge-field fluctuations, which originates from the emergent one-dimensional dynamics of spinons, does not cause any renormalization effects to the effective dynamics of both critical charge fluctuations and spinon excitations. However, it turns out that the coupling between U(1) gauge-field fluctuations and both matter-field excitations still persists at this Mott transition, which results in novel mean-field dynamics to explain the nature of the spin-liquid Mott quantum criticality. We discuss physical implications of effective one-dimensional spin dynamics and extremely overdamped gauge dynamics at the Mott quantum criticality.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Neural Network for Principle of Least Action

The principle of least action is the cornerstone of classical mechanics, theory of relativity, quantum mechanics, and thermodynamics. Here, we describe how a neural network (NN) learns to find the trajectory for a Lennard-Jones (LJ) system that maintains balance in minimizing the Onsager–Machlup (OM) action and maintaining the energy conservation. The phase-space trajectory thus calculated is in excellent agreement with the corresponding results from the “ground-truth” molecular dynamics (MD) simulation. Furthermore, we show that the NN can easily find structural transformation pathways for LJ clusters, for example, the basin-hopping transformation of an LJ 38 from an incomplete Mackay icosahedron to a truncated face-centered cubic octahedron. Unlike MD, the NN computes atomic trajectories over the entire temporal domain in one fell swoop, and the NN time step is a factor of 20 larger than the MD time step. The NN approach to OM action is quite general and can be adapted to model morphometrics in a variety of applications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Automatic Generation of Chemical Mechanisms for Electrochemical Systems: Solid Electrolyte Interphase Formation in Lithium Batteries

Electrolytes in many lithium ion batteries decompose at the low potentials near the anode. The decomposition products form a layer termed the solid electrolyte interphase (SEI). The composition and growth of the SEI layer significantly affect both the capacity fade and safety of lithium ion batteries. However, SEI formation and growth kinetics are not well understood. In this work, we present an extension of the Reaction Mechanism Generator (RMG) software to automatically generate mechanisms for SEI formation. We extend RMG’s solvation correction framework to account for kinetic solvent effects and demonstrate the accuracy of our technique. We calculate thermochemical parameters for 252 species and rate coefficients for 69 reactions, most with associated solvation corrections. This and additional quantum chemistry data are used to extend RMG’s thermodynamic group additivity and solute parameter estimation schemes to handle lithiated species and add 14 new reaction families to RMG. RMG is additionally extended to simulate electrocatalytic systems. Lastly, we demonstrate RMG on the decomposition of acetonitrile and ethylene carbonate near a battery anode. Furthermore, while this framework does not yet resolve individual ions, as appropriate thermochemistry estimators are not available, and thus, cannot yet resolve more complex electrochemical pathways, RMG is able to generate reasonable pathways for SEI formation that match literature pathways and products. In particular, RMG identifies a new important reaction pathway that is not present in literature.

Chemical reactions↗
DOE OSTI2020

Integrating Electrical Resistivity Tomography with Vadose Zone Characterization - 20246

The environmental management mission at the U.S. Department of Energy Hanford Site includes addressing subsurface contaminants that are present in the vadose zone because of historical waste discharges. Characterization is needed to support evaluation of whether contaminants at specific disposal sites will migrate to and contaminate the groundwater in the future. Drilling of boreholes and analysis of retrieved sediment samples is a standard approach to evaluating contaminant distribution in the vadose zone for the types of inorganic and radionuclide contaminants at the Hanford Site. However, borehole sampling has limitations for understanding the volumetric distribution of contaminants. The Hanford Site is now applying electrical resistivity tomography to characterize the path and volumetric extent of contaminant migration in the vadose zone. This information aids selection of appropriate locations for obtaining borehole samples and provides three-dimensional contaminant distribution information for incorporation into conceptual site models and to help configure fate and transport models. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2021

MechElastic: A Python library for analysis of mechanical and elastic properties of bulk and 2D materials

We report the MechElastic Python package evaluates the mechanical and elastic properties of bulk and 2D materials using the elastic coefficient matrix ( C ij ) obtained from any ab-initio density-functional theory (DFT) code. The current version of this package reads the output of VASP, ABINIT, and Quantum Espresso codes (but it can be easily generalized to any other DFT code) and performs the appropriate post-processing of elastic constants as per the requirement of the user. This program can also detect the input structure's crystal symmetry and test the mechanical stability of all crystal classes using the Born-Huang criteria. Various useful material-specific properties such as elastic moduli, longitudinal and transverse elastic wave velocities, Debye temperature, elastic anisotropy, 2D layer modulus, hardness, Pugh's ratio, Cauchy's pressure, Kleinman's parameter, and Lame's coefficients, can be estimated using this program. Another existing feature of this program is to employ the ELATE package (2016) [29] and plot the spatial variation of several elastic properties such as Poisson's ratio, linear compressibility, shear modulus, and Young's modulus in three dimensions. Further, the MechElastic package can plot the equation of state (EOS) curves for energy and pressure for a variety of EOS models such as Murnaghan, Birch, Birch-Murnaghan, and Vinet, by reading the inputted energy/pressure versus volume data obtained via numerical calculations or experiments. This package is particularly useful for the high-throughput analysis of elastic and mechanical properties of materials.

2D materials↗
DOE OSTI2026

Field-tailoring quantum materials via magneto-synthesis: metastable metallic and magnetically suppressed phases in a trimer iridate

We demonstrate that applying modest magnetic fields (< 0.1 T) during high-temperature crystal growth can profoundly alter the structure and ground state of a spin-orbit-coupled, antiferromagnetic trimer lattice. Using BaIrO₃ as a model system, whose ground state is intricately dictated by the trimer lattice, we show that magneto-synthesis , a field-assisted synthesis approach, stabilizes a structurally compressed, metastable metallic and magnetically suppressed phases inaccessible via conventional methods. These effects include a 0.85% reduction in unit cell, 4-order-of-magnitude decrease in resistivity, a 10-fold enhancement of the Sommerfeld coefficient, and the collapse of long-range magnetic order -- all intrinsic and bulk in origin. First-principles calculations confirm that the field-stabilized structure lies substantially above the ground state in energy, highlighting its metastable character. These large, coherent and correlated changes across multiple bulk properties, unlike those caused by dilute impurities, defects or off-stoichiometry, point to an intrinsic field-induced mechanism. The findings establish magneto-synthesis as a powerful new pathway for accessing non-equilibrium quantum phases in strongly correlated materials.

magneto-synthesis↗