Search NASA⌕ Search

SEARCH · Search NASA

Results for “Reaction rate constants”

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

Assessment of the Second-Ionization Potential of Lawrencium: Investigating the End of the Actinide Series with a One-Atom-at-a-Time Gas-Phase Ion Chemistry Technique

Experiments were performed at the Lawrence Berkeley National Laboratory 88-Inch Cyclotron facility to investigate the electron-transfer reduction reaction of dipositive Lr ( Z = 103) with O 2 gas. Ions of 255 Lr were produced in the fusion-evaporation reaction 209 Bi( 48 Ca,2n) 255 Lr and were studied with a novel gas-phase ion chemistry technique. The produced 255 Lr 2+ ions were trapped and O 2 gas was introduced, such that the charge-exchange reaction to reduce 255 Lr 2+ to 255 Lr 1+ was observed and the reaction rate constant was determined to be k = 1.5(7) × 10 -10 cm 3 /mol/s. The observation that this reaction proceeds establishes the lower limit on the second ionization potential of Lr to be 13.3(3) eV. This gives further support that the actinide series terminates with Lr. Additionally, this result can be used to better interpret the situation concerning the placement of Lu and Lr on the periodic table within the current framework of the actinide hypothesis. The success of this experimental approach now identifies unique opportunities for future gas-phase reaction studies on actinide and super heavy elements.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Formation of phenanthrenyl radicals via the reaction of acenaphthyl with acetylene

Ab initio G3(MP2,CC)/B3LYP/6–31G** calculations of the potential energy surface for the reaction of 3-acenaphthyl radical with acetylene were combined with Rice-Ramsperger-Kassel-Marcus Master Equation calculations of the temperature- and pressure-dependent reaction rate constants and product branching ratios with the goal to examine the transformation of an edge five-membered ring into a six-membered one. The reaction mechanism has been shown to contain two stages: (1) the formation of 3-ethynylacenaphthylene + H in the 3-acenaphthyl + C 2 H 2 reaction via acetylene addition to the radical site and an H loss, (2) the formation of the 1-phenanthrenyl radical in the 3-ethynylacenaphthylene + H reaction via H addition to the five-membered ring followed by the ring flip. The results of kinetic modeling show this pathway being faster than the pathway proposed earlier featuring the CCH chemisorption on the six-membered ring without involvement of the five-membered ring. A kinetic test demonstrated the dominance of the former pathway for the formation of pyrene from naphthalene. Here, it is expected that the new mechanism plays a major role in the transformation of five-membered rings on edges of PAH structures.

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

Bipolar Membranes for Ion Management in (Photo)Electrochemical Energy Conversion

