Search NASA⌕ Search

SEARCH · Search NASA

Results for “Discrete Element Method”

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.

700 records · Page 31
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Tunable Conversion of Poly(Acrylic Acid) to Vinyl Methyl Ether Copolymers via Electrochemical Decarboxylation

Recognizing the potential of polymer-to-polymer transformations for plastics valorization, a new electrochemical strategy for the decarboxylative functionalization of polyacrylates has been developed. In contrast to most electrochemical approaches to plastic waste valorization that use small-molecule mediators or only form deconstruction products, this report describes one of the first direct electrochemical redox upconversions of a commodity polymer and involves activation of poly(acrylic acid), a material commonly found in disposable diapers and other personal hygiene products. Via electrochemical decarboxylation, poly(acrylic acid) is converted to methyl ether-containing copolymers. Although methods to prepare the 1:1 regularly alternating poly(vinyl methyl ether-alt-maleic anhydride) are known, routes to generate copolymers with variable levels of vinyl methyl ether and acrylate moieties are unknown, but such polymers may have utility in similar applications, such as dentifrices. Tunable conversions with up to 62 mol% vinyl methyl ether content were achieved via controlled-current electrolyses at >100 mA. This approach enables efficient access to copolymers with broad range of comonomer content from a common waste polymer starting material.

