Search NASA⌕ Search

SEARCH · Search NASA

Results for “solvates”

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 397 records · Page 22
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Localized High Concentration Electrolyte and Its Effects on Polysulfide Structure in Solution

Lithium–sulfur batteries (LiSB) are a promising next-generation lithium energy storage technology that offers a multifold improvement over the traditional lithium-ion battery. However, the LiSB still faces the unresolved issue of the polysulfide shuttle effect. This phenomenon arises from the dissolution of sulfur intermediate reduction products into the electrolyte, which then causes a cascade of issues throughout the battery. Many mitigation strategies have been proposed to counteract this effect including the use of novel electrolyte compositions. Recently, there has been increased focus on the use of hydrofluorinated ethers as a major constituent of an electrolyte for lithium batteries. Previous studies have reported that the presence of these species create localized high concentration electrolytes (LHCE) which can have many advantages resulting in improved battery performance. Here, we report on how the inclusion of a hydrofluorinated ether BTFE (bis(2,2,2-trifluoroethyl) ether) modifies the general properties of a LiSB electrolyte and the structure of dissolved polysulfide species in the electrolyte. Here, we found that the inclusion of BTFE does not modify the primary solvation structure of Li + directly but actively participates in secondary solvation shells. With a high concentration of BTFE, the LHCE formation is observed by the presence of clusters of non-BTFE molecules. The structure of polysulfide species in solution was modified by the BTFE in the same way. Much of the primary solvation structure was kept, but the presence of BTFE increased polysulfide–polysulfide clustering. These results indicate that the polysulfide solubility will be limited due to the promotion of clustering.

25 ENERGY STORAGE↗
DOE OSTI2025

Epoxide Alcoholysis over M-BEA Zeolites: Effects of Alcohol Chain Length on Rates and Regioselectivities

The structures of nucleophilic reactants affect their coordination behavior among solvent molecules and kinetics of reactions with surface intermediates within the confines of fluidfilled pores of zeolites and other microporous materials. Consequently, rates and regioselectivities of diverse chemistries may depend sensitively on nucleophile identity in manners not observed for classic fluid phase reactions. Here, we examine the impact of varying the primary alcohol (ROH) chain length on the kinetics of 1,2-epoxybutane (C 4 H 8 O) ring-opening within Brønsted (Al-BEA) and Lewis acid (Zr-BEA) zeolites. Turnover rates increase by factors of ∼6 (Al-BEA) and 4-fold (Zr-BEA) between methanol and 1-hexanol, yet the reaction mechanisms remain comparable. Despite modest rate differences, apparent activation enthalpies calculated from rates and activities of solvated reactants decrease linearly by 12 (Al-BEA) to 33 kJ mol −1 (Zr- BEA) with increased proton affinity, which suggests bond formation energies for the nucleophile strongly influence rate increases. The molecular interpretation of these trends demonstrates, however, that the solvation of ring-opening transition states by zeolite pore structures and solvent molecules also governs rates. The impact of local solvating interactions appears most directly as changes in regioselectivities, which tend to enhance terminal alcohol formation with increasing ROH chain length. Regioselectivities largely do not vary with differences in fluid composition for a given ROH. The addition of H 2 O increases the number of hydrogen bonds among reactive species, and trends in regioselectivities imply that the decreased hydrogen bonding ability of longer chain ROH, and not the nucleophile strength or steric bulk, determines the regioselectivities of the resulting products. This work provides direct experimental evidence that nucleophilicity and hydrogen bonding influence reaction barriers and regioselectivities in zeolitecatalyzed epoxide ring-opening, offering pathways to better control reaction kinetics.

acidic zeolites↗
DOE OSTI2025

Locally Confined Polysulfide-Reactive Electrolytes for Shuttle-Free Sodium–Sulfur Batteries

