Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “electronic structure methods”

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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Hybrid gausslet/Gaussian basis sets

We introduce hybrid gausslet/Gaussian basis sets, where a standard Gaussian basis is added to a gausslet basis in order to increase accuracy near the nuclei while keeping the spacing of the grid of gausslets relatively large. The Gaussians are orthogonalized to the gausslets, which are already orthonormal, and approximations are introduced to maintain the diagonal property of the two electron part of the Hamiltonian so that it continues to scale as the second power of the number of basis functions rather than the fourth. We introduce several corrections to the Hamiltonian designed to enforce certain exact properties, such as the values of certain two-electron integrals. We also introduce a simple universal energy correction that compensates for the incompleteness of the basis stemming from the electron–electron cusps based on the measured double occupancy of each basis function. We perform a number of Hartree Fock and full configuration interaction (full-CI) test calculations on two electron systems and Hartree Fock on a ten-atom hydrogen chain to benchmark these techniques. The inclusion of the cusp correction allows us to obtain complete basis set full-CI results for the two electron cases at the level of several micro-Hartrees, and we see similar apparent accuracy for Hartree Fock on the ten-atom hydrogen chain.

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

Excitons: Energetics and spatiotemporal dynamics

The concept of an exciton as a quasiparticle that represents collective excited states was originally adapted from solid-state physics and has been successfully applied to molecular aggregates by relying on the well-established limits of the Wannier exciton and the Frenkel exciton. However, the study of excitons in more complex chemical systems and solid materials over the past two decades has made it clear that simple concepts based on Wannier or Frenkel excitons are not sufficient to describe detailed excitonic behavior, especially in nano-structured solid materials, multichromophoric macromolecules, and complex molecular aggregates. In addition, important effects such as vibronic coupling, the influence of charge-transfer (CT) components, spin-state interconversion, and electronic correlation, which had long been studied but not fully understood, have turned out to play a central role in many systems. This has motivated new experimental approaches and theoretical studies of increasing sophistication. This article provides an overview of works addressing these issues that were published for A Special Topic of the Journal of Chemical Physics on ``Excitons: Energetics and spatio-temporal dynamics" and discusses their implications.

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

Real-time equation-of-motion CC cumulant and CC Green’s function simulations of photoemission spectra of water and water dimer

Newly developed coupled-cluster (CC) methods enable simulations of ionization potentials and spectral functions of molecular systems in a wide range of energy scales ranging from core-binding to valence. This paper discusses the results obtained with the real-time equation-of-motion CC cumulant (RT-EOM-CC) approach and CC Green’s function (CCGF) approaches in applications to the water and water dimer molecules. Here, we compare the ionization potentials obtained with these methods for the valence region with the results obtained with the coupled-cluster with singles, doubles, and perturbative triples formulation as a difference of energies for N and N – 1 electron systems. All methods show good agreement with each other. They also agree well with the experiment with errors usually below 0.1 eV for the ionization potentials. We also analyze unique features of the spectral functions, associated with the position of satellite peaks, obtained with the RT-EOM-CC and CCGF methods employing single and double excitations, as a function of the monomer OH bond length and the proton transfer coordinate in the dimer. Finally, we analyze the impact of the basis set effects on the quality of calculated ionization potentials and find that the basis set effects are less pronounced for the augmented-type sets.

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

Mass-selected ion–molecule cluster beam apparatus for ultrafast photofragmentation studies

We describe an apparatus for investigating the excited-state dissociation dynamics of mass-selected ion–molecule clusters by mass-resolving and detecting photofragment-ions and neutrals, in coincidence, using an ultrafast laser operating at high repetition rates. The apparatus comprises a source that generates ion–molecule clusters, a time-of-flight spectrometer, and a mass filter that selects the desired anions, and a linear-plus-quadratic reflectron mass spectrometer that discriminates the fragment anions after the femtosecond laser excites the clusters. The fragment neutrals and anions are then captured by two channeltron detectors. The apparatus performance is tested by measuring the photofragments: I – , CF 3 I – , and neutrals from photoexcitation of the ion–molecule cluster CF 3 I·I – using femtosecond UV laser pulses with a wavelength of 266 nm. Here, the experimental results are compared with our ground state and excited state electronic structure calculations as well as the existing results and calculations, with particular attention to the generation mechanism of the anion fragments and dissociation channels of the ion–molecule cluster CF 3 I·I – in the charge-transfer excited state.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2024

Hybrid programming-model strategies for GPU offloading of electronic structure calculation kernels

