Search NASA⌕ Search

SEARCH · Search NASA

Results for “gradient calculation”

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 523 records · Page 29
All sourcesNASA NTRSDOE OSTI
DOE OSTI2014

The multifacet graphically contracted function method. I. Formulation and implementation

The basic formulation for the multifacet generalization of the graphically contracted function (MFGCF) electronic structure method is presented. The analysis includes the discussion of linear dependency and redundancy of the arc factor parameters, the computation of reduced density matrices, Hamiltonian matrix construction, spin-density matrix construction, the computation of optimization gradients for single-state and state-averaged calculations, graphical wave function analysis, and the efficient computation of configuration state function and Slater determinant expansion coefficients. Timings are given for Hamiltonian matrix element and analytic optimization gradient computations for a range of model problems for full-CI Shavitt graphs, and it is observed that both the energy and the gradient computation scale as O(N{sup 2}n{sup 4}) for N electrons and n orbitals. The important arithmetic operations are within dense matrix-matrix product computational kernels, resulting in a computationally efficient procedure. An initial implementation of the method is used to present applications to several challenging chemical systems, including N{sub 2} dissociation, cubic H{sub 8} dissociation, the symmetric dissociation of H{sub 2}O, and the insertion of Be into H{sub 2}. The results are compared to the exact full-CI values and also to those of the previous single-facet GCF expansion form.

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

Materials for Ultra‐Coherent, Mobile, Electron‐Spin Qubits

This research project has had the goal of gaining a better understanding of the physics of electrons bound to the surface of superfluid helium from both experimental and theoretical perspectives. It has particularly been aimed at two areas which had not been well studied: the relaxation and decoherence of the spin of the electrons on the helium surface and how the properties of underlying metallic layers affect the behavior of the electrons when the helium covering the metal is thin. This work is motivated in part by interest in using the spin of these electrons as a quantum bit, or qubit. Low levels of decoherence are advantageous for qubits, and moving the electrons, as one might do in a quantum processor, will be easiest if thin helium films can be employed. It had been suggested that spin decoherence should be very weak for electrons bound to superfluid He, but before this work there have been no quantitative studies of spin relaxation and decoherence. It is especially important to know how moving the electrons across the helium surface would affect their spin coherence. Calculations performed as part of this project show that the Rashba effective magnetic field, the mechanism which limits the spin coherence of mobile electrons in silicon-based devices (an actively pursued qubit technology), is exceptionally weak for electrons bound to helium. This project has identified other decoherence mechanisms which are stronger, but still weak compared to analogous silicon-based structures. Calculated spin coherence times for mobile electrons approach one day, as compared to microseconds in silicon. With coherence times of this magnitude, the spin qubit errors on helium will be completely dominated by errors in the quantum gates. In related work, the possibility of using an artificial spin-orbit interaction (a gradient magnetic field) for quantum operations on the electrons spins was considered. The calculations show that a moderate gradient field, small enough to be generated by a narrow superconducting wire, will enable high-fidelity quantum operations on electrons held in lithographically-defined quantum dots by driving them with a microwave electric field. The spin and motional coherence of the electrons is sufficient to allow high-fidelity 2-qubit quantum operations between electrons in neighboring quantum dots. As an outgrowth of experiments aiming to measure electron spin coherence it was discovered that very high densities of electrons can be stably supported on thin helium films coating ultra-smooth amorphous metallic layers. The measured densities are high enough that the electron system has almost certainly transitioned from an ordered array of electrons, known as a Wigner crystal (ordered by the electrons’ mutual repulsion), to a quantum fluid known as a Fermi liquid. This transition has been a subject of intense interest for over 40 years, since the electron Wigner crystal was first observed with electrons bound to superfluid helium, but it has never been unambiguously observed. Experiments are still underway in these new structures to definitively determine whether true quantum melting of the Wigner crystal has been demonstrated. This work has also catalyzed the development of a new approach for measuring the transport of electrons across very thin helium films, as will be needed for some of the quantum computing applications. The high electron density experiments as well as experiments with electrons bound in quantum dots have led to new techniques which may enable spin coherence measurements.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2016

CORONAL DENSITY STRUCTURE AND ITS ROLE IN WAVE DAMPING IN LOOPS