Conspectus (Photo)electrochemical energy conversion is important in the development of a carbon-neutral energy economy because it can provide a pathway for mitigating the intermittency of renewable energy sources such as wind and solar. In order to operate efficiently, these technologies, which include photoelectrochemical cells, water and CO 2 electrolyzers, fuel cells, and redox flow batteries, require fast charge transfer kinetics at the electrode/electrolyte interface as well as robust ion management in the electrolyte. In conventional electrolyzers and fuel cells, the electrolyte is strongly acidic or basic and ionic current is carried by H + or OH – ions. In contrast, photoelectrodes and electrocatalysts for water splitting are often studied in buffered solutions. The question of ion balance in these systems led us to analyze the polarization losses due to ion concentration gradients in cells that employed various buffer–membrane combinations. Continuously driving the buffer ions across an ionomer membrane not only lowers the buffer capacity of an aqueous electrolyte but also introduces pH gradients that result in significant energy losses. To address the problem, we and other groups have studied the use of reverse-biased bipolar membranes (BPMs) in (photo)electrolytic cells. BPMs consist of an anion exchange layer (AEL) laminated with a cation exchange layer (CEL) and are usually equipped with a catalytic layer in between to accelerate the water dissociation reaction. At the AEL/CEL interface, water dissociates into protons and hydroxide ions, which replenish those consumed at the cathode and anode. Compared to conventional water electrolyzers with proton/anion exchange membranes (PEM/AEM), BPM electrolyzers provide the unique advantage of continuously operating the cathode and anode under different pH conditions, which is desirable when the two electrode reactions have different pH requirements. BPMs also enable the use of buffered electrolytes at pH values that are optimized for electrode stability and product selectivity in applications such as CO 2 electrolysis. Product crossover losses and CO 2 pumping can be dramatically reduced in BPM-based CO 2 electrolyzers, relative to conventional alkaline membranes, by electrostatic repulsion (of anionic products) and electroosmotic drag (for neutral products). BPM-based gas fed CO 2 electrolyzers can achieve high current density, but they suffer from low Faradaic efficiency (FE) due to the acidic local environment of the CEL. This problem can be mitigated by adding an aqueous buffering layer or by creating a weak acid cation exchange film on the CEL face of the membrane. The use of BPMs in fuel cells and redox flow batteries offers some interesting advantages. Configurations with both reverse and forward bias have been studied, but forward bias has been favored due to material compatibility, reaction kinetics, and thermodynamic considerations. The net reaction at the AEL/CEL interface is the acid–base neutralization reaction, which has a high inherent reaction rate constant, but in the BPM is limited to a nanometric space-charge layer and requires efficient catalysis to achieve high current density. Understanding the mechanism of the acid–base neutralization and the opposite process, the water dissociation reaction, will be essential for improving the performance of forward-biased BPMs. In conclusion, this Account reviews our current understanding of the working mechanisms of BPMs and discusses how we can use them to effectively manage ions for various (photo)electrochemical applications.

Yan, Zhifei↗
DOE OSTI2023

Hydroxyl radical‐driven transformations of bisphenol A and 2,4‐dinitroanisole: Experimental and computational analysis

Abstract This study used experimental and computational analysis to investigate the advanced oxidation of bisphenol A (BPA) and 2,4‐dinitroanisole (DNAN). The pseudo first‐order reaction rate constants depended on the molar peroxide ratio and were between 0.13 and 0.28 min −1 for BPA and between 0.018 and 0.032 min −1 for DNAN. The kinetic differences appear to be due in part to the energy requirements for oxidation, which depended on the reaction mechanism but were typically lower for BPA than they were for DNAN. Density functional theory (DFT) was used to develop transformation pathways that included experimentally‐detected byproducts. The most energetically favored pathway for BPA oxidation begins with the formation of hydroxylated derivatives, while for DNAN, the most energetically favorable degradation pathway begins with the substitution of the methoxy group. Overall, these findings demonstrate the power of combining experimental and computational tools to reveal transformation mechanisms during water treatment. Practitioner Points Advanced oxidation transformations for two emerging water pollutants, bisphenol A and dinitroanisole, was investigated. The observed reaction kinetics depended on molar peroxide ratio in a manner that is in keeping with previous findings. Density functional theory‐based analysis revealed reaction energy requirements and degradation pathways.

Engineering↗
DOE OSTI2025

Developing the Science Basis for Understanding Polymer Encapsulant Degradation Mechanisms: DuraMAT 2.0 Final Project Report

Polymeric encapsulants are essential materials in photovoltaic modules, protecting sensitive electronics from the environment while providing mechanical integrity to the multilayered assembly. However, these polymeric materials are susceptible to degradation processes driven by the ingress of environmental species, ultraviolet radiation, thermal stresses, and mechanical loading. In this study, we employ a combined atomistic simulation and accelerated aging experimental approach to study the molecular-scale mechanisms of encapsulant degradation. Classical molecular dynamics simulations quantify the diffusion of environmental and degradation species through the polymer matrix, producing composition-specific diffusion coefficients. Reactive simulations characterize activation energy barriers and reaction rate constants for key chemical pathways. In parallel, thermal-desorption analyses coupled with mass spectrometry monitor the emergence and concentration profiles of degradation products under controlled stressor conditions. By integrating simulation and experiment, we establish quantitative correlations between polymer composition, species diffusivity, and chemical reactivity. We anticipate that these relations and quantitative values could serve as high-fidelity inputs to reaction-diffusion models, enabling physics-informed lifetime predictions and guiding the design of more durable encapsulant materials for solar energy applications.

