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

Investigating Inherent Numerical Stabilization for the Moist, Compressible, Non‐Hydrostatic Euler Equations on Collocated Grids

This study investigates inherent numerical dissipation due to upwind fluxes and reconstruction strategies for collocated Finite-Volume integration of the Euler equations. Idealized supercell simulations are used without any explicit dissipation. Flux terms are split into: mass flux, pressure, and advected quantities. They are computed with the following upwind strategies: central, advectively upwind, and acoustically upwind. This is performed for third and ninth-order-accurate reconstructions with and without Weighted Essentially Non-Oscillatory limiting. Acoustic-only upwinding for pressure and mass flux terms and advective-only upwinding for advected quantities is the most flexible simplification found. It reduces data movement and computations. Assuming a constant speed of sound in acoustic upwinding gives similar results to using the true speed of sound. Dissipation from upwind adapts automatically to grid spacing, time step, reconstruction accuracy, and flow smoothness. While stability is maintained even at 21st-order spatial accuracy, there is a limit to the spatial order of accuracy for which upwinding alone can create a realizable solution in the conditions of this study. Convex combinations of upwind and central solutions for flux terms also reduced dissipation, but as the central proportion grows, solutions become physically unrealizable. The range of length scales of the kinetic energy spectra can be extended along k −5/3 to smaller spatial scales by reducing dissipation either with higher-order reconstructions or using convex combinations of upwind and central fluxes. However, not all extensions of the length scale range along k −5/3 exhibit physically realizable solutions, even though the spectra appear to be physical.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Paramagnetic electronic structure of CrSBr: Comparison between ab initio G W theory and angle-resolved photoemission spectroscopy

Here, we explore the electronic structure of paramagnetic CrSBr by comparative first-principles calculations and angle-resolved photoemission spectroscopy. We theoretically approximate the paramagnetic phase using a supercell hosting spin configurations with broken long-range order and applying quasiparticle self-consistent GW theory, without and with the inclusion of excitonic vertex corrections to the screened Coulomb interaction (QS GW and QS $G\widehat{W}$, respectively). Comparing the quasiparticle band-structure calculations to angle-resolved photoemission data collected at 200 K results in excellent agreement. This allows us to qualitatively explain the significant broadening of some bands as arising from the broken magnetic long-range order and/or electronic dispersion perpendicular to the quasi-two-dimensional layers of the crystal structure. The experimental band gap at 200 K is found to be at least 1.51 eV at 200 K. At lower temperature, no photoemission data can be collected as a result of charging effects, pointing towards a significantly larger gap, which is consistent with the calculated band gap of approximately 2.1 eV.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

rmc-discord : reverse Monte Carlo refinement of diffuse scattering and correlated disorder from single crystals

A user-friendly program has been developed to analyze diffuse scattering from single crystals with the reverse Monte Carlo method. The approach allows for refinement of correlated disorder from atomistic supercells with magnetic or structural (occupational and/or displacive) disorder. The program is written in Python and optimized for performance and efficiency. Refinements of two user cases obtained with legacy neutron-scattering data demonstrate the effectiveness of the approach and the developed program. It is shown with bixbyite, a naturally occurring magnetic mineral, that the calculated three-dimensional spin-pair correlations are resolved with finer real-space resolution compared with the pair distribution function calculated directly from the reciprocal-space pattern. With the triangular lattice Ba 3 Co 2 O 6 (CO 3 ) 0.7 , refinements of occupational and displacive disorder are combined to extract the one-dimensional intra-chain correlations of carbonate molecules that move toward neighboring vacant sites to accommodate strain induced by electrostatic interactions. The program is packaged with a graphical user interface and extensible to serve the needs of single-crystal diffractometer instruments that collect diffuse-scattering data.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

PyFaults: A Python Toolkit for Stacking Fault Screening

PyFaults is an open-source Python library designed to model stacking fault disorder in crystalline materials and qualitatively assess the characteristic selective broadening effects in powder X-ray diffraction (PXRD). Here, the main capabilities of PyFaults are presented, including unit cell and supercell model construction, PXRD pattern calculation, assessment against experimental PXRD, and methods for rapid screening of candidate models within a set of possible stacking vectors and fault occurrence probabilities. This program aims to serve as a computationally inexpensive tool for identifying and screening potential stacking fault models in materials with planar disorder. Three diverse case studies, involving GaN, Li2MnO3 and Li3YCl6, are presented to illustrate the program functionality across a range of structure types and stacking fault modalities.

MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Urbanization may enhance tornado potential: A single case report

Tornadoes pose a risk of catastrophic economic loss and casualty in the United States. Modification of land use by urbanization alters the meteorological conditions that may impact tornado formation and intensification processes. Here we explored the simulated impact of Kansas City urbanization on the tornado potential of a supercell storm. In this studied case, we found that urbanization might enhance tornado potential by a) strengthening the low-level streamwise vorticity in the storm inflow region, thus forming stronger rotating updrafts; and b) intensifying near-surface horizontal vorticity near the boundary of the forward-flank cold pool which increases the ingestion, tilting, and stretching of streamwise horizontal vorticity into vertical vorticity. The former results from the stronger east-to-west pressure perturbation gradient due to the faster, stronger outflow boundary, and the latter is mainly a result of stronger cold fronts and a better alignment of storm-relative inflow with the horizontal vorticity vector. We emphasize that our conclusions only represent one possibility of how urbanization would affect tornado potential, and a more systematic examination is needed to achieve a more general conclusion.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Phase‐Dependent Band Gap Engineering in Alloys of Metal‐Semiconductor Transition Metal Dichalcogenides

Abstract Bandgap engineering plays a critical role in optimizing the electrical, optical and (photo)‐electrochemical applications of semiconductors. Alloying has been a historically successful way of tuning bandgaps by making solid solutions of two isovalent semiconductors. In this work, a novel form of bandgap engineering involving alloying non‐isovalent cations in a 2D transition metal dichalcogenide (TMDC) is presented. By alloying semiconducting MoSe 2 with metallic NbSe 2 , two structural phases of Mo 0.5 Nb 0.5 Se 2 , the 1T and 2H phases, are produced each with emergent electronic structure. At room temperature, it is observed that the 1T and 2H phases are semiconducting and metallic, respectively. For the 1T structure, scanning tunneling microscopy/spectroscopy (STM/STS) is used to measure band gaps in the range of 0.42–0.58 at 77 K. Electron diffraction patterns of the 1T structure obtained at room temperature show the presence of a nearly commensurate charge density wave (NCCDW) phase with periodic lattice distortions that result in an uncommon 4 × 4 supercell, rotated approximately 4° from the lattice. Density‐functional‐theory calculations confirm that local distortions, such as those in a NCCDW, can open up a band gap in 1T ‐Mo 0.5 Nb 0.5 Se 2 , but not in the 2H phase. This work expands the boundaries of alloy‐based bandgap engineering by introducing a novel technique that facilitates CDW phases through alloying.

Wang, Shuxi↗
DOE OSTI2022

Dative Epitaxy of Commensurate Monocrystalline Covalent van der Waals Moiré Supercrystal

Realizing van der Waals (vdW) epitaxy in the 1980s represents a breakthrough that circumvents the stringent lattice matching and processing compatibility requirements in conventional covalent heteroepitaxy. However, due to the weak vdW interactions, there is little control over film qualities by the substrate. Typically, discrete domains with a spread of misorientation angles are formed, limiting the applicability of vdW epitaxy. In this study, the epitaxial growth of monocrystalline, covalent Cr 5 Te 8 2D crystals on monolayer vdW WSe 2 by chemical vapor deposition is reported, driven by interfacial dative bond formation. The lattice of Cr 5 Te 8 , with a lateral dimension of a few tens of micrometers, is fully commensurate with that of WSe 2 via 3 × 3 (Cr 5 Te 8 )/7 × 7 (WSe 2 ) supercell matching, forming a single-crystalline moiré superlattice. This work establishes a conceptually distinct paradigm of thin-film epitaxy, termed “dative epitaxy”, which takes full advantage of covalent epitaxy with chemical bonding for fixing the atomic registry and crystal orientation, while circumventing its stringent lattice matching and processing compatibility requirements; conversely, it ensures the full flexibility of vdW epitaxy, while avoiding its poor orientation control. Cr 5 Te 8 2D crystals grown by dative epitaxy exhibit square magnetic hysteresis, suggesting minimized interfacial defects that can serve as pinning sites.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Dirac Fermion Cloning, Moiré Flat Bands, and Magic Lattice Constants in Epitaxial Monolayer Graphene

