SEARCH · Engineering Papers

Results for “electronic structure method”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 235 records · Page 13

All sources NASA NTRS DOE OSTI

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 OSTI2026

Machine learning approach for vibronically renormalized electronic band structures

Here, we present a machine learning (ML) method for efficient computation of vibrational thermal expectation values of physical properties from first principles. Our approach is based on the nonperturbative frozen phonon formulation in which stochastic Monte Carlo algorithm is employed to sample configurations of nuclei in a supercell at finite temperatures based on a first-principles phonon model. A deep-learning neural network is trained to accurately predict physical properties associated with sampled phonon configurations, thus bypassing the time-consuming ab initio calculations. To incorporate the point-group symmetry of the electronic system into the ML model, group-theoretical methods are used to develop a symmetry-invariant descriptor for phonon configurations in the supercell. We apply our ML approach to compute the temperature dependent electronic energy gap of silicon based on density functional theory (DFT). We show that, with less than a hundred DFT calculations for training the neural network model, an order of magnitude larger number of sampling can be achieved for the computation of the vibrational thermal expectation values. Our work highlights the promising potential of ML techniques for finite temperature first-principles electronic structure methods.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2023

Nanoscale imaging of He-ion irradiation effects on amorphous TaO x toward electroforming-free neuromorphic functions

Resistive switching in thin films has been widely studied in a broad range of materials. Yet, the mechanisms behind electroresistive switching have been persistently difficult to decipher and control, in part due to their non-equilibrium nature. Here, we demonstrate new experimental approaches that can probe resistive switching phenomena, utilizing amorphous TaO x as a model material system. Specifically, we applied scanning microwave impedance microscopy and cathodoluminescence (CL) microscopy as direct probes of conductance and electronic structure, respectively. These methods provide direct evidence of the electronic state of TaO x despite its amorphous nature. For example, CL identifies characteristic impurity levels in TaO x , in agreement with first principles calculations. We applied these methods to investigate He-ion-beam irradiation as a path to activate conductivity of materials and enable electroforming-free control over resistive switching. However, we find that even though He-ions begin to modify the nature of bonds even at the lowest doses, the films' conductive properties exhibit remarkable stability with large displacement damage and they are driven to metallic states only at the limit of structural decomposition. Finally, we show that electroforming in a nanoscale junction can be carried out with a dissipated power of <20 nW, a much smaller value compared to earlier studies and one that minimizes irreversible structural modifications of the films. Importantly, the multimodal approach described here provides a new framework toward the theory/experiment guided design and optimization of electroresistive materials.

36 MATERIALS SCIENCE↗

NASA NTRS2004

Curing Composite Materials Using Lower-Energy Electron Beams

In an improved method of fabricating composite-material structures by laying up prepreg tapes (tapes of fiber reinforcement impregnated by uncured matrix materials) and then curing them, one cures the layups by use of beams of electrons having kinetic energies in the range of 200 to 300 keV. In contrast, in a prior method, one used electron beams characterized by kinetic energies up to 20 MeV. The improved method was first suggested by an Italian group in 1993, but had not been demonstrated until recently. With respect to both the prior method and the present improved method, the impetus for the use of electron- beam curing is a desire to avoid the high costs of autoclaves large enough to effect thermal curing of large composite-material structures. Unfortunately, in the prior method, the advantages of electron-beam curing are offset by the need for special walls and ceilings on curing chambers to shield personnel from x rays generated by impacts of energetic electrons. These shields must be thick [typically 2 to 3 ft (about 0.6 to 0.9 m) if made of concrete] and are therefore expensive. They also make it difficult to bring large structures into and out of the curing chambers. Currently, all major companies that fabricate composite-material spacecraft and aircraft structures form their layups by use of automated tape placement (ATP) machines. In the present improved method, an electron-beam gun is attached to an ATP head and used to irradiate the tape as it is pressed onto the workpiece. The electron kinetic energy between 200 and 300 keV is sufficient for penetration of the ply being laid plus one or two of the plies underneath it. Provided that the electron-beam gun is properly positioned, it is possible to administer the required electron dose and, at the same time, to protect personnel with less shielding than is needed in the prior method. Adequate shielding can be provided by concrete walls 6 ft (approximately equal to 1.8 m) high and 16 in. (approximately equal to 41 cm) thick, without a ceiling. The success of the present method depends on the use of a cationic epoxy as the matrix material in the prepreg tape, heating the prepreg tape to a temperature of 50 C immediately prior to layup, and exposing the workpiece to an electron-beam dose of approximately 2 Mrad. Experiments have shown that structures fabricated by the present method have the same mechanical properties as those of nominally identical structures fabricated by the prior method with electron beams of 3 to 4 MeV.

