Search NASA⌕ Search

SEARCH · Search NASA

Results for “deformation function”

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 1,243 records · Page 69

Ultrafast Formation of Jahn–Teller Polarons Revealed by State-Selective Excitation in Correlated Spinel Co 3 O 4

Jahn–Teller polarons are quasiparticles that stem from symmetry breaking and strong local electron–phonon coupling. They originate from an excess charge carrier being dressed by a local lattice distortion, caused by the Jahn–Teller effect, and they critically impact electrical, structural, and magnetic properties in transition metal oxides. The observation of the microscopic steps involved in their formation is essential for enabling control over material properties through the targeted activation of local, site-specific modifications with light pulses. While Jahn–Teller distortion associated with polaron formation was predicted to contribute significantly to changes in electronic band gap and optical properties in Co 3 O 4 , its experimental observation remains elusive, requiring signatures of local symmetry reduction. In this work, we demonstrate Jahn–Teller polaron formation in spinel Co 3 O 4 . By exciting electronic transitions at 3.10 eV and 1.55 eV, we target either the Oh Cobalt(III) or the Td Cobalt(II) ions, and drive the subsequent coherent responses of the system through two different pathways. For the former, we demonstrate that ligand-to-metal charge transfer leads to Jahn–Teller polaron formation, which is linked to the deformation potential and magnetoelastic coupling. For the latter, we identify the coherent excitation of a T 2g phonon mode launched by on-site d-d electronic transitions. Key to our observations is the ability to target site-specific electronic excitations in spinel Co 3 O 4 using ultrafast optical pulses, while monitoring the ensuing low-energy collective modes through the coherent time-domain response of the material. Our approach, which combines the comparative analysis of experimental fingerprints with the support from density functional theory calculations, is broadly applicable to systems in which Jahn–Teller polaron physics has been theoretically predicted but remains experimentally unverified, and it underscores the potential of electronic-state targeting as a route to selectively excite and probe quasiparticle dynamics in solids.

Lattices↗

Linking electronic structure calculations to generalized stacking fault energies in multicomponent alloys

Abstract The generalized stacking fault energy is a key ingredient to mesoscale models of dislocations. Here we develop an approach to quantify the dependence of generalized stacking fault energies on the degree of chemical disorder in multicomponent alloys. We introduce the notion of a “configurationally-resolved planar fault” (CRPF) energy and extend the cluster expansion method from alloy theory to express the CRPF as a function of chemical occupation variables of sites surrounding the fault. We apply the approach to explore the composition and temperature dependence of the unstable stacking fault energy (USF) in binary Mo–Nb alloys. First-principles calculations are used to parameterize a formation energy and CRPF cluster expansion. Monte Carlo simulations show that the distribution of USF energies is significantly affected by chemical composition and temperature. The formalism is broadly applicable to arbitrary crystal structures and alloy chemistries and will enable the development of rigorous models for deformation mechanisms in high-entropy alloys.

Natarajan, Anirudh Raju (ORCID:0000000202529480)↗

CFT and lattice correlators near an RG domain wall between minimal models

Conformal interfaces separating two conformal field theories (CFTs) provide maps between different CFTs, and naturally exist in nature as domain walls between different phases. One particularly interesting construction of a conformal interface is the renormalization group (RG) domain wall between CFTs. For a given Virasoro minimal model \mathcal{M}_{k+3,k+2} ℳ k + 3 , k + 2 , an RG domain wall can be generated by a specific deformation which triggers an RG flow towards its adjacent Virasoro minimal model \mathcal{M}_{k+2,k+1} ℳ k + 2 , k + 1 with the deformation turned on over part of the space. An algebraic construction of this domain wall was proposed by Gaiotto in [J. High Energy Phys. 12, 103 (2012)]. In this paper, we will provide a study of this RG domain wall for the minimal case k=2 k = 2 , which can be thought of as a nonperturbative check of the construction. In this case the wall is separating the Tricritical Ising Model (TIM) CFT and the Ising Model (IM) CFT. We will check the analytical results of correlation functions from the RG brane construction with the numerical density matrix renormalization group (DMRG) calculation using a lattice model proposed in [arXiv.1206.1332, Science 344, 280 (2014)], and find a perfect agreement. We comment on possible experimental realizations of this RG domain wall.