It has long been established that gradients in the Alfvén speed, and in particular the plasma density, are an essential part of the damping of waves in the magnetically closed solar corona by mechanisms such as resonant absorption and phase mixing. While models of wave damping often assume a fixed density gradient, in this paper the self-consistency of such calculations is assessed by examining the temporal evolution of the coronal density. It is shown conceptually that for some coronal structures, density gradients can evolve in a way that the wave-damping processes are inhibited. For the case of phase mixing we argue that (a) wave heating cannot sustain the assumed density structure and (b) inclusion of feedback of the heating on the density gradient can lead to a highly structured density, although on long timescales. In addition, transport coefficients well in excess of classical are required to maintain the observed coronal density. Hence, the heating of closed coronal structures by global oscillations may face problems arising from the assumption of a fixed density gradient, and the rapid damping of oscillations may have to be accompanied by a separate (non-wave-based) heating mechanism to sustain the required density structuring.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Density functional theory calculations of the thermodynamic and kinetic properties of point defects in β -U

Density functional theory (DFT) calculations of the thermodynamic and kinetic properties of point defects in the β phase of uranium are reported. Defect energies and entropies were calculated using 2 x 2 x 2 supercells and the Generalized Gradient Approximation (GGA) of Perdew-Burke-Ernzerhof (PBE) for the exchange-correlation potential. Due to computational cost, calculations of the vibrational properties governing entropies were performed by only displacing atoms within (roughly) the 3rd or 4th nearest neighbor shell of the defect, which implicitly assumes that atoms beyond this distance are unaffected by the defect. Migration barriers were estimated by nudged elastic band (NEB) calculations. The low symmetry of the β-U phase (the unit cell is tetragonal and contains 30 atoms) results in many point defect configurations and even more migration pathways. A connectivity map, starting from the most stable point defects, was developed in order to identify the rate-limiting step controlling the net diffusion rate in each crystallographic direction. The uranium self-diffusivity tensor was calculated by combining the defect formation energies, entropies, migration barriers and attempt frequencies. The fastest diffusion rate was determined to be a vacancy mechanism in the z crystallographic direction. The predicted uranium self-diffusivity for this mechanism agrees well with available experimental data. Finally, the diffusion mechanisms and rates identified in this study will inform fuel performance models of swelling and gas evolution.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Structural and electronic properties of rare-earth chromites: A computational and experimental study

Here, in this work, the structural, optical, and electronic properties of rare-earth perovskites of the general formula R⁢ CrO 3 , where R⁢ represents the rare-earth Gd, Tb, Dy, Ho, Er, and Tm, have been studied in detail. These compounds were synthesized through a facile citrate route. X-ray diffraction, Raman spectroscopy, and UV-Visible spectroscopy were utilized to reveal the structural evolutions in R⁢ ⁢CrO 3 . The lattice parameters, Cr 3+ –O 2 –Cr 3+ bond angle, and CrO 6 octahedral distortions were found to strongly depend on the ionic radii of R⁢ . First-principles calculations based on density-functional theory within the generalized gradient approximation (GGA) of Perdew-Burke-Ernzerhof (PBE) and strongly constrained and appropriately normed (SCAN) meta-GGA were also employed to calculate the structural and electronic properties of R⁢ ⁢CrO 3 . The ground-state energy, lattice constants, electronic structures, and density of states of R⁢ ⁢CrO 3 were calculated. These provide some insights into the electronic characteristics of the R⁢ ⁢CrO 3 compounds. The calculated values of lattice parameters and band gaps with Hubbard U correction (SCAN+ U ) agree well with values measured experimentally and show more accuracy in predicting the ground-state crystal structure and band structure compared to PBE+ U approximation. The band gap of R⁢ ⁢CrO 3 is found to be independent of the ionic radii of R⁢ from both experiments and calculations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Origin of Unusual Acidity and Li+ Diffusivity in a Series of Water-in-Salt Electrolytes.

