Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “STATISTICAL MECHANICS”

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 OSTI2023

Uncertainty quantification for high explosive reactant and product equations of state

Equations of state (EOSs) are typically represented as physics-informed models with tunable parameters that are adjusted to replicate calibration data as closely as possible. Uncertainty quantification (UQ) allows for the development of an ensemble of EOS parameters that are consistent with the calibration data instead of a single EOS. In this work, we perform UQ for the reactant and product EOSs for a variety of high explosives (HEs). In doing so, we demonstrate a strategy for dealing with heterogeneous (both experimental and calculated) data. We also use a statistical distance metric to quantify the differences between the various HEs using the UQ results.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Mechanical Performance and Mesostructure Analysis of Proton-Irradiated Fused Filament Fabrication Acrylonitrile Butadiene Styrene Material

The use of fused filament fabrication (FFF) acrylonitrile butadiene styrene (ABS) and other thermoplastics in radiation environments is beginning to be studied as possible replacements for traditionally manufactured parts and tools. Interlayer adhesion within the mesostructure has been shown in published literature to be an integral component in the strength of FFF ABS. Research of irradiated 3D printed polymers has primarily focused on the influence of gamma irradiation on the mechanical properties of FFF ABS samples, without evaluating its impact on the mesostructure of the samples. The purpose of this paper is to understand the mechanical damage caused by proton radiation on FFF ABS samples through the evaluation of the mesostructure of the ABS samples. To achieve this objective, proton radiation at 40 MeV was applied to FFF ABS samples at radiation doses up to 1.0 MGy at a dose rate of 1 MGy/hr. Following the irradiation, tensile testing was performed on the samples. Here, the fractured surfaces of the tested samples were subsequently observed using a scanning electron microscope. The mechanical test results show an increase in the ultimate tensile strength (UTS) and a decrease in the ductility of the irradiated samples. Statistical analysis on the results shows that there is a statistically significant difference in the UTS and the ductility of unirradiated and 1.0 MGy irradiated samples and between irradiated samples and 1.0 MGy samples. In addition, the percentage elongation at break has statistically significant differences in the means between irradiated samples and 1.0 MGy samples. The UTS has statistically significant differences in the means between 0.1 MGy and 1.0 MGy. The difference in the means between the unirradiated and the 1.0 MGy samples is the most significant for both the UTS and elongation at break. The scanning electron microscopy (SEM) results indicated that interlayer adhesion improved as a function of radiation dose corresponding to the increase in tensile strength. The SEM results also showed that crazing and plastic deformation were reduced; aligning with the loss in ductility observed in the tensile tests results. The proton radiation is causing these mechanical and physical changes through two mechanisms: (1) dose rate effects on ionizing radiation-induced oxidative degradation; and (2) radiation heating effects of high energy (>1 MeV) charged particles.

electron microscopy↗
DOE OSTI2025

Examining the Impact of Local Constraint Violations on Energy Computations in DFT

ABSTRACT This work examines the impact of locally imposed constraints in Density Functional Theory (DFT). Using a metric referred to as the extent of violation index (EVI), we quantify how well exchange‐correlation functionals adhere to local constraints. Applying EVIs to a diverse set of molecules for GGA functionals reveals constraint violations, particularly for semi‐empirical functionals. We leverage EVIs to explore potential connections between these violations and errors in chemical properties. While no correlation is observed for atomization energies, a significant statistical correlation emerges between EVIs and total energies. Similarly, the analysis of reaction energies suggests weak positive correlations for specific constraints. However, definitive conclusions about error cancellation mechanisms cannot be made at this time. These observations revealed by EVIs may be useful for consideration when designing future generations of semilocal functionals.

