Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “covalency”

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 1,351 records · Page 75
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Interplay between electron localization, magnetic order, and Jahn-Teller distortion dictates LiMnO2 phase stability

The development of manganese (Mn)-rich cathodes for Li-ion batteries promises to alleviate potential supply chain bottlenecks in battery manufacturing. Fundamental challenges in Mn-rich cathodes arise from phenomena such as structural changes due to cooperative Jahn-Teller (JT) distortions of in octahedral environments, Mn migration, and phase transformations to spinel-like order, all of which affect the electrochemical performance. These physically complex phenomena motivate an re-examination of the Li-Mn-O rock-salt space, with a focus on the thermodynamics of the prototypical, polymorphs. It is found that the generalized gradient approximation (GGA-PBEsol) and meta-GGA ( ) density functionals with empirically fitted on-site Hubbard corrections yield spurious stable phases for , such as predicting a phase with -like order ( ) to be the ground state instead of the orthorhombic (Pmmn) phase, which is the experimentally known ground state. Accounting for antiferromagnetic order in each structure is shown to have a substantial effect on the total energies and resulting phase stability. By using hybrid-GGA (HSE06) and GGA with self-consistent Hubbard parameters (on-site and inter-site ) calculated from linear response theory, the experimentally observed phase stability trends are recovered. The calculated on-site between Mn- states in the experimentally observed orthorhombic, layered, and spinel phases are significantly smaller than in and disordered layered structures, by within GGA. The smaller values of are shown to be correlated with a collinear ordering of JT distortions, in which all orbitals are oriented in the same direction. This cooperative JT effect can lead to greater electron delocalization from Mn along the states due to increased Mn-O covalency, which contributes to the greater electronic stability compared to the phases with noncollinear JT arrangements. The structures with collinear ordering of JT distortions also generate greater vibrational entropy, which helps stabilize these phases at high temperature. These phases are shown to be strongly insulating with large calculated band gaps , which are computed using HSE06 and .

Kam, Ronald L↗
DOE OSTI2025

Carbon-carbon interactions in warm dense titanium carbide

Encouraged by experimental reports of diamond precipitation in carbon-containing materials under warm dense matter conditions and in shock compression experiments, we examine the behavior of carbon in TiC using density functional theory-molecular dynamics simulations. Two polymorphs of TiC are considered, the ambient-pressure B1 ($Fm$$\overline{3}$$m$) structure in which C-C interactions are prevented by geometric constraints, and a high-pressure Cmcm structure in which C atoms condense into zigzag chains. Chemistry-inspired bonding analyses confirm the covalently bound nature of these chains and further illuminate important interatomic interactions in both structures. Upon melting of B1 TiC, new short-range C-C interactions develop in the liquid, while the pre-existing C-C interactions in Cmcm TiC persist in the liquid. Here, the resulting carbon networks in the melt provide a promising environment for eventual diamond nucleation.

Ab initio molecular dynamics↗
DOE OSTI2020

Topological spin texture in Janus monolayers of the chromium trihalides Cr(I, X ) 3

Topological magnetic states are promising for ultradense memory and logic devices. Recent progress in two-dimensional magnets encourages the idea to realize topological states, such as skyrmions and merons, in freestanding monolayers. However, monolayers such as CrI 3 lack Dzyaloshinskii-Moriya interactions (DMIs) and thus do not naturally exhibit skyrmions/merons but rather a ferromagnetic state. Here we propose the fabrication of Cr(I, $X)_3$ Janus monolayers, in which the Cr atoms are covalently bonded to the underlying I ions and top-layer Br or Cl atoms. By performing first-principles calculations and Monte Carlo simulations, we identify strong enough DMIs, which leads to not only helical cycloid phases, but also to topologically nontrivial states, such as the intrinsic domain wall skyrmions in Cr(I, Br) 3 and the magnetic-field-induced bimerons in Cr(I, Cl) 3 . Finally, microscopic origins of such spin textures are revealed as well.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Temperature-induced structural evolution in liquid Ag-Ga alloys

