Search NASA⌕ Search

SEARCH · Search NASA

Results for “growth kinetics”

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

Electroconvection near an ion-selective surface with Butler–Volmer kinetics

We study the effects of interfacial kinetics on the electro-hydrodynamics of ion transport near an ion-selective surface using a combination of linear stability analysis and numerical simulation. The finite kinetics of the electrolyte–electrode interface affects the ion transfer and electroconvection in many ways. On a surface of fixed topography, such as a metal surface of slow and stable ion deposition or covered by a polymer membrane, the finite kinetics reduces the current in one-dimensional ion diffusion/migration, increases the critical voltage for the onset of the electroconvective instability, changes the dynamics of the electroconvection and the overlimiting current, and enhances the lateral ion diffusion within the interfacial layer. The first three effects are indirectly caused by the reaction kinetics and can be characterized by an effective voltage difference across the liquid electrolyte. In comparison, the last effect is controlled by a direct interplay between kinetics and nonlinear electroconvection. Scaling laws for ion transport and features of electroconvection are proposed. We also analyse the linear stability of a surface which evolves under ion deposition and find that the finite kinetics decreases the growth rate of both electroconvective and morphological instabilities and therefore modifies the wavenumber of the most unstable mode.

Mechanics↗
DOE OSTI2020

Semiconducting High-Entropy Chalcogenide Alloys with Ambi-ionic Entropy Stabilization and Ambipolar Doping

Entropy stabilization is a novel materials-design paradigm to realize new compounds with widely tunable properties. However, almost all entropy-stabilized materials so far are either conducting metals or insulating ceramics, with a clear dearth in the semiconducting regime. Here, a new class of the multi-cationic and -anionic entropy-stabilized chalcogenide alloys based on the (Ge,Sn,Pb)(S,Se,Te) formula are synthesized and characterized experimentally. The con-figurational entropy from the disorder of both the anion and the cation sublattices reaches a record value of ~2.2 R mol -1 for the equimolar composition and stabilizes the single-phase solid solution. Theoretical calculations and experiments both show that the synthesized alloys are thermodynamically stable at the growth temperature and kinetically metastable at room temperature, segregating by spinodal decomposition at moderate temperatures. In conclusion, doping and electronic transport measurements verify that the synthesized materials are ambipolarly dopable semiconductors, which paves the way for the wider adoption of entropy-stabilized chalcogenide alloys in functional applications.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Size-Dependent Adsorption and Adhesion Energetics of Ag Nanoparticles on Graphene Films on Ni(111) by Calorimetry

Interest in the use of carbon supports for late transition metal nanoparticle catalysts has expanded rapidly due to the increasing importance of electrocatalysts for clean energy and environmental technologies and the use and storage of renewable electricity. Compared to oxide supports, almost nothing is known about the effect of metal nanoparticle size on the energies of the metal atoms within carbon-supported nanoparticles, yet these energies are crucial for understanding their surface reactivity and sintering kinetics. Here, the growth morphology and adsorption energetics of vapor-deposited Ag onto clean graphene/Ni(111) surfaces have been studied using a combination of single-crystal adsorption calorimetry (SCAC) and He + low-energy ion scattering (LEIS). The differential heat of Ag adsorption is 207 kJ/mol for making ~30 atom Ag particles on graphene terraces at 100 K and 16 kJ/mol higher for making ~9 atom Ag clusters at defect sites at the same temperature. The heat of adsorption increases rapidly with Ag coverage as 3D Ag nanoparticles nucleate and grow in size, asymptotically reaching within 5 kJ/mol of the bulk Ag sublimation enthalpy (285 kJ/mol) by 2 ML. The heats of adsorption and Ag nanoparticle densities from LEIS (~10 16 /m 2 ) were combined to provide the Ag/graphene adhesion energy (E adh = 1.8 J/m 2 in the large-particle limit) and the Ag chemical potential (μ) versus effective particle diameter (D). The Ag chemical potential was well-fitted by μ(D) = (3γ v /M – E adh )(1 + (1.5 nm)/D)(2V m /D), where γ v /M is the surface energy of bulk Ag and V m is its molar volume. The same equation is known to fit similar data for late transition metals on clean surfaces of metal oxide single crystals. The adhesion energy of Ag measured here on graphene falls within the wide range measured for Ag on those oxide surfaces and is almost as large as on the oxide that binds Ag particles most strongly, namely CeO 2 (111), which is well-known to be very effective at resisting catalyst deactivation by metal sintering. Furthermore, these results imply that carbon supports will be effective at resisting sintering and that Ag particles smaller than 6 nm on graphene will bind small adsorbed reaction intermediates more weakly than supports with weaker adhesion to Ag, like MgO(100).

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

