Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Long duration energy storage”

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.

847 records · Page 12
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Thermally induced mimicry of quantum cluster excitations and implications for the magnetic transition in FePSe 3

In two dimensional magnets, the interplay of thermal fluctuations and spin anisotropy control the existence of long-range magnetic order. In the van der Waals antiferromagnets FePX 3 , orbital degeneracy in the 𝑡 2⁢𝑔 levels of the Fe 2+ ions in octahedral coordination yields strong uniaxial anisotropy, which stabilizes magnetic order up to 𝑇 ≈ 100 K. Recent inelastic neutron scattering measurements around the magnetic ordering transition have shown the existence of a broad spectrum of magnetic fluctuations with nontrivial momentum dependence, which has been interpreted as evidence for localized entangled cluster excitations. In this paper, we offer an alternative interpretation using classical nonlinear spin dynamics simulations. We present stochastic Landau Lifshitz dynamics simulations that reproduce the neutron scattering measurements of Chen et al. [npj Quantum Mater. 9, 40 (2024)] on FePSe 3 . These calculations faithfully explain the dynamical structure factor's momentum and energy dependence and point to a classical origin for the excitations observed in neutron spectroscopy and that the order-disorder transition can be understood in terms of thermal fluctuations overcoming the anisotropy energy.

Landau-Lifschitz-Gilbert equation↗
DOE OSTI2023

Effect of the Nature of Both Cation and Anion Substitution on the Structural Symmetry of Li‐Rich 3 d ‐Metal Chalcogenide Electrodes

Abstract Li‐rich layered chalcogenides have recently led to better understanding of the anionic redox process and its associated high capacity while providing ways to overcome its practical limitations of voltage fade and irreversibility. This study reports on the feasibility of triggering anionic activity in Li 2 TiS 3 , through anionic substitution (Se for S) or cationic substitution (Fe for Ti). Herein, the chalcogenide chemical space is further explored to prepare mono‐substituted Li 1.7 Ti 0.85 Mn 0.45 Ch 3 (Ch = S/Se) and doubly substituted cationic and anionic phases (Li 1.7 Ti 0.85 Fe 0.45 S 3‐z Se z ) which crystallize either in the O3‐ or O1‐type structures depending upon substituents. All series show a bell‐shape capacity variation as function of the transition metal (TM) substitution degree with values up to 240 mAh g −1 . For specific compositions, a structural O3 to O1 phase transition is observed upon Li removal, which is not reversible upon Li re‐insertion due to kinetic limitations and negatively affects long‐term cycling performance. Density functional theory (DFT) calculations confirm the O3/O1 relative stability along the different series and point subtle electronic differences in the TM‐doping, rationalizing the structural and electrochemical behaviors of these phases upon cycling. These findings provide further insights into the link between structural and electronic stability, which is of key importance for designing chalcogenide‐based anionic redox compounds.

Chemistry↗
DOE OSTI2025

Interaction of magnetic fields with spinons in a fractionalized state

The 4d-electron trimer lattice Ba₄Nb₁₋ₓRu₃₊ₓO₁₂ exhibits either a quantum spin liquid (QSL) or a heavy-fermion strange metal (HFSM) phase, depending on Nb content. In the QSL state, itinerant spinons act as effective heat carriers, enhancing thermal conductivity. Strikingly, applying a magnetic field up to 14 T causes an abrupt, up-to-5000% increase in heat capacity below 150 mK, disrupting the linear temperature dependence typical of both phases. Meanwhile, AC susceptibility and electrical resistivity remain nearly unchanged, while thermal conductivity drops by up to 40% below 4 K. These results suggest spinons, despite being charge-neutral, are highly sensitive to magnetic fields at low temperatures. We propose that the magnetic field could induce Anderson localization of spinons, creating emergent non-magnetic two-level systems responsible for the Schottky-like anomaly in heat capacity. These findings point to a previously unexplored regime of spinon dynamics, potentially governed by field-induced localization and distinct from conventional magnetic or transport signatures.

Condensed-matter physics↗
DOE OSTI2026

Heavy Element Spectroscopy in the Gas Phase