Byrne, Catherine A.↗

DOE OSTI2025

Convergent Protocols for Computing Protein–Ligand Interaction Energies Using Fragment-Based Quantum Chemistry

Fragment-based quantum chemistry methods offer a way to sidestep the steep nonlinear scaling of electronic structure calculations so that large molecular systems can be investigated using high-level methods. Here, we use fragmentation to compute protein–ligand interaction energies in systems with several thousand atoms, using a new software platform for managing fragment-based calculations that implements a screened many-body expansion. Convergence tests using a minimal-basis semiempirical method (HF-3c) indicate that two-body calculations, with single-residue fragments and simple hydrogen caps, are sufficient to reproduce interaction energies obtained using conventional supramolecular electronic structure calculations, to within 1 kcal/mol at about 1% of the computational cost. We also demonstrate that the HF-3c results are illustrative of trends obtained with density functional theory in basis sets up to augmented quadruple-ζ quality. Strategic deployment of fragmentation facilitates the use of converged biomolecular model systems alongside high-quality electronic structure methods and basis sets, bringing ab initio quantum chemistry to systems of hitherto unimaginable size. This will be useful for generation of high-quality training data for machine learning applications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Density-functional theory of material design: fundamentals and applications-I

This article is part-I of a review of density-functional theory (DFT) that is the most widely used method for calculating electronic structure of materials. The accuracy and ease of numerical implementation of DFT methods has resulted in its extensive use for materials design and discovery and has thus ushered in the new field of computational material science. In this article, we start with an introduction to Schrödinger equation and methods of its solutions. After presenting exact results for some well-known systems, difficulties encountered in solving the equation for interacting electrons are described. How these difficulties are handled using the variational principle for the energy to obtain approximate solutions of the Schrödinger equation is discussed. The resulting Hartree and Hartree–Fock theories are presented along with results they give for atomic and solid-state systems. We then describe Thomas–Fermi theory and its extensions which were the initial attempts to formulate many-electron problem in terms of electronic density of a system. Having described these theories, we introduce modern DFT by discussing Hohenberg–Kohn theorems that form its foundations. We then go on to discuss Kohn–Sham (KS) formulation of DFT in its exact form. Next, local density approximation (LDA) is introduced and solutions of KS equation for some representative systems, obtained using the LDA, are presented. We end part-I of the review describing the contents of part-II.

36 MATERIALS SCIENCE↗

DOE OSTI2020

Toward a systematic improvement of the fixed-node approximation in diffusion Monte Carlo for solids—A case study in diamond

While Diffusion Monte Carlo (DMC) is in principle an exact stochastic method for ab initio electronic structure calculations, in practice, the fermionic sign problem necessitates the use of the fixed-node approximation and trial wavefunctions with approximate nodes (or zeros). This approximation introduces a variational error in the energy that potentially can be tested and systematically improved. Here, we present a computational method that produces trial wavefunctions with systematically improvable nodes for DMC calculations of periodic solids. These trial wavefunctions are efficiently generated with the configuration interaction using a perturbative selection made iteratively (CIPSI) method. A simple protocol in which both exact and approximate results for finite supercells are used to extrapolate to the thermodynamic limit is introduced. This approach is illustrated in the case of the carbon diamond using Slater–Jastrow trial wavefunctions including up to one million Slater determinants. Fixed-node DMC energies obtained with such large expansions are much improved, and the fixed-node error is found to decrease monotonically and smoothly as a function of the number of determinants in the trial wavefunction, a property opening the way to a better control of this error. The cohesive energy extrapolated to the thermodynamic limit is in close agreement with the estimated experimental value. Interestingly, this is also the case at the single-determinant level, thus, indicating a very good error cancellation in carbon diamond between the bulk and atomic total fixed-node energies when using single-determinant nodes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS1983

Studies in electron phenomena in MOS structures: The pulsed C-V method

