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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Magnetic Response of Excitons and Excitonic Complexes in Defective Hexyl Ammonium Lead Iodide Self-Assembled Quantum Wells

Correlating the magnetic behavior of low-dimensional semiconductors with the chemistry used to form these excitonic materials remains crucial to the development of devices applied to quantum technologies. In this study, we apply large magnetic fields during the collection of low-temperature photoluminescence (PL) to assess the magnetic properties of excitons in hexyl ammonium lead iodide (HA 2 PbI 4 ) self-assembled quantum wells (SAQWs) formed at liquid−liquid interfaces. The effect of incident laser power and temperature was used to assign lower energy features in these samples’ PL spectra to both excitons trapped at defect sites and trions. Measured peaks shifts allow us to estimate coupling between the applied magnetic fields and the charged defect excitons. Our conclusions are supported by density functional theory calculations on supercells of the proposed defective HA 2 PbI 4 SAQW structure. Additionally, we find an anomalous magnetic response from trion states that resembles the behavior of interacting excitonic complexes in similar, quantum-confined materials at low temperatures. These results highlight the crucial role that chemical conditions can play in the magnetic response of 2D semiconductors applied in quantum technologies.

defect chemistry↗
DOE OSTI2021

Temperature-Dependent Properties of Molten Li 2 BeF 4 Salt Using Ab Initio Molecular Dynamics

Molten lithium tetrafluoroberyllate (Li 2 BeF 4 ) salt, also known as FLiBe, with a 2:1 mixture of LiF and BeF 2 is being proposed as a coolant and solvent in advanced nuclear reactor designs, such as the molten salt reactor or the fluoride salt cooled high-temperature reactor. We present the results on the structure and properties of FLiBe over a wide range of temperatures, 0–2000 K, from high-throughput ab initio molecular dynamics simulation using a supercell model of 504 atoms. The variations in the local structures of solid and liquid FLiBe with temperature are discussed in terms of a pair distribution function, coordination number, and bond angle distribution. The temperature-dependent electronic structure and optical and mechanical properties of FLiBe are calculated. The optical and mechanical property results are reported for the first time. The results above and below the melting temperature (~732 K) are compared with the experimental data and with data for crystalline FLiBe. The electronic structure and interatomic bonding results are discussed in correlation with the mechanical strength. A novel concept of total bond order density (TBOD), an important quantum mechanical parameter, is used to characterize the internal cohesion and strength in the simulated models. The results show a variation in the rate of change in properties in solid and liquid phases with anomalous behavior across the melting region. The observed trend is the decrease in mechanical strength, band gap, and TBOD in a nonlinear fashion as a function of temperature. The refractive index shows a surprising minimum at 850 K, among the tested temperatures, which lies above the melting point. These findings provide a new platform to understand the interplay between the temperature-dependent structures and properties of FLiBe salt.

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

The Subchalcogenides Ir 2 In 8 Q (Q = S, Se, Te): Dirac Semimetal Candidates with Re-entrant Structural Modulation

Subchalcogenides are uncommon compounds where the metal atoms are in unusually low formal oxidation states. They bridge the gap between intermetallics and semiconductors and can have unexpected structures and properties because of the exotic nature of their chemical bonding as they contain both metal–metal and metal–main group (e.g., halide, chalcogenide) interactions. Finding new members of this class of materials presents synthetic challenges as attempts to make them often result in phase separation into binary compounds. We overcome this difficulty by utilizing indium as a metal flux to synthesize large (millimeter scale) single crystals of novel subchalcogenide materials. Herein, we report two new compounds Ir 2 In 8 Q (Q = Se, Te) and compare their structural and electrical properties to the previously reported Ir 2 In 8 S analogue. Ir 2 In 8 Se and Ir 2 In 8 Te crystallize in the $P_{4_2}$/mnm space group and are isostructural to Ir2In8S, but also have commensurately modulated (with q vectors q = 1/6a* + 1/6b* and q = 1/10a* + 1/10b* for Ir 2 In 8 Se and Ir 2 In 8 Te, respectively) low-temperature phase transitions, where the chalcogenide anions in the channels experience a distortion in the form of In–Q bond alternation along the ab plane. Both compounds display re-entrant structural behavior, where the supercells appear on cooling but revert to the original subcell below 100 K, suggesting competing structural and electronic interactions dictate the overall structure. Notably, these materials are topological semimetal candidates with symmetry-protected Dirac crossings near the Fermi level and exhibit high electron mobilities (~1500 cm 2 V –1 s –1 at 1.8 K) and moderate carrier concentrations (~10 20 cm –3 ) from charge transport measurements. This work highlights metal flux as a synthetic route to high quality single crystals of novel intermetallic subchalcogenides with Dirac semimetal behavior.

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

Expandable Li Percolation Network: The Effects of Site Distortion in Cation-Disordered Rock-Salt Cathode Material

