Search NASA⌕ Search

SEARCH · Search NASA

Results for “parallel clustering”

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 541 records · Page 30
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

Results from the CBC3 readout ASIC for CMS 2S-modules

The CBC3 is the latest version of the CMS Binary Chip for readout of the outer radial region of the upgraded CMS Tracker at the High Luminosity LHC. This 254-channel, 130 nm CMOS ASIC is designed to be bump-bonded to a substrate to which sensors will be wire-bonded. It will instrument double-layer 2S-modules, containing two overlaid silicon microstrip sensors, aligned with a parallel orientation. On-chip logic identifies Level-1 trigger primitives from high transverse-momentum tracks by selecting correlated clusters in the two sensors. The CBC3 was delivered in late 2016; wafer probing and performance tests have been carried out. Several prototype modules using the CBC3 have been produced and tested in the lab and in different beams. The results show that the CBC3 satisfies CMS requirements and only small corrections are needed for the final version of the chip for production.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2015

Mechanistic Details and Reactivity Descriptors in Oxidation and Acid Catalysis of Methanol

Acid and redox reaction rates of CH₃OH-O₂ mixtures on polyoxometalate (POM) clusters, together with isotopic, spectroscopic, and theoretical assessments of catalyst properties and reaction pathways, were used to define rigorous descriptors of reactivity and to probe the compositional effects for oxidative dehydrogenation (ODH) and dehydration reactions. ³¹P-MAS NMR, transmission electron microscopy and titrations of protons with di-tert-butylpyridine during catalysis showed that POM clusters retained their Keggin structure upon dispersion on SiO₂ and after use in CH₃OH reactions. The effects of CH₃OH and O₂ pressures and of D-substitution on ODH rates show that C-H activation in molecularly adsorbed CH₃OH is the sole kinetically relevant step and leads to reduced centers as intermediates present at low coverages; their concentrations, measured from UV-vis spectra obtained during catalysis, are consistent with the effects of CH₃OH/O₂ ratios predicted from the elementary steps proposed. First-order ODH rate constants depend strongly on the addenda atoms (Mo vs W) but weakly on the central atom (P vs Si) in POM clusters, because C-H activation steps inject electrons into the lowest unoccupied molecular orbitals (LUMO) of the clusters, which are the d-orbitals at Mo⁶⁺ and W⁶⁺ centers. H-atom addition energies (HAE) at O-atoms in POM clusters represent the relevant theoretical probe of the LUMO energies and of ODH reactivity. The calculated energies of ODH transition states at each O-atom depend linearly on their HAE values with slopes near unity, as predicted for late transition states in which electron transfer and C-H cleavage are essentially complete. HAE values averaged over all accessible O-atoms in POM clusters provide the appropriate reactivity descriptor for oxides whose known structures allow accurate HAE calculations. CH₃OH dehydration proceeds via parallel pathways mediated by late carbenium-ion transition states; effects of composition on dehydration reactivity reflect changes in charge reorganizations and electrostatic forces that stabilize protons at Brønsted acid sites.

Environmental Molecular Sciences Laboratory↗
DOE OSTI2011

The microstructure of dislocation clusters in industrial directionally solidified multicrystalline silicon

The microstructure of commonly occurring dislocation patterns in industrial directionally solidified multicrystalline silicon has been systematically studied by light microscopy, electron backscatter diffraction, and transmission electron microscopy. The work has been focused on dislocation clusters on wafers near the top of cast blocks. In near {l_brace}111{r_brace} grain surface, dislocation arrays parallel to {l_brace}110{r_brace} plane traces are lying in parallel rows of {l_brace}111{r_brace} planes inclined to the surface, in mainly <112>30 deg. orientation. The dislocation configuration suggests that the microstructure may result from a recovery process. The dislocations formed during crystal growth and cooling have undergone transformations at high temperature in order to achieve low energy configurations for minimization of dislocation and crystal energy.

36 MATERIALS SCIENCE↗
DOE OSTI2012

SNP Discovery with EST and NextGen Sequencing in Switchgrass ( Panicum virgatum L. )

