Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “linear”

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 163 records · Page 9

Breaking the Linear Scaling Relations for the Oxygen Reduction Reaction with a Dual‐Atom Catalyst Composed of a MnFe‐Porphyrrole Aerogel

Bimetallic catalysts offer enhanced catalytic performance through synergistic interactions between the two metals, allowing them to break the linear scaling relations and reach high electrocatalytic activity. This study presents bimetallic aerogel-based catalyst synthesized as a covalent, three-dimensional framework containing neighboring iron and manganese sites. The aerogel structure provides a high surface area and porosity, facilitating an ultra-high active site density and efficient mass transport. The MnFe porphyrrole's unique structure is obtained by alternately linking Mn-porphyrin and Fe-corrole complexes. It exhibited outstanding performance with an onset potential of 0.99 V RHE . Comparative studies with a free-base Fe porphyrrole catalyst (E onset 0.97 V RHE ) revealed that while Mn incorporation led to only a slight improvement in half-cell performance, it resulted in significantly enhanced performance in anion exchange membrane fuel cell. The MnFe catalyst achieved an OCV of 0.97 V and a peak power density of 0.27 W cm −2 , outperforming the free-base Fe counterpart. Using density functional theory calculations, we show that the higher ORR activity of MnFe-porphyrrole is due to charge transfer between Mn and Fe atoms, which is absent in the reference free-base Fe-porphyrrole. These findings underscore the advantages of bimetallic catalysts in improving ORR activity and fuel cell efficiency by leveraging synergistic effects.

Aerogel↗

Linear Dichroism of the Optical Properties of SnS and SnSe Van der Waals Crystals

Abstract Tin monochalcogenides SnS and SnSe, belonging to a family of Van der Waals crystals isoelectronic to black phosphorus, are known as environmentally friendly materials promising for thermoelectric conversion applications. However, they exhibit other desired functionalities, such as intrinsic linear dichroism of the optical and electronic properties originating from strongly anisotropic orthorhombic crystal structures. This property makes them perfect candidates for polarization‐sensitive photodetectors working in near‐infrared spectral range. A comprehensive study of the SnS and SnSe crystals is presented, performed by means of optical spectroscopy and photoemission spectroscopy, supported by ab initio calculations. The studies reveal the high sensitivity of the optical response of both materials to the incident light polarization, which is interpreted in terms of the electronic band dispersion and orbital composition of the electronic bands, dictating the selection rules. From the photoemission investigation the ionization potential, electron affinity and work function are determined, which are parameters crucial for the design of devices based on semiconductor heterostructures.

Chemistry↗

Efficient BiVO 4 /CoFeO x H y photoanodes using controlled annealing and conformal linear-sweep electrocatalyst photodeposition

Although monoclinic bismuth vanadate (BiVO 4 ) is a promising photoanode for solar water splitting, its practical use is hindered by imperfect photocurrent generation/collection, low photovoltage compared to the bandgap, and corrosion side reactions that limit durability. Here, we introduce a controlled-annealing sol–gel process for BiVO 4 thin-film photoanodes along with an optimized linear-sweep-voltammetry photodeposition of CoFeO x H y cocatalysts. The resulting BiVO 4 films annealed at 550 °C exhibited a photocurrent density of 4.1 mA/cm 2 at 1.23 V RHE under 1 sun AM 1.5G solar simulation and a low onset potential of 0.26 V RHE due to high majority carrier conductivity, a crystalline bulk with reduced defects as evidenced by x-ray photoelectron spectroscopy and photoluminescence lifetime analysis, and thus enhanced photocarrier collection. However, significant degradation in performance was found due to interfacial photocorrosion. To protect the surface and speed the oxygen-evolution reaction CoFeO x H y cocatalyst layers were deposited. By varying the number of consecutive sweeps and adjusting the applied bias range, an ultra-thin (~15 nm) CoFeO x H y cocatalyst layer was uniformly grown deposited over 30 cycles on the BiVO 4 surface. The resulting BiVO 4 /CoFeO x H y yielded 4.03 mA/cm 2 at 1.23 V RHE and onset potential of 0.24 V RHE , with stable operation (~15 % loss in photocurrent at 1.23 V RHE relative to ~60 % loss in the uncatalyzed control sample). These conformal CoFeO x H y catalytic layers function simultaneously to selectively collect photoexcited holes from the BiVO 4 , catalyze the water-oxidation reaction, and protect the BiVO 4 from photodegradation.

