Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “STRENGTH CALCULATIONS”

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 523 records · Page 29

Electronic structure of Co 3 d states in the Kitaev material candidate honeycomb cobaltate Na 3 Co 2 SbO 6 probed with x-ray dichroism

The recent prediction that honeycomb lattices of Co 2+ (3d 7 ) ions could host dominant Kitaev interactions provides an exciting direction for exploration of new routes to stabilizing Kitaev’s quantum spin liquid in real materials. Na 3 Co 2 SbO 6 has been singled out as a potential material candidate provided that spin and orbital moments couple into a J eff = $\frac{1}{2}$ ground state, and that the relative strength of trigonal crystal field and spin-orbit coupling acting on Co ions can be tailored. Using x-ray linear dichroism (XLD) and x-ray magnetic circular dichroism (XMCD) experiments, alongside configuration interaction calculations, we confirm the counterintuitive positive sign of the trigonal crystal field acting on Co 2+ ions and test the validity of the J eff = $\frac{1}{2}$ description of the electronic ground state. Furthermore, the results lend experimental support to recent theoretical predictions that a compression (elongation) of CoO 6 octahedra along (perpendicular to) the trigonal axis would drive this cobaltate toward the Kitaev limit, assuming the J eff = $\frac{1}{2}$ character of the electronic ground state is preserved.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Low spin spectroscopy of neutron-rich Cl 43 , 44 , 45 via β - and β n decay

In this article, β - decay studies of neutron-rich isotopes 43,45 S performed at the National Superconducting Cyclotron Laboratory are reported. β -delayed $\gamma$ transitions were detected by an array of 16 clover detectors surrounding the Beta Counting Station, which consists of a 40×40 double-sided silicon strip detector followed by a single-sided silicon strip detector. β -decay half-lives were extracted for 43,45 S by correlating implants and decays in the pixelated implant detector and by considering further coincidences with $\gamma$ transitions in the daughter nuclei. Further, the level structures of 43,45 Cl are expanded by the addition of 20 and 8 new $\gamma$ transitions in 43 Cl and 45 Cl respectively, and core-excited negative-parity states were observed in both nuclei for the first time. For 45 S, a large fraction of the β -decay strength was observed feeding neutron-unbound states in 45 Cl, which, decaying by delayed neutrons, populated excited states in the β ⁢1⁢ n daughter, 44 Cl. Experimental observations were compared to detailed shell-model calculations using the SDPFSDG-MU interaction to highlight the role of the diminished N = 28 neutron shell gap and the near degeneracy of the proton s 1/2 and d 3/2 orbitals in the structure of the neutron-rich Cl isotopes. The current work also provides further support to a ground-state spin-parity assignment of 3/2 + in 45 Cl.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Phenomenology of Spectral Functions in Disordered Spin Chains at Infinite Temperature

Studies of disordered spin chains have recently experienced a renewed interest, inspired by the question to which extent the exact numerical calculations comply with the existence of a many-body localization phase transition. For the paradigmatic random field Heisenberg spin chains, many intriguing features were observed when the disorder is considerable compared to the spin interaction strength. Here, we introduce a phenomenological theory that may explain some of those features. The theory is based on the proximity to the noninteracting limit, in which the system is an Anderson insulator. Taking the spin imbalance as an exemplary observable, we demonstrate that the proximity to the local integrals of motion of the Anderson insulator determines the dynamics of the observable at infinite temperature. In finite interacting systems our theory quantitatively describes its integrated spectral function for a wide range of disorders.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Repurposing the θ (Al 2 Cu) phase to simultaneously increase the strength and ductility of an additively manufactured Al–Cu alloy