Although yield trials for switchgrass (Panicum virgatum L.), a potentially high value biofuel feedstock crop, are currently underway throughout North America, the genetic tools for crop improvement in this species are still in the early stages of development. Identification of high-density molecular markers, such as single nucleotide polymorphisms (SNPs), that are amenable to high-throughput genotyping approaches, is the first step in a quantitative genetics study of this model biofuel crop species. We generated and sequenced expressed sequence tag (EST) libraries from thirteen diverse switchgrass cultivars representing both upland and lowland ecotypes, as well as tetraploid and octoploid genomes. We followed this with reduced genomic library preparation and massively parallel sequencing of the same samples using the Illumina Genome Analyzer technology platform. EST libraries were used to generate unigene clusters and establish a gene-space reference sequence, thus providing a framework for assembly of the short sequence reads. SNPs were identified utilizing these scaffolds. We used a custom software program for alignment and SNP detection and identified over 149,000 SNPs across the 13 short-read sequencing libraries (SRSLs). Approximately 25,000 additional SNPs were identified from the entire EST collection available for the species. This sequencing effort generated data that are suitable for marker development and for estimation of population genetic parameters, such as nucleotide diversity and linkage disequilibrium. Based on these data, we assessed the feasibility of genome wide association mapping and genomic selection applications in switchgrass. Overall, the SNP markers discovered in this study will help facilitate quantitative genetics experiments and greatly enhance breeding efforts that target improvement of key biofuel traits and development of new switchgrass cultivars.

09 BIOMASS FUELS↗
NASA NTRS1981

Prime agricultural land monitoring and assessment component of the California Integrated Remote Sensing System

The use of digital LANDSAT techniques for monitoring agricultural land use conversions was studied. Two study areas were investigated: one in Ventura County and the other in Fresno County (California). Ventura test site investigations included the use of three dates of LANDSAT data to improve classification performance beyond that previously obtained using single data techniques. The 9% improvement is considered highly significant. Also developed and demonstrated using Ventura County data is an automated cluster labeling procedure, considered a useful example of vertical data integration. Fresno County results for a single data LANDSAT classification paralleled those found in Ventura, demonstrating that the urban/rural fringe zone of most interest is a difficult environment to classify using LANDSAT data. A general raster to vector conversion program was developed to allow LANDSAT classification products to be transferred to an operational county level geographic information system in Fresno.

Estes, J. E.↗
NASA NTRS2006

Bragg Curve, Biological Bragg Curve and Biological Issues in Space Radiation Protection with Shielding

The space environment consists of a varying field of radiation particles including high-energy ions, with spacecraft shielding material providing the major protection to astronauts from harmful exposure. Unlike low-LET gamma or X-rays, the presence of shielding does not always reduce the radiation risks for energetic charged particle exposure. Since the dose delivered by the charged particle increases sharply as the particle approaches the end of its range, a position known as the Bragg peak, the Bragg curve does not necessarily represent the biological damage along the particle traversal since biological effects are influenced by the track structure of both primary and secondary particles. Therefore, the biological Bragg curve is dependent on the energy and the type of the primary particle, and may vary for different biological endpoints. To achieve a Bragg curve distribution, we exposed cells to energetic heavy ions with the beam geometry parallel to a monolayer of fibroblasts. Qualitative analyses of gamma-H2AX fluorescence, a known marker of DSBs, indicated increased clustering of DNA damage before the Bragg peak, enhanced homogenous distribution at the peak, and provided visual evidence of high linear energy transfer (LET) particle traversal of cells beyond the Bragg peak. A quantitative biological response curve generated for micronuclei (MN) induction across the Bragg curve did not reveal an increased yield of MN at the location of the Bragg peak. However, the ratio of mono-to bi-nucleated cells, which indicates inhibition in cell progression, increased at the Bragg peak location. These results, along with other biological concerns, show that space radiation protection with shielding can be a complicated issue.

Honglu, Wu↗
NASA NTRS2019

Leveraging STARE for Co-aligned Data Locality with netCDF and Python MPI