42 ENGINEERING↗

Hempseed cell wall polysaccharides are dominated by linear xylans and cellulose: Comprehensive structural profiling of ten cultivars of industrial hemp, Cannabis sativa L .

Hempseed is a rich source of dietary fiber; however, there has been limited research on the variability of carbohydrate composition in hempseed cell walls. The primary aim of this study was to conduct a comprehensive chemical and structural analysis of the cell wall polysaccharides in ten hempseed cultivars. Water-soluble polysaccharides (WSP) and water-insoluble residues (WIR) were isolated and subsequently analyzed for their monosaccharide composition using HPAEC-PAD, glycosyl linkage analysis using GC–MS, and structural characterization via NMR spectroscopy. All hempseed cultivars contained a high proportion of insoluble fibers and smaller amounts of soluble polysaccharides. Glucose and xylose were the most abundant components of the WIR fractions, while the WSP fractions contained abundant amounts of galactose, galacturonic acid, arabinose, rhamnose, and mannose. The results of linkage and spectroscopic analysis were consistent with the compositional analysis, identifying cellulose and acetylated linear xylans as primary components of WIR, and arabinogalactans, rhamnogalacturonans, heteromannans, xyloglucans, and arabinan as predominant in WSP. Altogether, the study revealed a comparable cell wall structure among the analyzed hemp seed varieties. The high fiber content of whole hempseed-based ingredients presents significant potential for food manufacturers seeking to develop products with enhanced dietary fiber content, offering both functional and nutritional benefits for consumers.

59 BASIC BIOLOGICAL SCIENCES↗

Upcycling linear low-density polyethylene waste to turbostratic graphene for high mass loading supercapacitors

Linear low-density polyethylene (LLDPE) waste is difficult to upcycle into more valuable carbon materials because it tends to completely decompose into small molecules during thermal processing. In this work, LLDPE is upcycled into a high quality turbostratic graphene using a pre-treatment step to oxidatively crosslink the polymer with the assistance of solid additives (KCl and K 2 CO 3 ) that improve crosslinking by increasing the effective surface area of the polymer melt during processing. After this pretreatment step, the crosslinked polymer could then be carbonized and catalytically graphenized between 400–950°C without decomposition of the polymer feedstock. The LLDPE derived graphene (LLDPE-G) obtained from this process has a Brunauer–Emmett–Teller (BET) specific surface area, up to 1800m 2 g -1 and average Raman I D /I G and I 2D /I G ratios of 0.85 and 0.57, respectively, indicating high quality graphene. When used as an electrode material in symmetric supercapacitors, LLDPE-G possesses a specific capacitance up to 175Fg -1 at a mass loading of 20mgcm -2 , which is two times the commercial requirement, yielding an areal capacitance of 3.5 Fcm -2 . Moreover, LLDPE-G exhibits exceptional cycling stability with a capacitance retention of 95.8 % after 100,000 cycles at a current density of 4.0 Ag -1 . Additionally, the KCl and K 2 CO 3 solids are recycled and reused over 3 complete reaction cycles to make new LLDPE-G with the material quality and electrocapacitive performance retained and verified after each cycle. Our approach creates new opportunities for upcycling waste LLDPE and other varieties of polyethylene into a higher value graphene used for electrochemical energy storage applications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Optimizing district energy systems by integrating Borehole Thermal Energy Storage Using a Mixed-Integer Linear Programming g-function framework with a Multi-Timescale Rolling Horizon method

