Search NASA⌕ Search

SEARCH · Search NASA

Results for “Time dependent density functional theory”

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

Defect-Limited Mobility and Defect Thermochemistry in Mixed A-Cation Tin Perovskites: (CH3NH3)1-xCsxSnBr3

Hybrid organic-inorganic semiconductors crystallizing in the perovskite structure present a significant opportunity for realizing defect-tolerant semiconductors. In this work, we examine the solid solution, (CH3NH3)1-xCsxSnBr3, and identify the thermochemistry dictating the intrinsic carrier concentrations and how local structural distortions influences this electronic behavior. This family of compounds exhibits the expected systematic trend in decreasing optical gap with the cesium to methylammonium A-site mixing ratio, x, in the visible region. However, the carrier mobility, as determined from time-resolved microwave conductivity measurements, trends opposite to that expected from first-principles calculations combined with Boltzmann scattering theory calculations of the carrier mobility. We propose that this is a result of increasing carrier scattering with x in (CH3NH3)1-xCsxSnBr3 due to a significant increase in the carrier density with x. By examining the dependence of the carrier density as a function of x, we infer the compositional-dependence of the average enthalpy and nonconfigurational entropy per defect. While diffraction reveals a cubic aristotypic perovskite structure as a function of x at room temperature, the pair distribution functions obtained from synchrotron X-ray total scattering from these materials are better described by symmetry-adapted displacement modes of the Pm3m crystal structure, which we attribute to a large degree of anharmonic dynamics of the atom positions. This analysis shows that CH3NH3+-rich compositions retain mostly linear Sn-Br-Sn bonding environments through the displacements, while Cs+-rich compositions lead to more significantly bent Sn-Br-Sn environments. We propose that these bent bonding arrangements in Cs-rich compositions yield a higher propensity for defect formation. This also provides a rationale for carrier trapping that gives rise to anomalous microwave transients. Together, these results provide insight into the structure-dynamics-properties relationships in this highly anharmonic system with high amplitude atomic motions and low defect formation energies.

defect thermochemistry↗

Mixed Valence {Ni 2+ Ni 1+ } Clusters as Models of Acetyl Coenzyme A Synthase Intermediates

Acetyl coenzyme A synthase (ACS) catalyzes the formation and deconstruction of the key biological metabolite, acetyl coenzyme A (acetyl-CoA). The active site of ACS features a {NiNi} cluster bridged to a [Fe4S4] n+ cubane known as the A-cluster. The mechanism by which the A-cluster functions is debated, with few model complexes able to replicate the oxidation states, coordination features, or reactivity proposed in the catalytic cycle. In this work, we isolate the first bimetallic models of two hypothesized intermediates on the paramagnetic pathway of the ACS function. The heteroligated {Ni 2+ Ni 1+ } cluster, [K(12-crown-4) 2 ][1], effectively replicates the coordination number and oxidation state of the proposed “A red ” state of the A-cluster. Addition of carbon monoxide to [1] - allows for isolation of a dinuclear {Ni 2+ Ni 1+ (CO)} complex, [K(12-crown-2) n ][2] (n = 1–2), which bears similarity to the “A NiFeC ” enzyme intermediate. Structural and electronic properties of each cluster are elucidated by X-ray diffraction, nuclear magnetic resonance, cyclic voltammetry, and UV/vis and electron paramagnetic resonance spectroscopies, which are supplemented by density functional theory (DFT) calculations. Calculations indicate that the pseudo-T-shaped geometry of the three-coordinate nickel in [1] – is more stable than the Y-conformation by 22 kcal mol –1 , and that binding of CO to Ni 1+ is barrierless and exergonic by 6 kcal mol –1 . UV/vis absorption spectroscopy on [2] - in conjunction with time-dependent DFT calculations indicates that the square-planar nickel site is involved in electron transfer to the CO π*-orbital. Further, we demonstrate that [2] - promotes thioester synthesis in a reaction analogous to the production of acetyl coenzyme A by ACS.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Vacancy-Dependent Diffusion Mechanism in Oxygen-Defective SrFeO 3 Perovskite Materials: First-Principles Density Functional Theory and Experimental Approach

