Search NASA⌕ Search

SEARCH · Search NASA

Results for “energy function”

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 415 records · Page 23
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

Correlation functions from a unified variational principle: Trial Lie groups

Time-dependent expectation values and correlation functions for many-body quantum systems are evaluated by means of a unified variational principle. It optimizes a generating functional depending on sources associated with the observables of interest. It is built by imposing through Lagrange multipliers constraints that account for the initial state (at equilibrium or off equilibrium) and for the backward Heisenberg evolution of the observables. The trial objects are respectively akin to a density operator and to an operator involving the observables of interest and the sources. We work out here the case where trial spaces constitute Lie groups. This choice reduces the original degrees of freedom to those of the underlying Lie algebra, consisting of simple observables; the resulting objects are labeled by the indices of a basis of this algebra. Explicit results are obtained by expanding in powers of the sources. Zeroth and first orders provide thermodynamic quantities and expectation values in the form of mean-field approximations, with dynamical equations having a classical Lie–Poisson structure. At second order, the variational expression for two-time correlation functions separates–as does its exact counterpart–the approximate dynamics of the observables from the approximate correlations in the initial state. Two building blocks are involved: (i) a commutation matrix which stems from the structure constants of the Lie algebra; and (ii) the second-derivative matrix of a free-energy function. The diagonalization of both matrices, required for practical calculations, is worked out, in a way analogous to the standard RPA. The ensuing structure of the variational formulae is the same as for a system of non-interacting bosons (or of harmonic oscillators) plus, at non-zero temperature, classical Gaussian variables. This property is explained by mapping the original Lie algebra onto a simpler Lie algebra. The results, valid for any trial Lie group, fulfill consistency properties and encompass several special cases: linear responses, static and time-dependent fluctuations, zero- and high-temperature limits, static and dynamic stability of small deviations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Orbital relaxation effects on Kohn–Sham frontier orbital energies in density functional theory

We explore effects of orbital relaxation on Kohn–Sham frontier orbital energies in density functional theory by using a nonempirical scaling correction approach developed in Zheng et al. [J. Chem. Phys. 138, 174105 (2013)]. Relaxation of Kohn–Sham orbitals upon addition/removal of a fractional number of electrons to/from a finite system is determined by a systematic perturbative treatment. The information of orbital relaxation is then used to improve the accuracy of predicted Kohn–Sham frontier orbital energies by Hartree–Fock, local density approximation, and generalized gradient approximation methods. The results clearly highlight the significance of capturing the orbital relaxation effects. Moreover, the proposed scaling correction approach provides a useful way of computing derivative gaps and Fukui quantities of N-electron finite systems (N is an integer), without the need to perform self-consistent-field calculations for (N ± 1)-electron systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2013

Bridging single and multireference coupled cluster theories with universal state selective formalism

The universal state selective (USS) multireference approach is used to construct new energy functionals which offers a unique possibility of bridging single and multireference coupled cluster theories (SR/MRCC). These functionals, which can be used to develop iterative and non-iterative approaches, utilize a special form of the trial wavefunctions, which assure additive separability (or size-consistency) of the USS energies in the non-interacting subsystem limit. When the USS formalism is combined with approximate SRCC theories, the resulting formalism can be viewed as a size-consistent version of the method of moments of coupled cluster equations (MMCC) employing a MRCC trial wavefunction. Special cases of the USS formulations, which utilize single reference state specific CC (V.V. Ivanov, D.I. Lyakh, L. Adamowicz, Phys. Chem. Chem. Phys. 11, 2355 (2009)) and tailored CC (T. Kinoshita, O. Hino, R.J. Bartlett, J. Chem. Phys. 123, 074106 (2005)) expansions are also discussed.

coupled cluster theories↗
DOE OSTI2013

NAMD - The Engine for Large-Scale Classical MD Simulations of Biomolecular Systems Based on a Polarizable Force Field (ALCF-2 Early Science Program Technical Report)

Large-scale Molecular dynamics (MD) simulation based on atomic models provide a powerful tool to understand the structure-dynamics-function relationships of important biological systems. Recent advances in computing power, especially in supercomputer, make MD simulation more and more popular in modern scientific community. However, the power of classical MD simulation has been limited mainly by the accuracy of the potential energy function and the efficiency of the dynamic algorithm enabling the adequate configurational sampling. This Early Science Project is aimed to make MD simulation method go beyond current limit using the leadership supercomputer Blue Gene/Q Mira. To address the issue of potential energy accuracy, classical Drude oscillator is developed to take into account the electronic polarizability in molecular systems. In order to develop and test the Drude polarizable force field for biomolecules, extensive MD simulations of typical polypeptides, lipids are needed. This requires efficient dynamic sampling algorithm that allows fast sampling. The sampling issue can be addressed by using advanced strategies based on multiple copies in order to enhance the sampling efficiency of brute-force MD. One such method is called replica-exchange MD (REMD).