We have leveraged STARE indexing to package partitioned data chunks from diverse datasets into netCDF files, distributed them on a cluster of 16 lightweight nodes with their placements spatiotemporally co-aligned, and demonstrated a few integrative analyses using netCDF parallel I/O and Python MPI, with single-user performance and scalability comparable to, or even better than, that of a parallel array database management system (ADBMS) such as SciDB. However, records of the node location and STARE index ranges for each data chunk, similar to the chunk maps of SciDB, must be maintained and consulted by the I/O and analysis code for coordinating the analytic operations in parallel, in order to achieve the good performance and scalability.

Kwo-Sen Kuo↗
DOE OSTI2020

Scalable All-pairs Shortest Paths for Huge Graphs on Multi-GPU Clusters

We present an optimized Floyd-Warshall (Floyd-Warshall) algorithm that computes the All-pairs shortest path (APSP) for GPU accelerated clusters. The Floyd-Warshall algorithm due to its structural similarities to matrix-multiplication is well suited for highly parallel GPU architectures. To achieve high parallel efficiency, we address two key algorithmic challenges: reducing high communication overhead and addressing limited GPU memory. To reduce high communication costs, we redesign the parallel (a) to expose more parallelism, (b) aggressively overlap communication and computation with pipelined and asynchronous scheduling of operations, and (c) tailored MPI-collective. To cope with limited GPU memory, we employ an offload model, where the data resides on the host and is transferred to GPU on-demand. The proposed optimizations are supported with detailed performance models for tuning. Our optimized parallel Floyd-Warshall implementation is up to 5x faster than a strong baseline and achieves 8.1 PetaFLOPS/sec on 256~nodes of the Summit supercomputer at Oak Ridge National Laboratory. This performance represents 70% of the theoretical peak and 80% parallel efficiency. The offload algorithm can handle 2.5x larger graphs with a 20% increase in overall running time.

Sao, Piyush↗
DOE OSTI2020

Unidirectional growth of graphene nano-islands from carbon cluster seeds on Ge(1 1 0)

The anisotropic twofold symmetry of Ge(1 1 0) makes it a unique substrate for the growth of single-crystalline graphene. However, the underlying mechanism during the initial stage of growth of graphene on Ge(1 1 0) surface is not well understood. Here, we use in-situ cryogenic scanning tunneling microscopy (STM) and spectroscopy (STS), to study the initial growth properties of graphene synthesized on Ge(1 1 0) surface with ethylene precursor gases. The STM results reveal that the unidirectional growth of the graphene nanoribbons (GNRs) initiates from carbon cluster seeds on terraces during the early stage of the growth. The orientations of the GNRs show the same directions which are nearly parallel to the direction of the Ge(1 1 0) surface. This orientation is the same as that of the wafer-scale graphene grown on Ge(1 1 0). Subsequent growth, transformed the GNRs to graphene nanoislands (GNIs), which eventually coalesced to form single crystalline monolayer graphene. STS measurements demonstrated that the GNRs have small bandgaps induced by the confinement effect of graphene. Furthermore, this study provides an in-depth understanding of the growth mechanism of graphene on Ge(1 1 0) surface in various synthesis conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Structure sensitivity in the nonscalable regime explored via catalysed ethylene hydrogenation on supported platinum nanoclusters

The sensitivity, or insensitivity, of catalysed reactions to catalyst structure is a commonly employed fundamental concept. Here we report on the nature of nano-catalysed ethylene hydrogenation, investigated through experiments on size-selected Pt n (n=8-15) clusters soft-landed on magnesia and first-principles simulations, yielding benchmark information about the validity of structure sensitivity/insensitivity at the bottom of the catalyst size range. Both ethylene-hydrogenation-to-ethane and the parallel hydrogenation–dehydrogenation ethylidyne-producing route are considered, uncovering that at the <1 nm size-scale the reaction exhibits characteristics consistent with structure sensitivity, in contrast to structure insensitivity found for larger particles. The onset of catalysed hydrogenation occurs for Pt n (n≥10) clusters at T>150 K, with maximum room temperature reactivity observed for Pt 13 . Structure insensitivity, inherent for specific cluster sizes, is induced in the more active Pt 13 by a temperature increase up to 400 K leading to ethylidyne formation. As a result, control of sub-nanometre particle size may be used for tuning catalysed hydrogenation activity and selectivity.

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

