Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ab initio phasing”

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

Ab Initio Simulations of Tritium Diffusion in Al 2 O and Intermetallic Al 12 (TM) 2.34 Aluminide Coating Phases

Density functional theory simulations have been carried out to investigate the diffusion of interstitial tritium in Al 12 (TM) 2.34 and Al 2 O bulk phases. In Al 12 (TM) 2.34 the transition metal (TM) sites are occupied on average by 58.93 at.% Fe, 18.52 at.% Cr, and 22.54 at.% Ni. While the insertion of interstitial tritium in Al 2 O lead to a strong disordering of the structure, we only investigated interstitial tritium in Al 12 (TM) 2.34 and its iron end-member (i.e., Al 12 Fe 2.34 ). Nine diffusion pathways have been investigated along the a-, b-, and c -axis of Al 12 (TM) 2.34 . The direction of fastest diffusion for interstitial tritium is found to be along the b-axis, with a calculated diffusion coefficient of ≈10 -10 m 2 .s -1 at 600 K, which is at least one order of magnitude faster than those previously calculated for interstitial tritium in other Al-rich phases such as Fe2Alx, Fe4Al13, and FeNiAl5 for which D T =10 -11 m 2 .s -1 , ≈10 -12 m 2 .s -1 , and D T ≈10 -13 m 2 .s -1 respectively. The comparison of the energy landscape between Al 12 (TM) 2.34 and Al 12 Fe 2.34 for the fastest diffusion pathways in each axis direction, found that some pathways are not sensitive to transition metal mixing, while other are more affected, leading to higher energy barrier in Al 12 Fe 2.34 in part due to a more energetically favorable formation of Fe—T bond compared to Ni—T. However, we found that both materials have similar diffusion coefficients as the fastest diffusion pathways occurs along pathways that are not very sensitive to transition metal mixing.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Diverse high-pressure chemistry in Y-NH 3 BH 3 and Y–paraffin oil systems

The yttrium-hydrogen system has gained attention because of near-ambient temperature superconductivity reports in yttrium hydrides at high pressures. We conducted a study using synchrotron single-crystal x-ray diffraction (SCXRD) at 87 to 171 GPa, resulting in the discovery of known (two YH3 phases) and five previously unknown yttrium hydrides. These were synthesized in diamond anvil cells by laser heating yttrium with hydrogen-rich precursors—ammonia borane or paraffin oil. The arrangements of yttrium atoms in the crystal structures of new phases were determined on the basis of SCXRD, and the hydrogen content estimations based on empirical relations and ab initio calculations revealed the following compounds: Y 3 H 11 , Y 2 H 9 , Y 4 H 23 , Y 13 H 75 , and Y 4 H 25 . The study also uncovered a carbide (YC 2 ) and two yttrium allotropes. Complex phase diversity, variable hydrogen content in yttrium hydrides, and their metallic nature, as revealed by ab initio calculations, underline the challenges in identifying superconducting phases and understanding electronic transitions in high-pressure synthesized materials.

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

Dynamic signatures of electronically nonadiabatic coupling in sodium hydride: a rigorous test for the symmetric quasi-classical model applied to realistic, ab initio electronic states in the adiabatic representation

Sodium hydride (NaH) in the gas phase presents a seemingly simple electronic structure making it a potentially tractable system for the detailed investigation of nonadiabatic molecular dynamics from both computational and experimental standpoints. The single vibrational degree of freedom, as well as the strong nonadiabatic coupling that arises from the excited electronic states taking on considerable ionic character, provides a realistic chemical system to test the accuracy of quasi-classical methods to model population dynamics where the results are directly comparable against quantum mechanical benchmarks. Here, using a simulated pump–probe type experiment, this work presents computational predictions of population transfer through the avoided crossings of NaH via symmetric quasi-classical Meyer–Miller (SQC/MM), Ehrenfest, and exact quantum dynamics on realistic, ab initio potential energy surfaces. The main driving force for population transfer arises from the ground vibrational level of the D 1 Σ + adiabatic state that is embedded in the manifold of near-dissociation C 1 Σ + vibrational states. When coupled through a sharply localized first-order derivative coupling most of the population transfers between t = 15 and t = 30 fs depending on the initially excited vibronic wavepacket. While quantum mechanical effects are expected due to the reduced mass of NaH, predictions of the population dynamics from both the SQC/MM and Ehrenfest models perform remarkably well against the quantum dynamics benchmark. Additionally, an analysis of the vibronic structure in the nonadiabatically coupled regime is presented using a variational eigensolver methodology.

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