36 MATERIALS SCIENCE↗
DOE OSTI2023

An exploratory steady–state redox model of photosynthetic linear electron transport for use in complete modelling of photosynthesis for broad applications

A photochemical model of photosynthetic electron transport (PET) is needed to integrate photophysics, photochemistry, and biochemistry to determine redox conditions of electron carriers and enzymes for plant stress assessment and mechanistically link sun-induced chlorophyll fluorescence to carbon assimilation for remotely sensing photosynthesis. Towards this goal, we derived photochemical equations governing the states and redox reactions of complexes and electron carriers along the PET chain. These equations allow the redox conditions of the mobile plastoquinone pool and the cytochrome b 6 f complex (Cyt) to be inferred with typical fluorometry. The equations agreed well with fluorometry measurements from diverse C 3 /C 4 species across environments in the relationship between the PET rate and fraction of open photosystem II reaction centres. We found the oxidation of plastoquinol by Cyt is the bottleneck of PET, and genetically improving the oxidation of plastoquinol by Cyt may enhance the efficiency of PET and photosynthesis across species. Redox reactions and photochemical and biochemical interactions are highly redundant in their complex controls of PET. Although individual reaction rate constants cannot be resolved, they appear in parameter groups which can be collectively inferred with fluorometry measurements for broad applications. The new photochemical model developed enables advances in different fronts of photosynthesis research.

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

Single-Atom Doping at the Molecule–Metal Interface: How Rh Affects Surface Explosion Kinetics

Controlling chemical reactions at interfaces is central to many fields, including atmospheric, environmental, and biological chemistry as well as catalysis and corrosion. We performed a fundamental study of how catalytically active Rh atoms placed at and above a molecule-metal interface influence the surface reaction kinetics of tartaric acid (TA) decomposition. Specifically, TA decomposition on Cu(100) exhibits autocatalytic decomposition kinetics involving a slow initiation step followed by fast decomposition of the molecular layer. These so-called “surface explosions” are extremely sensitive to initial conditions, making it fundamentally interesting to study how the placement of reactive atoms affects the reaction rate. Temperature-programmed reaction (TPR) experiments reveal that Rh atoms embedded in the Cu(100) surface beneath TA (TA/Rh/Cu) or atop the TA layer (Rh/TA/Cu) reduce the decomposition temperature and modify the reaction rate constants. Isothermal TPR analysis reveals that Rh enhances both the initiation rate constant, k i , and the explosion rate constant, k e , facilitating explosive decomposition at lower temperatures. This result provides evidence that both the initiation and explosion steps occur at the metal-molecule interface and are accelerated by the presence of Rh at this interface. This study illustrates how small amounts of reactive elements may be used to control nonlinear kinetic processes at interfaces.

desorption↗
DOE OSTI2026

Catalytic resonance theory for the kinetics of photon-promoted catalysis

The illumination of catalytic surfaces with a continuous or pulsed stream of photons dynamically modulates surface chemistry for faster rates, non-equilibrium conversion, or product selectivity control. To establish fundamental principles of dynamic photon-modulated catalysis, the photocatalytic conversion of a generic surface reaction was simulated using the kinetic Monte Carlo method to understand the kinetic implications of an independent stream of photons that promotes surface product desorption. The time-averaged photocatalytic rate at differential conditions for varying photon flux and temperatures indicated three kinetic regimes described by product thermal desorption control, surface reaction control, and an intermediate kinetic regime with a zero slope Arrhenius plot, consistent with a degree of rate control dominated by the photon arrival frequency (i.e., per-site photon flux). Here, the maximum photocatalytic rate occurred orders of magnitude above the Sabatier limit at the resonance frequency, identified as the photon arrival frequency matching the surface reaction rate constant.

