Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “DFT Calculations”

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 OSTI2021

Two-bond 13 C– 13 C spin-coupling constants in saccharides: dependencies on exocyclic hydroxyl group conformation

Seven doubly 13 C-labeled isotopomers of methyl β-D-glucopyranoside, methyl β-D-xylopyranoside, methyl β-D-galactopyranoside, methyl β-D-galactopyranosyl-(1→4)-β-D-glucopyranoside and methyl β-D-galactopyranosyl-(1→4)-β-D-xylopyranoside were prepared, crystallized, and studied by single-crystal X-ray crystallography and solid-state 13 C NMR spectroscopy to determine experimentally the dependence of 2 J C1,C3 values in aldopyranosyl rings on the C1–C2–O2–H torsion angle, θ 2 , involving the C2 carbon of the C1–C2–C3 coupling pathway. Using X-ray crystal structures to determine θ 2 in crystalline samples and by selecting compounds that exhibit a relatively wide range of θ 2 values in the crystalline state, 2 J C1,C3 values measured in crystalline samples were plotted against θ 2 and the resulting plot compared to that obtained from density functional theory (DFT) calculations. For θ 2 values ranging from ~90° to ~240°, very good agreement was observed between the experimental and theoretical plots, providing strong validation of DFT-calculated spin-coupling dependencies on exocyclic C–O bond conformation involving the central carbon of geminal C–C–C coupling pathways. we report that these findings provide new experimental evidence supporting the use of 2 J CCC values as non-conventional spin-coupling constraints in MA'AT conformational modeling of saccharides in solution, and the use of NMR spin-couplings not involving coupled hydroxyl hydrogens as indirect probes of C–O bond conformation. Solvomorphism was observed in crystalline βGal-(1→4)-βGlcOCH 3 wherein the previously-reported methanol solvate form was found to spontaneously convert to a monohydrate upon air-drying, leading to small but discernible conformational changes in, and a new crystalline form of, this disaccharide.

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

Dy adsorption on and intercalation under graphene on 6 H -SiC(0001) surface from first-principles calculations

Previous experimental observations motivate clarification of configuration stabilities and kinetic processes for intercalation of guest atoms into a layered van der Waals material such as a graphene-SiC system. From our first-principles density functional theory (DFT) calculations, we analyze Dy adsorption and intercalation for graphene on a 6H-SiC(0001) surface, where the system includes two single-atom-thick graphene layers: the top-layer graphene (TLG) and the underling buffer-layer graphene (BLG) above the terminal Si layer. Our chemical potential analysis shows that intercalation of a single Dy atom into the gallery between TLG and BLG is more favorable than adsorption on TLG but that intercalation into the gallery underneath BLG is highly unfavorable. We obtain diffusion barriers of ~0.45 and 0.54 eV for a Dy atom diffusing on and under TLG, respectively. We find that the direct penetration of a Dy atom from the graphene top into the gallery under TLG is almost inhibited below a temperature of ~1400 K due to a large global barrier of at least ~3.5 eV. Instead, we find that a single Dy atom on TLG can easily intercalate by crossing a TLG step (e.g., a zigzag step presaturated by a Dy chain or a reconstructed zigzag step zz57). Additionally, we also perform DFT calculations for different Dy coverages to demonstrate how the favorability of Dy intercalation, as well as the corresponding interlayer spacings, depend on the coverage. Consequently, we can provide general insight and guidance for extensively studied systems involving intercalation of foreign atoms into graphene on a SiC substrate.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Hidden Local Symmetry Breaking in a Kagome-Lattice Magnetic Weyl Semimetal

Exploring the relationship between intriguing physical properties and structural complexity is a central topic in studying modern functional materials. Co 3 Sn 2 S 2 , a newly discovered kagome-lattice magnetic Weyl semimetal, has triggered intense interest owing to the intimate coupling between topological semimetallic states and peculiar magnetic properties. However, the origins of the magnetic phase separation and spin glass state below T C in this ordered compound are two unresolved yet important puzzles in understanding its magnetism. Here, we report the discovery of local symmetry breaking surprisingly co-emerges with the onset of ferromagnetic order in Co 3 Sn 2 S 2 , by a combined use of neutron total scattering and half-polarized neutron diffraction. An anisotropic distortion of the cobalt kagome lattice at the atomic/nano level is also found, with distinct distortion directions among the two Co1 and four Co2 atoms. The mismatch of local and average symmetries occurs below T C , indicating that Co 3 Sn 2 S 2 evolves to an intrinsically lattice disordered system when the ferromagnetic order is established. The local symmetry breaking with intrinsic lattice disorder provides new understanding of the puzzling magnetic properties. Our density functional theory (DFT) calculation indicates that the local symmetry breaking is expected to reorient local ferromagnetic moments, unveiling the existence of the ferromagnetic instability associated with the lattice instability. Furthermore, DFT calculation unveils that the local symmetry breaking could affect the Weyl property by breaking the mirror plane. In conclusion, our findings highlight the fundamentally important role that the local symmetry breaking plays in advancing our understanding on the magnetic and topological properties in Co 3 Sn 2 S 2 , which may draw attention to explore the overlooked local symmetry breaking in Co 3 Sn 2 S 2 , its derivatives and more broadly in other topological Dirac/Weyl semimetals and kagome-lattice magnets.

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