Temperature-dependent atomic structural evolutions of liquid Ag 60 Ga 40 and Ag 70 Ga 30 alloys have been studied by in situ high-energy x-ray-diffraction (HEXRD) experiments combined with ab initio molecular-dynamics simulations. The experimental data show a reversible structural crossover at about 1050 ~ 1100 (± 50) K in both liquid Ag 60 Ga 40 and Ag 70 Ga 30 alloys. Obvious changes of the electrical resistivity, absolute thermoelectric power, and atomic diffusivity around the similar temperature range for both Ag-Ga liquids strongly support the HEXRD results. The origin of the liquid-to-liquid crossover in both Ag-Ga liquids was suggested to link with the rearrangements of Ag and Ga atoms, i.e., Ag and Ga atoms prefer to associate with themselves in the higher temperature range above 1100 K, consistent with the accelerated increase of the strong covalently bonded Ga–Ga dimers in both Ag-Ga liquids. In addition, more studies from the energy aspect are still desirable to understand the rearrangements of Ag and Ga atoms in the higher temperature range in both Ag-Ga liquids.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Empty perovskites as Coulomb floppy networks: Entropic elasticity and negative thermal expansion

Floppy networks (FNs) provide valuable insight into the origin of anomalous mechanical and thermal properties in soft matter systems, from polymers, rubber, and biomolecules to glasses and granular materials. Here, we use the same FN concept to construct a quantitative microscopic theory of empty perovskites, a family of crystals with ReO 3 structure, which exhibit a number of unusual properties. One remarkable example is ScF 3 , which shows a near-zero-temperature structural instability and large negative thermal expansion (NTE). Furthermore, we trace these effects to an FN-like crystalline architecture formed by strong nearest-neighbor bonds, which is stabilized by net electrostatic repulsion that plays a role similar to osmotic pressure in polymeric gels. NTE in these crystalline solids, which we conceptualize as Coulomb floppy networks, emerges from the tension effect of Coulomb repulsion combined with the FN's entropic elasticity and has the same physical origin as in gels and rubber. Our theory provides an accurate, quantitative description of phonons, thermal expansion, compressibility, and structural phase diagram, all in excellent agreement with experiments. The entropic stabilization of critical soft modes, which play only a secondary role in NTE, explains the observed phase diagram. Significant entropic elasticity resolves the puzzle of a marked, ≈50% discrepancy between the experimentally observed bulk modulus and ab initio calculations. The Coulomb FN approach is potentially applicable to other important materials with markedly covalent bonds, from perovskite oxides to iron chalcogenides, whose anomalous vibrational and structural properties are still poorly understood.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

First-principles study of superconductivity in α and β gallium

Elemental gallium can exist in several phases under ambient pressure. The stable α phase has a superconducting transition temperature, T c , of 0.9 K. By contrast, the T c of the metastable β phase is around 6 K. To understand the significant improvement in T c in the β phase, we first calculate the electronic structure, phonon dispersion, and the electron-phonon coupling of gallium in the α and β phases. Next, we solve the Eliashberg equations to obtain the superconducting gaps and the transition temperatures. Additionally, using these results, we relate the increased T c in the β phase to structural differences between the phases that affect the electronic and phonon properties. The structure motif of the α phase is Ga 2 dimers, which form strong covalent bonds leading to bonding and antibonding states that reduce the density of states at the Fermi level. The β-Ga structure consists of arrays of Ga chains that favor strong coupling between the lattice vibrations and the electronic states near the Fermi level. The increased density of states and strong coupling to the phonons for the β-Ga chains compared to the αGa 2 dimers enhance superconductivity in the β-Ga phase.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Topological invariants of a filling-enforced quantum band insulator

Traditional ionic and covalent compound insulators arise from a commensuration between electron count and system volume. On the other hand, conventional topological insulators, outside of quantum Hall effect systems, do not typically display such a commensuration. Instead, they can undergo a phase transition to a trivial insulator that preserves the electron filling. Nevertheless, in some crystalline insulators, termed filling-enforced quantum band insulators (feQBIs), electron filling can dictate nontrivial topology in the insulating ground state. Currently, little is known about the relation between feQBIs and conventional topological invariants. In this work, we study such relations for a particularly interesting example of a half-filling feQBI that is realized in space group 106 with spinless electrons. In this work, we prove that any four-band feQBI in space group 106 with filling 2 must have a nontrivial topological invariant, namely, the $\mathbb{Z}$ 2 glide invariant, and thus must have a quantized magnetoelectric polarizability θ = π. We thus have found a three-dimensional example where electron filling and band topology are tied. Such a locking raises intriguing questions about the generality of the band-inversion paradigm in describing the transition between trivial and topological phases.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Electronic structure of the frustrated diamond lattice magnet NiRh 2 O 4