To address the challenge of performance portability and facilitate the implementation of electronic structure solvers, we developed the basic matrix library (BML) and Parallel, Rapid O(N), and Graph-based Recursive Electronic Structure Solver (PROGRESS) library. The BML implements linear algebra operations necessary for electronic structure kernels using a unified user interface for various matrix formats (dense and sparse) and architectures (CPUs and GPUs). Focusing on density functional theory and tight-binding models, PROGRESS implements several solvers for computing the single-particle density matrix and relies on BML. In this paper, we describe the general strategies used for these implementations on various computer architectures, using OpenMP target functionalities on GPUs, in conjunction with third-party libraries to handle performance critical numerical kernels. In this study, we demonstrate the portability of this approach and its performance in benchmark problems.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Exploring direct photodetachment and photodissociation–photodetachment dynamics of platinum iodide anions (PtI n - , n = 2–5) using cryogenic photoelectron spectroscopy

The direct photodetachment and two-photon photodissociation–photodetachment processes of a series of PtI n - (n = 2–5) anions were systematically studied using cryogenic anion photoelectron spectroscopy and first-principles electronic structure calculations. The adiabatic/vertical detachment energies (ADEs/VDEs) of these anions were determined from their 193 nm photoelectron (PE) spectra, i.e., 3.54/3.63, 4.04/4.09, 4.33/4.36, and 4.37/4.41 eV for n = 2–5, respectively, and well reproduced by B3LYP-D3(BJ)/aug-cc-pVTZ-pp calculations. As the coordination number increases, the electron affinity (EA) of PtI n • (n = 2–5) neutrals (equivalent to the corresponding anion’s ADE) gradually increases, exceeding the EA of Cl at n = 3 and exhibiting superhalogen characteristics for n ≥ 3. Meanwhile, the ground state transition contributed from detaching electrons in the highest occupied molecular orbital gradually evolves from the central metal Pt to the iodine ligands. For the PtI 3 - anion, besides one-photon direct detachment, four distinct two-photon photodissociation–photodetachment channels were identified, and the competition between them was discussed.

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

Collision integrals within the Chapman–Enskog theory for a generalized Lennard-Jones potential

Here we report the values of the collision integrals, needed for the calculation of the macroscopic transport properties such as viscosity (η) and diffusion coefficient (D) of gases within the Chapman–Enskog kinetic gas theory, for a generalized Lennard-Jones potential (gLJ), a more general potential with an adjustable long range 1/r dependence that can describe a wide range of intermolecular interactions.

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

Symmetry breaking as predicted by a phase space Hamiltonian with a spin Coriolis potential

Here, we perform electronic structure calculations for a set of molecules with degenerate spin-dependent ground states ( 3 CH 2 , 2 CH$^{•}_{3}$, 3 O 2 ) going beyond the Born–Oppenheimer approximation and accounting for nuclear motion. According to a phase space approach that parameterizes electronic states (|Φ⟩) and electronic energies (E) by nuclear position and momentum [i.e., |Φ(R, P)⟩ and E(R, P)], we find that the presence of degenerate spin degrees of freedom leads to broken symmetry ground states. More precisely, rather than a single degenerate minimum at (R, P) = (R min , 0), the ground state energy has two minima at (R,P)=(R' min ,±P min ) (where R' min is close to R min ), dramatically contradicting the notion that the total energy of the system can be written in separable form as E = $\frac{P^2}{2M}$ + V el . Although we find that the broken symmetry solutions have small barriers between them for the small molecules, we hypothesize that the barriers should be macroscopically large for metallic solids, thus offering up a new phase-space potential energy surface for simulating the Einstein–de Haas effect.

Berry connection↗
DOE OSTI2021

Variational quantum eigensolver for approximate diagonalization of downfolded Hamiltonians using generalized unitary coupled cluster ansatz

In this paper we discuss the utilization of Variational Quantum Solver (VQE) and recently introduced Generalized Unitary Coupled Cluster (GUCC) formalism for the diagonalization of downfolded/effective Hamiltonians in active spaces. In addition to effective Hamiltonians defined by the downfolding of a subset of virtual orbitals we also consider their form defined by freezing core orbitals, which enables us to deal with larger systems. Furthermore, we also consider various solvers to identify solutions of the GUCC equations. We use N 2 , H 2 O, and C 2 H 4 , and benchmark systems to illustrate the performance of the combined framework.

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

Correlated purification for restoring 𝑁-representability in quantum simulation

Experimentally measured reduced density matrices (RDMs) often violate constraints that ensure they represent N-electron states—known as N-representability conditions—because of statistical and hardware noise. In this work, we present a correlated purification framework based on semidefinite programming to restore the accuracy of a noisy, unphysical two-electron RDM (2-RDM). The method performs a bi-objective optimization that minimizes both the many-electron energy and the nuclear norm of the correction to the measured 2-RDM. The nuclear norm, often employed in matrix completion, promotes low-rank corrections, while the energy term acts as a regularization term that can improve the purity of the ground state. While the method is particularly effective for ground states, it can also be applied to excited and nonstationary states by decreasing the weight of the energy relative to the error norm. In an application to fermionic shadow tomography of large hydrogen chains, correlated purification yields substantial reductions in both energy and 2-RDM error, achieving chemical accuracy across dissociation curves. This framework provides a robust strategy for tomography in many-body quantum simulations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Gauge invariance of light-matter interactions in first-principle tight-binding models