Cogburn, Cameron V.↗

Oceanic lithosphere and asthenosphere - Thermal and mechanical structure

A coupled thermomechanical subsolidus model of the oceanic lithosphere and asthenosphere is developed which includes vertical heat conduction, a temperature-dependent thermal conductivity, heat advection by a horizontal and vertical mass flow that depends on depth and age, contributions of viscous dissipation or shear heating, a linear or nonlinear deformation law relating shear stress and strain rate, as well as a temperature- and pressure-dependent viscosity. The model requires a constant horizontal velocity and temperature at the surface, but zero horizontal velocity and constant temperature at great depths. The depth- and age-dependent temperature, horizontal and vertical velocities, and viscosity structure of the lithosphere and asthenosphere are determined along with the age-dependent shear stress in those two zones. The ocean-floor topography, oceanic heat flow, and lithosphere thickness are deduced as functions of ocean-floor age; seismic velocity profiles which exhibit a marked low-velocity zone are constructed from the age-dependent geotherms and assumed values of the elastic parameters. It is found that simple boundary-layer cooling determines the thermal structure at young ages, while effects of viscous dissipation become more important at older ages.

Schubert, G.↗

Demonstration of a Segment Alignment Maintenance System on a Seven-Segment Sub-Array of the Hobby-Eberly Telescope

NASA's Marshall Space Flight Center, in collaboration with Blue Line Engineering of Colorado Springs, Colorado, is developing a Segment Alignment Maintenance System (SAMS) for McDonald Observatory's Hobby-Eberly Telescope (HET). The SAMS shall sense motions of the 91 primary mirror segments and send corrections to HET's primary mirror controller as the mirror segments misalign due to thermo-elastic deformations of the mirror support structure. The SAMS consists of inductive edge sensors supplemented by inclinometers for global radius of curvature sensing. All measurements are sent to the SAMS computer where mirror motion corrections are calculated. In October 2000, a prototype SAMS was installed on a seven-segment cluster of the HET. Subsequent testing has shown that the SAMS concept and architecture are a viable practical approach to maintaining HET's primary mirror figure, or the figure of any large segmented telescope. This paper gives a functional description of the SAMS sub-array components and presents test data to characterize the performance of the sub-array SAMS.

Rakoczy, John↗

Quantitative study of the effect of grain boundary parameters on the slip system level Hall-Petch slope for basal slip system in Mg-4Al

Several theoretical studies have reported that the geometry and structure of grain boundaries in polycrystalline materials could impose a significant effect on the Hall-Petch slope. However, experimental observations are primarily limited by the ability of the techniques to accurately quantify the grain boundary strength and validate these theoretical models. Using high-resolution electron backscatter diffraction (HR-EBSD), the local stress tensor ahead of a slip band blocked by a grain boundary was quantified and coupled with a continuum dislocation pile-up model to assess the barrier strength of specific grain boundaries to specific slip systems, referred to as micro-Hall-Petch coefficient. For basal slip system in a deformed Mg-4Al alloy, the micro-Hall-Petch coefficient ($k^{basal}_{μ}$) varied significantly, from 0.054 to 0.184 MPa - m 1/2 for nine different grain boundaries. These results were correlated with geometric descriptors of the respective grain boundaries, with three-dimensional GB profile additionally measured via focused ion beam milling. It was found that the angle between the two slip plane traces on the grain boundary plane was the most sensitive parameter affecting $k^{basal}_{μ}$, followed by the angle between the slip directions. A functional form for calculation of $k^{basal}_{μ}$ depending on these two angles is proposed to augment crystal plasticity constitutive models with slip resistance dependent on some measure of the grain size. The method allows a new pathway to calibrate grain size strengthening parameters in crystal plasticity models, allowing further computational investigations of the interrelationship between texture, grain morphology, and the Hall Petch effect.

36 MATERIALS SCIENCE↗

Restraint of Liquid Jets by Surface Tension in Microgravity Modeled