Shallow geothermal has gained increasing attention in recent years; however, a reliable framework for its accurate incorporation into large-scale energy system optimization remains lacking. This study proposes a Mixed-Integer Linear Programming (MILP) framework combined with the g-function approach to integrate Borehole Thermal Energy Storage (BTES) technology into energy system optimization. Validation against a Modelica-based reservoir network simulation demonstrates that the proposed framework effectively captures the ground thermal response under varying energy loads and accurately estimates the borefield energy supply. To enhance scalability, a Rolling Horizon with Multi-Timescale (RH-MTS) method is further introduced, reducing computational time by 73 % for the 1-year optimization model with only minor loss of optimality. The framework is demonstrated through the case study of the UC Berkeley campus. Results indicate that BTES is a cost-effective and low-carbon solution: two borefields comprising 382 boreholes can meet 8.0 % and 6.6 % of the total campus heating and cooling demand, respectively, at an average energy rate of 0.70–0.77 USD/kWh and carbon intensity of 0.54 kg-CO2/kWh. Short-term analysis reveals a 35%–65% decline in BTES energy flow after 3–6 months of continuous heating/cooling operation, while long-term simulation shows that annual energy production of BTES can vary by up to 12.0 % after four years before stabilizing. Overall, this study develops a novel optimization framework that couples physics-based g-function method with MILP optimization framework, thereby advancing methodological development for shallow-geothermal integration and providing actionable guidance for BTES deployment in district-energy systems.

Yang, Jiahui↗

MITgcm-AD v2: Open source tangent linear and adjoint modeling framework for the oceans and atmosphere enabled by the Automatic Differentiation tool Tapenade

The Massachusetts Institute of Technology General Circulation Model (MITgcm) is widely used by the climate science community to simulate planetary atmosphere and ocean circulations. A defining feature of the MITgcm is that it has been developed to be compatible with an algorithmic differentiation (AD) tool, TAF, enabling the generation of tangent-linear and adjoint models. These provide gradient information which enables dynamics-based sensitivity and attribution studies, state and parameter estimation, and rigorous uncertainty quantification. Importantly, gradient information is essential for computing comprehensive sensitivities and performing efficient large-scale data assimilation, ensuring that observations collected from satellites and in-situ measuring instruments can be effectively used to optimize a large uncertain control space. As a result, the MITgcm forms the dynamical core of a key data assimilation product employed by the physical oceanography research community: Estimating the Circulation and Climate of the Ocean (ECCO) state estimate. Although MITgcm and ECCO are used extensively within the research community, the AD tool TAF is proprietary and hence inaccessible to a large proportion of these users. The new version 2 (MITgcm-AD v2) framework introduced here is based on the source-to-source AD tool Tapenade, which has recently been open-sourced. Another feature of Tapenade is that it stores required variables by default (instead of recomputing them) which simplifies the implementation of efficient, AD-compatible code. The framework has been integrated with the MITgcm model’s main branch and is now freely available.

Adjoints↗

Non-linear elastic-plastic behavior of the invert glass lithium phosphorous oxynitride (LiPON)

Micro and nano-scale mechanical behavior of network/ionic glasses is dictated by their composition and the number of constraints per structural unit. From this perspective, the glass LiPON presents the opportunity for enhancement of the microscale ductility due to its rather unconstrained orthophosphate structure. Here, we use instrumented nanoindentation with different tip geometries to investigate the mechanical response of LiPON glass. The results reveal that the elastic modulus of LiPON is not constant and depends on pressure. With the method utilizing spherical nanoindentation and continuous stiffness measurement (CSM) we determine the yield point of LiPON. We propose the Drucker-Prager type of yield criterion for LiPON and estimate its yield stress in compression as 2.4 GPa. There exists another stress threshold, however, around 490 MPa, at which the elastic deformation becomes non-linear, and this can be mistaken for the yield stress.

Constitutive model↗

Multi-parametric analysis for mixed integer linear programming: An application to transmission upgrade and congestion management

