Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Structural variation”

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

Modeling Short-Range Order in Disordered Rocksalt Cathodes by Pair Distribution Function Analysis

Pair distribution function (PDF) analysis is a powerful technique for the characterization of short-range order (SRO) in disordered materials. Accurate interpretation of experimental PDF data is critically reliant on the development of structural models that can account for local variations in site occupancies and bond lengths. To this end, we outline an approach to model SRO using first-principles calculations based on the cluster-expansion formalism. These methods are validated on neutron scattering data from two disordered rocksalt oxyfluorides, Li 1.3 Mn 0.4 Ti 0.3 O 1.7 F 0.3 and Li 1.3 Mn 0.4 Nb 0.2 Ti 0.1 O 1.7 F 0.3 . For each composition, we demonstrate that an average structure without any SRO fails to reproduce several key features in the experimental PDF. To pinpoint the origin of the suspected SRO in these materials, configurational and displacive effects were separately investigated using two disparate models. Special quasi-random structures were relaxed using density functional theory to account for local changes in bond lengths while maintaining a near-random ionic configuration. This leads to slightly improved accuracy but still misrepresents asymmetry in the first few peaks of the PDF. Monte Carlo simulations were performed to model configurational SRO on a fixed lattice, which by itself is shown to have a minimal influence on the PDF. Instead, we find that it is the bond length relaxations within environments created by SRO which controls the details of the PDF, thereby highlighting the subtle but important coupling between configurational and displacive SRO in disordered materials.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Ag-doped As–S–Se chalcogenide glasses: a correlative study of structural and dielectrical properties

