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

An atomic cluster expansion potential for twisted multilayer graphene

Twisted multilayer graphene, characterized by its moiré patterns arising from inter-layer rotational misalignment, serves as a rich platform for exploring quantum phenomena. Machine learning interatomic potentials (MLIPs) are a promising approach to model such systems. Our work develops a method to generate training and test datasets for fitting MLIPs that capture all possible misalignments but remain small-scale to facilitate efficient data generation and parameter estimation. To achieve this, we generate configurations with periodic boundary conditions suitable for density functional theory calculations, and then introduce an internal twist and shift within those supercell structures. Using this technique, supplemented with an active learning workflow, we fit an Atomic Cluster Expansion potential for simulating twisted multilayer graphene and test it for accuracy and robustness on a range of simulation tasks.

2D materials↗
DOE OSTI2025

Carrier Localization and Spontaneous Formation of Two-Dimensional Polarization Domain in Halide Perovskites

Halide perovskites are known for their rich phase diagram and superior performance in diverse optoelectronics applications. The latter property is often attributed to the long electron-hole recombination time, whose underlying physical mechanism has been a long-standing controversy. In this Letter, we investigate the transport and localization properties of electron and hole carriers in a prototypical halide perovskite (CsPbBr 3 ), through ab initio tight-binding nonadiabatic dynamics approach for large-scale (tens of nm size) supercell calculations, to simulate electron and ion dynamics on the same footing. Here, we found distinct structural, lattice polarization, and electron-phonon coupling properties at low (below 100 K) and high temperatures, consistent with experimental observations. In particular, at low temperature we find spontaneous formation of polar grain boundaries in the nonpolar bulk systems, which result in two-dimensional polarization patterns that serve to localize and separate electrons and holes. We reveal phonon-assisted variable-range hopping mostly responsible for low-temperature transport, and their characteristic frequency correlates with temperature-dependent phonon power spectrum and energy oscillation frequency in nonadiabatic dynamics. We answer the critical questions of long electron-hole recombination lifetime at low temperature and offer the correlation among polarization domains, electron-phonon couplings, and photocarrier dynamics.

Carrier generation & recombination↗
DOE OSTI2026

Lattice vacancy migration barriers in Fe-Ni alloys, and an indication as to why Ni atoms diffuse slowly: A first-principles study

Lattice vacancy migration barriers in ferromagnetic Fe 𝑥 ⁢Ni 1−𝑥 alloys (0.4 ≤ 𝑥 ≤ 0.6) are accurately quantified within the framework of ab initio electronic structure calculations using the nudged elastic band (NEB) method. Both the atomically disordered (A1) fcc phase, as well as the atomically ordered, tetragonal L⁢1 0 phase—which is under consideration as a material for a rare-earth-free gap magnet for advanced engineering applications—are investigated. Across an ensemble of NEB calculations performed on supercell configurations spanning a range of compositions and containing disordered, partially ordered, and fully ordered structures, we find that Ni-vacancy interchanges encounter significantly higher energetic barriers than do Fe-vacancy interchanges. We contend that this aspect is a key factor in determining the differences in mobility between Fe and Ni atoms in this ferromagnetic alloy. Moreover, we are able to interpret these findings in terms of the ferromagnetic alloy's underlying spin-polarized electronic structure. Specifically, we report a coupling between the size of local lattice distortions and the magnitude of the local electronic spin polarization around vacancies. This causes Fe atoms to relax into lattice vacancies, while Ni atoms remain rigidly fixed to their original lattice positions. These results give atomic-scale insight into the longstanding experimental observation that Ni exhibits remarkably slow atomic diffusion in Fe-Ni alloys.

density functional theory↗
DOE OSTI2025

Quantitative surface studies of moiré intercalated Pb structures under Gr/SiC