Multiscale modeling of metal-hydride interphases—quantification of decoupled chemo-mechanical energies

Abstract The quantification of interphase properties between metals and their corresponding hydrides is crucial for modeling the thermodynamics and kinetics of the hydrogenation processes in solid-state hydrogen storage materials. In particular, interphase boundary energies assume a pivotal role in determining the kinetics of nucleation, growth, and coarsening of hydrides, alongside accompanying morphological evolution during hydrogenation. The total interphase energy arises from both chemical bonding and mechanical strains in these solid-state systems. Since these contributions are usually coupled, it is challenging to distinguish via conventional computational approaches. Here, a comprehensive atomistic modeling methodology is developed to decouple chemical and mechanical energy contributions using first-principles calculations, of which feasibility is demonstrated by quantifying chemical and elastic strain energies of key interfaces within the FeTi metal-hydride system. Derived materials parameters are then employed for mesoscopic micromechanical analysis, predicting crystallographic orientations in line with experimental observations. The multiscale approach outlined verifies the importance of the chemo-mechanical interplay in the morphological evolution of growing hydride phases, and can be generalized to investigate other systems. In addition, it can streamline the design of atomistic models for the quantitative evaluation of interphase properties between dissimilar phases and allow for efficient predictions of their preferred phase boundary orientations.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Dynamics of homogeneous cavitation with pressure feedback

Theoretical studies of homogeneous cavitation have largely been based on the classical nucleation theory. However, existing cavitation models cannot adequately describe its dynamics at nanosecond timescale, which has been called for in other fields. We develop a model coupling nucleation kinetics with cavity growth and pressure feedback as saturation mechanisms. Our numerical studies reveal the exponential dependence of cavitation characteristics such as saturation cavity density and most probable cavity radius on model parameters: Tolman length and initial pressure. This work also sheds light on the possibility of accurately determining Tolman length, whose value has a large spread in the literature.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

Saturated random packing built of arbitrary polygons under random sequential adsorption protocol

Random packings and their properties are a popular and active field of research. Numerical algorithms that can efficiently generate them are useful tools in their study. This paper focuses on random packings produced according to the random sequential adsorption (RSA) protocol. Developing the idea presented by G. Zhang, where saturated random packings built of regular polygons were studied, we create an algorithm that generates strictly saturated packings built of any polygons. Then, the algorithm was used to determine the packing fractions for arbitrary triangles. In this work, the highest mean packing density, 0.552 814 ± 0.000 063, was observed for triangles of side lengths 0.63 : 1 : 1. Additionally, microstructural properties of such packings, kinetics of their growth, as well as distributions of saturated packing fractions and the number of RSA iterations needed to reach saturation were analyzed.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Mechanism-Based Energy Regularization in Computational Modeling of Quasibrittle Fracture

Abstract Quasibrittle materials are featured by a strain-softening constitutive behavior under many loading scenarios, which could eventually lead to localization instability. It has long been known that strain localization would result in spurious mesh sensitivity in finite element (FE) simulations. Previous studies have shown that, for the case of fully localized damage, the mesh sensitivity can be mitigated through energy regularization of the material constitutive law. However, depending on the loading configuration and structural geometry, quasibrittle structures could exhibit a complex damage process, which involves both localized and diffused damage patterns at different stages of loading. This study presents a generalized energy regularization method that considers the spatial and temporal evolution of damage pattern. The method introduces a localization parameter, which describes the local damage pattern. The localization parameter governs the energy regularization of the constitutive model, which captures the transition from diffused to localized damage during the failure process. The method is cast into an isotropic damage model, and is further extended to rate-dependent behavior. The energy regularization scheme is directly incorporated into the kinetics of damage growth. The model is applied to simulate static and dynamic failures of ceramic specimens. It is shown that the present model is able to effectively mitigate the spurious mesh sensitivity in FE simulations of both types of failure. The present analysis demonstrates the essential role of mechanism-based energy regularization of constitutive relation in FE simulations of quasibrittle fracture.

