Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “time dependent behavior”

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 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Deviations from the Isobaric Multiplet Mass Equation due to Threshold States

Recent studies have completed the A=16 isospin quintets for states with J π =0 + and 2 + . The dependence of their masses as a function of isospin projection shows evidence for deviations from quadratic behavior indicating isospin violation beyond the expectation from two-body forces. The deviation is most pronounced for the 2 + states. Predictions from the shell model embedded in the continuum (SMEC) allow us to explain that this isospin violation is associated with a modification of the nuclear structure due to the open-quantum-system nature of the proton-rich members of the quintet. In particular, the 0 + and 2 + states in 16 Ne and the 2 + state in 16 F are threshold resonances located just above a proton-decay threshold where s-wave coupling to the continuum is expected. The measured deviations of these threshold states from the quadratic behavior of the remaining members of the multiplets makes it possible to obtain information on the magnitude and the energy dependence of the continuum-coupling energy correction. Finally, continuum coupling is also indicated for the ground state of 8 C, but this time through p-wave coupling

coulomb energies & analogue states↗
DOE OSTI2020

Polymer Retention Evaluation in a Heavy Oil Sand for a Polymer Flooding Application on Alaska’s North Slope

For a polymer flooding field trial in a heavy oil reservoir on Alaska’s North Slope, polymer retention is a key parameter. Because of the economic impact of retention, this parameter was extensively studied using field core material and conditions. Several important findings were noted. Polymer retention based on effluent viscosity measurements can be overestimated unless the correct (non-linear) relation between polymer concentration and viscosity is used. Polymer degradation (either mechanical or oxidative) can also lead viscosity-based measurements to overestimate retention. Inaccessible pore volume (IAPV) can be overestimated if insufficient brine is flushed through the sand between polymer banks. Around 100 pore volumes (PV) of brine may be needed to displace mobile polymer to approach a true residual resistance factor and properly measure IAPV. Even for a sand pack with kwsor = 20 mD, IAPV was zero for HPAM with molecular weight (Mw) of 18 million g/mol. Fine-grain particles (<20 μm) strongly impacted polymer retention values. Native NB#1 sand with a significant component of particles <20 μm exhibited 290 μg/g, while the same sand exhibited 28 μg/g after these small particles were removed. Polymer retention did not necessarily correlate with mineral composition. The NB#1, NB#3, and OA sands had similar elemental and clay compositions, but the NB#1 sand exhibited ~10 times higher retention than the NB#3 sand. Polymer retention did not necessarily correlate with permeability. NB#1 sand exhibited much higher retention than OA sand, even though NB#1 sand was twice as permeable as OA sand. No evidence of chromatographic separation of HPAM molecular weights was found in our experiments. Although retention tended to be greater without a residual oil saturation (than at Sor), the effect was not strong. Aging a core (with high oil saturation) at 60°C reduced HPAM retention by a factor of two. As expected, under similar conditions, polymer retention was greater for a higher Mw HPAM (18 million g/mol) than for a lower Mw HPAM (10-12 million g/mol). In many cases with high polymer retention values (e.g., 240 μg/g), polymer arrival at the end of the core was relatively quick, but achieving the injected concentration occurred gradually over many pore volumes. This effect was not due chromatographic separation of polymer molecular weights. Results from modeling of this behavior was consistent with concentrationdependent polymer retention. The form assumed for the retention function in a simulator can have an important impact on the timing and magnitude of the oil response from a polymer flood. Field-based observations can underestimate polymer retention, depending on when the tracer and polymer concentrations were measured and the assumptions made about reservoir heterogeneity.

Wang, Dongmei↗
DOE OSTI2022

Connection between coherent phonons and electron-phonon coupling in Sb (111)