An isolated water droplet in the aqueous solution of a supramolecular tetrahedral cage

Water under nanoconfinement at ambient conditions has exhibited low-dimensional ice formation and liquid–solid phase transitions, but with structural and dynamical signatures that map onto known regions of water’s phase diagram. Using terahertz (THz) absorption spectroscopy and ab initio molecular dynamics, we have investigated the ambient water confined in a supramolecular tetrahedral assembly, and determined that a dynamically distinct network of 9 ± 1 water molecules is present within the nanocavity of the host. The low-frequency absorption spectrum and theoretical analysis of the water in the Ga 4 L 6 12– host demonstrate that the structure and dynamics of the encapsulated droplet is distinct from any known phase of water. A further inference is that the release of the highly unusual encapsulated water droplet creates a strong thermodynamic driver for the high-affinity binding of guests in aqueous solution for the Ga 4 L 6 12– supramolecular construct.

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

Molecular dynamics simulation of metallic Al–Ce liquids using a neural network machine learning interatomic potential

Al-rich Al-Ce alloys have the possibility of replacing heavier steel and cast-irons for use in high-temperature applications. Knowledge about the structures and properties of Al-Ce alloys at liquid state is vital for optimizing the manufacture process to produce desired allows. However, reliable molecular dynamics simulation of Al-Ce alloy systems remains a great challenge due to the lack of accurate Al-Ce interatomic potential. In this work, an artificial neural network (ANN) deep machine learning (ML) method is used to develop a reliable interatomic potential for Al-Ce alloy. Ab initio molecular dynamics (AIMD) simulation data on Al-Ce liquid with small unit cell (~200 atoms) and on the known Al-Ce crystalline compounds are collected to train the interatomic potential using ANN-ML. The obtained ANN-ML model reproduces well the energies, forces, and atomic structure of Al 90 Ce 10 liquid and crystalline phases of Al-Ce compounds in comparison with ab initio results. The developed ANN-ML potential is applied in molecular dynamics simulations to study the structures and properties of metallic Al 90 Ce 10 liquid, which would provide useful insight for guiding experimental process to produce desired Al-Ce allows.

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

Accurate and efficient parameterization of an atomic cluster expansion (ACE) potential for ammonia under extreme conditions

We present a machine learning interatomic potential for ammonia designed to capture its complex multiphase behavior, including both molecular and superionic phases. The potential is based on the atomic cluster expansion (ACE) formulation and has been parameterized to facilitate high-fidelity molecular dynamics simulations of ammonia under extreme conditions, for pressures up to 100 GPa and for temperatures above 500 K and up to 6000 K. A diverse range of configurations was generated through high-quality ab initio molecular dynamics simulations, covering insulating and superionic ice phases, liquid ammonia, molecular nitrogen (N 2 ) and hydrogen (H 2 ), and metastable compounds that form upon dissociation, including $NH^{+}_{4}$, $H^{+}_{3}$, N 2 H 4 , and N 3 H. We demonstrate that the ammonia ACE potential accurately reproduces experimental and density functional theory predicted isotherms and Hugoniots. Crucially, the potential is able to capture the intricate phase behavior of ammonia, including the transition from insulating molecular fluid to the superionic phase. This work provides a robust interatomic potential that can be used for large-scale, accurate simulations of ammonia under extreme thermodynamic conditions, offering a powerful tool for investigating its behavior in various phases and applications.

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

Thermodynamic and Kinetic Roles of H2 in Structure Evolution of Urchin-like Co: A density functional theory study