Tension in Microgravity Modeled Microgravity poses many challenges to the designer of spacecraft tanks. Chief among these are the lack of phase separation and the need to supply vapor-free liquid or liquidfree vapor to the spacecraft processes that require fluid. One of the principal problems of phase separation is the creation of liquid jets. A jet can be created by liquid filling, settling of the fluid to one end of the tank, or even closing a valve to stop the liquid flow. Anyone who has seen a fountain knows that jets occur in normal gravity also. However, in normal gravity, the gravity controls and restricts the jet flow. In microgravity, with gravity largely absent, jets must be contained by surface tension forces. Recent NASA experiments in microgravity (Tank Pressure Control Experiment, TPCE, and Vented Tank Pressure Experiment, VTRE) resulted in a wealth of data about jet behavior in microgravity. VTRE was surprising in that, although it contained a complex geometry of baffles and vanes, the limit on liquid inflow was the emergence of a liquid jet from the top of the vane structure. Clearly understanding the restraint of liquid jets by surface tension is key to managing fluids in low gravity. To model this phenomenon, we need a numerical method that can track the fluid motion and the surface tension forces. The fluid motion is modeled with the Navier-Stokes equation formulated for low-speed incompressible flows. The quantities of velocity and pressure are placed on a staggered grid, with velocity being tracked at cell faces and pressure at cell centers. The free surface is tracked via the introduction of a color function that tracks liquid as 1/2 and gas as -1/2. A phase model developed by Jacqmin is used. This model converts the discrete surface tension force into a barrier function that peaks at the free surface and decays rapidly. Previous attempts at this formulation have been criticized for smearing the interface. However, by sharpening the phase function, double gridding the fluid function, and using a higher order solution for the fluid function, interface smearing is avoided. These equations can be rewritten as two coupled Poisson equations that also include the velocity. The method of solution is as follows: first, the phase equations are solved from this solution, a velocity field is generated, then a successive overrelaxation scheme is used to solve for a pressure field consistent with the velocity solution. After the code was implemented in axisymmetric form and verified by several test cases, the drop tower runs of Aydelott were modeled. The model handed the free-surface deformation quite nicely, even to the point of modeling geyser growth in the regime where the free surface was no longer restrained. A representative run is shown.

Chato, David J.↗

Formation and Composition-Dependent Properties of Alloys of Cubic Halide Perovskites

Distinct shortcomings of individual halide perovskites for solar applications, such as restricted range of band gaps, propensity of $ABX_3$ to decompose into $AX + BX_2$, or oxidation of 2$ABX_3$ into $A_2BX_6$, have led to the need to consider alloys of individual perovskites such as [FA,Cs][Pb,Sn][Br,I] 3 . This proposition creates a nontrivial material selection problem associated with a six-component structure, spanning a continuum of three sets of compositions (one for each sublattice) and requiring control of phase separation or ordering in each alloyed subfield. Not surprisingly, material and structure choices were made thus far mostly via trial-and-error explorations among a large number of arrangements. In this work, we use ideas from solid-state theory of semiconductor alloys to analyze the behaviors of the canonical [FA,Cs][Pb,Sn]I 3 alloy system, where FA is formamidinium. Density functional calculations utilizing specially constructed supercells are used to calculate the composition dependence of band gaps, energy of decomposition, and alloy mixing enthalpies. A number of clear trends are observed for A-site alloys [Cs,FA]SnI 3 and [Cs,FA]PbI 3 , as well as for B-site alloys Cs[Sn,Pb]I 3 and FA[Sn,Pb]I 3 . To understand the physical reasons that control these trends, we decompose the alloy properties into distinct physical terms: (i) the energies associated with removing the octahedral deformations (tilting, rotations, B-site displacements) of the individual components, (ii) the energies for compressing the larger component and expanding the smaller one to the alloy volume V(x), (iii) the charge-transfer energies associated with placing the alloyed units onto a common lattice, and finally, (iv) structural relaxation of all bonds within the cells. This analysis clarifies the origin of the observed trends in band gaps, decomposition energies, and mixing enthalpies. Unlike a number of previous calculations, we find that the proper description of alloy physics requires that even the pure, nonalloyed, end-point compounds need to be allowed to develop local environment-dependent octahedral deformation that lowers significantly the total energy and raises their band gaps.

36 MATERIALS SCIENCE↗

Multiscale formulation of frictional contact mechanics at the pore scale