Tuning interactions between Dirac states in graphene has attracted enormous interest because it can modify the electronic spectrum of the 2D material, enhance electron correlations, and give rise to novel condensed-matter phases such as superconductors, Mott insulators, Wigner crystals, and quantum anomalous Hall insulators. Previous works predominantly focus on the flat band dispersion of coupled Dirac states from different twisted graphene layers. Here in this work, a new route to realizing flat band physics in monolayer graphene under a periodic modulation from substrates is proposed. Graphene/SiC heterostructure is taken as a prototypical example and it is demonstrated experimentally that the substrate modulation leads to Dirac fermion cloning and, consequently, the proximity of the two Dirac cones of monolayer graphene in momentum space. Theoretical modeling captures the cloning mechanism of the Dirac states and indicates that moiré flat bands can emerge at certain magic lattice constants of the substrate, specifically when the period of modulation becomes nearly commensurate with the ($\sqrt{3}$ x $\sqrt{3}$) R30° supercell of graphene. Our results show that epitaxial single monolayer graphene on suitable substrates is a promising platform for exploring exotic many-body quantum phases arising from interactions between Dirac electrons.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Insight into Pyroelectricity and Phase Transitions in Ferroelectrics from Nonequilibrium Approach: The Case of PbTiO 3

Computational methodology is proposed to simulate time dependent temperature variations via thermostat approach. Here, the methodology is applied to study pyroelectricity under nonequilibrium conditions in prototypical ferroelectric PbTiO 3 . It is found that at room temperature the nonequilibrium effects in pyroelectricity begin manifesting above 1.0 THz and originate from the overlap with the intrinsic soft mode frequency. Therefore, the proposed computational methodology can be used to predict equilibrium pyroelectric response by choosing a frequency of temperature variation below the intrinsic polar mode frequencies for the material. Probing the dynamics of ferroelectric phase transition with the proposed methodology revealed that it takes about a nanosecond for the computational supercell to undergo phase transition. It is concluded that computational techniques capable of simulating picoseconds can be used to compute pyroelectric response from nonequilibrium dynamical approach while those that can reach nanoseconds are suitable for simulating phase transitions in ferroelectrics.

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

Effect of XC functionals and dispersion corrections on the DFT‐computed structural and vibrational properties of SrCl 2 –NaCl and ZrF 4 –LiF

Density functional theory (DFT) calculations were performed to examine the impact of exchange–correlation (XC) functionals and van der Waals corrections (specifically the D3 method) on the structural and vibrational properties of the SrCl 2 –NaCl and ZrF 4 –LiF salt systems. Multiple XC functionals, including the local density approximation (LDA), the generalized gradient approximation using the Perdew–Burke–Ernzerhof (PBE) model, and its modified form suitable for solids (PBEsol), the dispersion-corrected PBE-D3 and PBEsol-D3, were considered. Of these functionals, LDA was found to exhibit the highest degree of error, while PBEsol and PBE-D3 displayed the least error. Underestimated lattice parameters compared with experimental values were observed to result in higher force constants, leading to an overprediction of vibrational frequencies. Conversely, an overestimation of lattice parameters was associated with lower vibrational frequencies. The methodology presented in this study yielded results that are in good agreement with experiment, irrespective of the method (finite differences vs. density functional perturbation theory) employed for calculating infrared and Raman spectra. It was further demonstrated that for alkali halides with weak Raman scattering, utilizing a supercell constructed from primitive cells better predicts Raman features than does the use of conventional cells.

Raman↗
DOE OSTI2022

Electronic Structure of InAs and InSb Surfaces: Density Functional Theory and Angle-Resolved Photoemission Spectroscopy

The electronic structure of surfaces plays a key role in the properties of quantum devices. However, surfaces are also the most challenging to simulate and engineer. Here, in this work, the electronic structure of InAs(001), InAs(111), and InSb(110) surfaces is studied using a combination of density functional theory (DFT) and angle-resolved photoemission spectroscopy (ARPES). Large-scale first principles simulations are enabled by using DFT calculations with a machine-learned Hubbard U correction [npj Comput. Mater. 6, 180 (2020)]. To facilitate direct comparison with ARPES results, a “bulk unfolding” scheme is implemented by projecting the calculated band structure of a supercell surface slab model onto the bulk primitive cell. For all three surfaces, a good agreement is found between DFT calculations and ARPES. For InAs(001), the simulations clarify the effect of the surface reconstruction. Different reconstructions are found to produce distinctive surface states, which may be detected by ARPES with low photon energies. For InAs(111) and InSb(110), the simulations help elucidate the effect of oxidation. Owing to larger charge transfer from As to O than from Sb to O, oxidation of InAs(111) leads to significant band bending and produces an electron pocket, whereas oxidation of InSb(110) does not. The combined theoretical and experimental results may inform the design of quantum devices based on InAs and InSb semiconductors, for example, topological qubits utilizing the Majorana zero modes.