This work is focused on investigating the effect of Ag doping on the dielectric and structural properties of the As–S–Se glass system and their suitability for potential applications in optoelectronics, such as electronic and switching components. A series of glasses from the system Agx(As40S30Se30)100-x (x = 5 at.% Ag) were prepared with melt-quenching technique. Morphological and compositional analysis was performed by scanning electron microscopy and energy-dispersive X-ray spectroscopy. Raman spectroscopy measurements have shown structural modifications of the glass network upon Ag doping, with creation of new Ag–(S,Se)–As structures, resulting in increased number of defect states. Variations in dielectric permittivity (e') and dielectric loss (e") were investigated in a wide frequency and temperature range, and explained by dominant dipolar and electronic polarization, as well as correlated with the observed structural changes. Improvement in dielectric properties is observed upon doping with Ag. Furthermore, the low dielectric tangent loss at high frequencies for all glassy samples points to a beneficial ability of the material to absorb energy of the external electric field. Good thermal stability of these glasses is confirmed by the calculated temperature coefficient of the dielectric permittivity. Density of localized states in all glasses was determined from AC conductivity and the dielectric tangent loss using a correlated barrier hopping model, and found to increase with Ag concentration. The overall results suggest that Ag doping enhances the electrical and optical quality of this type of glasses, making them suitable as nonlinear optical materials and their applications in devices.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Application of electron backscatter diffraction techniques to quantify effects of aging on sub-grain and spatial heterogeneity in NMC cathodes

Identification and evaluation of structural heterogeneity in spent cathode materials is crucial to developing appropriate remediation strategies for novel recycling processes. Native heterogeneities may be exacerbated during the cell's operational lifetime, as sub-particle-scale variations induce anisotropic expansion and contraction upon cycling. Structural transformations resulting from repeated cycling and calendar aging predominantly occur at the secondary particle surface and at the grain boundaries (GBs) between primary particles. However, the diffusion and stress build up around and across GBs are poorly understood. In this study, electron backscatter diffraction (EBSD) is employed to track sub-grain lattice structure across a statistically relevant number of Li(Ni0.33Mn0.33Co0.33)O2 (NMC-111) particles. Specifically, differences in lattice misorientation – measured as the deviation from the grain's average orientation – are tracked as a function of position within the electrode (near current collector, middle, and near separator) and as a function of electrochemical cycling. Further, a novel method of structural analysis is developed, offering insight into sub-grain diffusion behavior by comparing lattice misorientation near the grain boundary versus in the grain bulk. The present results suggest that electrode-scale spatial heterogeneity in lattice structure is induced by initial manufacturing conditions, and that radial gradients in lattice misorientation evolve at the primary particle scale with repeated electrochemical cycling.

25 ENERGY STORAGE↗
DOE OSTI2021

Correlative analysis of structure and chemistry of Li x FePO 4 platelets using 4D-STEM and X-ray ptychography

Lithium iron phosphate (Li x FePO 4 ), a cathode material used in rechargeable Li-ion batteries, phase separates upon de/lithiation under equilibrium. The interfacial structure and chemistry within these cathode materials affects Li-ion transport, and therefore battery performance. Here, correlative imaging of Li x FePO 4 was performed using four-dimensional scanning transmission electron microscopy (4D-STEM), scanning transmission X-ray microscopy (STXM), and X-ray ptychography in order to analyze the local structure and chemistry of the same particle set. Over 50,000 diffraction patterns from 10 particles provided measurements of both structure and chemistry at a nanoscale spatial resolution (16.6–49.5 nm) over wide (several micron) fields-of-view with statistical robustness. Li x FePO 4 particles at varying stages of delithiation were measured to examine the evolution of structure and chemistry as a function of delithiation. In lithiated and delithiated particles, local variations were observed in the degree of lithiation even while local lattice structures remained comparatively constant, and calculation of linear coefficients of chemical expansion suggest pinning of the lattice structures in these populations. Partially delithiated particles displayed broadly core–shell-like structures, however, with highly variable behavior both locally and per individual particle that exhibited distinctive intermediate regions at the interface between phases, and pockets within the lithiated core that correspond to FePO 4 in structure and chemistry. The results provide insight into the Li x FePO 4 system, subtleties in the scope and applicability of Vegard’s law (linear lattice parameter-composition behavior) under local versus global measurements, and demonstrate a powerful new combination of experimental and analytical modalities for bridging the crucial gap between local and statistical characterization.

25 ENERGY STORAGE↗
DOE OSTI2022

Fire, insect and disease‐caused tree mortalities increased in forests of greater structural diversity during drought

Abstract Structural diversity is an emerging dimension of biodiversity that accounts for size variations in organs among individuals in a community. Previous studies show significant effects of structural diversity on forest growth, but its effects on forest mortality are not known, particularly at a large scale. To address this knowledge gap, we quantified structural diversity using stem structural diversity (SSD) based on both tree diameter and height. We obtained U.S. Forest Service Forest Inventory and Analysis (FIA) data from over 2400 plots across southcentral U.S. forests that have suffered a recent drought. Using data from multiple sampling times, we calculated SSD and compared the relative importance of SSD, species diversity, functional diversity and other stand attributes in determining tree mortalities caused by fire, insects and diseases. We also used FIRETEC, a physics‐based fire model, to test the effect of SSD on canopy consumption by fire. Our results showed that (1) SSD was positively associated with tree mortalities caused by all three disturbances; (2) species richness was negatively associated with insect‐ and disease‐caused mortalities; (3) functional diversity was negatively associated with fire‐ and disease‐caused mortalities and (4) more phylogenetically related species had more similar mortality rates by insect and disease but not fire. Moreover, the FIRETEC model showed increasing canopy consumption by fire in stands with greater SSD. Together, the different tree mortalities during drought associated with SSD more consistently than the other biodiversity metrics were evaluated. Synthesis . Our results suggest that SSD could be considered in modelling forest dynamics and planning management to sustain forest health under disturbances.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

How Do Ionic Liquids “Fold” Ionenes? Computational and Experimental Analysis of Imidazolium Polymers Based on Ether and Alkyl Chain Variations Dissolved in an Ionic Liquid

Molecular dynamics (MD) simulations and complementary experiments are used here to understand the inter- and intramolecular structure and conformational properties of ionenes (i.e., cationic polymers formed from condensation reactions) dissolved in an ionic liquid (IL), 1-ethyl-3-methylimidazolium bistriflimide [C 2 mim + ][Tf 2 N - ]. The simulated structural properties are benchmarked against experimental analyses of these same polymers, mainly using dynamic light scattering (DLS) experiments. Four different imidazolium ionenes are considered, corresponding to variations in the chemistry and structure of the repeat units including: (1) poly (decylimidazolium) (PD 10 ); (2) poly (tetraethyleneglycolimidazolium) (PE 10 ); (3) alternating copolymer (P(ED) 5 ); and (4) a block copolymer (PE 5 D 5 ). Detailed computational analyses of the polymer structure and conformational properties were performed, including the radius of gyration, end-to-end distance, torsional distributions, and site-site and spatial distribution functions. Overall, there is a competition between intramolecular associations between the imidazolium groups and the ether sites in the polymers versus the intermolecular interactions of the polymer imidazolium groups with the surrounding anion molecules. The polymer with only ether linkages (PE 10 ) results in the most intramolecular interactions, leading to significant coiling behavior and chain contraction. These strong interactions reduce chain flexibility, but they also result in much more linear chain configurations and alignment of the imidazolium groups. The polymer structure is strongly affected as the concentration of alkyl groups is increased, and we find that the polymer architecture (alternating versus block copolymer) also has an important influence. Furthermore, these detailed observations are important for understanding the connection between the molecular design of imidazolium polymers and their emergent structural properties, which could lead to unique opportunities for creating polymer composites.

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

Testing unified theories for ozone response in C 4 species

Abstract There is tremendous interspecific variability in O 3 sensitivity among C 3 species, but variation among C 4 species has been less clearly documented. It is also unclear whether stomatal conductance and leaf structure such as leaf mass per area (LMA) determine the variation in sensitivity to O 3 across species. In this study, we investigated leaf morphological, chemical, and photosynthetic responses of 22 genotypes of four C 4 bioenergy species (switchgrass, sorghum, maize, and miscanthus) to elevated O 3 in side‐by‐side field experiments using free‐air O 3 concentration enrichment (FACE). The C 4 species varied largely in leaf morphology, physiology, and nutrient composition. Elevated O 3 did not alter leaf morphology, nutrient content, stomatal conductance, chlorophyll fluorescence, and respiration in most genotypes but reduced net CO 2 assimilation in maize and photosynthetic capacity in sorghum and maize. Species with lower LMA and higher stomatal conductance tended to show greater losses in photosynthetic rate and capacity in elevated O 3 compared with species with higher LMA and lower stomatal conductance. Stomatal conductance was the strongest determinant of leaf photosynthetic rate and capacity. The response of both area‐ and mass‐based leaf photosynthetic rate and capacity to elevated O 3 were not affected by LMA directly but negatively influenced by LMA indirectly through stomatal conductance. These results demonstrate that there is significant variation in O 3 sensitivity among C 4 species with maize and sorghum showing greater sensitivity of photosynthesis to O 3 than switchgrass and miscanthus. Interspecific variation in O 3 sensitivity was determined by direct effects of stomatal conductance and indirect effects of LMA. This is the first study to provide a test of unifying theories explaining variation in O 3 sensitivity in C 4 bioenergy grasses. These findings advance understanding of O 3 tolerance in C 4 grasses and could aid in optimal placement of diverse C 4 bioenergy feedstock across a polluted landscape.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Non-Solvent Induced Phase Separation Enables Designer Redox Flow Battery Electrodes

Porous carbonaceous electrodes are performance-defining components in redox flow batteries (RFBs), where their properties impact the efficiency, cost, and durability of the system. The overarching challenge is to simultaneously fulfill multiple seemingly contradictory requirements—i.e., high surface area, low pressure drop, and facile mass transport—without sacrificing scalability or manufacturability. Here, non-solvent induced phase separation (NIPS) is proposed as a versatile method to synthesize tunable porous structures suitable for use as RFB electrodes. The variation of the relative concentration of scaffold-forming polyacrylonitrile to pore-forming poly(vinylpyrrolidone) is demonstrated to result in electrodes with distinct microstructure and porosity. Tomographic microscopy, porosimetry, and spectroscopy are used to characterize the 3D structure and surface chemistry. Flow cell studies with two common redox species (i.e., all-vanadium and Fe2+/3+) reveal that the novel electrodes can outperform traditional carbon fiber electrodes. It is posited that the bimodal porous structure, with interconnected large (>50 µm) macrovoids in the through-plane direction and smaller (<5 µm) pores throughout, provides a favorable balance between offsetting traits. Although nascent, the NIPS synthesis approach has the potential to serve as a technology platform for the development of porous electrodes specifically designed to enable electrochemical flow technologies.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Tuning the structure and superconductivity of SrNi 2⁢ P 2 by Rh substitution

The compound SrNi 2⁢ P 2 is unique among the ThCr 2 ⁢Si 2 class since it exhibits a temperature-induced transition upon cooling from an uncollapsed tetragonal (ucT) state to a one-third-collapsed orthorhombic (tcO) state where one out of every three P-rows bond across the Sr layers. This compound is also known for exhibiting bulk superconductivity below 1.4 K at ambient pressure. Here, in this paper, we report on the effects of Rh substitution in Sr⁢(Ni 1-x ⁢Rh x ) 2 ⁢P 2 on the structural and superconducting properties. We studied the variation of the nearest P-P distances as a function of the Rh fraction at room temperature, as well as its temperature dependence for selected compositions. We find that increasing the Rh fraction leads to a decrease in the transition temperature between the ucT and tcO states, until a full suppression of the tcO state for x ≥ 0.166. The superconducting transition first remains nearly insensitive to the Rh fraction, and then it increases to 2.3 K after the tcO state is fully suppressed. These results are summarized in a phase diagram, built upon the characterization by energy dispersive x-ray spectroscopy, x-ray diffraction, resistance, magnetization, and specific heat measurements done on crystalline samples with varying Rh content. The relationship between band structure, crystal structure, and superconductivity is discussed based on previously reported band structure calculations on SrRh 2 ⁢P 2 . Moreover, the effect of Rh fraction on the stress-induced structural transitions is also addressed by means of strain-stress studies done by uniaxial compression of single-crystalline micropillars of Sr⁢(Ni 1-x ⁢Rh x ) 2 ⁢P 2 .

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Effects of thermophoresis on high-pressure binary-species boundary layers with uniform and non-uniform compositions

Direct numerical simulations of high-pressure binary-species temporal boundary layers are performed to investigate the flow physics for three situations: (1) uniform and equal composition, (2) uniform but unequal compositions and (3) non-uniform composition. Both colder- and hotter-wall situations compared with the free stream are simulated. The working fluid is a nitrogen/methane mixture. The analysis is performed at a case-specific self-similar state. Even when the initial composition is uniform, the methane mean mass fraction decreases near the colder wall, whereas it increases near the hotter wall and the mass fraction fluctuates in the entire boundary layer. Analysis of the species-mass diffusion balance and flow structures reveal that both mass-fraction variation and fluctuations are induced by the Soret effect. When the initial composition is non-uniform and the wall is colder, the methane mean mass fraction monotonically increases from the wall akin to its initial profile. However, when the wall is hotter the mean mass fraction decreases near the wall in contrast to its initial profile, a fact traced through the species-mass diffusion balance to the Soret effect being large and enriching methane near the wall. In contrast, the direction of the Soret flux is opposite for the colder wall, thus keeping the methane concentration small. Although the initial magnitude of the difference between the wall and free-stream temperature is the same in all cases, the situation is not symmetric between colder- and hotter-wall cases; the flow structure exhibits much smaller scales when the wall is hotter than when the wall is colder.

Mechanics↗
DOE OSTI2026

Lattice Relaxations in Bixbyite‐Structured Rare‐Earth Sesquioxides

This report examines atomic relaxations in rare‐earth sesquioxides with the cubic bixbyite crystal structure. In particular, the experimentally determined unit cell atom positions in bixbyite compounds are compared to “ideal” bixbyite atom positions, where the concept of ideality is based on similarities between the bixbyite and fluorite crystal structures. Analyses indicate that all rare‐earth sesquioxide compounds deviate from ideality in similar ways and that there is a trend wherein the magnitudes of these atomic deviations decrease with increasing cationic size. The deviations from an ideal, fluorite‐like structure produce highly distorted coordination polyhedra and variations in metal‐oxygen bond lengths. A modified bond valence sum analysis revealed that the first nearest neighbor cation‐anion bond length variations lead to large (20%) variations in cation‐anion bond strength. This report also provides a qualitative analysis regarding the accuracy of atomic relaxed deviation parameters, based on x‐ray diffraction measurements. The analysis indicates that uncertainties in the accuracy of measured oxygen relaxed deviation parameters in bixbyite sesquioxides are likely much greater than the published uncertainties.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Definitive Assessment of the Accuracy, Variationality, and Convergence of Relativistic Coupled Cluster and Density Matrix Renormalization Group in 100-Orbital Space

Accuracy, variationality, and convergence underpin the reliability of modern electronic structure methods, yet definitive benchmarks in the relativistic regime remain elusive due to the absence of numerically exact full configuration interaction (CI) references. Recent algorithmic advances in the CI framework, enabled by the small-tensor-product (STP) decomposition approach, have dramatically extended the tractable size of the configuration space, making numerically exact CI calculations feasible in large active spaces previously beyond reach. In this paper, we employ the recently developed STP-CI framework to perform large-scale numerically exact CI calculations and directly benchmark relativistic coupled cluster and density matrix renormalization group methods. Definitive benchmarking of approximate relativistic electronic structure methods is ensured through the application of the gap theorem, which provides rigorous error bounds on the CI reference and establishes a controlled standard for assessing accuracy, variationality, and convergence.

Chemical calculations↗
DOE OSTI2023

Molecular-dynamics simulation methods for macromolecular crystallography

It is investigated whether molecular-dynamics (MD) simulations can be used to enhance macromolecular crystallography (MX) studies. Historically, protein crystal structures have been described using a single set of atomic coordinates. Because conformational variation is important for protein function, researchers now often build models that contain multiple structures. Methods for building such models can fail, however, in regions where the crystallographic density is difficult to interpret, for example at the protein–solvent interface. To address this limitation, a set of MD–MX methods that combine MD simulations of protein crystals with conventional modeling and refinement tools have been developed. In an application to a cyclic adenosine monophosphate-dependent protein kinase at room temperature, the procedure improved the interpretation of ambiguous density, yielding an alternative water model and a revised protein model including multiple conformations. The revised model provides mechanistic insights into the catalytic and regulatory interactions of the enzyme. The same methods may be used in other MX studies to seek mechanistic insights.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Constraining level densities through quantitative correlations with cross-section data

The adopted level densities (LD) for the nuclei produced through different reaction mechanisms significantly impact the accurate calculation of cross sections for the different reaction channels. Many common LD models make simplified assumptions regarding the overall behavior of the total LD and the intrinsic spin and parity distributions of the excited states. However, very few experimental constraints are taken into account in these models: LD at neutron separation energy coming from average spacings of s- and p-wave resonances ( D 0 and D 1 , respectively) whenever they have been previously measured, and the sometimes subjective extrapolation of discrete levels. These, however, constrain the LD only in very specific regions of excitation energy, and for specific spins and parities. This work aims to establish additional experimental constraints on LD through quantitative correlations between cross sections and LD. This allows for fitting and the determination of detailed structures in LD. For this we use the microscopic Hartree-Fock-Bogoliubov (HFB) LD as a starting point as the HFB LD provide a more realistic spin and parity distributions than phenomenological models such as Gilbert-Cameron (GC). We then associate variations predicted by the HFB model with the structure observed in double-differential cross sections at low outgoing neutron energy, a region that is dominated by the LD input. We also use (n, p) on 56 Fe , as an example case where angle-integrated cross sections are extremely sensitive to LD. For comparison purposes we also perform calculations with the GC model. With this approach we are able to perform fits of the LD based on actual experimental data, constraining the model and ensuring its consistency. This approach can be particularly useful in extrapolating the LD to nuclei for which high-excited discrete levels and/or values of D 0 or D 1 are unknown. Finally, it also predicts neutron-induced inelastic γ cross sections that in some cases can differ significantly from more standard phenomenological LD models such as GC.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

The influence of cooling rates on strain phase diagrams and domain structures of ferroelectric thin films: A case study of PbTiO 3

Strain engineering has been established as an effective approach to control phase equilibria, domain configurations, and functional properties of ferroelectric thin films. Temperature-strain phase diagrams have been used as powerful tools for providing insights into strain engineering. However, almost all existing phase diagrams established using the phase-field approach assume quenching conditions without considering actual cooling rates during the post-deposition annealing process of ferroelectric thin films. Within this work, we systematically investigate the influence of cooling rates on domain structures and the strain-phase diagram of ferroelectric thin films using phase-field simulations, taking PbTiO 3 thin films as a model system. We found that both the position of phase boundaries in the strain phase diagrams and the domain morphology are significantly influenced by the cooling rates. It is revealed that while the paraelectric-ferroelectric phase boundary remains invariant, the phase boundaries between single-phase and multi-phase regions tend to shift toward the corresponding multi-phase region as the cool rate reduces. Slow cooling generally leads to more ordered domain structures with increased domain size. Using the obtained equilibrium domain structures, we calculated effective thermal conductivities and found significant variations that can be tuned by the cooling rates. In conclusion, this work reveals an underexplored yet critical impact of cooling rates on phase equilibria and domain structures in ferroelectric thin films, which may inspire further fine-tuning of domains and domain walls in low-dimensional ferroelectrics for multifunctional applications.

Domain structures↗
DOE OSTI2022

Simulation of a PCM Integrated Heat Pump Using Time-of-Use Utility Structure based Control Strategy for Demand Response

To mitigate the variation in demand on the electric grid, thermal energy storage (TES) is an alternative to electric batteries and constructing new peaking power plants. The aggressive carbon reduction goals of the United States are accelerating the electrification of building equipment. If space heating and cooling in the US were electrified using conventional heat pump technologies, the winter electric grid peak would approximately double. Integrating TES with heat pump (HP) enables electrification of space cooling and heating devices without overtaxing the grid.This paper evaluates the benefits of a novel phase change material (PCM) integrated heat pump configuration via simulation. A one-dimensional PCM heat exchanger model which discretizes the PCM tank and refrigerant tubes into small control volumes is developed. Each control volume can have different PCM temperatures, PCM properties, and heat transfer coefficients. The PCM tank is charged by a wrapped tank condenser and discharged by an internal refrigerant coil. The PCM heat exchanger model is integrated into DOE/ORNL Heat Pump Design Model for heat pump system simulation.To demonstrate the performance of the PCM integrated heat pump, a case study in Chicago was performed. A Time-of-Use utility structure-based control strategy is implemented to schedule the PCM tank charging and discharging mode switching. Compared with a conventional electric heat pump, the PCM integrated heat pump shows superior performance on load shifting and utility cost reduction. As a result, the proposed system demonstrates 24.6% utility saving for cooling application and 25.8% utility saving for heating application.

Li, Zhenning↗
DOE OSTI2021

Determining the Three-Dimensional Structures of Silica-Supported Metal Complexes from the Ground Up

The immobilization of molecularly precise metal complexes to substrates, such as silica, provides an attractive platform for the design of active sites in heterogeneous catalysts. Specific steric and electronic variations of the ligand environment enable the development of structure–activity relationships and the knowledge-driven design of catalysts. At present, however, the three-dimensional environment of the precatalyst, much less the active site, is generally not known for heterogeneous single-site catalysts. We explored the degree to which NMR-based surface-to-complex interatomic distances could be used to solve the three-dimensional structures of three silica-supported metal complexes. Furthermore, the structure solution revealed unexpected features related to the environment around the metal that would be difficult to discern otherwise. This approach appears to be highly robust and, due to its simplicity, is readily applied to most single-site catalysts with little extra effort.

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

Predicting interface structure using the minima hopping method

Here, we adapt the minima hopping method (MHM) to the problem of interfacial structure prediction and apply it to study a canonical problem, the tilt grain boundaries in SrTiO 3 . Our method employs a hybrid approach by first exploring the potential energy surface (PES) of different grain boundary samplings with an empirical force field, among which the fifteen candidates with lower energies are then refined using ab initio density functional theory (DFT) calculations. During the exploratory stage, we bias the search using a local order parameter to primarily sample various reconstructions in the vicinity of the interface, while preserving the crystallinity of the bulk regions. We further enhance the search by incorporating initial structures with rigid body displacements to account for translational variations between bulk phases, enabling the MHM to effectively generate both stoichiometric and nonstoichiometric SrTiO 3 Σ⁢3(111)[110] and Σ⁢3(112)[110] grain boundaries. From an algorithmic standpoint, MHM outperforms earlier studies based on genetic algorithms (GA) by identifying more stable interfacial structures of several SrTiO 3 grain boundaries. The performance of the present implementation of the MHM approach is primarily limited by exploring an approximate description of the PES with a rather simple Buckingham potential. This limitation leads to variations in performance when compared to approaches utilizing more advanced surrogate PES models, such as direct DFT-PES sampling or GA with the embedded atom method (EAM). Despite the present limitations, the MHM approach is able to yield interfacial structures with comparable or lower interfacial energies in specific cases, such as Σ⁢3(111)[110] Γ=1, ±0.5 and Σ⁢3(112)[110] Γ= ±1, −2, underscoring the robustness of the MHM approach even with a simple approximation of the DFT PES. The MHM interfacial structure prediction method thus offers an efficient approach to understanding the grain boundaries and heterointerfaces at the atomic scale, providing an important prerequisite for effective materials design.

density functional theory↗