The Al–Cu–Mn–Zr (ACMZ) cast family of alloys offers unique properties and value propositions for higher strength, higher temperature lightweight components of future vehicles. Earlier research has demonstrated trade-offs in the selection of the alloy chemistry in which an increase in Cu content from 6 up to 9 wt% improves hot tear resistance but lowers ductility. However, a recent study has demonstrated that higher-Cu (9Cu) ACMZ fabricated with laser powder bed fusion additive manufacturing (AM) results in an increase in both ductility and strength when compared to as-aged microstructure of cast 9Cu alloys. The mechanisms of differing mechanical performance of the cast and AM ACMZ alloys are elucidated in the current paper through the utilization of in situ high energy x-ray diffraction (HEXRD) tensile testing wherein lattice strains of different phases are calculated and correlated to their stresses. The ACMZ alloys consisted of theta (θ) and theta prime (θ') phases (both Al 2 Cu in nominal composition) within an aluminum matrix. The larger micron-size θ phase which decorated the grain boundaries in 9Cu ACMZ cast alloys recorded small lattice strains, while the submicron, homogeneously distributed θ phase in the 9Cu ACMZ AM alloy recorded considerably higher lattice strains. The maximum stress reached in the θ phase for the cast 9Cu alloy was found to be ~280 MPa, which was lower than the AM 9Cu alloy which registered a maximum stress of ~1.4 GPa. These measurements indicate that delayed fracture of the finer intermetallic phases simultaneously improves the ductility and strength of AM 9Cu alloy relative to the cast 9Cu alloy, which exhibits early fracture of the larger intermetallic particles.

36 MATERIALS SCIENCE↗

Hybridization effect on the x-ray absorption spectra for actinide materials: Application to PuB 4

Studying the local moment and 5 f -electron occupations sheds insight into the electronic behavior in actinide materials. X-ray absorption spectroscopy (XAS) has been a powerful tool to reveal the valence electronic structure when assisted with theoretical calculations. However, the analysis currently taken in the community on the branching ratio of the XAS spectra generally does not account for the hybridization effects between local f orbitals and conduction states. In this paper, we discuss an approach which employs the density functional theory plus Gutzwiller rotationally invariant slave boson method to obtain a local Hamiltonian for the single-impurity Anderson model, and calculates the XAS spectra by the exact diagonalization (ED) method. A customized numerical routine was implemented for the ED XAS part of the calculation. By applying this technique to the recently discovered 5 f -electron topological Kondo insulator Pu B 4 , we determined the signature of 5 f -electronic correlation effects in the theoretical x-ray spectra. Furthermore, we found that the Pu 5 f - 6 d hybridization effect provides an extra channel to mix the j = 5 / 2 and 7 / 2 orbitals in the 5 f valence. As a consequence, the resulting electron occupation number and spin-orbit coupling strength deviate from the intermediate-coupling regime.

36 MATERIALS SCIENCE↗

Reversible Modification of Rashba States in Topological Insulators at Room Temperature by Edge Functionalization

Quantum materials with novel spin textures from strong spin-orbit coupling (SOC) are essential components for a wide array of proposed spintronic devices. Topological insulators have a necessary strong SOC that imposes a unique spin texture on topological states and Rashba states that arise on the boundary, but there is no established methodology to control the spin texture reversibly. Here, it is demonstrated that functionalizing Bi 2 Se 3 films by altering the step-edge termination directly changes the strength of SOC and thereby modifies the Rashba strength of 1D edge states. Scanning tunneling microscopy/spectroscopy shows that these Rashba edge states arise and subsequently vanish through the Se functionalization and reduction process of the step edges. The observations are corroborated by density functional theory calculations, which show that a subtle chemical change of edge termination fundamentally alters the underlying electronic structure. Importantly, fully reversible and repeatable switching of Rashba edge states across multiple cycles at room temperature is experimentally demonstrated. The results imply Se functionalization as a practical method to control SOC and spin texture of quantum states in topological insulators.

Rashba edge states↗

BISON validation of FeCrAl cladding mechanical failure during simulated reactivity-initiated accident conditions*

Here, a reactivity-initiated accident (RIA) is a postulated design basis accident in light water reactors (LWRs) in which a rapid reactivity insertion induces a fission rate increase and a fuel pellet temperature rise. During RIA, the fuel pellet thermally expands and may cause pellet-cladding mechanical interaction (PCMI). Separate effects PCMI tests were performed on C26M FeCrAl cladding tube samples, introducing biaxial stress via a well understood modified burst test (MBT) system. A high-speed camera in the MBT system captured the projections, covering a 360° view of the cladding deformation and enabling a digital image correlation (DIC) method to quantify the surface strains with high fidelity. Representative hot zero-power RIA mechanical loading conditions were applied to the sample, and the test duration ranged from 20 to 500 ms at an average temperature of 573 K. BISON finite element–based fuel performance code modeling was performed against the high-fidelity DIC data produced from the MBTs. Validation calculations were conducted with 2D models and systematically compared with test data of the burst time, burst pressure, burst hoop strains, and hoop strain rates. Based on the behaviors from the separate effects test, BISON calculations satisfactorily predicted the cladding deformation behaviors. Sensitivity analysis was conducted to identify highly influential mechanical properties responsible for the cladding failure behavior during the MBT experiments. The results highlight the significance of the cladding's mechanical strength in governing cladding strain, followed by the significance of the mechanical interactions between the pellet and the cladding.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Structure of high-lying levels populated in the 96 Y → 96 Zr β decay

The nature of the high-lying final levels of the 96 Y gs β decay, one of the three most important contributors to the high-energy reactor antineutrino spectrum, has been investigated in high-resolution γ-ray spectroscopy following the β decay as well as in a campaign of inelastic photon scattering experiments. The comprehensive approach establishes 1 - levels associated with the Pygmy Dipole Resonance as high-lying final levels in the β decay. Branching ratios extracted from β decay complement photon scattering and allow the absolute E1 excitation strength to be determined for levels populated in both reactions. The combined data represents a comprehensive approach to the wavefunction of the 1 - levels below the Qβ value, which are investigated in the Quasiparticle Phonon Model. The calculations reveal that the components populated in β decay contribute only with small amplitudes to the complex wavefunction of these 1 - levels. A comparison of the β decay results to data from total absorption γ -ray spectroscopy demonstrates a good agreement between both measurements.

Pygmy dipole resonance↗

Luther-Emery liquid and dominant singlet superconductivity in the hole-doped Haldane spin-1 chain

We investigate the pairing tendencies in the hole-doped Haldane spin-1 chain. To allow for doping, we extend the original spin chain Hamiltonian into a fermionic model involving a two-orbital Hubbard chain at intermediate or strong repulsive interaction strengths U and for degenerate orbitals. At half filling and large U, the ferromagnetic Hund's coupling, J H , generates effective spin-1 moments, with antiferromagnetic correlations between sites. Using large-scale density matrix renormalization group calculations, we accurately study the system's behavior under light hole-doping. For U = 1.6 in units of the noninteracting bandwidth and for J H /U ≳ 0.275, we find that singlet pairing dominates the long-distance physics, establishing this system as a promising platform for repulsively mediated superconductivity. We provide concrete examples of materials that could realize the physics described here. We also provide evidence that the system approaches a Luther-Emery liquid state at large system sizes, reminiscent of the behavior of doped one-orbital two-leg ladders at weak coupling, which also have superconducting tendencies. The numerically calculated central charge approaches one in the thermodynamic limit, indicating a single gapless mode as is expected for the Luther-Emery state. Exponents characterizing the power-law decays of singlet pair-pair and charge density-density correlations are determined, and found to approximately satisfy the Luther-Emery identity.

1-dimensional systems↗

Shell model analysis of the B ( E 2,2 + → 0 + ) values in the A = 70, T = 1 triplet 70 Kr, 70 Br, and 70 Se

The B(E2, 2 + → 0 + ) transition strengths of the T=1 isobaric triplet 70 Kr, 70 Br, 70 Se, recently measured at the RIKEN Radioactive Isotope Beam Factory (RIBF), are discussed in terms of state-of-the-art large scale shell model calculations using the JUN45 and JUN45 + LNPS plus Coulomb interactions. In this Letter we argue that, depending on the effective charges used, the calculations are either in line with the experimental data within statistical uncertainties, or the anomaly happens in 70 Br, rather than 70 Kr. In the latter case, we suggest that it can be due to the presence of a hitherto undetected 1 + T = 0 state below the yrast 2 + T = 1 state. Furthermore, our results do not support a shape change of 70 Kr with respect to the other members of the isobaric multiplet.

59 ≤ A ≤ 89↗

Rapid Bayesian High Entropy Alloy Designs Fabricated via Wire Arc Additive Manufacturing

Purpose: This project seeks to demonstrate a new high-throughput (rapid) alloy design technique applied to creating new high entropy alloys (HEAs) for extreme environments. High entropy alloys shift the design paradigm from being focused on a single principal element (e.g. nickel-based alloys) to target alloys that include high atomic fractions (X >10%) of multiple elements. These HEA materials can exhibit sluggish diffusion and enhanced corrosion resistance, ideal for potential applications in advanced ultra supercritical (A-USC) steam cycles for power generation. Scope: The addition of multiple elements in high atomic fractions creates an enormous design space that cannot easily be investigated by traditional material design strategies such as designed of experiments (DOE). This project utilizes a Bayesian machine learning algorithm that has been modified to work with calculation of phase diagrams (CALPHAD) software. This Bayesian algorithm reduces manual inputs and increase the likelihood of achieving an optimal solution. Compositional inputs to this algorithm will be assessed using existing material property models for high temperature strength and corrosion resistance. The target for alloy performance will be a 15% (~100 ⁰C) increase in allowable service temperature beyond heat-resistant stainless steels while maintaining or improving alloy cost and corrosion resistance. Haynes 230 was selected as a baseline, which is 57 wt% Ni with 22 wt% Cr 14 wt% W, and 2 wt% Mo as solid solution strengtheners. In addition to rapid design via Bayesian machine learning, the alloys were rapidly fabricated using a multi-wire arc additive manufacturing (mWAAM) technique which allows for precise control of alloy composition and assessing of alloy design “windows” to study composition effects. Build speeds for wire-arc additive processes are among the highest for additive technologies enabling rapid and reliable sample fabrication when compared to conventional methods such as arc button melting. The mWAAM samples will be rapidly characterized via instrumented indentation for room temperature modulus and strength and for elevated temperature strength via hot hardness tests. After being screened with hardness testing, potential alloys will be further evaluated with conventional microscopy techniques including scanning electron microscopy (SEM) and transmission electron microscopy (TEM) to assess agreement with modeling results. The most promising compositions will also be evaluated by printing full sized tensile specimens for mechanical behavior tests at elevated temperatures. Results: Bayesian machine learning of a single performance function was initially used to optimize five performance metrics: 1) single phase stability, 2) yield strength, 3) creep resistance (low diffusion coefficient), 4) freezing range (weldability), and 5) material cost. The single performance function was suboptimal as assumptions had to be made about the results while formulating the optimization. A goal-oriented Bayesian optimization strategy (Hanaoka, 2021) was implemented with CALPHAD for use with the five metrics above. This multi-objective Bayesian optimization (MOBO) enabled the design of NiCrCoFe alloys with V and W additions. A base composition of NiCoCr was selected as Ni provides a stable FCC matrix, Cr aids corrosion/oxidation resistance, and Co is a solid-solutions strengthener that also improves creep by increasing the activation energy. Fe helps reduce diffusion coefficients and cost. Finally, V and W were selected for their reasonable solubility and high atomic misfit to aid in solid solution strengthening. Cracking of the mWAAM specimens was an early issue, and the Easton solidification cracking model (Easton et al., 2014a) was selected for addition to the MOBO function. High performing alloys fabricated by mWAAM included Ni 28 Cr 25 Co 26 Fe 15 V 8 and Ni 62 Cr 18 Co 1 Fe 3 W 15 . It was observed that even after adapting the mWAAM process for W, the W did not fully dissolve. To fully evaluate the Ni 62 Cr 18 Co 1 Fe 3 W 15 composition, a cored wire (80-20 NiCr sheath/powder core) was manufactured and printed via WAAM, and HIP’ing was utilized to homogenize and densify the printed alloy. The V and W alloys produced met metrics 1 (solid solution), 4 (solidification cracking), and 5 (cost). However, an unmodeled mechanism of thermal stress cracking was identified in the WAAM produced materials, perhaps exacerbated by the lack of grain boundary strengthening elements (B, C). Conclusions & Recommendations: A high-throughput (rapid) alloy design technique was applied to designing and manufacturing new high entropy alloys (HEAs) for extreme environments utilizing MOBO and mWAAM. The developed process was rapid and effective in addressing the mechanisms included in the model. The lack of grain boundary strengthening element additions (e.g., B, C) was a simplification that likely produced thermal stress cracking that turned into a large part of the investigation. Additions on the order of 0.005 wt% B and 0.05 wt% C likely would have minimized thermal stress grain boundary cracking. Overall, the high throughput design strategy is promising for rapid design of metrics-driven alloys for advanced ultra supercritical (A-USC) steam cycles for power generation. The MOBO and mWAAM process could be commercialized to accelerate metrics-driven alloy design. In addition, the cored-wire process utilized for scale-up is a promising high-volume process for WAAM alloy development and scale-up.

