Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “supercell”

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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Revealing intrinsic domains and fluctuations of moiré magnetism by a wide-field quantum microscope

Moiré magnetism featured by stacking engineered atomic registry and lattice interactions has recently emerged as an appealing quantum state of matter at the forefront of condensed matter physics research. Nanoscale imaging of moiré magnets is highly desirable and serves as a prerequisite to investigate a broad range of intriguing physics underlying the interplay between topology, electronic correlations, and unconventional nanomagnetism. Here we report spin defect-based wide-field imaging of magnetic domains and spin fluctuations in twisted double trilayer (tDT) chromium triiodide CrI 3 . We explicitly show that intrinsic moiré domains of opposite magnetizations appear over arrays of moiré supercells in low-twist-angle tDT CrI 3 . In contrast, spin fluctuations measured in tDT CrI 3 manifest little spatial variations on the same mesoscopic length scale due to the dominant driving force of intralayer exchange interaction. Our results enrich the current understanding of exotic magnetic phases sustained by moiré magnetism and highlight the opportunities provided by quantum spin sensors in probing microscopic spin related phenomena on two-dimensional flatland.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Electron–phonon physics from first principles using the EPW code

EPW is an open-source software for ab initio calculations of electron–phonon interactions and related materials properties. The code combines density functional perturbation theory and maximally localized Wannier functions to efficiently compute electron–phonon coupling matrix elements, and to perform predictive calculations of temperature-dependent properties and phonon-assisted quantum processes in bulk solids and low-dimensional materials. Here, we report on significant developments in the code since 2016, namely: a transport module for the calculation of charge carrier mobility under electric and magnetic fields using the Boltzmann transport equation; a superconductivity module for calculations of phonon-mediated superconductors using the anisotropic multi-band Eliashberg theory; an optics module for calculations of phonon-assisted indirect transitions; a module for the calculation of small and large polarons without supercells; and a module for calculating band structure renormalization and temperature-dependent optical spectra using the special displacement method. For each capability, we outline the methodology and implementation and provide example calculations.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Trajectory sampling and finite-size effects in first-principles stopping power calculations

Abstract Real-time time-dependent density functional theory (TDDFT) is presently the most accurate available method for computing electronic stopping powers from first principles. However, obtaining application-relevant results often involves either costly averages over multiple calculations or ad hoc selection of a representative ion trajectory. We consider a broadly applicable, quantitative metric for evaluating and optimizing trajectories in this context. This methodology enables rigorous analysis of the failure modes of various common trajectory choices in crystalline materials. Although randomly selecting trajectories is common practice in stopping power calculations in solids, we show that nearly 30% of random trajectories in an FCC aluminum crystal will not representatively sample the material over the time and length scales feasibly simulated with TDDFT, and unrepresentative choices incur errors of up to 60%. We also show that finite-size effects depend on ion trajectory via “ouroboros” effects beyond the prevailing plasmon-based interpretation, and we propose a cost-reducing scheme to obtain converged results even when expensive core-electron contributions preclude large supercells. This work helps to mitigate poorly controlled approximations in first-principles stopping power calculations, allowing 1–2 order of magnitude cost reductions for obtaining representatively averaged and converged results.

36 MATERIALS SCIENCE↗
DOE OSTI2024

A high-throughput framework for lattice dynamics

We develop an automated high-throughput workflow for calculating lattice dynamical properties from first principles including those dictated by anharmonicity. The pipeline automatically computes interatomic force constants (IFCs) up to 4th order from perturbed training supercells, and uses the IFCs to calculate lattice thermal conductivity, coefficient of thermal expansion, and vibrational free energy and entropy. It performs phonon renormalization for dynamically unstable compounds to obtain real effective phonon spectra at finite temperatures and calculates the associated free energy corrections. The methods and parameters are chosen to balance computational efficiency and result accuracy, assessed through convergence testing and comparisons with experimental measurements. Deployment of this workflow at a large scale would facilitate materials discovery efforts toward functionalities including thermoelectrics, contact materials, ferroelectrics, aerospace components, as well as general phase diagram construction.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

Optical properties of a diamond NV color center from capped embedded multiconfigurational correlated wavefunction theory