Small gas molecules working as good capping agents play important roles in controlling the morphologies and surface structures of the metal nanocrystals. In present work, the thermodynamic and kinetic roles of H2 molecule played on the morphology of Co nanocrystal were systematically studied based on density functional theory (DFT). The Gibbs surface free energies (?) of Co(100), Co(110) and Co(111) under different hydrogen surface coverages were investigated by ab initio thermodynamics, based on which i) the phase diagrams of stable H coverage on each plane were obtained; ii) the morphology evolutions of the Co nanocrystals with various surface hydrogen coverages were further predicted by Wulff construction. It was revealed that addition of H2 could change the facet stability generating diverse morphological Co nucleus. The kinetic role of H2 was further explored by DFT during adatom Co surface diffusions under different H coverages (?H), which suggested that except for Co(100) at ?H of 0.56 ML, the surface H would hinder Co surface diffusions. The projected density of state (PDOS) gave a deeper insight that the electronic structures of the Co adatom could be alerted by addition of the surface H, which thereby affected its surface diffusion ability. X. Wang, N. Liu, Q. Zhang, X. Liang, and B. Chen were financially supported by the National Natural Science Foundation of China (No. 21476012, 21571012, and 91534201). D. Mei was supported by the US Department of Energy (DOE), Office of Science, Office of Basic Energy Sciences, Division of Chemical Sciences, Geosciences & Biosciences. Pacific Northwest National Laboratory (PNNL) is a multiprogram national laboratory operated for DOE by Battelle. Computing time was granted by the facility user proposal grand challenge of computational catalysis of the William R. Wiley Environmental Molecular Sciences Laboratory (EMSL). EMSL is a national scientific user facility located at Pacific Northwest National Laboratory (PNNL) and sponsored by DOE’s Office of Biological and Environmental Research.

Nanocrystal Co, Surface free energy, Surface diffu↗
DOE OSTI2024

Complete separation of benzene-cyclohexene-cyclohexane mixtures via temperature-dependent molecular sieving by a flexible chain-like coordination polymer

Abstract The separation and purification of C6 cyclic hydrocarbons (benzene, cyclohexene, cyclohexane) represent a critically important but energy intensive process. Developing adsorptive separation technique to replace thermally driven distillation processes holds great promise to significantly reduce energy consumption. Here we report a flexible one-dimensional coordination polymer as an efficient adsorbent to discriminate ternary C6 cyclic hydrocarbons via an ideal molecular sieving mechanism. The compound undergoes fully reversible structural transformation associated with removal/re-coordination of water molecules and between activated and hydrocarbon-loaded forms. It exhibits distinct temperature- and adsorbate-dependent adsorption behavior which facilitates the complete separation of benzene, cyclohexene and cyclohexane from their binary and ternary mixtures, with the record-high uptake ratios for C 6 H 6 /C 6 H 12 and C 6 H 10 /C 6 H 12 in vapor phase and highest binary and ternary selectivities in liquid phase. In situ infrared spectroscopic analysis and ab initio calculations provide insight into the host-guest interactions and their effect on the preferential adsorption and structural transformation.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Elucidating the role of reactive nitrogen intermediates in hetero-cyclization during hydrothermal liquefaction of food waste

Hydrothermal liquefaction (HTL) has emerged as a promising strategy for converting abundant, water-rich organic streams into an energy-dense, sustainable, biocrude. Despite major strides in improving oil yields and process intensification, a key underlying challenge persists in relating the composition of the feedstock to the ultimate fate of nitrogen, which is especially important for biocrude quality. To elucidate how nitrogen partitions into the respective gas, aqueous, oil and solid hydrothermal products, food waste was treated under HTL conditions (15 wt% solids, 575 K reaction temperature, and 1 hour reaction time) and elemental, chromatographic and mass spectroscopy analyses were performed on the products to establish dominant product classes and to close mass, carbon, and nitrogen balances. Liquid products (aqueous and biocrude) were found to contain nearly 60% of the nitrogen, with a majority of the biocrude-phase nitrogen in the form of amides. Ab initio thermodynamic and kinetic simulations on model reactants were used to evaluate potential reaction pathways involving reactive nitrogen intermediates. A subsequent kinetic model was evaluated at reaction temperatures ranging 300–1000 K and for a range of feedstock compositions to identify primary reaction pathways. The Maillard reaction is revealed to be the dominant pathway for converting reactive nitrogen intermediates into observed products, especially those in which nitrogen is present as a heterocyclic aromatic. Subsequent product analysis by Fourier Transform Ion Cyclotron Resonance Mass Spectrometry (FT-ICR MS) confirms strong agreement among the experimental data, kinetic model and proposed Maillard pathways.