We report the A -site spinel NiRh 2 O 4 is the only known realization of a spin-1 diamond lattice magnet and is predicted to host unconventional magnetic phenomena driven by frustrated nearest and next-nearest neighbor exchange as well as orbital degeneracy. Previous works found no sign of magnetic order but found a gapped dispersive magnetic excitation indicating a possible valence bond magnetic ground state. However, the presence of many competing low energy degrees of freedom and limited empirical microscopic constraints complicates further analysis. Here we carry out resonant inelastic x-ray scattering (RIXS), x-ray absorption spectroscopy (XAS), and inelastic neutron scattering (INS) to characterize the local electronic structure and lattice dynamics of NiRh 2 O 4 . The RIXS data can be partly described by a single-ion model for tetrahedrally coordinated Ni 2+ and indicates a tetragonal distortion Δt 2 = 70 meV that splits the t 2 orbitals into a high energy orbital singlet and lower energy orbital doublet. We identify features of the RIXS spectra that are consistent with a Rh-Ni two-site excitation indicating strong metal-metal hybridization mediated by oxygen in NiRh 2 O 4 . We also identify signatures of electron-phonon coupling through the appearance of phonon sidebands that dress crystal field excitations. These results establish the key energy scales relevant to the magnetism in NiRh 2 O 4 and further demonstrate that covalency and lattice dynamics play essential roles in controlling the magnetic ground states of A -site spinels.

36 MATERIALS SCIENCE↗
DOE OSTI2022

A -type antiferromagnetic order in the Zintl-phase insulator EuZn 2 P 2

Zintl phases, containing strongly covalently bonded frameworks with separate ionically bonded ions, have emerged as a critical materials family in which to couple magnetism and strong spin-orbit coupling to drive diverse topological phases of matter. Here we report the single-crystal synthesis, magnetic, thermodynamic, transport, and theoretical properties of the Zintl compound EuZn 2 P 2 that crystallizes in the anti-La 2 O 3 (CaAl 2 Si 2 ) P-3m1 structure, containing triangular layers of Eu 2+ ions. In-plane resistivity measurements reveal insulating behavior with an estimated activation energy of E g = 0.11eV. Specific heat and magnetization measurements indicate antiferromagnetic ordering at T N = 23K. Curie-Weiss analysis of in-plane and out of plane magnetic susceptibility from T = 150 to 300 K yields p eff = 8.61 for μ 0 H⊥c and p eff = 7.74 for μ 0 H//c, close to the expected values for the 4f 7 J = S = 7/2 Eu 2+ ion and indicative of weak anisotropy. Below T N , a significant anisotropy of χ ⊥ /χ // ≈ 2.3 develops, consistent with A-type magnetic order as observed in isostructural analogs and as predicted by the density functional theory calculations reported herein. The positive Weiss temperatures of θ W =19.2K for μ 0 H⊥c and θ W =41.9K for μ 0 H//c show a similar anisotropy and suggest competing ferromagnetic and antiferromagnetic interactions. Comparing Eu magnetic ordering temperatures across trigonal EuM 2 X 2 (M= divalent metal, X= pnictide) shows that EuZn 2 P 2 exhibits the highest ordering temperature, with variations in T N correlating with changes in expected dipolar interaction strengths within and between layers and independent of the magnitude of electrical conductivity. These results provide experimental validation of the crystochemical intuition that the cation Eu 2+ layers and the anionic (M 2 X 2 ) 2– framework can be treated as electronically distinct subunits, enabling further predictive materials design.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Lattice vibrations and energy landscape of the isoelectronic semiconductor series CuBr, ZnSe, GaAs, and Ge: The special case of CuBr and its $\mathcal{d}$-level chemistry

Here we have examined the lattice vibrations and the energy landscape of the isoelectronic diamond and zincblende semiconductor series CuBr, ZnSe, GaAs, and Ge. Vibrations are found to be an increasing function of ionicity, with the cation sublattice always vibrating more strongly than the anion sublattice. These findings are consistent with density functional theory (DFT) calculations of the energy landscape and temperature-dependent molecular dynamics simulations of the atomic-position fluctuations. For CuBr, inclusion of the Cu 3d Hubbard U term is found necessary to stabilize the zincblende structure and to bring its vibrational amplitudes into agreement with experiment. In addition, vibrations are found to strongly affect the CuBr near-edge x-ray absorption fine structure that we have successfully modeled by including displacements in our theoretical Bethe-Salpeter equation calculations. Reverse Monte Carlo structural refinements using large atomic configurations to simultaneously fit x-ray absorption and x-ray total-scattering data support these conclusions, and they reveal strong Cu-Br first-neighbor correlations and asymmetric distributions of interatomic distances in the temperature ranges of both negative and positive thermal expansion. Delineation of the CuBr valence band photoelectron spectrum into its Cu 3d and Br 4p states uniquely reveals their covalent mixing and further supports the DFT results.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Breaking up CO 2 through pressure-induced redox reaction with rhenium