Roux, Benoit↗
NASA NTRS1986

Simulation studies for surfaces and materials strength

During this reporting period three investigations were carried out. The first area of research concerned the analysis of the structure-energy relationship in small clusters. This study is very closely related to the improvement of the potential energy functions which are suitable and simple enough to be used in atomistic simulation studies. Parameters obtained from ab initio calculations for dimers and trimers of Al were used to estimate energetics and global minimum energy structures of clusters continuing up to 15 Al atoms. The second research topic addressed modeling of the collision process for atoms impinging on surfaces. In this simulation study qualitative aspects of the O atom collision with a graphite surface were analyzed. Four different O/graphite systems were considered and the aftermath of the impact was analyzed. The final area of investigation was related to the simulation of thin amorphous Si films on crystalline Si substrates. Parameters obtained in an earlier study were used to model an exposed amorphous Si surface and an a-Si/c-Si interface. Structural details for various film thicknesses were investigated at an atomistic level.

Halicioglu, T.↗
NASA NTRS1993

High-Pressure Equation of State for Partially Ionic Solids

Recently, we showed that the cohesive energy of partially ionic solids may be characterized by a two-term energy relationship consisting of a Coulomb term arising from the valence-charge transfer delta Z between the atoms, and a scaled universal energy function E(sup *)(a(sup *)), which accounts for the partially covalent character of the bond and for the repulsion between the atomic cores for small R; a(sup *) is a scaled length. Normalized cohesive-energy curves of alkali halide crystals and of Ti and Ag halide crystals were obtained, and the cohesive-energy-curve parameters were used to generate theoretical equation-of-state (EOS) curves for the Li, Na, K, Cs, and Ag halides. Good agreement was obtained with the experimental isothermal compression curves over a wide pressure range (0-90 kbar). In this paper we verify that the cohesive-energy relationship is valid for divalent partially ionic solids; physically reasonable charge-transfer values (1.80 less than delta Z less than 2.0) are obtained for MgO, CaO, and CaS. Next, EOS curves for LiF, NaF, Nal, CsCl, Csl, MgO, CaO, and CaS are generated in terms of the cohesive-energy parameters. These EOS's yield excellent fits to experimental isothermal-compression data and to shock-wave data to very high pressures (P(sub max)= 250-1350 kbar).

Schlosser, Herbert↗
DOE OSTI2016

A well-scaling natural orbital theory

Here, we introduce an energy functional for ground-state electronic structure calculations. Its variables are the natural spin-orbitals of singlet many-body wave functions and their joint occupation probabilities deriving from controlled approximations to the two-particle density matrix that yield algebraic scaling in general, and Hartree–Fock scaling in its seniority-zero version. Results from the latter version for small molecular systems are compared with those of highly accurate quantum-chemical computations. The energies lie above full configuration interaction calculations, close to doubly occupied configuration interaction calculations. Their accuracy is considerably greater than that obtained from current density-functional theory approximations and from current functionals of the oneparticle density matrix.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2018

Spatial and Spin Symmetry Breaking in Semidefinite-Programming-Based Hartree–Fock Theory

The Hartree–Fock problem was recently recast as a semidefinite optimization over the space of rank-constrained two-body reduced-density matrices (RDMs) [ Phys. Rev. A 2014, 89, 010502(R)]. This formulation of the problem transfers the nonconvexity of the Hartree–Fock energy functional to the rank constraint on the two-body RDM. We consider an equivalent optimization over the space of positive semidefinite one-electron RDMs (1-RDMs) that retains the nonconvexity of the Hartree–Fock energy expression. The optimized 1-RDM satisfies ensemble N-representability conditions, and ensemble spin-state conditions may be imposed as well. The spin-state conditions place additional linear and nonlinear constraints on the 1-RDM. We apply this RDM-based approach to several molecular systems and explore its spatial (point group) and spin (Ŝ 2 and Ŝ 3 ) symmetry breaking properties. When imposing Ŝ 2 and Ŝ 3 symmetry but relaxing point group symmetry, the procedure often locates spatial-symmetry-broken solutions that are difficult to identify using standard algorithms. For example, the RDM-based approach yields a smooth, spatial-symmetry-broken potential energy curve for the well-known Be–H 2 insertion pathway. We also demonstrate numerically that, upon relaxation of Ŝ 2 and Ŝ 3 symmetry constraints, the RDM-based approach is equivalent to real-valued generalized Hartree–Fock theory.

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