09 BIOMASS FUELS↗
DOE OSTI2022

Experimentally informed structure optimization of amorphous TiO 2 films grown by atomic layer deposition

Amorphous titanium dioxide TiO 2 (a-TiO 2 ) has been widely studied, particularly as a protective coating layer on semiconductors to prevent corrosion and promote electron–hole conduction in photoelectrochemical reactions. The stability and longevity of a-TiO 2 is strongly affected by the thickness and structural heterogeneity, implying that understanding the structure properties of a-TiO 2 is crucial for improving the performance. This study characterized the structural and electronic properties of a-TiO 2 thin films (~17 nm) grown on Si by atomic layer deposition (ALD). Fluctuation spectra V(k) and angular correlation functions were determined with 4-dimensional scanning transmission electron microscopy (4D-STEM), which revealed the distinctive medium-range ordering in the a-TiO 2 film. A realistic atomic model of a-TiO 2 was established guided by the medium-range ordering and the previously reported short-range ordering of a-TiO 2 film, as well as the interatomic potential. The structure was optimized by the StructOpt code using a genetic algorithm that simultaneously minimizes energy and maximizes the match to experimental short- and medium-range ordering. The StructOpt a-TiO 2 model presents improved agreements with the medium-range ordering and the k-space location of the dominant 2-fold angular correlations compared with a traditional melt-quenched model. The electronic structure of the StructOpt a-TiO 2 model was studied by ab initio calculations and compared to the crystalline phases and experimental results. Finally, this work uncovered the medium-range ordering in a-TiO 2 thin films and provided a realistic a-TiO 2 structure model for further investigation of structure–property relationships and materials design. In addition, the improved multi-objective optimization package StructOpt was provided for structure determination of complex materials guided by experiments and simulations.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Shock compression of crystalline TeO 2 to the high-pressure fluid regime: Insights from ab initio molecular dynamics simulations

The shock response of fully-dense and porous crystalline tellurium dioxide (TeO 2 ⁠) to the high-pressure and high-temperature fluid regime was investigated within the framework of density functional theory with Mermin’s generalization to finite temperatures. The principal and porous shock Hugoniot curves were predicted from canonical ab initio molecular dynamics (AIMD) simulations, with the phase space sampled along isotherms up to 80 000 K, for densities ranging from ρ = 3 to 17 g/cm 3 . The polymorphs investigated are α-TeO 2 paratellurite (⁠P4 1 2 1 2), TeO 2 cotunnite (⁠Pnma⁠), and TeO 2 post-cotunnite (⁠P2 1 /m⁠). Based on the discontinuity found in the calculated U s – u p slope of TeO 2 post-cotunnite at a shock velocity of U s ≃ 8.35km/s and a particle velocity of u p ≃ 3.64 km/s, the shock melting temperature and pressure are predicted to be ≃ 6500 K and ≃ 170 GPa. Results from the AIMD simulations are in line with the static compression data of TeO 2 paratellurite and cotunnite, and with the recent shock Hugoniot data for single-crystal α- TeO 2 for pressures up to 85 GPa, obtained using the inclined-mirror method and the velocity interferometer system for any reflector combined with powder gun and two-stage light-gas gun.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

Ab-initio vibrational properties of WTe2 and Li0.5WTe2

We present ab-initio calculations of the Td and 1T’ phases of WTe2 and the Td’ phase of Li 0.5 WTe 2 using density functional theory (DFT). These include lattice relaxations, total energy calculations, phonon calculations, and calculations of elastic constants. We also provide a detailed view of the transformation from Td to Td’ using a non-primitive representation of the Td lattice that more closely resembles the Td’ structure. The details presented here are intended as a supplement to a forthcoming experimental-computational paper on lithiated WTe 2 .