Mechanics↗
DOE OSTI2025

Nickel hydroxide–nickel carbonate competitive growth on carbonate surfaces

The thermodynamic and kinetic factors controlling the competitive heterogeneous nucleation and growth of ubiquitous metal carbonate and hydroxide phases are poorly understood. In this work, calcite (CaCO 3 ) and magnesite (MgCO 3 ) powders were reacted with NiCl 2 (0–600 μM) for 7 days at 22 °C and 5 °C. The reacted powders were analyzed with X-ray photoelectron spectroscopy (XPS), scanning electron microscopy, and energy-dispersive X-ray spectroscopy to characterize the Ni surface precipitates formed. Evidence from these techniques pointed to the formation of mixed Ni carbonate-Ni hydroxide amorphous surface precipitates. On calcite, XPS detected primarily Ni(OH) 2 despite the initial solutions being more supersaturated with respect to gaspéite (NiCO 3 ) than to theophrastite (Ni(OH) 2 ) by a factor of 17–18. In contrast, NiCO 3 was the dominant component detected by XPS on magnesite in the same conditions. Decreasing the temperature had the effect of increasing the proportion of NiCO 3 to the detriment of Ni(OH) 2 . The experimental observations were consistent with low lattice/cation size mismatch favoring NiCO 3 nucleation and temperature most influencing Ni(OH) 2 nucleation. Comparison to previous work on Co-reacted powders indicated the differences in lattice/cation size mismatch and/or water exchange rate impacted the composition of the surface precipitates more than the relative thermodynamic stabilities of the competing minerals. This work explored the heterogeneous growth regime of Ni carbonate and hydroxide phases on carbonate surfaces and shed light on the factors that control the competition between surface precipitates when mineral surfaces are in contact with aqueous solutions supersaturated with respect to multiple mineral phases. In conclusion, these results contribute to geochemists’ efforts toward interpreting data from geochemical systems with elevated Ni concentrations, improving Ni environmental remediation and recovery strategies, and predicting the fate and transport of Ni in geochemical systems.

Heterogeneous nucleation↗
DOE OSTI2026

Time-Evolved Hetero-Alkali Interphases Enable Long-Life Sulfide-Based Anode-Free Solid-State Batteries

Sulfide-based anode-free solid-state batteries (AFSSBs) offer compelling advantages in terms of energy density and safety, yet their practical implementation is severely hindered by undesirable interfacial reactions between sulfide solid electrolytes (SEs) and freshly plated lithium (Li), as well as non-uniform Li plating/stripping behavior. Herein, an effective interfacial stabilization strategy by incorporating sodium bis(fluorosulfonyl)imide (NaFSI) additive into the Li5.4PS4.4Cl1.6 (LPSC) is investigated. Unlike conventional Li-based additives that form static passivation layers, NaFSI introduces a transient hetero-alkali chemistry that kinetically governs interphase evolution during fresh Li plating. NaFSI induces a timesequenced interphase evolution: an initial NaF/LiF-rich layer that suppresses early sulfide reduction, followed by a LiF/Li3N-rich layer that optimizes Li⁺ transport during repeated anode-free cycling. This evolved robust and fast ion conducting layer mitigates interfacial impedance growth, enhances Li + transport kinetics, and suppresses localized Li growth and filamentary shorting. As a result, the anode-free full cell with NaFSI modified LPSC as the interlayer exhibits an excellent cycling stability over 500 cycles at 0.2 C with a capacity retention of 77.6%, whereas the cell with bare LPSC suffers from rapid capacity decay after 100 cycles, retaining only 32.1% of its initial capacity. This work establishes dynamic heteroalkali additive chemistry as a general strategy to kinetically program solid-solid interphases, guiding the interface design in anode-free solid-state batteries.