Understanding oxygen diffusion at the atomic scale in SrFeO 3−δ perovskites is crucial for developing oxygen storage materials with optimal performance. Such materials are required to have high stability, corrosion resistance, and acceptable oxygen storage capacity at moderate operating temperatures and pressures. Here, in this study, we used first-principles density functional theory and thermogravimetric analysis to study the vacancy-dependent oxygen diffusion in oxygen-deficient SrFeO 3−δ (δ = 0, 0.065, 0.125, 0.25, 0.5) perovskites. The electronic structures, including the partial- and spin-resolved density of states, for different SrFeO 3−δ phases were calculated and compared with available experimental and theoretical results. By mapping the migration pathways, we investigated diffusion mechanisms and calculated the energy barriers for oxygen diffusion in cubic, orthorhombic, and brownmillerite phases of SrFeO 3−δ perovskites. Using the calculated energy barriers, we deduced the diffusion time scales and diffusion coefficients within SrFeO 3−δ . A diffusion coefficient on the order of 10 –8 m 2 /s was obtained for SrFeO 2.875 . We experimentally investigated the roles of temperature and oxygen partial pressures on the redox kinetics and deduced the kinetics rate and diffusion density, which agreed well with the calculated values for the density of diffusing oxygen vacancy in the lattice. Our results showed that the energy barrier tends to reduce at higher oxygen concentrations. Our results serve as an important guideline for designing oxygen storage materials with optimal redox kinetics.

chemical looping with oxygen uncoupling (CLOU)↗

Coexistence of Interacting Charge Density Waves in a Layered Semiconductor

Coexisting orders are key features of strongly correlated materials and underlie many intriguing phenomena from unconventional superconductivity to topological orders. Here, we report the coexistence of two interacting charge-density-wave (CDW) orders in EuTe 4 , a layered crystal that has drawn considerable attention owing to its anomalous thermal hysteresis and a semiconducting CDW state despite the absence of perfect FS nesting. By accessing unoccupied conduction bands with time- and angle-resolved photoemission measurements, we find that mono- and bi-layers of Te in the unit cell host different CDWs that are associated with distinct energy gaps. The two gaps display dichotomous evolutions following photoexcitation, where the larger bilayer CDW gap exhibits less renormalization and faster recovery. Surprisingly, the CDW in the Te monolayer displays an additional momentum-dependent gap renormalization that cannot be captured by density-functional theory calculations. This phenomenon is attributed to interlayer interactions between the two CDW orders, which account for the semiconducting nature of the equilibrium state. Finally, our findings not only offer microscopic insights into the correlated ground state of EuTe 4 but also provide a general non-equilibrium approach to understand coexisting, layer-dependent orders in a complex system.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Photoionization of Benzene and Small Polycyclic Aromatic Hydrocarbons in Ultraviolet-Processed Astrophysical Ices: A Computational Study

We employed density functional theory (DFT) calculations to model the photoionization behavior of benzene and small polycyclic aromatic hydrocarbons when they are embedded in a matrix of water ice in order to investigate issues raised by recent experimental work by Gudipati and Allamandola. The ionization energies of benzene, naphthalene, anthracene, and pyrene were found to be lowered by 1.5-2.1 eV in water ice. Low-lying vertical electronic excitation energies were computed with time-dependent DFT for both neutral and ionized species and are found in both cases to be remarkably unaffected by the ice matrix. Chemical behavior in ultraviolet-photoprocessed ices is also discussed, with a focus on electron recombination and pathways leading to phenol and analogous products.

Woon, D. E.↗

Ripening of Rh Nanoparticle Catalysts in Reverse Water–Gas Shift via a Data-Driven Model Combining Physics, Theory, and Experiment