36 MATERIALS SCIENCE↗
DOE OSTI2020

An ab initio molecular dynamics exploration of associates in Ba-Bi liquid with strong ordering trends

Fictive associates are widely used to describe and model liquid phases with strong ordering trends. However, little evidence is known about the assumed associates in most cases. In the present work, an ab initio molecular dynamics (AIMD) study is employed to investigate the characters of the Ba-Bi liquid, in which associates have been assumed in existing thermodynamic modeling. It is found that in the Ba rich melt, the Bi atoms are almost completely surrounded by Ba atoms. The Bi-centered coordination polyhedrons are strongly associated to crystalline structures of Ba 5 Bi 3 and Ba 4 Bi 3 with a longer lifetime than other polyhedrons during the AIMD simulations. In addition, these Bi-centered polyhedrons in Ba rich melt connect with each other through vertex, edge, face, and/or bipyramid sharing to form medium range orders (MRO). In the Bi rich melt, the Ba-centered polyhedrons also form MROs, but they are both structurally and compositionally diverse with a shorter lifetime. These findings from AIMD study provide evidences that there exist a strongly ordering Ba 4 Bi 3 associate and a weakly ordering BaBi 3 associate in the Ba-Bi liquid. Here, the predicted enthalpy of mixing in the liquid agrees well with the results by the CALPHAD modeling in the literature.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Modeling the Solvation and Acidity of Carboxylic Acids Using an Ab Initio Deep Neural Network Potential

Formic and acetic acid constitute the simplest of carboxylic acids, yet they exhibit fascinating chemistry in the condensed phase such as proton transfer and dimerization. The go-to method of choice for modeling these rare events have been accurate but expensive ab-initio molecular dynamics simulations. Here, we present a deep neural network potential trained using accurate ab-initio data that can be used in tandem with enhanced-sampling methods to perform an efficient exploration of the free-energy surface of aqueous solutions of weak carboxylic acids. In particular, we show that our model captures proton dissociation and provides a good estimate of the pK a , as well as the dimerization of formic and acetic acid. This provides a suitable starting point for applications in different research areas where computational efficiency coupled with the accuracy of ab-initio methods is required.

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

Investigating a novel magnetic MAX phase nitride and its (001)-surfaces

We report magnetic MAX phases including their surfaces exhibit promising functional properties for magneto-electronic devices and self-monitoring smart coatings. By employing an integrated ab-initio approach, here we investigate a new member of the MAX phase, Mn 2 AuN. This compound satisfies the chemical, mechanical and dynamical stability criteria, leading to a possibility of its synthesis, and exhibits electronic and elastic anisotropy. The identified ferromagnetic configuration (with high Curie temperature) is the most energetically favorable among the five different spin configurations i.e., non-magnetic, ferromagnetic, and three anti-ferromagnetic. From the coatings perspective, the surface properties of Mn 2 AuN(001) terminations are investigated considering the four possible surface termination models. Remarkably, the magnetism does not get vanished even after cleaving the bulk unit cell into the (001)-surfaces. While evaluating the surface energetics in different chemical potentials, the N-001 terminated surface comes out to be the most stable termination contrasting with carbide MAX phases and is the only termination that exhibits the magnetic characteristics.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Switching of Hybrid Improper Ferroelectricity in Oxide Double Perovskites