25 ENERGY STORAGE↗
DOE OSTI2020

A quantitative phase-field model for gas bubble evolution in UO2

Due to the large formation energy of vacancies and noble gas atoms in the form of interstitials or substitutional atoms in nuclear fuel (UO2), the thermodynamic equilibrium concentrations of these species are very low in the nuclear fuel matrix even at very high temperatures, which imposes difficulties upon the quantitative study of bubble evolution via the phase-field method. In this study, a quantitative phase-field model is proposed to deal with this problem. The system’s free energy density is derived according to the principles of thermodynamics, with consideration of the elastic effect and with the use of real material parameters from experiments. The model is useful for the study of the kinetics of gas bubble growth with very dilute concentrations of vacancy and gas atoms in the matrix. This model is applied to study single bubble growth and multiple bubble growth under various concentrations of vacancy and gas atoms and at various temperatures. The elastic effect and the effects of the generation rate of vacancies and gas atoms on bubble growth are analyzed.

low defect concentration, nuclear fuel, gas bubble↗
DOE OSTI2021

Stability Analysis of Parameter Varying Genetic Toggle Switches Using Koopman Operators

The genetic toggle switch is a well known model in synthetic biology that represents the dynamic interactions between two genes that repress each other. The mathematical models for the genetic toggle switch that currently exist have been useful in describing circuit dynamics in rapidly dividing cells, assuming fixed or time-invariant kinetic rates. There is a growing interest in being able to model and extend synthetic biological function for growth conditions such as stationary phase or during nutrient starvation. As cells transition from one growth phase to another, kinetic rates become time-varying parameters. In this paper, we propose a novel class of parameter varying nonlinear models that can be used to describe the dynamics of genetic circuits, including the toggle switch, as they transition from different phases of growth. We show that there exists unique solutions for this class of systems, as well as for a class of systems that incorporates the microbial phenomena of quorum sensing. Further, we show that the domain of these systems, which is the positive orthant, is positively invariant. We also showcase a theoretical control strategy for these systems that would grant asymptotic monostability of a desired fixed point. We then take the general form of these systems and analyze their stability properties through the framework of time-varying Koopman operator theory. A necessary condition for asymptotic stability is also provided as well as a sufficient condition for instability. A Koopman control strategy for the system is also proposed, as well as an analogous discrete time-varying Koopman framework for applications with regularly sampled measurements.

Harrison, Jamiree (ORCID:0000000265351579)↗
DOE OSTI2023

High-Throughput Automated Exploration of Phase Growth Behaviors in Quasi-2D Formamidinium Metal Halide Perovskites

Quasi-2D metal halide perovskites (MHPs) are an emerging material platform for sustainable functional optoelectronics, but the uncontrollable, broad phase distribution remains a critical challenge for applications. Nevertheless, the basic principles for controlling phases in quasi-2D MHPs remain poorly understood, due to the rapid crystallization kinetics during the conventional thin-film fabrication process. In this work, a high-throughput automated synthesis-characterization-analysis workflow is implemented to accelerate material exploration in formamidinium (FA)-based quasi-2D MHP compositional space, revealing the early-stage phase growth behaviors fundamentally determining the phase distributions. Upon comprehensive exploration with varying synthesis conditions including 2D:3D composition ratios, antisolvent injection rates, and temperatures in an automated synthesis-characterization platform, it is observed that the prominent n = 2 2D phase restricts the growth kinetics of 3D-like phases—α-FAPbI 3 MHPs with spacer-coordinated surface—across the MHP compositions. Thermal annealing is a critical step for proper phase growth, although it can lead to the emergence of unwanted local PbI 2 crystallites. Additionally, fundamental insights into the precursor chemistry associated with spacer-solvent interaction determining the quasi-2D MHP morphologies and microstructures are demonstrated. The high-throughput study provides comprehensive insights into the fundamental principles in quasi-2D MHP phase control, enabling new control of the functionalities in complex materials systems for sustainable device applications.

2D perovskites↗
DOE OSTI2024