Canavan, Jesse R. [University of Minnesota, Minnea↗
DOE OSTI2021

Analytical modeling for redox flow battery design

Deeper market penetration of redox flow batteries requires optimization of the cell performance. Though important for optimization, detailed analytical solutions have not been developed for electrolyte flow, mass and charge transport, and reaction kinetics within redox flow batteries. To this end, here we present analytical solutions to active species concentration and over-potential based on advection-diffusion transport for ions and Bulter-Volmer model for interface reaction kinetics. The solutions were validated with results from a finite element model. These solutions were then applied to investigate the relationship between over-potential and state of charge, current density, reaction rate constant, flow velocity, diffusivity, total active species concentration, and electrode structure. Explicit formulas were identified for minimum activation over-potential and limiting current density as well as their dependence on electrolyte properties, operation conditions, and electrode structure. With our new mathematical formulas, this work provides a theoretical framework for flow battery design.

25 ENERGY STORAGE↗
DOE OSTI2025

Catalytic Resonance Theory: Experimental and Kinetic Interpretation of Programmable Catalysis

The complexity of programmable catalysts that change in the physical or electronic state on the time scale of a catalytic turnover results in uncertainty in the kinetic outcome of common measurements with dynamic catalyst experiments. To provide an interpretation of utility to experimentalists, a general programmable catalytic reaction exhibiting high turnover efficiency was simulated to understand the kinetic behavior of catalysts with varying frequency, binding energy, oscillation amplitude, and minimum oscillation binding energy at different temperatures and reactant gas pressures. The time-averaged catalytic rate for varying temperature in the form of an Arrhenius plot exhibited three distinct kinetic regimes, with the intermediate temperatures or applied frequencies exhibiting zero slope, indicative of barrier-free kinetic control. Transitions in experimentally measurable kinetic regimes with both temperature and applied catalyst oscillation frequency were associated with transitions in the sensitivity of reaction rate constants and degrees of rate control. Furthermore, these distinct kinetic macroscopic features specific to programmable catalysts were identified for observation in experimental dynamic kinetic measurements.

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

Applying a Phase-Separation Parameterization in Modeling Secondary Organic Aerosol Formation from Acid-Driven Reactive Uptake of Isoprene Epoxydiols under Humid Conditions

Secondary organic aerosol (SOA) from acid-driven reactive uptake of isoprene epoxydiols (IEPOX) contributes up to 40% of organic aerosol (OA) mass in fine particulate matter. Our previous work showed substantial conversions of particulate inorganic sulfates to surface-active organosulfates (OSs) by IEPOX decreases aerosol acidity and creates a viscous organic-rich shell that poses as a diffusion barrier, inhibiting additional reactive uptake of IEPOX. In this study, to account for this “self-limiting” effect, a phase-separation box model was developed to evaluate parameterizations of IEPOX reactive uptake against time-resolved chamber measurements of IEPOX-SOA tracers, including 2-methyltetrols (2-MT) and methyltetrol sulfates (MTS), at ~ 50% relative humidity. The phase-separation model was most sensitive to mass accommodation coefficient, IEPOX diffusivity in the organic shell, and ratio of the third-order reaction rate constants forming 2-MT and MTS (k MT /k MTS ). In particular, k MT /k MTS had to be lower than 0.1 to bring model predictions of 2-MT and MTS in closer agreement with chamber measurements, while prior studies reported values larger than 0.71. The model-derived rate constants favor more particulate MTS formation due to 2-MT likely off-gassing at ambient-relevant OA loadings. Incorporating this parametrization into chemical transport models is expected to predict lower IEPOX-SOA mass and volatility due to the predominance of OSs.

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

Development of an Energy-Efficient and High-Productivity Ammonia Recovery and Removal Process Using Resin-Wafer Electrodeionization

As efforts to develop various energy resources, ammonia energy is emerging as a promising carbon-free fuel alternative. Recovering high-concentration ammonia and ammonium from wastewater potentially offers significant environmental and economic benefits. However, research on recovery technologies for industries such as semiconductors remains limited. This study develops an efficient, energy-saving ammonia recovery technology using Resin Wafer Electrodeionization (RW-EDI), specifically for the semiconductor industry. RW-EDI shows promise for recovering ammonia from high-concentration wastewater. By optimizing parameters, such as voltage and initial concentration, a balance between productivity and energy consumption is achieved. Results indicate that ammonium and fluoride ion transport kinetics are similar, with minimal competition and selectivity values between 0.95 and 1.08. As the initial ammonia/ammonium concentration increases from 500 to 8000 ppm, reaction rate constants and the overall mass transfer coefficient decrease. Increasing the voltage can enhance mass transfer and eliminate barriers. Ion transport primarily occurs in the RW solid phase, accounting for 93.89% of the total current. Additionally, RW-EDI shows superior specific energy consumption within a concentration range of 500 ppm to 8000 ppm, outperforming technologies like electrodialysis by reducing specific energy consumption from 8−15 kWh/kg-NH 4 + to 1.2−2.5 kWh/kg-NH 4 + . This study highlights RW-EDI’s potential for ammonia recovery, providing valuable insights for future applications in wastewater treatment.

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

On the mechanism of soot nucleation. II. E-bridge formation at the PAH bay

A recently proposed mechanism of soot nucleation based upon the formation of a rotationally-activated dimer in the collision of an aromatic molecule and a radical leading to a stable, doubly-bonded E-bridge between them, rooted in the existence of a five-membered ring on the molecule edge, has been further investigated with a focus on the 5-6 E-bridge forming in the reaction of a PAH cyclopenta group with a bay site of another PAH. As a prototype reaction of this kind, we examined the reaction between 4-phenanthrenyl and acenaphthylene and, to project these results to larger aromatic structures, we also explored the size effect of the E-bridge forming reactions by computing the 1-naphthyl + acenaphthylene system and comparing these results with the prior results for pyrenyl + acepyrene. The two systems have been studied through high-level G3(MP2,CC)//B3LYP/6-311G(d,p) ab initio calculations of their potential energy surfaces combined with the RRKM/Master Equation calculations of reaction rate constants. With PAH monomers of similar sizes involved, the formation of E-bridge structures at the bay radical sites appeared to be faster and lead to increased nucleation rates as compared to the zigzag-forming ones. Furthermore, a model combining both the bay and zigzag rotationally-induced formation of E-bridges successfully reaches the level of nucleation fluxes comparable to those of the irreversible pyrene dimerization, thus affirming the rotationally-activated dimerization as a feasible mechanism for soot particle nucleation.

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

Deciphering the morphology of transition metal carbonate cathode precursors

The performance and life of Li-ion battery cathode materials are determined by both the composition (crystal structure and transition metal ratio) and the morphology (particle size, size distribution, and surface area). Careful control of these two aspects is the key to long lasting, high-energy batteries that can undergo fast charge. Developing such cathodes requires manipulation of the synthesis conditions, namely the coprecipitation process to develop the precursor and a calcination step to lithiate and convert it to a transition metal oxide. In this paper, we utilize a combination of controlled synthesis, microscopic and spectroscopic characterization, and multi-scale mathematical modeling to shed light on the synthesis of cathode precursors. The complex interplay between the various chemical reactions in the co-precipitation process is studied to provide experimentalists with guidance on achieving composition control during synthesis. Further, the formation of a variety of morphologies of the primary particles and the driving force for agglomeration is mathematically described, for the first time, based on an energy minimization approach. Results suggest that the presence of Ni and/or Co significantly lowers the reaction rate constant compared to Mn, resulting in agglomerated growth in the former and single crystal growth in the latter. Modeling studies are used to provide a phase map describing the synthesis conditions needed to control the secondary particle size and corresponding size distribution. Finally, this paper represents an important step in developing a computationally guided approach to the synthesis of battery cathode materials.

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

2,4–dinitroanisole transformation during advanced oxidation with ultraviolet light emitting diodes and hydrogen peroxide

This study investigated the kinetics and reaction mechanisms associated with advanced oxidation of 2,4-dinitroanisole (DNAN) with light emitting diodes (LEDs) in the UV/H2O2 process. The pseudo first-order reaction rate constants for DNAN oxidation depended on the H 2 O 2 /DNAN molar ratio and were between 0.017 and 0.026 min –1 . Hydroxyl radical scavenging was strongly implicated by the experimental results. High mass spectrophotometric chromatograms showed the presence of several by-products, including three with ecotoxicity that is greater than that of DNAN. Transformation pathways were analysed using Density Functional Theory (DFT), which helped to predict the presence of experimentally detected DNAN by-products by considering Mulliken charges. Furthermore, this study is the first to elucidate DNAN transformation mechanisms using both experimental studies and computational chemistry.

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

Models to Incorporate Reaction Mechanisms into DG-OSPREY: Fixed-Bed Simulations for Organic Iodides Capture Using Ag 0 Z

As one of the most potent radioisotopes released during spent nuclear fuel reprocessing, 129 I is strictly regulated and must be removed before discharge. Organic iodides (primarily alkyl iodides with different chain lengths, i.e., CH 3 I, C 4 H 9 I, and C 12 H 25 I) comprise ~2% of the total iodine in the reprocessing off-gases and are primarily present in vessel off-gas (VOG). Reduced silver mordenite (Ag 0 Z) is predominantly considered for the removal of radioiodine; however, its capture performance and underlying interaction processes with long-chain organic iodides are not fully understood. Two major tasks were accomplished in this study. First, to improve upon the previous experimental studies where Ag 0 Z was used to capture CH 3 I, C 4 H 9 I, and C 12 H 25 I at different concentrations, we comprehensively investigated the corresponding capture mechanisms by characterizing fully loaded Ag 0 Z samples. Second, computational codes were implemented to perform fixed-bed simulations that account for transport and reaction mechanisms. Scanning electron microscopy with energy dispersive X-ray analysis (SEM-EDX), powder X-ray diffraction (PXRD), UV-visible diffuse reflectance spectroscopy (UV-vis DRS), and thermogravimetric analysis (TGA) were conducted on Ag 0 Z samples that are saturated with I 2 , CH 3 I, C 4 H 9 I, and C 12 H 25 I), respectively. Results indicate that AgI is the predominant adsorption product regardless of the adsorbed iodine species, yet alkyl iodides with different carbon chain lengths may have different compositions of α- and γ-AgI. Synchrotron pair distribution function (PDF) measurements and TGA coupled with a Fourier transformed infrared detector (TGA FTIR) have been performed, and experimental data are currently analyzed. Results are expected to provide further insights into the adsorption mechanisms. The fixed-bed capture performance for CH 3 I was successfully simulated by solving mathematical equations that describe the underlying transport processes and adsorption reactions. The computational framework, Catalytic After Treatment System (CATS), that was originated in our research group, was used to solve the governing equations. Kinetic parameters, including the pore diffusivity and reaction rate constant were obtained by optimization techniques using data from thin-bed experiments performed at Oak Ridge National Laboratory. The performance of a deep bed predicted using the optimized parameters showed promising agreement with the experimental data.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2023

Kinetics and Mechanism of Tetravalent Cerium Extraction by Di-2-Ethylhexyl Isobutyramide (DEHiBA)

DEHiBA is a fully-incinerable, CHON extractant with a high selectivity for hexavalent actinides. However, as a recently-developed ligand, much remains unknown about its physical and chemical properties, including its extraction and phase disengagement rates. Here, we present data on the phase transfer kinetics of Ce(IV) extraction and stripping by both DEHiBA and TBP collected using microfluidic devices. Ce(IV) was selected as a surrogate for Pu(IV) and Np(IV) in the initial tests of the microfluidics system prior to testing with the actinides. Pseudo first-order reaction rate constants are determined for both systems, and the corresponding Damköhler numbers (Da) are calculated to establish whether the extraction and stripping processes are kinetically or diffusion controlled. The stoichiometry of the extracted Ce(IV)-DEHiBA complex is determined, and the phase disengagement properties under extraction conditions are assessed using the dimensionless dispersion number. A possible explanation for the favorable extraction of Ce(IV) by DEHiBA, which is known to poorly extract many tetravalent cations, is presented.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