Sodium-sulfur batteries promise high-energy-density and sustainable electrochemical energy storage but suffer from uncontrolled polysulfides dissolution and high sodium reactivity. These challenges fundamentally originate from poor electrolyte-electrode compatibility. Current electrolyte research inadequately addresses the trade-off between minimal polysulfides solvation and stabilizing sodium interfaces. Here, we present a locally-confined polysulfide-reactive electrolyte strategy that mediates the polysulfide dissolution dynamics and sodium stability by leveraging an electrophilic solvating species with a localized high-concentration electrolyte. This design enables shuttle-free cell operation by synergistically restricting the global solvating power of the electrolyte through intermolecular interactions and locally scavenging sparingly dissolved polysulfides via electrolyte electrophilicity. The precisely confined surface reaction facilitates a protective cathode-electrolyte interface, realizing a quasi-solid-state sulfur conversion in our liquid ether-based electrolyte, which crucially avoids crossover-induced catastrophic sodium-metal degradation. The proposed electrolyte demonstrates long-term cycling of high-mass-loading sulfur cathodes (> 3 mg S cm −2 with commercial carbon host and 70 wt% sulfur content), which afford 710 mA h g −1 over 400 cycles in coin cells and steady pouch cell operation over 180 cycles. Furthermore, this work establishes a scalable electrolyte design protocol that regulates the reaction chemistry of highly reactive electrodes, offering a pathway toward sustainable renewable energy storage.

25 ENERGY STORAGE↗
DOE OSTI2024

Development and application of hybrid AIMD/cDFT simulations for atomic-to-mesoscale chemistry

Many important chemical processes involve reactivity and dynamics in complex solutions. Gaining a fundamental understanding of these reaction mechanisms is a challenging goal that requires advanced computational and experimental approaches. However, important techniques such as molecular simulation have limitations in terms of scales of time, length, and system complexity. Furthermore, among the currently available solvation models, there are very few designed to describe the interaction between the molecular scale and the mesoscale. To help address this challenge, here, we establish a novel hybrid approach that couples first-principles plane-wave density functional theory with classical density functional theory (cDFT). In this approach, a region of interest described by ab initio molecular dynamics (AIMD) interacts with the surrounding medium described using cDFT to arrive at a self-consistent ground state. cDFT is a robust but efficient mesoscopic approach to accurate thermodynamics of bulk electrolyte solutions over a wide concentration range (up to 2M concentrations). Benchmarking against commonly used continuum models of solvation, such as SMD, as well as experiments, demonstrates that our hybrid AIMD–cDFT method is able to produce reasonable solvation energies for a variety of molecules and ions. With this model, we also examined the solvent effects on a prototype S N 2 reaction of the nucleophilic attack of a chloride ion on methyl chloride in the solution. The resulting reaction pathway profile and the solution phase barrier agree well with experiment, showing that our AIMD/cDFT hybrid approach can provide insight into the specific role of the solvent on the reaction coordinate.

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

Electric-Field-Induced Spatially Dynamic Heterogeneity of Solvent Motion and Cation Transference in Electrolytes

While electric fields primarily result in migration of charged species in electrolytic solutions, the solutions are dynamically heterogeneous. Solvent molecules within the solvation shells of the cation will be dragged by the field while free solvent molecules will not. Here we combine electrophoretic NMR measurements of ion and solvent velocities under applied electric fields with molecular dynamics simulations to interrogate different solvation motifs in a model liquid electrolyte. Measured values of the cation transference number $(t^0_+)$ agree quantitatively with simulation-based predictions over a range of electrolyte concentrations. Solvent-cation interactions strongly influence the concentration-dependent behavior of $t^0_+$. We identify a critical concentration at which most of the solvent molecules lie within solvation shells of the cations. The dynamic heterogeneity of solvent molecules is minimized at this concentration where $t^0_+$ is approximately equal to zero.

42 ENGINEERING↗
DOE OSTI2022

Non‐Flammable Ester Electrolyte with Boosted Stability Against Li for High‐Performance Li Metal Batteries

Abstract In traditional non‐flammable electrolytes a trade‐off always exists between non‐flammability and battery performance. Previous research focused on reducing free solvents and forming anion‐derived solid‐electrolyte interphase. However, the contribution of solvated anions in boosting the stability of electrolyte has been overlooked. Here, we resolve this via introducing anions into Li + solvation sheaths using anions with similar Gutmann donor number (DN) to that of solvents. Taking trimethyl phosphate fire‐retardant (DN=23.0 kcal mol −1 ) and NO 3 − (DN=22.2 kcal mol −1 ) as an example, NO 3 − is readily involved in the Li + solvation sheath and reduces the polarity of solvent. This results in boosted stability of electrolyte against Li. The developed non‐flammable electrolyte has low viscosity, high ionic conductivity and is low cost. The reversibility of Li‐Cu cell was improved to 99.49 % and the lifespan of practical LMBs was extended by >100 %.