Superconcentrated aqueous electrolytes ("water-in-salt" electrolytes, or WiSEs) enable various aqueous battery chemistries beyond the voltage limits imposed by the Pourbaix diagram of water. However, their detailed structural and transport properties remain unexplored and could be better understood through added studies. Here, we report on our observations of strong acidity (pH 2.4) induced by lithium bis(trifluoromethane sulfonyl)imide (LiTFSI) at superconcentration (at 20 mol/kg). Multiple nuclear magnetic resonance (NMR) and pulsed-field gradient (PFG) diffusion NMR experiments, density functional theory (DFT) calculations, and molecular dynamics (MD) simulations reveal that such acidity originates from the formation of nanometric ion-rich structures. The experimental and simulation results indicate the separation of water-rich and ion-rich domains at salt concentrations >= 5 m and the acidity arising therefrom is due to deprotonation of water molecules in the ion-rich domains. As such, the ion-rich domain is composed of hydrophobic -CF3 (of TFSI-) and hydrophilic hydroxyl (OH-) groups. At 20 m concentration, the tortuosity and radius of water diffusion channels are estimated to be similar to 10 and similar to 1 nm, respectively, which are close to values obtained from hydrated Nafion membranes that also have hydrophobic polytetrafluoroethylene (PTFE) backbones and hydrophilic channels consisting of SO3- ion cluster networks providing for the transport of ions and water. Thus, we have discovered the structural similarity between WiSE and hydrated Nafion membranes on the nanometer scale.

Han, Kee Sung↗
DOE OSTI2018

Strong Correlation Effect in Ferrimagnetic Half-Metallic V{sub 2}CoAl and V{sub 2}CoGa Heusler Compounds

We present first-principles calculations based on the full-potential linearized augmented plane wave (FP-LAPW) method to investigate the structural, elastic, electronic, and magnetic properties of the Heusler compounds V{sub 2}CoAl and V{sub 2}CoGa. We explore different magnetic configurations in both regular and inverse Heusler structures: the Cu{sub 2}MnAl-type and Hg{sub 2}CuTi-type structures respectively. Our compounds are found to be ferrimagnets in the inverse Heusler structure. To take into account the effects of the strong correlation among the localized d-states, we perform a calculation using the GGA+U method (the general gradient approximation with the Hubbard on-site Coulomb correction U). Then we show that both compounds exhibit a half-metallic character and have a magnetic moment of 2.00 μ{sub B} per formula unit. This latter statement is coherent with the generalized Slater-Pauling rule (SP). In terms of the band structure calculations, we explain the formation of the indirect energy band gaps in the spin-majority channels. We also use the mean-field approximation (MFA) to discuss the magnetic exchange couplings and predict the Curie temperature for both compounds.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Modeling and simulation of oxygen transport in high burnup LWR fuel

We have developed a formulation for oxygen transport in uranium dioxide nuclear fuel that accounts for the effects of irradiation. The overall simulation combines the evolving isotopic composition, thermochemistry, and oxygen transport in irradiated fuel. The driving forces for oxygen transport are computed from local thermodynamic equilibrium calculations and account for the effects of temperature gradients and composition, including fission products. The proposed method provides a mechanism for including complex thermodynamic models of nuclear fuel in modeling of mass redistribution, and alleviates difficulties associated with the common thermodiffusion formulation. Here, the transport model has been implemented within the nuclear fuel performance code BISON utilizing the thermochemistry code Thermochimica, with burnup calculations provided by the ORIGEN isotopic transmutation code.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2019

Effect of Fe Doping and O Vacancies on the Magnetic Properties of Rutile TiO{sub 2}

Magnetism sources of Fe doping and O vacancies, which coexist in the presence of rutile TiO{sub 2}, are controversial and require a resolution. The effects of Fe doping and point defect on the magnetism of rutile TiO{sub 2} were studied using geometry optimization and energy calculation based on the first principle generalized gradient approximation + U method (GGA+U) of density functional theory. Fe/Vo ratio is 1:1 in rutile TiO{sub 2}. The next closest distance between Fe and Vo is 0.04082 nm. The formation energy of the doping system is the smallest. This doping system has the highest stability. The coexistence of Fe doping and O vacancy achieved ferromagnetic long-range order. The different ratios of Fe/Vo are 1:1, 1:2, 2:1, and 2:2. The magnetic properties of the doping system are significantly different. At a ratio of 2:1, the magnetic moment is the largest among all doping systems. Moreover, the Curie temperature of the doping system can be higher than room temperature. These findings can obtain dilute magnetic semiconductors with potential application value. At a Fe/V{sub O} ratio of 2:2, the Ti{sub 14}Fe{sub 2}O{sub 30} system produces a half-metallic property of 100% electron spin polarizability. This is the most valuable application for designing and preparing novel spintronic injection source dilute magnetic semiconductors. The magnetic source is mainly mediated by the holes generated by Fe doping and O vacancy complexes, thereby causing the spin polarization double exchange effect on the electron of Ti–3d orbital atoms near the O vacancy, O–2p electron orbital, and Fe–3d electron orbital. This feature is consistent with the following theories: the mean field approximation and the double exchange mechanism theories.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2019