Structural heterogeneity in non-crystalline Te x Se1−x thin films

Rapid crystallization behavior of amorphous TexSe1−x thin films limits the use of these alloys as coatings and in optoelectronic devices. Understanding the short- and medium-range ordering of the amorphous structure and the fundamental physics governing the crystallization of the films is crucial. Although the lack of long range crystalline order restricts the characterization of the amorphous films, electron microscopy offers a way to extract information about the nanoscale ordering. In this paper, the local ordering of amorphous TexSe1−x thin films with x=0.22, 0.61, 0.70, 0.90, and 1 grown by thermal evaporation is investigated using radial distribution function (RDF) and fluctuation electron microscopy (FEM) analysis. RDF results show that the nearest-neighbor distances of selenium (Se) and tellurium (Te) in their crystalline structure are preserved, and their bond lengths increase with the addition of Te. Density functional theory (DFT) calculations predict structures with interatomic distances similar to those measured experimentally. Additionally, fluctuations in atomic coordination are analyzed. Medium range order (MRO) analysis obtained from FEM and DFT calculations suggests that there are at least two populations within the chain network structure, which are close to the Se–Se and Te–Te intrachain distances. For the binary alloy with x > 0.61, TexSe1−x, Te–Te like populations increase and Te fragments might form, suggesting that the glass forming ability decreases rapidly.

Sari, Bengisu (ORCID:0000000244217098)↗
DOE OSTI2023

Machine Learning for First Principles Calculations of Material Properties for Ferromagnetic Materials

The investigation of finite temperature properties using Monte-Carlo (MC) methods requires a large number of evaluations of the system’s Hamiltonian to sample the phase space needed to obtain physical observables as function of temperature. DFT calculations can provide accurate evaluations of the energies, but they are too computationally expensive for routine simulations. To circumvent this problem, machine-learning (ML) based surrogate models have been developed and implemented on high-performance computing (HPC) architectures. In this paper, we describe two ML methods (linear mixing model and HydraGNN) as surrogates for first principles density functional theory (DFT) calculations with classical MC simulations. These two surrogate models are used to learn the dependence of target physical properties from complex compositions and interactions of their constituents. We present the predictive performance of these two surrogate models with respect to their complexity while avoiding the danger of overfitting the model. An important aspect of our approach is the periodic retraining with newly generated first principles data based on the progressive exploration of the system’s phase space by the MC simulation. The numerical results show that HydraGNN model attains superior predictive performance compared to the linear mixing model for magnetic alloy materials.

Eisenbach, Markus↗
DOE OSTI2025

Hybrid Doping Strategy with High‐Entropy Cu/Fe Surface Modification and Zr Bulk Incorporation for Ni‐Rich Cathodes

A hybrid doping strategy combining Zr 4+ bulk doping with high-entropy Cu 2+ /Fe 3+ surface doping is developed to enhance the structural and interfacial stability of Ni-rich layered oxide cathodes. Cu and Fe are selectively introduced at the particle surface via a surface-selective ion-exchange process, forming a ≈15 nm Fe-rich layer while preserving the layered framework. Compared to the pristine cathode, the hybrid sample exhibits significantly improved electrochemical performance in both half-cell and full-cell configurations. In half-cells, the hybrid retains 88.5% and 90.2% after 100 cycles at 1C under 4.6 and 4.5 V, respectively. During high-voltage full-cell cycling, the hybrid cathode maintains over 80% capacity retention, whereas the pristine counterpart retains less than 10% under identical conditions over the same cycling period. XPS, EELS, and DEMS analyses confirm improved oxygen retention, suppressed gas evolution, and stable surface chemistry, while DFT calculations indicate enhanced Me–O bonding in the selected Fe 0.75 Cu 0.25 (Mn 1/16 Co 2/16 Ni 13/16 )O 2 surface composition, which is identified through DFT-calculated mixing energy reaching a minimum at this ratio, indicating the most thermodynamically favorable configuration. In conclusion, these results demonstrate the effectiveness of this hybrid doping strategy in mitigating coupled degradation pathways in Ni-rich cathodes.