Diamond defects are among the most promising qubits. Modeling their properties through accurate quantum mechanical simulations can further their development into robust units of information. We use the recently developed capped density functional embedding theory (capped-DFET) with the multiconfigurational n-electron valence second-order perturbation theory to characterize the electronic excitation energies for different spin manifolds of the well-characterized negatively charged substitutional N defect adjacent to a vacancy (V C ) in diamond (N C V C − ). We successfully reproduce vertical excitation energies for both triplet and singlet states of N C V C − with errors < 0.1 eV. Unlike other embedding methods, capped-DFET exhibits robust predictions that are approximately independent of the embedded cluster size: it only requires a cluster to contain the defect atoms and their nearest neighbors (as small as a 40-atom capped cluster). Furthermore, our method is free from slowly converging Coulomb interactions between charged defects, and thus also only weakly dependent on supercell size.

Chemistry↗
DOE OSTI2023

Geometric frustration of Jahn–Teller order in the infinite-layer lattice

The Jahn–Teller effect, in which electronic configurations with energetically degenerate orbitals induce lattice distortions to lift this degeneracy, has a key role in many symmetry-lowering crystal deformations. Lattices of Jahn–Teller ions can induce a cooperative distortion, as exemplified by LaMnO 3 . Although many examples occur in octahedrally or tetrahedrally coordinated transition metal oxides due to their high orbital degeneracy, this effect has yet to be manifested for square-planar anion coordination, as found in infinite-layer copper, nickel, iron and manganese oxides. Here, for this work, we synthesize single-crystal CaCoO 2 thin films by topotactic reduction of the brownmillerite CaCoO 2.5 phase. We observe a markedly distorted infinite-layer structure, with ångström-scale displacements of the cations from their high-symmetry positions. This can be understood to originate from the Jahn–Teller degeneracy of the d xz and d yz orbitals in the d 7 electronic configuration along with substantial ligand–transition metal mixing. A complex pattern of distortions arises in a $2\sqrt2$ x $2\sqrt2$ x $1$ tetragonal supercell, reflecting the competition between an ordered Jahn–Teller effect on the CoO 2 sublattice and the geometric frustration of the associated displacements of the Ca sublattice, which are strongly coupled in the absence of apical oxygen. As a result of this competition, the CaCoO 2 structure forms an extended two-in–two-out type of Co distortion following ‘ice rules’.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Accelerating computational modeling and design of high-entropy alloys

High-entropy alloys, with N elements and compositions {$c_{ν = 1,N}$} in competing crystal structures, have large design spaces for unique chemical and mechanical properties. In this work, to enable computational design, we use a metaheuristic hybrid Cuckoo search (CS) to construct alloy configurational models on the fly that have targeted atomic site and pair probabilities on arbitrary crystal lattices, given by supercell random approximates (SCRAPs) with S sites. Our Hybrid CS permits efficient global solutions for large, discrete combinatorial optimization that scale linearly in a number of parallel processors, and linearly in sites S for SCRAPs. For example, a four-element, 128-site SCRAP is found in seconds—a more than 13,000-fold reduction over current strategies. Our method thus enables computational alloy design that is currently impractical. We qualify the models and showcase application to real alloys with targeted atomic short-range order. Being problem-agnostic, our Hybrid CS offers potential applications in diverse fields.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Defect graph neural networks for materials discovery in high-temperature clean-energy applications

We present a graph neural network approach that fully automates the prediction of defect formation enthalpies for any crystallographic site from the ideal crystal structure, without the need to create defected atomic structure models as input. Here we used density functional theory reference data for vacancy defects in oxides, to train a defect graph neural network (dGNN) model that replaces the density functional theory supercell relaxations otherwise required for each symmetrically unique crystal site. Interfaced with thermodynamic calculations of reduction entropies and associated free energies, the dGNN model is applied to the screening of oxides in the Materials Project database, connecting the zero-kelvin defect enthalpies to high-temperature process conditions relevant for solar thermochemical hydrogen production and other energy applications. The dGNN approach is applicable to arbitrary structures with an accuracy limited principally by the amount and diversity of the training data, and it is generalizable to other defect types and advanced graph convolution architectures. In conclusion, it will help to tackle future materials discovery problems in clean energy and beyond.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Molecular dynamics simulations for glass transition temperature predictions of polyhydroxyalkanoate biopolymers

