Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Single crystal model”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

A phase-field model for non-isothermal phase transformation and plasticity in polycrystalline yttria-stabilized tetragonal zirconia

Here we propose an elastoplastic phase-field (PF) model to investigate the mechanics of tetragonal-to-monoclinic phase transformation (TMPT) and elastoplastic deformation of polycrystalline yttria-stabilized tetragonal zirconia (YSTZ). A Landau polynomial with non-vanishing chemical energy at the equilibrium temperature is introduced to account for the actual formation energies of the phases. The effects of different grain orientations, latent heat, and temperature on TMPT and deformation mechanisms are considered. The suppressive transformation effects of the grain boundaries (GBs) is modeled using an inhomogeneous kinetic coefficient in the bulk and GBs. The simulation results for single crystals demonstrate the capability of the model to reproduce the orientation-dependent compressive deformation of YSTZ similar to atomistic simulations and micropillar experiments. The single crystal with [100] crystallographic orientation along the loading direction (SC[100]) displays both TMPT and plasticity, SC[101] experiences only phase transformation, while SC[001] undergoes only plastic yielding. The TMPT induced by compressive loading exhibits shape memory effect (SME) below the equilibrium transformation temperature and pseudoelasticity above it, while the critical transformation stress increases with increasing loading temperature. The irrecoverable plastic strain is found to trap a part of the monoclinic phase, which prevents a complete reverse transformation. The polycrystalline cases also display SME and PE at low and high temperatures, respectively. Due to the orientation differences between grains and the stress concentrations at geometric nonlinearities, plastic deformation occurs in polycrystalline YSTZ for an applied load less than the yield stress. The results suggest a possible limitation of plasticity and an improvement of the shape recovery of YSTZ if one can control the orientation of the grains and/or increase the density of stacking faults at the GBs during material processing.

36 MATERIALS SCIENCE↗
DOE OSTI2020

One Structure, Two Elements—LuGe 2 Superconductor vs Ordinary Metallic Conductor LuSn 2 . A Case Study on How Site-Selective Germanium for Tin Atom Substitution Leads to Modulating of the Charge Distribution

The substitution of chemically similar elements in a given crystal structure is an effective way to enhance physical properties, but the understanding on such improvements is usually impeded because the substitutions are random, and the roles of the different atoms cannot be distinguished by crystallographic symmetry. Herein, we provide a detailed crystallographic analysis and property measurements for the continuous solid solutions LuGe x Sn 2–x (0 < x < 2). The results show that there is no apparent change of the global symmetry, with the end-members LuGe 2 and LuSn 2 , as well as the intermediate LuGe x Sn 2–x compositions adopting the ZrSi 2 type structure (space group Cmcm, Pearson index oC12). Yet, the refinements of the crystal structures from single-crystal X-ray diffraction data show that Ge–Sn atom substitutions are not random, but occur preferentially at the zig-zag chain. The patterned distribution of two group 14 elements leads to a significant variation in chemical bonding and charge ordering within the other structural fragment, the 2D square nets, thereby resulting in tuned electron transport. The enhancement is greater than the typical Bloch-Gruneisen model, and more akin to the parallel-resistor model. Magnetization measurements on single crystals show bulk superconductivity in all LuGe x Sn 2–x samples with shielding fractions as high as 90%. Specific heat data confirm the effect to originate from residual metallic tin in the material, indicating that Sn atom substitutions in the 2D square nets cause disruptions of the hypervalent bonding and local anisotropy, which ultimately leads to vanishing of the superconducting state in the end-member LuGe 2 . This work sheds light on how the complexity in chemical interactions by two different carbon congeners leads to changes in the physical properties and how they can be correlated with the induced charge distribution. Furthermore, these studies also provide a general approach to modulation of charge density, and thus, of emerging physical properties in other classes of intermetallic systems based on the main-group elements of groups 13 to 15.

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

Machine learning inversion of interatomic force constants from single-crystal inelastic neutron scattering