Degradation via sintering is an ongoing challenge that impedes the broad commercial success of supported metallic nanoparticle catalysts. To mitigate degradation via informed catalyst design and process operations, here we aim to disambiguate the underlying mechanisms of sintering by combining theory and experiment in a quantitative framework. While mechanistic sintering models exist, they only model a single sintering pathway, even though multiple sintering mechanisms can occur simultaneously or dominate at different stages of the process. Data-driven machine learning models have emerged as a means to represent complex processes through data regression. However, machine learning models have very large data needs and lack mechanistic insights due to their black-box encoding. To develop an interpretive model of catalyst degradation via sintering, we constructed a hybrid model combining mechanistic “physics-based” models and data-driven methods to obtain both reliable predictions and mechanistic insights regarding experimentally observed sintering phenomena. Focusing on nanoparticle sintering in the Rh–TiO 2 catalyst for the reverse water–gas shift (RWGS) reaction, the hybrid model couples a mechanistic term for Ostwald ripening with energy values calculated via density functional theory (DFT) with a parametric, data-driven discrepancy function term for unmodeled mechanisms. The hybrid model is trained using Bayesian inference with data collected from small-angle X-ray scattering (SAXS) in situ experiments wherein average nanoparticle diameter versus time was measured at three relevant operating temperatures. The calibrated hybrid model results show that an Ostwald ripening-only model parameterized with fixed DFT energies does not fully capture the time and temperature dependence of the SAXS-observed sintering kinetics, and that an additional functional contribution, or DFT energy calibration, is required to reconcile simulation and experiment. Analysis of the hybrid-model error confirms that the hybrid model outperforms both the purely mechanistic and purely data-driven alternatives in terms of expected predictive accuracy for time-evolving average particle sizes. Furthermore, the results support the hypothesis that the Ostwald ripening mechanism is less important for explaining the sintering phenomena as operating temperature increases under an assumed fixed DFT parameterization. This could be explained in one of two ways: either latent, unmodeled sintering mechanisms dominate at higher temperatures, or the DFT uncertainty increases with temperature. The proposed modeling approach directly links theory to experiments and simulations via a statistical hybrid modeling framework and can be extended to other catalytic systems to improve predictive models and mechanistic understanding.

Bayesian hybrid modeling↗

Speciation and diffusive dynamics in hydrated grain boundaries of complex oxide Gd2Ti2O7

Abstract Grain boundaries in polycrystalline materials significantly affect their properties, such as ionic transport, corrosion, and chemical durability. The pyrochlore compound (Gd 2 Ti 2 O 7 ) is employed as a model for complex oxides and is known for its diverse applications, including nuclear waste immobilization. Density functional theory-based first-principles molecular dynamics simulations were performed at different temperatures on the hydrated grain boundary system. The results show extensive transformations within the grain boundaries among hydrous water species (OH − , H 2 O, and H 3 O + ). The temperature dependence of self-diffusion coefficients follows Arrhenius behavior, with an activation energy of 35.9 kJ/mol for hydrogen and 46.3 kJ/mol for oxygen. The lifetime of OH − is about three to four times longer than that of H 2 O at temperatures from 800 to 2100 K, suggesting the greater stability of OH − over H 2 O, a unique characteristic of the grain boundaries. The estimated lifetime of the hydrous species decreases as the temperature increases, with an activation energy of 9.9 kJ/mol for OH − and 13.4 kJ/mol for H 2 O. While Gd 3 + is more mobile than Ti 4+ , both the Gd 3 + and Ti 4+ cations are orders of magnitude less mobile than the water species. The results suggest that water species are much more mobile within grain boundaries than in the bulk crystal and have the potential to penetrate deep into polycrystalline materials through grain boundaries, leading to grain boundary degradation and dissolution. The different mobilities of cations in complex oxides can lead to leaching of certain cations and incongruent dissolution during the chemical weathering of Earth and industrial materials.

B. Ghosh, Dipta↗

Utilizing local phase transformation strengthening for nickel-based superalloys

Almost 75 years of research has been devoted to producing superalloys capable of higher operating temperatures in jet turbine engines, and there is an ongoing need to increase operating temperature further. Here, a new disk Nickel-base superalloy is designed to take advantage of strengthening atomic-scale dynamic complexions. This local phase transformation strengthening provides the alloy with a three times improvement in creep strength over similar disk superalloys and comparable strength to a single crystal blade alloy at 760 °C. Ultra-high-resolution chemical mapping reveals that the improvement in creep strength is a result of atomic-scale η (D024) and χ (D019) formation along superlattice stacking faults. To understand these results, the energy differences between the L12 and competing D024 and D019 stacking fault structures and their dependence on composition are computed by density functional theory. This study can help guide researchers to further optimize local phase transformation strengthening mechanisms for alloy development.

Timothy M. Smith↗