In ABO 3 -type perovskite oxides with Pnma symmetry, rotation (Q R+ , a 0 a 0 c + ) and tilt (Q T , a – a – c 0 ) of BO 6 octahedra are the two primary order parameters. These order parameters establish an inherent trilinear coupling with anti-ferroelectric A-site displacement (Q AFE ) to form the low-symmetry phase. The symmetry is further lowered in double perovskite oxides (DPOs) due to A/A' cation ordering. It in turn makes these systems polar via hybrid improper ferroelectric mechanism, primarily driven by Q R+ and Q T . Naturally, it has been believed that functionalities such as polarization can also be switched by tuning these primary order parameters. However, mystery around finding switching mechanism still remains. Our study based on density functional theory calculations combined with finite-temperature molecular dynamics simulations shows that the polarization switching is a two-step process, driven by out-of-phase rotation (Q R– , a 0 a 0 c – when Q T = 0 or, a – a – b – when Q T ≠ 0). A series of polar DPOs such as KLnFeOsO 6 [Ln = Sm, Gd, Dy, Tm (lanthanides) and Y (rare earth)], all belonging to P2 1 symmetry, are considered in this investigation. The polarization switching P ($\overrightarrow{P}$) occurs at a very high temperature of ~1150 K through a phase transition, from a polar (P2 1 ) phase with $\overrightarrow{P}$(+) to $\overrightarrow{P}$(-) via a non-polar P4/n phase. The switching itself is metastable in nature. The switching (both polarization and spin state) is only observed for a very short period of time (~23 ps) that poses limitation on using such a mechanism in memory device realization. We demonstrate a concurrent heating–cooling procedure to overcome such shortcoming. In conclusion, simulations conducted at 600 K further imply that long lasting switching can be achieved, at least for 1.2 ns for 600 K, and ideally for an infinite time, if the material is heated just above the T c followed by rapid cooling to a temperature below T c .

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

Static Subspace Approximation for Random Phase Approximation Correlation Energies: Applications to Materials for Catalysis and Electrochemistry

Modeling complex materials using high-fidelity, ab initio methods at low cost is a fundamental goal for quantum chemical software packages. The GW approximation and random phase approximation (RPA) provide a unified description of both electronic structure and total energies using the same physics in a many-body perturbative approach that can be more accurate than generalized-gradient density functional theory (DFT) methods. However, GW/RPA implementations have historically been limited to either specific materials classes or application toward small chemical systems. Here, the static subspace approximation allows for reduced cost full-frequency GW/RPA calculations and has previously been benchmarked thoroughly for GW calculations. Here, we describe our approach to including partial occupations of electronic orbitals in full-frequency GW and RPA calculations for the study of electrocatalysts. We benchmarked RPA total energy calculations using the subspace approximation across a diverse test suite of materials for a variety of computational parameters. The benchmarking quantifies the impact of different extrapolation procedures for representing the static polarizability at infinite screened cutoff, and shows that using screened cutoffs above 20-25 Ryd result in diminishing accuracy returns for predicting RPA total energies. Additionally, for moderately sized electrocatalytic models, 2-3 times fewer computational resources are used to compute RPA total energies by representing the static polarizability with 20-30% of the static subspace basis, with an error of approximately 0.01 eV or better in RPA adsorption energy calculations. Finally, we show that for these electrochemical models RPA can shift DFT adsorption energy shifts by up to 0.5 eV and that GW can frequently shift DFT eigenvalues of surface and adsorbate states by approximately 0.5-1 eV.

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

Topological electride phase of sodium at high pressures and temperatures

Sodium, a textbook example of a nearly-free-electron metal, exhibits unforeseen pressure-induced behavior including the stabilization of numerous polymorphs—some possessing extremely complex unit cells—as well as a metal-to-insulator transition to the iconic hP4 phase. However, until recently, most of the experimental and theoretical studies on solid sodium have been restricted to the low-temperature regime. In this work, using ab initio evolutionary structure searches coupled with quasiharmonic calculations, we discover seven new phases of sodium that are more stable than the known hP4 phase at pressure-temperature conditions that were recently attained in ramp-compression experiments. From these, our calculations suggest that P6 3 /m Na is the ground state between ≈250 GPa at 710 K and ≈350 GPa at 1900 K. Electronic structure calculations show that this phase is a topological semimetal with a Dirac nodal surface that is protected by a nonsymmorphic symmetry S 2z and an electride owing to its non-nuclear charge localized within one-dimensional honeycomb channels and zero-dimensional cages. Our results highlight the complexity of dense sodium's free-energy landscape and intricate electronic structure that emerges at finite temperatures and conditions where ionic cores overlap.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