Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ab initio data”

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

Lattice dynamics and thermodynamics for δ -plutonium from density functional theory

Here, we present results from density functional theory (DFT) calculations of the lattice dynamics (phonons) and thermodynamics for δ-phase plutonium. The fully relativistic electronic structure is calculated assuming a three-dimensional noncollinear magnetic structure in conjunction with DFT and the general gradient approximation for the electron exchange and correlation interactions. The electronic-structure model is further enhanced by addressing strong orbital-orbital coupling via the conventional orbital-polarization (OP) scheme as has been successfully done for plutonium. The temperature dependence of the phonons is calculated within the self-consistent ab initio lattice dynamics approach. The obtained phonons compare very well with measurements although a modest overestimation of the transverse L-point [ξξξ] phonon is acknowledged. Calculated thermal vibration amplitudes and the associated Debye-Waller temperatures are close to experiments. Lattice, electronic, and magnetic contributions to the heat capacity are predicted and consistent to a few percent with that deduced from experimental data. Good agreement is only achieved when a magnetic contribution to the specific heat is recognized. The parameter-free DFT+OP electronic model is thus capable of predicting phonon properties and thermodynamic behavior of δ-phase plutonium rather accurately.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2018

Band and Impurity States in Dimond with the (MV){sup –} (M = Si, Ge, Sn) Centers Based on ab Initio Calculations

Ab initio calculations were used to determine the structure, electronic zone dispersion, partial (atomic and orbital) densities of states, and levels of impurity–vacancy cluster (MV){sup –} (M = Si, Ge, Sn) in diamond. Wavefunctions were determined in the density functional theory approximation using hybrid exchange-correlation functional, localized basis set of Gaussians, and a supercell consisting of 192 atoms. Results are compared with data on undoped diamond and with experimental results. The dependence of the direct and indirect energy gap on the cluster type is discussed. Calculations are carried out taking into account the spin–orbit interaction to identify impurity levels in the valence band. The effect of electron–electron interaction on the position of impurity levels is discussed. The dominant contribution of the d orbitals to the density of states at the impurity level in all clusters is detected.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2019

Structure Formation of Hexagonal Diamond: Ab Initio Calculations

The formation of structure of hexagonal diamond from graphite and cubic diamond is simulated with the density functional theory. Orthorhombic AB graphite transforms into hexagonal diamond under uniaxial compression at pressures above 61 GPa. The pressure should be higher than 57 GPa to form hexagonal diamond from AA graphite. Another possible way to obtain hexagonal diamond is uniaxial compression of cubic diamond under pressures from 300 to 380 GPa. The theoretical calculations were used to interpret the diffraction data obtained after X-ray diffraction and electron microscopic analysis of natural and synthetic diamond-like materials. We found that there are no crystals with ideal cubic and hexagonal structures in diamond-like carbon materials, and their structure is characterized by random packing of molecular layers.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Theoretical studies of chemical reactions related to the formation and growth of polycyclic aromatic hydrocarbons (PAH) and molecular properties of their key intermediates (Final Progress Report)