We report time- and angle-resolved photoemission spectroscopy (trARPES) measurements on the Sb(111) surface. We observe band- and momentum-dependent binding-energy oscillations in the bulk and surface bands driven by $A_{1g}$ and $E_{g}$ coherent phonons. While the bulk band shows simultaneous $A_{1g}$ and $E_{g}$ oscillations, the surface bands show either $A_{1g}$ or $E_{g}$ oscillations. The observed behavior is reproduced by frozen-phonon calculations based on density-functional theory. Furthermore, this evidences the connection between electron-phonon coupling and coherent binding energy dynamics tied to lattice vibration and confirms that band-, momentum-, and mode-dependent electron-phonon coupling can indeed be probed by trARPES in the low fluence limit.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Scaling and adiabaticity in a rapidly expanding gluon plasma

In this work we aim to gain qualitative insight on the far-from-equilibrium behavior of the gluon plasma produced in the early stages of a heavy-ion collision. It was recently discovered [1] that the distribution functions of quarks and gluons in QCD effective kinetic theory (EKT) exhibit self-similar “scaling” evolution with time-dependent scaling exponents long before those exponents reach their pre-hydrodynamic fixed-point values. In this work we shed light on the origin of this time-dependent scaling phenomenon in the small-angle approximation to the Boltzmann equation. We first solve the Boltzmann equation numerically and find that time-dependent scaling is a feature of this kinetic theory, and that it captures key qualitative features of the scaling of hard gluons in QCD EKT. We then proceed to study scaling analytically and semi-analytically in this equation. We find that an appropriate momentum rescaling allows the scaling distribution to be identified as the instantaneous ground state of the operator describing the evolution of the distribution function, and the approach to the scaling function is described by the decay of the excited states. That is to say, there is a frame in which the system evolves adiabatically. Furthermore, from the conditions for adiabaticity we can derive evolution equations for the time-dependent scaling exponents. In addition to the known free-streaming and BMSS fixed points, we identify a new “dilute” fixed point when the number density becomes small before hydrodynamization. Corrections to the fixed point exponents in the small-angle approximation agree quantitatively with those found previously in QCD EKT and arise from the evolution of the ratio between hard and soft scales.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

The synthesis and electrical transport properties of carbon/Cr 2 GaC MAX phase composite microwires

While MAX phases offer an exotic combination of ceramic and metallic properties, rendering them a unique class of materials, their applications remain virtually hypothetical. To overcome this shortcoming, a sol–gel based route is introduced that allows access to microwires in the range of tens of micrometers. Thorough structural characterization through XRD, SEM, EDS, and AFM demonstrates a successful synthesis of carbonaceous Cr2GaC wires, and advanced low temperature electronic transport measurements revealed resistivity behavior dominated by amorphous carbon. Here, the tunability of electronic behavior of the obtained microwires is shown by a halide post-synthesis treatment, allowing purposeful engineering of the microwires’ electrical conductivity. Raman studies revealed the polyanionic nature of the intercalated halides and a slow decrease in halide concentration was concluded from time-dependent conductivity measurements. Based on these findings, the process is considered a viable candidate for fabricating chemiresistive halogen gas sensors.

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

Excess-iron driven spin glass phase in Fe 1 + y Te 1 – x Se x

The iron-chalcogenide superconductor FeTe 1– x Se x displays a variety of exotic features distinct from iron pnictides. Although much effort has been devoted to understanding the interplay between magnetism and superconductivity near x = 0.5, the existence of a spin glass phase with short-range magnetic order in the doping range ( x ~ 0.1–0.3) has rarely been studied. Here, we use DC/AC magnetization and (quasi) elastic neutron scattering to confirm the spin-glass nature of the short-range magnetic order in a Fe 1.07 Te 0.8 Se 0.2 sample. The AC-frequency dependent spin-freezing temperature T f generates a frequency sensitivity Δ T f ( ω )/[ T f ( ω )Δlog 10 ω ] ≈ 0.028 and the description of the critical slowing down with τ = τ 0 ( T f / T SG – 1) - z v gives T SG ≈ 22 K and zv ≈ 10, comparable to that of a classical spin-glass system. We have also extended the frequency-dependent T f to the smaller time scale using energy-resolution-dependent neutron diffraction measurements, in which the T N of the short-range magnetic order increases systematically with increasing energy resolution. We find that by removing the excess iron through annealing in oxygen, the spin-freezing behavior disappears, and bulk superconductivity is realized. Thus, the excess Fe is the driving force for the formation of the spin-glass phase detrimental to bulk superconductivity.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Exact solution of two simple non-equilibrium electron-phonon and electron-electron coupled systems: The atomic limit of the Holstein-Hubbard model and the generalized Hatsugai-Komoto model

