Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “structural modeling”

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 451 records · Page 25

rcsb-api : Python Toolkit for Streamlining Access to RCSB Protein Data Bank APIs

The Protein Data Bank (PDB) was founded in 1971 as the first open-access digital data resource in biology to serve as the single global archive for three-dimensional (3D) macromolecular structure data. Current PDB holdings exceed 230,000 experimentally determined structures of proteins, nucleic acids, viruses, and macromolecular machines. The RCSB Protein Data Bank RCSB.org research-focused web portal facilitates search, analyses, and visualization of every PDB structure along with more than one million Computed Structure Models from AlphaFold DB and the ModelArchive. It is powered by a set of publicly available Application Programming Interfaces (APIs) that both support RCSB.org users and provide programmatic access to PDB data. Given the breadth and levels of granularity encompassed in this rich data collection, efficiently accessing the information programmatically may be challenging for new users. RCSB PDB has developed a Python software package, rcsb-api , that facilitates easy and efficient use of RCSB PDB APIs within a Python environment. This software tool is designed to streamline access to the extensive corpus of data housed within the PDB, enabling researchers to search, retrieve, and analyze 3D biostructure data seamlessly. Its use will accelerate research in structural biology, molecular biology and biochemistry, drug discovery, and bioinformatics by providing more efficient tools for data integration and analysis. The new toolkit is available on GitHub (github.com/rcsb/py-rcsb-api) and published to the public Python package repository (PyPI) to foster wider usage and support basic and applied research in fundamental biology, biomedicine, and the energy sciences.

FAIR principles↗

Understanding and modeling current and future coastal wetland methane dynamics (Final Report)

The coastal terrestrial-aquatic interface (TAI) is a highly dynamic component of the Earth system that plays a critical role in biogeochemical cycling. Due to its dynamic nature, the processes that regulate decomposition and methane (CH 4 ) emissions are of greater significance at the TAI than in upland systems. Despite this, we have limited mechanistic understanding of how climate stressors interact to regulate the electron acceptors and donors that determine decomposition pathways within TAIs, including the generation of hot spots and hot moments. Accurately modeling these processes is critical for incorporating the coastal TAI into Earth systems models, such as DOE’s E3SM. With previous DOE support, we adapted an aerobic terrestrial representation of decomposition using PFLOTRAN, a reactive flow and transport model, and added anerobic decomposition pathways, salinity, and oxygen (O 2 ) that fluctuates independently of water table level. However, because this model (PFLOTRAN TAI ) is based on decomposition rates and organic matter carbon to nitrogen ratios from terrestrial systems, its performance in TAI systems is currently limited by the lack of empirical data to properly parameterize variables. In addition, while PFLOTRAN TAI can simulate movement of O 2 into sediments, it is not currently capable of tracking the movement of CH 4 gas through plant tissues due to both current model structure and lack of available data. Our overall objective of the project was the increase our mechanistic understanding of CH 4 dynamics in response to environmental change, such that we can improve the representation of these dynamics in PFLOTRAN. We installed automated flux chambers in a new field-scale active soil warming experiment and set up marsh organs (mesocosms) to test effects of warming, flooding, and salinity. This resulted in a new dataset consisting of chamber-level CH 4 flux measurements across multiple sites, ecological conditions, and timeframes, as well as corresponding measurements on porewater chemistry, soil carbon quality, plant biomass, and redox reaction rates. Using these data, we improved PFLOTRAN TAI to more accurately model CH 4 dynamics and successfully tested our hypotheses. This grant contributed to the professional development of 3 postdocs, 4 undergraduate interns, 4 teacher interns, 17 technicians, and 15 participatory scientists. The automated chamber technology designed for this grant has also been shared with multiple new projects.

54 ENVIRONMENTAL SCIENCES↗

Process-based balance of system cost modeling for offshore wind power plants in the United States

This paper describes the development of a process-based and open-source balance of system cost model that provides the capability to evaluate both existing and novel offshore wind technologies. Individual design and installation steps are represented with bottom-up engineering models that compute times and costs associated with the process; furthermore, operational constraints are assigned to each process so that delays caused by weather and presence of marine mammals may be accounted for in the overall project timeline. The model structure, assumptions, inputs, and results are vetted with industry partners and compared against actual projects for validation. Installation times show reasonable agreement with real data. Project cost sensitivities are investigated to compute the system-level impact of different design choices. First, individual vessel efficiencies are computed for varying numbers of installation vessels and weather time series to show the diminishing returns of more than two feeder barges. Then, array cable capital costs and installation times are determined for a representative project with different turbine sizes. These values quantify the cost-benefit tradeoffs and show a net-cost savings of decreasing numbers of turbines, increased turbine spacing, and fewer turbine terminations. These results demonstrate that the balance of system model features the accuracy, functionality, and accessibility to serve as the foundation for a wide range of analyses to identify cost reduction potentials for offshore wind energy in the United States.