The pulse hysteresis capacitance voltage (C-V) provides a straight forward technique for measuring the change of various charges in MOS structures and a tool for investigating the kinetics of various electron phenomena is developed and described. The method can be used for measuring the energy distribution and kinetics of surface states with the resolution of about 1/5 x 10 to the -9 power cm eV. Some transients in an MOS structure, particularly, the thermal generation of minority charge carriers via surface states and the relaxation of minority charge carriers supplied from the inversion layer outside the MOS structure are theoretically investigated. Analytical expressions which clearly present the physics of those electron phenomena are derived.

Kaplan, G.↗

DOE OSTI2025

Surface preparation method for investigating the three-dimensional electronic structure of perovskite nickelates

The investigation of the electronic structures on perovskite oxides using surface-sensitive spectroscopy techniques is often hindered by their “uncleavable” nature, typically requiring expensive and complex in situ experimental setups that integrate the capabilities of sample synthesis and spectroscopy measurement under ultrahigh vacuum condition. Here, we address this challenge by developing an ozone-annealing process that yields atomically flat surfaces on perovskite oxide thin films, making them suitable for high-resolution angle-resolved photoemission spectroscopy measurements. Using this method, we present a three-dimensional electronic structure study of Nd 1−𝑥 ⁢Sr 𝑥 ⁢NiO 3 (𝑥=0 and 0.175) thin films with unprecedented accuracy. The experimentally determined low-energy fermiology exhibits quantitative agreements with two-band tight-binding simulations, which is further validated by first-principles calculations considering the material's actual crystal structure. This work provides an accessible approach for ex situ ARPES measurements on perovskite oxides and other strongly correlated oxides, including the recently discovered high-𝑇 c nickelates.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

NASA NTRS1988

Symmetry and equivalence restrictions in electronic structure calculations

A simple method for obtaining MCSCF orbitals and CI natural orbitals adapted to degenerate point groups, with full symmetry and equivalnece restrictions, is described. Among several advantages accruing from this method are the ability to perform atomic SCF calculations on states for which the SCF energy expression cannot be written in terms of Coulomb and exchange integrals over real orbitals, and the generation of symmetry-adapted atomic natural orbitals for use in a recently proposed method for basis set contraction.

Bauschlicher, Charles W., Jr.↗

DOE OSTI2021

Assessment of DFT methods for the prediction of detachment energies and electronic structures of complex and multiply charged anions

Here, in this work, we assess the ability of different density functional theory methods to reproduce the vertical/adiabatic detachment energies (VDEs/ADEs) of a series of 40 anions, including monoanion versus multiply charged anions (MCAs) and covalent versus noncovalent interaction. The statistical analysis of errors is firstly performed by comparing the theoretical values with the experimental benchmark data obtained from the negative ion photoelectron spectroscopy (NIPES). The proposed optimally tuned range-separated (OTRS) functionals are proved to not only well reproduce the experimental VDEs/ADEs, but also simulate experimental NIPES. The result implies a more serious issue of delocalization error for MCAs and a balanced description of electron-rich/- deficient region of anions is necessary. The radius of the spherically symmetric average electron localization function (ELF) region is demonstrated as a useful descriptor for the characterization of electronic structure of anionic systems and a quasilinear fitting model is proposed to efficiently obtain their OTRS parameters.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Overdestabilization vs Overstabilization in the Theoretical Analysis of f-Orbital Covalency

The complex nature of the f-orbital electronic structures and their interaction with the chemical environment pose significant computational challenges. Advanced computational techniques that variationally include scalar relativities and spin–orbit coupling directly at the molecular orbital level have been developed to address this complexity. Among these, variational relativistic multiconfigurational multireference methods stand out for their high accuracy and systematic improvement in studies of f-block complexes. Additionally, these advanced methods offer the potential for calibrating low-scaling electronic structure methods such as density functional theory. However, studies on the Cl K-edge X-ray absorption spectra of the [Ce(III)Cl 6 ] 3– and [Ce(IV)Cl 6 ] 2– complexes show that time-dependent density functional theory with approximate exchange–correlation kernels can lead to inaccuracies, resulting in an overstabilization of 4f orbitals and incorrect assessments of covalency. In contrast, approaches utilizing small active space wave function methods may understate the stability of these orbitals. The results herein demonstrate the need for large active space, multireference, and variational relativistic methods in studying f-block complexes.

chemical calculations↗

DOE OSTI2024

Shubnikov–de Haas oscillations of biaxial-strain-tuned superconductors in pulsed magnetic field up to 60 T