Quantitative growth studies using surface diffraction and scanning tunneling microscopy are carried out to determine the intercalation conditions for Pb under Gr/SiC, with initial phase a mixture of buffer-layer and single-layer graphene. Intercalation is a kinetically complex, multistage process with intercalated Pb moving between different subsurface locations with annealing, eventually to the location of higher stability. The transfer is determined by monitoring the difference of a given spot intensity from its clean value, as a function of growth parameters. In addition, a broad, low-intensity (10×10) gr moiré phase is observed exclusively originating from intercalated Pb. Well-defined growth conditions are identified for the moiré spot to be observed reversibly. The (10×10) gr phase follows changes of the Pb(10) spot during growth and annealing experiments, instead of changes of the other Gr/SiC spots indicating preferred nucleation of Pb islands on the (10×10) gr domains. These structural results are discussed in connection to different models proposed for the (10×10) gr supercell and the subsurface Pb coverage. In conclusion, mass transport of subsurface Pb is unusually high, possibly because of the very high compression of the intercalated Pb, when compared to metallic Pb(111).

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

First-principles computations of the Stark shift of a defect-bound exciton: The case of the T center in silicon

The T center in silicon has recently drawn a lot of attention for its potential in quantum information science. The sensitivity of the zero-phonon line (ZPL) to electrical field was recently investigated by a combination of different experimental methods but there are still few first principles studies on the Stark shift of the T center. Dealing with the defect-bound exciton nature of the excited state is particularly challenging using density functional theory because of the large spatial delocalization associated with the wavefunction. Here, in this work, we tackle this issue by performing a convergence study over the supercell size. We obtain an exciton binding energy of 28.5 meV, in good agreement with experimental results. We then calculate the Stark shift through the dipole moment change of the ZPL transition of the T center using the modern theory of polarization formalism and find a modest linear coefficient of Δ⁢𝜇=0.79⁢D along X and Δ⁢𝜇=0.03⁢𝐷 along Y. We discuss our results in light of the recent experimental measurements of the Stark shift. Our analysis suggests that bound-exciton defects could be particularly sensitive to local field effect as a result of their large spatial extent.

color centers↗
DOE OSTI2025

Evolution of charge correlations in the hole-doped kagome superconductor CsV 3−𝑥⁢ Ti 𝑥 ⁢Sb 5

The interplay between superconductivity and charge correlations in the kagome metal CsV 3⁢ Sb 5 can be tuned by external perturbations such as doping or pressure. Here we present a study of charge correlations and superconductivity upon hole doping via Ti substitution on the V kagome sites in CsV 3−𝑥 ⁢Ti 𝑥 ⁢Sb 5 via synchrotron x-ray diffraction and scanning superconducting quantum interference device measurements. While the superconducting phase, as viewed via the vortex structure, remains conventional and unchanged across the phase diagram, the nature of charge correlations evolves as a function of hole doping from the first superconducting dome into the second superconducting dome. For Ti doping in the first superconducting dome, competing 2×2×2 and 2×2×4 supercells form within the charge density wave state and are suppressed rapidly with carrier substitution. In the second superconducting dome, no charge correlations are detected. Here, comparing these results to those observed for CsV 3 ⁢Sb 5−𝑥 ⁢Sn 𝑥 suggests important differences between hole doping via chemical substitution on the V and Sb sites, particularly in the disorder potential associated with each dopant.

Charge density waves↗
DOE OSTI2026

Structural modulation, physical properties, and electronic band structure of the kagome metal UCr 6 ⁢Ge 6

The chemical flexibility of the 𝑅⁢𝑀 6 ⁢𝑋 6 stoichiometry, where an 𝑓-block element is intercalated in the CoSn structure type, allows for the tuning of flatbands associated with kagome lattices to the Fermi level and for emergent phenomena due to interactions between the 𝑓- and 𝑑-electron lattices. Yet, 5⁢𝑓 members of the “166” compounds are underrepresented compared with 4⁢𝑓 members. In this work, we report single-crystal growth of UCr 6 ⁢Ge 6 , which crystallizes in a monoclinically distorted Y 0.5⁢ Co 3 ⁢Ge 3 -type structure. The real-space character of the modulation, which is unique within the 𝑅⁢𝑀 6 ⁢𝑋 6 family, is approximated by a 3×1×2 supercell of the average monoclinic cell. The compound has kagome-lattice flatbands near the Fermi level and a moderately enhanced electronic heat capacity, as evidenced by its low-temperature Sommerfeld coefficient (𝛾=86.5 mJ mol −1 K −2 ) paired with band structure calculations. The small, isotropic magnetization and featureless resistivity of UCr 6 ⁢Ge 6 suggest itinerant uranium 5⁢𝑓 electrons and Pauli paramagnetism. Angle-resolved photoemission spectroscopy results provide evidence for uranium 5⁢𝑓 weight at the Fermi level and for a flatband near the Fermi level associated with the chromium 3⁢𝑑 kagome lattice. The isotropic magnetic behavior of the uranium 5⁢𝑓 electrons starkly contrasts with localized behavior in other uranium 166 compounds, highlighting the high tunability of the magnetic ground state across the material family.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Chiral Topological Surface States on a Finite Square Photonic Crystal Bounded by Air

