Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “polymorphous configurations”

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.

14 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Stacking polymorphism of PtSe 2 : its implication to layer-dependent metal-insulator transitions

Using diffusion Monte Carlo (DMC) and density functional theory (DFT) calculations, we examine the structural stability and interlayer binding properties of PtSe 2 . Our DMC study reveals that AA and AB-r bilayer stacking modes are nearly degenerate, highlighting the significant role of interlayer hybridization in offsetting the energy cost due to larger interlayer separations in the AB-r mode. Additionally, our DMC-benchmarked DFT calculations with the r 2 SCAN+rVV10 functional uncover pronounced stacking polymorphism in few-layer PtSe 2 , driven by degenerate AA and AB-r interfaces, which leads to substantial band gap variations across different stacking configurations. This polymorphism, along with selenium vacancies, influences a layer-dependent metal-insulator transition observed in few-layer PtSe 2 . Our findings emphasize the importance of both van der Waals interactions and interlayer hybridization in determining the phase stability and electronic properties of PtSe 2 , advancing our understanding of its fundamental properties and refining theoretical models for practical applications in nanoelectronic devices.

Electronic materials↗
DOE OSTI2026

Accelerated Selectrion of Optimal Perovskite Alloys for Solar PV using a Combined Quantum and Machine Learning Hierachiral Approach

The project aims to: (i) accelerate the discovery of “Missing HP alloys” by combining quantum mechanics and artificial intelligence machine learning approaches, and (ii) analyze the stabilities of candidate alloys, including those that do not pass selection filters (and are thus expected to degrade over time) to decipher the nature of the instabilities to guide the development of durable solar cell materials. Successful candidates will be subjected to validation experiments at NREL's state-of-the-art facilities. The discoveries this effort will provide will be directly testable and implementable and will greatly impact U.S. progress in HP PV as they will provide clear direction and motivation for experimental studies including specific material synthetic targets, device optimization, and device stability protocols. A key advantage of this effort is the feedback and guidance provided by the Industry Collaborative Work Group that we established to coordinate academic and national lab research with industry needs. The proposed work will provide a basis for and direct the development of robust and reliable HP PV. It will also provide a timely, valuable and extensive roadmap to the experimental HP PV community to enable it to focus its efforts on improving and fine-tuning promising HP compositions that this effort predicts will likely be the best performers rather than wandering in the vast chemical space for decades spending enormous resources mostly evaluating unpromising candidate materials.

14 SOLAR ENERGY↗
DOE OSTI2025

Damage buildup in multilayer hexagonal boron nitride films under Ar ion bombardment

Ion irradiation can be used to modify properties of two-dimensional hexagonal boron nitride (2D-hBN). However, the behavior of multilayer 2D-hBN under ion bombardment remains essentially unexplored. Here, in this study, we systematically study the damage buildup in 40-nm-thick multilayer 2D-hBN films irradiated at room temperature with 500 keV Ar ions with doses up to 0.64 displacements per atom. Transmission electron microscopy, soft X-ray absorption spectroscopy, and Raman spectroscopy are employed to study effects of ion irradiation. Results reveal a monotonic buildup of damage that can be described by the phenomenological zero cascade overlap model. Compared to the behavior of other group-III nitride semiconductors, radiation damage in hBN is influenced by polymorphism, involving transitions between bonding configurations of hexagonal, cubic, and amorphous phases.