42 ENGINEERING↗
DOE OSTI2022

Novel self-assembled two-dimensional layered oxide structure incorporated with Au nanoinclusions towards multifunctionalities

Two-dimensional (2D) layered oxides have recently attracted wide attention owing to the strong coupling among charges, spins, lattice, and strain, which allows great flexibility and opportunities in structure designs as well as multifunctionality exploration. In parallel, plasmonic hybrid nanostructures exhibit exotic localized surface plasmon resonance (LSPR) providing a broad range of applications in nanophotonic devices and sensors. A hybrid material platform combining the unique multifunctional 2D layered oxides and plasmonic nanostructures brings optical tuning into the new level. Here, a novel self-assembled Bi 2 MoO 6 (BMO) 2D layered oxide incorporated with plasmonic Au nanoinclusions has been demonstrated via one-step pulsed laser deposition (PLD) technique. Comprehensive microstructural characterizations, including scanning transmission electron microscopy (STEM), differential phase contrast imaging (DPC), and STEM tomography, have demonstrated the high epitaxial quality and particle-in-matrix morphology of the BMO-Au nanocomposite film. DPC-STEM imaging clarifies the magnetic domain structures of BMO matrix. Three different BMO structures including layered supercell (LSC) and superlattices have been revealed which is attributed to the variable strain states throughout the BMO-Au film. Owing to the combination of plasmonic Au and layered structure of BMO, the nanocomposite film exhibits a typical LSPR in visible wavelength region and strong anisotropy in terms of its optical and ferromagnetic properties. This study opens a new avenue for developing novel 2D layered complex oxides incorporated with plasmonic metal or semiconductor phases showing great potential for applications in multifunctional nanoelectronics devices.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Order Parameter Engineering for Random Systems

The chemical short-range order (CSRO) in crystalline materials influences the properties, and its effect is significant in the context of multicomponent materials. Here, we propose a scheme for the CSRO parameter or Δ-parameter in terms of the number of like and unlike bonds in the multicomponent systems. The OPERA or Order Parameter Engineering for RAndom Systems scheme for semi-canonical and canonical ensembles is proposed. The proposed framework of Δ-parameter with OPERA framework can generate the single-phase supercell with desired CSRO without explicit energy calculations and provides a computationally efficient scheme for exploration of the CSRO. We demonstrate the applicability of the Δ-parameter as a scalar quantity for describing the CSRO in multicomponent alloys and oxides (FCC-CoCrNi, BCC-MoNbTaW, and (CoCuMgNiZn)O).

36 MATERIALS SCIENCE↗
DOE OSTI2021

A nodal transport method for coupled fast-thermal reactor analysis

We report that for efficient core analyses of coupled fast-thermal reactors, a new deterministic method has been developed based on the variational nodal transport method of VARIANT. A new multigroup cross section generation procedure was devised by combining Monte Carlo lattice calculations for thermal assemblies and the two-step procedure of the MC 2 -3 code for fast assemblies. To reduce assembly homogenization errors in nodal transport calculations, a new nodal equivalence method was developed based on the partial current discontinuity factor (PCDF) and incorporated in VARIANT. A practical procedure to calculate PCDFs with fixed source supercell calculations was also devised. The performance of the proposed method was investigated using a test problem derived from the versatile coupled test reactor (VCTR) design. The new procedure produced multigroup cross sections accurately for all assemblies. VARIANT transport calculations with PCDFs produced accurate multiplication factor and power distribution compared to reference Serpent-2 Monte Carlo solutions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Coverage-dependent structures and thermodynamic stability of intercalated Gd layers beneath buffer-layer graphene on SiC(0001)