Chiral, topologically protected, photonic surface states can be found at the boundary between gyrotropic photonic crystals where a changing magnetic field induces different topology across the interface. Typically, photonic crystals with either a suitable band structure on both sides of the interface to provide a band gap and evanescent decay of the surface states away from the interface, or an outer layer with engineered material properties is required. In this paper, we show the emergence of topological, unidirectional surface states at the termination of finite gyrotropic photonic crystals with a simple square lattice and C 4 rotational symmetry bounded by a vacuum, eliminating the need for an outside layer to enable chiral surface modes. Here, we start from an infinite, time-reversal-symmetry-breaking photonic crystal with a band gap associated with bands with nonzero Chern numbers, different from all-zero Chern numbers in air. We then modify the photonic crystal to move this band gap below the light line, while maintaining the Chern-number discontinuities. Band-structure calculations for a supercell approximating a photonic crystal finite in the direction normal to the surface demonstrate the existence, dispersion, and chirality of the surface mode. Extensive direct scattering calculations for a point source and spatial Fourier analysis further reveal a unidirectional free-space topological surface state, which propagates counterclockwise around the surface of a finite photonic crystal, providing a nearly foolproof way to cross-check the surface-mode band structure unaffected by backscattering from local defects. Additionally, scattering simulations allow an independent characterization of the state dispersion and unveil the robustness of the topological plasmonic mode propagation around the 90° bends of the structure, being due to only radiation leakage. In contrast to buried topological surface states, the observed surface modes at the photonic crystal–air interface have the advantage of being accessible to the outside world, allowing one to take advantage of the defect-tolerant backscattering-free surface modes to engineer emission from photonic crystal surfaces into arbitrary free-space beam shapes and directions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Effects of lithium intercalation in twisted bilayer graphene

In this work, we investigate the effects of lithium intercalation in twisted bilayers of graphene, using first-principles electronic structure calculations. To model this system we employ commensurate supercells that correspond to twist angles of 7.34° and 2.45°. From the energetics of lithium absorption we demonstrate that for low Li concentration the intercalants cluster in the AA regions with double the density of a uniform distribution. The charge donated by the Li atoms to the graphene layers results in modifications to the band structure that can be qualitatively captured using a continuum model with modified interlayer couplings in a region of parameter space that has yet to be explored either experimentally or theoretically. Thus, the combination of intercalation and twisted layers simultaneously provides the means for spatial control over material properties and an additional knob with which to tune moiré physics in twisted bilayers of graphene, with potential applications ranging from energy storage and conversion to quantum information.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Electronic band gaps from quantum Monte Carlo methods

We develop a method for calculating the fundamental electronic gap of semiconductors and insulators using grand canonical quantum Monte Carlo simulations. Further, we discuss the origin of the bias introduced by supercell calculations of finite size and show how to correct the leading and subleading finite size errors either based on observables accessible in the finite-sized simulations or from density-functional theory calculations. Our procedure is applied to carbon, silicon, and molecular hydrogen crystals, and compared to experiment for carbon and silicon.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Temperature dependence of electron-phonon interactions in vanadium

First-principles calculations were used to study the Fermi surface of body-centered cubic vanadium at elevated temperatures. Supercell calculations accounted for effects of thermal atom displacements on band energies, and band unfolding was used to project the spectral weight of the electron states into the Brillouin zone of a standard bcc unit cell. An electronic topological transition (ETT, or Lifshitz transition) occurred near the Γ point with increasing temperature, but the large thermal smearings from the atomic disorder and the Fermi-Dirac distribution reduced the effect of this ETT on the electron-phonon interactions. The phonon dispersions showed thermal stiffening of their Kohn anomalies near the Γ point and of the longitudinal N phonon mode. Finally, in general the effects of the ETT were overcome by the thermal smearing of the Fermi surface that reduces the spanning vector densities for anomalous phonon modes.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Many-body electronic structure of LaScO 3 by real-space quantum Monte Carlo