Upgrading the capacity of existing transmission lines is essential for meeting the growing energy demands, facilitating the integration of renewable energy, and ensuring the security of the transmission system. This study focuses on the selection of lines whose capacities and by how much should be expanded from the perspective of the Independent System Operators (ISOs) to minimize the total system cost. We employ advanced multi-parametric programming and an enhanced branch-and-bound algorithm to address complex mixed-integer linear programming (MILP) problems, considering multi-period time constraints and physical limitations of generators and transmission lines. To characterize the various decisions in transmission expansion, we model the increased capacity of existing lines as parameters within a specified range. This study first relaxes the binary variables to continuous variables and applies the Lagrange method and Karush-Kuhn-Tucker (KKT) conditions to obtain optimal solutions and identify critical regions associated with active and inactive constraints. Moreover, we extend the traditional branch-and-bound (B&B) method by determining the problem’s upper and lower bounds at each node of the B&B decision tree, helping to manage computational challenges in large-scale MILP problems. Here, we compare the difference between the upper and lower bounds to obtain an approximate optimal solution within the decision-makers’ tolerable error range. In addition, the first derivative of the objective function on the parameters of each line is used to inform the selection of lines for easing congestion and maximizing social welfare. Finally, the capacity upgrades are selected by weighing the reductions in system costs against the expense of upgrading line capacities. The findings are supported by numerical simulations and provide transmission-line planners with decision-making guidance.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Protocol to detect dilution cycles in chemostat experiments and estimate growth rate slopes with linear modeling with R software chemostat_regression

Chemostat growth chambers measure optical density over time and require manual calculation of growth rates. Here, we present chemostat_regression, R software that enables users to automatically identify chemostat cycles and estimate growth rate using a linear regression approach. We describe steps for creating requisite software environment(s), formatting input data, executing the software via command line/RStudio/R-Shiny, interpreting results, assessing the validity of results, and modifying input parameters.

59 BASIC BIOLOGICAL SCIENCES↗

Polymorphism and Negative Linear Compressibility in Pyrazine-d4

Pyrazine (1,4-diazine) is a widely used linker in coordination polymers, molecular magnets, and metal–organic frameworks, yet its solid-state polymorphism is controlled by weak C–H···N interactions. Using variable-temperature and high-pressure neutron diffraction, a previously unrecognized low-temperature transition to phase IV is identified, showing that the form stabilized under pressure can also be accessed below ∼90 K at ambient pressure in the perdeuterated compound. The transition from phase III → IV involves a symmetry-breaking lattice distortion driven by the reorganization of C–H···N hydrogen bonds, accompanied by negative linear compressibility, which is primarily accommodated through the reduction in void-space within the structure. Structural analysis reveals that while the overall monoclinic framework of phase IV is preserved under both temperature and pressure routes, the evolution of hydrogen-bonded network is distinct depending on the pathway taken. These results establish pyrazine as a model molecular solid for studying pressure–temperature polymorphism and demonstrate how weak directional interactions can produce anomalous mechanical responses in organic crystals.