Khanna, Vaibhav [Department of Chemistry Universit↗
DOE OSTI2024

The metamorphosis of semi-classical mechanisms of confinement: from monopoles on ℝ 3 × S 1 to center-vortices on ℝ 2 × T 2

There are two distinct regimes of Yang-Mills theory where we can demonstrate confinement, the existence of a mass gap, and the multi-branch structure of the effective potential as a function of the theta angle using a reliable semi-classical calculation. The two regimes are deformed Yang-Mills theory on ℝ 3 × S 1 , and Yang-Mills theory on ℝ 2 × T 2 where the torus is threaded by a ’t Hooft flux. The weak coupling regime is ensured by the small size of the circle or torus. In the first case the confinement mechanism is related to self-dual monopoles, whereas in the second case self-dual center-vortices play a crucial role. These two topological objects are distinct. In particular, they have different mutual statistics with Wilson loops. On the other hand, they carry the same topological charge and action. We consider the theory on ℝ × T 2 × S 1 and extrapolate both the monopole and vortex regimes to a quantum mechanical domain, where a cross-over takes place. Both sides of the cross-over are described by a deformed ℤ N TQFT. On ℝ 2 × S 1 × S 1 , we derive an effective field theory (EFT) of vortices from the EFT of monopoles in the presence of a ’t Hooft flux. This construction is based on a two-stage Higgs mechanism, reducing SU(N) to U(1) N−1 in 3d first, followed by reduction to a ℤ N EFT in 2d in the second step. This result shows how monopoles transmute into center-vortices, and suggests adiabatic continuity between the two confinement mechanisms. The basic mechanism is flux fractionalization: the magnetic flux of the monopoles splits up and is collimated in such a way that 2d Wilson loops detect it as a center vortex.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Towards understanding the formation of internal fragments generated by collisionally activated dissociation for top-down mass spectrometry

Top-down mass spectrometry (TD-MS) generates fragment ions that returns information on the polypeptide amino acid sequence. In addition to terminal fragments, internal fragments that result from multiple cleavage events can also be formed. Traditionally, internal fragments are largely ignored due to a lack of available software to reliably assign them, mainly caused by a poor understanding of their formation mechanism. To accurately assign internal fragments, their formation process needs to be better understood. Here we applied a statistical method to compare fragmentation patterns of internal and terminal fragments of peptides and proteins generated by collisionally activated dissociation (CAD). Internal fragments share similar fragmentation propensities with terminal fragments (e.g., enhanced cleavages N-terminal to proline and C-terminal to acidic residues), suggesting that their formation follows conventional CAD pathways. Internal fragments should be generated by subsequent cleavages of terminal fragments and their formation can be explained by the well-known mobile proton model. Additionally, internal fragments can be coupled with terminal fragments to form complementary product ions that span the entire protein sequence. These enhance our understanding of internal fragment formation and can help improve sequencing algorithms to accurately assign internal fragments, which will ultimately lead to more efficient and comprehensive TD-MS analysis of proteins and proteoforms.

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

Nonlinear coupling of electromagnetic and electron acoustic waves in multi-species degenerate astrophysical plasma

Nonlinear wave-coupling is studied in a multispecies degenerate astrophysical plasma consisting of two electron species (at different temperatures): a highly degenerate main component plus a smaller classical relativistic flow immersed in a static neutralizing ion background. It is shown that the high frequency electromagnetic waves through their strong nonlinear interactions with the electron-acoustic waves [sustained by a multielectron component (degenerate) plasma surrounding a compact astrophysical object] can scatter to lower frequencies so that the radiation observed faraway will be spectrally shifted downward. It is also shown that, under definite conditions, the electromagnetic waves could settle into stationary solitonic states. It is expected that the effects of such structures may persist as detectable signatures in forms of modulated micropulses in the radiation observed far away from the accreting compact object. Both these effects will advance our abilities to interpret the radiation coming out of the compact objects.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

The quantum p -spin glass model: a user manual for holographers

We study a large-N bosonic quantum mechanical sigma-model with a spherical target space subject to disordered interactions, more colloquially known as the p-spin spherical model. Replica symmetry is broken at low temperatures and for sufficiently weak quantum fluctuations, which drives the system into a spin glass (SG) phase. The first half of this paper is dedicated to a discussion of this model's thermodynamics, with particular emphasis on the marginally stable SG. This phase exhibits an emergent conformal symmetry in the strong coupling regime, which dictates its thermodynamic properties. It is associated with an extensive number of nearby states in the free energy landscape. We discuss in detail an elegant approximate solution to the SG equations, which interpolates between the conformal regime and an ultraviolet-complete short distance solution. In the second half of this paper we explore the real-time dynamics of the model and uncover quantum chaos as measured by out-of-time-order four-point functions, both numerically and analytically. We find exponential Lyapunov growth, which intricately depends on the model's couplings and becomes strongest in the quantum critical regime. We emphasize that the SG phase also exhibits quantum chaos, albeit with parametrically smaller Lyapunov exponent than in the replica symmetric phase. An analytical calculation in the marginal SG phase suggests that this Lyapunov exponent vanishes in a particular infinite coupling limit. Finally, we comment on the potential meaning of these observations from the perspective of holography.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Gauge-symmetrization method for energy-momentum tensors in high-order electromagnetic field theories

For electromagnetic field theories, canonical energy-momentum conservation laws can be derived from the underpinning spacetime translation symmetry according to the Noether procedure. However, the canonical energy-momentum tensors (EMTs) are neither symmetric nor gauge-symmetric (gauge invariant). The Belinfante-Rosenfeld (BR) method is a well-known procedure to symmetrize the EMTs, which also renders them gauge symmetric for first-order field theories. High-order electromagnetic field theories appear in the study of gyrokinetic systems for magnetized plasmas and the Podolsky system for the radiation reaction of classical charged particles. For these high-order field theories, gauge-symmetric EMTs are not necessarily symmetric and vice versa. In the present study, we develop a new gauge-symmetrization method for EMTs in high-order electromagnetic field theories. The Noether procedure is carried out using the Faraday tensor $F_{μν}$, instead of the 4-potential $A_{μ}$, to derive a canonical EMT $T^{μν}_{N}$. We show that the gauge-dependent part of $T^{μν}_{N}$ can be removed using the displacement-potential tensor $F^{σμν}$ ≡ $D^{σμ}A^{ν}/4π$, where $D^{σμ}$ is the antisymmetric electric displacement tensor. This method gauge-symmetrizes the EMT without necessarily making it symmetric, which is adequate for applications not involving general relativity. For first-order electromagnetic field theories, such as the standard Maxwell system, $F^{σμν}$ reduces to the familiar BR superpotential $S^{σμν}$, and the method developed can be used as a simpler procedure to calculate $S^{σμν}$ without employing the angular momentum tensor in 4D spacetime. When the electromagnetic system is coupled to classical charged particles, the gauge-symmetrization method for EMTs is shown to be effective as well.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Lieb-Robinson Light Cone for Power-Law Interactions

The Lieb-Robinson theorem states that information propagates with a finite velocity in quantum systems on a lattice with nearest-neighbor interactions. What are the speed limits on information propagation in quantum systems with power-law interactions, which decay as 1 / r α at distance r ? Here, we present a definitive answer to this question for all exponents α > 2 d and all spatial dimensions d . Schematically, information takes time at least r min { 1 , α - 2 d } to propagate a distance r . As recent state transfer protocols saturate this bound, our work closes a decades-long hunt for optimal Lieb-Robinson bounds on quantum information dynamics with power-law interactions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Sterilization of N95 Respirators via Gamma Radiation: Comparison of Post-sterilization Efficacy

This study evaluated gamma irradiation for sterilization and reuse of two models of N95 respirators after gamma radiation sterilization as a method to increase availability of N95 respirators during a shortage. The Sandia National Laboratories Gamma Irradiation Facility was used to irradiate two different models of N95 filtering facepiece respirators at doses ranging from 0 kGy(tissue) to 50 kGy(tissue). The following tests were used to determine the efficacy of the respirator after irradiation sterilization: Ambient Aerosol Condensation Nuclei Counter Quantitative Fit Test, tensile test, strain cycling, oscillatory dynamic mechanical analysis, microscopic image analysis of fiber layers, and electrostatic field measurements. Both of the respirator models exhibited statistically significant changes after gamma irradiation as shown by the Quantitative Fit Test, electrostatic testing and the aerosol testing. The change in electrostatic capability of the filter reduced the efficiency of challenging particles near the 200 nm size by approximately 40-50%. Both tested respirators showed statistically significant changes associated with gamma sterilization. However, our results indicate that choices in materials and manufacturing methods to achieve N95 filtration lead to different magnitudes of damage when exposed to gamma radiation at sterilization relevant doses. This damage results in lower filtration performance. While our sample size (2 different types of respirators) was small, we did observe a change in electrostatic properties on a filter layer that coincided with the failure on the Quantitative Fit Test.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2020

Sterilization of N95 Respirators via Gamma Radiation: Comparison of Post-sterilization Efficacy

This study evaluated gamma irradiation for sterilization and reuse of two models of N95 respirators after gamma radiation sterilization as a method to increase availability of N95 respirators during a shortage. The Sandia National Laboratories Gamma Irradiation Facility was used to irradiate two different models of N95 filtering facepiece respirators at doses ranging from 0 kGy(tissue) to 50 kGy(tissue). The following tests were used to determine the efficacy of the respirator after irradiation sterilization: Ambient Aerosol Condensation Nuclei Counter Quantitative Fit Test, tensile test, strain cycling, oscillatory dynamic mechanical analysis, microscopic image analysis of fiber layers, and electrostatic field measurements. Both of the respirator models exhibited statistically significant changes after gamma irradiation as shown by the Quantitative Fit Test, electrostatic testing and the aerosol testing. The change in electrostatic charge of the filter was correlated with a reduction in capturing particles near the 200 nm size by approximately 40-50%. Both tested respirators showed statistically significant changes associated with gamma sterilization. However, our results indicate that choices in materials and manufacturing methods to achieve N95 filtration lead to different magnitudes of damage when exposed to gamma radiation at sterilization relevant doses. This damage results in lower filtration performance. While our sample size (2 different types of respirators) was small, we did observe a change in electrostatic properties on a filter layer that coincided with the failure on the Quantitative Fit Test and reduction in aerosol filtering efficiency. Key Words: N95 respirators, respirators, airborne transmission, pandemic prevention, COVID-19, gamma sterilization

COVID-19↗
DOE OSTI2022

Does a Free Electron Laser Exhibit Non-Standard Statistics? [Slides and Report]

The originally proposed work under this award sought to develop a theoretical basis for previous experimental observations in free-electron laser (FEL) physics involving certain non-classical effects in the light produced by these devices. A sound and valid theoretical interpretation of these experiments would represent a transformational understanding of free-electron lasers, and of electron-photon interactions in general. The specific objective of the completed work, as described in the attached progress report by Jeongwan Park, was to perform a detailed analysis of an experiment by Chen and Madey (Physical Review Letters; volume 86 number 26, 2001) concerning the purported sub-Poissonian photon statistics in an FEL. Since that published work contradicted FEL theory as presently understood, the purpose of the work under the present subcontract was to obtain a better understanding of the quantum nature of the FEL, to investigate the possibility of theoretical and experimental evidence for the non-standard photon statistics of the FEL, and to critically examine the experimental evidence through a re-evaluation of the data analysis of Chen and Madey. The results of this study showed that there were numerous experimental conditions under which the observed photon statistics could be explained by combining both the photon clustering property and the dead-time effect, and consequently, that the observed photon statistics in the Chen-Madey experiment did not require sub-Poissonian photon statistics for their explanation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Isoscaling and nuclear reaction dynamics

Isoscaling parameters α and β have been explored as a function of breakup angle in binary excited projectile-like fragment decays produced in collisions of 70 Zn+ 70 Zn and 64 Zn+ 64 Zn at 35 MeV per nucleon. In this analysis, focus was placed on isoscaling the second heaviest fragment with 4 ≤ Z L ≤ 8 emitted from the excited projectile-like fragment in events that contained a heavy fragment with Z H ≥ 12. The breakup orientation θ prox was defined as the angle between the heavy and light fragments’ center-of-mass velocity and the fragment pairs’ relative velocity. Breakups between 0° < θ prox ≤ 80° have been shown to be dominated by dynamical contributions, while break-up angles of θ prox > 100° are predominantly statistical. Historically, isoscaling has often been understood and applied in a statistical context, assuming that the fragments are produced after statistical equilibrium is achieved. Finally, studying isoscaling parameters as a function of θ prox reveals the sensitivity of α and β to the mechanism of fragment production.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Lagrangian statistics of a shock-driven turbulent dynamo in decaying turbulence

Small-scale fluctuating magnetic fields of order nG are observed in supernova shocks and galaxy clusters, where its amplification is likely caused by the Biermann battery mechanism. However, these fields cannot be amplified further without the turbulent dynamo, which generates magnetic energy through the stretch-twist-fold (STF) mechanism. Thus, we present here novel 3D magnetohydrodynamic (MHD) simulations of a laser-driven shock propagating into a stratified, multiphase medium, to investigate the post-shock turbulent magnetic field amplification via the turbulent dynamo. The configuration used here is currently being tested in the shock tunnel at the National Ignition Facility (NIF). In order to probe the statistical properties of the post-shock turbulent region, we use 384 × 512 × 384 tracers to track its evolution through the Lagrangian framework, thus providing a high-fidelity analysis of the shocked medium. Our simulations indicate that the growth of the magnetic field, which accompanies the near-Saffman kinetic energy decay (E kin ∝ t -1.15 ) without turbulence driving, exhibits slightly different characteristics as compared to periodic box simulations. Seemingly no distinct phases exist in its evolution because the shock passage and time to observe the magnetic field amplification during the turbulence decay are very short (~0.3 of a turbulent turnover time). Yet, the growth rate is still consistent with those expected for compressive (curl-free) turbulence driving in subsonic, compressible turbulence. Phenomenological understanding of the dynamics of the magnetic and velocity fields are also elucidated via Lagrangian frequency spectra, which are consistent with the expected inertial range scalings in the Eulerian–Lagrangian bridge.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Attribution of heterogeneous stress distributions in low-grain polycrystals under conditions leading to damage

In high-purity polycrystalline metallic materials, voids tend to favor grain boundaries as nucleation sites due to the elevated stress states produced by granular interactions and the weakened grain boundary from the relative atomic disorder. To quantify the key factors of this elevated stress state, simple compression of a small multi-grain cylinder of body-centered cubic tantalum was simulated using a single crystal plasticity model that incorporates non-Schmid effects. Four increasingly complex synthetic microstructures were created to tractably incorporate grain boundary interactions, and a statistically significant number of combinations were performed by varying the initial crystallographic orientations of the microstructure. Most of these simulations produce the maximum von Mises stress on a grain boundary and less frequently at the multi-grain junctions. To build a statistical model for the maximum von Mises stress at the grain boundary, physically based features that could contribute to the elevated stress state were selected. Then, a learning algorithm based on information theory was used to identify which of these features contributed the most information to the data set. The identified features include a grain’s propensity to accommodate both elastic and plastic deformations and their directional components. The misalignment of the direction of each grain’s mechanical response was found to be strongly correlated to the magnitude of the stress near the grain boundary. For all of the synthetic microstructures, the statistical models produce a residual distribution that is nearly Gaussian with a variance of, at most, 10% of the prior distribution. The successful performance of the statistical model implies the correct identification of the physical features that cause severe stress localization in polycrystalline materials. The statistical models constructed here can be used to formulate a physically motivated void nucleation model which is sensitive to a microstructure’s propensity to produce elevated stress states. As a result, these statistical models also enable the design of material microstructures, in which the crystallographic orientation is chosen to resist void nucleation.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Quantification of microtubule stutters: dynamic instability behaviors that are strongly associated with catastrophe

Microtubules (MTs) are cytoskeletal fibers that undergo dynamic instability (DI), a remarkable process involving phases of growth and shortening separated by stochastic transitions called catastrophe and rescue. Dissecting DI mechanism(s) requires first characterizing and quantifying these dynamics, a subjective process that often ignores complexity in MT behavior. We present a Statistical Tool for Automated Dynamic Instability Analysis (STADIA) that identifies and quantifies not only growth and shortening, but also a category of intermediate behaviors that we term “stutters.” During stutters, the rate of MT length change tends to be smaller in magnitude than during typical growth or shortening phases. Quantifying stutters and other behaviors with STADIA demonstrates that stutters precede most catastrophes in our in vitro experiments and dimer-scale MT simulations, suggesting that stutters are mechanistically involved in catastrophes. Related to this idea, we show that the anticatastrophe factor CLASP2γ works by promoting the return of stuttering MTs to growth. STADIA enables more comprehensive and data-driven analysis of MT dynamics compared with previous methods. The treatment of stutters as distinct and quantifiable DI behaviors provides new opportunities for analyzing mechanisms of MT dynamics and their regulation by binding proteins.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Onion-Like nanoparticles of the metal–organic framework UiO-66 synthesized by sequential spike crystal growth

Reproducible size and morphology of metal–organic framework (MOF) particles are essential for tuning these materials for applications like heterogeneous catalysis. Particle size measurements from electron micrographs provide information about the impact of changing synthetic parameters on MOF growth, but reported data are usually limited to averages and statistical distributions. Further elucidating synthetic control parameters could facilitate studies of facet-dependent and non-classical growth mechanisms. Here, in this work, phase-pure nanoparticles of the MOF UiO-66 were synthesized using alternating spikes of zirconium and hafnium-based nodes, forming an onion-like structure with alternating layers of Zr-UiO-66 and Hf-UiO-66. In conventional, bright-field transmission electron micrographs, Hf-UiO-66 and Zr-UiO-66 layers appeared darker and lighter, respectively. Alternating layers of Hf-UiO-66 and Zr-UiO-66 were also apparent in elemental maps obtained using scanning transmission electron microscopy with energy-dispersive X-ray spectroscopy. Average particle sizes increased monotonically with each spike, and no Zr-only or Hf-only particles were observed after the first spike, indicating that new material was consistently incorporated onto existing seeds. Surface area of the onion-like particles and their uptake of Ni 2+ from solution were both between the values measured for Zr-UiO-66 and Hf-UiO-66. Results demonstrate the use of sequential spikes of nodes in monitoring crystal growth in MOFs.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Bottom-Up Simulation, Reconstruction, and Quantification of Macromolecule Sequences from Experimental Polymerizations

Motivated by the canonical sequence–structure–function paradigm, tools to characterize chemical patterning in natural biomacromolecules, from proteins to nucleic acids, have grown exponentially in recent years. However, analogous strategies for synthetic macromolecules remain in nascent stages, complicated by sequence polydispersity and analytical limitations. To address this, we have developed a comprehensive and open-source Python package, PRISM (polymer rate insights and sequence modeling), an end-to-end workflow that provides a path from experimental kinetics measurements to quantitative and qualitative metrics for describing chemical patterning in stochastic polymers. First, a numerical integration strategy was constructed to simulate and fit experimental data from reversible addition–fragmentation chain transfer (RAFT) polymerization kinetics, enabling the facile estimation of relevant reactivity ratios. These ratios were then used in a mechanism-specific stochastic kinetic simulation strategy to simulate sequence ensembles corresponding to model systems spanning experimental copolymers, classes of statistical polymers (e.g., alternating, block, and gradient), and multiblock copolymers. Lastly, inspired by sequence homology metrics from bioinformatics, we introduce visualization strategies and quantitative metrics to facilitate comparisons of different sequence ensembles. As the sequence–structure–function paradigm becomes increasingly central in de novo design of synthetic macromolecules, this toolkit provides a first step toward accurate and representative sequence description and featurization.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