The formation mechanisms of polycyclic aromatic hydrocarbons, (PAHs) – organic molecules carrying fused benzene rings – are of great interest to scientists and engineers due to their importance in combustion chemistry and astrochemistry. On Earth, PAHs are largely produced in incomplete combustion of fossil fuel and are considered as critical precursors to unwanted soot particles leading to combustion inefficiency and causing air pollution along with detrimental health effects. Simple PAH molecules initially formed in the gas phase, are further involved in a build-up process in combustion flames leading to larger PAH, bowl-shaped nanostructures, fullerenes, and solid-phase species including carbonaceous dust, graphene particles, and soot. In deep space, PAH and their derivatives are potential key intermediates and nucleation sites leading eventually to carbonaceous nanoparticles (“interstellar grains”). Therefore, the understanding of the key processes in the synthesis of PAHs along with their precursors and their degradation mechanisms in combustion systems and in interstellar, circumstellar, and planetary atmospheric environments will provide critical insights into how complex aromatic structures, carbonaceous nanoparticles, and fullerenes are formed and destroyed. Achieving this understanding is an important step in the development of the efficient combustion processes and of the ecofriendly devices with reduced environmental pollution as well as technological strategies for the production of hydrogen and solid carbon through thermal or plasma-assisted pyrolysis of natural gas and biomass. Also, the understanding of the key processes of PAH and soot growth will help in our comprehension of chemical evolution in the universe. Detailed information on the mechanisms and reliable rate constants of the key elementary chemical reactions involved in PAH formation and destruction processes and in inception of soot particles is often missing, with the main deficiencies being the absence of temperature- and pressure-dependent rate constants for the broad range of conditions occurring in various terrestrial and interstellar processes and the lack of data on the reaction products and their branching ratios. Complementary to experimental studies, these gaps in knowledge can be filled by using quantum chemical calculations of reaction potential energy surfaces providing us with accurate energies of reaction products, intermediates, and transition states, revealing the reaction mechanism, and giving the molecular properties required to compute rate constants for relevant reaction steps and product branching ratios using the RRKM-Master Equation (ME) method. Molecular dynamics (MD) simulations can be used in cases when a reaction rate cannot be properly described by statistical theories. During the terminal renewal project period we employed these ab initio/RRKM-ME and MD approaches to complete our studies on several key reactions relevant to the formation/growth of PAH and inception of soot particles including (1) the reaction mechanism and kinetics of the resonance stabilized fulvenallenyl radical with propargyl and C 3 H 4 isomers; (2) the reaction mechanism and kinetics for the C + indene and C 2 + styrene reactions producing naphthyl or azulenyl radicals in low-temperature environments; (3) the MD study of non-equilibrium dimerization of acepyrene and coronene and its radical. The information derived from our theoretical calculations contributed to a better fundamental understanding of the reaction mechanisms and provide missing critical kinetic data to improve combustion models of hydrocarbon fuels and astrochemical models of the growth of carbonaceous molecules and particles in cold molecular clouds, circumstellar envelopes, and planetary atmospheres.

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

The roles of three-nucleon force and continuum coupling in mirror symmetry breaking of oxygen mass region

With both three-nucleon force and continuum coupling included, we have developed a self-consistent ab initio Gamow shell model within the Gamow Hartree-Fock (GHF) basis obtained by the realistic interaction itself. With the chiral two-nucleon N 3 LO and three-nucleon N 2 LO interactions, the Gamow shell model has been applied to the mirror systems of Z = 8 neutron-rich isotopes and N = 8 proton-rich isotones, giving good agreements with data in binding energies, dripline positions and excitation spectra. The GHF calculated that the 0 ds/2 , 1 S1/2 and 1 p3/2 orbitals are resonances. The resonance states and their interplay with nonresonant continua play a crucial role in the descriptions of nuclei around driplines. Excitation spectra and Thomas-Ehrman shifts observed can be better described when both three-nucleon force and continuum coupling are considered in calculations. The three-nucleon force and continuum coupling produce a combined effect on the Thomas-Ehrman shift, e.g., for the 1/2 + resonance level of 19Na. The calculations help the understandings of related nuclear astrophysical processes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Probing the antiknock effect of anisole through an ignition, speciation and modeling study of its blends with isooctane

