Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Elastic energy”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Nanocrystal Geometry Governs Phase Transformation Pathways in Palladium Hydride

Pathways and structural dynamics of phase transformations impact performance of materials in energy and information storage technologies. Palladium hydride (PdH x ) nanocrystals are an ideal model system for studying solute-induced phase transformations, where elastic energy from lattice mismatch between α-PdH x and β-PdH x phases is often considered a key to determining the transformation pathways. α/β-PdH x interfacial elastic energy is affected by the confined geometry of a nanocrystal. However, how nanocrystal geometry influences phase transformation pathways is largely unknown. Using in situ liquid phase transmission electron microscopy, we directly visualize hydrogenation in Pd nanocrystals with two geometries, a nanocube and a hexagonal nanoplate. Both follow similar sequences of an initially curved nucleus, interface flattening, and reverse-stage nucleation; however, their evolving α/β-PdH x interfaces exhibit geometry-dependent crystallographic alignments. In nanocubes, {100}-aligned configurations conform to static elastic energy ordering, representing a pathway that maintains a local mechanical equilibrium, whereas nanoplates display both {110}- and {211}-aligned interfaces. Theoretical simulations show that geometry determines the accessibility of alternative phase transformation pathways as the system is driven far from equilibrium during hydrogenation. In conclusion, these findings identify geometry as a fundamental parameter for directing phase transformation pathways, offering design principles for accessing atypical configurations and improving properties of intercalation-based devices.

