Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “molecular replacement”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Linking NH$^+_4$ motion to magnetism in molecular multiferroic (NH 4 ) 2 ⁢[FeCl 5 ⁢(H 2 ⁢O)]: A neutron vibrational spectroscopy study

Here, we present a neutron vibrational spectroscopy study to investigate the influence of NH$^+_4$ motion on the magnetism in [(NH 4 ) 1–x K x ] 2 ⁢[FeCl 5 ⁢(H 2 ⁢O)]. The parent compounds, (NH 4 ) 2 ⁢[FeCl 5 ⁢(H 2 ⁢O)] (x = 0) and K 2 ⁡[FeCl 5 ⁢(H 2 ⁢O)]⁢(x = 1) are isostructural at room temperature, yet displaying drastically different magnetic and multiferroic behavior. K 2 ⁡[FeCl 5 ⁢(H 2 ⁢O)] is nonmultiferroic with type-A collinear antiferromagnetic structure below T N ≈ 14.06 K, whereas (NH 4 ) 2 ⁢[FeCl 5 ⁢(H 2 ⁢O)] is a type-II multiferroic with incommensurate cycloidal spin structure below T FE ≈ 6.8 K. A recent study of the dielectric, structure, and magnetic properties in the mixed [(NH 4 ) 1–x ⁢K x ] 2 ⁢[FeCl 5 ⁢(H 2 ⁢O)] shows that a small amount of potassium substitution to replace NH$^+_4$ transforms the spin structure from incommensurate cycloidal (x ≤ 0.06) into commensurate collinear antiferromagnetic (x ≥ 0.15), indicating NH$^+_4$ is essential to the emergent phenomena observed in this molecular multiferroic compound. Our vibrational spectroscopy study reveals that NH$^+_4$ libration and torsion motion exhibit substantial temperature dependence at low temperatures. The intensity of NH$^+_4$ libration and torsion modes increases slightly at 5 K in comparison with data at 25 K behaving like a magnon, indicating that they are coupled to the magnetism in (NH 4 ) 2 ⁢[FeCl 5 ⁢(H 2 ⁢O)]. Comparing data of x=0, 0.06, 0.09, and 0.15 samples further illustrates that the strength of the increased signal in NH$^+_4$ libration mode is very sensitive to potassium concentration. The signal diminishes quickly with increasing potassium concentration and vanishes in the x = 0.15 sample corresponding to the magnetic structure change for x ≥ 0.15. The results directly link the anomalous behavior in NH$^+_4$ libration motion to the magnetism in [(NH 4 ) 1–x ⁢K x ] 2 ⁢[FeCl 5 ⁢(H 2 ⁢O)], providing new insights into the crucial role NH$^+_4$ plays in the coupled phenomena in (NH 4 ) 2 ⁢[FeCl 5 ⁢(H 2 ⁢O)]. The unique information opens a new door to go through in searching for new multifunctional materials by incorporation of NH 4 via a material-by-design approach.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Molecular-based analysis of nanoparticle solvation: Classical density functional approach

Proper statistical mechanics understanding of nanoparticle solvation processes requires an accurate description of the molecular structure of the solvent. Achieving this goal with standard molecular dynamics (MD) simulation methods is challenging due to large length scales. An alternative approach to this problem can be formulated using classical density functional theory (cDFT), where a full configurational description of the positions of all the atoms is replaced by collective atomic site densities in the molecule. Using an example of the negatively charged silica-like system in an aqueous polar environment represented by a two-site water model, we demonstrate here that cDFT can reproduce MD data at a fraction of the computational cost. An important implication of this result is the ability to understand how the solvent molecular features may affect the system’s properties at the macroscopic scale. A concrete example highlighted in this work is the analysis of nanoparticle interactions with sizes of up to 100 nm in diameter.

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

Molecular structures of residual solvent in polyacrylonitrile based electrolytes: Implications for conductivity and stability