We present real space quantum Monte Carlo (QMC) calculations of the scandate LaSc0 3 that proved to be challenging for traditional electronic structure approaches due to strong correlation effects resulting in inaccurate band gaps from DFT and GW methods when compared with existing experimental data. Besides calculating an accurate QMC band gap corrected for supercell size biases and in agreement with numerous experiments, we also predict the cohesive energy of the crystal using the standard fixed-node QMC without any empirical or non-variational parameters. We show that promotion (optical) gap and fundamental gap agree with each other illustrating a clear absence of significant excitonic effects in the ideal crystal. We obtained these results in perfect consistency in two independent tracks that employ different basis sets (plane wave vs. localized gaussians), different codes for generating orbitals (QUANTUM ESPRESSO vs. CRYSTAL), different QMC codes (QMCPACK vs. QWALK) and different high-accuracy pseudopotentials (ccECPs vs. Troullier-Martins) presenting the maturity and consistency of QMC methodology and tools for studies of strongly correlated problems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Ab initio electronic density in solids by many-body plane-wave auxiliary-field quantum Monte Carlo calculations

We present accurate many-body results of the electronic densities in several solid materials, including Si, NaCl, and Cu. These results are obtained using the ab initio auxiliary-field quantum Monte Carlo (AFQMC) method working in a plane-wave basis with norm-conserving, multiple-projector pseudopotentials. AFQMC has been shown to be an excellent many-body total energy method. Computation of observables and correlation functions other than the ground-state energy requires back-propagation, whose adaption and implementation in the plane-wave basis AFQMC framework are discussed in the present paper. This development allows us to compute correlation functions, electronic densities, and interatomic forces, paving the way for geometry optimizations and calculations of thermodynamic properties in solids. Finite supercell size effects are considerably more subtle in the many-body framework than in independent-electron calculations. We analyze the convergence of the electronic density, and obtain best estimates for the thermodynamic limit. The densities from several typical density functionals are benchmarked against our near-exact results. The electronic densities we obtained can also be used to help construct improved density functionals.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Cohesion and excitations of diamond-structure silicon by quantum Monte Carlo: Benchmarks and control of systematic biases

We have carried out quantum Monte Carlo (QMC) calculations of silicon crystal focusing on the accuracy and systematic biases that affect the electronic structure characteristics. The results show that 64 and 216 atom supercells provide an excellent consistency for extrapolated energies per atom in the thermodynamic limit for ground, excited, and ionized states. We have calculated the ground state cohesion energy with both systematic and statistical errors below ≈ 0.05 eV. The ground state exhibits a fixed-node error of only 1.3(2)% of the correlation energy, suggesting an unusually high accuracy of the corresponding single-reference trial wave function. We obtain a very good agreement between optical and quasiparticle gaps that affirms the marginal impact of excitonic effects. Our most accurate results for band gaps differ from the experiments by about 0.2 eV. This difference is assigned to a combination of residual finite-size and fixed-node errors. In conclusion, we have estimated the crystal Fermi level referenced to vacuum that enabled us to calculate the edges of valence and conduction bands in agreement with experiments.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Enhanced orbital anisotropy through the proximity to a SrTiO 3 layer in the perovskite iridate superlattices