17 WIND ENERGY↗

Engineering Robust Metallic Zero-Mode States in Olympicene Graphene Nanoribbons

Metallic graphene nanoribbons (GNRs) represent a critical component in the toolbox of low-dimensional functional materials technology serving as 1D interconnects capable of both electronic and quantum information transport. The structural constraints imposed by on-surface bottom-up GNR synthesis protocols along with the limited control over orientation and sequence of asymmetric monomer building blocks during the radical step-growth polymerization have plagued the design and assembly of metallic GNRs. Here, we report the regioregular synthesis of GNRs hosting robust metallic states by embedding a symmetric zero-mode (ZM) superlattice along the backbone of a GNR. Tight-binding electronic structure models predict a strong nearest-neighbor electron hopping interaction between adjacent ZM states, resulting in a dispersive metallic band. First-principles density functional theory-local density approximation calculations confirm this prediction, and the robust, metallic ZM band of olympicene GNRs is experimentally corroborated by scanning tunneling spectroscopy.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Structure determination of ζ-N 2 from single-crystal X-ray diffraction and theoretical suggestion for the formation of amorphous nitrogen

The allotropy of solid molecular nitrogen is the consequence of a complex interplay between fundamental intermolecular as well as intramolecular interactions. Understanding the underlying physical mechanisms hinges on knowledge of the crystal structures of these molecular phases. That is especially true for ζ-N 2 , key to shed light on nitrogen’s polymerization. Here, we perform single-crystal X-ray diffraction on laser-heated N 2 samples at 54, 63, 70 and 86 GPa and solve and refine the hitherto unknown structure of ζ-N 2 . In its monoclinic unit cell (space group C 2/ c ), 16 N 2 molecules are arranged in a configuration similar to that of ε-N 2 . The structure model provides an explanation for the previously identified Raman and infrared lattice and vibrational modes of ζ-N 2 . Density functional theory calculations give an insight into the gradual delocalization of electronic density from intramolecular bonds to intermolecular space and suggest a possible pathway towards nitrogen’s polymerization.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The local ordering of polar solvents around crystalline carbon nitride nanosheets in solution

The crystalline graphitic carbon nitride, poly-triazine imide (PTI) is highly unusual among layered materials since it is spontaneously soluble in aprotic, polar solvents including dimethylformamide (DMF). The PTI material consists of layers of carbon nitride intercalated with LiBr. When dissolved, the resulting solutions consist of dissolved, luminescent single to multilayer nanosheets of around 60–125 nm in diameter and Li+ and Br- ions originating from the intercalating salt. To understand this unique solubility, the structure of these solutions has been investigated by high-energy X-ray and neutron diffraction. Although the diffraction patterns are dominated by inter-solvent correlations there are clear differences between the X-ray diffraction data of the PTI solution and the solvent in the 4–6 Å -1 range, with real space differences persisting to at least 10 Å. Structural modelling using both neutron and X-ray datasets as a constraint reveal the formation of distinct, dense solvation shells surrounding the nanoparticles with a layer of Br - close to the PTI-solvent interface. This solvent ordering provides a configuration that is energetically favourable underpinning thermodynamically driven PTI dissolution.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Gate-based quantum computing for protein design

Protein design is a technique to engineer proteins by permuting amino acids in the sequence to obtain novel functionalities. However, exploring all possible combinations of amino acids is generally impossible due to the exponential growth of possibilities with the number of designable sites. The present work introduces circuits implementing a pure quantum approach, Grover’s algorithm, to solve protein design problems. Our algorithms can adjust to implement any custom pair-wise energy tables and protein structure models. Moreover, the algorithm’s oracle is designed to consist of only adder functions. Quantum computer simulators validate the practicality of our circuits, containing up to 234 qubits. However, a smaller circuit is implemented on real quantum devices. Our results show that using iterations, the circuits find the correct results among all N possibilities, providing the expected quadratic speed up of Grover’s algorithm over classical methods (i.e.,).