Dauzvardis, Fabian [University of Delaware, Newark↗
DOE OSTI2025

Effect of Preparation Conditions of Fe@SiO2 Catalyst on Its Structure Using High-Pressure Activity Studies in a 3D-Printed SS Microreactor

Fischer–Tropsch synthesis (FTS) in a 3D-printed stainless steel (SS) microchannel microreactor was investigated using Fe@SiO2 catalysts. The catalysts were prepared by two different techniques: one pot (OP) and autoclave (AC). The mesoporous structure of the two catalysts, Fe@SiO2 (OP) and Fe@SiO2 (AC), ensured a large contact area between the reactants and the catalyst. They were characterized by N2 physisorption, H2 temperature-programmed reduction (H2-TPR), scanning electron microscopy (SEM), transmission electron microscopy (TEM), X-ray diffraction (XRD), Fourier-transform infrared spectroscopy (FTIR), X-ray photoelectron microscopy (XPS), and thermogravimetric analysis–differential scanning calorimetry (TGA-DSC) techniques. The AC catalyst had a clear core–shell structure and showed a much greater surface area than that prepared by the OP method. The activities of the catalysts in terms of FTS were studied in the 200–350 °C temperature range at 20-bar pressure with a H2/CO molar ratio of 2:1. The Fe@SiO2 (AC) catalyst showed higher selectivity and higher CO conversion to olefins than Fe@SiO2 (OP). Stability studies of both catalysts were carried out for 30 h at 320 °C at 20 bar with a feed gas molar ratio of 2:1. The Fe@SiO2 (AC) catalyst showed higher stability and yielded consistent CO conversion compared to the Fe@SiO2 (OP) catalyst.

Biochemistry & Molecular Biology↗
NASA NTRSDate not supplied

First-Principles Investigation of Phase Transitions in RuNb, RuTa, and Ru-Nb-Ta Ultrahigh Temperature Shape Memory Alloys

Ultrahigh temperature shape memory alloys (UHT-SMAs) have transition temperatures above 600 C, and they have found applications for sensing and actuating devices in aerospace industry. Among very few such UHT-SMAs currently known are Ru-based alloys such as RuNb and RuTa, whose martensite structures and phase transitions are totally different from those of NiTi-based SMAs and were poorly understood. In this work, we carried out a systematical study of RuNb using first-principles total energy calculations and molecular dynamics (MD) simulations. The transition paths and mechanisms in cubic → tetragonal →monoclinic transitions are revealed. and the transition sequence and martensitic transition temperatures are determined (MTTs) by evaluating the Gibbs free energies using thermodynamic integration. The calculated MTTs are in very good agreement with the experimental data. We found that the monoclinic phase at the second transition has the 𝑃21/m symmetry instead of experimentally identified 𝑃2/m. Our calculations demonstrate that RuTa has very similar phase transitions to those of RuTa. Furthermore, we studied the Ru0.5Nb0.25Ta0.25 ternary. Our results can explain the measured significant drop of MTT (~ 200 C) in the second transition for Ru0.5Nb0.25Ta0.25 compared with those of RuNb and RuTa, while in the first transition its MTT is between MTTs of RuNb and RuTa. The insights gained by this study and the verified ab initio methods for accurate MTT calculations can be applied to fast screen and quantitatively design novel UHT-SMAs having similar properties with desirable MTTs and much reduced cost.

Zhigang Wu↗
DOE OSTI2026

A Chemist’s Guide to Solid-State Qubits: Diamond Color Centers as Embedded Vacancy-Centered Molecules

Diamond point defects or color centers are at the crux of solid-state quantum technology applications, from computing to sensing. Modeling molecular qubits has advanced by leaps and bounds mainly due to the advanced quantum mechanical (QM) tools in the arsenal of quantum chemists for modeling (isolated) molecular systems, which should be easily translatable for analogous systems in the solid state. Here, we simulate a series of diamond defects: a carbon vacancy (V C ) and second row substitutions that are adjacent to a V C or D C V C , with D = B, N, and O, to evaluate their magnetic and optical properties. We applied molecular orbital principles, an embedded cluster method, and a multiconfigurational QM theory to computationally characterize the above-mentioned defects for a multitude of spin, charge, and excited states. Our computational approach delivered high-quality QM insights that compare well with experiments. Furthermore, with this framework, we confirm the spin-triplets B C V C – and O C V C 0 to be analogous to N C V C – , with calculated bright transitions that respectively correspond to UV (3.80 eV) and blue (2.65 eV) light absorption and emission lifetimes (τ) of 2.3 and 5.6 ns, which are higher in energy and more short-lived (and thus brighter) than that for N C V C – : green (2.26 eV), τ = 9.5 ns.

Martirez, John Mark P. [Princeton Plasma Physics L↗
DOE OSTI2026

Quantum Filtering and Analysis of Multiplicities in Eigenvalue Spectra

Fine-grained spectral properties of quantum Hamiltonians, including both eigenvalues and their multiplicities, provide useful information for characterizing many-body quantum systems as well as for understanding phenomena such as topological order. Extracting such information with small additive error is #BQP-complete in the worst case. In this work, we introduce QFAMES (quantum filtering and analysis of multiplicities in eigenvalue spectra), a quantum algorithm that efficiently identifies clusters of closely spaced dominant eigenvalues and determines their multiplicities under physically motivated assumptions, which allows us to bypass worst-case complexity barriers. QFAMES also enables the estimation of observable expectation values within targeted energy clusters, providing a powerful tool for studying quantum phase transitions and other physical properties. We validate the effectiveness of QFAMES through numerical demonstrations, including its applications to characterizing quantum phases in the transverse-field Ising model and estimating the ground-state degeneracy of a topologically ordered phase in the two-dimensional toric code model. We also generalize QFAMES to the setting of mixed initial states. Our approach offers rigorous theoretical guarantees and significant advantages over existing subspace-based quantum spectral analysis methods, particularly in terms of the sample complexity and the ability to resolve degeneracies.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Field-tailoring quantum materials via magneto-synthesis: metastable metallic and magnetically suppressed phases in a trimer iridate

We demonstrate that applying modest magnetic fields (< 0.1 T) during high-temperature crystal growth can profoundly alter the structure and ground state of a spin-orbit-coupled, antiferromagnetic trimer lattice. Using BaIrO₃ as a model system, whose ground state is intricately dictated by the trimer lattice, we show that magneto-synthesis , a field-assisted synthesis approach, stabilizes a structurally compressed, metastable metallic and magnetically suppressed phases inaccessible via conventional methods. These effects include a 0.85% reduction in unit cell, 4-order-of-magnitude decrease in resistivity, a 10-fold enhancement of the Sommerfeld coefficient, and the collapse of long-range magnetic order -- all intrinsic and bulk in origin. First-principles calculations confirm that the field-stabilized structure lies substantially above the ground state in energy, highlighting its metastable character. These large, coherent and correlated changes across multiple bulk properties, unlike those caused by dilute impurities, defects or off-stoichiometry, point to an intrinsic field-induced mechanism. The findings establish magneto-synthesis as a powerful new pathway for accessing non-equilibrium quantum phases in strongly correlated materials.

magneto-synthesis↗
DOE OSTI2026

Challenges and Opportunities for HVAC + Phase Change Material Thermal Energy Storage in Buildings

Phase-change material (PCM) thermal-energy storage (TES) integrated with HVAC and domestic hot water (DHW) can shift a large share of building thermal loads. By flattening and shifting loads, PCM TES reduces peak electricity use, eases stress on local and grid infrastructure, and lowers costs. It can also defer costly upgrades to service panels, distribution, and transmission. Higher energy density relative to chilled or hot water storage makes PCM TES practical for small, space limited, and retrofit projects, while packaged HVAC-integrated systems expand cost-effective load shifting to commercial buildings that previously lacked options. Despite this promise, deployment faces barriers. This paper presents challenges, opportunities, and lessons learned from lab and field integrations of PCM TES with packaged HVAC systems. Key challenges include misalignment between default heat pump controls tuned for direct-to-load operation and TES charge/discharge objectives, PCM properties that diverge from manufacturer claims and evolve with cycling, and high integration and deployment costs with off-the-shelf products. These studies underscore the need for factory built, integrated systems that arrive as prepackaged modules with proven controls, reducing design and installation effort and delivering predictable performance across projects. The paper outlines steps to enable viable and scalable PCM TES HVAC and DHW systems, including 1) validated methods that can be used to characterize PCM properties and TES system performance, 2) supervisory controls that optimize charge/discharge scheduling, and 3) factory integrated packaging that eliminates bespoke field engineering.