Force Field X: A computational microscope to study genetic variation and organic crystals using theory and experiment

Force Field X (FFX) is an open-source software package for atomic resolution modeling of genetic variants and organic crystals that leverages advanced potential energy functions and experimental data. FFX currently consists of nine modular packages with novel algorithms that include global optimization via a many-body expansion, acid–base chemistry using polarizable constant-pH molecular dynamics, estimation of free energy differences, generalized Kirkwood implicit solvent models, and many more. Applications of FFX focus on the use and development of a crystal structure prediction pipeline, biomolecular structure refinement against experimental datasets, and estimation of the thermodynamic effects of genetic variants on both proteins and nucleic acids. The use of Parallel Java and OpenMM combines to offer shared memory, message passing, and graphics processing unit parallelization for high performance simulations. Overall, the FFX platform serves as a computational microscope to study systems ranging from organic crystals to solvated biomolecular systems.

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

On the piecewise convex or concave nature of ground state energy as a function of fractional number of electrons for approximate density functionals

In this paper, we provide a rigorous proof that the Hartree Fock energy, as a function of the fractional electron number, E(N), is piecewise concave. Moreover, for semi-local density functionals, we show that the piecewise convexity of the E(N) curve, as stated in the literature, is not generally true for all fractions. By an analysis based on exchange-only local density approximation and careful examination of the E(N) curve, we find for some systems, there exists a very small concave region, corresponding to adding a small fraction of electrons to the integer system, while the remaining E(N) curve is convex. Several numerical examples are provided as verification. Although the E(N) curve is not convex everywhere in these systems, the previous conclusions on the consequence of the delocalization error in the commonly used density functional approximations, in particular, the underestimation of ionization potential, and the overestimation of electron affinity, and other related issues, remain unchanged. Finally, this suggests that instead of using the term convexity, a modified and more rigorous description for the delocalization error is that the E(N) curve lies below the straight line segment across the neighboring integer points for these approximate functionals.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2026

A new effective theory for stochastic relativistic hydrodynamics

Thermal fluctuations are a fundamental feature of dissipative systems that are essential for understanding physics near the expected critical point of QCD and in small systems. When such fluctuations are modeled naively in relativistic systems, strange features can appear such as negative self-correlation functions. We construct an effective theory for nonlinear stochastic relativistic hydrodynamics that ensure a well-posed mathematical formulation. Using Crooks fluctuation theorem, we derive a symmetry of the effective action that incorporates fluctuations through a suitable free energy functional. For divergence type theories, the action can then be fully specified using a single vector generating current. The equations of motion obtained using this procedure are guaranteed to be flux conservative and symmetric hyperbolic when the dynamics is causal. This ensures that these equations are well-posed (for suitable initial data) and are in a form that can easily be simulated, including with Metropolis techniques.