59 BASIC BIOLOGICAL SCIENCES↗

Numerical study of the effect of propped surface area and fracture conductivity on shale gas production: Application for multi-size proppant pumping schedule design

Hydraulic fracturing is a technique extensively used in the oil and gas industry, where water, proppant (sand) and additives are injected into unconventional reservoirs to enhance the recovery of shale hydrocarbon. Although some previous studies have developed pumping schedules that maximize gas production for a single-size proppant, there are very few studies that consider the effect of varying proppant diameters across pumping stages on shale gas production. Motivated by this, we carried out an extensive sensitivity analysis to determine the effect of different proppant diameters on the average fracture conductivity (FC), average propped surface area (PSA) and cumulative shale gas production volume. We found out that the cumulative shale gas production volume depends on both the average PSA and average FC. We also found out that small-diameter proppant resulted in higher average PSA and lower average FC, whereas large-diameter proppant resulted in lower average PSA and higher average FC. Hence, we designed a multi-size proppant pumping schedule considering both of these parameters into account for simultaneously propagating multiple fractures to maximize shale gas production from unconventional reservoirs. Since the size of injected proppant particles determines the average PSA and average FC for the propped hydraulic fractures, we developed a novel framework called Sequentially Interlinked Modeling Structure (SIMS) to predict the average PSA, average FC and cumulative shale gas production volume at the end of 10 years for a given pumping schedule. Then, we used this SIMS framework to obtain a multi-size proppant pumping schedule that maximizes shale gas production. In closing, we demonstrated that the obtained pumping schedule gives a cumulative shale gas production volume greater than the values obtained from the existing pumping schedules.

03 NATURAL GAS↗

3D Printed Parahydrophobic Surfaces as Multireaction Platforms

We report Parahydrophobic surfaces (PHSs) composed of arrays of cubic μ-pillars with a double scale of roughness and variable wettability were systematically obtained in one step and a widely accessible stereolithographic Formlabs 3D printer. The wettability control was achieved by combining the geometrical parameters (H = height and P = pitch) and the surface modification with fluoroalkyl silane compounds. Homogeneous distribution of F and Si atoms onto the pillars was observed by XPS and SEM-EDAX. A nano-roughness on the heads of the pillars was achieved without any post-treatment. The smallest P values lead to surfaces with static contact angles (CAs) >150° regardless of the H utilized. Interestingly, the relationship 0.6 ≤ H/P ≤ 2.6 obtained here was in good agreement with the H/P values reported for nano- and submicron pillars. Furthermore, experimental CAs, advancing and receding CAs, were consistent with the theoretical prediction from the Cassie–Baxter model. Structures covered with perfluorodecyltriethoxysilane with high H and short P lead to PHSs. Conversely, structures covered with perfluorodecyltrimethoxysilane exhibited a superhydrophobic behavior. Finally, several aqueous reactions, such as precipitation, coordination complex, and nanoparticle synthesis, were carried out by placing the reactive agents as microdroplets on the parahydrophobic pillars, demonstrating the potential application as chemical multi-reaction array platforms for a large variety of relevant fields in microdroplet manipulation, microfluidics systems, and health monitoring, among others.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Bivalent recognition of fatty acyl-CoA by a human integral membrane palmitoyltransferase

S-acylation, also known as palmitoylation, is the most abundant form of protein lipidation in humans. This reversible posttranslational modification, which targets thousands of proteins, is catalyzed by 23 members of the DHHC family of integral membrane enzymes. DHHC enzymes use fatty acyl-CoA as the ubiquitous fatty acyl donor and become autoacylated at a catalytic cysteine; this intermediate subsequently transfers the fatty acyl group to a cysteine in the target protein. Protein S-acylation intersects with almost all areas of human physiology, and several DHHC enzymes are considered as possible therapeutic targets against diseases such as cancer. These efforts would greatly benefit from a detailed understanding of the molecular basis for this crucial enzymatic reaction. Here, in this study, we combine X-ray crystallography with all-atom molecular dynamics simulations to elucidate the structure of the precatalytic complex of human DHHC20 in complex with palmitoyl CoA. The resulting structure reveals that the fatty acyl chain inserts into a hydrophobic pocket within the transmembrane spanning region of the protein, whereas the CoA headgroup is recognized by the cytosolic domain through polar and ionic interactions. Biochemical experiments corroborate the predictions from our structural model. We show, using both computational and experimental analyses, that palmitoyl CoA acts as a bivalent ligand where the interaction of the DHHC enzyme with both the fatty acyl chain and the CoA headgroup is important for catalytic chemistry to proceed. This bivalency explains how, in the presence of high concentrations of free CoA under physiological conditions, DHHC enzymes can efficiently use palmitoyl CoA as a substrate for autoacylation.