Dutton, Spencer↗
DOE OSTI2026

Robust Heat-Flux Sensors for Coal-Fired Boiler Extreme Environments

In this project, robust heat-flux measurement systems were developed. The heat-flux sensors utilize thermoelectric effects to directly transduce the heat-flux inputs to analog electrical voltage signals. They were constructed from dedicated materials that can withstand temperatures of at least 1000°C and maintain adequate performance at these conditions for prolonged periods of time. The proposed approaches took into account numerous considerations, including system cost, sensor head resilience, sensor footprint, data accuracy, response time, and maintenance requirements. Through modern thermoelectric materials design, methodical materials selection and rigorous testing in materials characterization labs and medium-scale fire research facilities, we have demonstrated functioning laboratory prototypes, upon which one could base industrial heat-flux sensing platforms capable of operating in the challenging high-temperature, corrosive environments of the boilers of coal-fired power plants. A distributed sensor array for heat-flux measurements throughout the furnace water-wall, the superheater area and the economizer coils can provide critical data for the power plant control systems to increase efficiency, improve safety and reduce down times. For example, the combined heat-flux sensor/control systems can contribute to the optimization of burner and boiler operations under flexible loads, the optimization of heat-exchange conditions and overall reduction of heat rate and emissions, the prediction of imminent overheating conditions, and the optimization of the soot-blowing protocols.

20 FOSSIL-FUELED POWER PLANTS↗
DOE OSTI2026

HTGR Multiphysics Application Drivers FY26 Updates

This report summarizes FY26 progress under the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program's high-temperature gas-cooled reactor (HTGR) application driver work, covering a wide range of activities such as code validation and multi-physics code assessment. 1) A detailed SAM model of the High-Temperature Engineering Test Reactor (HTTR) was developed using a unique-block grouping approach, with an extended parallel thermal network method to capture block-to-block conduction and radiation heat transfer, and applied to steady-state simulations of the HTTR 30~MW and 9~MW cases. 2) In another activity, SAM's newly implemented multi-component gas flow model was validated against the Natural convection Shutdown heat removal Test Facility (NSTF) argon ingress experiment, correctly capturing the density-driven suppression and thermal recovery of natural circulation observed when argon is introduced into the air-cooled Reactor Cavity Cooling System (RCCS) loop. 3) For the OECD/NEA High Temperature Test Facility (HTTF) benchmark, we co-led the international benchmark activities as well as the OECD/NEA final benchmark report to be released at the end of this year. 4) Finally, the coupled Griffin-SAM modeling capability for pebble-bed HTGRs was advanced by verifying the Griffin neutronics solution against Serpent Monte Carlo for a realistic non-uniform temperature distribution, resolving several deficiencies in the SAM-to-Griffin temperature transfer scheme, and enabling distinct fuel kernel, moderator, and coolant temperatures for cross section feedback. These new features were demonstrated in a PBR load-following transient.