Mullins, Nicki [University of Illinois at Urbana-C↗
DOE OSTI2020

Wannier–Koopmans method calculations for transition metal oxide band gaps

The widely used density functional theory (DFT) has a major drawback of underestimating the band gaps of materials. Wannier–Koopmans method (WKM) was recently developed for band gap calculations with accuracy on a par with more complicated methods. WKM has been tested for main group covalent semiconductors, alkali halides, 2D materials, and organic crystals. Here we apply the WKM to another interesting type of material system: the transition metal (TM) oxides. TM oxides can be classified as either with d 0 or d 10 closed shell occupancy or partially occupied open shell configuration, and the latter is known to be strongly correlated Mott insulators. We found that, while WKM provides adequate band gaps for the d0 and d10 TM oxides, it fails to provide correct band gaps for the group with partially occupied d states. This issue is also found in other mean-field approaches like the GW calculations. We believe that the problem comes from a strong interaction between the occupied and unoccupied d-state Wannier functions in a partially occupied d-state system. We also found that, for pseudopotential calculations including deep core levels, it is necessary to remove the electron densities of these deep core levels in the Hartree and exchange–correlation energy functional when calculating the WKM correction parameters for the d-state Wannier functions.

36 MATERIALS SCIENCE↗
NASA NTRS1989

Significance of norms and completeness in variational based methods

By means of a simple structural problem, an important requirement often overlooked in practice on the basis functions used in Rayleigh-Ritz-Galerkin type methods is brought into focus. The problem of the static deformation of a uniformly loaded beam is solved variationally by expanding the beam displacement in a Fourier Cosine series. The potential energy functional is rendered stationary subject to the geometric boundary conditions. It is demonstrated that the variational approach does not converge to the true solution. The object is to resolve this paradox, and in so doing, indicate the practical implications of norms and completeness in an appropriate inner product space.

Storch, Joel A.↗
NASA NTRS1993

Global expression for representing cohesive-energy curves. II

Schlosser et al. (1991) showed that the R dependence of the cohesive energy of partially ionic solids may be characterized by a two-term energy relationship consisting of a Coulomb term arising from the charge transfer, delta-Z, and a scaled universal energy function, E*(a *), which accounts for the partially covalent character of the bond and for repulsion between the atomic cores for small R; a* is a scaled length. In the paper by Schlosser et al., the normalized cohesive-energy curves of NaCl-structure alkali-halide crystals were generated with this expression. In this paper we generate the cohesive-energy curves of several families of partially ionic solids with different crystal structures and differing degrees of ionicity. These include the CsCl-structure Cs halides, and the Tl and Ag halides, which have weaker ionic bonding than the alkali halides, and which have the CsCl and NaCl structures, respectively. The cohesive-energy-curve parameters are then used to generate theoretical isothermal compression curves for the Li, Na, K, Cs, and Ag halides. We find good agreement with the available experimental compression data.

Schlosser, Herbert↗
DOE OSTI2015

Impact of geometric, thermal and tunneling effects on nano-transistors

Electronic transistors are fundamental building blocks of large scale integrated circuits in modern advanced electronic equipments, and their sizes have been down-scaled to nanometers. Modeling and simulations in the framework of quantum dynamics have emerged as important tools to study functional characteristics of these nano-devices. This work explores the effects of geometric shapes of semiconductor–insulator interfaces, phonon–electron interactions, and quantum tunneling of three-dimensional (3D) nano-transistors. First, we propose a two-scale energy functional to describe the electron dynamics in a dielectric continuum of device material. Coupled governing equations, i.e., Poisson–Kohn–Sham (PKS) equations, are derived by the variational principle. Additionally, it is found that at a given channel cross section area and gate voltage, the geometry that has the smallest perimeter of the channel cross section offers the largest channel current, which indicates that ultra-thin nanotransistors may not be very efficient in practical applications. Moreover, we introduce a new method to evaluate quantum tunneling effects in nanotransistors without invoking the comparison of classical and quantum predictions. It is found that at a given channel cross section area and gate voltage, the geometry that has the smallest perimeter of the channel cross section has the smallest quantum tunneling ratio, which indicates that geometric defects can lead to higher geometric confinement and larger quantum tunneling effect. Furthermore, although an increase in the phonon–electron interaction strength reduces channel current, it does not have much impact to the quantum tunneling ratio. Finally, advanced numerical techniques, including second order elliptic interface methods, have been applied to ensure computational accuracy and reliability of the present PKS simulation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

Using the electron localization function to correct for confinement physics in semi-local density functional theory

We have previously proposed that further improved functionals for density functional theory can be constructed based on the Armiento-Mattsson subsystem functional scheme if, in addition to the uniform electron gas and surface models used in the Armiento-Mattsson 2005 functional, a model for the strongly confined electron gas is also added. However, of central importance for this scheme is an index that identifies regions in space where the correction provided by the confined electron gas should be applied. The electron localization function (ELF) is a well-known indicator of strongly localized electrons. We use a model of a confined electron gas based on the harmonic oscillator to show that regions with high ELF directly coincide with regions where common exchange energy functionals have large errors. This suggests that the harmonic oscillator model together with an index based on the ELF provides the crucial ingredients for future improved semi-local functionals. For a practical illustration of how the proposed scheme is intended to work for a physical system we discuss monoclinic cupric oxide, CuO. A thorough discussion of this system leads us to promote the cell geometry of CuO as a useful benchmark for future semi-local functionals. Very high ELF values are found in a shell around the O ions, and take its maximum value along the Cu–O directions. An estimate of the exchange functional error from the effect of electron confinement in these regions suggests a magnitude and sign that could account for the error in cell geometry.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1985

Simulation studies for surfaces and materials strength

From intermolecular force studies, it is now known that the overall non-additive contribution to the lattice enegy is positive so that analysis based on only pairwise additivity suggests a shallower intermolecular potential than the true value. Two body contributions alone are also known to be categorically unable to even qualitatively describe some configurations of molecular clusters in the gas phase or the general relaxation and reconstruction of fcc crystal surfaces. In addition, the many-body contribution was shown to play a key role in the stability of certain crystal structures. In these recent analyses, a relatively simple potential energy function (PEF), comprising only a two-body Mie-type potential plus a three-body Axilrod-Teller-type potential, was found to be extremely effective. This same parametric PEF is applied to describe the bulk stability and surface energy for the diamond cubic structure. To test the stability condition, the FCC, BCC, diamond cubic, graphite and beta-tin structures were considered.

Halicioglu, T.↗