Electronic properties of two-dimensional (2D) materials are strongly influenced by their atomic arrangements, making the theoretically-aided characterization of experimentally-synthesized 2D structures crucial. Using first-principles density functional theory, we analyze nearly 200 configurations of intercalated Gd layers beneath buffer-layer graphene on SiC(0001) over a Gd coverage range of 0.01 < θ < 1.2. By fully relaxing selectively-constructed configurations at each coverage within a large, low-strain supercell, we determine the coverage dependence of the chemical potential for intercalated Gd structures. Thermodynamically-preferred configurations below θ ≈ 0.8 form single-atom-thick monolayers, while 3D-like or multilayer structures emerge beyond θ ≈ 0.9. Most structures are amorphous-like, including the configuration at the chemical potential minimum around θ ≈ 0.4. In contrast, a strongly stretched Gd(0001)-like monolayer at θ = 1/3 and a nearly perfect Gd(0001) monolayer at θ = 1 are significantly less favorable with 0.16 eV and 0.82 eV higher chemical potentials above the minimum, respectively. Furthermore, the graphene layer decoupled by intercalated Gd near the chemical potential minimum is significantly flatter compared to its morphology above intercalated 3D structures at higher coverages and nearly isolated Gd atoms in the lowest coverage region. In conclusion, these findings align with our experimental results and underscore the need for further research on this unique intercalated system, which holds significant potential for diverse applications.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Rapid generation of optimal generalized Monkhorst-Pack grids

In this report computational modeling of the properties of crystalline materials has become an increasingly important aspect of materials research, consuming hundreds of millions of CPU-hours at scientific computing centers around the world each year, if not more. A routine operation in such calculations is the evaluation of integrals over the Brillouin zone. We have previously demonstrated that performing such integrals using generalized Monkhorst-Pack k-point grids can roughly double the speed of these calculations relative to the widely-used traditional Monkhorst-Pack grids. However the generation of optimal generalized Monkhorst-Pack grids is not implemented in most software packages due to the computational cost and difficulty of identifying the best grids. To address this problem, we present new algorithms that allow rapid generation of optimal generalized Monkhorst-Pack grids on the fly. We demonstrate that the grids generated by these algorithms are on average significantly more efficient than those generated using existing algorithms across a range of grid densities. For grids that correspond to a real-space supercell with at least 50 Å between lattice points, which is sufficient to converge density functional theory calculations within 1 meV/atom for nearly all materials, our algorithm finds optimized grids in an average of 0.19 s on a single processing core. To facilitate the widespread adoption of this approach, we present new open-source tools including a library designed for integration with third-party software packages.

36 MATERIALS SCIENCE↗
DOE OSTI2021

An efficient computational framework for charge density estimation in twisted bilayer graphene

Electronic properties such as band structure and Fermi velocity in low-angle twisted bilayer graphene (TBG) are intrinsically dependent on the atomic structure. Rigid rotation between individual graphene layers provides an approximate description of the bilayer symmetry. Upon relaxation, in-plane displacement of the atoms in low angle TBG causes a change in the symmetry through the enlargement of the AB stacking regions and the reduction in size of AA and SP stacking regions. However, the effect of this in-plane relaxation on the charge density remains unexplored, because the necessary electronic structure calculations of such large supercells of low twist angle TBG are computationally infeasible. Therefore, we develop a computationally efficient framework that enables the exploration of the charge density symmetry of the low twist angle TBG. This framework is based on the Fourier representation of the charge density which presents high intensity Bragg peaks. Here we find that with the decrease of twist angle, low intensity satellite peaks also become apparent. Our framework incorporates these satellite peaks which reveals transformation of symmetry in the charge density distribution from high to low twist angle TBG. One striking outcome is the demonstration of the electron localization in the AA region of low twist angle TBG. Our framework helps to explain the effect of the atomistic relaxation on the charge density distribution and thus, it provides information about exotic electronic properties of low twist angle TBG at a low computational expense.

36 MATERIALS SCIENCE↗
DOE OSTI2022

First-principles study of Al/Al 3 Ni interfaces

Al-Ni alloys have shown promise for high-temperature applications due to the strengthening of Al 3 Ni fibers resistant to coarsening and spheroidization up to 400°C. While the interface between Al and Al 3 Ni phases affects the coarsening rate of Al 3 Ni at elevated temperatures, its characteristics are largely unknown to date. Here, we have constructed various supercells to model this interface and performed a first-principles study based on density functional theory (DFT). We have considered three groups of Al/Al 3 Ni interfaces: experimentally reported orientation relationships from the solidification studies, crystallographically similar Fe-Fe 3 C pearlite interfaces, and the family of low-index (100) termination planes. We have analyzed the correlation between the DFT Al/Al 3 Ni interfacial energies and characteristic features, e.g., excess free volume and the number of broken bonds. We outline the further experimental and computational analysis required to improve the interface modeling of Al/Al 3 Ni.

36 MATERIALS SCIENCE↗