SEARCH · Engineering Papers

Results for “reversible computation”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 955 records · Page 53

All sources NASA NTRS DOE OSTI

NASA NTRS2011

Planetary Suit Hip Bearing Model for Predicting Design vs. Performance

Designing a planetary suit is very complex and often requires difficult trade-offs between performance, cost, mass, and system complexity. In order to verifying that new suit designs meet requirements, full prototypes must eventually be built and tested with human subjects. Using computer models early in the design phase of new hardware development can be advantageous, allowing virtual prototyping to take place. Having easily modifiable models of the suit hard sections may reduce the time it takes to make changes to the hardware designs and then to understand their impact on suit and human performance. A virtual design environment gives designers the ability to think outside the box and exhaust design possibilities before building and testing physical prototypes with human subjects. Reductions in prototyping and testing may eventually reduce development costs. This study is an attempt to develop computer models of the hard components of the suit with known physical characteristics, supplemented with human subject performance data. Objectives: The primary objective was to develop an articulating solid model of the Mark III hip bearings to be used for evaluating suit design performance of the hip joint. Methods: Solid models of a planetary prototype (Mark III) suit s hip bearings and brief section were reverse-engineered from the prototype. The performance of the models was then compared by evaluating the mobility performance differences between the nominal hardware configuration and hardware modifications. This was accomplished by gathering data from specific suited tasks. Subjects performed maximum flexion and abduction tasks while in a nominal suit bearing configuration and in three off-nominal configurations. Performance data for the hip were recorded using state-of-the-art motion capture technology. Results: The results demonstrate that solid models of planetary suit hard segments for use as a performance design tool is feasible. From a general trend perspective, the suited performance trends were comparable between the model and the suited subjects. With the three off-nominal bearing configurations compared to the nominal bearing configurations, human subjects showed decreases in hip flexion of 64%, 6%, and 13% and in hip abduction of 59%, 2%, and 20%. Likewise the solid model showed decreases in hip flexion of 58%, 1%, and 25% and in hip abduction of 56%, 0%, and 30%, under the same condition changes from the nominal configuration. Differences seen between the model predictions and the human subject performance data could be attributed to the model lacking dynamic elements and performing kinematic analysis only, the level of fit of the subjects with the suit, the levels of the subject s suit experience.

Cowley, Matthew S.↗

NASA NTRS2015

Comparison of Numerically Simulated and Experimentally Measured Performance of a Rotating Detonation Engine

A quasi-two-dimensional, computational fluid dynamic (CFD) simulation of a rotating detonation engine (RDE) is described. The simulation operates in the detonation frame of reference and utilizes a relatively coarse grid such that only the essential primary flow field structure is captured. This construction and other simplifications yield rapidly converging, steady solutions. Viscous effects, and heat transfer effects are modeled using source terms. The effects of potential inlet flow reversals are modeled using boundary conditions. Results from the simulation are compared to measured data from an experimental RDE rig with a converging-diverging nozzle added. The comparison is favorable for the two operating points examined. The utility of the code as a performance optimization tool and a diagnostic tool are discussed.

combustion↗

DOE OSTI2014

Changes in Functional Lung Regions During the Course of Radiation Therapy and Their Potential Impact on Lung Dosimetry for Non-Small Cell Lung Cancer

Purpose: To study changes in functional activity on ventilation (V)/perfusion (Q) single-photon emission computed tomography (SPECT) during radiation therapy (RT) and explore the impact of such changes on lung dosimetry in patients with non-small cell lung cancer (NSCLC). Methods and Materials: Fifteen NSCLC patients with centrally located tumors were enrolled. All patients were treated with definitive RT dose of ≥60 Gy. V/Q SPECT-CT scans were performed prior to and after delivery of 45 Gy of fractionated RT. SPECT images were used to define temporarily dysfunctional regions of lung caused by tumor or other potentially reversible conditions as B3. The functional lung (FL) was defined on SPECT by 2 separate approaches: FL1, a threshold of 30% of the maximum uptake of the patient's lung; and FL2, FL1 plus B3 region. The impact of changes in FL between initiation of RT and delivery of 45 Gy on lung dosimetry were analyzed. Results: Fourteen patients (93%) had larger FL2 volumes than FL1 pre-RT (P<.001). Dysfunctional lung became functional in 11 patients (73%) on V SPECT and in 10 patients (67%) on Q SPECT. The dosimetric parameters generated from CT-based anatomical lung had significantly lower values in FL1 than FL2, with a median reduction in the volume of lung receiving a dose of at least 20 Gy (V{sub 20}) of 3%, 5.6%, and mean lung dose of 0.95 and 1.55 on V and Q SPECT respectively. Conclusions: Regional ventilation and perfusion function improve significantly during RT in centrally located NSCLC. Lung dosimetry values vary notably between different definitions of functional lung.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