Direct numerical simulation (DNS) yields the highest fidelity predictions of mechanical deformation at the pore scale, but is prohibitively expensive for analyzing large or many samples. Discrete element methods (DEM) are an efficient alternative, but are limited to granular media and incapable of estimating or controlling prediction errors. We present a pore-level multiscale method (PLMM) that approximates DNS efficiently and with controllable accuracy. We focus on the linear elastic response of a consolidated geologic porous medium with arbitrary microstructure, heterogeneous mineralogy, containing cracks or defects. PLMM decomposes the solid phase into non-overlapping subdomains, on which local basis functions are constructed. The bases are then coupled with a global interface problem that accounts for slip or stick contact conditions between the subdomains. PLMM produces an initial, but accurate, approximation to DNS that can be iteratively improved. It is amenable to parallelism and allows for different mesh, models, and physics in each subdomain. An algebraic interpretation of PLMM as a preconditioner is also presented to allow non-intrusive implementation into existing solvers. Lastly, this work extends previous developments of PLMM in fluid dynamics to solid mechanics and enables future extensions towards modeling coupled flow and mechanics problems.

42 ENGINEERING↗

Three-dimensional explicit finite element formulation for shear localization with global tracking of embedded weak discontinuities

Shear localization is a damage mechanism characterized by the concentration of plastic deformation within thin bands of materials, often observed as a precursor to failure in metals. In this study, we present a three-dimensional explicit finite element formulation with embedded weak discontinuities for the treatment of shear localization under dynamic loading conditions. In this formulation, the onset of localization is detected via a material stability analysis suitable for rate-sensitive materials, and the continuity of propagating shear bands is ensured using a global tracking strategy involving a heat-conduction type boundary value problem, solved for a scalar level set function over the global domain of the problem. In addition, we propose a modified quadrature rule that is used to compute the contribution of an individual element to the global finite element arrays, taking into account the position of the embedded shear band within that element.

42 ENGINEERING↗

Correlation of finite element free vibration predictions using random vibration test data

Finite element analysis is regularly used during the engineering cycle of mechanical systems to predict the response to static, thermal, and dynamic loads. The finite element model (FEM) used to represent the system is often correlated with physical test results to determine the validity of analytical results provided. Results from dynamic testing provide one means for performing this correlation. One of the most common methods of measuring accuracy is by classical modal testing, whereby vibratory mode shapes are compared to mode shapes provided by finite element analysis. The degree of correlation between the test and analytical mode shapes can be shown mathematically using the cross orthogonality check. A great deal of time and effort can be exhausted in generating the set of test acquired mode shapes needed for the cross orthogonality check. In most situations response data from vibration tests are digitally processed to generate the mode shapes from a combination of modal parameters, forcing functions, and recorded response data. An alternate method is proposed in which the same correlation of analytical and test acquired mode shapes can be achieved without conducting the modal survey. Instead a procedure is detailed in which a minimum of test information, specifically the acceleration response data from a random vibration test, is used to generate a set of equivalent local accelerations to be applied to the reduced analytical model at discrete points corresponding to the test measurement locations. The static solution of the analytical model then produces a set of deformations that once normalized can be used to represent the test acquired mode shapes in the cross orthogonality relation. The method proposed has been shown to provide accurate results for both a simple analytical model as well as a complex space flight structure.

Chambers, Jeffrey A.↗

Influence of Material Distribution on Impact Resistance of Hybrid Composites

Impact events occur in a wide variety of circumstances. A typical example is a bullet impacting a target made of composite material. These impact events produce time-varying loads on a structure that can result in damage. As a first step to understanding the damage resistance issue in composite laminates, an accurate prediction of the transient response during an impact event is necessary. The analysis of dynamic loadings on laminated composite plates has undergone considerable development recently. Rayleigh-Ritz energy method was used to determine the impact response of laminated plates. The impact response of composite plates using shear deformation plate theory was analyzed. In recent work a closed-form solution was obtained for a rectangular plate with four edges simply supported subjected to a center impact load using classical plate theory. The problem was further investigated and the analysis results compared of both classical plate theory and shear deformation theory, and found that classical plate theory predicts very accurate results for the range of small deformations considered. In this study, the influence of cross sectional material distribution on the comparative impact responses of hybrid metal laminates subjected to low and medium velocity impacts is investigated. A simple linear model to evaluate the magnitude of the impact load is proposed first, and it establishes a relation between the impact velocity and the impact force. Then a closed-form solution for impact problem is presented. The results were compared with the finite element analysis results. For an 11 layer-hybrid laminate, the impact response as a function of material distribution in cross-section is presented. With equal areal weight, the effect of the number of laminate layers on the impact resistance is also investigated. Finally, the significance of the presented results is discussed.