Phase field modeling of hydrogen release in nitrile-butadiene rubber composites after high-pressure hydrogen exposure

Experimental data showed that a simple diffusion model doesn’t suffice to explain the hydrogen (H 2 ) release kinetics in nitrile butadiene rubber (NBR) composites (one with no filler or plasticizer, one with plasticizer only, one with fillers and plasticizer, and one with fillers only) after exposed to hydrogen gas at a pressure of 27.6 MPa. In this work, we developed a phase field model that considers the effect of H 2 gas bubble evolution on H 2 release kinetics. With the model, we simulated the effect of thermodynamic and kinetics properties on the nucleation, growth, and shrinkage of gas bubbles and H 2 release kinetics during and after decompression. Here, the results demonstrated that 1) the model built upon certain thermodynamic and kinetics properties of a given material can well describe its H 2 release kinetics measured in experiments, and 2) the model can predict the gas bubble evolution, which is associated with the material property degradation and failures, with accurate thermodynamic and kinetics properties or H 2 release data from experiments.

08 HYDROGEN↗
DOE OSTI2023

An Integrated Simulation of Multiple-Pass U-10Mo Alloy Hot Rolling and Static Recrystallization

To achieve a desired microstructure and minimize the thickness variation in rolled foils, researchers must understand the effects of foil fabrication process variables on microstructure evolution. We developed an integrated simulation of deformation and recrystallization that employs the finite element method (FEM) and the kinetic Monte Carlo (KMC) Potts model, respectively, to investigate microstructure evolution during multiple-pass hot rolling and heat treatment in polycrystalline U-10Mo fuel. Scanning electron microscopy and electron backscatter diffraction images of microstructures were directly used as input in FEM calculation of deformation, and the calculated strains were used to determine the driving force of nucleation and growth of recrystallized grains in the Potts model. Grain structures predicted by the Potts model were used to update the grain structure and material properties for FEM. Simulation alternated between FEM and the Potts model to simulate grain structure evolution during multiple rolling and heat treatments. The initial model parameters were determined by benchmarking the recrystallization kinetics against experimental data. Then, the model was applied to predict the grain structure evolution. Results showed that our model can capture the coupling between deformation and recrystallization and can quantitatively reproduce the observed U-10Mo recrystallization and grain growth kinetics. The simulation results demonstrated that the developed model can predict U-10Mo grain structures as a function of initial microstructure and foil fabrication parameters.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Phase-field model of freeze casting

Directional solidification of water-based solutions has emerged as a versatile technique for templating hierarchical porous materials. However, the underlying mechanisms of pattern formation remain incompletely understood. In this work, we present a detailed derivation and analysis of a quantitative phase-field model for simulating this nonequilibrium process. The phase-field model extends the thin-interface formulation of dilute binary alloy solidification with antitrapping to incorporate the highly anisotropic energetic and kinetic properties of the partially faceted ice-water interface. This interface is faceted in the basal plane normal to the ⟨0001⟩ directions and atomically rough in other directions within the basal plane. On the basal plane, the model reproduces a linear or nonlinear relationship between the interface growth rate and the kinetic undercooling that can be linked to experimental measurements. In both cases, spontaneous parity breaking of the solidification front is observed when the preferred growth direction is aligned with the temperature gradient. This phenomenon leads to the formation of partially faceted ice lamellae that drift laterally in one of the ⟨0001⟩ directions. Here, we demonstrate that the drifting velocity of the ice lamellae is controlled by the kinetics on the basal plane and converges as the thickness of the diffuse solid-liquid interface decreases. Furthermore, we examine the effect of the form of the kinetic anisotropy, which is chosen here such that the inverse of the kinetic coefficient varies linearly from a finite value in the ⟨0001⟩ directions to zero in all other directions within the basal plane, consistent with the assumption that the interface grows in local thermodynamic equilibrium in this plane. Our results indicate that the drifting velocity of ice lamellae is not affected by the slope of this linear relation, and the radius and undercooling at the tip of an ice lamella converge at relatively small slope values. Consequently, the phase-field simulations remain quantitative with computationally tractable choices of both the interface thickness and the slope assumed in the form of the kinetic anisotropy.