NASA NTRS1993

New techniques for reversible compression of multispectral data

While spatial correlations are adequately exploited by standard lossless image compression techniques, little success has been attained in exploiting spectral correlations when dealing with multispectral image data. In this paper, we present some new lossless image compression techniques that capture spectral correlations as well as spatial correlation in a simple and elegant manner. The schemes are based on the notion of a prediction tree, which defines a non-causal prediction model for an image. We present a backward adaptive technique and a forward adaptive technique. We then give a computationally efficient way of approximating the backward adaptive technique. The approximation gives good results and is extremely easy to compute. Simulation results show that for high spectral resolution images, significant savings can be made by using spectral correlations in addition to spatial correlations. Furthermore, the increase in complexity incurred in order to make these gains is minimal.

Memon, N. D.↗

DOE OSTI2024

Coupling Anionic Oxygen Redox with Selenium for Stable High‐Voltage Sodium Layered Oxide Cathodes

Utilizing anion redox reaction is crucial for developing the next generation of high-energy density, low-cost sodium-ion batteries. However, the irreversible oxygen redox reaction in Na-ion layered cathodes, which leads to voltage fading and reduced overall lifespan, has hindered their practical application. In this study, selenium is incorporated as a synergistic redox active center of oxygen to improve the stability of Na-ion cathodes. The redesigned cathode maintains stable voltage by demonstrating reversible oxygen redox while significantly suppressing the redox activity of manganese. The anionic redox contribution capacity of the selenium-doped Na 0.6 Li 0.2 Mn 0.8 O 2 cathode remains as high as 84% after 50 cycles, while the pristine Na 0.6 Li 0.2 Mn 0.8 O 2 cathode experiences a reduction to 39% of its initial capacity. The X-ray photoelectron spectroscopy data and computational analysis further revealed that selenium doping participates in redox as Se +4/5 which stabilizes the charged state and increases the energy step for O─O dimerization, thus improving the stability and lifespan of Na 0.6 Li 0.2 Mn 0.8 O 2 cathodes. In conclusion, the findings highlight the potential of redox coupling design to address the issue of voltage fade caused by irreversible anionic redox.

25 ENERGY STORAGE↗

DOE OSTI2026

Photoresponsive Rotaxanes Switch Lipid Bilayer Neuromorphic Behavior with Light

A rotaxane consisting of a macrocycle ring with two azobenzene units mechanically interlocked onto a bolaamphiphilic axle was incorporated into droplet interface bilayers (DIBs). The azobenzene groups on the ring underwent quasi-reversible, photoisomerization-induced cycling between 1-E and 1-Z configurations when irradiated with 370 and 467 nm light, respectively, enabling programmable access to different history-dependent electrical behaviors from the same membrane. In the 1-E configuration, bilayers exhibited type-IIactive memristance that coincided with increasingly elevated ionic conduction, associated with progressively enhanced bilayer permeability during voltage cycling. In the 1-Z configuration, bilayers displayed type-I, passive memcapacitive behavior, reflecting tighter lipid packing and reduced ionic permeability. Photoswitching also yielded a nonvolatile, photoresponsive memcapacitor that could be modulated repetitively with negligible loss, likely via reversible changes in membrane thickness. Concurrent ohmic leakage currents across the membrane were less than 0.3%. These results agree with previous studies with increasing membrane permeability using photoswitchable rotaxanes and provide new insights into the coupling between volatile and nonvolatile memcapacitance during photoisomerization. More broadly, they demonstrate a new strategy for the manipulation of neuromorphic behaviors in soft materials using light, with implications for brain-inspired computation and sensing.

droplet interface bilayers↗

DOE OSTI2023

Spy in the GPU-box: Covert and Side Channel Attacks on Multi-GPU System