59 BASIC BIOLOGICAL SCIENCES↗

CO2 response screen in grass Brachypodium reveals the key role of a MAP kinase in CO2-triggered stomatal closure

Abstract Plants respond to increased CO2 concentrations through stomatal closure, which can contribute to increased water use efficiency. Grasses display faster stomatal responses than eudicots due to dumbbell-shaped guard cells flanked by subsidiary cells working in opposition. However, forward genetic screening for stomatal CO2 signal transduction mutants in grasses has yet to be reported. The grass model Brachypodium distachyon is closely related to agronomically important cereal crops, sharing largely collinear genomes. To gain insights into CO2 control mechanisms of stomatal movements in grasses, we developed an unbiased forward genetic screen with an EMS-mutagenized B. distachyon M5 generation population using infrared imaging to identify plants with altered leaf temperatures at elevated CO2. Among isolated mutants, a “chill1” mutant exhibited cooler leaf temperatures than wild-type Bd21-3 parent control plants after exposure to increased CO2. chill1 plants showed strongly impaired high CO2-induced stomatal closure despite retaining a robust abscisic acid-induced stomatal closing response. Through bulked segregant whole-genome sequencing analyses followed by analyses of further backcrossed F4 generation plants and generation and characterization of sodium azide and CRISPR-cas9 mutants, chill1 was mapped to a protein kinase, Mitogen-Activated Protein Kinase 5 (BdMPK5). The chill1 mutation impaired BdMPK5 protein-mediated CO2/HCO3− sensing together with the High Temperature 1 (HT1) Raf-like kinase in vitro. Furthermore, AlphaFold2-directed structural modeling predicted that the identified BdMPK5-D90N chill1 mutant residue is located at the interface of BdMPK5 with the BdHT1 Raf-like kinase. BdMPK5 is a key signaling component that mediates CO2-induced stomatal movements and is proposed to function as a component of the primary CO2 sensor in grasses.

Lopez, Bryn N. K. (ORCID:0009000937288216)↗

Progress on Demonstration of a MOOSE-Based Coupled Capability for Hot Channel Factors in Fast Reactors

Hot channel factors (HCFs) are computed values that account for the impact on predicted peak fuel, cladding, and coolant temperatures due to uncertainties in the as-built reactor’s material properties and geometry as well as uncertainties due to modeling approximations. Reduction in computed HCF values via reduction or elimination of modeling approximations may translate to significant economic savings if the reactor power can be raised due to the extra temperature margin gained. While limited historical datasets exist for sodium-cooled fast reactors (SFRs), there are no available HCF data for lead-cooled fast reactors (LFRs) outside of work generated previously within NEAMS. The computation of HCFs involves insights from reactor physics, thermal fluids and heat conduction calculations to determine how the peak temperatures respond to various uncertainties in the design. Due to the significant advantages for multi-physics coupling offered by the MOOSE framework, Griffin (MOOSE-based reactor physics code), MOOSE Heat Conduction Module, and Cardinal (MOOSE-wrapped multi-physics application which includes the NekRS thermal fluids code) are being coupled together using the MOOSE MultiApp System to develop a highfidelity multi-physics modeling capability for HCF simulations. This high-fidelity coupling workflow may also be beneficial for other fast reactor applications in the future. In previous work, Griffin and NekRS were individually assessed to ensure the necessary capabilities were in place. This work describes initial efforts to couple the codes (including folding in the MOOSE Heat Conduction Module) and determining the workflow for the perturbed calculations which will leverage the Stochastic Tools Module (STM). To our knowledge, this is the first coupling of Griffin and NekRS as well as the first exploratory use of Stochastic Tools Module for Cardinal. In this report, the neutronics code Griffin, the heat conduction solver in MOOSE, and the MOOSE-wrapped application containing NekRS (Cardinal) are linked together to demonstrate the coupled capability. Griffin and Cardinal are linked dynamically by specifying shared libraries. Different coupling hierarchies are tested for selecting the most appropriate coupling strategy. A coupling scheme is selected based on the efficiency of calculation and ease of data communication. Multiple tests are performed to choose suitable mesh structure, model configurations, scheme setup and boundary conditions to avoid loss of energy due to data interpolation between different modules or weak imposition of fluxes in finite element codes. Computational experiments are performed to study the tolerance control of each type of iteration to avoid false convergence. The coupled capability is demonstrated in both single pin and 7-pin models based on LFR materials and geometry. The study finds that the use of too large a time step size in the heat conduction module can lead to temperature oscillation even though the heat conduction equation does not have a time-derivative kernel, but only the time-dependent boundary condition. A 7-pin model without duct region achieved good convergence in the coupled calculation while a 7 pin model with duct region experienced data communication issues which need to be resolved.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Probing the Local Site Disorder and Distortion in Pyrochlore High-Entropy Oxides