15 GEOTHERMAL ENERGY↗
DOE OSTI2023

Minimum conditions for accurate modeling of urea production via co-electrolysis

Co-electrolysis of carbon oxides and nitrogen oxides promise to simultaneously help restore the balance of the C and N cycles while producing valuable chemicals such as urea. However, co-electrolysis processes are still largely inefficient and numerous knowledge voids persist. Here, we provide a solid thermodynamic basis for modelling urea production via co-electrolysis. First, we determine the energetics of aqueous urea produced under electrochemical conditions based on experimental data, which enables an accurate assessment of equilibrium potentials and overpotentials. Next, we use density functional theory (DFT) calculations to model various co-electrolysis reactions producing urea. The calculated reaction free energies deviate significantly from experimental values for well-known GGA, meta-GGA and hybrid functionals. These deviations stem from errors in the DFT-calculated energies of molecular reactants and products. In particular, the error for urea is approximately -0.25 ± 0.10 eV. Finally, we show that all these errors introduce large inconsistencies in the calculated free-energy diagrams of urea production via co-electrolysis, such that gas-phase corrections are strongly advised.

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

Path Integral Monte Carlo Simulations of Iron Plasmas (Final Technical Report)

This documents is the final technical report for our grant entitled "Path Integral Monte Carlo Simulations of Iron Plasmas" that focused on developing path integral Monte Carlo (PIMC) computer simulations. This techniques will be developed to study plasmas composed of heavier elements including iron and other third row elements. Equations of state (EOS) and transport properties will be derived in the regime of warm dense matter (WDM) and dense plasmas where existing first-principles methods cannot be applied. While standard density functional theory (DFT) has been used to accurately predict the structure of many solids and liquids up to temperatures on the order of 100,000 K, this method is not applicable at much higher temperature because the number of partially occupied electronic orbitals reaches intractably large numbers or the use of finite-temperature free energy functionals in orbital-free DFT introduces an uncontrolled approximation. Here we focus on PIMC methods that become more and more efficient with increasing temperatures and still include all electronic correlation effects. In this approach, electronic excitations increase the efficiency rather than reduce it. While it had commonly been assumed this method could only be applied to elements without core electrons, we showed that PIMC with free-particle nodes works well for first-row elements (PRL 108 (2012) 115502). Most recently, we extended the applicability range of all-electron PIMC to second-row elements by adopting localized nodal surfaces (PRL 115 (2015) 176403). To simulate third-row elements efficiently under WDM conditions, we propose a new method to remove core electrons by introducing pseudo-nodes. We explain our approach step by step and present preliminary results. We focus our method development on getting PIMC simulations of iron to work because of its fundamental importance for WDM and astrophysics. Then we move on to krypton and copper-doped beryllium, a ICF ablator material. We plan to continue working on key second-row material such as Na, Mg, MgO, Al, silica, and silicon-doped plastic ablators. Our collaborators at LLNL, will use our PIMC EOS data both as comparisons to existing semi-empirical, EOS-generating schemes, and as input for continuum radiation hydrodynamics simulations. We will establish an efficient pipeline from PIMC to macroscopic continuum studies of materials response. An emphasis will be placed on benchmarking such methods for plasmas of heavy elements at the very high temperatures (~100 eV) and low densities that are generated when Hohlraum radiation heats the ablator material in indirect drive laser experiments. Results from changes to the EOS will be of immeasurable importance to the designers at the National Ignition Facility (NIF) and at other facilities. Starting with our EOS of Cu-doped Be, our second collaborator at LLE, will perform real-time simulations of laser fusion experiments at the Omega laser and at the NIF to determine how sensitive the compression path depends on the ablator EOS. Since our collaborator also has experience in performing orbital-free DFT calculations, we propose to compare predictions from this method with PIMC results. In joint publications, we plan to analyze the accuracy of different free-energy functionals in order to understand why existing orbital-free DFT calculations do not predict compression peaks along the shock Hugoniot curve that we see with PIMC. The peaks are caused by the ionization of various electron shells. Their accurate characterization is important to compare with experimental results. We will break new ground by developing PIMC techniques that can simulate iron and all other third row elements in the plasma and WDM regimes. We introduce the concept of pseudo-nodes for the efficient treatment of the core-electrons. The EOS and transport properties will be derived and published online in the form of a new WDM database. Our PIMC EOS calculations will benchmark and possibly replace semi-analytical EOS tables like QEOS or SESAME, which will impact the hydrocode simulation community and will affect the design of NIF targets. Our PIMC results will help to improve the accuracy of orbital-free DFT simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Doping Effects on the Ductility of a Lightweight Refractory High-Entropy Alloy: Grain Boundary and Bulk Lattice Aspects