In order to unravel the reaction pathways relevant to anisole co-oxidation within a fuel blend, a detailed study of isooctane/anisole blends was performed with the ULille RCM. Ignition delays as well as mole fraction profiles were measured during a two-stage ignition delay using sampling and GC techniques. These results are used to validate a kinetic model developed from ab initio calculations for the most relevant rate constants which included H-atom abstraction reactions from anisole, and reactions on the potential energy surfaces of methoxyphenyl + O 2 and anisyl + O 2 . Pressure dependent rate constants were computed for the methoxyphenyl + O 2 and anisyl + O 2 reactive systems using master equation code analysis. The new kinetic model shows good agreement with the experimental data. Dual brute-force sensitivity analysis was performed, on both first- and second-stages of ignition, allowing the identification of the most important reactions in the prediction of both ignition delays. It was observed that while pure anisole does not show NTC behavior, a 60/40 isooctane/anisole blend displays such behavior, as well as two-stage ignition. This suggests anisole addition may not be as beneficial to knock resistance as expected from its high octane number. The kinetic modeling results demonstrate the importance of H-abstraction reactions both from the methoxy group and from the aryl ring in ortho-position and the addition of the resultant radicals to O2, mostly leading to the formation of polar or non-aromatic products.

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

Learning Grain-Boundary Segregation: From First Principles to Polycrystals

The segregation of solute atoms at grain boundaries (GBs) can strongly impact the structural and functional properties of polycrystals. Yet, due to the limited availability of simulation tools to study polycrystals at the atomistic scale (i.e., interatomic potentials), there is a minimal understanding of the variation of solute segregation tendencies across the very complex space of GB microenvironments and the large range of alloys in which it can occur. Here, we develop an algorithmic framework that can directly learn the full spectrum of segregation energies for a metal solute atom in a metal polycrystal from ab initio methods, bypassing the need for alloy interatomic potentials. This framework offers a pathway to a comprehensive catalog of GB solute segregation with quantum accuracy, for the entire alloy space. As an initial demonstration in this pursuit, we build an extensive GB segregation database for aluminum-based alloys across the periodic table, including dozens of alloys for which there are substantially no prior data.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1981

Theoretical Treatment of the X(sup 1)Sigma(sup +), A(sup 1)Sigma(sup +), and B(sup 1)Pi states of LiH