Abatan, Ayu↗

Optical Performance Modeling of FUSE Telescope Mirror

We describe the Metrology Data Processor (METDAT), the Optical Surface Analysis Code (OSAC), and their application to the image evaluation of the Far Ultraviolet Spectroscopic Explorer (FUSE) mirrors. The FUSE instrument - designed and developed by the Johns Hopkins University and launched in June 1999 is an astrophysics satellite which provides high resolution spectra (lambda/Delta(lambda) = 20,000 - 25,000) in the wavelength region from 90.5 to 118.7 nm The FUSE instrument is comprised of four co-aligned, normal incidence, off-axis parabolic mirrors, four Rowland circle spectrograph channels with holographic gratings, and delay line microchannel plate detectors. The OSAC code provides a comprehensive analysis of optical system performance, including the effects of optical surface misalignments, low spatial frequency deformations described by discrete polynomial terms, mid- and high-spatial frequency deformations (surface roughness), and diffraction due to the finite size of the aperture. Both normal incidence (traditionally infrared, visible, and near ultraviolet mirror systems) and grazing incidence (x-ray mirror systems) systems can be analyzed. The code also properly accounts for reflectance losses on the mirror surfaces. Low frequency surface errors are described in OSAC by using Zernike polynomials for normal incidence mirrors and Legendre-Fourier polynomials for grazing incidence mirrors. The scatter analysis of the mirror is based on scalar scatter theory. The program accepts simple autocovariance (ACV) function models or power spectral density (PSD) models derived from mirror surface metrology data as input to the scatter calculation. The end product of the program is a user-defined pixel array containing the system Point Spread Function (PSF). The METDAT routine is used in conjunction with the OSAC program. This code reads in laboratory metrology data in a normalized format. The code then fits the data using Zernike polynomials for normal incidence systems or Legendre-Fourier polynomials for grazing incidence systems. It removes low order terms from the metrology data, calculates statistical ACV or PSD functions, and fits these data to OSAC models for the scatter analysis. In this paper we briefly describe the laboratory image testing of FUSE spare mirror performed in the near and vacuum ultraviolet at John Hopkins University and OSAC modeling of the test setup performed at NASA/GSFC. The test setup is a double-pass configuration consisting of a Hg discharge source, the FUSE off-axis parabolic mirror under test, an autocollimating flat mirror, and a tomographic imaging detector. Two additional, small fold flats are used in the optical train to accommodate the light source and the detector. The modeling is based on Zernike fitting and PSD analysis of surface metrology data measured by both the mirror vendor (Tinsley) and JHU. The results of our models agree well with the laboratory imaging data, thus validating our theoretical model. Finally, we predict the imaging performance of FUSE mirrors in their flight configuration at far-ultraviolet wavelengths.

Saha, Timo T.↗

Stretchable electronic systems with containment chambers

The present invention provides electronic systems, including device arrays, comprising functional device(s) and/or device component(s) at least partially enclosed via one or more containment chambers, such that the device(s) and/or device component(s) are at least partially, and optionally entirely, immersed in a containment fluid. Useful containment fluids for use in containment chambers of electronic devices of the invention include lubricants, electrolytes and/or electronically resistive fluids. In some embodiments, for example, electronic systems of the invention comprise one or more electronic devices and/or device components provided in free-standing and/or tethered configurations that decouple forces originating upon deformation, stretching or compression of a supporting substrate from the free standing or tethered device or device component.

Rogers, John A.↗

Spin Signature of the C 60 Fullerene Anion: A Combined X- and D-Band EPR and DFT Study