Lithium-ion batteries increasingly play significant roles in modern technologies; however, increased energy density also raises concerns about electrolyte safety. Traditional electrolytes that use volatile organic solvents face risks of thermal runaways and fires from electrode shorting. In response, polymer-based solid electrolytes have been developed for replacement. Polyacrylonitrile (PAN) is a promising fire-resistant component for electrolyte fabrication, but its limited solubility necessitates using low-volatility solvents, which are notoriously difficult to remove in subsequent drying processes. Here, we use femtosecond two-dimensional infrared spectroscopy to provide an in-depth understanding of how residual solvent from processing affects the molecular structures and dynamics within a polymer electrolyte. To this end, linear and nonlinear infrared spectroscopies are employed to interrogate the molecular interactions in PAN-based electrolytes containing various contents of N,N-dimethylformamide (DMF). We show that the amount of DMF within the PAN electrolyte affects the Li+ structure. Further, the coordination can proceed through the carbonyl group and/or the amide nitrogen to form antiparallel structures with the nitrile groups of PAN through dipole–dipole interactions. The free motion of DMF is drastically inhibited upon interaction with Li+ and PAN, which decreases the ionic conductivity and potentially affects the stability (resistance toward removal and chemical decomposition). These findings have implications for the design and processing of solid polymer electrolytes.

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

Increasing thermostability of the key photorespiratory enzyme glycerate 3‐kinase by structure‐based recombination

As global temperatures rise, improving crop yields will require enhancing the thermotolerance of crops. One approach for improving thermotolerance is using bioengineering to increase the thermostability of enzymes catalysing essential biological processes. Photorespiration is an essential recycling process in plants that is integral to photosynthesis and crop growth. The enzymes of photorespiration are targets for enhancing plant thermotolerance as this pathway limits carbon fixation at elevated temperatures. We explored the effects of temperature on the activity of the photorespiratory enzyme glycerate kinase (GLYK) from various organisms and the homologue from the thermophilic alga Cyanidioschyzon merolae was more thermotolerant than those from mesophilic plants, including Arabidopsis thaliana. To understand enzyme features underlying the thermotolerance of C. merolae GLYK (CmGLYK), we performed molecular dynamics simulations using AlphaFold-predicted structures, which revealed greater movement of loop regions of mesophilic plant GLYKs at higher temperatures compared to CmGLYK. Based on these simulations, hybrid proteins were produced and analysed. These hybrid enzymes contained loop regions from CmGLYK replacing the most mobile corresponding loops of AtGLYK. Two of these hybrid enzymes had enhanced thermostability, with melting temperatures increased by 6 °C. One hybrid with three grafted loops maintained higher activity at elevated temperatures. Whilst this hybrid enzyme exhibited enhanced thermostability and a similar Km for ATP compared to AtGLYK, its Km for glycerate increased threefold. This study demonstrates that molecular dynamics simulation-guided structure-based recombination offers a promising strategy for enhancing the thermostability of other plant enzymes with possible application to increasing the thermotolerance of plants under warming climates.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Assessment of software methods for estimating protein-protein relative binding affinities