Ab initio self-consistent-field plus configuration-interaction calculations are reported for the X(sup 1)Sigma(sup +), A(sup 1)Sigma(sup +), and B(sup 1)Pi states of LiH using a 22(sugma)12(pi)7(delta) function Slater basis set. The resulting dissociation energies, with the experimental values in parentheses, are D(sub e)(X(sup 1)Sigma(sup +)) = 19 972 (20288)/ cm, and D(sub e)(A(sup 1)Sigma(sup +)) = 9042 (8682)/ cm, and D(sup e)(B(sup 1)Pi) = 284 (288) /cm. This is the first ab initio treatment to quantitatively account for the binding in the B(sup 1)Pi state. Calculated dipole moments and electronic transition moments for the X(sup 1)Sigma(sup+)-A(sup 1)Sigma(sup +), X(sup 1)Sigma(sup +)-B((sup 1)Pi, and A(sup 1)Sigma(sup +)- B(sup 1)PI band systems are in excellent agreement with existing theoretical and experimental data. Radiative transition probabilities and lifetimes, including both the bound-bound and bound-free contributions, are computed for all vibrational levels of the A(sup 1)Sigma(sup+) and B(sup 1)Pi states. Comparison with previous results using experimentally based potentials provides insight into the sensitivity of the radiative lifetimes to the detailed nature of the uppermost region of the potentials. Our calculated lifetimes for the lower vibrational levels of the A(sup 1)Sigma(sup +) state are within the experimental uncertainty. Our calculated lifetimes for the three vibrational levels of the B(sup 1)Pi state are in excellent agreement with those of Zemke and Stwalley (values in parentheses), increasing with (upsilon)' from 11.3 (11.3) nsec at (upsilon)' = 0, to 17.0 (17.0) nsec at (upsilon)' = 1, and then to 23.5 (24.0) nsec at (upsilon)' = 2.

Partridge, Harry↗
DOE OSTI2018

Large-scale exact diagonalizations reveal low-momentum scales of nuclei

Ab initio methods aim to solve the nuclear many-body problem with controlled approximations. Virtually exact numerical solutions for realistic interactions can only be obtained for certain special cases such as few-nucleon systems. In this paper, we extend the reach of exact diagonalization methods to handle model spaces with dimension exceeding ${10}^{10}$ on a single compute node. This allows us to perform no-core shell model (NCSM) calculations for $^{6}\mathrm{Li}$ in model spaces up to ${N}_{\mathrm{max}}=22$ and to reveal the $^{4}\mathrm{He}$+d halo structure of this nucleus. Still, the use of a finite harmonic-oscillator basis implies truncations in both infrared (IR) and ultraviolet (UV) length scales. These truncations impose finite-size corrections on observables computed in this basis. We perform IR extrapolations of energies and radii computed in the NCSM and with the coupled-cluster method at several fixed UV cutoffs. It is shown that this strategy enables information gain also from data that is not fully UV converged. IR extrapolations improve the accuracy of relevant bound-state observables for a range of UV cutoffs, thus making them profitable tools. We relate the momentum scale that governs the exponential IR convergence to the threshold energy for the first open decay channel. Finally, using large-scale NCSM calculations we numerically verify this small-momentum scale of finite nuclei.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Calibration and validation of the foundation for a multiphase strength model for tin

In this work, the Common Model of Multi-phase Strength and Equation of State (CMMP) model was applied to tin. Specifically, calibrations of the strength-specific elements of the CMMP foundation were developed with a combination of experiments and theory, and then the model was validated experimentally. The first element of the foundation is a multi-phase analytic treatment of the melt temperature and the shear modulus for the solid phases. These models were parameterized for each phase based on ab initio calculations using the software VASP (Vienna Ab initio Simulations Package) based on density functional theory. The shear modulus model for the ambient phase was validated with ultrasonic sound speed measurements as a function of pressure and temperature. The second element of the foundation is a viscoplastic strength model for the phase, upon which strength for inaccessible higher-pressure phases can be scaled as necessary. The stress–strain response of tin was measured at strain rates of 10 -3 to 3 x 10 3 s -1 and temperatures ranging from 87 to 373 K. The Preston–Tonks–Wallace (PTW) strength model was fit to that data using Bayesian model calibration. For validation, six forward and two reverse Taylor impact experiments were performed at different velocities to measure large plastic deformation of tin at strain rates up to ⁠10 5 s -1 . The PTW model accurately predicted the deformed shapes of the cylinders, with modest discrepancies attributed to the inability of PTW to capture the effects of twinning and dynamic recrystallization. Some material in the simulations of higher velocity Taylor cylinders reached the melting temperature, thus testing the multiphase model because of the presence of a second phase, the liquid. In simulations using a traditional modeling approach, the abrupt reduction of strength upon melt resulted in poor predictions of the deformed shape and non-physical temperatures. With CMMP, the most deformed material points evolved gradually to a mixed solid–liquid but never a fully liquid state, never fully lost strength, stayed at the melt temperature as the latent heat of fusion was absorbed, and predicted the deformed shape well.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Precision Mass Measurements Reveal Low Neutron Pairing in Tin beyond 𝑁=82 and Its Impact on Stellar Nucleosynthesis

We present a study on neutron-rich tin (𝑍 =50) isotopes beyond the doubly closed shell of 𝑁 = 82 through high-precision mass measurements, including the first-ever measurements of the masses of 136 Sn, 137 Sn, and 138 Sn isotopes. These measurements enhance our understanding of the nuclear structure and astrophysical nucleosynthesis in this previously unexplored region. The new mass data are used for evaluation of the final abundances of mass numbers 𝐴 =135 and 137 in 𝑟-process network calculations. Our findings reveal a notable change in the empirical pairing gap for tin isotopes beyond the 𝑁 = 82 closed shell and a shift in the two-neutron-separation energy slope compared to heavier elements above the shell closure. A new set of ab initio calculations effectively describes these observed trends.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Randomized low-rank decompositions of nuclear three-body interactions

First-principles simulations of many-fermion systems are commonly limited by the computational requirements of processing large data objects. As a remedy, we propose the use of low-rank approximations of three-body interactions, which are the dominant such limitation in nuclear physics. We introduce a randomized decomposition technique to handle the excessively large matrix dimensions and study the sensitivity of low-rank properties to interaction details. The developed low-rank three-nucleon interactions are benchmarked in ab initio simulations of few- and many-body systems. Exploiting low-rank properties provides a promising route to extend the microscopic description of atomic nuclei to large systems where storage requirements exceed the computational capacities of the most advanced high-performance computing facilities.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

NbTaV_BCC_SolidSolution_128atoms_VASP6

We performed density functional theory (DFT) calculations for body-centered-cubic (BCC) structures with 128 lattices sites of solid solution ternary alloys niobium-tantalum-vanadium (Nb-Ta-V). The first-principle code that has been used to run the calculations is the closed-source Vienna Ab-Initio Simulation Package (VASP). The calculations have been collected by sampling chemical compositions across the entire compositional range. The chemical compositions have been sampled by progressively changing the number of atoms per constituent by 8. For each chemical composition of binaries and ternaries, the first-principle calculations have been run for 100 randomized arrangements of the constituents on the BCC lattice sites. We collected data for a total of 10,500 randomized atomic structures over 105 chemical compositions. The calculations have been collected on NERSC-Perlmutter using the VASP 6.3.2. The VASP calculations for every atomic structure have been performed in 2 main steps: 1. Starting from an ideal body-centered-cubic (BCC) structure, geometry optimization with low precision has been executed to perform a preliminary optimization of the atomic structure. The output for this calculations is available in the files 0.CONTCAR, 0.OUTCAR, rlx1.out. 2. Using the atomic structure resulting from the preliminary geometry optimization, a second geometry optimization has been performed using normal precision. The output for this calculations is available in the files CONTCAR, OUTCAR, rlx2.out, vaspout.h5, and vasprun.xml. Every chemical composition sampled across the composition range in the dataset has its own directory. The convention used to name the directories for ternary alloys is AXBYCZ, where A, B, and C refer to the constituents, and X, Y, and Z are positive integers that represent the number of atoms for each constituent and their values still sum up to 128. Each atomic structure associated with a specific chemical composition has its own sub-directory within the directory of the corresponding chemical composition. The sub-directories for each atomic structure for each chemical composition are named 'case-*', where * is a positive integer that spans all the values from 1 through 100, extremes included. The files contained in each sub-directory 'case-*' for each atomic structure are as follows: 1. INCAR: input file that contains various parameters and settings for controlling the behavior of the electronic structure calculations 2. 0.POSCAR: input file that defines the atomic structure of a system 3. 0.CONTCAR: output file that provides the atomic positions and cell parameters after the first geometry optimization has been run with the precision variable set to PREC=Low in the INCAR file 4. 0.OUTCAR: output file that contains detailed information about the progress of a calculation after the first geometry optimization has been run with the precision variable set to PREC=Low in the INCAR file 5. rlx1.out: file with diagnostic information about the execution of the first geometry optimization with precision variable set to PREC=Low in the INCAR file 6. POSCAR: input file that defines the atomic structure of a system after the first geometry optimization has been run at low precision. This represents the input for the second geometry optimization run with the precision variable set to PREC=Normal in the INCAR file 7. CONTCAR: output file that provides the atomic positions and cell parameters after the second geometry optimization has been run with the precision variable set to PREC=Normal in the INCAR file 8. OUTCARL: output file that contains detailed information about the progress of a calculation after the second geometry optimization has been run with the precision variable set to PREC=Normal in the INCAR file 9. rlx2.out: file with diagnostic information about the execution of the second geometry optimization with precision variable set to PREC=Normal in the INCAR file 10. vaspout.h5: hierarchical HDF5 file containing the inputs and outputs of a VASP calculation. To analyze the data in this file we recommend using py4vasp. This file is only produced if the VASP version used is compiled with HDF5 support 11. vasprun.xml: contains similar information to OUTCAR, but in an xml format. This research is sponsored by the Artificial Intelligence Initiative as part of the Laboratory Directed Research and Development (LDRD) Program of Oak Ridge National Laboratory, managed by UT-Battelle, LLC, for the US Department of Energy under contract DE-AC05-00OR22725. This work used resources of the Oak Ridge Leadership Computing Facility, which is supported by the Office of Science of the U.S. Department of Energy under Contract No. DE-AC05-00OR22725, under Directorate Discretionary awards MAT025 (Materials Science) and LRN026 (Machine Learning), and INCITE award MAT201. This work also used resources of the National Energy Research Scientific Computing Center, which is supported by the Office of Science of the U.S. Department of Energy under Contract No. DE-AC02-05CH11231, under award ERCAP0025216.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Systematics of the Third Row Transition Metal Melting: The HCP Metals Rhenium and Osmium

The melting curves of rhenium and osmium to megabar pressures are obtained from an extensive suite of ab initio quantum molecular dynamics (QMD) simulations using the Z method. In addition, for Re, we combine QMD simulations with total free energy calculations to obtain its phase diagram. Our results indicate that Re, which generally assumes a hexagonal close-packed (hcp) structure, melts from a face-centered cubic (fcc) structure in the pressure range 20–240 GPa. We conclude that the recent DAC data on Re to 50 GPa in fact encompass both the true melting curve and the low-slope hcp-fcc phase boundary above a triple point at (20 GPa, 4240 K). A linear fit to the Re diamond anvil cell (DAC) data then results in a slope that is 2.3 times smaller than that of the actual melting curve. The phase diagram of Re is topologically equivalent to that of Pt calculated by us earlier on. Regularities in the melting curves of Re, Os, and five other 3rd-row transition metals (Ta, W, Ir, Pt, Au) form the 3rd-row transition metal melting systematics. We demonstrate how this systematics can be used to estimate the currently unknown melting curve of the eighth 3rd-row transition metal Hf.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Role of Nematicity in Controlling Spin Fluctuations and Superconducting Tc in Bulk FeSe

FeSe undergoes a transition from a tetragonal to a slightly orthorhombic phase at 90 K and becomes a superconductor below 8 K. The orthorhombic phase is sometimes called a nematic phase because quantum oscillation, neutron, and other measurements detect a significant asymmetry in x and y. How nematicity affects superconductivity has recently become a matter of intense speculation. Here, we employ an advanced ab initio Green's function description of superconductivity and show that bulk tetragonal FeSe would, in principle, superconduct with almost the same critical temperature Tc as the nematic phase. The mechanism driving the observed nematicity is not yet understood. Since the present theory underestimates it, we simulate the full nematic asymmetry by artificially enhancing the orthorhombic distortion. For benchmarking, we compare theoretical spin susceptibilities against experimentally observed data over all energies and relevant momenta. When the orthorhombic distortion is adjusted to correlate with observed nematicity in spin susceptibility, the enhanced nematicity causes spectral weight redistribution in the Fe-3dxz and Fe-dyz orbitals, but it leads to at most a 10-15% increment in Tc. This is because the dxy orbital always remains the most strongly correlated and provides most of the source of the superconducting glue. Nematicity suppresses the density of states at the Fermi level; nevertheless, Tc increases, in contradiction to both BCS theory and the theory of Bose-Einstein condensation. We show how the increase is connected to the structure of the particle-particle vertex. Our results suggest that while nematicity may be an intrinsic property of bulk FeSe, it is not the primary force driving the superconducting pairing.

bulk FeSe↗
DOE OSTI2016

Anisotropic type-I superconductivity and anomalous superfluid density in OsB 2

Here, we present a microscopic study of superconductivity in OsB 2 , and discuss the origin and characteristic length scales of the superconducting state. From first-principles we show that OsB 2 is characterized by three different Fermi sheets, and we prove that this fermiology complies with recent quantum-oscillation experiments. Using the found microscopic properties, and experimental data from the literature, we employ Ginzburg-Landau relations to reveal that OsB 2 is a distinctly type-I superconductor with a very low Ginzburg-Landau parameter κ—a rare property among compound materials. We show that the found coherence length and penetration depth corroborate the measured thermodynamic critical field. Moreover, our calculation of the superconducting gap structure using anisotropic Eliashberg theory and ab initio calculated electron-phonon interaction as input reveals a single but anisotropic gap. The calculated gap spectrum is shown to give an excellent account for the unconventional behavior of the superfluid density of OsB 2 measured in experiments as a function of temperature. This reveals that gap anisotropy can explain such behavior, observed in several compounds, which was previously attributed solely to a two-gap nature of superconductivity.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Coordination Sphere of Lanthanide Aqua Ions Resolved with Ab Initio Molecular Dynamics and X-ray Absorption Spectroscopy

To resolve the fleeting structures of lanthanide Ln 3+ aqua ions in solution, we (i) performed the first ab initio molecular dynamics (AIMD) simulations of the entire series of Ln 3+ aqua ions in explicit water solvent using pseudopotentials and basis sets recently optimized for lanthanides and (ii) measured the symmetry of the hydrating waters about Ln 3+ ions (Nd 3+ , Dy 3+ , Er 3+ , Lu 3+ ) for the first time with extended X-ray absorption fine structure (EXAFS). EXAFS spectra were measured experimentally and generated from AIMD trajectories to directly compare simulation, which concurrently considers the electronic structure and the atomic dynamics in solution, with experiment. In this study, we performed a comprehensive evaluation of EXAFS multiple-scattering analysis (up to 6.5 Å) to measure Ln–O distances and angular correlations (i.e., symmetry) and elucidate the molecular geometry of the first hydration shell. This evaluation, in combination with symmetry-dependent L 3 - and L 1 -edge spectral analysis, shows that the AIMD simulations remarkably reproduces the experimental EXAFS data. The error in the predicted Ln–O distances is less than 0.07 Å for the later lanthanides, while we observed excellent agreement with predicted distances within experimental uncertainty for the early lanthanides. Our analysis revealed a dynamic, symmetrically disordered first coordination shell, which does not conform to a single molecular geometry for most lanthanides. This work sheds critical light on the highly elusive coordination geometry of the Ln 3+ aqua ions.

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

Organo‐Functionalized Lacunary Double Cubane‐Type Oxometallates: Synthesis, Structure, and Properties of [(M II Cl) 2 (V IV O) 2 {((HOCH 2 CH 2 )(H)N(CH 2 CH 2 O))(HN(CH 2 CH 2 O) 2 )} 2 ] (M=Co, Zn)

Abstract Organofunctionalized tetranuclear clusters [(M II Cl) 2 (V IV O) 2 {((HOCH 2 CH 2 )(H)N(CH 2 CH 2 O))(HN(CH 2 CH 2 O) 2 )} 2 ] (1, M=Co,2: M=Zn) containing an unprecedented oxometallacyclic {M 2 V 2 Cl 2 N 4 O 8 } (M=Co, Zn) framework have been prepared by solvothermal reactions. The new oxo‐alkoxide compounds were fully characterized by spectroscopic methods, magnetic susceptibility measurement, DFT and ab initio computational methods, and complete single‐crystal X‐ray diffraction structure analysis. The isostructural clusters are formed of edge‐sharing octahedral {VO 5 N} and trigonal bipyramidal {MO 3 NCl} units. Diethanolamine ligates the bimetallic lacunary double cubane core of1and2in an unusual two‐mode fashion, unobserved previously. In the crystalline state, the clusters of1and2are joined by hydrogen bonds to form a three‐dimensional network structure. Magnetic susceptibility data indicate weakly antiferromagnetic interactions between the vanadium centers [J iso (V IV −V IV )=−5.4(1); −3.9(2) cm −1 ], and inequivalent antiferromagnetic interactions between the cobalt and vanadium centers [J iso (V IV −Co II )=−12.6 and −7.5 cm −1 ] contained in1.

Chemistry↗