Lee, Alvin↗
NASA NTRS2026

Spatially Refined Satellite Gravimetry Captures Human Signatures in Global Terrestrial Water Storage Trends

Human activities have directly altered the water cycle through water management, aquifer pumping, agricultural irrigation, and land use change. Although satellite gravimetry has transformed global hydrological research, its coarse resolution limits attribution of freshwater change to human activities at many management-relevant scales. Here we assessed global terrestrial water storage (TWS) trends from April 2002 to November 2025 using “stacked” regression of Level-1B intersatellite ranging data, which leverages temporal information and variability to dramatically improve effective spatial resolution relative to standard approaches. We combined this refined product with rigorous uncertainty analysis, autocorrelation-robust geostatistical methods, and literature assessment to evaluate TWS trend associations with land and water use, climate variability, and glacial mass loss. We identified TWS trend hotspots exhibiting significant spatial associations with anthropogenic 40 drivers, including groundwater and surface-water irrigation, rainfed agriculture, deforestation, and reservoir impoundment. Across these regions, cumulative TWS losses (3,122 Gt) substantially exceeded gains (2,432 Gt). Compared with traditional regression of monthly mascons, our approach yielded regional trend magnitudes that are on average 33% larger, revealing that global freshwater depletion, particularly from groundwater pumping, is considerably more acute than previously estimated. Multivariate regression models show that humans account for a significant share of the spatial variability in TWS trends on every non-polar continent except Australia. We detected localized TWS gains linked to rainfed agriculture, surface water irrigation, and reservoir filling that were unresolved in earlier gravimetric studies. The methodology provides a foundation for future gravity missions to independently track decadal freshwater change with unprecedented spatial fidelity.

groundwater↗
DOE OSTI2026

Investigation Into LiCl-LiF+Li 2 O as a Low Temperature Electrolyte in Direct Electrolytic Reduction of UO 2

Direct electrolytic reduction (DER) of UO 2 in a binary LiCl-Li2O (1–3 wt%) at 650 °C has been widely reported. The process temperature can be reduced to 550 °C in theory by using a near eutectic ternary LiCl-LiF+Li 2O (∼2 wt%) electrolyte, which will result in substantial reduction in rate of salt vaporization. The feasibility of using the ternary electrolyte for the process was tested via material compatibility testing, electrochemical polarization testing, and controlled potential or current DER experiments. MgO, Pt, and 625 nickel samples were each contacted with the ternary salt for 168 h. Negligible changes to dimensions and mass of the samples were measured, and concentrations of corrosion products in the salt were 102 ppm or less. O2−oxidation was observed to occur at a lower potential than Cl-or F-on a Pt anode. DER of UO2was performed in six trials. Process variables included type of reference electrode, electrochemical control method, type of cathode basket, and electrolyte composition. UO 2conversion ranged from 28.2 to 99.9 % , and current efficiency ranged from 38 to 52%. Overall, results indicate that the use of the ternary salt is feasible for DER.

36 MATERIALS SCIENCE↗
NASA NTRS2026

Hybrid Modeling Study on Grain Evolution in the Metal Welding Process and Its Potential Lunar Application