36 MATERIALS SCIENCE↗

Continuum contribution to charged-current absorption of low-energy $ν_e$ on $^{40}$Ar

Accurate modeling of the absorption of tens-of-MeV $ν_e$ on $^{40}$Ar is needed to enable measurements of astrophysical neutrinos using large liquid argon time projection chamber (LArTPC) detectors, such as those planned for the Deep Underground Neutrino Experiment (DUNE). We revisit the MARLEY neutrino interaction model used in present estimates of DUNE sensitivity to supernova and solar neutrino signals. Multiple theoretical refinements are pursued, especially in the unbound continuum region of nuclear excitation energy. Inclusive charged-current neutrino-argon cross sections are calculated using a hybrid strategy. Nuclear transitions to unbound states are treated using a Hartree-Fock Continuum Random Phase Approximation (HF-CRPA) model, including forbidden contributions. Allowed transitions to low-lying discrete levels are also included using indirect measurements and approximate corrections for the momentum transfer dependence. Exclusive predictions are obtained by coupling these calculations with a statistical nuclear de-excitation model. The impact on observables of interest for DUNE and similar experiments is examined in terms of both total and differential cross sections. Our refined calculations predict a lower allowed portion of the cross section relative to the prior MARLEY model. At neutrino energies appreciably below 100 MeV, the inclusion of forbidden transitions does not fully compensate for the loss of allowed strength. For a representative neutrino burst from a galactic core-collapse supernova, our results suggest that MARLEY 1.2.0 overestimates the event yield in a DUNE-like detector by approximately 20%. However, because this overestimation is more severe at backwards angles, use of the charged-current $ν_e$-$^{40}$Ar reaction for supernova pointing may be more feasible than previously expected.