Polyhydroxyalkanoates (PHAs) represent an emerging class of biosynthetic and biodegradable polyesters that exhibit considerable potential to replace petroleum-based plastics towards a sustainable future. Despite the promise, general structure–property mappings within this class of polymers remain largely unexplored. An efficient exploration of this vast chemical space calls for the development and validation of predictive methods for accurate estimation of a diverse range of properties for PHA-based polymers. Towards this aim, we present and validate in this work the results of our molecular dynamics (MD) simulation based approach aimed at predicting glass transition temperatures (T g ) of PHA-based polymers. Since generally available and widely used polymer forcefields exhibit a relatively poor performance for T g predictions, we have developed a new forcefield by modifying the polymer consistent force field (PCFF) via refining a selected set of torsion potentials of the polymer backbone using accurate density functional theory (DFT) computations. After carefully assessing the dependence of critical simulation parameters, such as, polymer chain length, number of polymer chains, supercell size, and thermal quenching rate used in the simulation, the applicability and transferability of the modified PCFF (mPCFF) is demonstrated by directly comparing the computed T g predictions of various polymers with different chemistries, polymer side chain lengths and functional groups forming the polymer side chains against the respective experimentally measured values. Furthermore, the transport properties such as self-diffusion coefficient and viscosity are computationally determined and their well-known correlation with the target properties is demonstrated. Lastly, we have employed the developed approach to predict T g values for a number of yet-to-be-synthesized PHA-based polymers with a diverse set of functional groups in the polymer side chains. The results are further rationalized by correlating the predicted T g values with the inter-chain H-bond formation tendencies of the different side chain functional groups. This work represents an important first step towards computationally guided design of PHA-based functional polymers and opens up new directions for a systematic investigation of composition- and configuration-dependent structure–property relationships in more complex binary and ternary copolymer systems.

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

Computational investigations of Dienes defect- and vacancy-induced changes in the electronic and vibrational properties of carbon fiber structural units

Carbon fiber (CF) is a promising lightweight alternative to steel and is of significant interest for energy applications. As CF continues to find new uses and is exposed to new external conditions, a noninvasive method of monitoring its structural integrity is critical. Raman spectroscopy is a commonly used method for this monitoring; however, it is highly inferential, and the interpretation of the data is not always straightforward. In this work, we perform density functional theory (DFT) calculations to investigate changes in the vibrational properties of CF structural units (i.e., graphene and graphite) caused by monovacancy and Dienes defects as a foundation for modeling more complex defects that move our model toward that of realistic CF. Using large computational supercells, we can understand how these defects change the electronic structure and vibrational properties of graphene and graphite for interdefect distances near those of the lower experimental limit. The monovacancy opens an electronic bandgap at the K point. Although no such electronic gap is opened by the Dienes defect, both defects introduce flat defect bands near the Fermi energy. The Dienes defect creates long-range deviations of the phonons, leading to substantial broadening of the highest frequency optical modes in the band structure compared to that of the pristine material. In contrast, the phonon changes caused by the monovacancy are short range, and only minor changes in the band structure or phonon density of states were observed. These findings can assist in the interpretation of experimental results by providing atomic-scale insight into key electronic and vibrational features.

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

Evaluation of polymorphism and charge transport in a BaO–CaO–Ta 2 O 5 perovskite phase diagram using TOF-neutron and synchrotron X-ray diffraction, the bond-valence method and impedance spectroscopy

In the present work, we develop a comprehensive functional phase diagram for the Ba–Ca–Ta–O quaternary system Ba 3 Ca 1+ x Ta 2− x O 9−3 x /2 (0 ≤ x ≤ 0.36) between 1000 and 1550 °C, coupled with theoretical calculations of the cationic ordering in supercells.

Singh, Kalpana↗
DOE OSTI2024

GeSn-on-GaAs with photoconductive carrier lifetime >400 ns: role of substrate orientation and atomistic simulation

Group IV GeSn quantum material finds application in electronics and silicon-compatible photonics. Synthesizing these materials with low defect density and high carrier lifetime is a potential challenge due to lattice mismatch induced defects and tin segregation at higher growth temperature. Recent advancements in the growth of these GeSn materials on Si, Ge, GaAs, and with substrate orientations, demonstrated different properties using epitaxial and chemical deposition methods. This article addresses the effect of GaAs substrate orientation and misorientation on the materials’ properties and carrier lifetimes in epitaxial Ge 0.94 Sn 0.06 layers. With starting GaAs substrates of (100)/2°, (100)/6°, (110) and (111)A orientations, Ge 0.94 Sn 0.06 epitaxial layers were grown with an intermediate Ge buffer layer by molecular beam epitaxy and analyzed by several analytical tools. X-ray analysis displayed good crystalline quality, and Raman spectroscopy measurements showed blue shifts in phonon wavenumber due to biaxial compressive strain in Ge 0.94 Sn 0.06 epilayers. Cross-sectional transmission electron microscopy analysis confirmed the defect-free heterointerface of Ge 0.94 Sn 0.06 /Ge/GaAs heterostructure. Minority carrier lifetimes of the unintentionally doped n-type Ge 0.94 Sn 0.06 epilayers displayed photoconductive carrier lifetimes of >400 ns on (100)/6°, 319 ns on (100)/2°, and 434 ns on (110) GaAs substrate at 1500 nm excitation wavelength. On the other hand, Ge 0.94 Sn 0.06 layer showed poor carrier lifetime on (111)A GaAs substrate. The observed differences in carrier lifetimes were correlated with the formation energy of the Ge on (100)/6° and (100)/2° GaAs heterointerface using Stillinger–Weber interatomic potential model-based atomistic simulation with different heterointerfacial bonding by Synopsys QuantumATK tool. Total energy computation of 6280-atom Ge/GaAs supercell on (100)/6° leads to lower formation energy than (100)/2°, making it more thermodynamically stable. Hence, the growth of the GeSn/III–V material system using misoriented (100) substrates that are more thermodynamically stable will enhance the performances of optoelectronic devices.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Porosity modeling in a TiNbTaZrMo high-entropy alloy for biomedical applications