Lee, Daewon [Lawrence Berkeley National Laboratory↗
DOE OSTI2024

Novel ceramic capacitors with ultrahigh energy density and efficiency (Final Technical Report)

Antiferroelectric ceramics are a special class of material that have shown great potential as the dielectric in electrical capacitors due to their high energy- and power-density. During each charge-discharge cycle, the ceramic undergoes transformation to a ferroelectric phase and resumes its antiferroelectric phase. The hysteresis associated with the transitions leads to a mediocre energy efficiency and service lifetime of antiferroelectric capacitors and, hence, their almost absence in commercial products. Under the support of this research project, we first formulated a universal lattice-compatibility theory that included electrostatic polarization energy along with elastic energy and thermal energy to understand the origin of the hysteresis in antiferroelectric oxides. Guided by this compatibility theory, we conducted high-throughput density functional theory (DFT) calculations to assess chemical modifiers and their effect on crystal structures of 400+ PbZrO 3 -based compositions. Down-selected compositions were experimentally validated for their suppressed hysteresis and higher energy efficiency. The verified low-hysteresis compositions were then expanded to an antiferroelectric ceramic library with nearly 500 new compositions (more than 1,500 samples) using high-throughput experiments involving ceramic synthesis and property screening. The large quantity of data generated (theory and experimental) in these tasks were processed by machine-learning techniques and identified trends were fed to the next iteration. In the end, we successfully discovered four compositions with near-zero hysteresis, yielding a world-record energy efficiency of 98.2% at an energy density of 3.0 J/cm 3 . Furthermore, our antiferroelectric ceramic capacitor reaches 79.5 million charge-discharge cycles lifetime, a factor of 80 enhancement over previous antiferroelectric ceramics with large hysteresis. These research accomplishments have not only met the milestones set in the SOPO, but also led to two patent filings, three journal publications (one of them was in Advanced Materials, impact factor 29.4), and nine oral presentations at various venues. Through the course of the project, three postdocs, four Ph.D. students, and one M.S. student were trained. In short, our project established a new methodology in searching next-generation functional ceramics on the fundamental side and discovered several high-efficiency antiferroelectric compositions for capacitors on the applied side. Once fabricated into the multilayer form for commercial applications, these ceramic capacitors can potentially enable the high temperature high power density DC-link capacitors that are critical for the next generation inverters in electric vehicles. The project also significantly contributed to the nation’s workforce development in the STEM fields.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Towards a Deeper Fundamental Understanding of (Al,Sc)N Ferroelectric Nitrides

Density functional theory (DFT) calculations, within the virtual crystal alloy approximation, are performed, along with the development of a Landau-type model employing a symmetry-allowed analytical expression of the internal energy and having parameters determined from first principles, to investigate properties and energetics of Al1-xScxN ferroelectric nitrides in their hexagonal forms. These DFT computations and this model predict the existence of two different types of minima, namely, the fourfold-coordinated wurtzite (WZ) polar structure and a five-fold coordinated paraelectric hexagonal phase (denoted as H5), for any Sc composition up to 40%. The H5 minimum progressively becomes the lowest-energy state within hexagonal symmetry as the Sc concentration increases from 0 to 0.4. Furthermore, the model points to several key findings. Examples include the crucial role of the coupling between polarization and strains to create the WZ minimum, in addition to polar and elastic energies, and that the origin of the H5 state overcoming the WZ phase as the global minimum within hexagonal symmetry when increasing the Sc composition mostly lies in the compositional dependency of only two parameters-one linked to the polarization and another one being purely elastic in nature. Other examples are that forcing Al1-xScxN systems to have no or a weak change in lattice parameters when heating them allows us to reproduce their finite-temperature polar properties well and that a value of the axial ratio close to that of the ideal WZ structure implies a large polarization at low temperatures but not necessarily at high temperatures because of the ordered-disordered character of the temperature-induced formation of the WZ state. Such findings should allow for a better fundamental understanding of (Al,Sc)N ferroelectric nitrides, which may be used to design efficient devices having, e.g., low operating voltages.

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

Reference lattice, sound, stiffness, and magnetic transitions of Ising monolayers

A reference lattice, away from which elastic distortions induced by the spin texturing of two-dimensional (2D) magnets take hold, is motivated from a picture of pairwise Biot-Savart interactions among identical solenoids that either elongate or compress a (“zero-current”) spring lattice. Applied to a paradigmatic CrSiTe 3 monolayer, the reference is given by the average between the atomic positions of ferromagnetic (FM) and Néel antiferromagnetic (AFM) lattices; such an atomic disposition permits understanding structural distortions and elastic energies due to magnetism readily. Furthermore, the anisotropic speed of sound in the magnetic ground state explains an observed anisotropy of vibrational frequencies on similar magnets. Elastic stiffness constants are reported, too. Magnetic energies in four Ising structural configurations were calculated and the strain needed for those 2D magnets to undergo an AFM to FM quantum phase transition was determined as well.

density functional theory↗
DOE OSTI2025

Phase-field modeling of orientation-dependent crack growth in ductile single crystals with anisotropic elasticity

Crack growth in ductile single crystals (DuSCs) is orientation dependent due to the anisotropies of crystal plasticity and elastic tensor. This study develops a phase-field model incorporating both crystal plasticity and crack growth and proposes a general method to decompose the elastic energy into compressive and tensile parts to prevent crack growth under compression in the phase-field description. The phase-field model, in combination with three Euler angles, is employed to simulate orientation-dependent crack growth in DuSCs. The contributions from crystal plasticity and anisotropic elasticity are compared, and the former is found to dominate in the anisotropy of crack growth in copper single crystals. Furthermore, the simulation results demonstrate that crystal orientation strongly affects the heterogeneous distribution of plastic strain and the interaction between plastic strain and crack growth. High-throughput phase-field simulations are performed with exhaustive crystal orientations, and the results are explained based on the anisotropy of the Taylor factor.

Computational Solid Mechanics↗
DOE OSTI2026

Computationally efficient method for determining limiting velocities of edge dislocations in anisotropic crystals

The continuum-limit theory of dislocations in crystals predicts divergences in the elastic energy at crystal-geometry dependent limiting velocities vL, which separate subsonic, transsonic, and supersonic dislocation glide regimes and are therefore import for material strength models at high strain rates. Although it is known how to calculate those limiting velocities, there is one special case - edge dislocations with reflection symmetry, but non-vanishing elastic constants c16 or c26 - where previous methods have been notoriously slow. In this letter, we address this deficiency by deriving a computationally efficient method for determining the limiting velocities of edge dislocations with reflection symmetry which is two orders of magnitude faster than the previous method.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Realizing mechanical frustration at the nanoscale using DNA origami

Structural designs inspired by physical and biological systems have been previously utilized to develop mechanical metamaterials with enhanced properties based on clever geometric arrangement of constituent building blocks. Here, we use the DNA origami method to realize a nanoscale metastructure exhibiting mechanical frustration, a counterpart of the well-known phenomenon of magnetic frustration. By selectively actuating reconfigurable struts, it adopts either frustrated or non-frustrated states, each characterized by distinct free energy profiles. While the non-frustrated state distributes the strain homogeneously, the frustrated mode concentrates it at a specific location. Molecular dynamics simulations reconcile the contrasting behaviors and provide insights into underlying mechanics. We explore the design space further by tailoring responses through structural modifications. Our work combines programmable DNA self-assembly with mechanical design principles to overcome engineering limitations encountered at the macroscale to design dynamic, deformable nanostructures with potential applications in elastic energy storage, nanomechanical computation, and allosteric mechanisms in DNA-based nanomachinery.

DNA nanostructures↗
DOE OSTI2026

Scattering Processes from Quantum Simulation Algorithms for Scalar Field Theories

We provide practical simulation methods for scalar field theories on a quantum computer that yield improved asymptotics as well as concrete gate estimates for the simulation and physical qubit estimates using the surface code. We achieve these improvements through two optimizations. First, we consider a finite volume approach for estimating the elements of the S-matrix. This approach is appropriate in general for 1+1D and for certain low-energy elastic collisions in higher dimensions. Second, we implement our approach using a series of different fault-tolerant simulation algorithms for Hamiltonians formulated both in the field occupation basis and field amplitude basis. Our algorithms are based on either second-order Trotterization or qubitization. The cost of Trotterization in occupation basis scales as O ( λ N 7 | Ω | 3 / ( M 5 / 2 ϵ 3 / 2 ) ) where λ is the coupling strength, N is the occupation cutoff, | Ω | is the volume of the spatial lattice, M is the mass of the particles and ϵ is the uncertainty in the energy calculation used for the S -matrix determination. Qubitization in the field basis scales as O ( | Ω | 2 ( k 2 Λ + k M 2 ) / ϵ ) , where k is the cutoff in the field and Λ is a scaled coupling constant. We find in both cases that the bounds suggest physically meaningful simulations can be performed using on the order of 4 × 10 6 physical qubits and 10 12 T -gates which corresponds to roughly one day on a superconducting quantum computer with surface code and a cycle time of 100 ns. This places the simulation of scalar field theory within striking distance of the gate counts for the best available chemistry simulation results.

Hardy, Andrew [Toronto U.] (ORCID:0000000235817382↗
DOE OSTI2025

Phase-field modeling and experiments of dynamic fracture in single crystal quartz

Predicting the onset and characteristics of brittle fracture is important for a wide range of engineering and geological material applications. In this paper, we study important aspects of brittle fracture in α-quartz by phase-field modeling and experiments using a top-down approach. In the modeling framework, the work term in the Griffith energy balance is replaced with internal energy contributions that represent surface energy, thermal energy, and elastic strain energy stored in defects. This allows parametrization of individual energy contributions in terms of internal state variables and keeps track of energy partitioning after the onset of fracture. The path and history dependence of fracture is included in evolution laws for internal state variables, e.g., entropy evolution, while the energy remains a true potential. In the experimental part, dynamic compression experiments coupled with X-ray phase contrast imaging are performed on cube-like samples with a hole. In the top-down analysis, dynamic compression and three point bending experiments from the literature are simulated with the developed phase-field damage model. In conclusion, the fitted model highlights the strain rate, size, and stress state dependence of damage nucleation and evolution in single crystal α-quartz.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Temperature and dose rate effects on dislocation loop formation and decoration in self-ion irradiated high-purity Fe-Cr model alloys

Body-centered cubic (BCC) Fe-based alloys are promising candidate materials for advanced nuclear reactors. A widespread perception is that the addition of Cr favors the formation of ½ <111> dislocation loops. However, the Cr effects on dislocation loop formation/evolution remain unclear due to the limited experimental studies. In this study, we irradiate ultra-high purity bulk BCC Fe-xCr (x = 3, 5, and 8 wt%) with 6.7-8 MeV Fe ions over a wide range of temperatures (Tirr = 250-500 °C) and dose rates (10⁻⁵ to 10⁻³ dpa/s) to two midrange doses (0.35 and 3.5 dpa). The major irradiation-induced microstructures are loop decoration along non-screw dislocation lines at low temperatures and high dose rates as well as self-assembled loops at elevated temperatures. Post-irradiation Burgers vector analysis via g•b method on loops and electron tomography reveals that <001> preponderance is observable in all irradiated Fe-Cr samples at elevated temperature, resulting from stronger elastic energy originating from <001> edge.

Li, Yao [University of Tennessee (UT)]↗
DOE OSTI2024

Insulating moiré homobilayers lack a threefold symmetric second-harmonic generation

Atoms within moiré bilayers relax in plane to minimize elastic energy; such relaxation brings their space group symmetries down to P1. Here, the ab initio second harmonic generation (SHG) of twisted and atomistically optimized hBN bilayers was determined at four twist angles (θ = 38.21°, 60.00°, 73.17°, and 98.21°) and for three displacements τ measured away from the ground state AA' configuration. All moiré bilayers have a P1 space symmetry after structural optimization. This situation is quite different to monolayers with hexagonal lattices, which retain a threefold symmetry. We point out that the actual symmetries of the SHG reported for hBN bilayers on two experimental works do not coincide with the sixfold symmetric theoretical profiles they provide [either sin 2 ⁡(3⁢Φ) or cos 2 ⁡(3⁢Φ)], and show that the intrinsic low structural symmetry of (atomically optimized) hBN bilayer moirés can in fact be read out from experimental SHG intensity profiles—which are tunable by θ and by the frequency of light ω: The SHG is most definitely not sixfold symmetric because moirés do not retain a threefold symmetry. Furthermore, an extrinsic twofold symmetry of the SHG emission is realized by tilting the pump by an angle α away from the 2D material's normal, regardless of θ and ω. Furthermore, the design of in-plane and ultrathin sources of SHG with low symmetry could be useful for the eventual creation of entanglement sources from 2D materials.

2-dimensional systems↗
DOE OSTI2025

Enhancing stability, magnetic anisotropy, and coercivity of manganese aluminum: Machine learning, ab initio , and micromagnetic modeling

The binary manganese aluminum (MnAl) alloy with L⁢1 0 crystal structure is a promising rare earth (RE) element-free permanent magnetic material because of its exceptional magnetic properties. However, experimentally synthesizing it in a stable bulk form is extremely challenging. Here, in this study, an alternative method of stabilizing the material, a pathway for experimental synthesis and validation, is proposed and theoretically verified. This is done by partially substituting Mn and Al sites with Fe and Ni and identifying its enhanced phase stability, saturation magnetization density, magnetic anisotropy, and coercivity from density functional theory (DFT), machine learning (ML) crystal graph convolution neural network (CGCNN), and micro-magnetic modeling. When considering a fixed 50% Ni, the magnetic anisotropy increases with the increasing Fe content but decreases the formation energy. The calculated formation energies, elastic constants, and phonon frequencies demonstrate that the binary and quaternary compositions are stable. Most importantly, in 50% Fe and Ni-substituted-equiatomic phase, magnetic anisotropy constants and saturation magnetization density increase by 56% and 23% as compared to the MnAl. Further, the coercivity of the equiatomic phase predicted with micro-magnetic modeling is higher by 17% than the parent compound.

Bhandari, Churna [Ames National Laboratory, and Io↗
DOE OSTI2025

Niobium hydride formation in superconducting qubit thin films

The formation of nonsuperconducting hydrides in 160–170-nm-thick films of niobium is examined. We identify six elastically distinct orientation relationships ɛ−Nb 4 ⁢H 3 takes within the matrix, solid-solution 𝛼−NbH 𝑥 phase. We employ a phase field model to assess the impact of elastic energy induced by the strain of phase transformation on the morphology and transformation dynamics of ɛ precipitates within a thin film. We consider the dimensions of the thin film, crystallographic growth direction, and diffusion rates to predict the timescale of hydride evolution. Leveraging the finite element method, we predict two-dimensional and three-dimensional equilibrium shapes of ɛ-hydrides within a bulk sample and in a thin film that has a traction-free surface. Our results suggest that niobium hydrides migrate to the free surface of the film. Precipitates which reach the free surface coarsen, while precipitates within the film dissolve. Precipitates in both two dimensions and three dimensions experience a repulsive interaction force at the free surface, that is attractive in the bulk, shown in experiment and theory of previous studies.

36 - MATERIALS SCIENCE↗
DOE OSTI2025

Machine learning interatomic potential for predicting the thermal properties of uranium nitride

We present a combined computational and experimental investigation of the thermal properties of uranium nitride (UN), focusing on the development of a machine learning interatomic potential (MLIP) using the moment tensor potential framework. The MLIP was trained on density functional theory (DFT) data and validated against various quantities including energies, forces, elastic constants, phonon dispersion, and defect formation energies, achieving excellent agreement with DFT calculations, prior experimental results, and our thermal conductivity measurement. The potential was then employed in molecular dynamics simulations to predict key thermal properties such as melting point, thermal expansion, specific heat, and lattice thermal conductivity. To further assess model accuracy, we fabricated a UN sample and performed new thermal conductivity measurements representative of single-crystal properties, which showed strong agreement with the MLIP predictions. This work confirms the reliability and predictive capability of the developed potential for determining the thermal properties of UN.

36 - MATERIALS SCIENCE↗
DOE OSTI2025

Coupled phase field damage and crystal plasticity analysis of intragranular fracture: The role of crystallographic orientation and voids

Damage evolution in engineering metal alloys at the grain scale exhibits significant microstructural heterogeneity and anisotropy. These heterogeneities create local hotspots for stress and strain localization, leading to void nucleation. Crystal orientation influences the active slip systems around voids, affecting lattice rotation and potentially forming discontinuities. At low triaxiality, voids may change shape due to lower stress, rotation, elongation, and coalescence. At high triaxiality, the correlation between crystal orientation and void growth rate becomes stronger, resembling the behavior observed in isolated single crystals. Therefore, understanding the effects of crystal orientation, heterogeneous strain, and defect evolution is crucial for single crystal fracture characterization. Here, in this work, a coupled phase-field damage (PFD) and crystal plasticity (CP) model is implemented within a finite element framework to analyze crystal deformation and failure. The CP method employs a dislocation density-based constitutive model, while intragranular failure is modeled using an anisotropic PFD method. The PFD model considers both the stored energy due to elastic stretching and the energy release due to defect formation and crack formation. A single crystal Al2219 with an intracrystalline spherical void is chosen to analyze fracture. The study finds that fracture propagation is strongly correlated with crystal orientations. This coupled CP-PFD model provides accurate failure prediction in crystalline materials by incorporating the effects of crystal orientations and existing voids. This study demonstrates how the local microstructure and defects influence plastic deformation and failure mechanisms in metal alloys.

Aluminum alloy↗
DOE OSTI2024

Multiscale modeling of metal-hydride interphases—quantification of decoupled chemo-mechanical energies

Abstract The quantification of interphase properties between metals and their corresponding hydrides is crucial for modeling the thermodynamics and kinetics of the hydrogenation processes in solid-state hydrogen storage materials. In particular, interphase boundary energies assume a pivotal role in determining the kinetics of nucleation, growth, and coarsening of hydrides, alongside accompanying morphological evolution during hydrogenation. The total interphase energy arises from both chemical bonding and mechanical strains in these solid-state systems. Since these contributions are usually coupled, it is challenging to distinguish via conventional computational approaches. Here, a comprehensive atomistic modeling methodology is developed to decouple chemical and mechanical energy contributions using first-principles calculations, of which feasibility is demonstrated by quantifying chemical and elastic strain energies of key interfaces within the FeTi metal-hydride system. Derived materials parameters are then employed for mesoscopic micromechanical analysis, predicting crystallographic orientations in line with experimental observations. The multiscale approach outlined verifies the importance of the chemo-mechanical interplay in the morphological evolution of growing hydride phases, and can be generalized to investigate other systems. In addition, it can streamline the design of atomistic models for the quantitative evaluation of interphase properties between dissimilar phases and allow for efficient predictions of their preferred phase boundary orientations.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Controlled Shifts of X-ray Emission Lines Measured with Transition Edge Sensors at the Advanced Photon Source

The measurement of shifts in the energy of X-ray emission lines is important for understanding the electronic structure and physical properties of materials. In this study, we demonstrate a method using a synchrotron source to introduce controlled eV-scale shifts of a narrow line in between fixed-energy fluorescence lines. We use this to characterize the ability of a hard X-ray superconducting Transition Edge Sensor (TES) array to measure line shifts. Fixed fluorescence lines excited by higher harmonics of the monochromatic X-ray beam are used for online energy calibration, while elastic scattering from the primary harmonic acts as the variable energy emission line under study. We use this method to demonstrate the ability to track shifts in the energy of the elastic scattering line of magnitude smaller than the TES energy resolution, and find we are ultimately limited by our calibration procedure. The method can be applied over a wide X-ray energy range and provides a robust approach for the characterization of the ability of high-resolution detectors to detect X-ray emission line shifts, and the quantitative comparison of energy calibration procedures.

Guruswamy, Tejas (ORCID:0000000241656765)↗
DOE OSTI2025

Stability, electronic quantum states, and magnetic interactions of Er 3+ ions in Ga 2 ⁢O 3

Here, we report an ab initio study of phase stability, defect formation, electronic structure, and multiple magnetic, Dzyaloshinskii-Moriya, optical, hyperfine, and crystal field interactions in erbium (Er)-doped wide band gap 𝛼- and 𝛽-gallium oxides (Ga 2 ⁢O 3 ), critically important to make a foundation for both optoelectronic and quantum information applications. The chemical, structural, mechanical, and dynamical stabilities of the pristine phases are confirmed from respective negative formation energies, negative cohesive energies, favorable elastic constants, and positive phonon frequencies. The phonon dispersions indicate that the Ga-O bonds are uniform in the 𝛼-phase, while they vary in the 𝛽-phase due to the anisotropic polyhedral movement. The defect formation energy analysis confirms that both Er-doped 𝛼- and 𝛽−Ga 2 ⁢O 3 prefer Er 3+ (neutral) state. The underestimated band gaps of the pristine phases from standard density functional theory (DFT) calculations as compared to experimental values are corrected by employing the hybrid functional calculations, resulting in the indirect band gaps of 5.21 eV in 𝛼−Ga 2 ⁢O 3 and 4.94 eV in 𝛽−Ga 2 ⁢O 3 . The site preference energy analysis indicates partial occupation of Er in the octahedral site of Ga. The anisotropic nature of hyperfine tensor coefficients of Er is similar in both phases, which may be due to the occupation of Er in the same octahedral Ga site. On the other hand, the calculated magnetic exchange interaction between two Er dopants is negative for 𝛼 and positive for 𝛽, indicating an antiferromagnetic (AFM) ground state in the former and a ferromagnetic (FM) ground state in the latter. Large values of Dzyaloshinskii-Moriya interactions (DMIs) are obtained along the 𝑥 direction in the 𝛼 and along the 𝑦 direction in the 𝛽. The large DMI may support exotic magnetic textures, a promising direction for spintronic applications. The analysis of dielectric constants and refractive indices of both pristine and Er-doped phases shows a good agreement with available experimental values. The calculated optical anisotropy is slightly higher in 𝛽 than those in 𝛼, which is due to the involvement of lower symmetry in 𝛽. The crystal field coefficients (CFCs) calculated from DFT are used to analyze 4⁢𝑓 multiplets and 4⁢𝑓 −4⁢𝑓 transitions. Thus calculated lowest energy level of the first excited state to the lowest energy level of the ground state is about 1.53 µ⁢m, which is in a good agreement with available experiments, and it falls within the quantum telecommunication wavelength range.

3-dimensional systems↗