Assessing the concept of structure sensitivity or insensitivity for sub-nanometer catalyst materials

In this report the nature of the nano-catalyzed hydrogenation of ethylene, yielding benchmark information pertaining to the concept of structure sensitivity/insensitivity and its applicability at the bottom of the catalyst particle size-range, is explored with experiments on size-selected Ptn (n = 7–40) clusters soft-landed on MgO, in conjunction with first-principles simulations. As in the case of larger particles both the direct ethylene hydrogenation channel and the parallel hydrogenation–dehydrogenation ethylidyne-producing route must be considered, with the fundamental uncovering that at the < 1 nm size-scale the reaction exhibits characteristics consistent with structure sensitivity, in contrast to the structure insensitivity found for larger particles. In this size-regime, the chemical properties can be modulated and tuned by a single atom, reflected by the onset of low temperature hydrogenation at T > 150 K catalyzed by Pt n (n ≥ 10) clusters, with maximum room temperature reactivity observed for Pt 13 using a pulsed molecular beam technique. Structure insensitive behavior, inherent for specific cluster sizes at ambient temperatures, can be induced in the more active sizes, e.g. Pt 13 , by a temperature increase, up to 400 K, which opens dehydrogenation channels leading to ethylidyne formation. This reaction channel was, however found to be attenuated on Pt 20 , as catalyst activity remained elevated after the 400 K step. Pt 30 displayed behavior which can be understood from extrapolating bulk properties to this size range; in particular the calculated d-band center. In the non-scalable sub-nanometer size regime, however, precise control of particle size may be used for atom-by-atom tuning and manipulation of catalyzed hydrogenation activity and selectivity.

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

Symmetry and Topology Code of the Cluster Self-Assembly of Intermetallic Compounds A{sub 2}{sup [16]}B{sub 4}{sup [12]} of the Friauf Families Mg{sub 2}Cu{sub 4} and Mg{sub 2}Zn{sub 4}

A geometrical and topological analysis of intermetallic compounds A{sub 2}{sup [16]}B{sub 4}{sup [12]} of the Friauf families Mg{sub 2}Cu{sub 4} (formed in 273 binary and 775 ternary systems) and Mg{sub 2}Zn{sub 4} (formed in 67 binary and 121 ternary systems) has been performed. Basic 3D nets are determined in the form of a graph, whose sites correspond to the centers of tetrahedral cluster precursors В{sub 4}: Cu-type net with a coordination number (CN) of 12 for Mg{sub 2}Cu{sub 4} and a BN-type net with a CN of 8 for Mg{sub 2}Zn{sub 4}. The symmetry and topology code is determined for the crystal structure self-assembly in the following form: primary chain S{sub 3}{sup 1} → microlayer S{sub 3}{sup 2} → microframework S{sub 3}{sup 3}. An identical mechanism of primary chain S{sub 3}{sup 1} self-assembly from B{sub 4} clusters with participation of large template atoms A{sup [16]} for both families and a mechanism of S{sub 3}{sup 2} microlayer self-assembly from parallel S{sub 3}{sup 1} chains in Mg{sub 2}Cu{sub 4} (basic 2D net 4{sup 4}) and shifted chains in Mg{sub 2}Zn{sub 4} (basic 2D net 3{sup 6}) have been determined for the first time. The constant chemical composition of intermetallic compounds A{sub 2}{sup [16]}B{sub 4}{sup [12]} is related to the functional difference between the B atoms forming 3D frameworks and template atoms А. The numbers of chemically different В{sub 4} clusters for the Mg{sub 2}Cu{sub 4} and Mg{sub 2}Zn{sub 4} families are found to be 25 and 19, respectively.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Transitioning NWChem to the Next Generation of Manycore Machines