Descriptors for Electrolyte-Renormalized Oxidative Stability of Solvents in Lithium-Ion Batteries

Electrolyte stability against oxidation is an important factor limiting the development of high energy density batteries. The highest occupied molecular orbital (HOMO) level of solvent molecules has been traditionally used for understanding trends in their oxidative stability, but this assumes a noninteracting environment. However, solvent HOMO levels are renormalized because of molecules in their solvation shells. In this work, we first demonstrate an inexpensive and accurate method to determine the HOMO level of the solvent followed by simple descriptors for renormalization of the HOMO level due to different electrolyte components. The descriptors are based on Gutmann donor and acceptor numbers of the solvent and other components. The method uses fast generalized gradient approximation-level density functional theory calculations compared to previously used expensive, experimental data-dependent methods. This method can be used to screen for unexplored stable solvents among the large number of known organic compounds to design novel high-voltage stable electrolytes.

25 ENERGY STORAGE↗
DOE OSTI2017

Atomistic non-adiabatic dynamics of the LH2 complex with a GPU-accelerated ab initio exciton model

We recently outlined an efficient multi-tiered parallel ab initio excitonic framework that utilizes time dependent density functional theory (TDDFT) to calculate ground and excited state energies and gradients of large supramolecular complexes in atomistic detail – enabling us to undertake non-adiabatic simulations which explicitly account for the coupled anharmonic vibrational motion of all the constituent atoms in a supramolecular system. Here we apply that framework to the 27 coupled bacterio-chlorophyll-a chromophores which make up the LH2 complex, using it to compute an on-the-fly nonadiabatic surface-hopping (SH) trajectory of electronically excited LH2. Part one of this article is focussed on calibrating our ab initio exciton Hamiltonian using two key parameters: a shift δ, which corrects for the error in TDDFT vertical excitation energies; and an effective dielectric constant ε, which describes the average screening of the transition-dipole coupling between chromophores. Using snapshots obtained from equilibrium molecular dynamics simulations (MD) of LH2, we tune the values of both δ and ε through fitting to the thermally broadened experimental absorption spectrum, giving a linear absorption spectrum that agrees reasonably well with experiment. In part two of this article, we construct a time-resolved picture of the coupled vibrational and excitation energy transfer (EET) dynamics in the sub-picosecond regime following photo-excitation. Assuming Franck–Condon excitation of a narrow eigenstate band centred at 800 nm, we use surface hopping to follow a single nonadiabatic dynamics trajectory within the full eigenstate manifold. Consistent with experimental data, this trajectory gives timescales for B800→B850 population transfer (τ B800→B850 ) between 650–1050 fs, and B800 population decay (τ 800→ ) between 10–50 fs. The dynamical picture that emerges is one of rapidly fluctuating LH2 eigenstates that are delocalized over multiple chromophores and undergo frequent crossing on a femtosecond timescale as a result of the atomic vibrations of the constituent chromophores. The eigenstate fluctuations arise from disorder that is driven by vibrational dynamics with multiple characteristic timescales. The scalability of our ab initio excitonic computational framework across massively parallel architectures opens up the possibility of addressing a wide range of questions, including how specific dynamical motions impact both the pathways and efficiency of electronic energy-transfer within large supramolecular systems.

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

Perturbational treatment of spin-orbit coupling for generally applicable high-level multi-reference methods

An efficient perturbational treatment of spin-orbit coupling within the framework of high-level multi-reference techniques has been implemented in the most recent version of the COLUMBUS quantum chemistry package, extending the existing fully variational two-component (2c) multi-reference configuration interaction singles and doubles (MRCISD) method. The proposed scheme follows related implementations of quasi-degenerate perturbation theory (QDPT) model space techniques. Our model space is built either from uncontracted, large-scale scalar relativistic MRCISD wavefunctions or based on the scalar-relativistic solutions of the linear-response-theory-based multi-configurational averaged quadratic coupled cluster method (LRT-MRAQCC). The latter approach allows for a consistent, approximatively size-consistent and size-extensive treatment of spin-orbit coupling. The approach is described in detail and compared to a number of related techniques. The inherent accuracy of the QDPT approach is validated by comparing cuts of the potential energy surfaces of acrolein and its S, Se, and Te analoga with the corresponding data obtained from matching fully variational spin-orbit MRCISD calculations. The conceptual availability of approximate analytic gradients with respect to geometrical displacements is an attractive feature of the 2c-QDPT-MRCISD and 2c-QDPT-LRT-MRAQCC methods for structure optimization and ab inito molecular dynamics simulations.

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