The deep learning revolution has been enabled in large part by GPUs, and more recently accelerators, which make it possible to carry out computationally demanding training and inference in acceptable times. As the size of machine learning networks and workloads continues to increase, multi-GPU machines have emerged as an important platform offered on High Performance Computing and cloud data centers. Since these machines are shared among multiple users, it becomes increasingly important to protect applications against potential attacks. In this paper, we explore the vulnerability of Nvidia's DGX multi-GPU machines to covert and side channel attacks. These machines consist of a number of discrete GPUs that are interconnected through a combination of custom interconnect (NVLink) and PCIe connections. We reverse engineer the interconnected cache hierarchy and show that it is possible for an attacker on one GPU to cause contention on the L2 cache of another GPU. We use this observation to first develop a covert channel attack across two GPUs, achieving the best bandwidth of around 4 MB/s. We also develop a prime and probe attack on a remote GPU allowing an attacker to recover the cache access pattern of another workload. This access pattern can be used in any number of side channel attacks: we demonstrate a proof of concept attack that fingerprints the application running on the remote GPU, with high accuracy. We also develop a proof of concept attack to extract hyperparameters of a machine learning workload. Our work establishes for the first time the vulnerability of these machines to microarchitectural attacks and can guide future research to improve their security.

Dutta, Sankha↗

DOE OSTI2019

Evaluation of Acoustic Channel Capacity for Complex Piping Topology

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2019

Final Report for Transmission of Information by Acoustic Communication along Metal Pathways in Nuclear Facilities

Transmission of information using elastic waves on existing metallic pipes provides an alternative communication option for a nuclear facility. The advantages of this approach consist of transmitting information through barriers, such as the containment building wall. A viable candidate for acoustic communication channel is a chemical volume control system (CVCS) stainless steel pipe, which penetrates through the containment building wall. A laboratory bench-scale system consisting of a nuclear grade CVCS-like pipe and ultrasonic transducers was assembled for a preliminary communication system analysis. Because of low bandwidth and spectral dispersion of ultrasonic transducers, on off keying (OOK) protocol was chosen for data communication. Operation of the acoustic communication during normal and post-accident conditions involves data transmission on the CVCS pipe, which according to literature guidelines, will be at 50°C and 150°C, respectively. Preliminary communication system development was performed using paintbrush PZT transducers capable of continuous operation at temperatures up to 80°C. Next, high-temperature LiNbO 3 acoustic transducers, originally developed for flow metering applications of EBRII reactor at Argonne, were utilized to develop a high-temperature acoustic communication system. A Gaussian pulse shape filter was introduced to suppress ringing and thus reduce inter-symbol interference. This resulted in enhancement of data transmission bitrate. Heating tapes, temperature controllers, and thermal insulation were installed on the laboratory pipe to study acoustic communication at elevated temperature. This included development of custom insertion heating elements for controlling temperature at pipe ends. Demonstrations of communication at high temperature included transmission of images and text files. Main achievements thus far include demonstration of transmission of 90KB image at the bitrate of 10Kbps with bit error rate smaller than 10 -3 across the pipe heated to 50°C and 150°C. Final phase of project development involved investigating strategies for deployment of the ultrasonic communication system in a representative environment. Mechanisms Engineering Test Loop (METL) facility at ANL was chosen as a viable candidate for communication system demonstration. Analysis of piping manifolds at METL has shown that in a typical scenario, ultrasonic signals have to be transmitted over straight and bent piping sections. This report describes preliminary studies on transmission of information with ultrasonic shear waves across bent piping sections. A test article consisting of a stainless steel pipe bent at 90° was developed for laboratory analysis of ultrasonic signal transmission. The bent piping test article developed by welding two straight pipes to an elbow. For consistency, the diameter and wall thickness of the bent pipe is the same as that of the straight pipe utilized in prior studies. COMSOL computer simulations were performed to study ultrasonic refracted shear wave coupling and transmission across bent piping. Preliminary evaluations of ultrasonic shear wave transmission were conducted with piezoelectric (PZT) and LiNbO 3 transducers, with the signals recorded with previously developed LabView interface. Preliminary results have shown significant dispersion of transmitted signals, most likely due to reflections and scattering by piping welds. Signal distortion compensation algorithm based on time reversal modulation (TRM) was initially investigated for transmission through solids, and demonstrated for signals transmitted over the bent pipe. Using TRM, it was shown that communication bitrates of 100Kbps were achievable, in principle.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

NASA NTRS2018

Hybrid Reduced Order Model for N2-N2 Interactions for Application to Dissociation and Energy Transfer Processes