Doping elements in small amounts often segregate to grain boundaries (GBs) in alloys and can significantly impact mechanical properties and performance. Refractory high-entropy alloys (RHEAs) are known for their poor ductility, especially at low temperatures. Promoting GB cohesion through segregation can be an effective approach to mitigate embrittlement. Here, in this study, first-principles density functional theory (DFT) calculations were performed to examine the effects of important interstitial dopants (O, B, C, and N) and substitutional dopants (Cr, Y, La, and Ce) on the Σ5(310) [001] tilt GB of a lightweight RHEA Nb 32.5 Ti 27.5 Mo 22.5 Ta 12.5 Hf 2.5 Zr 2.5 . The DFT calculations reveal that certain dopants, such as B, C, Cr, N, and Y, exhibit favorable GB strengthening effects by improving bonding interactions with the bulk alloy. The impact of doping on the ductility parameter of the bulk lattice, defined as the ratio of surface energy to unstable stacking fault energy for the {110} <111> slip system, was also studied; and the results show that doping reduces the intrinsic ductility of the alloy, decreasing the D-parameter from 2.90 to 2.55, depending on the specific dopant. The present findings provide a foundational understanding at atomic level of the effect of representative dopants on mechanical properties of RHEAs and can be used to guide future alloy design for improved mechanical properties.

deformation charge density↗
DOE OSTI2022

Atomic-Level Structure of Mesoporous Hexagonal Boron Nitride Determined by High-Resolution Solid-State Multinuclear Magnetic Resonance Spectroscopy and Density Functional Theory Calculations

Mesoporous hexagonal boron nitride (p-BN) has received significant attention over the last decade as a promising candidate for water cleaning/pollutant removal and hydrogen storage applications. In this work, high-resolution solid-state NMR spectroscopy and plane-wave density-functional theory (DFT) calculations are used to obtain an atomic-level description of p-BN. 1 H– 15 N or 1 H– 14 N heteronuclear (HETCOR) correlation experiments recorded with either conventional NMR at room temperature or dynamic nuclear polarization surface-enhanced spectroscopy (DNP-SENS) at ca. 100 K reveal NB 2 H, NBH 2 , NBH 3 + species residing on the edges of BN sheets. Ultra-high field 35.2 T 11 B NMR spectroscopy was used to resolve 11 B NMR signals from BN 3 , BN 2 O x (OH) 1–x (x = 0–1), BNO x (OH) 2–x (x = 0–2), BO x (OH) 3–x (x = 0–3), and BO x (OH) 4–x – (x = 0–4). Importantly, 2D 11 B dipolar double-quantum–single-quantum homonuclear correlation spectra reveal that many pore/defect sites are composed of boron oxide/hydroxide clusters connected to the BN framework through BN 2 O units. 1D and 2D 11 B{ 15 N} HETCOR NMR experiments, in addition to plane-wave DFT calculations of nine different structural models, further confirm the assignment of all NMR signals. The detailed structure determination of the pore and edge/defect sites within p-BN should further enable the rational design and development of next-generation p-BN-based materials. In addition, the techniques outlined here should be applicable to determine structure within other porous and/or boron-based materials.

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

Photochemistry of Hypervalent Iodoazide Derivatives

The photochemical properties and reaction mechanisms of a series of hypervalent iodoazide compounds (R–IN3) were investigated, with substituents (−CH 3 , −H, −CF 3 ) tuning the electronic density of the phenyl ring. With the help of UV irradiation, ultrafast time-resolved spectroscopy, and density functional theory (DFT) calculations, we elucidated the mechanisms of azide radical (N 3 • ) release and its subsequent reactivity. UV–vis spectroscopy reveals that the photoconversion rates follow the trend CH3 > H > CF3, aligning with DFT-calculated ΔG values for ring-opening transitions. Homolytic cleavage of the I–N bond is identified as the dominant pathway for N3• generation upon UV irradiation, occurring within 400 fs. The released azide radicals react with the solvent molecules, while the iodo radical R–I C • fragments undergo a thermodynamically uphill lactone ring-opening (RO) reaction and subsequent hydrogen atom abstraction from the solvent to form carboxylic acids R–I–COOH, as validated by NMR and IR spectroscopy. The study also highlights the role of substituents in influencing reaction kinetics and intermediate stability, with electron-donating groups accelerating the release of the N 3 • species. This work bridges experimental observations with computational predictions, offering a foundation for future advancements in azide-based reactions and materials.

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