High-entropy alloys (HEAs) have attracted great attention for many biomedical applications. However, the nature of interatomic interactions in this class of complex multicomponent alloys is not fully understood. We report, for the first time, the results of theoretical modeling for porosity in a large biocompatible HEA TiNbTaZrMo using an atomistic supercell of 1024 atoms that provides new insights and understanding. Our results demonstrated the deficiency of using the valence electron count, quantification of large lattice distortion, validation of mechanical properties with available experimental data to reduce Young's modulus. We utilized the novel concepts of the total bond order density (TBOD) and partial bond order density (PBOD) via ab initio quantum mechanical calculations as an effective theoretical means to chart a road map for the rational design of complex multicomponent HEAs for biomedical applications.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Resolving local ordering and structure in Mn x Ge 1- x Te alloys through thermodynamic ensembles of pair distribution functions

Characterizing local bonding environments in complex materials is essential for understanding and optimizing their properties. Equally as important is the ability to predict local motifs as a function of synthesis conditions, enhancing chemists’ ability to design properties into materials. In this study, we present an approach to leverage statistical mechanics to generate temperature- and energy-informed ensemble averaged pair distribution functions (PDFs). This method, which we have named Thermodynamic Ensemble Averages of PDFs for Ordering and Transformations (TEAPOT), utilizes density functional theory (DFT) to relax supercells while incorporating energetic penalties for local order, enabling accurate and computationally efficient analysis of local structure. We apply this method to the neutron PDF measurements of the pseudobinary MnTe–GeTe (MGT) alloy, demonstrating its capability to resolve complex local distortions and chemical ordering. Our results reveal detailed insights into phase transformations and local distortions driven by Mn substitution. For compositions that globally present as rock salt, our analysis reveals that Ge coordination geometry is heavily impacted by synthesis temperature. We propose that high temperature synthesis conditions promote a lowered Ge polyhedra distortion, promoting high charge carrier mobility due to the alignment of local and global structure. Incorporating statistical mechanics and computation into experimental analysis thus guides synthesis of tailored local structure.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Re-evaluation of the TSL for Yttrium Hydride

Yttrium hydride (YH x ) is of interest as a high-temperature moderator material in advanced nuclear reactor systems because of its superior ability to retain hydrogen at elevated temperatures. Thermal neutron scattering laws (TSL) for hydrogen bound in yttrium hydride (H-YH 2 ) and yttrium bound in yttrium hydride (Y-YH 2 ) were previously evaluated by Naval Nuclear Laboratory using the ab initio approach and released in ENDF/B-VIII.0. In that work, density functional theory, incorporating the generalized gradient approximation (GGA) for the exchange-correlation energy, was used to simulate the face-centered cubic structure of YH 2 and calculate the interatomic Hellmann-Feynman forces for a 2×2×2 supercell containing 96 atoms. Lattice dynamics calculations using PHONON were used to determine the phonon density of states. The calculated phonon density of states for H and Y in YH 2 were then used to prepare H-YH 2 and Y-YH 2 TSL evaluations, in the incoherent approximation, using the LEAPR module of NJOY. In addition, elastic scattering was assumed to be incoherent for both H and Y. While the incoherent elastic scattering approximation is appropriate for H-YH 2 , it introduces an undesirable approximation for Y-YH 2 . In this work, we re-evaluate the TSL for Y-YH 2 using FLASSH (Full Law Analysis Scattering System Hub). Y-YH 2 is evaluated using the FLASSH generalized coherent elastic scattering capability in order to capture the Bragg peaks associated with the YH 2 crystal structure which were neglected in the prior NJOY-based evaluation due to limitations in LEAPR. An experimental approach to validate the Y-YH 2 TSL using neutron transmission measurements is discussed.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