One of the challenges in many-body physics is determining the effects of phonons on strongly correlated electrons. The difficulty arises from strong correlations at differing energy scales—for band metals, Migdal-Eliashberg theory accurately determines electron-phonon coupling effects due to the absence of vertex corrections, but strongly correlated electrons require a more complex description, and the standard Migdal-Eliashberg approach does not necessarily apply. In this work, we solve for the atomic limit Green's function of the Holstein-Hubbard model with both time-dependent electron-electron and electron-phonon couplings. Then, we examine the photoemission spectra of this model in and out of equilibrium. Next, we use similar methods to exactly solve an extended version of the Hatsugai-Komoto model and examine its behavior in and out of equilibrium. These calculations lead us to propose using the first moment of the photoemission spectra to signal non-equilibrium changes in electron-electron and electron-phonon couplings.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Polarity effects on wake–up behavior of Al 0.94 B 0.06 N ferroelectrics

Wurtzite ferroelectric materials are promising candidates for energy-efficient memory technologies, particularly for applications requiring high operating temperatures. Asymmetric wake-up behaviors, in which the polarization reversal depends both on polarity and cycle number for the first few dozen cycles, must be better understood for reliable device operation. Here, the detailed analysis of the asymmetric wake-up behavior of thin film Al 0.94 B 0.06 N was performed combining time-resolved switching measurements with Rayleigh analysis, piezoelectric measurements, and etching experiments of progressively switched samples. Furthermore, the analysis shows that the gradual opening of the polarization hysteresis loops associated with wake-up is driven by a gradual increase in the domain-wall density and/or domain-wall mobility with electric field cycle to the polarity opposite to the growth polarity. The insights of this discovery will help to guide interface and polarity design in the eventual deployment of reliable devices based on these materials.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Spatiotemporal Plasma-Particle Characterization of Dry Aerosols Using Nanosecond, Femtosecond, and Filament Laser-Produced Plasmas

The ability to rapidly characterize dry aerosols in air using laser-induced breakdown spectroscopy (LIBS) with femtosecond laser pulses promises advancement towards real-time atmospheric sampling and standoff capabilities. Of particular interest is the ability to apply LIBS in the context of low-particle loaded environments where discrete particle interactions must be observed within the sampling volume of the laser-produced plasma (LPP). In this study, dry nanoparticles in suspension are generated from a standard solution and sampled in air using Q-switched nanosecond (ns-) pulses, short-focus (SF) femtosecond (fs-) pulses, and filaments. Short time-gated plasma images are captured to observe spatially and temporally varying discrete plasma-particle interactions, which is shown to influence early air breakdown behavior and subsequent plasma evolution. Along with images, photo-multiplier tube (PMT) measurements are conducted where strong spatiotemporal dependencies are exhibited by the collected emission signal on particle proximity and plasma expansion behavior. Finally, conditional analysis is performed on LIBS measurements to determine associated sampling probabilities and filter out spectra with poor or absent emission peaks with an adaptive threshold algorithm.

42 ENGINEERING↗
DOE OSTI2022

Behavior and mechanisms of Doppler wind lidar error in varying stability regimes