The NorthWest Chemistry (NWChem) modeling software is a popular molecular chemistry simulation software that was designed from the start to work on massively parallel processing supercomputers[6, 28, 49]. It contains an umbrella of modules that today includes Self Consistent Field (SCF), second order Mller-Plesset perturbation theory (MP2), Coupled Cluster, multi-conguration selfconsistent eld (MCSCF), selected conguration interaction (CI), tensor contraction engine (TCE) many body methods, density functional theory (DFT), time-dependent density functional theory (TDDFT), real time time-dependent density functional theory, pseudopotential plane-wave density functional theory (PSPW), band structure (BAND), ab initio molecular dynamics, Car-Parrinello molecular dynamics, classical molecular dynamics (MD), QM/MM, AIMD/MM, GIAO NMR, COSMO, COSMO-SMD, and RISM solvation models, free energy simulations, reaction path optimization, parallel in time, among other capabilities[ 22]. Moreover new capabilities continue to be added with each new release.

Environmental Molecular Sciences Laboratory↗
DOE OSTI2017

Divide and Conquer (DC) BLAST: fast and easy BLAST execution within HPC environments

Bioinformatics is currently faced with very large-scale data sets that lead to computational jobs, especially sequence similarity searches, that can take absurdly long times to run. For example, the National Center for Biotechnology Information (NCBI) Basic Local Alignment Search Tool (BLAST and BLAST+) suite, which is by far the most widely used tool for rapid similarity searching among nucleic acid or amino acid sequences, is highly central processing unit (CPU) intensive. While the BLAST suite of programs perform searches very rapidly, they have the potential to be accelerated. In recent years, distributed computing environments have become more widely accessible and used due to the increasing availability of high-performance computing (HPC) systems. Therefore, simple solutions for data parallelization are needed to expedite BLAST and other sequence analysis tools. However, existing software for parallel sequence similarity searches often requires extensive computational experience and skill on the part of the user. In order to accelerate BLAST and other sequence analysis tools, Divide and Conquer BLAST (DCBLAST) was developed to perform NCBI BLAST searches within a cluster, grid, or HPC environment by using a query sequence distribution approach. Scaling from one (1) to 256 CPU cores resulted in significant improvements in processing speed. Thus, DCBLAST dramatically accelerates the execution of BLAST searches using a simple, accessible, robust, and parallel approach. DCBLAST works across multiple nodes automatically and it overcomes the speed limitation of single-node BLAST programs. DCBLAST can be used on any HPC system, can take advantage of hundreds of nodes, and has no output limitations. Thus, this freely available tool simplifies distributed computation pipelines to facilitate the rapid discovery of sequence similarities between very large data sets.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2017

Large-scale atomistic calculations of clusters in intense x-ray pulses

Here, we present the methodology of our recently developed Monte-Carlo/ Molecular-Dynamics method for studying the fundamental ultrafast dynamics induced by high-fluence, high-intensity x-ray free electron laser (XFEL) pulses in clusters. The quantum nature of the initiating ionization process is accounted for by a Monte Carlo method to calculate probabilities of electronic transitions, including photo absorption, inner-shell relaxation, photon scattering, electron collision and recombination dynamics, and thus track the transient electronic configurations explicitly. The freed electrons and ions are followed by classical particle trajectories using a molecular dynamics algorithm. These calculations reveal the surprising role of electron-ion recombination processes that lead to the development of nonuniform spatial charge density profiles in x-ray excited clusters over femtosecond timescales. In the high-intensity limit, it is important to include the recombination dynamics in the calculated scattering response even for a 2- fs pulse. We also demonstrate that our numerical codes and algorithms can make e!cient use of the computational power of massively parallel supercomputers to investigate the intense-field dynamics in systems with increasing complexity and size at the ultrafast timescale and in non-linear x-ray interaction regimes. In particular, picosecond trajectories of XFEL clusters with attosecond time resolution containing millions of particles can be e!ciently computed on upwards of 262,144 processes.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

Parallel performance of algebraic multigrid domain decomposition