Cation-disordered rock-salt (DRX) materials receive intensive attention as a new class of cathode candidates for high-capacity lithium-ion batteries (LIBs). Unlike traditional layered cathode materials, DRX materials have a three-dimensional (3D) percolation network for Li + transportation. The disordered structure poses a grand challenge to a thorough understanding of the percolation network due to its multiscale complexity. In this work, we introduce the large supercell modeling for DRX material Li 1.16 Ti 0.37 Ni 0.37 Nb 0.10 O 2 (LTNNO) via the reverse Monte Carlo (RMC) method combined with neutron total scattering. Here, through a quantitative statistical analysis of the material’s local atomic environment, we experimentally verified the existence of short-range ordering (SRO) and uncovered an element-dependent behavior of transition metal (TM) site distortion. A displacement from the original octahedral site for Ti 4+ cations is pervasive throughout the DRX lattice. Density functional theory (DFT) calculations revealed that site distortions quantified by the centroid offsets could alter the migration barrier for Li + diffusion through the tetrahedral channels, which can expand the previously proposed theoretical percolating network of Li. The estimated accessible Li content is highly consistent with the observed charging capacity. The newly developed characterization method here uncovers the expandable nature of the Li percolation network in DRX materials, which may provide valuable guidelines for the design of superior DRX materials.

25 ENERGY STORAGE↗
DOE OSTI2024

Correlated Anion Disorder in Heteroanionic Cubic TiOF 2

Resolving anion configurations in heteroanionic materials is crucial for understanding and controlling their properties. For anion-disordered oxyfluorides, conventional Bragg diffraction cannot fully resolve the anionic structure, necessitating alternative structure determination methods. We have investigated the anionic structure of anion-disordered cubic (ReO 3 -type) TiOF 2 using X-ray pair distribution function (PDF), 19 F MAS NMR analysis, density functional theory (DFT), cluster expansion modeling, and genetic-algorithm structure prediction. Our computational data predict short-range anion ordering in TiOF 2 , characterized by predominant cis-[O 2 F 4 ] titanium coordination, resulting in correlated anion disorder at longer ranges. To validate our predictions, we generated partially disordered supercells using genetic-algorithm structure prediction and computed simulated X-ray PDF data and 19 F MAS NMR spectra, which we compared directly to experimental data. To construct our simulated 19 F NMR spectra, we derived new transformation functions for mapping calculated magnetic shieldings to predicted magnetic chemical shifts in titanium (oxy)fluorides, obtained by fitting DFT-calculated magnetic shieldings to previously published experimental chemical shift data for TiF 4 . We find good agreement between our simulated and experimental data, which supports our computationally predicted structural model and demonstrates the effectiveness of complementary experimental and computational techniques in resolving anionic structure in anion-disordered oxyfluorides. From additional DFT calculations, we predict that increasing anion disorder makes lithium intercalation more favorable by, on average, up to 2 eV, highlighting the significant effect of variations in short-range order on the intercalation properties of anion-disordered materials.

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

Higher Dimensionality in the Mg–Co–B System: Synthesis and Structure of Incommensurate Composite Mg 1+ε Co 4 B 4

Guided by high-temperature in situ X-ray diffraction, the discovery and synthesis of Mg 1+ε Co 4 B 4 (ε ≈ 0.272) using a MgH 2 hydride precursor is reported, along with a detailed crystal structure description and measurement of magnetic properties. The mismatch in lattice periodicities between Mg and Co–B substructures places Mg 1+ε Co 4 B 4 in the family of incommensurate composite crystals and prompted structural refinement in a (3 + 1)-dimensional model. The structure of Mg 1+ε Co 4 B 4 (P4 2 /ncm(00γ)s00s, a = 6.75847(7) Å, c = 3.94007(8) Å, q = (0, 0, 1.2721(3))) was refined from neutron powder diffraction and high-resolution powder X-ray diffraction data and confirmed by scanning transmission electron microscopy and electron diffraction. Mg 1+ε Co 4 B 4 is isostructural to Nd 1+ε Fe 4 B 4 and several related ternary borides with 0.07 ≤ ε ≤ 0.17, with Mg occupying the rare-earth site. Satellite reflections in the electron diffraction patterns hinted at positional modulation of the transition metal–boron substructure by Mg atoms, but this could not be refined from the neutron or X-ray diffraction data. Low-temperature magnetic measurements show no indications of long-range magnetic ordering or superconductivity down to 5 K. DFT calculations confirmed the absence of a magnetically ordered ground state and the stability of a 5:4 supercell (ε = 0.25) relative to the fully commensurate structure. Neutron diffraction and synthesis from elemental Mg demonstrated that Mg 1+ε Co 4 B 4 is not a hydrogen-stabilized phase. Mg 1+ε Co 4 B 4 represents the second compound reported in the Mg–Co–B system and the first superspace symmetry model of a Nd 1+ε Fe 4 B 4 -type incommensurate composite compound refined from powder diffraction data.

chemical structure↗
DOE OSTI2025

Multiscale Modeling Framework Using Element‐Based Galerkin Methods for Moist Atmospheric Limited‐Area Simulations

This paper presents a multiscale modeling framework (MMF) to model moist atmospheric limited-area weather. The MMF resolves large-scale convection using a coarse grid while simultaneously resolving local features through numerous fine local grids and coupling them seamlessly. Both large- and small-scale processes are modeled using the compressible Navier-Stokes equations within the Nonhydrostatic Unified Model of the Atmosphere (NUMA), and are discretized using a continuous element-based Galerkin method (spectral elements) with high-order basis functions. Consequently, the large-scale and small-scale models share the same dynamical core but have the flexibility to be adjusted individually. The proposed MMF method is tested in 2D and 3D idealized limited-area weather problems involving storm clouds produced by squall line and supercell simulations. Numerical results from the MMF showed enhanced representation of cloud processes compared to the coarse model.

Kang, Soonpil [Naval Postgraduate School, Monterey↗
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↗