Abstract. Wind lidars are widespread and important tools in atmospheric observations. An intrinsic part of lidar measurement error is due to atmospheric variability in the remote-sensing scan volume. This study describes and quantifies the distribution of measurement error due to turbulence in varying atmospheric stability. While the lidar error model is general, we demonstrate the approach using large ensembles of virtual WindCube V2 lidar performing a profiling Doppler-beam-swinging scan in quasi-stationary large-eddy simulations (LESs) of convective and stable boundary layers. Error trends vary with the stability regime, time averaging of results, and observation height. A systematic analysis of the observation error explains dominant mechanisms and supports the findings of the empirical results. Treating the error under a random variable framework allows for informed predictions about the effect of different configurations or conditions on lidar performance. Convective conditions are most prone to large errors (up to 1.5 m s−1 in 1 Hz wind speed in strong convection), driven by the large vertical velocity variances in convective conditions and the high elevation angle of the scanning beams (62∘). Range-gate weighting induces a negative bias into the horizontal wind speeds near the surface shear layer (−0.2 m s−1 in the stable test case). Errors in the horizontal wind speed and direction computed from the wind components are sensitive to the background wind speed but have negligible dependence on the relative orientation of the instrument. Especially during low winds and in the presence of large errors in the horizontal velocity estimates, the reported wind speed is subject to a systematic positive bias (up to 0.4 m s−1 in 1 Hz measurements in strong convection). Vector time-averaged measurements can improve the behavior of the error distributions (reducing the 10 min wind speed error standard deviation to <0.3 m s−1 and the bias to <0.1 m s−1 in strong convection) with a predictable effectiveness related to the number of decorrelated samples in the time window. Hybrid schemes weighting the 10 min scalar- and vector-averaged lidar measurements are shown to be effective at reducing the wind speed biases compared to cup measurements in most of the simulated conditions, with time averages longer than 10 min recommended for best use in some unstable conditions. The approach in decomposing the error mechanisms with the help of the LES flow field could be extended to more complex measurement scenarios and scans.

17 WIND ENERGY↗
DOE OSTI2025

Nanoscale Magnetic Ordering Dynamics in a High Curie Temperature Ferromagnet

Thermally driven transitions between ferromagnetic and paramagnetic phases are characterized by critical behavior with divergent susceptibilities, long-range correlations, and spin dynamics that can span kHz to GHz scales as the material approaches the critical temperature T c , but it has proven technically challenging to probe the relevant length and time scales with most conventional measurement techniques. In this study, we employ scanning nitrogen-vacancy center based magnetometry and relaxometry to reveal the critical behavior of a high-T c ferromagnetic oxide near its Curie temperature. Cluster analysis of the measured temperature-dependent nanoscale magnetic textures points to a 3D universality class with a correlation length that diverges near T c . Meanwhile, the temperature-dependent spin dynamics, measured through all optical relaxometry suggest that the phase transition is in the XY universality class. In conclusion, our results capture both static and dynamic aspects of critical behavior, providing insights into universal properties that govern phase transitions in magnetic materials.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Organismic materials for beyond von Neumann machines

The elementary basis of intelligence in organisms with a central nervous system includes neurons and synapses and their complex interconnections forming neural circuits. In non-neural organisms such as slime mold with gel-like media, viscosity modulation enables adaptation to changing environments. At a larger scale, collective intelligence emerges via social interactions and feedback in animal colonies. Learning and memory are therefore multi-scale features that evolve as a result of constant interactions with the environment. There is growing interest in emulating such features of intelligence in computing machines and autonomous systems. Materials that can respond to their environment in a manner similar to organisms (referred to as "organismic materials") therefore may be of interest as hardware components in artificial intelligence machines. In this brief review, we present a class of semiconductors called correlated oxides as candidates for learning machines. The term "correlated" refers to the fact that electrons in such lattices strongly interact and the ground state is not what is predicted by classical band theory. Such materials can undergo insulator-metal transitions at near ambient conditions under external stimuli such as thermal or electrical fields, strain, and chemical doping. Depending on the mechanism driving the transition, intermediate states can be metastable with different volatilities, and the time scales of phase change can be controlled over many orders of magnitude. The change in electronic properties can be sharp or gradual, leading to digital or analog behavior. These properties enable the realization of artificial neurons and synapses and emulate the associative and non-associative learning characteristics found in various organisms. We examine microscopic properties concerning electronic and structural transitions leading to collective behavior and theoretical treatments of the ground state and dynamical response, showcasing VO 2 as a model system. Next, we briefly review algorithms designed from the plasticity demonstrated by phase changing systems. Finally, we conclude the brief review with suggestions for future research toward realizing non-von Neumann machines.