We present results from density-functional-theory (DFT) based molecular dynamics simulations of fluid mixtures of carbon dioxide (CO 2 ) and rhenium (Re) metal at pressures up to 45 GPa and a temperature of 4000 K. Covalently-bonded carbon chains form as a result of a redox reaction between CO 2 and Re, whereas a pure CO 2 fluid is stable against decomposition at the same conditions. Furthermore, these findings have implications in reducing CO 2 by pressure-induced redox reactions with transition metals, which may play a role in CO 2 destruction and carbon storage as related to possible solutions to the current climate crisis.

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

Temperature-induced gradual polymorphic structural crossover in liquid indium, tin, and antimony

Temperature-induced polymorphic structural crossover in metallic liquid indium (In), tin (Sn), and antimony (Sb) is investigated by in situ high-energy x-ray diffraction and ab initio molecular dynamics simulations. The results demonstrate the existence of a temperature-induced reversible and gradual structural crossover in these three liquids, reflecting from both “static” structures, e.g., peak positions of the structure factor and pair distribution function, bond angle distribution, number of tetrahedra, cluster connection, free volume, and dynamical behaviors, including the relaxation time, diffusion coefficient, and non-Gaussian parameter. Further, the main difference from liquid In, Sn, to Sb is that not only the large-sized free volume but also the anharmonicity and the dynamical heterogeneity increase significantly with the appearance of strong covalent bonds. These findings will deepen the understanding of liquid structure and the liquid-liquid transition in metallic liquids.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Viscosity bounds in liquids with different structure and bonding types

Recently, it was realized that liquid viscosity has a lower bound which is nearly constant for all liquids and is governed by fundamental physical constants. This was supported by experimental data in noble and molecular liquids. Here, we perform large-scale molecular dynamics simulations to ascertain this bound in two other important liquid types: the ionic molten salt system LiF and metallic Pb. We find that these ionic and metallic systems similarly have lower viscosity bounds corresponding to the minimum of kinematic viscosity of ∼ 10 − 7 m 2 / s . We show that this agrees with experimental data in other systems with different structures and bonding types, including noble, molecular, metallic, and covalent liquids. This expands the universality of viscosity bounds into the main system types known. Published by the American Physical Society 2024

Materials Science↗
DOE OSTI2024

Unraveling the magnetic ground state in the alkali-metal lanthanide oxide Na 2 ⁢Pr⁡O 3

Here, a comprehensive set of muon spin spectroscopy and neutron scattering measurements supported by ab initio and model Hamiltonian simulations have been used to investigate the magnetic ground state of Na 2 ⁢PrO 3 . μSR reveals a Néel antiferromagnetic order below T N ~ 4.9K, with a small static magnetic moment m static ≤ 0.22 μ B /Pr collinearly aligned along the c axis. Inelastic neutron measurements reveal the full spectrum of crystal field excitations and confirm that the Pr 4+ ground-state wave function deviates significantly from the Γ 7 limit that is relevant to the Kitaev model. Single- and two-magnon excitations are observed in the ordered state below T N =4.6K and are well described by nonlinear spin wave theory from the Néel state using a magnetic Hamiltonian with Heisenberg exchange J=1 meV and symmetric anisotropic exchange Γ/J=0.1, corresponding to an XY model. Intense two magnon excitations are accounted for by g-factor anisotropy g z /g ± = 1.29. A fluctuating moment δm 2 = 0.57 ⁢(22) ⁢μ$^2_B$/Pr extracted from the energy and momentum integrated inelastic neutron signal is reduced from expectations for a local J = 1/2 moment with average g factor g avg ≈ 1.1. Together, the results demonstrate that the small moment in Na 2 ⁢PrO 3 arises from crystal field and covalency effects and the material does not exhibit significant quantum fluctuations.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Phase Diagram of a Deep Potential Water Model