Atomic vibrations govern many macroscopic properties of materials, but experiments to comprehensively probe them remain challenging. Inelastic neutron scattering (INS) is a powerful technique to map phonon dispersions in crystals, especially when leveraging modern time-of-flight (ToF) spectrometers with large detectors. However, efficiently and robustly extracting interatomic force constants (FCs) parameterizing phonon dynamics from experimental spectra remains a bottleneck due to the complexity and high dimensionality of ToF INS datasets. Here, we present a machine learning approach for the direct inversion of FCs from single-crystal INS measurements. The framework leverages synthetic training data generated using universal machine-learned force fields and an efficient physics-based forward model. We benchmark two neural architectures–one emphasizing structured latent representation learning and the other direct, supervised spectral regression–across simulated datasets for two materials under idealized and noisy conditions. The latent-representation model is subsequently applied to experimental single-crystal INS data on germanium. The model is shown to reproduce FCs derived from both first-principles simulations and from iterative optimization, and furthermore achieves reliable inference even from sparse, single-orientation measurements representing short data acquisitions. Analysis of the learned latent space reveals semantically continuous and physically interpretable encodings that support strong cross-domain generalization. By bridging theoretical and experimental domains, we establish a path toward rapid inversion of experimental spectra and data-driven interpretation of temperature-dependent lattice dynamics.

42 ENGINEERING↗
DOE OSTI2022

Accurate quantification of lattice temperature dynamics from ultrafast electron diffraction of single-crystal films using dynamical scattering simulations

In ultrafast electron diffraction (UED) experiments, accurate retrieval of time-resolved structural parameters, such as atomic coordinates and thermal displacement parameters, requires an accurate scattering model. Unfortunately, kinematical models are often inaccurate even for relativistic electron probes, especially for dense, oriented single crystals where strong channeling and multiple scattering effects are present. This article introduces and demonstrates dynamical scattering models tailored for quantitative analysis of UED experiments performed on single-crystal films. As a case study, we examine ultrafast laser heating of single-crystal gold films. Comparison of kinematical and dynamical models reveals the strong effects of dynamical scattering within nm-scale films and their dependence on sample topography and probe kinetic energy. Applying to UED experiments on an 11 nm thick film using 750 keV electron probe pulses, the dynamical models provide a tenfold improvement over a comparable kinematical model in matching the measured UED patterns. Also, the retrieved lattice temperature rise is in very good agreement with predictions based on previously measured optical constants of gold, whereas fitting the Debye–Waller factor retrieves values that are more than three times lower. Altogether, these results show the importance of a dynamical scattering theory for quantitative analysis of UED and demonstrate models that can be practically applied to single-crystal materials and heterostructures.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Asymmetry of the Ferroelectric Phase Transition in BaTiO 3

Symmetry changes during phase transformations fundamentally determine the behavior of thermodynamic systems, governing phenomena as diverse as water evaporation, fermion condensation, and epidemic spreading. Phase transitions are conventionally divided into two classes, and this classification is typically assumed to remain invariant upon reversing the transition. Here, an asymmetric phase transformation is uncovered in the ferroelectric–paraelectric transition of single-crystal BaTiO 3 , a model system for first-order transitions. Under slow temperature variation (≤0.1 °C min −1 ), thermodynamic, dielectric, and domain-structure measurements reveal that the ferroelectric-to-paraelectric transition exhibits latent heat, phase coexistence, and a discontinuous order parameter, while these signatures are absent upon cooling. Complementary phase-field simulations demonstrate similar behavior, attributing it to distinct elastic strain energy accumulation and release during heating and cooling. These findings reveal a first-order character upon heating but a second-order-like behavior upon cooling, challenging the conventional paradigm of symmetric phase-transition classification and suggesting new possibilities for ferroelectric-based energy storage.