Actinides are inherently unstable and undergo nuclear decay processes with a concurrent release of energy. Consequently, they are used for nuclear power generation, nuclear weapons, and nuclear medicine. However, the radioactive decay processes also pose significant technological problems for the safe treatment and storage of spent nuclear materials. Cost-effective extraction of the actinides is the key first step in the remediation of nuclear waste, but the appropriate chemical means have yet to be determined. Our present understanding of the chemistry of actinides is limited, with the role of the 5f electrons posing a set of particularly challenging questions. The work reported here is focused on the use of electronic spectroscopy to probe the bonding of small molecules in the gas phase that contains thorium or uranium. Analyses of these data, carried out within the framework of ligand field theory, reveal clear evidence that the 5f electrons are spectators that retain their atomic metal ion character.

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

Emergence of low-energy spin waves in superconducting electron-doped cuprates

In order to fully utilize the technological potential of unconventional superconductors, an enhanced understanding of the superconducting mechanism is necessary. In the best performing superconductors, the cuprates, superconductivity is intimately linked with magnetism, although the details of this coupling remain elusive. Here, we address this gap by studying the electron-doped cuprate Nd 1.85 Ce 0.15 CuO 4−δ that has an antiferromagnetic ground state when synthesized and only becomes superconducting after a reductive annealing process. Using neutron spectroscopy, we show that the as-grown crystal exhibits a large spin pseudogap in the magnetic fluctuation spectrum. Annealing removes defects introduced by the commonly employed synthesis method and significantly reduces the spin pseudogap. While the spin pseudogap in the annealed sample likely arises from superconductivity, in the as-grown sample it results from the absence of long-wavelength spin waves. These results reveal a direct connection between defects, magnetism, and superconductivity, offering new insight into the mechanisms underlying high-temperature superconductivity and guiding the design of improved superconducting materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

Atomic Evolution of Hydrogen Intercalation Wave Dynamics in Palladium Nanocrystals Revealed by Liquid-Phase Transmission Electron Microscopy

Solute-intercalation-induced phase separation creates spatial heterogeneities in host materials, a phenomenon ubiquitous in batteries, hydrogen storage, and other energy devices. Despite many efforts, probing intercalation processes at the atomic scale has been a significant challenge. By utilizing liquid-phase transmission electron microscopy (TEM), we study hydrogen (de)intercalation in palladium nanocrystals as a model system and have achieved unprecedented atomic-resolution imaging of hydrogen intercalation wave dynamics. Our observations reveal that intercalation wave mechanisms, instead of shrinking-core mechanisms, prevail at ambient temperature for palladium nanocubes ranging from ∼60 nm down to ∼10 nm. Systematic image analysis uncovers the atomic evolution of the hydrogen intercalation wave, transitioning from nonplanar and inclined boundaries to those closely aligned with {100} planes. Our kinetic Monte Carlo simulations demonstrate that the observed intercalation wave dynamics correspond to sorption pathways minimizing the lattice mismatch strain at the phase boundary. In conclusion, unveiling the atomic intercalation pathways holds profound implications for engineering intercalation-mediated devices and advancements in energy sciences.