Ligand Controls Excited Charge Carrier Dynamics in Metal-Rich CdSe Quantum Dots: Computational Insights

Small metal-rich semiconducting quantum dots (QDs) are promising for solid-state lighting and single-photon emission due to their highly tunable yet narrow emission line widths. Nonetheless, the anionic ligands commonly employed to passivate these QDs exert a substantial influence on the optoelectronic characteristics, primarily owing to strong electron–phonon interactions. In this work, we combine time-domain density functional theory and nonadiabatic molecular dynamics to investigate the excited charge carrier dynamics of Cd 28 Se 17 X 22 QDs (X = HCOO – , OH – , Cl – , and SH – ) at ambient conditions. These chemically distinct but regularly used molecular groups influence the dynamic surface-ligand interfacial interactions in Cd-rich QDs, drastically modifying their vibrational characteristics. The strong electron–phonon coupling leads to substantial transient variations at the band edge states. The strength of these interactions closely depends on the physicochemical characteristics of passivating ligands. Consequently, the ligands largely control the nonradiative recombination rates and emission characteristics in these QDs. Our simulations indicate that Cd 28 Se 17 (OH) 22 has the fastest nonradiative recombination rate due to the strongest electron–phonon interactions. Conversely, QDs passivated with thiolate or chloride exhibit considerably longer carrier lifetimes and suppressed nonradiative processes. The ligand-controlled electron–phonon interactions further give rise to the broadest and narrowest intrinsic optical line widths for OH and Cl-passivated single QDs, respectively. Finally, obtained computational insights lay the groundwork for designing appropriate passivating ligands on metal-rich QDs, making them suitable for a wide range of applications, from blue LEDs to quantum emitters.

36 MATERIALS SCIENCE↗

Eigenfunctions of the magnetospheric radial-diffusion operator

The applicability of eigenfunction methods to the theory of magnetospheric radial diffusion is investigated analytically, with a focus on the operator (Lambda) for radiation-belt and ring-current particles. It is shown that the eigenfunctions of Lambda can be expressed in terms of a linear combination of ordinary Bessel functions of the first and second kinds, and that the solution for the drift-averaged phase-space density can be expressed as the sum of a quasi-static solution and a time-dependent superposition of the eigenfunctions of Lambda. Results for illustrative and practical sample problems are presented in extensive tables and graphs and discussed in detail.

Schulz, Michael↗

The total neutron cross section of liquid and solid ammonia

Ammonia is a material of interest for future neutron moderators at high-power sources due to its high hydrogen density, low melting point, and resistance to polymerization in an intense radiation field. Its performance in such applications cannot currently be calculated due to the absence of suitable computer models for the interaction of neutrons with ammonia under relevant conditions. In an effort to develop suitable scattering kernels for computer simulations of moderator performance, we have conducted a series of Density Functional Theory and Molecular Dynamics calculations of the molecular-level thermal properties of ammonia at various temperatures within both the solid and liquid phases. In this paper, we compare computer calculations for the energy-dependent total neutron cross section of ammonia, based on these models, to experimental measurements of those cross sections at temperatures of 221 K, 180 K, and 35 K. The experimental data were collected over an energy range from 0.1 meV to 10 eV using time-of-flight techniques at the Low Energy Neutron Source (LENS) facility at Indiana University. This comparison provides a first validation in the development of thermal scattering libraries for Monte Carlo source design simulations based on liquid and solid ammonia. In conclusion, we also provide some insights into where additional development of tools for creating such models may be needed.

Ammonia↗

Fe‐Triazolate Metal–Organic Frameworks as Water Oxidation Catalysts with Dual Photoanode Functionality