42 ENGINEERING↗
DOE OSTI2020

Numerical methodology for treating static and dynamic dislocation problems near a free surface

Abstract Simulation of dislocation dynamics enables researchers and scientists to explore the plastic behavior of crystalline materials under loading. Analytic solutions for the stress field due to a linear dislocation segment near a free surface are case-specific, e.g. dealing with either a horizontal segment or a vertical segment, and therefore hard to implement in time-dependent dislocation dynamics simulations as different dislocation segments could be differently oriented. This article presents a generalized numerical framework to find the stress field beneath a free surface due to the presence of a dislocation segment. The framework can be expanded to non-flat surfaces. Also, three-dimensional discrete dislocation dynamics simulations are performed here, which clearly show the effect of free surfaces on the flow stress of a material.

Siddique, Abu Bakar↗
DOE OSTI2025

Symmetry instability induced by topological phase transitions

The symmetry-topology interplay dictates how to define order parameters and classify material ordered phases. However, current understanding of this interplay has been predominately approached from a one-sided perspective, with topological states being classified within the constraints imposed by specific fixed symmetries. Here, in this study, we complete this full circle by demonstrating spontaneous symmetry breaking that results from a periodic alteration of topological phases induced by light in a centrosymmetric Dirac material ZrTe 5 . The distinguishing feature is the observation of robust correlation and striking anomalies in the fluence and temperature dependence of key transport parameters. First, both shift current 𝐽 𝑠 and displacement current 𝐽 𝑑 , arising from interband transition and infrared phonon driving, respectively, along with charge carrier pumping, exhibit similar behaviors. Second, they all peak at similar low pump fluence, followed by a subsequent reduction as the fluence further increases. This behavior cannot be explained by conventional energetically allowed, direct excitations. Third, all the three observables exhibit anomalies when they approach the topological phase transition temperature. These results highlight the unique low-energy pumping behaviors in ZrTe 5 , characterized by reversible fluence dependence and a “hinge-like” interaction that connects various electronic and lattice observables, including phonons, charge carriers, and currents. Our findings, supported by model analysis, provide key insights into the fragility of crystalline (inversion) and time-reversal symmetries during the dynamics of topological phase transitions. This fragility drives spontaneous symmetry breaking, evidenced by the synchronized emergence of off-resonant infrared phonons and broken-symmetry photocurrents.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

High-Fidelity Modeling of Fuel-To-Coolant Thermomechanical Transport Behaviors Under Transient Conditions

This report summarizes the work completed under NEUP project number 21-24006. The objectives of this project are to advance the high-fidelity modeling capabilities and important phenomena that is important for high-burnup UO 2 and accident tolerant fuels (ATF) during transient conditions. Accurate modeling of the time-dependent phenomena that impact material performance must be used to determine the figures of merit and safety margin. Phenomena such as fuel fragmentation, cladding oxidation, pellet-clad interaction, clad ballooning, and clad rupture are examples that pose challenges to modeling during these transients. This project focused on the development of high-fidelity tightly coupled multiphysics models that can capture the time-dependent material response and associated thermal hydraulic conditions during these events. These models can then be validated against existing separate effects tests and in-pile integral experiments and will be used to model Transient Reactor Test facility (TREAT) loss-of-coolant accidents (LOCA) experiments. To achieve the project objective, we used a combination of NEAMS and NRC codes to model various LOCA test sets for the separate effects and in-pile integral experiments. BlueCRAB tool set, which can accurately predict material response at a sub-fuel pin level, as well as modeling the entire reactor system response to these events. Fuel performance was modeled using BISON (various versions) and FAST (version 1.2.1). BISON and FAST can model on a sub-fuel pin level the fuel performance under transient conditions. Both have simplified thermal hydraulic models that are capable of providing basic coolant boundary conditions. To better capture the thermomechanical interaction between the fuel, clad, and coolant, more sophisticated thermal hydraulic models are necessary.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

A Theoretical Analysis on the Dynamic Behavior and Electromechanical Response of PZT