Wang, Zhijie↗
DOE OSTI2022

Non‐Flammable Ester Electrolyte with Boosted Stability Against Li for High‐Performance Li Metal Batteries

Abstract In traditional non‐flammable electrolytes a trade‐off always exists between non‐flammability and battery performance. Previous research focused on reducing free solvents and forming anion‐derived solid‐electrolyte interphase. However, the contribution of solvated anions in boosting the stability of electrolyte has been overlooked. Here, we resolve this via introducing anions into Li + solvation sheaths using anions with similar Gutmann donor number (DN) to that of solvents. Taking trimethyl phosphate fire‐retardant (DN=23.0 kcal mol −1 ) and NO 3 − (DN=22.2 kcal mol −1 ) as an example, NO 3 − is readily involved in the Li + solvation sheath and reduces the polarity of solvent. This results in boosted stability of electrolyte against Li. The developed non‐flammable electrolyte has low viscosity, high ionic conductivity and is low cost. The reversibility of Li‐Cu cell was improved to 99.49 % and the lifespan of practical LMBs was extended by >100 %.

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

Comparative Study of Solvatomorphs of Stryker's Reagent Using MicroED and Quantum Mechanics

Abstract The atomic position of hydrogen atoms in metal hydrides has been a long‐standing structural question in inorganic chemistry given that hydride delivery is integral to diverse chemical reactions. Microcrystal electron diffraction (microED), with it's increased sensitivity toward hydrogen atoms relative to X‐ray diffraction, offers a potential path to addressing this challenge. Herein, the first microED study of Stryker's reagent is reported, resulting in the structure of a new benzene solvate. Improved accuracy for hydrogen atom positions was obtained via a quantum crystallography (QCr) approach, Hirshfeld atom refinement (HAR). Structural and topological analysis supports edge bridging hydrides in the microED structure of a THF solvate form, consistent with previous diffraction studies. Interestingly, analysis of a new benzene solvate, discovered in this study, is consistent with mixed edge‐ and face‐bridging hydrides.