Algebraic multigrid (AMG) is a widely used scalable solver and preconditioner for large-scale linear systems resulting from the discretization of a wide class of elliptic PDEs. While AMG has optimal computational complexity, the cost of communication has become a significant bottleneck that limits its scalability as processor counts continue to grow on modern machines. This article examines the design, implementation, and parallel performance of a novel algorithm, algebraic multigrid domain decomposition (AMG-DD), designed specifically to limit communication. The goal of AMG-DD is to provide a low-communication alternative to standard AMG V-cycles by trading some additional computational overhead for a significant reduction in communication cost. Numerical results show that AMG-DD achieves superior accuracy per communication cost compared with AMG, and speedup over AMG is demonstrated on a large GPU cluster.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2017

Parallel Fokker–Planck-DSMC algorithm for rarefied gas flow simulation in complex domains at all Knudsen numbers

A major challenge for the conventional Direct Simulation Monte Carlo (DSMC) technique lies in the fact that its computational cost becomes prohibitive in the near continuum regime, where the Knudsen number (Kn)—characterizing the degree of rarefaction—becomes small. In contrast, the Fokker–Planck (FP) based particle Monte Carlo scheme allows for computationally efficient simulations of rarefied gas flows in the low and intermediate Kn regime. The Fokker–Planck collision operator—instead of performing binary collisions employed by the DSMC method—integrates continuous stochastic processes for the phase space evolution in time. This allows for time step and grid cell sizes larger than the respective collisional scales required by DSMC. Dynamically switching between the FP and the DSMC collision operators in each computational cell is the basis of the combined FP-DSMC method, which has been proven successful in simulating flows covering the whole Kn range. Until recently, this algorithm had only been applied to two-dimensional test cases. In this contribution, we present the first general purpose implementation of the combined FP-DSMC method. Utilizing both shared- and distributed-memory parallelization, this implementation provides the capability for simulations involving many particles and complex geometries by exploiting state of the art computer cluster technologies.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Local Distortions and Metal–Semiconductor–Metal Transition in Quasi-One-Dimensional Nanowire Compounds AV 3 Q 3 O δ (A = K, Rb, Cs and Q = Se, Te)

Metal cluster compounds have garnered renewed interest in the search for novel superconductors and topological semimetals owing to structural instabilities of metal-cluster geometries and broken symmetries. Here we synthesized needle-like crystals of the V-cluster-based quasi-one-dimensional (Q1D) materials AV 3 Q 3 O δ (A = K, Rb, Cs, Q= Se, Te) which can also be viewed as being composed of parallel nanowires. We examine how changes in their average and local structure control their electronic properties. All compounds crystallize in the TIFe 3 Te 3 -type structure (P6 3 /m space group) with infinite (V 3 )Q 3 ) - double-walled columnar chains separated by A + cations. Our single-crystal and synchrotron powder diffraction studies indicate oxygen atoms partially occupy the center site of the V 6 octahedral metal cluster cages in KV 3 Te 3 O 0.33 , RbV 3 Te 3 O 0.32 , and CsV 3 Te 3 O 0.35 , whereas KV 3 Se 3 is structurally oxygen-free. Our synchrotron X-ray pair distribution function (PDF) analyses indicate that the oxygen-free V 6 cluster octahedra in KV 3 Se 3 are highly distorted perpendicular to the chain direction even at room temperature, reducing the symmetry of the average structure from hexagonal P6 3 /m to monoclinic P2 1 /m. Our theoretical calculation supports this P2 1 /m distortion and suggests the structure further distorts to P2 1 or P2 1 /c at lower temperatures. In contrast, the oxygen-centered V-cluster in KV 3 Te 3 O 0.33 exhibits a V-3-triangle-trimerization along the chain direction. This feature is discernible from the local PDF and is consistent with lattice dynamical calculations based on density functional theory. Resistivity measurements indicate that KV 3 Se 3 exhibits metallic behavior, whereas a dramatic metal-semiconductor-metal transition emerges in KV 3 Te 3 O 6.33 , RbV 3 Te 3 O 0.32 , and CsV 3 Te 3 O 0.35 because of oxygen disorder and changes in local structure captured from our electronic structure analyses of the Fermi surface. Our investigation of the AV 3 Q 3 O 6 family demonstrates the importance of understanding local changes in structure driven by electronic instabilities, which can guide the search for new quantum materials in other low-dimensional cluster-compound materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