Artificial photosynthesis is an emerging technology that achieves renewable fuels, such as hydrogen, from sunlight. Its realization depends on finding highly active and stable catalysts of water splitting and photoactive materials for light absorption. To be scalable, these should contain only abundant elements. Here, for the first time, Fe-triazolate (Fe(ta) 2 ) and its metal substituted derivatives (Fe-Metal(ta) 2 ) Metal-organic frameworks (MOFs) are characterized as new dual-function materials for photo-absorption and water oxidation catalysis in acidic media. The materials were studied by a range of structural, spectroscopic, and computational density functional theory (DFT) techniques. Fe(ta) 2 and Fe-Mn(ta) 2 were found to be highly active and stable in chemical and photochemical water oxidation, and in addition function as photoanodes, with photo-electrocatalytic currents (∼2.00 x 10 −3 Acm −2 at + 1.4 V vs. Ag/AgCl) at pH = 1. The possibility of a unique catalytic mechanism where O─O bond formation is possible from the coupling of two adjacent Fe IV = O fragments was demonstrated by DFT analysis. Thus, Fe-triazolate MOF has been established as a new, stable, scalable, versatile, and efficient platform for sustainable energy conversion in the realm of artificial photosynthesis.

Artificial photosynthesis↗

Nature of molybdenum carbide surfaces for catalytic hydrogen dissociation using machine-learned potentials: an ensemble-averaged perspective

Molybdenum carbides with an electronic structure similar to noble metals have gained attention as a promising low-cost catalyst for biomass valorization and the hydrogen evolution reaction. However, our fundamental understanding of the catalyst surface and how different phases of these catalysts behave at varying reaction conditions is limited to ground state density functional theory calculations as ab initio molecular dynamics (AIMD) is computationally prohibitive at relevant length and time scales. Here, in this work, we train a multi-atomic cluster expansion (MACE) machine-learned interatomic potentials (MLIP) to study hydrogen dissociation and dynamics over Mo, δ-MoC, α-Mo 2 C, and β-Mo 2 C surfaces at varying temperatures and hydrogen partial pressures. Our simulations identify unique and different molecular and atomic hydrogen adsorption sites on different surfaces that do not depend on the temperature. At low hydrogen pressures, the surface coverage is monolayer, which transitions to two-layer adsorption at higher pressures. We find that atomic hydrogen diffusion and recombinations are preferred over molybdenum atom hollow sites, while the diffusion over carbon-terminated facets was negligible, signifying particularly strong C–H interactions. In contrast, molecular hydrogen adsorption occurs mostly atop Mo or the bridging sites. At a comparable hydrogen loading, β-Mo 2 C (001) is the most active surface for hydrogen dissociation reaction. This work provides insights into the dynamic nature of the hydrogen dissociation chemistry and the diversity of hydrogen adsorption sites on molybdenum carbides.

08 HYDROGEN↗

Influence of Markovianity and self-consistency on time-resolved spectral functions of driven quantum systems

We present a systematic comparison of the real-time Dyson expansion (RTDE) with established nonequilibrium Green's function (GF) approaches for simulating driven, interacting quantum systems. Focusing on density matrix dynamics, time-off-diagonal GFs, and time-resolved photoemission spectra, we benchmark RTDE against fully self-consistent Kadanoff-Baym equation (KBE) calculations, the generalized Kadanoff-Baym ansatz, and exact diagonalization for small systems using second-order many-body perturbation theory. Using a driven two-band Hubbard model, we show that mean-field single-particle density matrix trajectories provide a reliable baseline for RTDE across a broad range of interaction strengths and excited-carrier populations. Further, RTDE accurately captures correlation effects in the GFs, including long-lived oscillations and revivals that are strongly suppressed by the overdamping inherent to self-consistent KBE schemes. As a consequence, RTDE resolves rich nonequilibrium spectral structure in time-resolved photoemission, such as interaction- and population-dependent quasiparticle splittings and band gap renormalization, which are largely washed out in self-consistent approaches yet are present in exact solutions. Furthermore, our results demonstrate that RTDE bridges the gap between mean-field propagation and full two-time KBE simulations, retaining favorable linear scaling while capturing essential dynamical correlations relevant for ultrafast spectroscopy.

Electronic structure↗

The TDHF code Sky3D version 1.2

The Sky3D code has been widely used to describe nuclear ground states, collective vibrational excitations, and heavy-ion collisions. The approach is based on Skyrme forces or related energy density functionals. The static and dynamic equations are solved on a three-dimensional grid, and pairing is been implemented in the BCS approximation. This updated version of the code aims to facilitate the calculation of nuclear strength functions in the regime of linear response theory, while retaining all existing functionality and use cases. The strength functions are benchmarked against available RPA codes, and the user has the freedom of choice when selecting the nature of external excitation (from monopole to hexadecapole and more). Some utility programs are also provided that calculate the strength function from the time-dependent output of the dynamic calculations of the Sky3D code.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Spontaneous emergence of straintronics effects and striped stacking domains in untwisted three-layer epitaxial graphene