Recent work in the aerothermodynamics community has focused on the development of reduced-order models for thermo-chemical non-equilibrium which avoid the restrictive assumptions of multi-temperature models and the prohibitive cost associated with State-to-State (StS) models. In the present work, this is accomplished by lumping energy states together and assuming groups of states are roughly in equilibrium. As a result, the non-equilibrium behavior of a gas can be captured at a reduced computational cost from a full StS simulation. In this work, we present a hybrid grouping model for studying energy transfer and dissociation in a mixture of nitrogen molecules due to N2-N2 reactions. This is accomplished by making use of a grouping strategy informed by data from the N2-N StS kinetic data. However, due to the massive computational cost associated with constructing StS data for the N2-N2 system, the kinetic data for the hybrid grouping model are calculated using the quasi-classical trajectory (QCT) method by sampling states for trajectory within the groups. This general framework is called the Maximum-Entropy Quasi-Classical Trajectory (ME-QCT) method. The primary challenge associated with this method is that rates for reverse grouped reactions cannot be obtained through detailed balance at a group level, due to the variation of group internal temperatures. To construct the full model for N2-N2 grouped kinetics using the ME-QCT method, detailed balance is invoked at the microscopic level, allowing for the calculation of the full kinetic data from QCT. Results will be presented using the full ME-QCT model for the N2-N2 system in an isothermal and isochoric reactor simulation. In addition, simple CFD test cases for a one-dimensional standing shock and for a quasi-one-dimensional nozzle will be used for demonstration of the ME-QCT method. This method allows for the calculation of non-equilibrium behavior for the N2-N2 system without the prohibitive cost of a full StS simulation. Moreover, it enables the construction of a unified model for the dissociating and recombining non-equilibrium flows.

nonequilibrium modeling↗

DOE OSTI2020

TBASS: A Robust Adaptation of Bayesian Adaptive Spline Surfaces

The R package TBASS is an extension of the BASS package created by Francom and Sansó (2019). The package is used to fit a Bayesian multivariate adaptive spline to a dataset that either follows a Student’s t-distribution or has outliers. Much of the framework for TBASS is adapted from the concepts of Bayesian Multivariate Adaptive Regression Splines (BMARS), specifically the work done by Denison, Mallick, and Smith (1998). The spline function is fit using a Reversible-Jump Markov Chain Monte Carlo algorithm,. By including this more robust generalization, a dataset with outliers can be accurately fit using the BMARS model, without the possibility of overfitting or variance inflation.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS2004

Fast Query-Optimized Kernel-Machine Classification

A recently developed algorithm performs kernel-machine classification via incremental approximate nearest support vectors. The algorithm implements support-vector machines (SVMs) at speeds 10 to 100 times those attainable by use of conventional SVM algorithms. The algorithm offers potential benefits for classification of images, recognition of speech, recognition of handwriting, and diverse other applications in which there are requirements to discern patterns in large sets of data. SVMs constitute a subset of kernel machines (KMs), which have become popular as models for machine learning and, more specifically, for automated classification of input data on the basis of labeled training data. While similar in many ways to k-nearest-neighbors (k-NN) models and artificial neural networks (ANNs), SVMs tend to be more accurate. Using representations that scale only linearly in the numbers of training examples, while exploring nonlinear (kernelized) feature spaces that are exponentially larger than the original input dimensionality, KMs elegantly and practically overcome the classic curse of dimensionality. However, the price that one must pay for the power of KMs is that query-time complexity scales linearly with the number of training examples, making KMs often orders of magnitude more computationally expensive than are ANNs, decision trees, and other popular machine learning alternatives. The present algorithm treats an SVM classifier as a special form of a k-NN. The algorithm is based partly on an empirical observation that one can often achieve the same classification as that of an exact KM by using only small fraction of the nearest support vectors (SVs) of a query. The exact KM output is a weighted sum over the kernel values between the query and the SVs. In this algorithm, the KM output is approximated with a k-NN classifier, the output of which is a weighted sum only over the kernel values involving k selected SVs. Before query time, there are gathered statistics about how misleading the output of the k-NN model can be, relative to the outputs of the exact KM for a representative set of examples, for each possible k from 1 to the total number of SVs. From these statistics, there are derived upper and lower thresholds for each step k. These thresholds identify output levels for which the particular variant of the k-NN model already leans so strongly positively or negatively that a reversal in sign is unlikely, given the weaker SV neighbors still remaining. At query time, the partial output of each query is incrementally updated, stopping as soon as it exceeds the predetermined statistical thresholds of the current step. For an easy query, stopping can occur as early as step k = 1. For more difficult queries, stopping might not occur until nearly all SVs are touched. A key empirical observation is that this approach can tolerate very approximate nearest-neighbor orderings. In experiments, SVs and queries were projected to a subspace comprising the top few principal- component dimensions and neighbor orderings were computed in that subspace. This approach ensured that the overhead of the nearest-neighbor computations was insignificant, relative to that of the exact KM computation.