In this work, we study the different ways of introducing light-matter interaction in first-principle tight-binding (TB) models. The standard way of describing optical properties is the velocity gauge, defined by linear coupling to the vector potential. In finite systems a transformation to represent the electromagnetic radiation by the electric field instead is possible, albeit subtleties arise in periodic systems. The resulting dipole gauge is a multi-orbital generalization of the Peierls substitution. In this work we investigate the accuracy for both pathways, with particular emphasis on gauge invariance, for TB models constructed from maximally localized Wannier functions. This approach accurately captures the light-matter interaction close to the Fermi level. Focusing on paradigmatic two-dimensional materials, we construct first-principle models and calculate the response to electromagnetic fields in linear response and for strong excitations. Benchmarks against fully converged first-principle calculations allow for ascertaining the accuracy of the TB models. We find that the dipole gauge provides a more accurate description than the velocity gauge in all cases. The main deficiency of the velocity gauge is an imperfect cancellation of paramagnetic and diamagnetic current. Formulating a corresponding sum rule however provides a way to explicitly enforce this cancellation. This procedure corrects the TB models in the velocity gauge, yielding excellent agreement with dipole gauge and thus restoring gauge invariance.

2-dimensional systems↗
DOE OSTI2021

Twisted symmetric trilayer graphene. II. Projected Hartree-Fock study

The Hamiltonian of the magic-angle twisted symmetric trilayer graphene (TSTG) can be decomposed into a twisted-bilayer-graphene- (TBG-) like flat band Hamiltonian and a high-velocity Dirac fermion Hamiltonian. We use Hartree-Fock mean field approach to study the projected Coulomb interacting Hamiltonian of TSTG developed in Calugaru et al. [Phys. Rev. B 103, 195411 (2021)] at integer fillings ν = -3, -2, -1, and 0 measured from charge neutrality. We study the phase diagram with w 0 /w 1 , the ratio of AA and AB interlayer hoppings, and the displacement field, which introduces an interlayer potential U and hybridizes the TBG-like bands with the Dirac bands. At small U, we find the ground states at all fillings ν are in the same phases as the tensor products of a Dirac semimetal with the filling ν TBG insulator ground states, which are spin-valley polarized at ν = -3, and fully (partially) intervalley coherent at ν = -2, 0(ν = -1) in the flat bands. An exception is ν = -3 with w 0 /w 1 ≳ 0.7 , which possibly becomes a metal with competing orders at small U due to charge transfers between the Dirac and flat bands. At strong U where the bandwidths exceed interactions, all the fillings ν enter a metal phase with small or zero valley polarization and intervalley coherence. Lastly, at intermediate U, semimetal or insulator phases with zero intervalley coherence may arise for ν = -2, -1, 0. Finally, our results provide a simple picture for the electron interactions in TSTG systems, and reveal the connection between the TSTG and TBG ground states.

2-dimensional systems↗
DOE OSTI2022

Interaction-expansion inchworm Monte Carlo solver for lattice and impurity models

Multiorbital quantum impurity models with general interaction and hybridization terms appear in a wide range of applications, including embedding, quantum transport, and nanoscience. However, most quantum impurity solvers are restricted to a few impurity orbitals, discretized baths, diagonal hybridizations, or density-density interactions. We generalize the inchworm quantum Monte Carlo method to the interaction expansion, and we explore its application to typical single- and multiorbital problems encountered in investigations of impurity and lattice models. Our implementation generically outperforms bare and bold-line quantum Monte Carlo algorithms in the interaction expansion. For the systems studied here, our implementation remains inferior to the more specialized hybridization expansion and auxiliary field algorithms. So, the problem of convergence to unphysical fixed points, which hampers so-called bold-line methods, is not encountered in inchworm Monte Carlo.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Adaptive time stepping for the two-time integro-differential Kadanoff-Baym equations

The nonequilibrium Green's function gives access to one-body observables for quantum systems. Of particular interest are quantities such as density, currents, and absorption spectra which are important for interpreting experimental results in quantum transport and spectroscopy. We present an integration scheme for the Green's function's equations of motion, the Kadanoff-Baym equations (KBE), which is both adaptive in the time integrator step size and method order as well as the history integration order. We analyze the importance of solving the KBE self-consistently and show that adapting the order of history integral evaluation is important for obtaining accurate results. To examine the efficiency of our method, we compare runtimes to a state-of-the-art fixed time step integrator for several test systems and show an order of magnitude speedup at similar levels of accuracy. Published by the American Physical Society 2024

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Investigating the magnetic and atomic interface configuration for a model Fe/CrN bilayer system