Lee, Daewon [Lawrence Berkeley National Laboratory↗
DOE OSTI2024

Fluorophilic Sigma(σ)‐Lock Self‐Healable Copolymers

Although F-Containing molecules and macromolecules are often used in molecular biology to increase the binding with Lewis acidic groups by introducing favorable C−F dipoles, there is virtually no experimental evidence and limited understanding of the nature of these interactions, especially their role in synthetic polymeric materials. These studies elucidate the molecular origin of inter- and intra-Chain interactions responsible for self-healing of F-Containing copolymers composed of pentafluorostyrene and n-butyl acrylate units (p(PFS/nBA). Guided by dynamic surface oscillating force (SOF) and spectroscopic measurements supported by molecular dynamics (MD) simulations, these studies show that the reformation of σ-σ orbitals in −C−F of PFS and CH 3 CH 2 − of nBA units enables the recovery of entropic energy via fluorophilic-σ-lock van der Waals forces when PFS/nBA molar ratios are ~50/50. The strength of these interactions determined experimentally for self-healable PFS/nBA compositions is in the order ~0.3 kcal/mol which primarily comes from fluorophilic-σ-lock (~70 %) contributions. These interactions are significantly diminished for non-self-healable counterparts. Strongly polarized −C−F σ orbitals create lateral dipolar forces enhancing the affinity towards −C−H orbitals, facilitating energetically favorable interactions. Entropic recovery driven by non-Covalent bonding offers a valuable tool in designing materials with unique functionalities, particularly self-healable batteries and energy storage devices.

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

Equation of state for Hf, Ta, W, Re, Os, Ir, Pt, and Au to multi-terapascal pressures from density-functional theory

We present the zero-temperature equation of state (pressure dependence of compression) and phase stability predictions for the 5d-transition metals obtained from all-electron density-functional theory (DFT) calculations. The results compare favorably with experiments but extend beyond current experimental capabilities to 10 TPa. Our study reveals phase changes that are explained from the calculated electronic structure. The cubic face-centered and body-centered structures (fcc and bcc), together with two-, three-, and four-layered hexagonal structures, play major roles under compression. The results’ dependence on the electron exchange and correlation in the DFT approach is investigated, and it is shown that the impact of the choice, while significant at lower pressures, diminishes in the terapascal regime. We further illustrate that the normal parabolic trends in atomic volume and bulk modulus with atomic number, due to the occupation of bonding and anti-bonding 5d states, break down at TPa pressures, suggesting drastically different chemical bonding at these extreme conditions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Time-Resolved Stochastic Dynamics of Quantum Thermal Machines

Steady-state quantum thermal machines are typically characterized by a continuous flow of heat between different reservoirs. However, at the level of discrete stochastic realizations, heat flow is unraveled as a series of abrupt quantum jumps, each representing an exchange of finite quanta with the environment. Here, in this work, we present a framework that resolves the dynamics of quantum thermal machines into cycles classified as enginelike, coolinglike, or idle. We analyze the statistics of individual cycle types and their durations, enabling us to determine both the fraction of cycles useful for thermodynamic tasks and the average waiting time between cycles of a given type. Central to our analysis is the notion of intermittency, which captures the operational consistency of the machine by assessing the frequency and distribution of idle cycles. Our framework offers a novel approach to characterizing thermal machines, with significant relevance to experiments involving mesoscopic transport through quantum dots.

full counting statistics↗
DOE OSTI2025

Highly textured metal anodes for stable aqueous batteries: Fabrication and characterization

We report a purely mechanical “cold-compression flow” method for fabricating Zn, Sn, and In substrates with tunable crystallographic textures. Using textured Zn as a model system, we investigate Zn electrocrystallization and demonstrate correlated growth of crystalline films with correlation lengths from tens to hundreds of micrometers. At 5 milliamperes per square centimeter (mA/cm 2 ), capacities between 20 and 82 milliampere hours per square centimeter (mA·hour/cm 2 ) are achieved depending on substrate texture level. At higher currents (40 mA/cm 2 ), capacities reach up to 604 mA·hour/cm 2 . Rotating disk electrode studies show that dominantly (002) textured Zn substrates exhibit enhanced corrosion resistance and reduced interphase passivation. We introduce an effective Damköhler number (Da*) to concisely describe morphological evolution during electrocrystallization across substrates with different textures. High-texture (002) Zn substrates substantially enhance performance in high-capacity (~20 mA·hour/cm 2 ) symmetric Zn||Zn cells and full cells (Zn||δ-MnO 2 and Zn||I 2 ), enabling fast-charging and prolonged energy storage in coin and pouch rechargeable Zn battery formats.

Science & Technology - Other Topics↗
DOE OSTI1995

Nuclear Safety [Vol. 36, No. 2, July-December 1995]

Nuclear Safety is a journal that covers significant issues in the field of nuclear safety. Its primary scope is safety in the design, construction, operation, and decommissioning of nuclear power reactors worldwide and the research and analysis activities that promote this goal, but it also encompasses the safety aspects of the entire nuclear fuel cycle, including fuel fabrication, spent-fuel processing and handling, and nuclear waste disposal, the handling of fissionable materials and radioisotopes, and the environmental effects of all these activities. Table of Contents for this issue follows. THE CHORNOBYL ACCIDENT: 195 The Chornobyl Accident Revisited, Part III: Chernobyl Source Term Release Dynamics and Reconstruction of Events During the Active Phase, A. R. Sich; GENERAL SAFETY CONSIDERATIONS: 218 Second ANS Workshop on the Safety of Soviet-Designed Nuclear Power Plants, R. A. Bari; 234 Elements of a Nuclear Criticality Safety Program, C. M. Hopper; 243 Rickover, Excellence, and Criticality Safety Programs, R. E. Wilson; ACCIDENT ANALYSIS: 249 Transient Analysis of the PIUS Advanced Reactor Design with the TRAC-PF1/MOD2 Code, B. E. Boyack, J. L Steiner, S. C. Harmony, H. J. Stumpf, and J. F. Lime; 278 The Hierarchy-By-Interval Approach to Identifying Important Models that Need Improvement in Severe-Accident Simulation Codes, T. J. Heames, M. Khatib-Rahbar, J. E. Kelly, R. P. Jenks-Johnson, and Y.-S. Chen; 290 RELAP5/MOD3 Code Coupling Model, R. P. Martin; 299 Missiles Caused by Severe Pressurized-Water Reactor Accidents, R. Krieg; DESIGN FEATURES: 310 Validation of COMMIX with Westinghouse AP-600 PCCS Test Data, J. G. Sun, T. H. Chien, J. Ding, and W T. Sha; ENVIRONMENTAL EFFECTS: 321 Spent Nuclear Fuel Characterization for a Bounding Reference Assembly for the Receiving Basin for Off-Site Fuel, S. D. Kahook, R. L. Garrett, L. R. Canas, and M J. Beckum; OPERATING EXPERIENCES: 332 Reactor Shutdown Experience, Compiled by J. W. Cletcher; U.S. NUCLEAR REGULATORY COMMISSION INFORMATION AND ANALYSES: 335 Reactor Coolant System Blowdown at Wolf Creek on September 17, 1994, J. V. Kauffman and S. L. Israel; RECENT DEVELOPMENTS: 344 Reports, Standards, and Safety Guides, D. S. Oueener; 349 Proposed Rule Changes as of June 30, 1995; ANNOUNCEMENTS: 320 Symposium on Acceptability of Risk From Radiation—Application to Manned Space Flight; 320 24th DOE/NRC Nuclear Air Cleaning and Treatment Conference; 361 Radiation Biology and Radiation Protection— Modern Developments and Tendencies in Radiation Biology; 362 1997 IEEE Sixth Conference on Human Factors and Power Plants; 354 The Authors; 360 Reviewers of Nuclear Safety, Vol. 36.

05 NUCLEAR FUELS↗
DOE OSTI2023

The active CGCG 077-102 NED02 galaxy within the Abell 2063 galaxy cluster

Context.Within the framework of investigating the link between the central super massive black holes in the cores of galaxies and the galaxies themselves, we detected a variable X-ray source in the center of CGCG 077-102 NED02, which is a member of the CGCG 077-102 galaxy pair within the Abell 2063 cluster of galaxies. Aims.Our goal is to combine X-ray and optical data to demonstrate that this object harbors an active super massive black hole in its core, and to relate this to the dynamical status of the galaxy pair within the Abell 2063 cluster. Methods.We usedChandraandXMM-Newtonarchival data to derive the X-ray spectral shape and variability. We also obtained optical spectroscopy to detect the expected emission lines that are typically found in active galactic nuclei. Finally, we used public ZTF imaging data to investigate the optical variability. Results.There is no evidence of multiple X-ray sources or extended components within CGCG 077-102 NED02. Single X-ray spectral models fit the source well. We detect significant, nonrandom inter-observation 0.5–10 keV X-ray flux variabilities, for observations separated by ∼4 days for short-term variations and by up to ∼700 days for long-term variations. Optical spectroscopy points toward a passive galaxy for CGCG 077-102 NED01 and a Seyfert for CGCG 077-102 NED02. The classification of CGCG 077-102 NED02 is also consistent with its X-ray luminosity of over 10 42 erg s −1 . We do not detect short-term variability in the optical ZTF light curves. However, we find a significant long-term stochastic variability in theg-band that can be well described by the damped random walk model with a best-fit characteristic damping timescale ofτ DRW = 30 −12 +28 days. Finally, the CGCG 077-102 galaxy pair is deeply embedded within the Abell 2063 potential, with a long enough history within this massive structure to have been affected by the influence of this cluster for a long time. Conclusions.Our observations point toward a moderately massive black hole in the center of CGCG 077-102 NED02 of ∼10 6 M ⊙ . As compared to another similar pair in the literature, CGCG 077-102 NED02 is not heavily obscured, perhaps because of the surrounding intracluster medium ram-pressure stripping.

Astronomy & Astrophysics↗
DOE OSTI2025

Decoupling Failure Pathways in Li-O 2 Cells Operated Under Lean Electrolyte Conditions

The operation of lithium-oxygen (Li-O 2 ) batteries under lean electrolyte conditions offers higher energy density but leads to rapid degradation and short cycle life. Although cathode passivation and electrolyte decomposition occur in all regimes, we show that under lean electrolyte conditions, failure is primarily driven by progressive electrolyte consumption at the lithium/solid electrolyte interphase (SEI), rather than by irreversible cathode passivation. Poor wetting of the lithium surface results in heterogeneous SEI growth and high local current densities, which accelerate electrolyte loss and cell failure. Strategies aimed at improving interfacial stability, including optimized wetting and SEI forming additives, significantly extend cycle life without compromising energy density. Our results establish anode-electrolyte interactions as the dominant degradation mechanism under lean electrolyte conditions, and emphasize the need to engineer a stable Li/SEI interface for long-lasting Li-O 2 batteries.

Córdoba, Daniel [Argonne National Laboratory (ANL)↗
DOE OSTI2026

Thermal Reservoir Networks for Modularly Expandable Thermal Microgrids

The Department of Defense (DoD) faces the substantial challenge of cost-effectively retrofitting one to two installations per month, each comprising approximately 1,000 buildings, to improve resilience, reduce energy consumption, and enhance energy supply security. Achieving these objectives requires optimal system selection and effective risk mitigation during system integration. To address this need, we introduce Platform-Based Design (PBD), a structured, hierarchical methodology adapted from other industrial sectors to the domain of energy system retrofits. We demonstrate the effectiveness of PBD through a techno-economic feasibility study comparing geothermal-coupled thermal energy networks (TENs) with conventional energy systems for heating, cooling, and powering 17 buildings at Joint Base Andrews (JBA) in Maryland. Our analysis illustrates that the PBD approach enables rigorous, data-driven, sequential decision making, resulting in a family of Pareto-optimal systems, among which the TEN emerged as the most promising solution. The selected TEN design integrates geothermal borefields, heat recovery heat pumps, photovoltaic (PV) arrays, and battery storage. Compared to the baseline system – gas heating combined with air-source chillers – the proposed TEN reduces annual imported energy by 74% and peak electricity demand by 45%, achieves a levelized cost of energy of $\$0.210$/kWh, and substantially enhances resilience. Life-cycle costs increase by approximately 6%, and initial investment costs are about 2.5 times higher than the baseline. However, if central plant infrastructure, district loops, and utility-scale PV and battery systems are privately funded and operated, the initial investment would fall below the baseline system cost. Critical to achieving these significant performance improvements were detailed nonlinear dynamic simulations coupling geothermal heat transfer, energy system operation, and realistic feedback control logic. These simulations identified essential design modifications and control strategy refinements that substantially reduced energy use, peak demand, and compressor shortcycling, thereby improving durability and reliability—issues that would have been significantly more expensive to resolve during operation. Additionally, the verification step highlighted sensitivities to key design parameters that could reduce initial investment by approximately $\$2$ million and reduce annual life-cycle costs more than $\$300,000$. We recommend adopting the PBD methodology for future feasibility studies and TEN pilot projects to gain valuable operational experience. Furthermore, we recommend that DoD invest in transferring and scaling the PBD methodology to other installations. This entails developing standardized computational frameworks and component libraries as well as training industry in conducting PBD. Such investments would enable rapid, robust, reliable, and cost-effective retrofits, supporting DoD’s ambitious energy system modernization goals.

Wetter, Michael [Lawrence Berkeley National Labora↗
DOE OSTI2025

Time evolution of a pumped molecular magnet—A time-resolved inelastic neutron scattering study

Introducing an experimental technique of time-resolved inelastic neutron scattering (TRINS), we explore the time-dependent effects of resonant pulsed microwaves on the molecular magnet Cr 8 F 8 Piv 16 . The octagonal rings of magnetic Cr 3+ atoms with antiferromagnetic interactions form a singlet ground state with a weakly split triplet of excitations at 0.8 meV. A 4.6 tesla field was applied to tune the splitting between two members of the triplet excited level |1$\rangle$ ↔ |2$\rangle$ to resonance with 105 GHz (0.434 meV) microwaves. The time-dependent occupations of the ground state |0$\rangle$, lower lying levels |1$\rangle$ and |2$\rangle$, and higher energy states |λ ≥ 3$\rangle$ were extracted during and after 20 s long microwave pulses incident along the (101) direction of a Cr 8 F 8 Piv 16 crystal held at 1.9 K. At significantly elevated spin temperatures, we found underpopulation relative to thermal equilibrium of |2$\rangle$ and spin-lattice thermalization time scales ranging from 1.6(2) s to 5.7(2) s depending on the power level. This contrasts with the relaxation time τ 1 (T → 0) = 27(5) 𝛍 s inferred for |2$\rangle$ from in situ Electron Spin Resonance measurements. By probing a broad range of excited states during intense microwave pumping, TRINS thus provides a first view of long lived excited states in a molecular antiferromagnet.

Science & Technology - Other Topics↗
DOE OSTI2026

Adaptation of virtual synchronous generators to dynamic conditions in power grids

Virtual synchronous generators (VSGs) are widely adopted as grid-forming controls for inverter-based resources. However, when grid conditions vary significantly as characterized by changes in short-circuit ratio (SCR) and the reactance-to-resistance (X/R) ratio, fixed-gain designs and the commonly used P–Q decoupling assumption can become inaccurate. Such conditions can degrade transient power performance, leading to oscillations, prolonged settling, and overshoot, particularly in stiff-grid operating points. This paper quantifies how grid strength and impedance-dependent coupling affect the active–reactive power dynamics of a conventional VSG over a broad range of SCR and X/R values. An adaptive VSG tuning framework is then developed by combining (i) a coupling-explicit, impedance-parameterized state-space model to enable systematic controller synthesis, (ii) a full-state-feedback law designed via pole placement to meet prescribed damping and settling-time specifications, and (iii) a physics-informed neural network (PINN)–based online grid-impedance estimator that updates controller gains in real time as grid conditions vary. Offline simulations in MATLAB/Simulink and real-time validation on an OPAL-RT platform show that the proposed method preserves consistent damping and settling behavior with reduced overshoot across wide SCR and X/R ranges, compared with fixed-gain VSG baselines.

Adaptive control↗
DOE OSTI2023

Hybrid Data‐Driven Discovery of High‐Performance Silver Selenide‐Based Thermoelectric Composites

Optimizing material compositions often enhances thermoelectric performances. However, the large selection of possible base elements and dopants results in a vast composition design space that is too large to systematically search using solely domain knowledge. To address this challenge, a hybrid data-driven strategy that integrates Bayesian optimization (BO) and Gaussian process regression (GPR) is proposed to optimize the composition of five elements (Ag, Se, S, Cu, and Te) in AgSe-based thermoelectric materials. Data is collected from the literature to provide prior knowledge for the initial GPR model, which is updated by actively collected experimental data during the iteration between BO and experiments. Within seven iterations, the optimized AgSe-based materials prepared using a simple high-throughput ink mixing and blade coating method deliver a high power factor of 2100 µW m −1 K −2 , which is a 75% improvement from the baseline composite (nominal composition of Ag 2 Se 1 ). In conclusion, the success of this study provides opportunities to generalize the demonstrated active machine learning technique to accelerate the development and optimization of a wide range of material systems with reduced experimental trials.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Quasicrystal stability and nucleation kinetics from density functional theory

The aperiodic order of quasicrystals bridges the amorphous and crystalline regime, so it has remained unclear whether quasicrystals are metastable or stable phases of matter. Density functional theory is often used to evaluate thermodynamic stability, but quasicrystals are long-range aperiodic and their energies cannot be calculated using conventional ab initio methods. Here, in this work, we perform first-principles calculations on quasicrystal nanoparticles of increasing size, from which we can directly extrapolate their bulk and surface energies. Using this technique, we determine with high confidence that the icosahedral quasicrystals ScZn 7.33 and YbCd 5.7 are ground-state phases, thus revealing that translational symmetry is not a necessary condition for the zero-temperature stability of inorganic solids. Although we found the ScZn 7.33 quasicrystal to be thermodynamically stable, we show on a mixed thermodynamic and kinetic phase diagram that its solidification from the melt is limited by nucleation, which illustrates why even stable materials may be kinetically challenging to grow. Our techniques broadly open the door to first-principles investigations into the structure–bonding–stability relationships of aperiodic materials.

density functional theory↗