Seo, M. [Lawrence Livermore National Laboratory (L↗
DOE OSTI2025

Comparison of nested geometry treatments within GPU-based Monte Carlo neutron transport simulations of fission reactors

Monte Carlo (MC) neutron transport provides detailed estimates of radiological quantities within fission reactors. This involves tracking individual neutrons through a computational geometry. CPU-based MC codes use multiple polymorphic tracker types with different tracking algorithms to exploit the repeated configurations of reactors, but virtual function calls have high overhead on the GPU. The Shift MC code was modified to support GPU-based tracking with three strategies: dynamic polymorphism with virtual functions, static polymorphism, and a single tracker type with tree-based acceleration. On the Frontier supercomputer these methods achieve 77.8%, 91.2%, and 83.4%, respectively, of the tracking rate obtained using a specialized tracker optimized for rectilinear-grid-based reactors. This indicates that all three methods are suitable for typical reactor problems in which tracking does not dominate runtime. The flexibility of the single tracker method is highlighted with a hexagonal-grid microreactor problem, performed without hexagonal-grid-specific tracking routines, providing a 2.19× speedup over CPU execution.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Neural network potentials with effective charge separation for non-equilibrium dynamics of ionic solids: a ZnO case study

Developing neural network potentials (NNPs) accurate under non-equilibrium dynamics is challenging, as such systems require extensive sampling beyond equilibrium phases. Here we construct high-fidelity NNPs for zinc oxide (ZnO), a polymorphic ionic solid, using density functional theory (DFT) reference data. To efficiently capture transitional configurations, we combine enhanced-sampling molecular dynamics with empirical potentials, data distillation, and pretraining on short-range atomic energies (A-Train), followed by transfer learning with DFT-relabeled datasets. This hierarchical approach improves transferability across polymorphs and stress states. We further introduce effective charge separation, treating long-range Coulombic terms analytically while short-range residual interactions are learned by the NNP. The optimal effective charges fall in the range 0.5–1.0 q e , consistent with dielectric-screened values derived from formal charges but distinct from Bader estimates. Motivated by this observation, we propose a simple data-driven protocol in which effective charges are optimized by comparing DFT reference energies with explicit Coulomb calculations, without additional NNP training. This strategy improves accuracy and transferability in DFT-level predictions of energies, forces, and stress. Together, these results provide a practical charge-selection framework for robust NNP development in ionic solids, enabling reliable simulation of polymorphic phase transformations and non-equilibrium dynamics.

Chemistry↗
DOE OSTI2024

Coherent-Precipitation-Stabilized Phase Formation in Over-Stoichiometric Rocksalt-Type Li Superionic Conductors

Rationalizing synthetic pathways is crucial for material design and property optimization, especially for polymorphic and metastable phases. Over-stoichiometric rocksalt (ORX) compounds, characterized by their face-sharing configurations, are a promising group of materials with unique properties; however, their development is significantly hindered by challenges in synthesizability. Here, taking the recently identified Li superionic conductor, over-stoichiometric rocksalt Li–In–Sn–O (o-LISO) material as a prototypical ORX compound, the mechanisms of phase formation are systematically investigated. It is revealed that the spinel-like phase with unconventional stoichiometry forms as coherent precipitate from the high-temperature-stabilized cation-disordered rocksalt phase upon fast cooling. This process prevents direct phase decomposition and kinetically locks the system in a metastable state with the desired face-sharing Li configurations. This insight enables us to enhance the ionic conductivity of o-LISO to be >1 mS cm -1 at room temperature through low-temperature post-annealing. This work offers insights into the synthesis of ORX materials and highlights important opportunities in this new class of materials.

36 MATERIALS SCIENCE↗
DOE OSTI2025

LeaPP: Learning Pathways to Polymorphs through Machine Learning Analysis of Atomic Trajectories

Understanding the mechanisms underlying crystal nucleation and growth is crucial for many technological applications. Due to the short length and time scales involved, crystal nucleation is often studied using molecular simulations. Most existing approaches to extract the nucleation mechanism from simulations focus on the analysis of static snapshots of the configurations, potentially overlooking subtle local fluctuations and the history of the particles involved in the formation of solid nuclei. Here, in this work, we propose a novel methodology called LeaPP that categorizes nucleation trajectories based on the temporal information of their constituent particles. We leverage the time evolution of the local environment of the crystallizing particles to encapsulate the relationship between the structure and dynamics and distinguish between different evolving particle paths. Identification of the distinct particle paths further enables characterizing the nucleation trajectories into different pathways. Collectively, LeaPP provides a more nuanced understanding of nucleation through an unsupervised approach with lesser dependence on traditional order parameters. Furthermore, the pathways identified by LeaPP are predictive of the resulting polymorph. We demonstrate LeaPP on three different systems─Lennard-Jones-like particles, Ni 3 Al, and water on surfaces. The general methodology underlying LeaPP─considering the time evolution of the building blocks─applies to a wide range of self-assembly problems.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Programmable Phase Selection between Altermagnetic and Noncentrosymmetric Polymorphs of MnTe on InP via Molecular Beam Epitaxy

This dataset contains DFT input and output files supporting the theoretical modeling in the associated publication (ACS Appl. Mater. Interfaces 2026, 18, 15654-15664). The calculations model the interfacial energetics of two MnTe polymorphs — NiAs-MnTe (hexagonal, alpha phase) and ZnS-MnTe (cubic, gamma phase) — on InP(111) substrates with two surface terminations: In-terminated InP(111)A and P-terminated InP(111)B. This gives four interface configurations: NiAs on In-terminated (experimentally observed), NiAs on P-terminated (computed for comparison), ZnS on In-terminated (computed for comparison), and ZnS on P-terminated (experimentally observed). The dataset is organized into four calculation types, each covering all four polymorph/termination combinations: (i) Slabs: Pristine MnTe/InP heterostructure slabs used to compute total energies and interface energy densities (Eint) for all four configurations, as reported in Fig. 6 of the main text. (ii) Disorder: Same slab geometries with a P_Te + Te_P antisite defect pair introduced near the interface, used to assess chemical intermixing effects on interface stability (Fig. S8, SI). (iii) Strain: Pristine slab calculations with in-plane lattice parameters strained by -1% and +1% relative to the InP lattice constant, used to evaluate strain-dependent interface energetics (Fig. S9, SI). (iv) Charge_Density: Single-point calculations on the full heterostructure, the isolated InP slab, and the isolated MnTe slab at fixed geometry, used to compute differential charge density plots showing interfacial charge accumulation and depletion as a function of surface termination (Fig. S10, SI). Each calculation folder contains INCAR, KPOINTS, POSCAR, CONTCAR, OUTCAR, and POTCAR_info.txt (PAW potential information, excluding the full POTCAR due to VASP licensing restrictions). The calculations were performed using VASP 6.4.3 with PBE exchange-correlation, PAW potentials, a Hubbard correction of Ueff = 5 eV on Mn d-states, and A-type AFM spin initialization.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Multi-scale structural analysis of swift heavy ion-irradiated ZrO 2 and HfO 2

The radiation-induced monoclinic-to-tetragonal phase transition in ZrO 2 and HfO 2 has been the subject of many investigations, but the transformation pathways and underlying structural mechanisms are still not well understood. In this study, microcrystalline powder samples of ZrO 2 and HfO 2 were irradiated with 946 MeV and 1470 MeV Au ions to a wide fluence range up to 3 × 10 13 ions/cm 2 . To characterize beam-induced structural modifications across all spatial length scales, complementary experimental techniques such as synchrotron X-ray diffraction and spallation neutron total scattering were used. The phase evolution of the tetragonal polymorph with increasing ion fluence is accurately described by a heterogeneous track-overlap model that incorporates both direct- and double-impact processes. These damage accumulation processes are an expression of a core-shell ion track morphology that depends on irradiation conditions and target material. Neutron pair distribution function analysis revealed that ion-beam-induced tetragonal ZrO 2 is merely a configurational average of short-range orthorhombic (Pbcn) domains stabilized by a dense network of domain walls. Furthermore, this knowledge is critical for a better understanding of how crystalline-to-crystalline phase transformations proceed at the atomic scale under extreme conditions.

Monoclinic-to-tetragonal phase transformation↗
DOE OSTI2025

Unraveling the transformation pathway of the 𝛽 to 𝛾 phase transition in Ga 2 ⁢O 3 from atomistic simulations

Defect spinel 𝛾−Ga 2 ⁢O 3 is the least stable polymorph of Ga 2 ⁢O 3 , so its frequent appearance as a structural defect within or on the surface of monoclinic 𝛽−Ga 2 ⁢O 3 remains a mystery. Through first-principles calculations, we explore potential pathways for the phase transition from 𝛽−Ga 2⁢ O 3 to 𝛾−Ga 2 ⁢O 3 , and examine two key driving forces: tensile strain and Ga deficiency. When configurational entropy contributions to phase energies are included, the 𝛾 phase becomes energetically competitive with the 𝛽 phase, with the free energy difference between these phases diminishing even further under Ga-deficient conditions. Notably, a stability crossover occurs at room temperature at high vacancy concentrations ([V$^{3−}_{Ga}$]>3%) . A simple model 𝛽 → 𝛾 transformation pathway is identified, comprising two primary reactions, that enables the formation of the 𝛾 phase via simultaneous migration of Ga atoms from tetrahedral lattice sites to octahedral interstitial positions. The transformation barriers are prohibitively large in pristine Ga 2 ⁢O 3 , but can be substantially reduced by: (1) the presence of Ga vacancies, (2) elongational strains along the crystallographic 𝑎-axis, and (3) when volumetric relaxations are possible during transformation. These results elucidate prior experimental observations, where 𝛾−Ga 2⁢ O 3 is seen on damaged surfaces or in highly 𝑛-type 𝛽−Ga 2⁢ O 3 environments, which support Ga deficiency and mechanical strain. The insights into the driving forces and mechanisms of 𝛾−Ga 2⁢ O 3 formation enhance understanding of how localized strain and nonequilibrium defect concentrations may facilitate its formation from the 𝛽 phase.

Defects↗
DOE OSTI2025

Atomic-scale structure of ZrO 2 : Formation of metastable polymorphs

Metastable phases can exist within local minima in the potential energy landscape when they are kinetically “trapped” by various processing routes, such as thermal treatment, grain size reduction, chemical doping, interfacial stress, or irradiation. Despite the importance of metastable materials for many technological applications, little is known about the underlying structural mechanisms of the stabilization process and atomic-scale nature of the resulting defective metastable phase. Investigating ion-irradiated and nanocrystalline zirconia with neutron total scattering experiments, we show that metastable tetragonal ZrO 2 consists of an underlying structure of ferroelastic, orthorhombic nanoscale domains stabilized by a network of domain walls. The apparent long-range tetragonal structure that can be recovered to ambient conditions is only the configurational ensemble average of the underlying orthorhombic domains. This structural heterogeneity with a distinct short-range order is more broadly applicable to other nonequilibrium materials and provides insight into the synthesis and recovery of functional metastable phases with unique physical and chemical properties.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2025

Proximity ferroelectricity in wurtzite heterostructures

Proximity ferroelectricity is an interface-associated phenomenon in electric-field-driven polarization reversal in a non-ferroelectric polar material induced by one or more adjacent ferroelectric materials. Here, in this paper, we report proximity ferroelectricity in wurtzite ferroelectric heterostructures. In the present case, the non-ferroelectric layers are AlN and ZnO, whereas the ferroelectric layers are Al 1−x B x N, Al 1−x Sc x N and Zn 1−x Mg x O. The layered structures include nitride–nitride, oxide–oxide and nitride–oxide stacks that feature two-layer (asymmetric) and three-layer (symmetric) configurations. Ferroelectric switching in both layers is validated by multimodal characterization methods, including polarization hysteresis, anisotropic chemical etching, second harmonic generation, piezo response force microscopy, electromechanical testing and atomic resolution polarization orientation imaging in real space by scanning transmission electron microscopy. We present a physical switching model in which antipolar nuclei originate in the ferroelectric layer and propagate towards the internal non-ferroelectric interface. The domain wall leading edge produces elastic and electric fields that extend beyond the interface at close proximity, reducing the switching barrier in the non-ferroelectric layer, and allowing complete domain propagation without breakdown. Density functional theory calculations of polymorph energies, reversal barriers and domain wall energies support this model. Proximity ferroelectricity enables polarization reversal in wurtzites without the chemical or structural disorder that accompanies elemental substitution, opening new questions and opportunities regarding interface-based ferroelectricity.

36 MATERIALS SCIENCE↗
DOE OSTI2026

The Molecular Structures of Liquid and Glassy Nifedipine and Felodipine and Their Incorporation into PVP

Background: Amorphous drug formulations are commonly used to improve the solubility and bioavailability of poorly soluble molecular pharmaceuticals, yet less is known about their molecular conformations and local bonding interactions than their crystalline phases. Methods: High-energy X-ray diffraction structure factor measurements have been made on liquid and glassy nifedipine (NIF), felodipine (FEL), NIF 1:3 polyvinylpyrrolidone (PVP), and FEL 1:3 PVP wt.% mixtures. The corresponding X-ray pair distribution functions have been interpreted using empirical potential structure refinement using different models and density functional theory conformer calculations. Results: In both NIF and FEL, the NH···O inter-molecular hydrogen bonds between the pyridyl nitrogen and ester carbonyls are found to be considerably weaker than those observed in the crystalline polymorphs. For nifedipine, it is proposed that either inter-molecular NH…ON nitro bonds are present and/or a fraction (<20%) of conformational changes, with the aryl ring flipped, occur in the liquid state. For felodipine, the models indicate significant disorder associated with the methyl and ethyl side chains in the liquid state, with the main peak intensity at 3.0 Å arising from intra-molecular Cl-Cl atom pairs. When nifedipine molecules are incorporated into PVP, our models show they possess stronger NH···O bonds to the PVP polymer than felodipine molecules, which have stronger affinity for bonding to the polymer than to other felodipine molecules. Conclusions: The amorphous forms of both NIF and FEL show much weaker hydrogen bonding than found in their crystalline phases. Liquid NIF also exhibits configurations which are not observed in the crystal phases.

Amorphous Solid Dispersion↗
DOE OSTI2025

Unraveling the Atomic Mechanism of the Crystalline Phase‐Dependent Structural Features and Special Spectral Design of α‐, β‐, and Ɛ‐Ga₂O₃

Atomic‐scale phase transformations profoundly influence the functional properties of Ga₂O₃ polymorphs. By combining irradiation experiments with microstructure characterization and theoretical approaches, phase‐specific energy‐dissipation pathways in α‐, β‐, and ε‐Ga₂O₃ are uncovered and strategies for targeted property design are outlined. Competing antiphase boundaries (APBs) and twin domain boundaries (TDBs) promote irreversible α→ε interconversion through domain fragmentation. In β‐Ga₂O₃, defect‐induced stress gradients drive two distinct local transformations: surface Ga‐aggregated β→δ that stabilizes transient states, and latent‐track‐confined β→κ phase transition with recoverable distortions via cation reordering. Under electronic excitation, β‐Ga₂O₃ forms nanohillocks via robust GaO₆ octahedra (high density/strong Ga─O bonds), while α/ε‐Ga₂O₃ generates nanopores from tetrahedral Ga looseness (low bonding energy), highlighting phase‐dependent surface dynamics shaped by atomic packing and bonding anisotropy. Defect‐regulated recombination suppresses visible photoluminescence in α/β‐Ga₂O₃, whereas in ε‐Ga₂O₃ bandgap narrowing of ΔE: 0.30 eV is observed, enhancing emission. Linking phase‐dependent defect‐carrier interactions and metastable‐phase engineering in Ga₂O₃ enables property optimization for power‐electronics and optoelectronics devices.

electronic state configuration↗