A growing number of computational tools have been developed to accurately and rapidly predict the impact of amino acid mutations on protein-protein relative binding affinities. Such tools have many applications, for example, designing new drugs and studying evolutionary mechanisms. In the search for accuracy, many of these methods employ expensive yet rigorous molecular dynamics simulations. By contrast, non-rigorous methods use less exhaustive statistical mechanics, allowing for more efficient calculations. However, it is unclear if such methods retain enough accuracy to replace rigorous methods in binding affinity calculations. This trade-off between accuracy and computational expense makes it difficult to determine the best method for a particular system or study. Here, eight non-rigorous computational methods were assessed using eight antibody-antigen and eight non-antibody-antigen complexes for their ability to accurately predict relative binding affinities (ΔΔG) for 654 single mutations. In addition to assessing accuracy, we analyzed the CPU cost and performance for each method using a variety of physico-chemical structural features. This allowed us to posit scenarios in which each method may be best utilized. Most methods performed worse when applied to antibody-antigen complexes compared to non-antibody-antigen complexes. Rosetta-based JayZ and EasyE methods classified mutations as destabilizing (ΔΔG < -0.5 kcal/mol) with high (83–98%) accuracy and a relatively low computational cost for non-antibody-antigen complexes. Some of the most accurate results for antibody-antigen systems came from combining molecular dynamics with FoldX with a correlation coefficient (r) of 0.46, but this was also the most computationally expensive method. Overall, our results suggest these methods can be used to quickly and accurately predict stabilizing versus destabilizing mutations but are less accurate at predicting actual binding affinities. This study highlights the need for continued development of reliable, accessible, and reproducible methods for predicting binding affinities in antibody-antigen proteins and provides a recipe for using current methods.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Composite Correlated Molecular Orbital Theory Calculations of Ring Strain for Use in Predicting Polymerization Reactions

Strained ring systems play an important role in synthesis and can be characterized by the ring strain energy (RSE). The RSE of 3, 4, 5, and 6 membered saturated and unsaturated ring systems containing N, O, P, and S heteroatoms and H, F, SiMe 3 , and SO 2 Me substituents were calculated at the G3(MP2) composite correlated molecular orbital theory level using up to 5 models to predict the RSE. Generally, the RSE decreased as ring size increased with a substantial decrease from 4 to 5 membered rings. Replacement of a ring CH 2 with P or S reduced the RSE, consistent with less angle strain. The RSE for unsaturated systems were generally greater than for saturated systems due to increased angle strain. No general trends were found with respect to substituent effects. In conclusion, the RSE values suggest that 3-pyrroline and 2-pyrroline and their derivatives may be able to support ring opening metathesis polymerization and warrant further study.

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

Subsurface spectroscopy of heterogeneous materials using optical wavefront shaping

Plastic-bonded explosives and propellants consist of energetic molecular crystals embedded in a polymer matrix that also includes additives, such as taggants, plasticizers, grit, antioxidants, etc. Furthermore, this heterogeneous composition renders these materials optically opaque, limiting optical characterization techniques (e.g., Raman spectroscopy) to probing chemical reactions at the sur-face. However, many reactions of interest are believed to occur inside the material where current optical techniques can not probe. To address this challenge, we have developed a new optical technique that utilizes feedback assisted wavefront shaping (using spatial light modulators) to focus probe light inside a heterogeneous material, such that we can perform subsurface spectroscopy. Recently, we demonstrated this technique by monitoring subsurface photodegradation and thermal degradation of Eu-doped molecular crystals using both fluorescence and Raman spectroscopy. Based off these successes we are now looking into improvements to this technique to allow for greater time resolution, including replacing our spatial light modulators with fast digital mirror devices, using ns optical phase conjugation, and self-healing optical beams.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Tuning the melting point of selected ionic liquids through adjustment of the cation’s dipole moment.

In prior research with thermally robust salts [Cassity et al., PCCP, 2017, 19, 31560] it was noted that an increase in the dipole moment of the cation generally led to a decrease in the melting point. Molecular dynamics simulations of the liquid state revealed that an increased dipole moment reduces cation-cation repulsions through dipole-dipole alignment. This was believed to reduce the liquid phase enthalpy, which would tend to lower the melting point of the IL. Here, we also test this principle by replacing hydrogen atoms with fluorine atoms at selected positions within the cation. This allows us to alter the electrostatics of the cation without substantially affecting the sterics. Furthermore, the strength of the dipole moment can be controlled by choosing different positions within the cation for replacement. We studied variants of four different parent cations paired with bistriflimide and determined their melting points, and enthalpies and entropies of fusion through DSC experiments. The decreases in the melting point were determined to be enthalpically driven. We found that the dipole moment of the cation, as determined by quantum chemical calculations, is inversely correlated with the melting point of the given compound. Molecular dynamics simulations of the crystalline and solid states of two isomers showed differences in their enthalpies of fusion that closely matched those seen experimentally. Moreover, this reduction in the enthalpy of fusion was determined to be caused by an increase in the enthalpy of the crystalline state. We provide evidence that dipole-dipole interactions between cations leads to the formation of cationic domains in the crystalline state. These cationic associations partially block favourable cation–anion interactions, which are recovered upon melting. If, however, the dipole- dipole interactions between cations is too strong they have a tendency to form glasses. This study provides a design rule for lowering the melting point of structurally similar ILs by altering their dipole moment.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2026