Understanding Trigger Linkage Dynamics in Energetic Materials Using Mixed Picramide Nitrate Ester Explosives

The ability to predict the handling sensitivity of new organic energetic materials has been a longstanding goal. We report the synthesis and characterization of six new nitropicramide energetic materials with mixed functional groups that mimic known explosives such as nitroglycerin, erythritol tetranitrate (ETN), and pentaerythritol tetranitrate (PETN). The molecules have been studied theoretically using quantum molecular dynamics (QMD) simulations and density functional theory (DFT) calculations to identify the weakest bond in the reactants - the trigger-linkages - which control handling sensitivity, and to quantify their specific enthalpies of explosion. In good accord with the drop weight impact sensitivity data, our calculations predict that the sensitivities of the molecules are very similar owing to the small variations of the energy output and rates of trigger linkage rupture. In addition, both the QMD and DFT calculations point to the nitropicramide N–NO 2 bonds as the trigger linkages rather than the more typical O–NO 2 bonds. We propose that the switch of the trigger linkage from the nitrate esters to the nitramine groups arises from the strongly electron withdrawing character of the adjacent trinitrobenzene groups.

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

Equation of state and strength of diamond in high-pressure ramp loading

Diamond is used extensively as a component in high-energy-density experiments, but existing equation of state (EOS) models do not capture its observed response to dynamic loading. In particular, in contrast with first-principles theoretical models of equilibrium multiphase EOS, no solid-solid phase changes have been detected, and no general-purpose EOS models match the measured ambient isotherm. We have performed density functional theory (DFT) calculations of the diamond phase to ~10 TPa, well beyond its predicted range of thermodynamic stability, and used these results as the basis of a Mie-Grüneisen EOS. We also performed DFT calculations of the elastic moduli, and calibrated an algebraic elasticity model for use in simulations. We then estimated the flow stress of diamond by comparison with the stress-density relation measured experimentally in ramp-loading experiments. Furthermore, the resulting constitutive model allows us to place a constraint on the Taylor-Quinney factor (the fraction of plastic work converted to heat) from the observation that diamond does not melt on ramp compression.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Revealing Robust Room Temperature Ferromagnetism in Gd‐Doped Few‐Layered MoS 2 Thin Films

2D MoS 2 holds great promise for spintronics, yet is limited by intrinsic diamagnetism. This study demonstrates inducing ferromagnetic behavior in MoS 2 films doped with 0.47% Gd, achieving an ultrahigh saturation magnetization of 454 emu/cm 3 in a few‐layered film over 11‐times higher than bulk films (40 nm). Raman spectroscopy, X‐ray photoelectron spectroscopy, X‐ray magnetic circular dichroism, and density functional theory (DFT) calculations reveal an interplay between Gd dopants and Mo, S vacancies (V 1Mo+2S ), leading to the formation of bound magnetic polarons (BMPs) that drive ferromagnetic ordering. H 2 S annealing and DFT calculations reveal that defect healing reduces the saturation magnetization by 83%. High sulfur migration barrier in few‐layered films helps preserve BMPs, thereby sustaining ferromagnetism, whereas lower migration barriers in bulk films lead to suppression. These findings highlight the synergy between Gd doping and defect engineering in achieving ultrahigh room‐temperature ferromagnetism, offering a scalable strategy for developing high‐performance 2D magnetic materials for spintronic applications.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Atomic Layer Deposition of Pt on the Surface Deactivated by Fluorocarbon Implantation: Investigation of the Growth Mechanism