Two-dimensional (2D) materials have gained increasing prominence not only in fundamental research but also in daily applications. However, to fully harness their potential, it is crucial to optimize their properties with an external parameter and track the electronic structure simultaneously. Magnetotransport over a wide magnetic field range is a powerful method to probe the electronic structure and, for metallic 2D materials, quantum oscillations superimposed on the transport signals encode Fermi surface parameters. In this manuscript, we utilize biaxial strain as an external tuning parameter and investigate the effects of strain on the electronic properties of two quasi-2D superconductors, MoTe 2 and RbV 3 Sb 5 , by measuring their magnetoresistance in pulsed magnetic fields up to 60 T. With a careful selection of insulating substrates, we demonstrate the possibility of both the compressive and tensile biaxial strains imposed on MoTe 2 and RbV 3 Sb 5 , respectively. For both systems, the applied strain has led to superconducting critical temperature enhancement compared to their free-standing counterparts, proving the effectiveness of this biaxial strain method at cryogenic temperatures. Clear quantum oscillations in the magnetoresistance—the Shubnikov–de Haas (SdH) effect—are obtained in both samples. In strained MoTe 2 , the magnetoresistance exhibits a nearly quadratic dependence on the magnetic field and remains non-saturating even at the highest field, whereas in strained RbV 3 Sb 5 , two SdH frequencies showed a substantial enhancement in effective mass values, hinting at a possible enhancement of charge fluctuations. Our results demonstrate that combining biaxial strain and pulsed magnetic field paves the way for studying 2D materials under unprecedented conditions.

2D materials↗

DOE OSTI2024

Shubnikov-de Haas oscillations of biaxial-strain-tuned superconductors in pulsed magnetic field up to 60 T

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2021

Understanding flavin electronic structure and spectra

Abstract Flavins have emerged as central to electron bifurcation, signaling, and countless enzymatic reactions. In bifurcation, two electrons acquired as a pair are separated in coupled transfers wherein the energy of both is concentrated on one of the two. This enables organisms to drive demanding reactions based on abundant low‐grade chemical fuel. To enable incorporation of this and other flavin capabilities into designed materials and devices, it is essential to understand fundamental principles of flavin electronic structure that make flavins so reactive and tunable by interactions with protein. Emerging computational tools can now replicate spectra of flavins and are gaining capacity to explain reactivity at atomistic resolution, based on electronic structures. Such fundamental understanding can moreover be transferrable to other chemical systems. A variety of computational innovations have been critical in reproducing experimental properties of flavins including their electronic spectra, vibrational signatures, and nuclear magnetic resonance (NMR) chemical shifts. A computational toolbox for understanding flavin reactivity moreover must be able to treat all five oxidation and protonation states, in addition to excited states that participate in flavoprotein's light‐driven reactions. Therefore, we compare emerging hybrid strategies and their successes in replicating effects of hydrogen bonding, the surrounding dielectric, and local electrostatics. These contribute to the protein's ability to modulate flavin reactivity, so we conclude with a survey of methods for incorporating the effects of the protein residues explicitly, as well as local dynamics. Computation is poised to elucidate the factors that affect a bound flavin's ability to mediate stunningly diverse reactions, and make life possible. This article is categorized under: Structure and Mechanism > Computational Biochemistry and Biophysics Electronic Structure Theory > Combined QM/MM Methods Theoretical and Physical Chemistry > Spectroscopy

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2023

TAMM: Tensor algebra for many-body methods

Tensor algebra operations such as contractions in computational chemistry consume a significant fraction of the computing time on large-scale computing platforms. The widespread use of tensor contractions between large multi-dimensional tensors in describing electronic structure theory has motivated the development of multiple tensor algebra frameworks targeting heterogeneous computing platforms. In this paper, we present Tensor Algebra for Many-body Methods (TAMM), a framework for productive and performance-portable development of scalable computational chemistry methods. TAMM decouples the specification of the computation from the execution of these operations on available high-performance computing systems. With this design choice, the scientific application developers (domain scientists) can focus on the algorithmic requirements using the tensor algebra interface provided by TAMM, whereas high-performance computing developers can direct their attention to various optimizations on the underlying constructs, such as efficient data distribution, optimized scheduling algorithms, and efficient use of intra-node resources (e.g., graphics processing units). The modular structure of TAMM allows it to support different hardware architectures and incorporate new algorithmic advances. We describe the TAMM framework and our approach to the sustainable development of scalable ground- and excited-state electronic structure methods. We present case studies highlighting the ease of use, including the performance and productivity gains compared to other frameworks.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