Environment-imposed selection rules for nuclear-spin conversion of H 2 in molecular crystals

Nuclear-spin conversion in molecular hydrogen is governed by strict symmetry rules that typically require magnetic fields or catalytic surfaces to break. Here we demonstrate that the intrinsic tensor composition of a non-magnetic molecular crystal field can impose and relax these rules without external fields. High-resolution infrared spectra of H 2 in crystalline CO 2 reveal large rank-2 (quadrupolar) crystal-field splittings of the 𝑚 sublevels, while nuclear-spin conversion occurs only through 𝛥𝑚 = 0 channels. Replacing CO 2 with polar N 2 O introduces rank-1 (dipole) components that partially open 𝛥⁢𝑚 ≠ 0 pathways, while incorporation of paramagnetic NO 2 fully lifts the restriction. Furthermore, these results establish a direct correspondence between crystal-field tensor rank and nuclear-spin dynamics, introducing a general symmetry-based framework for designing and controlling spin-isomer populations and quantum-state connectivity in molecular solids.

McLane, Nathan [University of Maryland, College Pa↗
DOE OSTI2026

Investigating PVC polymer–plasticizer interactions with atomistic MD simulations and potential of mean force calculations

For this work, atomistic molecular dynamics (MD) simulations coupled with potential of mean force (PMF) calculations were employed to investigate the interactions between PVC polymer chains containing 6–20 repeating units and various plasticizers, with the goal of identifying potential replacements for the toxic plasticizer di(2-ethylhexyl)phthalate (DEHP) used in blood bags. The selected plasticizers belong to various chemical families, such as orthophthalates, citrates, adipates, and the terephthalate DEHT. Both the polymer and the plasticizers lack ionizable groups, and their interactions are primarily governed by van der Waals and electrostatic forces. A model correlating PMF profiles with interaction forces was developed and validated across polyvinyl chloride (PVC) polymers of different lengths and all investigated plasticizers. This model provides insight into how structural variations in plasticizers influence their respective PMF values. The study ranks the investigated plasticizers based on binding affinity, identifying TOTM (tris(2-ethylhexyl) trimellitate, TEHTM), BTHC (butyryl trihexyl citrate, Citroflex B-6), ATHC (acetyl trihexyl citrate, Citroflex A-6, CA-6), and DEHT (bis(2-ethylhexyl)terephthalate, DOTP/DEHTP) as promising alternatives for further investigation. This comprehensive study encompasses PVC polymers of four different lengths and 14 plasticizers, with all data averaged over 30 independent simulations. The approach provides a deeper understanding of molecular interactions, enabling the tailoring of polymer–plasticizer systems for diverse applications, including extractables and leachables, with relevance spanning materials science to biomedical engineering, and also serves as a basis for developing coarse-grained simulation protocols. Overall, this study provides valuable insights for designing safer and more efficient plasticizer substitutes and is well supported by other studies.

Shet, Sai Athmeeya G. [Sri Sathya Sai Institute of↗
DOE OSTI2023

Non-Empirical and Self-Interaction Corrections for DFTB: Towards Accurate Quantum Simulations for Large Mesoscale Systems (Final Report)

This project was comprised of two complementary (but parallel) thrusts: (1) implementing massively-parallelized computing hardware (with new computational hardware that may replace GPUs) in the density functional tight binding (DFTB) approach and (2) developing new capabilities in DFTB to calculate the electronic structure and dynamics of large chemical systems. While classical molecular dynamics can handle hundreds of thousands of atoms, it cannot provide a first-principles based description of chemical systems at the quantum level. At the other extreme, conventional Kohn-Sham DFT methods can probe the true quantum mechanical nature of chemical systems; however, these methods cannot tackle the large sizes and length scales relevant to dynamics simulations of realistic systems. The DFTB formalism utilized in this project provides a viable approach for probing these large systems at a quantum mechanical level of detail. However, to utilize the DFTB approach for accurate calculations of electronic properties, it is crucial to incorporate quantum-based corrections in DFTB since exchange-correlation effects can still remain very strong in these large systems. At the same time, enhancing the computational efficiency of DFTB is also essential since optimal computational performance is required for addressing the large size scales associated with realistic chemical systems. As such, the new non-empirical corrections and computing hardware enhancements implemented in this project will enable accurate and computationally efficient approaches to directly probe electronic properties in these large, complex systems.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Chemical Kinetic Method for Active-Site Quantification in Fe-N-C Catalysts and Correlation with Molecular Probe and Spectroscopic Site-Counting Methods

Mononuclear Fe ions ligated by nitrogen (FeN x ) dispersed on nitrogen-doped carbon (Fe-N-C) serve as active centers for electrocatalytic O 2 reduction and thermocatalytic aerobic oxidations. Despite their promise as replacements for precious metals in a variety of practical applications, such as fuel cells, the discovery of new Fe-N-C catalysts has relied primarily on empirical approaches. In this context, the development of quantitative structure–reactivity relationships and benchmarking of catalysts prepared by different synthetic routes and by different laboratories would be facilitated by the broader adoption of methods to quantify atomically dispersed FeN x active centers. In the present study, we develop a kinetic probe reaction method that uses the aerobic oxidation of a model hydroquinone substrate to quantify the density of FeN x centers in Fe- N-C catalysts. The kinetic method is compared with low-temperature Mössbauer spectroscopy, CO pulse chemisorption, and electrochemical reductive stripping of NO derived from NO 2 – on a suite of Fe-N-C catalysts prepared by diverse routes and featuring either the exclusive presence of Fe as FeN x sites or the coexistence of aggregated Fe species in addition to FeN x . The FeN x site densities derived from the kinetic method correlate well with those obtained from CO pulse chemisorption and Mössbauer spectroscopy. The broad survey of Fe-N-C materials also reveals the presence of outliers and challenges associated with each site quantification approach. Furthermore, the kinetic method developed here does not require pretreatments that may alter active-site distributions nor specialized equipment beyond reaction vessels and analytical instrumentation (e.g., NMR).

36 MATERIALS SCIENCE↗
DOE OSTI2020

Electrolyte design for Li metal-free Li batteries

Li metal, with the lowest thermodynamically achievable negative electrochemical potential and the highest specific capacity (3860 mAh g -1 ), is the ultimate anode choice for Li batteries. However, the highest reported Li plating/stripping Coulombic efficiency (CE) of 99.5% after extensive efforts is still too low for the Li metal-free (all the Li metal in cycling comes from cathode, without anode pre-lithiation) Li metal batteries. The low CE is attributed to both non-uniform Li plating/stripping on the lithiophobic Cu current collector and Li dendrite growth through lithiophilic organic–inorganic solid electrolyte interphase (SEI) formed in carbonate electrolytes. Here, we use a lithiophilic Bismuth graphite blend (Bi–Gr) substrate to replace lithiophobic Cu current collector to seed a uniform Li nucleation, and form a lithiophobic LiF-rich SEI rather than lithiophilic organic-rich SEI to suppress Li dendrite growth. Molecular dynamics simulations reveal the preferential reduction of anions in 2.0 M LiPF 6 in tetrahydrofuran/2-methyl tetrahydrofuran (2.0 M LiPF 6 –mixTHF) electrolyte to generate LiF-rich SEI on plated Li. Bi–Gr substrate and 2.0 M LiPF 6 –mixTHF electrolyte enable the Li anodes to achieve a record high CE of 99.83% at a high capacity of 1.0 mAh cm -2 and current of 0.5 mA cm -2 . The Bi particles serve as dispersed nucleation centers that promote uniform Li deposition with strong adhesion to the substrate to avoid dead Li, while the lithiophobic LiF-rich SEI promotes lateral Li growth and suppresses the vertical Li dendrite growth even at a high current density of 3.0 mA cm -2 and high areal capacities of 3.0 mAh cm -2 . The regulation of Li nucleation and growth enables the Li metal-free LiFePO 4 full cells to achieve 100 cycles at a practical areal capacity of >2.0 mAh cm -2 . This article highlights the benefits of simultaneous substrate design to improve Li nucleation and electrolyte design to promote lithiophobic SEI growth, enabling a promising and practical route Li metal-free Li metal batteries.

36 MATERIALS SCIENCE↗
DOE OSTI2021

High-level coupled-cluster energetics by merging moment expansions with selected configuration interaction

Inspired by our earlier semi-stochastic work aimed at converging high-level coupled-cluster (CC) energetics, we propose a novel form of the CC(P; Q) theory in which the stochastic Quantum Monte Carlo propagations, used to identify dominant higher-than-doubly excited determinants, are replaced by the selected configuration interaction (CI) approach using the perturbative selection made iteratively (CIPSI) algorithm. The advantages of the resulting CIPSI-driven CC(P; Q) methodology are illustrated by a few molecular examples, including the dissociation of F 2 and the automerization of cyclobutadiene, where we recover the electronic energies corresponding to the CC calculations with a full treatment of singles, doubles, and triples based on the information extracted from compact CI wave functions originating from relatively inexpensive Hamiltonian diagonalizations.

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

Interlayer pillaring influences the octahedral tilting and electrochemical capacity of tungsten oxides

One strategy to tune the electrochemical properties of layered inorganic materials for energy storage and conversion is to introduce molecular organic pillars. We investigated how pillaring tungsten oxides with alkylammonium cations influences their physical and electrochemical properties relative to the non-pillared host. While the presence of alkylammonium increased the interlayer spacing, we found that the electrochemical capacity of the materials decreased in non-aqueous electrolytes. In aqueous acidic electrolytes, replacement of the interlayer alkylammonium pillars with water molecules led to the recovery of the electrochemical capacity. We rationalize these results based on the influence of interlayer pillars on the degree of octahedral tilting within the inorganic tungsten oxide layers. The presence of alkylammonium in the interlayer leads to an increase in octahedral tilting relative to interlayer water. In turn, this leads to an increase in the band gap and corresponding decrease in the electrochemical capacity. Our findings emphasize the correlation between the chemical nature of the interlayer molecular pillars and the electronic structure, which in turn affects the electrochemical capacity.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Elucidate the molecular basis of ampholytic chitosan as a high-performance cryoprotectant to myosin denaturation: The importance of saccharide charges

The uses of charged poly/oligo-saccharides enable the retarding of protein denaturation against various environmental stresses during food storage and manufacturing. However, at subzero temperatures, the molecular basis of such stabilization behaviors, i.e., cryoprotections, remain less explored. Here, in this study, we introduced an ampholytic saccharide, carboxymethyl chitooligosaccharide (CMCO) that effectively inhibited the freezing-induced myosin denaturation. The in-depth cryoprotective mechanism was systematically investigated by using molecular dynamic simulation and multispectral characterizations. Results showed that CMCO may interact with myosin through hydrogen bonding and electrostatic interactions, which caused the expelling of water at protein surfaces and the reduced conformational flexibility of myosin molecules. Due to this water replacement event, both secondary and ternary structures of myosin became freezing-resistant, leading to the inhibited protein aggregations and retained functionalities, such as solubility, Ca 2+ -ATPase activity, and gelling properties. Moreover, cryoprotective behaviors of CMCO were charge-dependent. CMCO with a higher degree of carboxymethyl substitution (DS: 1.2) was inclined to bind and stabilize myosin molecules better than the low-DS (DS: 0.8) one, even though both outperformed other cryoprotective saccharides. Therefore, this investigation not only introducing a high-performance myosin cryoprotectant, but also elaborated the cryoprotective mechanism of ampholytic saccharides.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

High Pressure Infiltration/Expulsion of Aqueous NaCl in Planar Hydrophobic Nanopores

Pressure-driven permeation of water in a poorly wettable material results in a conversion of mechanical work into surface free energy representing a new form of energy storage, or absorption. When water is replaced by a concentrated electrolyte solution, the storage capacity of a nanoporous medium becomes comparable to high-end supercapacitors. The addition of salt can also reduce the hysteresis of the infiltration/expulsion cycle. Our molecular simulations provide a theoretical perspective into the mechanisms involved in the process, and underlying structures and interactions in compressed nanoconfined solutions. Specifically, we consider aqueous NaCl in planar confinements of widths of 1.0 nm and 1.64 nm and pressures of up to 3 kbar. Open ensemble Monte Carlo simulations utilizing fractional exchanges of molecules for efficient additions/removal of ions have been utilized in conjunction with pressure-dependent chemical potentials to model bulk phases under pressure. Confinements open to these pressurized bulk, aqueous electrolyte phases show reversibility at narrow pore sizes and irreversibility in wider ones, consistent with experiment, as well as strong hysteresis at both pore size. The addition of salt results in significant increases in the solid/liquid interfacial tension in narrower pores and associated infiltration and expulsion pressures. Furthermore, these changes are consistent with strong desalination effects at the lower pore size observed irrespective of external pressure and initial concentration.

25 ENERGY STORAGE↗
DOE OSTI2020

Analysis of an all-solid state nanobattery using molecular dynamics simulations under an external electric field

Present Li-ion battery (LIB) technology requires strong improvements in performance, energy capacity, charging-time, and cost to expand their application to e-mobility and grid storage. Li-metal is one of the most promising materials to replace commercial anodes such as graphite because of its 10 times higher specific capacity. However, Li-metal has high reactivity with commercial liquid electrolytes; thus, new solid materials are proposed to replace liquid electrolytes when Li-metal anodes are used. We present a theoretical analysis of the charging process in a full nanobattery, containing a LiCoO 2 cathode, a Li 7 P2S 8 I solid-state electrolyte (SSE), a Li-metal anode as well as Al and Cu collectors for the cathode and anode, respectively. In addition, we added a Li 3 P/Li 2 S film as a solid electrolyte interphase (SEI) layer between the Li-anode and SSE. Thus, we focus this study on the SEI and SSE. We simulated the charging of the nanobattery with an external voltage by applying an electric field. We estimated temperature profiles within the nanobattery and analyzed Li-ion transport through the SSE and SEI. Here, we observed a slight temperature rise at the SEI due to reactions forming $PS_{3}^{–}$ and $P_{2}S_{7}^{4}$$^{–}$ fragments at the interfaces; however, this temperature profile changes due to the charging current under the presence of the external electric field ε = 0.75 V Å –1 . Without the external field, the calculated open-circuit voltage (OCV) was 3.86 V for the battery, which is within the range of values of commercial cobalt-based LIBs. This voltage implies a spontaneous fall of available Li-ions from the anode to the cathode (during discharge). The charge of this nanobattery requires overcoming the OCV plus an additional voltage that determines the charging current. Thus, we applied an external potential able to neutralize the OCV, plus an additional 1.6 V to induce the transport of Li + from the cathode up to the anode. Several interesting details about Li + transport paths through the SSE and SEI are discussed.

25 ENERGY STORAGE↗