Materials science↗
DOE OSTI2021

Kinetically Controlled Synthesis of Rhodium Nanocrystals with Different Shapes and a Comparison Study of Their Thermal and Catalytic Properties

Here, we report the synthesis of Rh nanocrystals with different shapes by controlling the kinetics involved in the growth of preformed Rh cubic seeds. Specifically, Rh nanocrystals with cubic, cuboctahedral, and octahedral shapes can all be obtained from the same cubic seeds under suitable reduction kinetics for the precursor. The success of such a synthesis also relies on the use of a halide-free precursor to avoid oxidative etching, as well as the involvement of a sufficiently high temperature to remove Br – ions from the seeds while ensuring adequate surface diffusion. The availability of Rh nanocrystals with cubic and octahedral shapes allows for an evaluation of the facet dependences of their thermal and catalytic properties. The data from in situ electron microscopy studies indicate that the cubic and octahedral Rh nanocrystals can keep their original shapes up to 700 and 500 °C, respectively. When tested as catalysts for hydrazine decomposition, the octahedral nanocrystals exhibit almost 4-fold enhancement in terms of H 2 selectivity relative to the cubic counterpart. As for ethanol oxidation, the order is reversed, with the cubic nanocrystals being about three times more active than the octahedral sample.

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

High-accuracy method for modeling nucleation and growth of particles

State-of-the-art numerical models describing the kinetics of aerosol particle nucleation and growth from a cooling vapor primarily use a nodal method, in which particles that are smaller than the critical size are omitted from consideration because they are thermodynamically unfavorable. This omission is based on the assumption that most newly formed particles are above the critical size, so that subcritical-size particles are not important to take into account. Due to the nature of the nodal method, it suffers from numerical diffusion, which can cause an artificial broadening of the cluster size distribution leading to a significant overestimation of the number of large-size particles. To address these issues, we propose a more accurate numerical method that explicitly models particles of all sizes, and uses a special numerical scheme that substantially reduces the numerical diffusion and provides high solution accuracy and numerical stability. We extensively compare this novel method to the commonly used nodal solver of the general dynamic equation (GDE) for particle growth and demonstrate that it offers GDE solutions with higher accuracy with low numerical diffusion. Incorporating small subcritical clusters into the solution is crucial for: 1) more precise determination of the entire particle size distribution function and 2) wider applicability of the model to experimental studies with non-monotonic temperature variations leading to particle evaporation. The computational code implementing this numerical method in Python is available upon request.

42 ENGINEERING↗
DOE OSTI2020

Solvent manipulation of the pre-reduction metal–ligand complex and particle-ligand binding for controlled synthesis of Pd nanoparticles

Understanding how to control the nucleation and growth rates is crucial for designing nanoparticles with specific sizes and shapes. Herein, we show that the nucleation and growth rates are correlated with the thermodynamics of metal–ligand/solvent binding for the pre-reduction complex and the surface of the nanoparticle, respectively. To obtain these correlations, we measured the nucleation and growth rates by in situ small angle X-ray scattering during the synthesis of colloidal Pd nanoparticles in the presence of trioctylphosphine in solvents of varying coordinating ability. The results show that the nucleation rate decreased, while the growth rate increased in the following order, toluene, piperidine, 3,4-lutidine and pyridine, leading to a large increase in the final nanoparticle size (from 1.4 nm in toluene to 5.0 nm in pyridine). Using density functional theory (DFT), complemented by 31 P nuclear magnetic resonance and X-ray absorption spectroscopy, we calculated the reduction Gibbs free energies of the solvent-dependent dominant pre-reduction complex and the solvent-nanoparticle binding energy. The results indicate that lower nucleation rates originate from solvent coordination which stabilizes the pre-reduction complex and increases its reduction free energy. At the same time, DFT calculations suggest that the solvent coordination affects the effective capping of the surface where stronger binding solvents slow the nanoparticle growth by lowering the number of active sites (not already bound by trioctylphosphine). The findings represent a promising advancement towards understanding the microscopic connection between the metal–ligand thermodynamic interactions and the kinetics of nucleation and growth to control the size of colloidal metal nanoparticles.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