Analytic first derivatives for a spin-adapted open-shell coupled cluster theory: Evaluation of first-order electrical properties

An analytic scheme is presented for the evaluation of first derivatives of the energy for a unitary group based spin-adapted coupled cluster (CC) theory, namely, the combinatoric open-shell CC (COSCC) approach within the singles and doubles approximation. The widely used Lagrange multiplier approach is employed for the derivation of an analytical expression for the first derivative of the energy, which in combination with the well-established density-matrix formulation, is used for the computation of first-order electrical properties. Derivations of the spin-adapted lambda equations for determining the Lagrange multipliers and the expressions for the spin-free effective density matrices for the COSCC approach are presented. Orbital-relaxation effects due to the electric-field perturbation are treated via the Z-vector technique. We present calculations of the dipole moments for a number of doublet radicals in their ground states using restricted open-shell Hartree-Fock (ROHF) and quasi-restricted HF (QRHF) orbitals in order to demonstrate the applicability of our analytic scheme for computing energy derivatives. We also report calculations of the chlorine electric-field gradients and nuclear quadrupole-coupling constants for the CCl, CH{sub 2}Cl, ClO{sub 2}, and SiCl radicals.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2016

The electronic structure of As{sub 2}S{sub 3}

The electronic band structure and density of states of crystalline As{sub 2}S{sub 3} are calculated using SIESTA code within generalized-gradient approximation (GGA) of density functional theory (DFT). Energy band gap in As{sub 2}S{sub 3} is found to be indirect with a value 1.925 eV. Density of states and partial density of states are in agreement with available experimental data.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

GPU acceleration of local and semilocal density functional calculations in the SPARC electronic structure code

We present a Graphics Processing Unit (GPU)-accelerated version of the real-space SPARC electronic structure code for performing Kohn–Sham density functional theory calculations within the local density and generalized gradient approximations. In particular, we develop a modular math-kernel based implementation for NVIDIA architectures wherein the computationally expensive operations are carried out on the GPUs, with the remainder of the workload retained on the central processing units (CPUs). Here, using representative bulk and slab examples, we show that relative to CPU-only execution, GPUs enable speedups of up to 6× and 60× in node and core hours, respectively, bringing time to solution down to less than 30 s for a metallic system with over 14 000 electrons and enabling significant reductions in computational resources required for a given wall time.

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

FeRh ground state and martensitic transformation

Cubic B2 FeRh exhibits a metamagnetic transition [(111) antiferromagnet (AFM) to ferromagnet (FM)] around 353 K and remains structurally stable at higher temperatures. However, the calculated zero-Kelvin phonons of AFM FeRh exhibit imaginary modes at M points in the Brillouin zone, indicating a premartensitic instability, which is a precursor to a martensitic transformation at low temperatures. Combining electronic-structure calculations with ab initio molecular dynamics, conjugate gradient relaxation, and the solid-state nudged-elastic band methods, we predict that AFM B2 FeRh becomes unstable at ambient pressure and transforms without a barrier to an AFM(111) orthorhombic (martensitic) ground state below 90±10K. In conclusion, we also consider competing structures, in particular, a tetragonal AFM(100) phase that is not the global ground state, as proposed, but a constrained solution.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Observation of Nonlocal Heat Flux Using Thomson Scattering

Here, nonlocal heat flux was measured in laser-produced coronal plasmas using a novel Thomson scattering technique. The measured heat flux was smaller than the classical values inferred from the measured plasma conditions in regions with large temperature gradients and agreed with classical values for weak gradients. Vlasov–Fokker–Planck simulations self consistently calculated the electron distribution functions used to reproduce the measured Thomson scattering spectra and to determine the heat flux. Multigroup nonlocal simulations overestimated the measured heat flux.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