Using the Deep Potential methodology, we construct a model that reproduces accurately the potential energy surface of the SCAN approximation of density functional theory for water, from low temperature and pressure to about 2400 K and 50 GPa, excluding the vapor stability region. The computational efficiency of the model makes it possible to predict its phase diagram using molecular dynamics. Satisfactory overall agreement with experimental results is obtained. Here, the fluid phases, molecular and ionic, and all the stable ice polymorphs, ordered and disordered, are predicted correctly, with the exception of ice III and XV that are stable in experiments, but metastable in the model. The evolution of the atomic dynamics upon heating, as ice VII transforms first into ice VII" and then into an ionic fluid, reveals that molecular dissociation and breaking of the ice rules coexist with strong covalent fluctuations, explaining why only partial ionization was inferred in experiments.

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

Molecule-surface interaction from van der Waals-corrected semilocal density functionals: The example of thiophene on transition-metal surfaces

Semilocal density functional approximations are widely used. None of them can capture the long-range van der Waals (vdW) attraction between separated subsystems, but they differ remarkably in the extent to which they capture intermediate-range vdW effects responsible for equilibrium bonds between neighboring small closed-shell subsystems. The local density approximation (LDA) often overestimates this effect, while the Perdew-Burke-Ernzerhof (PBE) generalized gradient approximation (GGA) underestimates it. The strongly constrained and appropriately normed (SCAN) meta-GGA often estimates it well. All of these semilocal functionals require an additive nonlocal correction such as the revised Vydrov-Van Voorhis 2010 (rVV10) to capture the long-range part. This work reports adsorption energies and the corresponding geometry of aromatic thiophene (C 4 H 4 S) bound to transition-metal surfaces. The adsorption process requires a genuine interplay of covalent and weak binding and requires a simultaneously accurate description of surface and adsorption energies with the correct prediction of the adsorption site. All these quantities must come from well-balanced short- and long-range correlation effects for a universally applicable method for weak interactions with chemical accuracy. Here, our methods indicate that the correct interplay is not present in any combination of recent meta-GGA's and rVV10. The simple short-range damping of the vdW correction scheme that is practically successful in the combination of GGA's and vdW approximations is less transferable in SCAN+rVV10 or in the revised version, revSCAN+rVV10. In addition, we present accurate random-phase-approximation-quality adsorption energies from a model based on the one of Zaremba and Kohn.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Different bonding type along each crystallographic axis: Computational study of poly($\textit{p}$-phenylene terephthalamide)

Poly($\textit{p}$-phenylene terephthalamide) (PPTA) exhibits van der Waals (vdW) bonding along its $\textit{a}$ axis, hydrogen bonding along its $\textit{b}$ axis, and covalent bonding along its $\textit{c}$ axis. In this work, we explore the structural and mechanical properties of PPTA using density functional theory with various functionals including LDA, PBE, PBE+rVV10L, SCAN, and SCAN+rVV10, compared with available experiments. The hierarchy of nonempirical semilocal functionals LDA, PBE, and SCAN (not fitted to any multicenter bonded system) includes differing amounts of intermediate-range vdW interaction. rVV10 is the long-range vdW correction. (rVV10L differs from rVV10 only in the value of a range parameter.) Among the tested functionals, SCAN shows the best performance for the lattice parameters of PPTA along the two crystal directions involving vdW or hydrogen-bond interaction. The equilibrium lattice constants obtained by SCAN and PBE+rVV10L are closest to experimental data, while SCAN+rVV10 slightly overbinds the system. We study the mechanical response of PPTA by applying strain along three lattice directions. Due to the inclusion of vdW interaction, SCAN, PBE+rVV10L, and SCAN+rVV10 all exhibit correct bonding strain-energy curves. On the contrary, PBE strongly underestimates the vdW interaction needed to resist uniaxial stretching along the $\textit{a}$ axis. The Young's modulus and yield strength of PPTA are computed and compared with previous results. The experimental values are much smaller than the computed ones, mainly due to the fact that the PPTA fiber samples used for measurements are mechanically weaker than the perfect molecular crystals considered in the simulations. Interestingly, when a compressive uniaxial stress of 25 Kbar is applied along the $\textit{b}$ axis, a structural phase transition, in which the hydrogen bonds reform along one diagonal of the $\textit{ab}$ rectangle, is predicted by SCAN.

36 MATERIALS SCIENCE↗