Fullerenes attract much attention in various scientific fields, but their electronic properties are still not completely understood. Here we report on a combined EPR and DFT study of the fullerene anion C 60 – in solid glassy environment. DFT calculations were used to characterize its electronic structure through spin density distribution and magnetic resonance parameters. The electron spin density is not uniformly distributed throughout the C 60 – cage but shows a pattern similar to PC 61 BM – . EPR spectroscopy reveals a rhombic g-tensor sensitive to the environment in the frozen glassy solutions, which can be rationalized by deformation of the fullerenes along lowfrequency vibrational modes upon cooling. DFT modeling confirms that these deformations lead to variation in the C 60 – g values. Furthermore, the decrease in g-tensor anisotropy with sample annealing is related to the lessening of g-tensor strain upon temperature relaxation of the most distorted sites in the glassy state.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Extreme matter meets extreme gravity: Ultraheavy neutron stars with phase transitions

The speed of sound of the matter within neutron stars may contain nonsmooth structure related to first- or higher-order phase transitions. Here we investigate what are the observable consequences of structure in the speed of sound, such as bumps, spikes, step functions, plateaus, and kinks. One of the main consequences is the possibility of ultraheavy neutron stars (with masses larger than 2.5 solar masses), mass twins in heavy (with masses larger than 2 solar masses) and ultraheavy neutron stars. These stars pass all observational and theoretical constraints, including those imposed by recent LIGO/Virgo gravitational-wave observations and NICER x-ray observations. Here we thoroughly investigate other consequences of this structure in the speed of sound to develop an understanding of how nonsmooth features affect astrophysical observables, such as stellar radii, tidal deformability, moment of inertia, and Love number. Our results have important implications for future gravitational wave and x-ray observations of neutron stars and their impact in nuclear astrophysics.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Myco-Architecture off Planet: Growing Surface Structures at Destination

Our work focused on filling major select key technical knowledge gaps, and technical aspects to be addressed in a Phase II proposal were identified such as the secretion of glues or plastics by the fungi or bacteria to form biocomposites. Specifically, we raised the technical TRL by assessing the growth of in-house mycelial-producing strains on potential food substrates and analyzed the advantages and disadvantages to their use off planet. Growth of the mycelia on sawdust and powdered nutrients including dried, powdered algae, was used as a baseline. We measured growth as a function of temperature, tested the relevant material properties of the mycelia products produced, and measured the mass of input and output volume for production per volume of material. Mycotecture was produced in a bag simulating the proposed mission implementation. Post-production such as heat treatment (as is done with terrestrial mycotecture)was assessed.The material properties of the dried, frozen and baked mycelial outputs included density, strength,thermal insulation and flame retardation. Tensile/compression testing on biocomposite material specimens was performed to obtain three-dimensional elastic constants, as well as to examine deformation and fracture behavior. These mechanical tests were performed on both dormant and activated samples to understand the evolution of the structural material. In addition, experimental data obtained from the mechanical testing was used to build a failure prediction model that accounts for material anisotropy. This aided in future structural design through a quantitative understanding of the mechanical limits of the material. To examine the texture, defects, fracture surfaces of the material, high-resolution microstructural imaging was utilized before and after the mechanical tests. This microstructural analysis informed us about the macro-structural behavior and influenced structural design. Thermal decomposition analysis was also performed tounderstand the thermal limits.

Biocomposites↗

Design principles for the ultimate gas deliverable capacity material: nonporous to porous deformations without volume change

Understanding the fundamental limits of gas deliverable capacity in porous materials is of critical importance as it informs whether technical targets (e.g., for on-board vehicular storage) are feasible. High-throughput screening studies of rigid materials, for example, have shown they are not able to achieve the original ARPA-E methane storage targets, yet an interesting question remains: what is the upper limit of deliverable capacity in flexible materials? In this work we develop a statistical adsorption model that specifically probes the limit of deliverable capacity in intrinsically flexible materials. The resulting adsorption thermodynamics indicate that a perfectly designed, intrinsically flexible nanoporous material could achieve higher methane deliverable capacity than the best benchmark systems known to date with little to no total volume change. Density functional theory and grand canonical Monte Carlo simulations identify a known metal–organic framework (MOF) that validates key features of the model. Therefore, this work (1) motivates a continued, extensive effort to rationally design a porous material analogous to the adsorption model and (2) calls for continued discovery of additional high deliverable capacity materials that remain hidden from rigid structure screening studies due to nominal non-porosity.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