Gardiner, Steven [Fermilab]↗

Domain Separated Density Functional Theory for Reaction Energy Barriers and Optical Excitations

We recently proposed Domain Separated Density Functional Theory (DS-DFT), a framework that allows for the combination of different levels of theory for the computation of the electronic structure of molecules. Here, this work discusses the application of DS-DFT to the computation of transition-state energy barriers and optical absorption spectra. We considered several hydrogen abstraction reactions and optical spectra of molecule/metal cluster systems, including the absorption of individual species such as carbon monoxide, methane, and molecular hydrogen to a Li6 cluster. We present and discuss two domain-separated methods: i), the screened-density approximation (SDA), and, ii), linearly weighted exchange (LWE). We find that SDA, which is applied as a hybridization based on atomic domains, could be useful to computing energy barriers, whereas LWE is suited for the analysis of electronic properties such as ground-state gaps, excitation energies, and oscillator strengths.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Multiphase tin equation of state using density functional theory

In this work, we perform density functional theory (DFT) calculations of five solid phases and the liquid phase of tin. The calculations include cold curves of the five solid phases, phonon calculations in the quasiharmonic approximation over a range of volumes for each solid phase, and DFT-based molecular dynamics (DFT-MD) simulations of the liquid phase, including those of the melt curve using the Z method. Using the DFT results, we construct a tabular multiphase sesame equation of state for tin, referred to as sesame 2162. Comparisons to experimental data are made and show a high level of agreement in isobaric data, isothermal data, shock data, and phase boundary measurements, including measurements of the melt curve. The 2162 EOS will be useful for hydrodynamics simulations and has been designed with an eye toward hydrodynamics simulations that incorporate materials strength models and allow for modeling of the kinetics of phase transitions.