Mazzoni, Dominic↗

DOE OSTI2014

Investigation of the Redox Chemistry of Anthraquinone Derivatives Using Density Functional Theory

Application of density functional calculations to compute electrochemical properties such as redox windows, effect of substitution by electron donating and electron withdrawing groups on redox windows, and solvation free energies for ~50 anthraquinone (AQ) derivatives are presented because of their potential as anolytes in all-organic redox flow batteries. Computations suggest that lithium ions can increase (by ~0.4 V) the reduction potential of anthraquinone due to the lithium ion pairing by forming a Lewis base-Lewis acid complex. To design new redox active species, the substitution by electron donating groups are essential to improve the reduction window of AQ with adequate oxidative stability. For instance, a complete methylation of AQ can improve its reduction window by ~0.4 V. The quantum chemical studies of the ~50 AQ derivatives are used to derive a relationship that connects the computed LUMO energy and the reduction potential that can be applied as a descriptor for screening thousands of AQ derivatives. Our computations also suggest that incorporating oxy-methyl dioxolane substituents in the AQ framework can increase its interaction with non-aqueous solvent and improve its solubility. Thermochemical calculations for likely bond breaking decomposition reactions of un-substituted AQ anions suggest that the dianions are relatively stable in the solution. These studies provide ideal platform to perform further combined experimental and theoretical studies to understand the electrochemical reversibility and solubility of new quinone molecules as energy storage materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2019

Single-layer tensor network study of the Heisenberg model with chiral interactions on a kagome lattice

Here, we study the antiferromagnetic kagome Heisenberg model with an additional scalar-chiral interaction by using the infinite projected entangled-pair state (iPEPS) ansatz. We discuss in detail the implementation of the optimization algorithm in the framework of the single-layer tensor network based on the corner transfer matrix technique. Our benchmark based on the full-update algorithm shows that the single-layer algorithm is stable, which leads to the same level of accuracy as the double-layer ansatz but with much less computation time. We further apply this algorithm to study the nature of the kagome Heisenberg model with a scalar-chiral interaction by computing the bond dimension scaling of magnetization, bond energy difference, chiral order parameter, and correlation length. In particular, we find that for strong chiral coupling the correlation length, which is extracted from the transfer matrix, saturates to a finite value for large bond dimension, representing a gapped spin-liquid state. Further comparison with density matrix renormalization group results supports that our iPEPS faithfully represents the time-reversal symmetry breaking chiral state. Our iPEPS simulation results shed new light on constructing PEPS for describing gapped chiral topological states.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

NASA NTRS2005

Fully Associative, Nonisothermal, Potential-Based Unified Viscoplastic Model for Titanium-Based Matrices

A number of titanium matrix composite (TMC) systems are currently being investigated for high-temperature air frame and propulsion system applications. As a result, numerous computational methodologies for predicting both deformation and life for this class of materials are under development. An integral part of these methodologies is an accurate and computationally efficient constitutive model for the metallic matrix constituent. Furthermore, because these systems are designed to operate at elevated temperatures, the required constitutive models must account for both time-dependent and time-independent deformations. To accomplish this, the NASA Lewis Research Center is employing a recently developed, complete, potential-based framework. This framework, which utilizes internal state variables, was put forth for the derivation of reversible and irreversible constitutive equations. The framework, and consequently the resulting constitutive model, is termed complete because the existence of the total (integrated) form of the Gibbs complementary free energy and complementary dissipation potentials are assumed a priori. The specific forms selected here for both the Gibbs and complementary dissipation potentials result in a fully associative, multiaxial, nonisothermal, unified viscoplastic model with nonlinear kinematic hardening. This model constitutes one of many models in the Generalized Viscoplasticity with Potential Structure (GVIPS) class of inelastic constitutive equations.

Source record↗

DOE OSTI2025

Dehydroxylation kinetics of kaolinite and montmorillonite examined using isoconversional methods