Metal is most commonly used structural material in a wide range of spacecraft, and welding is the principal method for joining metal components into functional systems. However, conducting welding experiments under extreme environments—such as microgravity or vacuum conditions in space—is prohibitively expensive and experimentally challenging. To overcome these limitations, multi-physics computational welding models provide a cost-effective and versatile alternative. In this work, the authors have developed a coupled thermal (fluid) microstructure simulation framework to model metal welding under varying gravity conditions. The framework integrates a mixed-mode heat transfer formulation (conduction, convection, and radiation) with molten pool fluid dynamics, enabling accurate prediction of temperature fields and weld-pool geometry. A grain growth model is further incorporated to capture the spatial and temporal evolution of microstructure, including grain size distribution and morphological transitions during solidification. This approach provides detailed insight into molten pool evolution and grain-level microstructure development throughout the welding process. By explicitly parameterizing environmental conditions, the model supports extrapolation to off-Earth manufacturing scenarios such as welding on the lunar surface. Tantalum—chosen in this study due to its high melting point, oxidation resistance, and mechanical stability at elevated temperatures—serves as the material system for model demonstration. Beyond Tantalum, the integrated multi-physics framework offers broad applicability for predictive welding simulations of various structural and refractory metals or alloys used in extreme terrestrial or extraterrestrial environments.

kinetic Monte Carlo (SPPARKS)↗
DOE OSTI2024

Efficient Leaching of Metal Ions from Spent Li-Ion Battery Combined Electrode Coatings Using Hydroxy Acid Mixtures and Regeneration of Lithium Nickel Manganese Cobalt Oxide

Extensive use of Li-ion batteries in electric vehicles, electronics, and other energy storage applications has resulted in a need to recycle valuable metals Li, Mn, Ni, and Co in these devices. In this work, an aqueous mixture of glycolic and lactic acid is shown as an excellent leaching agent to recover these critical metals from spent Li-ion laptop batteries combined with cathode and anode coatings without adding hydrogen peroxide or other reducing agents. An aqueous acid mixture of 0.15 M in glycolic and 0.35 M in lactic acid showed the highest leaching efficiencies of 100, 100, 100, and 89% for Li, Ni, Mn, and Co, respectively, in an experiment at 120 °C for 6 h. Subsequently, the chelate solution was evaporated to give a mixed metal-hydroxy acid chelate gel. Pyrolysis of the dried chelate gel at 800 °C for 15 h could be used to burn off hydroxy acids, regenerating lithium nickel manganese cobalt oxide, and the novel method presented to avoid the precipitation of metals as hydroxide or carbonates. The Li, Ni, Mn, and Co ratio of regenerated lithium nickel manganese cobalt oxide is comparable to this metal ratio in pyrolyzed electrode coating and showed similar powder X-ray diffractograms, suggesting the suitability of α-hydroxy carboxylic acid mixtures as leaching agents and ligands in regeneration of mixed metal oxide via pyrolysis of the dried chelate gel.

Electrochemistry↗
DOE OSTI2026

Large Language Model for Validation, Optical Calibration, and Learning (VOCAL) Distributed Temperature Sensing Interface

Distributed temperature sensing (DTS) using fiber optic sensors (FOS) offers a promising method for temperature measurements in advanced reactors, such as sodium fast reactors and molten salt cooled reactors. To support the calibration and validation of DTS measurements, Argonne National Laboratory developed the Validation, Optical Calibration, and Learning (VOCAL) software package. This report describes the integration of a local large language model (LLM) with a retrieval-augmented generation (RAG) system into the VOCAL interface to serve as an interactive user assistant. The LLM framework enhances the VOCAL platform’s accessibility to users by explaining interface components, clarifying inputs and outputs, and answering user queries dynamically in real-time. The accuracy of the LLM assistant performance was evaluated with 20 queries regarding the interface and its parameters using experimental data from the Thermal Hydraulic Experimental Test Article (THETA) facility. Results demonstrate that the LLM achieved a 95% accuracy rate, with a BERTScore of 0.8816 and SBERT value of 0.7417. Furthermore, validation of the RAG system within the LLM framework showed optimal accuracy with k-values between 1 and 2 using the k-refinement convergence test. The prompt perturbation analysis demonstrated good initial consistency for the RAG system, exhibiting the highest accuracy under punctuation variations and the greatest sensitivity under query reordering. Notably, the model’s errors were limited to data retrieval failures rather than factual hallucinations, reinforcing its baseline reliability. The integration of LLM provides a highly accurate, userfriendly enhancement to the VOCAL platform without disrupting its core computational capabilities for FOS calibration and validation.

