SEARCH · Engineering Papers
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.
A new electron localization function is justified and compared with the laplacian of the electronic charge density. Some advantages of the new function and its novel predictions are discussed. Several slides and some video of visualizations using the new function in comparison to the laplacian will be shown.
Hydrogen‐based materials are able to possess extremely high superconducting critical temperatures, T c s , due to hydrogen's low atomic mass and strong electron–phonon interaction. Recently, a descriptor based on the Electron Localization Function (ELF) has enabled the rapid estimation of the T c of hydrogen‐containing compounds from electronic networking properties, but its applicability has been limited by the small size and homogeneity of the training dataset used. Herein, the model is re‐examined, compiling a publicly available combined dataset of 244 binary and ternary hydride superconductors. The analysis shows that though ELF‐based networking remains a valuable descriptor, its predictive power declines with increasing compositional complexity. However, by introducing the molecularity index, defined as the highest value of the ELF at which two hydrogen atoms connect, and applying symbolic regression, the accuracy of the predictions can be substantially enhanced. These results establish a more robust framework for assessing superconductivity in hydride materials, facilitating accelerated screening of novel candidates through integration with crystal structure prediction methods or high‐throughput searches.
It has long been assumed that all matter will adopt simple close-packed lattices and become metallic under pressure, in accordance with the Thomas–Fermi–Dirac (TFD) model. However, this model struggles to explain pressure-driven complex structural transitions that have been observed in elements, including sodium, challenging our conventional understanding of compressed matter. Moreover, in stark contrast to the TFD model, first-principles calculations suggest that various elements and compounds become electrides under pressure. Electrides, characterized by concentrations of charge density at interstitial regions, can be thought of as ionic compounds where electrons behave as the anions. Though ambient-pressure molecular electrides have been extensively studied via experiments and computations, high-pressure electrides (HPEs) are not well-understood. The identification and characterization of HPEs have been, to date, based purely on theory, including topological analysis of the electron density and the electron localization function. Here, we review these theoretical analysis tools and suggest guidelines that can be used to classify systems as electrides. Moreover, we describe models used to rationalize the electronic structure of HPEs, drawing parallels with ambient-pressure molecular systems, and encourage the development of experimental techniques that provide evidence for the theoretically calculated charge localization.
A combination of ab initio, atomistic and finite element methods (FEM) were used to investigate the structures, energetics and lattice thermal conductance of grain boundaries for the ultra high temperature ceramic ZrB2. Atomic models of idealized boundaries were relaxed using density functional theory. Information about bonding across the interfaces was determined from the electron localization function. The Kapitza conductance of larger scale versions of the boundary models were computed using non-equilibrium molecular dynamics. The interfacial thermal parameters together with single crystal thermal conductivities were used as parameters in microstructural computations. FEM meshes were constructed on top of microstructural images. From these computations, the effective thermal conductivity of the polycrystalline structure was determined.
While hydrogen-rich materials have been demonstrated to exhibit high T c superconductivity at high pressures, there is an ongoing search for ternary, quaternary, and more chemically complex hydrides that achieve such high critical temperatures at much lower pressures. First-principles searches are impeded by the computational complexity of solving the Eliashberg equations for large, complex crystal structures. Here, we adopt a simplified approach using electronic indicators previously established to be correlated with superconductivity in hydrides. This is used to study complex hydride structures, which are predicted to exhibit promisingly high critical temperatures for superconductivity. In particular, we propose three classes of hydrides inspired by the Fm$\overline{3}$ m RH 3 structures that exhibit strong hydrogen network connectivity, as defined through the electron localization function. The first class [RH 11 X 3 Y] is based on a Pm$\overline{3}$m structure showing moderately high T c , where the T c estimate from electronic properties is compared with direct Eliashberg calculations and found to be surprisingly accurate. The second class of structures [(RH 11 ) 2 X 6 YZ] improves on this with promisingly high density of states with dominant hydrogen character at the Fermi energy, typically enhancing T c . The third class [(R 1 H 11 )(R 2 H 11 )X 6 YZ] improves the strong hydrogen network connectivity by introducing anisotropy in the hydrogen network through a specific doping pattern. These design principles and associated model structures provide flexibility to optimize both T c and the structural stability of complex hydrides.
This study employs Density Functional Theory (DFT) and Molecular Dynamics (MD) simulations to investigate interactions between water molecules and Poly(N-isopropylacrylamide) (PNIPAM). DFT reveals preferential water binding sites, with enhanced binding energy observed in the linker zone. Quantum Theory of Atoms in Molecules (QTAIM) and electron localization function (ELF) analyses highlight the roles of hydrogen bonding and steric hindrance. MD simulations unveil temperature-dependent hydration dynamics, with structural transitions marked by changes in the radius of gyration (Rg) and the radial distribution function (RDF), aligning with DFT findings. Our work goes beyond prior studies by combining a DFT, QTAIM and MD simulations approach across different PNIPAM monomer-to-30mer structures. It introduces a systematic quantification of pseudo-saturation thresholds and explores water clustering dynamics with structural specificity, which have not been previously reported in the literature. These novel insights establish a more complete molecular-level picture of PNIPAM hydration behavior and temperature responsiveness, emphasizing the importance of amide hydrogen and carbonyl oxygen sites in hydrogen bonding, which weakens above the lower critical solution temperature (LCST), resulting in increased hydrophobicity and paving the way for understanding water sorption mechanisms, offering guidance for future applications such as dehumidification and atmospheric water harvesting.
This dataset contains DFT input and output files supporting the theoretical modeling in the associated publication (Chem. Mater. 2024, 36, 11306-11316). The calculations characterize MnRhBi3, an orthorhombic (Cmmm) van der Waals-layered intermetallic compound that cleaves easily between neighboring Bi layers. The dataset is organized into three calculation types: (i) Bulk: Structural relaxations of the periodic MnRhBi3 crystal in antiferromagnetic (AFM) and ferromagnetic (FM) configurations, using the vdW-DF-optB86b functional. These provide the equilibrium lattice constants, magnetic energy differences (AFM is 0.5 meV/f.u. lower than FM), and magnetic moments (4.4 µB/Mn, 0.17 µB/Rh, 0.18 µB/Bi) reported in Table 1 of the main text. (ii) Slab: Same magnetic configurations computed with an 18 Ang vacuum layer introduced between Bi layers, used to calculate the cleavage energy Ec = 0.56 J/m2 (AFM) and 0.57 J/m2 (FM), establishing MnRhBi3 as a van der Waals-layered material comparable to graphite, MoS2, and CrI3. (iii) ELF: Single-point calculation on the relaxed bulk AFM geometry with LELF=.TRUE., producing the ELFCAR file used to generate electron localization function isosurfaces and contour maps (Fig. 2, main text) showing Bi lone pairs directed into the van der Waals gaps. All folders contain CONTCAR, INCAR, KPOINTS, OUTCAR, and POSCAR. The ELF/ folder additionally contains ELFCAR. Calculations were performed using VASP 6.3.2 with PBE + vdW-DF-optB86b, PAW potentials, and an energy cutoff of 800 eV.
Four new d 0 transition metal tellurites, namely InVTeO 6 (1), β-InVTe 2 O 8 (2), FeVTe 2 O 8 (3), and Fe 2 MoTe 2 O 10 (4), were synthesized under hydrothermal conditions, and their crystal structures were determined using single-crystal X-ray diffraction. Compound 1 crystallizes in the noncentrosymmetric space group P2 1 2 1 2 1 (No. 19), with unit cell parameters of a = 5.0759(2) Å, b = 8.5030(3) Å, c = 11.6376(5) Å, V = 502.28(3) Å 3 , and Z = 4, while the other three compounds crystallize in centrosymmetric structures, with space group Cmca and unit cell parameters of a = 7.0633(7) Å, b = 8.962(1) Å, c = 20.162(2) Å, V = 1276.3(2) Å 3 , and Z = 8 for compound 2, space group P2 1 /n and unit cell parameters of a = 7.8901(2) Å, b = 4.9617(1) Å, c = 16.4290(3) Å, β = 93.6198(8)°, V = 641.88(2) Å 3 , and Z = 4 for compound 3, and space group Pnma and unit cell parameters of a = 8.7878(2) Å, b = 6.1327(2) Å, c = 15.2423(4) Å, V = 821.45(4) Å 3 , and Z = 4 for compound 4. Compound 3 adopts a two-dimensional (2D) sheet structure composed of FeO 6 octahedra, VO 4 tetrahedra, and TeO 4 polyhedra. The other three compounds exhibit three-dimensional (3D) framework structures: Compound 1 is built from InO 6 octahedra, VO 4 tetrahedra, TeO 4 polyhedra; Compound 2 consists of InO 6 octahedra, VO 5 square pyramids, and TeO 3 polyhedra; and Compound 4 is composed of FeO 6 octahedra, MoO 6 octahedra, and TeO 3 polyhedra. The band structures, density of states, and electron localization functions of all four compounds were calculated. The vibrational spectra were examined for compounds 1, 2, and 4. Furthermore, powder X-ray diffraction, element analysis, and thermogravimetric analysis were conducted for compound 4.
Phase-pure CuIn 2 Se 4 , a ternary metal chalcogenide that forms in a disordered stannite crystal structure, was synthesized to investigate the structure thermal property relationships as well as reveal the origin of the ultralow thermal conductivity this material possesses over a large temperature range. Modeling of the temperature-dependent heat capacity and thermal conductivity revealed distinctive thermal properties and large lattice anharmonicity. Electron localization function calculations highlight the asymmetric bonding inherent to CuIn 2 Se 4 , which together with lattice anharmonicity directly impacts the thermal properties. Our findings reveal the specific atomic arrangement and bonding governing the thermal properties of this ternary metal chalcogenide. Our findings underscore the specific atomic arrangement and bonding governing the thermal properties in this ternary chalcogenide. This study advances the fundamental understanding of stannites, and our findings can be applied to these and other multinary metal chalcogenides of interest for applications where low thermal conductivity is desirable.
Metal silicon phosphides composed of earth-abundant Si and P tend to exhibit semiconducting properties and adopt diverse crystal structures with relatively small additions of structure-directing elements. The potential of silicon phosphide materials in nonlinear optical applications has been hindered by the inability to systematically produce noncentrosymmetric structures with such a flexible framework. Here, in this work, two isostructural compounds with a novel noncentrosymmetric structure were made possible by the inclusion of elements with stereochemically active lone pairs (Sn 2+ and Pb 2+ ). The structures were determined through single-crystal and synchrotron powder X-ray diffraction. Analysis of chemical bonding in real space through the electron localization function revealed stereochemically active Pb 2+ and Sn 2+ species in a trigonal pyramidal coordination with {Pb/Sn}–P bonds. Such covalent bonding between Pb and P is quite uncommon in extended solids and has been reported in a few rare instances. Band structure calculations and linear optical measurements confirm the semiconducting nature of Cs X Si 15 P 21 ( X = Sn or Pb). The synthesis was optimized to yield high-purity polycrystalline samples. The nonlinear optical properties show promising second-harmonic generation (SHG) coefficients from the Kurtz–Perry method. First-principles calculations of the nonlinear optical properties support the experimentally determined SHG values and provide moderate values of birefringence, suggesting Cs X Si 15 P 21 could be phase-matchable and practical nonlinear optical materials in the mid-IR region.
We provide an exact evaluation of the out-of-time correlation (OTOC) functions for the localized 𝑓-particle states in the Falicov-Kimball model within dynamical mean-field theory. Different regimes of quantum chaos and quantum scrambling are distinguished by the winding numbers of the block Toeplitz matrices used in the calculation. The similarities of these fermionic OTOCs and their logarithmic derivatives for time evolution with the OTOCs for quantum spin models with disorder are also discussed.
The Topaz3 auroral sounding rocket made the following observations concerning the transfer of precipitating auroral electron energy to transverse ion acceleration in the topside auroral zone. During the course of the flight, the precipitating electron beam was modified to varying degrees by interaction with VLF hiss, at times changing the beam into a field-aligned plateau. The electron distribution functions throughout the flight are classified according to the extent of this modification, and correspondences with ion acceleration events are sought. The hiss power during most of this rocket flight apparently exceeded the threshold for collapse into solitary structures. At the times of plateaued electron distributions, the collapse of these structures was limited by Landau damping through the ambient ions, resulting in a velocity-dependent acceleration of both protons and oxygen. This initial acceleration is sufficient to supply the number flux of upflowing ions observed at satellite altitudes. The bursty ion acceleration was anticorrelated, on 1-s or smaller timescales, with dispersive bursts of precipitating field-aligned electrons, although on longer timescales the bursty ions and the bursty electrons are correlated.
The effect of electron-nuclear spin interactions on qubit operations is investigated for a qubit represented by the spin of an electron localized in an InGaAs self-assembled quantum dot. The localized electron wave function is evaluated within the atomistic tight-binding model. The electron Zeeman splitting induced by the electron-nuclear spin interaction is estimated in the presence of an inhomogeneous environment characterized by a random nuclear spin configuration, by the dot-size distribution, alloy disorder, and interface disorder. Due to these inhomogeneities, the electron Zeeman splitting varies from one qubit to another by the order of 10(-6), 10(-6), 10(-7), and 10(-9) eV, respectively. Such fluctuations cause errors in exchange operations due to the inequality of the Zeeman splitting between two qubits. However, the error can be made lower than the quantum error threshold if an exchange energy larger than 10(-4) eV is used for the operation. This result shows that the electron-nuclear spin interaction does not hinder quantum-dot based quantum computer architectures from being scalable even in the presence of inhomogeneous environments.
A detailed first principle kinetic theory for electrons which is neither a classical fluid treatment nor an exospheric calculation is presented. This theory illustrates the global and local properties of the solar wind expansion that shape the observed features of the electron distribution function, such as its bifurcation, its skewness and the differential temperatures of the thermal and suprathermal subpopulations. Coulomb collisions are substantial mediators of the interplanetary electron velocity distribution function and they place a zone for a bifurcation of the electron distribution function deep in the corona. The local cause and effect precept which permeates the physics of denser media is modified for electrons in the solar wind. The local form of transport laws and equations of state which apply to collision dominated plasmas are replaced with global relations that explicitly depend on the relative position of the observer to the boundaries of the system.
Here, we introduce the Wannier function decomposition of excitons (WFDX) method to quantify exciton localization in solids within the ab initio Bethe-Salpeter equation framework. By decomposing each Bloch exciton wave function into products of single-particle electron and hole maximally localized Wannier functions, this real-space approach provides well-defined orbital- and spatial-resolved measures of both Frenkel and charge-transfer excitons at low computational cost. We apply WFDX to excitons in acene crystals, quantifying how the number of rings, the exciton spin state, and the center-of-mass momentum affect spatial localization. Additionally, we show how this real-space representation reflects structural nonsymmorphic symmetries that are hidden in standard reciprocal-space descriptions. We demonstrate how the WFDX framework can be used to efficiently interpolate exciton expansion coefficients in reciprocal-space and outline how it may facilitate evaluation of observables involving position operators, highlighting its potential as a general tool for both analyzing and computing excitonic properties in solids.
The very high symmetry of icosahedral C60 suggests that, as a first approximation, an atom trapped inside C60 would be subject to a potential that is radially symmetric about the center. All-electron local-density-functional calculations of the total energy of a sodium ion as a function of radial displacement from the center along the fivefold axis of C60 serve to refine such a radial potential. In particular, the calculations suggest studying potentials that have minima displaced from the center. An analytic functional form for a radial potential having a positive cusp at the origin is proposed, and the s-wave radial solutions of the corresponding Schroedinger equation are examined.
Measurements from the Magnetospheric Multiscale (MMS) mission are reported to show distinct features of electron energization and mixing in the diffusion region of the terrestrial magnetopause reconnection. At the ion jet and magnetic field reversals, distribution functions exhibiting signatures of accelerated meandering electrons are observed at an electron out-of-plane flow peak. The meandering signatures manifested as triangular and crescent structures are established features of the electron diffusion region (EDR). Effects of meandering electrons on the electric field normal to the reconnection layer are detected. Parallel acceleration and mixing of the inflowing electrons with exhaust electrons shape the exhaust flow pattern. In the EDR vicinity, the measured distribution functions indicate that locally, the electron energization and mixing physics is captured by two-dimensional reconnection, yet to account for the simultaneous four-point measurements, translational invariant in the third dimension must be violated on the ion-skin-depth scale.
Viscosity coefficient curve fits for neutral gas species are available from many sources. Many do a good job of reproducing experimental and computational chemistry data. The curve fits are usually expressed as a function of temperature only. This is consistent with the governing equations used to derive an expression for the neutral species viscosity coefficient. Ionized species pose a more complicated problem. They are subject to electrostatic as well as intermolecular forces. The electrostatic forces are affected by a shielding phenomenon where electrons shield the electrostatic forces of positively charged ions beyond a certain distance. The viscosity coefficient for an ionized gas species is a function of both temperature and local electron number density. Currently available curve fits for ionized gas species, such as those presented by Gupta/Yos, are a function of temperature only. What they did was to assume an electron number density. The problem is that the electron number density they assumed was unrealistically high. The purpose of this paper is two-fold. First, the proper expression for determining the viscosity coefficient of an ionized species as a function of both temperature and electron number density will be presented. Then curve fit coefficients will be developed using the more realistic assumption of an equilibrium electron number density. The results will be compared against previous curve fits and against highly accurate computational chemistry data.