Accurate Nonempirical Range-Separated Hybrid van der Waals Density Functional for Complex Molecular Problems, Solids, and Surfaces
Not Available
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.
Not Available
The precise synthesis of high-purity materials is crucial in accelerating materials discovery. However, the lack of theoretical understanding and practical guidance poses challenges, particularly for materials with compositional and structural complexity. Here, we propose a feasible principle toward synthesizing complex inorganic solids. This principle involves the introduction of an inducer that induces crucial intermediates, which in turn guide the synthesis pathway toward the target materials through structural templating, named inducer-facilitated assembly through structural templating (i-FAST). We validate this principle with three distinct oxides: garnet Li6.5La3Zr1.5Ta0.5O12, perovskite BaCo0.8Sn0.2O3, and pyrochlore Gd1.5La0.5Zr2O7. This structural templating approach enables synthesis along predesigned pathways, forming intermediates that are thermodynamically favored for prior formation and kinetically preferred for the final product, resulting in precisely synthesizing high-purity target materials. This study not only represents a substantial advancement in comprehending the interplay between thermodynamics/kinetics and phase evolution in complex solid synthesis but also provides an effective strategy for guiding exploratory solid-state synthesis.
Halide solid electrolytes have emerged as promising candidates for solid-state batteries owing to their high oxidative stability and ionic conductivity. Among them, Li3InCl6 (LIC) has attracted significant attention. However, diffraction patterns of LIC synthesized via different methods exhibit distinct differences particularly at low-angle reflectionsindicative of underlying structural disorder. These variations are attributed to deviations from ideal crystallographic order, especially stacking faults, whose impact on structure and ion transport remains poorly understood. Here, we identify and quantify stacking faults in LIC samples prepared under different synthetic conditions. Using X-ray diffraction and time-of-flight neutron diffraction, we construct and refine stacking fault models that accurately reproduce the experimental diffraction features. LIC samples with higher degrees of stacking faults exhibit only negligible differences in ionic conductivities and activation energies. This indicates that stacking faults have a limited impact on altering the Li+ diffusion pathway along the c-axis, likely due to the high concentration of vacancies in the In layers, while Li+ diffusion remains nearly unchanged in the ab-plane. Our results account for the observed differences in diffraction patterns across samples and provide a quantitative assessment of faulting probabilities and stacking sequences. The insights gained from this study are expected to be broadly applicable to other layered halide solid electrolytes and contribute to a deeper understanding of the role of structural disorder in ion transport.
Here we present a fully-integrated lattice Boltzmann (LB) method for fluid–structure interaction (FSI) simulations that efficiently models deformable solids in complex suspensions and active systems. Our Eulerian method (LBRMT) couples finite-strain solids to the LB fluid on the same fixed computational grid with the reference map technique (RMT). An integral part of the LBRMT is a new LB boundary condition for moving deformable interfaces across different densities. With this fully Eulerian solid–fluid coupling, the LBRMT is well-suited for parallelization and simulating multi-body contact without remeshing or extra meshes. We validate its accuracy via a benchmark of a deformable solid in a lid-driven cavity, then showcase its versatility through examples of soft solids rotating and settling. The LBRMT achieves a spatial convergence rate between first-order and second-order for FSI simulations and is designed for low to intermediate Reynolds number flows with finite inertia at small Mach numbers. With simulations of complex suspensions mixing, we highlight the potential of the LBRMT for studying collective behavior in soft matter and biofluid dynamics.
Reaction of 4 equiv. of [Li(TMEDA)] 2 [1,2-S 2 C 6 H 4 ] with [ThCl 4 (DME) 2 ] or [UCl 4 (THF) 3 ] in THF results in formation of [Li(THF) 2 ] 4 [An(1,2-S 2 C 6 H 4 ) 4 ] (An = Th, 1; An = U, 2), whereas reaction of 4 equiv. of [Li(TMEDA)] 2 [1,2-S 2 C 6 H 4 ] with UCl 4 in Et 2 O results in formation of [Li(TMEDA)] 4 [U(1,2-S 2 C 6 H 4 ) 4 ] (3). Complexes 1–3 represent the first reported benzenedithiolate complexes of the actinides. Here, they were characterized by NMR spectroscopy and X-ray crystallography. In the solid state, complexes 1–3 exhibit triangular dodecahedral geometries about their actinide centers. Additionally, their Li + cations are bound by two sulfur atoms of adjacent [1,2-S 2 C 6 H 4 ] 2− ligands, in addition to two solvent donor atoms. In solution, complexes 2 and 3 exhibit spectral data consistent with S 4 symmetry (and non-exchanging Li + sites), whereas complex 1 exhibits spectral properties consistent with labile Li + cations.
Crystalline solids assembled from superatomic building blocks are attractive functional materials due to their hierarchical structure, multifunctionality, and tunability. An interesting example is Re 6 Se 8 Cl 2 , in which the Re 6 Se 8 building blocks are covalently linked into two-dimensional (2D) sheets that are stacked into a layered van der Waals solid. It is an indirect gap semiconductor that, when heavily doped, becomes a superconductor at low temperatures. Given the finite electronic bandwidths (300–400 meV), carrier properties in this material are expected to be strongly influenced by coupling to phonons. Here, we apply angle-resolved photoemission spectroscopy to probe the valence band edge (VBE) of Re 6 Se 8 Cl 2 . We find that dispersion of the VBE is a strong function of temperature. The bandwidth is W = 120 ± 30 meV at 70 K and decreases by one order of magnitude to W ~ 10 ± 20 meV as temperature is increased to 300 K. This observation reveals the dominant polaronic effects in Re 6 Se 8 Cl 2 , consistent with the Holstein polaron model commonly used to describe molecular solids.
Abstract A four‐coordinate cobalt–imido complex, ( t Bu mPNP)Co=NMes ( t Bu mPNP=modified PNP pincer ligand) has been synthesized from addition of 2,4,6‐trimethylphenylazide (Mes–N 3 ) to the corresponding dinitrogen complex. The solid‐state structure determined by X‐ray diffraction established a rare, idealized planar geometry with a Co=N bond distance of 1.716(2) Å. Magnetic measurements revealed an S =1 ground state with CAS‐SCF calculations supporting radical character on the imide nitrogen. Thermolysis of the cobalt–imido compound induced selective insertion of the imido group into a Co−P bond and yielded a three‐coordinate cobalt complex with a distorted T‐shaped geometry. Transition state analysis conducted with DFT calculations established the thermodynamic stability of the P–N coupled product and provided insight into the exclusive selectivity.
As described here, interstitial electron density ρ o is offered as a direct metric for maximum strength in metals, arising from universal properties derived from an electron-gas – ρ o sets the exchange-correlation parameter r s in density-functional theory (DFT). It holds also for maximum shear strength τ max in polycrystals [Phys. Rev. Lett. 124, 125501]. Elastic moduli and τ max for polycrystalline (amorphous) metals are linear with ρ o and melting T m (glass-transition T g ) temperature. ρ o or r s , even with rule-of-mixture estimate, predicts relative strength for rapid, reliable selection of high-strength alloys with ductility, as confirmed for elements to steels to complex solid-solutions, and validated experimentally.
The fundamental principles that govern antiferroelectric (AFE)–ferroelectric (FE) transitions are not well understood for many solid solutions of perovskite compounds. For example, crystal chemical considerations based on the average Goldschmidt tolerance factor or ionic polarizability do not precisely predict the boundary between the AFE and FE phases in dilute solid solutions of alkali niobates, such as K x Na 1– x NbO 3 ( x ≤ 0.02). Here, based on detailed structural analysis from neutron total scattering experiments, we provide insights about how the relative local distortions around the A- and B-sites of the ABO 3 perovskite structure affect the AFE/FE order of the average crystallographic phases in K x Na 1– x NbO 3 . Specifically, we show that a higher (lower) ratio of B-site-centered distortions over A-site-centered distortions drives transition toward a long-range FE (AFE) phase, which is based on a competition between the long-range polarizing field of the Nb–O dipoles and the disordering effect of local distortions around the A-site. Our study provides a predictive tool for designing complex solid-solution perovskites with tunable (anti)ferroelectric polarization properties, which can be of interest for various energy-related applications such as high-density energy storage and solid-state cooling.
Complex-oxide superlattices provide a pathway to numerous emergent phenomena because of the juxtaposition of disparate properties and the strong interfacial interactions in these unit-cell-precise structures. This is particularly true in superlattices of ferroelectric and dielectric materials, wherein new forms of ferroelectricity, exotic dipolar textures, and distinctive domain structures can be produced. Here, in this study, relaxor-like behavior, typically associated with the chemical inhomogeneity and complexity of solid solutions, is observed in (BaTiO 3 ) n /(SrTiO 3 ) n (n = 4–20 unit cells) superlattices. Dielectric studies and subsequent Vogel–Fulcher analysis show significant frequency dispersion of the dielectric maximum across a range of periodicities, with enhanced dielectric constant and more robust relaxor behavior for smaller period n. Bond-valence molecular-dynamics simulations predict the relaxor-like behavior observed experimentally, and interpretations of the polar patterns via 2D discrete-wavelet transforms in shorter-period superlattices suggest that the relaxor behavior arises from shape variations of the dipolar configurations, in contrast to frozen antipolar stripe domains in longer-period superlattices (n = 16). Moreover, the size and shape of the dipolar configurations are tuned by superlattice periodicity, thus providing a definitive design strategy to use superlattice layering to create relaxor-like behavior which may expand the ability to control desired properties in these complex systems.
Lithium-sulfur (Li-S) batteries involve complex solid-liquid-solid phase transformations during both discharging and charging processes, where cathode materials, formulation, and structure play a crucial role. Here, a design of experiments (DoE) methodology and an empirical model are developed to systematically explore the interactions and trade-offs among cathode factors and process variables, and to obtain generalizable effects estimates for the multivariate system. Compared to the conventional one-factor-at-a-time (OFAT) approach, this work demonstrates advantages in both efficiency and accuracy by allowing the data to guide future research and decisions. Further, an optimized cathode formulation and processing parameters are predicted and validated experimentally, achieving over 1000 mAh g -1 in discharge capacity and improved cycling under practical lean electrolyte (4 µL mg -1 S) and high S-loading cathodes (>4 mg cm -2 ) conditions. The optimized cathode was scaled up and assembled into Li-S pouch cells, achieving 316 Wh kg -1 in cell-level energy, proving that the comprehensive and rigorous framework for optimizing complex systems with DoE leads to improved performance in a practical pouch cell system.
High intensity, short-pulse laser interaction with a solid metal target produces broadband hard x-rays potentially for various applications of x-ray radiography. In this work, experimental benchmarking of numerical modelling for short-pulse laser-driven broadband x-ray radiography is presented. Angular dependent x-ray spectra are first calculated with a hybrid particle-in-cell code, Large Scale Plasma (LSP), using fast electron parameters inferred from an analysis of measured bremsstrahlung signals. Subsequently, a calculated x-ray spectrum in the direction of radiography is used in photon transport calculations using a Monte Carlo code, Particle and Heavy Ion Transport code System (PHITS), to simulate a radiographic image including a modelled 3D test object, an x-ray attenuation filter and an image plate detector. Simulated radiographic images are compared with measurements obtained in an experiment using a 50-TW Leopard short-pulse laser at the University of Nevada Reno. Results show that simulations reproduce the experimental images well for three different attenuation filters (plastic, aluminium, and brass), while 1D transmission profiles for the plastic and aluminium filters are quantitatively in good agreement. The modelling approach established in this work could be used as a predictive tool to simulate radiographic images of complex 3D solid objects at any arbitrary angular position or to optimize experimental components such as the source spectrum, x-ray attenuation filters and a detector type depending on a radiographic object without carrying out radiographic experiments.
Two type-I clathrates were synthesized by introducing Ga into the framework of the Ba 8 Cu 16 P 30 type-I clathrate. The introduction of minute amounts of Ga, 1.9% Ga/M total (where M total = Cu + Ga), resulted in the disturbance of the completely ordered Pbcn superstructure of Ba 8 Cu 16 P 30 . Ba 8 Cu 15.43(2) Ga 0.3 P 30.26(3) crystallizes in a partially ordered orthorhombic Pmna clathrate-I superstructure with five out of 15 framework sites being jointly occupied by metal+phosphorus. Increasing the Ga content resulted in all framework sites being occupied by metal + phosphorus in the archetype cubic Pm n clathrate-I crystal structure of Ba 8 Cu 14.5(3) Ga 1.3 P 30.2(4) with 8.2% Ga/M total . A combination of energy dispersive x-ray spectroscopy, inductively coupled plasma mass spectroscopy, and single crystal x-ray diffraction was used to determine the structures alongside the compositions. The positional disorder was verified by 31 P solid state NMR spectroscopy. Characterization of the transport properties indicated that the Ga-substituted samples exhibit higher Seebeck coefficients and electrical resistivities compared to its pristine counterpart, in line with the expected reduction of the hole concentration due to Ga/Cu substitution. Moderate improvements in the thermoelectric power factor and overall figure-of-merit were observed for samples with 6.9% and 3.8% Ga/M total as compared to those for the pristine Ba 8 Cu 16 P 30 clathrate. Finally, band structure calculations shed light on how Ga substitution affects the electronic structure and thermoelectric properties of studied clathrates.
A new algorithm based on a rigorous theorem and quantum data computationally mined from element 118 guarantees automated construction of initial Fermi–Löwdin-Orbital (FLO) starting points for all elements in the Periodic Table. It defines a means for constructing a small library of scalable FLOs for universal use in molecular and solid-state calculations. The method can be systematically improved for greater efficiency and for applications to excited states such as x-ray excitations and optically silent excitations. FLOs were introduced to recast the Perdew–Zunger self-interaction correction (PZSIC) into an explicit unitarily invariant form. The FLOs are generated from a set of N quasi-classical electron positions, referred to as Fermi-Orbital descriptors (FODs), and a set of N-orthonormal single-electron orbitals. FOD positions, when optimized, minimize the PZSIC total energy. However, creating sets of starting FODs that lead to a positive definite Fermi orbital overlap matrix has proven to be challenging for systems composed of open-shell atoms and ions. The proof herein guarantees the existence of a FLOSIC solution and further guarantees that if a solution for N electrons is found, it can be used to generate a minimum of N – 1 and a maximum of 2 N – 2 initial starting points for systems composed of a smaller number of electrons. As a result, applications to heavy and super-heavy atoms are presented. All starting solutions reported here were obtained from a solution for element 118, Oganesson.
High-temperature disordered multi-component alloys, including high-entropy alloys, experience either segregation or partially-ordered phases to reach low-temperature phases. For Ti 0.25 CrFeNiAl x (0 ≤ x ≤ 1), experiments suggest a partially-ordered B2 phase, whereas CALculation of PHAse Diagrams (CALPHAD) predicts a region of L21+B2 coexistence. We employ first-principles KKR-CPA to assess stability of phases with arbitrary order and a KKR-CPA linear-response theory to predict atomic short-range order (SRO) in the disorder phase that reveals the competing long-range ordered (LRO) phases in a given Bravais lattice. The favorable SRO provides a specific concentration-waves (site occupation probabilities and partially-ordered unit cells) and estimated energy gains that can then be assess directly by KKR-CPA formation enthalpies. Our results are in good agreement with experiments and CALPHAD in Al-poor regions (x ≤ 0.75) and with CALPHAD in Al-rich region (1≥ x >0.75). Our first-principles KKR-CPA and SRO-based concentration-wave analysis is shown to be a powerful and fast method to assess competing LRO phases in complex solid-solution alloys, and our results suggests more careful experiments in Al-rich region are needed.
This paper introduces a new microstructural template for high entropy alloys (HEAs), where the face centered cubic (FCC) complex concentrated solid solution is reinforced with a high density of strong, yet deformable, nanorods of an ordered multi-component intermetallic L12 compound. Thermodynamic modeling has been employed to design this HEA with a large L12 volume fraction. Thermo-mechanical processing by isothermal annealing of the conventionally processed bulk cold-rolled alloy directly at precipitation temperatures, has been applied to produce a high density of uniformly distributed L12 nanorods within refined FCC grains, resulting from concomitant recrystallization and discontinuous precipitation processes. The nanorod morphology of the discontinuous L12 product has been established from three-dimensional atom probe tomography. The refined grains result in a complete coverage of the microstructure with discontinuously precipitated intermetallic nanorods. This nanorod strengthened HEA exhibits an exceptionally high room temperature yield strength of ~1630 MPa, good tensile ductility of ~15%, and an ultimate tensile strength of ~1720 MPa. Furthermore, a single L12 phase alloy, melted based on the precipitate composition in the two-phase FCC + L12 HEA, exhibits very high compressive deformability and strain hardenability, unusual for ordered intermetallic compounds. These results open a new strategy for design of fine-grained microstructures strengthened via ordered intermetallic phases, exploiting the beneficial effects of discontinuous precipitation, for achieving very high room temperature tensile strengths while maintaining good ductility.
Here, a high-order accurate implicit-mesh discontinuous Galerkin framework for wave propagation in single-phase and bi-phase solids is presented. The framework belongs to the embedded-boundary techniques and its novelty regards the spatial discretization, which enables boundary and interface conditions to be enforced with high-order accuracy on curved embedded geometries. High-order accuracy is achieved via high-order quadrature rules for implicitly-defined domains and boundaries, whilst a cell-merging strategy addresses the presence of small cut cells. The framework is used to discretize the governing equations of elastodynamics, written using a first-order hyperbolic momentum-strain formulation, and an exact Riemann solver is employed to compute the numerical flux at the interface between dissimilar materials with general anisotropic properties. The space-discretized equations are then advanced in time using explicit high-order Runge–Kutta algorithms. Several two- and three-dimensional numerical tests including dynamic adaptive mesh refinement are presented to demonstrate the high-order accuracy and the capability of the method in the elastodynamic analysis of single- and bi-phases solids containing complex geometries.
Gas switching reforming for hydrogen production (GSR-H2) presents an efficient, low-carbon hydrogen production method that incorporates integrated carbon capture, offering efficiency gains over traditional methods such as proton exchange membrane (PEM) electrolysis, steam methane reforming (SMR) and the newer method of chemical looping reforming (CLR). GSR-H2 has been demonstrated in lab scale which operates as an exothermic process that eliminates the need for additional natural gas combustion, using its own waste heat to generate process steam and partially offset energy usage through electricity production. Beyond its thermal self-sufficiency, GSR-H2 advances upon CLR by integrating all reaction stages within a single reactor cluster, eliminating the complexities of solid circulation, reducing capital costs, and enhancing overall process efficiency. This streamlined design simplifies scale-up and enables inherent CO2 separation with minimal energy penalty, making GSR-H2 a highly competitive pathway for low-carbon hydrogen production. This study presents the first life cycle assessment (LCA) of GSR-H2, offering a novel evaluation of this new process’s environmental impacts across diverse energy scenarios. Key findings reveal that in the renewables-powered scenario, GSR-H2 achieves a GWP of 2.77 kg CO2 eq per kg H2, a substantial improvement over SMR’s 10.4 kg CO2 eq and close to the low emissions of CLR (1.84 kg CO2 eq) and PEM electrolysis (1.85 kg CO2 eq). These results demonstrate GSR-H2’s competitive advantage as a lower-emission alternative, combining design simplicity and efficiency gains, especially in renewable-integrated systems. These results establish GSR-H2 as a competitive, scalable option for hydrogen production, particularly in decarbonization efforts.