High-entropy oxides (HEOs) have attracted great interest in diverse fields because of their inherent opportunities to tailor and combine materials functionalities. The control of local order/disorder in the class is by extension a grand challenge toward realizing their vast potential. In this paper, we report the first examples of pyrochlore HEOs with five M-site cations, for Nd 2 M 2 O 7 , in which the local structure has been investigated by neutron diffraction and pair distribution function (PDF) analysis. The average structure of the pyrochlores is found to be orthorhombic Imma, in agreement with radius-ratio rules governing the structural archetype. The computed PDFs from density functional theory relaxed special quasirandom structure models are compared with real space PDFs in this work to evaluate M-site order/disorder. Reverse Monte Carlo combined with ab initio molecular dynamics and Metropolis Monte Carlo simulations demonstrates that Nd 2 (Ta 0.2 Sc 0.2 Sn 0.2 Hf 0.2 Zr 0.2 ) 2 O 7 is synthesized with its M-site local to nanoscale order highly randomized/disordered, while Nd 2 (Ti 0.2 Nb 0.2 Sn 0.2 Hf 0.2 Zr 0.2 ) 2 O 7+x exhibits a strong distortion of the TiO 6 octahedron and small degree of Ti chemical short-range order (SRO) on the subnanometer scale. Calculations suggest that this may be intrinsic, energetically favored SRO rather than due to sample processing. These results offer an important demonstration that the engineered variation of participating ions in HEOs, even among those with very similar radii, provides richly diverse opportunities to control local order/disorder motifs—and therefore materials properties for future designs. This work also hints at the exquisite level of detail that may be needed in computational and experimental data analysis to guide structure–property tuning in the emerging HEO materials class.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Reaching the potential of electron diffraction

Microcrystal electron diffraction (MicroED) is an emerging structural technique in which submicron crystals are used to generate diffraction data for structural studies. Structures allow for the study of molecular-level architecture and drive hypotheses about modes of action, mechanisms, dynamics, and interactions with other molecules. Combining cryoelectron microscopy (cryo-EM) instrumentation with crystallographic techniques, MicroED has led to three-dimensional structural models of small molecules, peptides, and proteins and has generated tremendous interest due to its ability to use vanishingly small crystals. In this perspective, we describe the current state of the field for MicroED methodologies, including making and detecting crystals of the appropriate size for the technique, as well as ways to best handle and characterize these crystals. Our perspective provides insight into ways to unlock the full range of potential for MicroED to access previously intractable samples and describes areas of future development.

3D ED↗

Interdependence in active mobility adoption: Joint modeling and motivational spillover in walking, cycling and bike-sharing

Active mobility offers an array of physical, emotional, and social well-being benefits. However, with the proliferation of the sharing economy, new nonmotorized means of transport are entering the fold, complementing some existing mobility options while competing with others. The purpose of this research study is to investigate the adoption of three active travel modes—namely walking, cycling, and bikesharing—in a joint modeling framework. Here, the analysis is based on an adaptation of the stages of change framework, which originates from the health behavior sciences. Multivariate ordered probit modeling drawing on U.S. survey data provides well-needed insights into individuals’ preparedness to adopt multiple active modes as a function of personal, neighborhood, and psychosocial factors. The research suggests three important findings. (1) The joint model structure confirms interdependence among different active mobility choices. The strongest complementarity is found for walking and cycling adoption. (2) Each mode has a distinctive adoption path with either three or four separate stages. We discuss the implications of derived stage-thresholds and plot adoption contours for selected scenarios. (3) Psychological and neighborhood variables generate more coupling among active modes than individual and household factors. Specifically, identifying strongly with active mobility aspirations, experiences with multimodal travel, possessing better navigational skills, along with supportive local community norms are the factors that appear to drive the joint adoption decisions. This study contributes to the understanding of how decisions within the same functional domain are related and help to design policies that promote active mobility by identifying positive spillovers and joint determinants.