3-dimensional systems↗

Assessment of errors in analytic modeling of permanent magnet electron spectrometers for laser-plasma accelerators

Here, we assess the error in three treatments of a critical electron beam diagnostic for laser-plasma acceleration (LPA) experiments: a permanent magnet electron spectrometer. Since LPA electron spectrometers are often difficult to calibrate due to the scarcity of well-characterized, tunable electron beam sources in the appropriate energy range and the mechanical complexity of electron spectrometers, the standard of calibration is Hall probe measurements. We first compare the electron spectrometer performance between SIMION calculations and the Hall probe measurements. We find up to a 60% (<8%) error in determining the absolute energy for electrons below (above) 80 MeV when using the SIMION-modeled magnetic field vs the measured magnetic field due to SIMION overpredicting fringe field strength. The difference in spectrometer energy resolution is within ∼10% between the modeled and measured fields. We then assess a simple block model, which is commonly used. We find that the simple block model presented here sufficiently captures the predicted as-built spectrometer performance for any application provided that the user can tolerate uncertainties on absolute energy determination up to 5% and errors on energy resolution up to 1%.

Electromagnetic optics↗

Electric dipole polarizability of 58 Ni

The electric dipole strength distribution in 58 Ni between 6 and 20 MeV has been determined from proton inelastic scattering experiments at very forward angles at RCNP, Osaka. The experimental data are rather well reproduced by quasiparticle random-phase approximation calculations including vibration coupling, despite a mild dependence on the adopted Skyrme interaction. They allow an estimate of the experimentally inaccessible high-energy contribution above 20 MeV, leading to an electric dipole polarizability α D ⁡ ( 58 Ni) = 3.48 (31)⁢ fm 3 . This serves as a test case for recent extensions of coupled-cluster calculations with chiral effective field theory interactions to nuclei with two nucleons on top of a closed-shell system.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Machine Learning Vacancy Formation Energy in Nickel-Based Superalloys