We have used angle-dependent soft x-ray absorption spectroscopy (XAS) at the O K edge and first-principles calculations to investigate the electronic structures of iridate-based superlattices ( SrIr O 3 ) m / ( SrTi O 3 ) ( m = 1 , 2, 3, and ∞). We focus on the pre-edge Ir 5 d t 2 g – O 2 p orbital hybridization feature in the XAS spectra. By varying the measurement geometry relative to the incident photon polarization, we are able to extract the dichroic contrast and observe the systematic increase in the anisotropy of Ir 5 d orbitals as m decreases. First-principles calculations elucidate the orbital anisotropy coming mainly from the enhanced out-of-plane compression of Ir O 6 octahedra in the SrIr O 3 layers that are adjacent to the inserted SrTi O 3 layers. As m decreases, the increased volume fraction of these interfacial SrIr O 3 layers and their contact with the SrTi O 3 layers within the ( SrIr O 3 ) m / ( SrTi O 3 ) supercell lead to enhanced orbital anisotropy. Furthermore, the tilt and rotation of Ir O 6 octahedra are shown to be essential to understand the subtle orbital anisotropy in these superlattices, and constraining these degrees of freedom will give an incorrect trend. Overall, our results demonstrate that the structural constraint from the inserted SrTi O 3 layers, in addition to other electronic means such as polar interface and charge transfer, can serve as a knob to control the orbital degree of freedom in these iridate-based superlattices.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Fast and robust all-electron density functional theory calculations in solids using orthogonalized enriched finite elements

Here, we present a computationally efficient approach to perform systematically convergent real-space all-electron Kohn-Sham density functional theory calculations for solids using an enriched finite element (FE) basis. The enriched FE basis is constructed by augmenting the classical FE basis with atom-centered numerical basis functions, comprising of atomic solutions to the Kohn-Sham problem. Notably, to improve the conditioning, we orthogonalize the enrichment functions with respect to the classical FE basis, without sacrificing the locality of the resultant basis. In addition to improved conditioning, this orthogonalization procedure also renders the overlap matrix block diagonal, greatly simplifying its inversion. Subsequently, we use a Chebyshev polynomial based filtering technique to efficiently compute the occupied eigenspace in each self-consistent field iteration. We demonstrate the accuracy and efficiency of the proposed approach on periodic unit cells and supercells. The benchmark studies show a staggering 130× speedup of the orthogonalized enriched FE basis over the classical FE basis. We also present a comparison of the orthogonalized enriched FE basis with the linearized augmented plane-wave + local orbitals basis, both in terms of accuracy and efficiency. Notably, we demonstrate that the orthogonalized enriched FE basis can handle large system sizes of ~10 000 electrons. Finally, we observe good parallel scalability of our implementation with 92% efficiency at 22× speedup for a system with 620 electrons

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Driven responses of periodically patterned superconducting films

Here, we simulate the motion of a commensurate vortex lattice in a periodic lattice of artificial circular pinning sites having different diameters, pinning strengths, and spacings using the time-dependent Ginzburg-Landau formalism. Above some critical DC current density J c , the vortices depin, and the resulting steady-state motion then induces an oscillatory electric field E (t) with a defect "hopping" frequency f 0 , which depends on the applied current density and the pinning landscape characteristics. The frequency generated can be locked to an applied AC current density over some range of frequencies, which depends on the amplitude of the DC as well as the AC current densities. Both synchronous and asynchronous collective hopping behaviors are studied as a function of the supercell size of the simulated system and the (asymptotic) synchronization threshold current densities determined.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Optical conductivity of overdoped cuprates from ab initio out-of-plane impurity potentials

Dopant impurity potentials determined by ab initio supercell density functional theory calculations are used to calculate the optical conductivity of overdoped La 2-x ⁢Sr x CuO 4 (LSCO) and Tl 2 ⁢Ba 2⁢ CuO 6+δ (Tl-2201) in the superconducting and normal states. Vertex corrections are included to account for the effect of forward scattering on two-particle properties. This approach was previously shown to provide good, semiquantitative agreement with measurements of superfluid density in LSCO. Here, in this study, we compare calculations of conductivity with measurements of terahertz conductivity on LSCO using identical impurity, band, and correlation parameters and find similarly good correspondence with experiment. In the process, we delineate the impact of the different disorder mechanisms on single-particle and transport relaxation processes. In particular, we reveal the critical role of apical oxygen vacancies in transport scattering and show that transport relaxation rates in LSCO are significantly reduced when apical oxygen vacancies are annealed out. These considerations are shown to be crucial for understanding the variability of experimental results for overdoped LSCO in samples of nominally identical doping but different types. Finally, we give predictions for Tl-2201 terahertz conductivity experiments.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