42 ENGINEERING↗

Exotic bonding interactions and coexistence of chemically distinct periodic lattice distortions in the charge density wave compound TaTe2

The superstructure of TaTe 2 is examined by total high-energy x-ray scattering and large-scale structure modeling over a broad temperature range. At room temperature, it features double zigzag chains of Ta atoms exhibiting metallic bonding. The chains are sandwiched between planes of weakly interacting Te atoms. Below 170 K, the chains appear broken to “butterflylike” clusters and single zigzag chains of covalently bonded Te atoms emerge while Ta-Te interactions remain largely ionic. This is a rare example of a charge density wave compound stabilized by bonding that is a mix of covalent, ionic, and metallic character giving rise to chemically distinct periodic lattice distortions. We argue that the formation of periodic lattice distortions in the Te sublattice in addition to those in the Ta sublattice favors charge delocalization, which may explain the unusually diminished resistivity and increased magnetic susceptibility exhibited by the low-temperature charge density wave phase of TaTe 2 . The effect may be common for 5d transition metal chalcogenides containing Te because of the extended electronic orbitals of the constituent atoms.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Resolving Nonequilibrium Shape Variations among Millions of Gold Nanoparticles

Nanoparticles, exhibiting functionally relevant structural heterogeneity, are at the forefront of cutting-edge research. Now, high-throughput single-particle imaging (SPI) with X-ray free-electron lasers (XFELs) creates opportunities for recovering the shape distributions of millions of particles that exhibit functionally relevant structural heterogeneity. To realize this potential, three challenges have to be overcome: (1) simultaneous parametrization of structural variability in real and reciprocal spaces; (2) efficiently inferring the latent parameters of each SPI measurement; (3) scaling up comparisons between 10 5 structural models and 10 6 XFEL-SPI measurements. Here, we describe how we overcame these three challenges to resolve the nonequilibrium shape distributions within millions of gold nanoparticles imaged at the European XFEL. These shape distributions allowed us to quantify the degree of asymmetry in these particles, discover a relatively stable “shape envelope” among nanoparticles, discern finite-size effects related to shape-controlling surfactants, and extrapolate nanoparticles’ shapes to their idealized thermodynamic limit. Ultimately, these demonstrations show that XFEL SPI can help transform nanoparticle shape characterization from anecdotally interesting to statistically meaningful.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Critical Role of Framework Flexibility and Disorder in Driving High Ionic Conductivity in LiNbOCl 4

Understanding Li-ion transport is key for the rational design of superionic solid electrolytes with exceptional ionic conductivities. LiNbOCl 4 is reported to be one of the most highly conducting materials in the recently realized new class of soft oxyhalide solid electrolytes, exhibiting an ionic conductivity of ~11 mS·cm -1 . Here, we apply X-ray/neutron diffraction and pair distribution function analysis - coupled with density functional theory/ab-initio molecular dynamics - to determine a structural model that provides a rationale for the high conductivity that we observe experimentally in this nanocrystalline solid. We show that it arises from unusually high framework flexibility at room temperature. This owes to isolated 1-D [NbOCl 4 ] - anionic chains which exhibit energetically favorable orientational disorder that is - in turn - correlated to multiple, disordered and equi-energetic Li + sites in the lattice. As the Li-ions sample the 3-D energy landscape with a fast predicted diffusion coefficient of 5.1 x 10 -7 cm 2 /s at room temperature (σ i calc = 17.4 mS·cm -1 ), the inorganic polymer chains can reorient or vice versa. The activation energy barrier for Li migration through the frustrated energy landscape is especially reduced by the elastic nature of the NbO 2 Cl 4 octahedra evident from very widely dispersed Cl-Nb-Cl bond angles in AIMD snapshots at 300 K. The phonon spectra are predominantly influenced by Cl vibrations in the low energy range, and there is strong overlap between the framework (Cl, Nb) and Li partial pDOS in the region between 1.2 - 4.0 THz. The framework flexibility is also reflected in a relatively low bulk modulus of 22 GPa. In conclusion, our findings pave the way for investigation of future “flex-ion” inorganic solids and open up a new direction for the design of high conductivity, soft solid electrolytes for all-solid-state batteries.

AIMD↗