The use of calcined clays as supplementary cementitious materials (SCMs) in concrete is a promising strategy towards decarbonizing the cement and concrete industry. This is especially relevant considering the ever-increasing demand for concrete. Comprehensive understanding of the kinetics of calcination is essential towards maximizing the potential reactivity of clay minerals while ensuring energy efficiency. In this study, the kinetics of the dehydroxylation of kaolinite and montmorillonite are investigated under non-isothermal conditions at constant heating rate. Activation energies ( E a ) are determined via Friedman differential and advanced Vyazovkin incremental methods over the isoconversional range; these are devoid of computational approximations, thus allowing kinetic analysis without assuming a specific reaction model. Kinetic equations—in the differential form as well as a combination of differential and integral forms are compared against the experimentally determined reaction models to identify the most probable dehydroxylation mechanism for kaolinite and montmorillonite. A reaction order mechanism is established for dehydroxylation of kaolinite, while montmorillonite is noted to undergo dehydroxylation via a single-step reversible diffusion-controlled process. Kinetic triplet—comprising activation energy, reaction model and pre-exponential factor—is used to predict isothermal calcination conditions, which is further verified using analytical techniques. Heat release rates of clay-portlandite blends from isothermal calorimetry are used within a thermodynamic framework to quantify reactivity of the calcined clays. Here, the study demonstrates a general approach based on isoconversional methods to predict calcination conditions for different clays that can be used in efficient and optimized production of blended cements or SCMs.

36 MATERIALS SCIENCE↗

DOE OSTI2019

Relaxation Dynamics of Hydrated Thymine, Thymidine, and Thymidine Monophosphate Probed by Liquid Jet Time-Resolved Photoelectron Spectroscopy

In this work, the relaxation dynamics of thymine and its derivatives thymidine and thymidine monophosphate are studied using time-resolved photoelectron spectroscopy applied to a water microjet. Two absorption bands are studied; the first is a bright ππ* state which is populated using tunable-ultraviolet light in the range 4.74–5.17 eV and probed using a 6.20 eV probe pulse. By reversing the order of these pulses, a band containing multiple ππ* states is populated by the 6.20 eV pulse and the lower energy pulse serves as the probe. The lower lying ππ* state is found to decay in ~400 fs in both thymine and thymidine independent of pump photon energy, while thymidine monophosphate decays vary from 670 to 840 fs with some pump energy dependence. The application of a computational quantum mechanical/molecular mechanical scheme at the XMS-CASPT2//CASSCF/AMBER level of theory suggests that conformational differences existing between thymidine and thymidine monophosphate in solution account for this difference. The higher lying ππ* band is found to decay in ~600 fs in all three cases, but it is only able to be characterized when the 5.17 eV probe pulse is used. Notably, no long-lived signal from an nπ* state can be identified in either experiment on any of the three molecules.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2015

Elucidating the structure of the magnesium aluminum chloride complex electrolyte for magnesium-ion batteries

Non-aqueous Mg-ion batteries offer a promising way to overcome safety, costs, and energy density limitations of state-of-the-art Li-ion battery technology. Here, we present a rigorous analysis of the magnesium aluminum chloride complex (MACC) in tetrahydrofuran (THF), one of the few electrolytes that can reversibly plate and strip Mg. We use ab initio calculations and classical molecular dynamics simulations to interrogate the MACC electrolyte composition with the goal of addressing two urgent questions that have puzzled battery researchers: (i) the functional species of the electrolyte, and (ii) the complex equilibria regulating the MACC speciation after prolonged electrochemical cycling, a process termed as conditioning, and after prolonged inactivity, a process called aging. A general computational strategy to untangle the complex structure of electrolytes, ionic liquids and other liquid media is presented. The analysis of formation energies and grand-potential phase diagrams of Mg–Al–Cl–THF suggests that the MACC electrolyte bears a simple chemical structure with few simple constituents, namely the electro-active species MgCl + and AlCl 4 - in equilibrium with MgCl 2 and AlCl 3 . Knowledge of the stable species of the MACC electrolyte allows us to determine the most important equilibria occurring during electrochemical cycling. We observe that Al deposition is always preferred to Mg deposition, explaining why freshly synthesized MACC cannot operate and needs to undergo preparatory conditioning. Similarly, we suggest that aluminum displacement and depletion from the solution upon electrolyte resting (along with continuous MgCl 2 regeneration) represents one of the causes of electrolyte aging. Finally, we compute the NMR shifts from shielding tensors of selected molecules and ions providing fingerprints to guide future experimental investigations.

25 ENERGY STORAGE↗