Creep performance plays a key role in nickel-based superalloys for high temeprature applications. Creep behavior depends on many parameters such as strength, dislocations, diffusivity, and microstructural stability in addition to temeprature, applied stress, and oxidation. This work focuses on predicting vacancy formation energy in nickel-based superalloys using machine learning approach. High-throughput density functional theory (DFT) calculations are performed on Ni-based alloys with the addition of various alloying elements to predict the vacancy formation energy and vacancy concentration. Machine learning is performed using various models including graph neural networks.

creep performance↗

Carbonic Anhydrase Enhanced UV-Crosslinked PEG-DA/PEO Extruded Hydrogel Flexible Filaments and Durable Grids for CO 2 Capture

In this study, poly (ethylene glycol) diacrylate/poly (ethylene oxide) (PEG-DA/PEO) interpenetrating polymer network hydrogels (IPNH) were extruded into 1D filaments and 2D grids. The suitability of this system for enzyme immobilization and CO 2 capture application was validated. IPNH chemical composition was verified spectroscopically using FTIR. The extruded filament had an average tensile strength of 6.5 MPa and elongation at break of 80%. IPNH filament can be twisted and bent and therefore is suitable for further processing using conventional textile fabrication methods. Initial activity recovery of the entrapped carbonic anhydrase (CA) calculated from esterase activity, showed a decrease with an increase in enzyme dose, while activity retention of high enzyme dose samples was over 87% after 150 days of repeated washing and testing. IPNH 2D grids that were assembled into spiral roll structured packings exhibited increased CO 2 capture efficiency with increasing enzyme dose. Long-term CO 2 capture performance of the CA immobilized IPNH structured packing was tested in a continuous solvent recirculation experiment for 1032 h, where 52% of the initial CO 2 capture performance and 34% of the enzyme contribution were retained. These results demonstrate the feasibility of using rapid UV-crosslinking to form enzyme-immobilized hydrogels by a geometrically-controllable extrusion process that uses analogous linear polymers for both viscosity enhancement and chain entanglement purposes, and achieves high activity retention and performance stability of the immobilized CA. Potential uses for this system extend to 3D printing inks and enzyme immobilization matrices for such diverse applications as biocatalytic reactors and biosensor fabrication.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