Ridley, Christopher J. [Oak Ridge National Laborat↗

Linear Free Energy Relationships Associated with Hydride Transfer From [(6,6'-R 2 -bpy)Re(CO) 3 H]: A Cautionary Tale in Identifying Hydrogen Bonding Effects in the Secondary Coordination Sphere

Here, six rhenium hydride complexes, [(6,6'-R 2 -bpy)Re(CO) 3 H] (bpy = 2,2'-bipyridine, R = OEt, OMe, NHMe, Me, F, Br), were synthesized. These complexes insert CO 2 to form rhenium formate complexes of the type [(6,6'-R 2 -bpy)Re(CO) 3 {OC(O)H}]. All the rhenium formate species were characterized using X-ray crystallography, which revealed that the bpy ligand is not coplanar with the metal coordination plane containing the two nitrogen donors of the bpy ligand but tilted. A solid-state structure of [(6,6'-Me 2 -bpy)Re(CO) 3 H] determined using MicroED also featured a tilted bpy ligand. The kinetics of CO 2 insertion into complexes of the type [(6,6'-R 2 -bpy)Re(CO) 3 H] were measured experimentally and the thermodynamic hydricities of [(6,6'-R 2 -bpy)Re(CO) 3 H] species were determined using theoretical calculations. A Brønsted plot constructed using the experimentally determined rate constants for CO 2 insertion and the calculated thermodynamic hydricities for [(6,6'-R 2 -bpy)Re(CO) 3 H] revealed a linear free energy relationship (LFER) between thermodynamic and kinetic hydricity. This LFER is different to the previously determined relationship for CO 2 insertion into complexes of the type [(4,4'-R 2 -bpy)Re(CO) 3 H]. At a given thermodynamic hydricity, CO 2 insertion is faster for complexes containing a 6,6'-substituted bpy ligand. This is likely in part due to the tilting observed for systems with 6,6'-substituted bpy ligands. Notably, the 6,6'-(NHMe) 2 -bpy ligand could in principle stabilize the transition state for CO 2 insertion via hydrogen bonding. This work shows that if only the rate of CO 2 insertion into [(6,6'-(NHMe) 2 -bpy)Re(CO) 3 H] is compared to [(4,4'-R 2 -bpy)Re(CO) 3 H] systems, the increase in rate could be easily attributed to hydrogen bonding, but in fact all 6,6'-substituted systems lead to faster than expected rates.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Linearized Pair-Density Functional Theory with Spin–Orbit Coupling

Here, we include spin–orbit coupling (SOC) effects in linearized pair-density functional theory (L-PDFT), which is a multistate extension of multiconfiguration pair-density functional theory (MC-PDFT). Both 1-electron and 2-electron SOC integrals are computed using Breit-Pauli and Douglas–Kroll–Hess Hamiltonians in the atomic mean-field approximation. SO-L-PDFT removes the unphysical J-symmetry breaking observed in MC-PDFT. The accuracy of SO-L-PDFT is validated by calculations of zero-field splittings, fine-structure excitation energies, and low-energy excited-state spectra for a diverse group of atoms and molecules spanning the whole range of the periodic table, including atoms of groups 3, 11, and 13–17, the Ce 3+ and U 5+ ions, group 16 monohydrides, group 17 monoxides, lanthanide hexachlorides ([CeCl 6 ] 3− , [PrCl 6 ] 3− , and [NdCl 6 ] 3− ), actinyl ions ([UO 2 ] + , [NpO 2 ] 2+ ), and tricarbonatoactinyl complexes ([UO 2 (CO 3 ) 3 ] 5− , [NpO 2 (CO 3 ) 3 ] 4− ). We also compare the results to new spin–orbit-inclusive calculations by single-state and multistate multireference perturbation theory.

Hamiltonians↗

Low-Lying Excited States of Linear All- Trans Polyenes: Insights from Analytic Gradient and Nonadiabatic Coupling Calculations Based on Multireference Configuration Interaction

Polyenes serve as a rigorous test for theoretical models and electronic structure methods, playing a key role in advancing computational and theoretical chemistry. Here, we present a high-level theoretical investigation of linear, all-trans polyenes using energy gradients and nonadiabatic coupling vectors based on an MR-CISD wave function to describe electronic transitions involving the ground state (1 1 A g – ) and three low-lying excited states (2 1 A g – , 1 1 B u + , and 2 1 B u – ) of hexatriene, octatetraene, and decapentaene. This approach enables accurate evaluation of both adiabatic and vertical excitation and emission energies, yielding results in excellent agreement with experiment, as well as locating minima on the crossing seam between adiabatic states. Our results show that vertical excitation energies to the 1 1 B u + state are blue-shifted by 0.2–0.3 eV relative to the experimental absorption maximum, whereas the vertical emission energy from the 2 1 A g – state is red-shifted by ∼0.2 eV relative to the experimental emission maximum. Upon relaxation from the Franck–Condon geometry, the 2 1 A g – state stabilizes by around 1 eV, compared to 0.2–0.3 eV for the 1 1 B u + state. An analysis of the S 1 /S 0 crossing seam in hexatriene shows that its minimum involves asymmetric backbone deformations and provides an efficient channel for ultrafast internal conversion to the ground state, consistent with the absence of detectable fluorescence in this molecule. These results demonstrate the power of analytic gradients and nonadiabatic coupling vectors based on an MR-CISD wave function for accurately characterizing the electronic structure and photophysics of polyenes.

Excited states↗

Linear-Scaling Asymmetric Triples Correction through the Solution of the DLPNO–CCSD Lambda Equations: DLPNO–CCSD(T) Λ

In this research, we derive equations for solving for the stationary points of the DLPNO–CCSD Lagrangian, in the t 1 -transformed formalism introduced earlier and as currently implemented in the P SI 4 quantum chemistry software package. These lambda equations in the local pair natural orbital basis allow for the evaluation of CCSD(T) Λ energetics with linear-scaling computational effort, also known as the asymmetric triples correction. This DLPNO–CCSD(T) Λ method allows for accurate triples contributions to be computed for larger molecules, especially in cases that CCSD(T) is known to be insufficient, such as with multireference systems and bond-breaking systems. We showcase the accuracy of our code on reaction energies, barrier heights, and noncovalent interaction energies. Also showcased are the capabilities of our code by evaluating DLPNO–CCSD(T) Λ energetics on large noncovalent dimers up to 112 atoms, as well as a rhodium catalyst complex containing 66 atoms.

Cluster chemistry↗

Quantifying the Decomposition Kinetics of Linear C 2 –C 4 Perfluoroalkanes

The pyrolysis of C 2 –C 4 linear perfluoroalkanes was investigated using laser schlieren densitometry (LS) in a diaphragmless shock tube. Mixtures of 1, 2, and 4% perfluoroethane (PFE), perfluoropropane (PFP), and perfluorobutane (PFB) in excess krypton were shock-heated to 1400–2500 K at pressures of 60, 120, and 240 Torr. The initiation step for all perfluoroalkanes involved carbon–carbon bond fission. The rate constants for unimolecular decomposition were computed using microcanonical rate theory, and these rate constants were combined into a detailed chemical kinetic mechanism. For PFE and PFP, there was one dissociation channel, while PFB dissociated via two competing bond-fission reactions. The perfluoroalkyl radicals subsequently recombined or underwent further unimolecular decomposition. In conclusion, the kinetic model accurately simulated the density gradient profiles of all perfluoroalkanes.

Guzman, Eduardo H. [Brown Univ., Providence, RI (U↗

Formation of Linear Plasmonic Heterotrimers Using Nanoparticle Docking to DNA Origami Cages

The fabrication of complex assemblies with interesting collective properties from plasmonic nanoparticles (NPs) is often challenging. While DNA-directed self-assembly has emerged as one of the most promising approaches to forming such complex assemblies, the resulting structures tend to have large variability in gap sizes and shapes, as the DNA strands used to organize these particles are flexible, and the polydispersity of the NPs leads to variability in these critical structural features. Here, we use a new strategy termed docking to DNA origami cages (D-DOC) to organize spherical NPs into a linear heterotrimer with a precisely defined geometrical arrangement. Instead of binding NPs to the exterior of the DNA templates, D-DOC binds the NPs to either the interior or the opening of a 3D cage, which significantly reduces the variability of critical structural features by incorporating multiple diametrically arranged capture strands to tether NPs. Additionally, such a spatial arrangement of the capture strand can work synergistically with shape complementarity to achieve tighter confinement. To assemble NPs via D-DOC, we developed a multistep assembly process that first encapsulates an NP inside a cage and then binds two other NPs to the openings. Microscopic characterization shows low variability in the bond angles and gap sizes. Both UV–vis absorption and surface-enhanced Raman scattering (SERS) measurements showed strong plasmonic coupling that aligned with predictions by electrodynamic simulations, further confirming the precision of the assembly. These results suggest D-DOC could open new opportunities in biomolecular sensing, SERS and fluorescence spectroscopies, and energy harvesting through the self-assembly of NPs into more complex 3D assemblies.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Mechanically Accelerated Depolymerization of Entangled Linear Polymer Melts

Mechanical forces can enhance the chemical depolymerization of synthetic polymers when shear flow accelerates chain scission. To quantify the extent of mechanically-accelerated scission, the effect of simple shear flow (duration and strength) with low Weissenberg and Deborah numbers was investigated by considering the impact of applied work in both simple shear and shear dominated mixed flows. Hydrogenated polyisoprene was chosen as a model linear, entangled system. The conditions (strain amplitude, frequency, and shearing time) necessary to increase chain scission were assessed in the rubbery melt. Shear flow accelerated chain scission at higher temperatures, suggesting an activated process. Isothermal scission versus work curves were superposed by applying shift factors a T,S , whose Arrhenius-like temperature dependence gave an apparent activation energy for chain scission of ~ 110 kJ/mol, which is likely a combination of the activation energy of viscosity and bond energy. This work provides a base for quantifying the impact of shear on depolymerization of polymer melts and highlight the connection between viscous dissipation and scission chemistry.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