Piezoelectric materials are used as a power source in stockpile components for safety and reliability assurances. The objective of this project is to gain insights into the fundamental pressure-induced behavior and electromechanical response of lead zirconate titanate (PZT). Specifically, to establish the basis for an accurate, physics-based material model for robust model-based design to rapidly optimize PZT-based materials and components for performance studies. Established ab-initio methods are used to interrogate and understand the dynamic behavior of PZT as a function of composition (50/50, 65/35, 80/20, 95/5) and dopants (La, Nb) to overcome the costly and time-consuming experimental methodologies. New cold curves for pure and doped single-crystal PZT are obtained as the reference equation-of-state (EOS). A negligible change in the pressure responses was observed for the systems and strains studied. Dielectric and piezoelectric responses of pure and doped single-crystal PZT were also calculated as a function of pressure. For undoped PZT, there is a clear and orderly pressure and composition dependence. For doped PZT, there is a significant increase in the responses, but the behavior is very disordered and inconclusive.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Elucidating Cycling Rate-Dependent Electrochemical Strains in Sodium Iron Phosphate Cathodes for Na-ion Batteries

Battery electrodes materials undergo significant mechanical instabilities which affects their longevity and exert rate-limitations during the cycling process. In this study, we investigate the rate-dependent mechanical response of sodium iron phosphate (NaFePO4, NFP) cathodes during Na intercalation via galvanostatic cycling at different rates by employing digital image correlation, electrochemical methods, and mathematical model. The mechanical behaviour of the electrode shows strong dependance on the applied scan rate. At slower rates, electrode shows asymmetrical strain generation between anodic and cathodic cycles, which is attributed to the formation of cathode-electrolyte interface layers. The electrode undergoes smaller strain generation when cycled at slower rates when the same amount of Na ions is removed or inserted into the electrode. A mathematical model was developed to predict strain evolution in the composite electrode as well as the concentration profile of the Na ions in the electrode particles. Rate-dependent and time dependent factors on the strain generation in the electrode are attributed to the capacity-dependent intercalation strains, rate-dependent mismatch strains, and time-dependent irreversible strains. The combination of in situ strain measurements with the analytical model provided new insight into the electrochemically induced mechanical deformations in Na-ion cathode electrodes.

Ozdogru, Bertan↗
DOE OSTI2025

(NH 3 (CH 2 ) 7 NH 3 ) 2 Sn 3 I 10 , a Vacancy-Ordered Three-Dimensional Tin(II) Perovskite-Derived Semiconductor

Ordering vacancies in hybrid Sn(II) halide semiconductors provides a strategy for preventing uncontrolled oxidation and formation of mobile holes. In this study, we report the structure and optical and electronic properties of (NH 3 (CH 2 ) 7 NH 3 ) 2 Sn 3 I 10 , a vacancy-ordered perovskite derivative with three-dimensional inorganic connectivity. The crystal structure resembles that of a Dion–Jacobson layered perovskite derivative, but with [SnI 5 ] square pyramids bridging the layers. UV–vis diffuse reflectance spectroscopy reveals a sharp onset of light absorption at 1.86(1) eV with the photoluminescence emission maximum at 1.90(1) eV. However, the maximum excitation occurs from 3.42 to 3.81 eV (325 to 370 nm), revealing a significant Stokes shift of 1.3 eV. The electronic properties determined from dark and time-resolved microwave conductivity measurements reveal a minimum carrier mobility of 4.3 × 10 –2 cm 2 V –1 s –1 and a maximum carrier density of 5.96 × 10 16 cm –3 , a uniquely low value for a hybrid Sn(II) halide semiconductor. The transport behavior in combination with first-principles calculations of the electronic band structure and dielectric permittivity suggest polaron-mediated electronic transport, yet the photogenerated carriers have a fast and fluence-dependent nonradiative recombination rate, suggestive of localized “defect-like” states at the band edge. The observed photoluminescence is most consistent with single-ion-like behavior of an asymmetric Sn(II) environment. Together, these results suggest that defect ordering presents a strategy for the reduction of mobile charge carriers at equilibrium.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