A bilayer of iron on chromium nitride (Fe/CrN) is an interesting system for exchange biasing and sensing applications as the Néel temperature of CrN is 280 K and the Curie temperature of Fe is 1043 K. In this paper, we study the crystal and magnetic structures of the Fe/CrN interface at the atomic level. High quality epitaxial Fe/CrN bilayers prepared by molecular beam epitaxy grow in 001 orientation on MgO(001) substrates with uniform layer thicknesses and sharp interfaces. Our data reveal the epitaxial correlation between Fe and CrN crystals and their magnetic structures at the interface. The magnetic anisotropy directions of Fe and CrN are found parallel to [110] MgO . We studied the electronic and magnetic properties of the interface by performing the first-principles total-energy calculations. Here, we present a model that combines the crystal and magnetic structures of the Fe/CrN bilayer and fully explains all results.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Targeted Dy intercalation under graphene/SiC for tuning its electronic band structure

Metal intercalation of graphene is a promising method to tune its electronic band structure and generate novel electronic and topological phases. The tuning depends critically on the ability to bond the intercalated atoms at predesigned, subsurface interlayer locations because the emerging band structure depends on metal location. In this work, we have studied Dy intercalation under single-layer graphene (SLG) on SiC using spot profile analysis–low-energy electron diffraction and scanning tunneling microscopy (STM). The experimental work is complemented with density-functional theory (DFT) analysis. Because different diffraction spots originate from different subsurface interlayer regions, it is possible to identify changes in the intercalation location by monitoring the spot intensity as a function of growth conditions. DFT calculations of the chemical potential as a function of intercalated Dy coverage support the variation of the stability of the intercalated phase at different intercalated locations. The preferred location is confirmed from STM studies showing the removal of the 6 × 6 moiré corrugation at the preferred location, observed at higher Dy coverage.

2-dimensional systems↗
DOE OSTI2024

Solution of the Schrödinger equation for quasi-one-dimensional materials using helical waves

We formulate and implement a spectral method for solving the Schrödinger equation, as it applies to quasi-one-dimensional materials and structures. This allows for computation of the electronic structure of important technological materials such as nanotubes (of arbitrary chirality), nanowires, nanoribbons, chiral nanoassemblies, nanosprings and nanocoils, in an accurate, efficient and systematic manner. Our work is motivated by the observation that one of the most successful methods for carrying out electronic structure calculations of bulk/crystalline systems — the plane-wave method — is a spectral method based on eigenfunction expansion. Our scheme avoids computationally onerous approximations involving periodic supercells often employed in conventional plane-wave calculations of quasi-one-dimensional materials, and also overcomes several limitations of other discretization strategies, e.g., those based on finite differences and atomic orbitals. The basis functions in our method — called helical waves (or twisted waves) — are eigenfunctions of the Laplacian with symmetry adapted boundary conditions, and are expressible in terms of plane waves and Bessel functions in helical coordinates. We describe the setup of fast transforms to carry out discretization of the governing equations using our basis set, and the use of matrix-free iterative diagonalization to obtain the electronic eigenstates. Miscellaneous computational details, including the choice of eigensolvers, use of a preconditioning scheme, evaluation of oscillatory radial integrals and the imposition of a kinetic energy cutoff are discussed. We have implemented these strategies into a computational package called HelicES (Helical Electronic Structure). We demonstrate the utility of our method in carrying out systematic electronic structure calculations of various quasi-one-dimensional materials through numerous examples involving nanotubes, nanoribbons and nanowires. We also explore the convergence properties of our method, and assess its accuracy and computational efficiency by comparison against reference finite difference, transfer matrix method and plane-wave results. We anticipate that our method will find applications in computational nanomechanics and multiscale modeling, for carrying out transport calculations of interest to the field of semiconductor devices, and for the discovery of novel chiral phases of matter that are of relevance to the burgeoning quantum hardware industry.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Nano-photocurrent Mapping of Local Electronic Structure in Twisted Bilayer Graphene

In this work, we report a combined nano-photocurrent and infrared nanoscopy study of twisted bilayer graphene (TBG) enabling access to the local electronic phenomena at length scales as short as 20 nm. We show that the photocurrent changes sign at carrier densities tracking the local superlattice density of states of TBG. We use this property to identify domains of varying local twist angle by local photothermoelectric effect. Consistent with the photocurrent study, infrared nanoimaging experiments reveal optical conductivity features dominated by twist-angle-dependent interband transitions. Our results provide a fast and robust method for mapping the electronic structure of TBG and suggest that similar methods can be broadly applied to probe electronic inhomogeneities of Moiré superlattices in other van der Waals heterostructures.

36 MATERIALS SCIENCE↗