Hershkovitz, Asaf [Technion-Israel Institute of Te↗
DOE OSTI2025

Constructing Water‐Stable Porous Organic Salts via Suppressed Proton Integration Using Fluorinated Tetrazole Tectons

Porous organic salts (POSs) are an emerging class of materials with ordered ionic architectures, offering excellent proton transfer and water uptake properties. However, conventional POS synthesis via strong acid–base neutralization (e.g., ─SO₃H and ─NH₂) leads to extensive hydrogen bonding with water, compromising stability in aqueous and water-lean environments. Here, we address this challenge by designing POSs with hydrophobic porous channels and minimal hydrogen bonding formation. Our key innovation is the use of fluorinated tetrazole as a weak acid tecton and a tetra-substituted imidazole precursor devoid of active protons as the base. Single-crystal analysis and computational modeling reveal that the structural integrity of the synthesized POSs arises primarily from cation–anion interactions, with water confined as clusters in the pores, independent of hydrogen bonding with the scaffold. Robustness of the POS structure under aqueous and water-lean conditions is confirmed by X-ray and neutron scattering, as well as computational modeling, confirming preserved packing and crystal structures. The stability of POS is further demonstrated in aqueous iodine capture, with imidazolium cations and C–F functionalizations serving as strong adsorption sites. As a result, the approach developed herein further pushes the boundary of POS materials to withstand both aqueous and water-lean conditions.

Fluorinated tecton↗
DOE OSTI2021

Role of crystallographic orientation on intragranular void growth in polycrystalline FCC materials

In this study, we study the effect of crystallographic orientation and applied triaxiality on the growth of intragranular voids. Two 3D full-field micromechanics methods are used, the dilatational visco-plastic fast-Fourier transform (DVP-FFT) and the crystal plasticity Finite Elements (CP-FE), both of which incorporate a combination of crystalline plasticity and dilatational plasticity. We demonstrate with several select cases that predictions of void growth from both formulations agree qualitatively. With the more computationally efficient DVP-FFT, additional effects of polycrystalline microstructure and the influence of nearest neighborhood are investigated. Crystals bearing a single intracrystalline void are studied in three types of 3D microstructural environments: isolated single crystals, individual equal-sized grains within a regular polycrystal, and individual variable sized grains within a polycrystal with grains and voids randomly located. We show that loading type plays a significant role. In strain-rate controlled conditions, voids in the hardest [111]-crystals grow the fastest in time, whereas in stress-controlled conditions, voids in the softest [100]-crystal grow the fastest in time. The analysis reveals that on average void growth is slower for the same starting orientation in the polycrystal than in the single crystal. We find that at the highest triaxiality tested that the correlation between crystal orientation and void growth rate in the polycrystal strengthens, drawing closer to that seen in the isolated single crystals. These results and model can help guide the microstructural design of polycrystalline materials with high strength and damage-tolerance in high-rate deformation.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Modifying Interface Solvation and Oxygen Reduction Electrocatalysis with Hydrophobic Species

Ionic liquid interlayers improve the oxygen reduction reaction (ORR) kinetics on bulk and nanostructured catalysts for both Pt and alloyed-Pt materials. Despite the demonstrated performance enhancement at the half-cell and membrane electrode assembly level, the mechanism of the improvement is not fully understood. In this work, we combine single-crystal experiments with microkinetic modeling to uncover the origin of the ORR kinetic improvement on Pt(111) in the presence of ionic liquids and hydrophobic cations. With the incorporation of a modified Frumkin isotherm, our model accurately simulates the disorder–order transition observed in hydroxyl and bisulfate adsorption on Pt(111) under acidic conditions. Voltametric analysis shows that ionic liquids impact solvation to break so-called scaling relations between the adsorption strength of OH ad and O ad , but these effects have little impact on ORR activity. Instead, destabilized OH ad reduces the overall hydroxyl (spectator) coverage, resulting in higher availability of active sites.

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

A Polycrystalline Pd Surface Studied by Two-Dimensional Surface Optical Reflectance during CO Oxidation: Bridging the Materials Gap

Industrial catalysts are complex materials systems operating in harsh environments. The active parts of the catalysts are nanoparticles that expose different facets with different surface orientations at which the catalytic reactions occur. However, these facets are close to impossible to study in detail under industrially relevant operating conditions. Instead, simpler model systems, such as single crystals with a well-defined surface orientation, have been successfully used to study gas–surface interactions such as adsorption and desorption, surface oxidation, and oxidation/reduction reactions. To more closely mimic the many facets exhibited by nanoparticles and thereby close the so-called materials gap, there has also been a recent move toward using polycrystalline surfaces and curved crystals. However, these studies are limited either by the pressure or spatial resolution at realistic pressures or by the number of surfaces studied simultaneously. In this work, we demonstrate the use of reflectance microscopy to study a vast number of catalytically active surfaces simultaneously under realistic and identical reaction conditions. As a proof of concept, we have conducted an operando experiment to study CO oxidation over a Pd polycrystal, where the polycrystalline surface acts as a collection of many single-crystal surfaces. Finally, we visualized the resulting data by plotting the reflectivity as a function of surface orientation. We think the techniques and visualization methods introduced in this work will be key toward bridging the materials gap in catalysis.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Elucidating the active phases of CoOx films on Au(111) in the CO oxidation reaction

Abstract Noble metals supported on reducible oxides, like CoO x and TiO x , exhibit superior activity in many chemical reactions, but the origin of the increased activity is not well understood. To answer this question we studied thin films of CoO x supported on an Au(111) single crystal surface as a model for the CO oxidation reaction. We show that three reaction regimes exist in response to chemical and topographic restructuring of the CoO x catalyst as a function of reactant gas phase CO/O 2 stoichiometry and temperature. Under oxygen-lean conditions and moderate temperatures (≤150 °C), partially oxidized films (CoO x<1 ) containing Co 0 were found to be efficient catalysts. In contrast, stoichiometric CoO films containing only Co 2+ form carbonates in the presence of CO that poison the reaction below 300 °C. Under oxygen-rich conditions a more oxidized catalyst phase (CoO x>1 ) forms containing Co 3+ species that are effective in a wide temperature range. Resonant photoemission spectroscopy (ResPES) revealed the unique role of Co 3+ sites in catalyzing the CO oxidation. Density function theory (DFT) calculations provided deeper insights into the pathway and free energy barriers for the reactions on these oxide phases. These findings in this work highlight the versatility of catalysts and their evolution to form different active phases, both topological and chemically, in response to reaction conditions exposing a new paradigm in the catalyst structure during operation.

Chen, Hao (ORCID:0000000337598352)↗
DOE OSTI2023

Light-driven hydrogen evolution via a novel pincer/no pincer mechanism including a possible concerted proton electron transfer

The road to efficient molecular catalysts for artificial photosynthesis requires a great deal of basic research to address the challenges of anthropogenic global warming. An unexpected result of such research is a novel pincer/no pincer mechanism discovered for the DFT-calculated catalytic cycle of PNP-C2 complexes. Structural parameters from theoretical modeling are supported by single crystal X-ray data for the corresponding complexes. Using the examples of [M(II)Cl(PNP-C2-R)] + (M = Pd 2+ , Ni 2+ ; R = CH 3 , H; PNP-C2 = N,N-bis{(di(2-methoxy-phenyl)phosphanyl)ethyl}-N-alkylamine) this pincer/no pincer mechanism is supported by differential pulse voltammetric, cyclic voltammetric, UV-Vis spectroscopic (combined with voltammetry) and 31 P nuclear magnetic resonance spectroscopic studies. The efficiency of hydrogen evolution during artificial photosynthesis depends significantly on the ratio of photosensitiser and water reduction catalysts. At first sight, the results from irradiation experiments (gas chromatography and mass spectrometry) seem obvious, but show a considerable impact of secondary and tertiary amines within the pincer structure. The best result obtained for the palladium pincer complex was with a secondary amine as a proton relay, delivering a turnover number of 2237 after only 2.5 hours. In contrast the analogous water reduction catalyst with a methylated (tertiary) amine instead of a secondary one gives only 972 turnovers under the same experimental conditions, over the same time. For all photocatalytic measurements the amount of hydrogen is approximately 2.5-times higher when a secondary amine is used as the proton relay. The influences of these different amines on the activity of water reduction catalysts during irradiation experiments in H 2 O or D 2 O lead to evidence of a kinetic isotope effect, supporting a possible concerted proton electron transfer. DFT calculations of reaction enthalpies, activation energies and intrinsic reaction coordinates of the catalytic cycle support the high efficiency of the presented catalysts. In particular, results from secondary amines suggest new candidates as examples of the rarely known concerted proton electron transfer.

14 SOLAR ENERGY↗
DOE OSTI2024

Electron mobility enhancement by electric field engineering of AlN/GaN/AlN quantum-well HEMTs on single-crystal AlN substrates

To enhance the electron mobility in quantum-well high-electron-mobility transistors (QW HEMTs), we investigate the transport properties in AlN/GaN/AlN heterostructures on Al-polar single-crystal AlN substrates. Theoretical modeling combined with experiment shows that interface roughness scattering due to high electric field in the quantum well limits mobility. Increasing the width of the quantum well to its relaxed form reduces the internal electric field and scattering, resulting in a binary QW HEMT with a high two-dimensional electron gas (2DEG) density of 3.68×1013 cm–2, a mobility of 823 cm2/Vs, and a record-low room temperature (RT) sheet resistance of 206 Ω/□. Further reduction of the quantum well electric field yields a 2DEG density of 2.53×1013 cm–2 and RT mobility > 1000 cm2/V s. These findings will enable future developments in high-voltage and high-power microwave applications on the ultrawide bandgap AlN substrate platform.

Physics↗
DOE OSTI2025

Boundary scattering in topological Kondo insulator SmB 6

We have studied the effects of phonon-boundary scattering on the thermal transport of topological Kondo insulator SmB6. The studies have been performed using the 3ω method across a temperature range 3–300 K. Our results indicate that the thermal conductivity of micro-sized SmB 6 is of an order of magnitude smaller than that of a bulk single crystal. Using the Callaway model, we analyzed the low-temperature lattice thermal conductivity of the microcrystal and demonstrated that phonon scattering at the sample boundaries is a major contributor to the thermal resistance in this topological material. Furthermore, our study reveals that the temperature dependence of the lattice thermal conductivity exhibits a double-peak structure, suggesting strong phonon–phonon or phonon–defect interactions in this material, characteristic of resonant scattering. Furthermore, these findings will help in a better understanding of thermal transport in advanced materials and devices at the micro scale.

Physics - Condensed matter physics↗
DOE OSTI2025

Beyond Magic Barrels: Digital manufacturing for crystallization, process development and optimization of explosive materials: Part II Resveratrol Exemplar

This SAND report summarizes work supported by an Engineering Sciences Research Foundation (ESRF) Lab Directed Research and Development (LDRD) project entitled “Beyond Magic Barrels: Digital manufacturing for crystallization, process development and optimization of explosive materials.” This SAND report is written in two parts with Part 1 discusses recrystallization of our explosive exemplar and Part 2 summarizing our work with recrystallization of resveratrol. We have studied resveratrol recrystallization with a multiscale approach combining experiments, modeling and simulation. At the single crystal scale, microscopy experiments illuminate crystal time-dependent growth rates using advanced image analysis. Bench scale experiments were carried out to look at growth of multiple particles in a small reactor creating thousands of particles and analyzing the results with microscopy and μCT. For the modeling we combine kinetic Monte Carlo (kMC) models with subscale information from density functional theory (DFT) or molecular dynamics. This work is discussed in Part 1 and can also be found in a paper from the project discussing a coarse-grained kMC model specifically developed for resveratrol. For well-mixed systems, we have population balance equations (PBE) linked with species mass conservation forming a set of ordinary differential equations that can be solved quickly. For more complicated geometries, such as the vat crystallization used throughout the complex, a coupled computational fluid dynamic (CFD)/PBE method was developed to account for gradients in temperature and concentration and differences in crystallization rates throughout the domain. These simulations are more complex and require high performance computing. We present results for two cases: 5% seed fast cool with parameters fit to the well-mixed case and 5% seed slow cool using the same parameters. We show reasonable agreement with experiments though are particles are significantly larger than the experiments.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Statistical evaluation of microscale stress conditions leading to void nucleation in the weak shock regime

Here, we investigate the heterogeneity of the stress state driven by anisotropic deformation response at the single crystal level through five statistical volume element (SVE) calculations of polycrystalline BCC tantalum. This work focuses on grain boundaries as a prominent material defect type prone to void nucleation based upon experimental observations of predominantly intergranular void nucleation in this material. The SVEs are constructed to be statistically representative of larger volumes of material and are meshed such that mean and standard deviation of grain size and orientation information is reconstructed. The computational meshes feature hexahedral (brick) elements and smooth conformal grain boundaries where significant stress concentration is known to occur, a tail effect of interest in the extreme events process of dynamic ductile damage. An existing micromechanical crystallographic plasticity model shown to capture the single crystal behavior of BCC tantalum well is used to perform the polycrystal calculations. The model includes representation of the non-Schmid effect of non-planar screw dislocation kinetics in tantalum. A three-dimensional stress state time profile predicted by damage modeling of a flyer plate impact experiment is applied as boundary conditions to each SVE. Resulting grain boundary stress state statistics are strongly non-Gaussian. Significant structural evolution is observed within the compressive hold before unloading into tension in the stress profile. Strong angular dependence of grain boundary traction magnitude with shock direction is observed. Non-Schmid effects continue to suggest their influence on propensity of microstructural defect types to nucleate voids. A general void nucleation criterion is proposed using probability theory. The general framework is specified to polycrystalline BCC tantalum in the weak shock regime to include the SVE calculations and literature molecular dynamics calculations of grain boundary void nucleation strength. Probability density functions (PDFs) are used to describe the interaction between the local stress state heterogeneity and the distributed grain boundary void nucleation strength state. A causation entropy maximization procedure removes the requirement for ad hoc selection of a PDF functional form and provides a rigorous procedure for data-based PDF determination. The resulting physically informed PDF describes the spatial appearance frequency of nucleated voids as a function of applied macroscale pressure. Lower length scale physics are thus packaged in a precise and computationally efficient way to provide computational plasticity insight to macroscale dynamic ductile damage models.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Inelastic deformation of diamond single crystals shock compressed to multimegabar stresses: Wave profile calculations

As the archetypal strong solid, the response of diamond shock compressed to multimegabar stresses is important for fundamental science and for numerical simulations of wave profiles for applications in high energy density physics experiments. Previous experiments and analysis have shown that the commonly used hydrodynamic assumption is invalid for diamond shock compressed to stresses below melt and an elastic–inelastic description is needed. Here, we present a phenomenological material model for calculating wave profiles in shock compressed diamond single crystals that incorporates this description. Also, to support the modeling effort, we carried out wave profile measurements on shock compressed diamond single crystals at the Sandia Z facility to augment previous measurements. Wave profiles for [100] and [111] diamond calculated using the material model provide a good match to the elastic–inelastic response (observed two-wave structure) measured at ∼325 and ∼360 GPa. Furthermore, the calculated peak stresses for single (overdriven) waves provide a good match to the measured Hugoniot states for stresses reaching ∼700 GPa, which is near melting conditions. The present results show that the diamond single crystal response at multimegabar shock stresses is characteristic of a brittle solid—pressure-dependent strength and strength loss due to inelastic deformation.

Deformation↗