Hong, Evan↗
DOE OSTI2026

Transition from Vehicular to Structural Ionic Transport in Electrified Alkali Aqueous Solutions

A molecular understanding of the solvation and dynamics of ions under static electric fields is crucial for modelling a wide range of natural and technological processes. Yet, traditional simulation methods suffer from a trade-off that has to be made between accuracy and statistical convergence. To bridge this gap, herein we extend our recently introduced Perturbed Neural Network Potential Molecular Dynamics (PNNP MD) approach to investigate the solvation structures and ionic transport mechanisms of electrified alkali cationic solutions. We obtain ionic conductivities for Li+, Na+ and Cs+ from the field dependence of the ionic current density in good agreement with experiment. Surprisingly, the migration mechanism is found to be strikingly different for the three ions despite their similar ionic conductivities. While Li+ conducts predominantly through vehicular migration of a stable 4-fold coordinated ion at all field strengths, Cs+ conducts strictly through a structural diffusion mechanism, where 9-12 transient first shell water coordination bonds are broken and renewed. Notably, aqueous Na+ emerges as a “Goldilocks” ion: its ion-water interactions are strong enough to maintain distinct 5-6 fold coordination shells at zero field (unlike Cs+) yet labile enough to be strongly perturbed by electric fields (unlike Li+). As a consequence, we observe an electric field-induced transition from vehicular to structural ionic transport for Na+ that is accompanied by a marked increase in the ionic current density. Our results imply that the conductance mechanism of ions with moderate ion-solvent interactions can be effectively tuned by external electric fields.

Joll, Kit↗
DOE OSTI2026

Integral Nuclear Data and Benchmarking Needs for Fusion Energy Systems

Fusion energy systems are currently being designed and optimized using radiation transport codes. To deal with the unique environment inside a fusion-based system, many of these designs incorporate novel materials able to withstand the high radiation fields, ensure adequate cooling and thermal protection, and produce tritium. Validation plays a vital role in building trust in the predictive power of these models and computational methods. Validation of a code consists of modeling documented real-world experiments and comparing the code-predicted response to the measured response. Adequate validation requires measured responses from real-world experiments, also known as integral data, that mimic the system being designed, including materials, impinging radiation, and temperature, among other variables. The most trusted integral data are experimental responses that have been through a rigorous benchmarking process that develops a recommended computational model and evaluates all experimental uncertainties. Finally, there are a few research groups around the world that have been producing integral data for fusion applications, but a substantial investment is needed to address the unique validation needs of the fusion community.

Fusion↗
DOE OSTI2026

Single-crystal NMC622 from combustion-assisted solid-state synthesis: Tunable particle size, voltage-dependent stability, and Li-ion diffusion

Cost-effective and scalable synthesis of single-crystal lithium–metal-oxide cathodes with controlled size and morphology remains challenging for advancing high-performance lithium-ion batteries. Here, we demonstrate a solution combustion–assisted solid-state method for producing single-crystal NMC622 with tunable particle size at various lithium stoichiometries. Structural and electronic characteristics were probed using synchrotron-XANES, EXAFS, XRD and XPS, confirming well-defined coordination environments and phase purity up to Li ≈ 1.3. The electrochemical behavior of LixNMC (x = 1.0 and 1.3) single crystals was evaluated across 2.8–4.3 V and 2.8–4.7 V windows to elucidate the effects of lithium content and operating voltage. At 2.8–4.3V, Li1.3NMC exhibits higher initial capacity due to its larger lithium inventory, while Li1.0NMC shows superior long-term retention driven by its larger crystal size and reduced structural distortion. Increasing the cutoff voltage to 4.7 V enhances the initial capacity of both compositions by 17–25%, but long-term capacity retention ultimately converges to values comparable to those at 4.3 V due to voltage-induced degradation, where Li+ diffusion constant remains lying ranging from 10−12 –10−11 cm2s−1. Overall, this study establishes a scalable combustion-assisted route for synthesizing tunable single-crystal NMC622 and highlights the interplay between lithium stoichiometry, particle size, and voltage window in governing cathode stability.

Roy, Subrata C [Jackson State University]↗
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↗