Emergent electronic phenomena, from superconductivity to ferroelectricity, magnetism, and correlated many-body band gaps, have been observed in domains created by stacking and twisting atomic layers of Van der Waals materials. In graphene, emergent properties have been observed in ABC stacking domains obtained by exfoliation followed by expert mechanical twisting and alignment with the desired orientation, a process very challenging and nonscalable. Here, conductive atomic force microscopy shows in untwisted epitaxial graphene grown on SiC the surprising presence of striped domains with dissimilar conductance, a contrast that demonstrates the presence of ABA and ABC domains since it matches exactly the conductivity difference observed in ABA/ABC domains in twisted exfoliated graphene and calculated by density functional theory. The size and geometry of the stacking domains depend on the interplay between strain, solitons crossing, and shape of the three-layer regions. Interestingly, we demonstrate the growth of three-layer regions in which the ABA/ABC stacking domains self-organize in stable stripes of a few tens of nanometers. The growth-controlled production of isolated and stripe-shaped ABA/ABC domains open the path to fabricate quantum devices on these domains. These findings on self-assembly formation of ABA/ABC epitaxial graphene stripes on SiC without the need of time-consuming and nonscalable graphene exfoliation, alignment, and twisting provide different potential applications of graphene in electronic devices.

Rejhon, Martin (ORCID:000000017775487X)↗

High-temperature structure, elasticity, and thermal expansion of ε-ZrH 1.8

Zirconium hydride is a promising candidate material for nuclear microreactor applications as a solid-state moderator component, owing to its favorable neutronics properties and good thermal stability over other metal hydrides. Here, in the present work, the crystal structure, thermal expansion, and elastic properties of the hydrogen-rich ε phase hydride were measured at elevated temperatures in the range 300–900 K. Samples were prepared by direct hydriding Zircaloy-4 metal – a nuclear-grade zirconium alloy. Room-temperature lattice parameters agree well with those reported from literature for unalloyed zirconium hydride and fall within an observed quadratic H-content dependence. The coefficients of thermal expansion, determined from lattice expansion and dilatometry, agree well within our work but were about 30 % lower than those reported by others for unalloyed hydrides. Density functional theory-based molecular dynamics simulations were used to compare with thermal expansion and elasticity measurements. Results showed lattice parameter temperature dependence and slope of thermal expansion align with those from measurements. Based on diffraction scans at select temperatures, ε phase remained stable in air up to at least 770 K. Likewise, dilatometry showed smooth thermal expansion up to the thermal decomposition temperature around 950 K. The precise decomposition temperature was not determined via diffraction due to sparse scanning. The complete elastic property measurements were gathered for ε-phase Ziracloy-4 hydride for the first time. Young's modulus was lower compared to the metal and δ hydride phases. High-temperature elasticity measurements were limited to <350 K due to acoustic dissipation effects.

36 MATERIALS SCIENCE↗

A topological Hund nodal line antiferromagnet

The interplay of topology, magnetism, and correlations gives rise to intriguing phases of matter. In this study, through state-of-the-art angle-resolved photoemission spectroscopy, density functional theory, and dynamical mean-field theory calculations, we visualize a fourfold degenerate Dirac nodal line at the boundary of the bulk Brillouin zone in the antiferromagnet YMn 2 Ge 2 . We further demonstrate that this gapless, antiferromagnetic Dirac nodal line is enforced by the combination of magnetism, space-time inversion symmetry, and nonsymmorphic lattice symmetry. The corresponding drumhead surface states traverse the whole surface Brillouin zone. YMn 2 Ge 2 thus serves as a platform to exhibit the interplay of multiple degenerate nodal physics and antiferromagnetism. Interestingly, the magnetic nodal line displays a d-orbital dependent renormalization along its trajectory in momentum space, thereby manifesting Hund's coupling. Our findings offer insights into the effect of electronic correlations on magnetic Dirac nodal lines, leading to an antiferromagnetic Hund nodal line.

36 MATERIALS SCIENCE↗