Jha, Kunal K. [Division of Chemistry &amp, Chemica↗
DOE OSTI2025

Mechanistic Insights into the Demethylation of Lignin‐Derived Structures Using Protic Ionic Liquids: A Density Functional Theory Study

Lignin valorization is restricted by the stability of its methoxy groups, creating a critical need for efficient demethylation strategies. Here, in this study, density functional theory (DFT) is employed to dissect the mechanistic pathways of demethylation in lignin model compounds, guaiacol and syringol, using protic ionic liquids (PILs) that act as both solvent and catalyst. Conductor‐like Screening Model for Real Solvents (COSMO‐RS) analysis identifies monoethanolammonium acetate ([MEOA][Ace]) as the most promising medium, attributed to its strong hydrogen bonding network and solvation ability. By integrating implicit and explicit solvation models, it is revealed that an acid‐catalyzed hydrolytic mechanism governs demethylation, with PILs stabilizing crucial transition states and intermediates. Complementary electronic structure evaluations, including highest occupied molecular orbital–lowest unoccupied molecular orbital (HOMO‐LUMO) gap analysis, charge distribution, and electrostatic potential mapping, demonstrate how PILs lower energetic barriers and enhance reactivity. To simulate realistic environments, this study is extended to lignin dimer complexes with varying water content, uncovering how water‐bridged solvation reverses demethylation preference from guaiacyl to syringyl units. This mechanistic shift aligns with experimental observations showing faster S‐unit reactivity in hydrated systems. Together, these findings provide atomic‐level insight into lignin demethylation dynamics and highlight how tuning acid concentration in PILs can accelerate kinetics, enabling a rational pathway toward next‐generation biomass conversion technologies.

COSMO-RS↗
DOE OSTI2025

Designing Advanced Electrolytes for High-Voltage High-Capacity Disordered Rocksalt Cathodes

Lithium (Li)-excess transition metal oxide materials which crystallize in the cation-disordered rock salt (DRX) structure are promising cathodes for realizing low-cost, high-energy-density Li batteries. However, the state-of-the-art electrolytes for Li-ion batteries cannot meet the high-voltage stability requirement for high-voltage DRX cathodes, thus new electrolytes are urgently demanded. It has been reported that the solvation structures and properties of the electrolytes critically influence the performance and stability of the batteries. In this study, the structure–property relationships of various electrolytes with different solvent-to-diluent ratios are systematically investigated through a combination of theoretical calculations and experimental tests and analyses. This approach guides the development of electrolytes with unique solvation structures and characteristics, exhibiting high voltage stability, and enhancing the formation of stable electrode/electrolyte interphases. These electrolytes enable the realization of Li||Li 1.094 Mn 0.676 Ti 0.228 O 2 (LMTO) DRX cells with improved performance compared to the conventional electrolyte. Specifically, Li||LMTO cells with the optimized advanced controlled-solvation electrolyte deliver higher specific capacity and longer cycle life compared to cells with the conventional electrolyte. Additionally, the investigation into the structure–property relationship provides a foundational basis for designing advanced electrolytes, which are crucial for the stable cycling of emerging high-voltage cathodes.

25 ENERGY STORAGE↗
DOE OSTI2021

Insight into the nanostructure of “water in salt” solutions: A SAXS/WAXS study on imide-based lithium salts aqueous solutions

"Water-in-salt " electrolyte (WISE) series have been studied for lithium-ion batteries and supercapacitor applications, but so far, most of the focus has been on the LiTFSI salt-based systems. Herein, we used small-angle X-ray scattering/wide-angle X-ray scattering (SAXS/WAXS) and molecular dynamics (MD) simulation to investigate the solvation structure of a series of lithium salts with four different symmetric anions: (bis(fluoro sulfonyl)imide (FSI), bis(trifluoromethane sulfonyl)imide (TFSI), bis(pentafluoroethane sulfonyl) imide (BETI) and bis(nonafluorobutane sulfonyl)imide (BNTI)), which have similar parent structures but different lengths of the fluorocarbon chains. Two competing structures were found in the four lithium salts aqueous solutions: anion solvated structure and anion network. The anion network plays a crucial role in obtaining a stable and enlarged electrochemical window. The d spacing of the anion solvated structure follows a linear correlation with the number of carbons in the fluorocarbon chains and an exponential correlation with the concentrations. Further, we also observed that the transition from one structure to another is predominantly controlled by the salt volume fraction for all imide-based aqueous solutions. This work provides a systematic study of the atomistic scale structures of a series of imide-based lithium salt aqueous solutions that could extend the knowledge of WIS electrolytes.

25 ENERGY STORAGE↗
DOE OSTI2026

Protonation Dynamics of Confined Ethanol–Water Mixtures in H-ZSM-5 from Machine Learning-Driven Metadynamics

Zeolites are indispensable heterogeneous catalysts in industrial chemical processes, valued for their strong Brønsted acidity, well-defined microporous frameworks, and tunable pore structures. Their catalytic activity arises primarily from Brønsted acid sites (BAS), typically present as bridging hydroxyl groups (Si–OH–Al). Under aqueous reaction conditions, these protons interact dynamically with water and alcohol molecules, leading to complex solvation and protonation behavior within confined pores. In this study, we investigate the protonation equilibrium occurring between ethanol and water at the BAS of acidic zeolites under varying hydration levels, i.e., C2H5OH–(H2O)n, n=1–4. Local structure was analyzed through an adaptive-learning global optimization algorithm, while enhanced sampling molecular dynamics simulations with Well-Tempered Metadynamics (WMetaD) and machine learning interatomic potentials (MLPs) provide free-energy surfaces (FES) at variable hydration levels. The results reveal a strong dependence of proton localization on the degree of hydration. At low hydration (1 water molecule), the proton resides predominantly on ethanol; with 2 water molecules, it shifts toward water, and at higher hydration (3 or more water molecules), it becomes extensively delocalized over the water cluster. These findings underscore the critical role of solvation in modulating acid site behavior and suggest that a minimum of three water molecules is necessary to fully stabilize the proton on water within the zeolite framework. This solvation threshold has significant implications for catalytic processes, particularly in biomass conversion reactions where alcohol protonation is a key step in dehydration mechanisms.

machine learning↗
DOE OSTI2023

Relationship of the Molecular Structure and Transport Properties of Imide-Based Lithium Salts of “Acetonitrile/Water-in-Salt” Electrolytes

“Water-in-salt” (WIS) electrolytes exhibit excellent safety and electrochemical performance. However, they possess high concentrations, relatively low diffusion coefficient, and high viscosity. “Acetonitrile/water in salt” (AWIS) electrolytes can overcome the disadvantages of WIS electrolytes. Under relatively low concentrations, AWIS electrolytes show good electrochemical performance comparable to WIS electrolytes and low conductivity. Herein, we investigate the relationship between the solvation structures and the transport properties using small-angle X-ray scattering and molecular dynamics simulation. We observed two solvation behaviors of AWIS: anions dissolved in acetonitrile forming small acetonitrile/anion clusters and additional water further dissolving the small acetonitrile/anion clusters. The introduction of acetonitrile weakens the water–solute interaction and enhances the cation–anion interaction, which results in an enhanced dynamical slowdown as the concentration increases. In conclusion, this work provides molecular-level understanding of the connection between two-stage solvation structures and transport properties for imide-based lithium salt solutions.

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

Correlated Ion Transport Governed by Dynamic Local Structure in High Concentration and Localized High Concentration Electrolytes

Understanding the dynamics of cluster formation and network percolation provides the mechanistic link between microscopic solvation structure and transport in concentrated electrolytes, including localized high-concentration electrolytes (LHCEs). Although recent studies have shown LHCEs to form micelle-like aggregates at specific compositions, a quantitative understanding of how solvation structures and transport properties depend on salt–solvent–diluent ratios remains limited. Here, we integrate molecular dynamics with Onsager transport analyses to chart the evolution of solvation microstructure and associated ionic transport in LiFSI/DMC/TTE, achieving good agreement with experimental conductivities across composition. We show that micelle formation arises from a dynamic instability of the cation–anion network, and that network percolation is the primary determinant of conductivity. Two compositional thresholds emerge: A critical network concentration (CNC) and a critical micelle concentration (CMC) that delineate transitions from extended percolating networks to micelle-like clusters and then to fragments. These structural transitions rationalize the nonintuitive conductivity decrease at intermediate dilution despite monotonically decreasing viscosity and provide composition-level design rules for LHCEs.

Mohanakrishnan, Rohith Srinivaas [University of Ca↗
DOE OSTI2022

Structure of Sulfuric Acid Solutions Using Pair Distribution Function Analysis

Solvation and mesoscale ordering of sulfuric acid and other strong acid solutions leads to suppressed freezing points and strong rheological changes with varying concentration. While the solid-state structures are well-understood, studies focused on the evolving solvation structure in the solution phase have probed a limited concentration range (~1–6 M). This study applies a total scattering approach in both the wide-angle X-ray scattering (WAXS) and pair distribution function (PDF) regimes to elucidate the evolving solvation structure over its full range of acid concentration (0–18 M). The emergence of a prepeak in the WAXS regime at intermediate concentrations indicates a transition from noninteracting sulfate molecules in the dilute limit to sterically limited sulfates at concentrations near its deep eutectic point. Fits to the PDF data quantify this trend, showing a transition from octahedrally hydrated sulfates up to 6–7 M concentrations, followed by gradual dehydration, and eventually reaching a solution structure similar to that of water-in-salt electrolyte systems at high acid concentrations.

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

Molecular Dynamics and Free Energy Calculations of Dicyclohexano-18-crown-6 Diastereoisomers with Sm 2+ , Eu 2+ , Dy 2+ , Yb 2+ , Cf 2+ , and Three Halide Salts in Tetrahydrofuran and Acetonitrile Using the AMOEBA Force Field

With the continual development of lanthanides (Ln) in current technological devices, an efficient separation process is needed that can recover greater amounts of these rare elements. Dicyclohexano-18-crown-6(DCH18C6) is a crown ether that may be a promising candidate for Ln separation, but additional research is required. As such, molecular dynamics(MD) simulations have been performed on four divalent lanthanide halide salts(Sm 2+ , Eu 2+ , Dy 2+ , and Yb 2+ ) and one divalent actinide halide salt (Cf 2+ ) bound to three diastereoisomers of DCH18C6. Dy 2+ , Yb 2+ , Cf 2+ , DCH18C6, and tetrahydrofuran (THF) solvent were parameterized for the AMOEBA polarizable force field for the first time, whereas existing parameters for Sm 2+ and Eu 2+ were utilized from our previous efforts. A coordination number (CN) of six for Ln 2+ /An 2+ –O solvated in THF indicated that the cations interacted almost entirely with the oxygens of the polyether ring. A CN of one for Ln 2+ /An 2+ -N solvated in acetonitrile for systems containing iodide suggested that theN atom of acetonitrile was competitive with I – for cation interactions. Fluctuation between five and six CNs for Dy 2+ and Yb 2+ suggested that although the cations remained in the polyether ring, the size of the ring may not be an ideal fit as these cations possess comparatively smaller ionic radii. Gibbs binding free energies of Sm 2+ in all DCH18C6 diastereoisomers solvated in THF were calculated. The binding free energy of the cis-syn-cis diastereoisomer was the most favorable, followed by cis-anti-cis, and then trans-anti-trans. Lastly, two major types of conformation were observed for each diastereoisomer that were related to the electrostatic interactions and charge density of the cations.

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

Modeling the Behavior of Complex Aqueous Electrolytes Using Machine Learning Interatomic Potentials: The Case of Sodium Sulfate

Understanding the structure and thermodynamics of solvated ions is essential for advancing applications in electrochemistry, water treatment, and energy storage. While ab initio molecular dynamics methods are highly accurate, they are limited by short accessible time and length scales whereas classical force fields struggle with accuracy. Herein, we explore the structure and thermodynamics of complex monovalent-divalent ion pairs using Na 2 SO 4 (aq) as a case study by applying a machine learning interatomic potential (MLIP) trained on density functional theory (DFT) data. Our MLIP-based approach reproduces key bulk properties such as density and radial distribution functions of water. We provide the hydration structure of the sodium and sulfate ions in the 0.1–2 M concentration range and the one-dimensional and two-dimensional potentials of mean force for the sodium–sulfate ion pairing at the low concentration limit (0.1 M), which are inaccessible to DFT. At low concentrations, the sulfate ion is strongly solvated, leading to the stabilization of solvent-separated ion pairs over contact ion pairs. Minimum energy pathway analysis revealed that coordinating two sodium ions with a sulfate ion is a multistep process whereby the sodium ions coordinate to the sulfate ion sequentially. Finally, we demonstrate that MLIPs allow the study of solvated ions beyond simple monovalent pairs with DFT-level accuracy in their low concentration limit (0.1 M) via statistically converged properties from ns-long simulations.

anions↗
DOE OSTI2020

Vibrational Spectroscopy Reveals Effects of Electron Push–Pull and Solvent Polarity on Electron Delocalization in Radical Anions of Donor–Acceptor Compounds

Herein, the nature of excess electrons has been studied in donor–acceptor (D–A) compounds based on substituted triarylamines and a nitrile-functionalized fluorene, by changing the substituents on the triarylamines and also the solvent polarity. We observed that both electron push-pull capability at the distant location in the amine donor unit and solvation in solvents of varying polarity significantly affect the nitrile $ν(C≡N)$ vibrations of the fluorene acceptor unit in radical anions of these D–A compounds. Quantum calculations show that the push-pull capability translates the position of an excess electron while keeping its width relatively constant. On the other hand, solvation changes both, making an electron more compact in polar media. The current study points to the idea that solvation plays a more significant role in controlling the nature of excess electrons, while synthetic modification that influences electron push-pull capability enables further tuning.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