Selective atomic layer deposition (ALD) using surface-controlled reactivity is attracting a great deal of attention as a simple bottom-up patterning process that can provide both continued dimensional scaling and accurate pattern placement for next-generation nanoelectronics. Previously, we reported topographically selective deposition through Pt ALD using a MeCpPtMe 3 precursor and an O 2 counter reactant on fluorocarbon (CF x )-modified surfaces; however, gradual loss of selectivity in the CF x -modified regions was observed during the Pt ALD process. This work develops a fundamental understanding of the microscopic growth mechanisms of Pt ALD on the CF x -modified surface using a combination of experimental analyses and theoretical methods. The Pt growth characteristics on the CF x surface are investigated within a temperature window from 225 to 350 °C, and the results show a sharp sensitivity to growth temperature, with significant Pt growth occurring at temperatures above 300 °C. Through density functional theory (DFT) calculations, the reaction energies for adsorption of oxygen and the MeCpPtMe 3 precursor as well as formation of reaction products of CF x degradation are determined. Based on experimental results in conjunction with the DFT calculations, we show that while lower temperature Pt ALD (<300 °C) on CF x -modified surfaces can be significantly retarded because of a lack of MeCpPtMe 3 chemisorption, the surface reaction pathways for Pt ALD on CF x -modified surfaces at temperatures above 300 °C can proceed through oxygen adsorption and CF 4 desorption followed by Pt nucleation in CF x -degraded regions.

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

KSSOLV 2.0: An efficient MATLAB toolbox for solving the Kohn-Sham equations with plane-wave basis set

KSSOLV (Kohn-Sham Solver) is a MATLAB toolbox for performing Kohn-Sham density functional theory (DFT) calculations with a plane-wave basis set. KSSOLV 2.0 preserves the design features of the original KSSOLV software to allow users and developers to easily set up a problem and perform ground-state calculations as well as to prototype and test new algorithms. Furthermore, it includes new functionalities such as new iterative diagonalization algorithms, k-point sampling for electron band structures, geometry optimization and advanced algorithms for performing DFT calculations with local, semi-local, and hybrid exchange-correlation functionals. It can be used to study the electronic structures of both molecules and solids. So, we describe these new capabilities in this work through a few use cases. We also demonstrate the numerical accuracy and computational efficiency of KSSOLV on a variety of examples.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Fermi surface and Berry phase analysis for Dirac nodal line semimetals: Cautionary tale from SrGa 2 and BaGa 2

A Berry phase of odd multiples of 𝜋 inferred from quantum oscillations (QOs) has often been treated as evidence for nontrivial reciprocal space topology. However, disentangling the Berry phase values from the Zeeman effect and the orbital magnetic moment is often challenging. In centrosymmetric compounds, the case is simpler as the orbital magnetic moment contribution is negligible. Although the Zeeman effect can be significant, it is usually overlooked in most studies of QOs in centrosymmetric compounds. Here, we present a detailed study on the nonmagnetic centrosymmetric SrGa 2 and BaGa 2 , which are predicted to be Dirac nodal line semimetals based on density functional theory (DFT) calculations. Evidence of the nontrivial topology is found in magnetotransport measurements. The Fermi surface topology and band structure are carefully studied through a combination of angle-dependent QOs, angle-resolved photoemission spectroscopy (ARPES), and DFT calculations, where the nodal line is observed in the vicinity of the Fermi level. Strong de Haas–van Alphen fundamental oscillations associated with higher harmonics are observed in both compounds, which are well fitted by the Lifshitz-Kosevich (LK) formula. However, even with the inclusion of higher harmonics in the fitting, we found that the Berry phases cannot be unambiguously determined when the Zeeman effect is included. We revisit the LK formula and analyze the phenomena and outcomes that were associated with the Zeeman effect in previous studies. Our experimental results confirm that SrGa 2 and BaGa 2 are Dirac nodal line semimetals. Additionally, we highlight the often overlooked role of spin-damping terms in Berry phase analysis.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Electronic signatures of successive itinerant, antiferromagnetic transitions in hexagonal La 2 Ni 7

We use high-resolution angle-resolved photoemission spectroscopy (ARPES) and density functional theory (DFT) calculations to study the electronic and magnetic properties of La 2 Ni 7 , an itinerant magnetic system with a series of three magnetic transition temperatures upon cooling, which end in a small moment antiferromagnetic (smAFM) ground state. Our APRES data reveal several electron and hole pockets that have hexagonal symmetry near the Γ point. We observe significant reconstruction of the band structure upon successive magnetic transitions at T 1 ~ 61 K, T 2 ~ 57 K and T 3 ~ 42 K. Several features observed in ARPES data were reasonably well reproduced by DFT calculations, while others were not. In particular, the flat band near E f predicted by DFT in AFM state, was seemingly absent in ARPES data. Our results detail the effects of magnetic ordering on the electronic structure in a Ni-based small moment antiferromagnet and highlight challenges of current computational methods.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