SEARCH · Engineering Papers

Results for “exemplar”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 109 records · Page 6

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Radiation-reduced diffusion of Nd in bcc Fe

Here, this study combines ab initio calculations and kinetic Monte Carlo simulations to investigate the bulk diffusion of Nd, an exemplar lanthanide element, in bcc Fe. Contrary to the conventional belief that irradiation accelerates diffusion processes, our study reveals that irradiation can unexpectedly reduce the diffusion of Nd due to the strong binding between substitutional Nd atoms and vacancies in bcc Fe. Combined with the high concentration of nonequilibrium vacancies, this promotes the abundant formation of Nd+3Va defect complexes containing a single Nd substitutional atom and three vacancies, accompanied by the depletion of Nd-vacancy pairs that mediate the diffusion of Nd.

Density functional theory (DFT)↗

DOE OSTI2017

Phase stability, physical properties and strengthening mechanisms of concentrated solid solution alloys

Here, we review recent research developments in a special class of multicomponent concentrated solid solution alloys (CSAs) – of which the recently discovered high entropy alloys (HEAs) are exemplars – that offer a new paradigm for the development of next generation structural materials. This review focuses on the role of inherent extreme chemical complexity on the phase stability, electronic, transport, and mechanical properties of this remarkable class of disordered solid solution alloys. Both experimental observations and theoretical models indicate that the phase stability of HEAs goes beyond the original conjecture that these alloys are stabilized by configurational/mixing entropy; rather, it results from competition between the homogeneously disordered phase and phase separation/intermetallic compound formation. Although the number of single-phase HEAs with equiatomic composition is limited, those that do exist often exhibit remarkable electronic, magnetic, transport, and mechanical properties. For the mechanical response, we discuss the solution strengthening mechanism which governs the strength and deformation behaviors of the CSAs, as well as the increasing evidence that low stacking fault energies (deformation twinning) plays an important role in the low temperature strength and ductility of CrMnFeCoNi related alloys. We also review the current understanding of the role of the number and type of alloy elements in determining the electronic, magnetic, and transport properties, in particular the dominant role of magnetic interactions in the properties of 3d-transition metal based alloys. Finally, we emphasize that, despite rapid progress in characterization and understanding of the phase stability and physical/mechanical responses of CSAs, there remain significant challenges to fully exploring the new paradigm that these alloys represent.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Solid face sheets enable lattice metamaterials to withstand high-amplitude impulsive loading without yielding

Owing to their ability to provide tunable mechanical responses, lattice materials are frequently studied to elucidate their response to static and dynamic loads. However, these roles are typically in opposition: static loads must be supported sufficiently far away from the onset of buckling or yielding, whereas dynamic loads are typically ameliorated by crushing of the lattice, which provides excellent energy-absorption due to the large plastic deformation accompanying densification. In contrast, this work considers the octet truss as an exemplar topology, in a structural role where it must simultaneously support static loads while enduring high-amplitude impulsive loads. This study focuses on the ability to withstand impulsive loads without yielding, an essential prerequisite to enduring dual loading. Computational studies using the ALE3D hydrocode were performed to examine the response of the octet truss under a short temporal width impulse shape associated with laser-driven shocks. A key finding was that covering the lattice with a solid face sheet and treating this face sheet thickness as a design variable allows the Taylor-like pulse to be attenuated prior to entering the weaker lattice, at the cost of added mass up front. Experimental validation was accomplished by laser-driven shock testing, using octet trusses printed out of Ti-5Al-5V-5Mo-3Cr. The results show that for a given quantity of mass, the attenuation is maximized when as much mass as possible is moved into the face sheet, leaving a more slender lattice structure. The effect of placing mass in the face sheet rather than lattice beams dominates the effect of relative density, to the point where a low-mass structure with most of the mass concentrated in the face sheet can outperform a high-mass structure with most of the mass in the lattice. Finally, by further understanding the propagation of short pulse width waves within under-dense structures, this study expand the domain of applicability of such structures, including lattice materials, to challenging dual-loading regimes spanning decades of strain rates.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Bayesian optimal experimental design for constitutive model calibration

Computational simulation is increasingly relied upon for high/consequence engineering decisions, which necessitates a high confidence in the calibration of and predictions from complex material models. However, the calibration and validation of material models is often a discrete, multi-stage process that is decoupled from material characterization activities, which means the data collected does not always align with the data that is needed. To address this issue, an integrated workflow for delivering an enhanced characterization and calibration procedure—Interlaced Characterization and Calibration (ICC)—is introduced and demonstrated. Further, this framework leverages Bayesian optimal experimental design (BOED), which creates a line of communication between model calibration needs and data collection capabilities in order to optimize the information content gathered from the experiments for model calibration. Eventually, the ICC framework will be used in quasi real-time to actively control experiments of complex specimens for the calibration of a high-fidelity material model. This work presents the critical first piece of algorithm development and a demonstration in determining the optimal load path of a cruciform specimen with simulated data. Calibration results, obtained via Bayesian inference, from the integrated ICC approach are compared to calibrations performed by choosing the load path a priori based on human intuition, as is traditionally done. The calibration results are communicated through parameter uncertainties which are propagated to the model output space (i.e. stress–strain). In these exemplar problems, data generated within the ICC framework resulted in calibrated model parameters with reduced measures of uncertainty compared to the traditional approaches.

42 ENGINEERING↗

DOE OSTI2026

Systematic analysis of melt pool dynamics in laser processing of mixed powder feedstocks

Functionally graded materials (FGMs) fabricated via additive manufacturing of blended powders offer the potential to spatially tailor properties for new technologies, such as fusion first-wall systems, turbine blades, and spacecraft. However, processing these materials is difficult due to the multiplicity of processing parameters to optimize, all of which must be changed as substrate material, powder feedstock compositions, and melt pool dynamics evolve. Here, this work systematically evaluates the qualitative and quantitative effects of these variables on the melt pool size, shape, composition, and particle distribution in an exemplar Ti-Ta system, and connects the experimental results to Marangoni flow behavior and phenomena observed in other systems. Increasing laser power linearly increases melt pool size and layer thickness, driving engineering considerations such as part/geometrical tolerances. Decreasing laser velocity changes the melt pool shape from lenticular to convex and reduces chemical homogeneity due to extreme thermal and compositional gradients between the melt pool center and boundaries. Thermophysical property differences between the powder feedstock and substrate material, as well as the directionality of the gradient, affect dilution and melt pool dynamics, which in turn affect the melt pool boundary characteristics, shape, and uniformity. Mixed powder feedstocks of intermediate compositions do not behave according to linear interpolations between single-material endpoints, instead building taller and wider melt pools. As such, it is recommended to quantify process maps for at least one intermediate composition in the FGM or multi-material system of interest to ensure optimized processing parameters, predictable melt pool sizes and shapes, and compositional and spatial precision.

Dissimilar↗

DOE OSTI2025

Development and characterization of an aerosol-generated multi-method analytical particle test material

Ceria (CeO 2 ) particles with low to ultra-low loading of nickel dopant were produced using an aerosol-based, droplet-to-particle synthesis via an in-line calcination technique. This aerosol-based synthesis method enables the production of particles with a monodisperse size distribution. These produced and well-characterized, multi-element, ceria-based particles demonstrate a material exemplar for multi-method analytical testing. They were prepared from a cerium nitrate feedstock where low loading nickel dopant was added at target Ni/(Ni + Ce) atomic percents of 1 %, 0.1 %, and 0.01 %, using a nickel nitrate spike. This methodology proved to produce ceria particles doped with a dynamic range of low to ultra-low loadings of nickel over a 24-h period, with consistent size distribution, morphology, and composition. The successful incorporation of nickel was demonstrated with bulk and single particle inductively coupled plasma mass spectroscopy and revealed notable particle-to-particle elemental homogeneity. X-ray photoelectron spectroscopy demonstrated the presence of a high concentration of nickel dopant incorporated preferentially toward the surface of the particles, and that this dopant aided oxidation of surface Ce(III) atoms to Ce(IV). These particle test materials were then validated through X-ray absorption near edge spectroscopy, comparing the ultra-low 0.01 % Ni and low 1 % Ni-doped ceria samples. This revealed a more-reduced oxidation state of the nickel with an increase in dopant concentration. Finally, this work demonstrates a synthesis and systematic characterization scheme to produce multi-method analytical test particulates.

36 MATERIALS SCIENCE↗

DOE OSTI2017

The Random Ray Method for neutral particle transport

A new approach to solving partial differential equations (PDEs) based on the method of characteristics (MOC) is presented. The Random Ray Method (TRRM) uses a stochastic rather than deterministic discretization of characteristic tracks to integrate the phase space of a problem. TRRM is potentially applicable in a number of transport simulation fields where long characteristic methods are used, such as neutron transport and gamma ray transport in reactor physics as well as radiative transfer in astrophysics. In this study, TRRM is developed and then tested on a series of exemplar reactor physics benchmark problems. The results show extreme improvements in memory efficiency compared to deterministic MOC methods, while also reducing algorithmic complexity, allowing for a sparser computational grid to be used while maintaining accuracy.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2025

Magnetic pair distribution function and half polarized neutron powder diffraction at the HB-2A powder diffractometer

Local magnetic order and anisotropy are often central for understanding fundamental behavior and emergent functional properties in quantum materials and beyond. Advances in neutron powder diffraction experiments and analysis tools now allow for quantitative determination. Here, we demonstrate this here with complementary total neutron scattering and polarized neutron measurements on the HB-2A neutron powder diffractometer at the High Flux Isotope Reactor (HFIR). In recent years, magnetic pair distribution function (mPDF) analysis has emerged as a powerful technique for probing local magnetic spin ordering of magnetic materials. This method can be broadly applied to any magnetic material but is particularly effective for studying systems with short-range magnetic order, such as materials with reduced dimensionality, geometrically frustrated magnets, thermoelectrics, multiferroics, and correlated paramagnets. Magnetic anisotropy often underpins the short-range order adopted. Half-polarized neutron powder diffraction (pNPD) can be used to determine the local susceptibility tensor on the magnetic sites to quantify the magnetic anisotropy. Combining the techniques of mPDF and pNPD can therefore provide valuable insights into local magnetic behavior. A series of measurements optimized for these techniques are presented as exemplar cases focused on frustrated materials where short-range order dominates, these include measurements to ultra-low temperature (<100 mK) not typically accessible for such experiments.

Half polarized neutron scattering↗

DOE OSTI2021

Development of small-scale plutonium electrorefining in molten CaCl 2

Molten salt extraction and electrorefining of plutonium were successfully conducted with only a 200-g initial charge. Due to this small scale, the electrorefining cell components were redesigned relative to conventional, multi-kg standards. The efficiency and yield of this small-scale process are comparable to those of larger-scale exemplars. Chemical analysis of the plutonium product and the anode heel is presented. Here, the product was 99.8% pure. The anode heel was able to be depleted beyond the expected freezing point without significant impurity carryover.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2024

Dual X-ray computed tomography-aided classification of melt pool boundaries and flaws in crept additively manufactured parts

In metal additive manufacturing (AM), understanding the process-structure-performance relationships requires a combination of multi-scale characterization techniques that allows for the measurement of the melt pool shape and boundary and classifying various defects and flaws in the AM parts. Such approaches can be destructive, only 2D in nature, or have a small field of view and can be complex to co-register and analyze. Here, in this work, we present a non-destructive 3D inspection technique that employs dual-energy X-ray computed tomography (XCT) along with a model-based iterative reconstruction (MBIR) and a new segmentation algorithm. The proposed approach and algorithm are not only capable of classifying and quantifying flaws such as pores, cracks, and inclusions, but they also allow for the extraction of microstructural features such as melt pool boundaries (MPB) and melt pool regions (MPR), that can help understand process-structure-performance relationships for alloys under study. As an exemplar application, we employed the method for characterization of an additively manufactured aluminum alloy crept under tensile stress at 300 °C for 1064 h. Our results demonstrate high quality segmentation and classification of various flaws and MPB and MPR, for the first time, using 3D X-ray CT inspection. The delineated MPB and MPR in the crept samples reveal the preferential growth paths of cracks that formed during creep deformation. The technique was used for successfully quantifying the characteristics (number of defects, their density, volume fraction, etc.) of the manufacturing-induced pores and creep-induced cracks, which is necessary to better understand the creep failure mechanisms of the material.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Automation and optimization of stopping and range of ions in matter simulation runtime

Prior to every ion implantation experiment a simulation of the ion range and other relevant parameters is performed using Monte-Carlo based codes. Although increasing computing power has improved the speed of these calculations, the demands on Monte-Carlo codes are also increasing, requiring evaluation of the optimal number of simulations while ensuring accuracy within threshold bounds. We evaluate the “Stopping and Range of Ions in Matter” (SRIM) code due to its widespread usage. We show how dividing simulations into multiple parallel simulations with different random seeds can lead to calculation speedup and find lower bounds for the required number of ion traces simulated based on an exemplar system of a Ga focused ion beam and a high energy C beam as used in high linear energy transfer testing. Here our results indicate simulations can yield results within the underlying data accuracy of SRIM at 10X and 100X shorter simulation time than the SRIM default values.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2018

Using strontium isotopes to evaluate the spatial variation of groundwater recharge

Recharge of alluvial aquifers is a key component in understanding the interaction between floodplain vadose zone biogeochemistry and groundwater quality. The Rifle Site (a former U-mill tailings site) adjacent to the Colorado River is a well-established field laboratory that has been used for over a decade for the study of biogeochemical processes in the vadose zone and aquifer. This site is considered an exemplar of both a riparian floodplain in a semiarid region and a post-remediation U-tailings site. In this paper we present Sr isotopic data for groundwater and vadose zone porewater samples collected in May and July 2013 to build a mixing model for the fractional contribution of vadose zone porewater (i.e. recharge) to the aquifer and its variation across the site. The vadose zone porewater contribution to the aquifer ranged systematically from 0% to 38% and appears to be controlled largely by the microtopography of the site. The area-weighted average contribution across the site was 8% corresponding to a net recharge of 7.5 cm. Given a groundwater transport time across the site of ~1.5 to 3 years, this translates to a recharge rate between 5 and 2.5 cm/yr, and with the average precipitation to the site implies a loss from the vadose zone due to evapotranspiration of 83% to 92%, both ranges are in good agreement with previously published results by independent methods. A uranium isotopic ( 234 U/ 238 U activity ratios) mixing model for groundwater and surface water samples indicates that a ditch across the site is hydraulically connected to the aquifer and locally significantly affects groundwater. Groundwater samples with high U concentrations attributed to natural bio-reduced zones have 234 U/ 238 U activity ratios near 1, suggesting that the U currently being released to the aquifer originated from the former U-mill tailings.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

Efficient prediction of concentrating solar power plant productivity using data clustering

Concentrating solar power (CSP) plants convert solar energy to electricity and can be deployed with a thermal storage capability to shift electricity generation from time periods with available solar resource to those with high electricity demand or electricity price. Rigorous optimization of plant design and operational strategies can improve the market-competitiveness and commercial viability; however, such optimization may require hundreds of annual performance simulations, each of which can be computationally expensive when including considerations such as optimization of dispatch scheduling, sub-hourly time resolution, and stochastic effects due to uncertain weather or electricity price forecasts. This paper proposes a methodology to reduce the computational burden associated with simulation of electricity yield and revenue for CSP plants over a single- or multi-year period. Data-clustering techniques are employed to select a small number of limited-duration time blocks for simulation that, when appropriately weighted, can reproduce generation and revenue over a single year or within each year of a multi-year period. After selection of appropriate data features and weighting factors defining similarity between time-series profiles, the methodology captured annual revenue within 2.3%, 1.7%, or 1.2% using simulation of 10, 30, or 50 three-day exemplar time blocks, respectively, for each of three single-year location/weather/market scenarios and five plant configurations ranging from low to high solar multiple and storage capacity. When applied to multi-year datasets, the proposed methodology can capture inter-year variability that is unavailable from typical meteorological year (TMY) datasets while simultaneously requiring simulation of less than a single year of data.

14 SOLAR ENERGY↗

DOE OSTI2023

Evaluating the pressure dependence of PZT structures using a virtual reality environment

Pb–Zr–Ti–O (PZT) perovskites span a large solid-solution range and have found widespread use due to their piezoelectric and ferroelectric properties that also span a large range. Crystal structure analysis via Rietveld refinement facilitates materials analysis via the extraction of the structural parameters. These parameters, often obtained as a function of an additional dimension (e.g., pressure), can help to diagnose materials response within a use environment. Often referred to as “ in-situ ” studies, these experiments provide an abundance of data. Viewing structural changes due to applied pressure conditions can give much-needed insight into materials performance. However, challenges exist for viewing/presenting results when the details are inherently three-dimensional (3D) in nature. For PZT perovskites, the use of polyhedra (e.g., Zr/Ti–O 6 octahedra) to view bonding/connectivity is beneficial; however, the visualization of the octahedra behavior with pressure dependence is less easily demonstrated due to the complexity of the added pressure dimension. We present a more intuitive visualization by projecting structural data into virtual reality (VR). We employ previously published structural data for Pb 0.99 (Zr 0.95 Ti 0.05 ) 0.98 Nb 0.02 O 3 as an exemplar for VR visualization of the PZT R3c crystal structure between ambient and 0.62 GPa pressure. This is accomplished via our in-house CAD2VR™ software platform and the new CrystalVR plugin. The use of the VR environment enables a more intuitive viewing experience, while enabling on-the-fly evaluation of crystal data, to form a detailed and comprehensive understanding of in-situ datasets. Discussion of methodology and tools for viewing are given, along with how recording results in video form can enable the viewing experience.

Rodriguez, Mark A. (ORCID:0000000269405044)↗

DOE OSTI2024

Characterization of MoS 2 films via simultaneous grazing incidence X-ray diffraction and grazing incidence X-ray fluorescence (GIXRD/GIXRF)

Physical vapor deposited (PVD) molybdenum disulfide (nominal composition MoS 2 ) is employed as a thin film solid lubricant for extreme environments where liquid lubricants are not viable. The tribological properties of MoS 2 are highly dependent on morphological attributes such as film thickness, orientation, crystallinity, film density, and stoichiometry. These structural characteristics are controlled by tuning the PVD process parameters, yet undesirable alterations in the structure often occur due to process variations between deposition runs. Nondestructive film diagnostics can enable improved yield and serve as a means of tuning a deposition process, thus enabling quality control and materials exploration. Grazing incidence X-ray diffraction (GIXRD) for MoS 2 film characterization provides valuable information about film density and grain orientation (texture). However, the determination of film stoichiometry can only be indirectly inferred via GIXRD. The combination of density and microstructure via GIXRD with chemical composition via grazing incidence X-ray fluorescence (GIXRF) enables the isolation and decoupling of film density, composition, and microstructure and their ultimate impact on film layer thickness, thereby improving coating thickness predictions via X-ray fluorescence. We have augmented an existing GIXRD instrument with an additional X-ray detector for the simultaneous measurement of energy-dispersive X-ray fluorescence spectra during the GIXRD analysis. This combined GIXRD/GIXRF analysis has proven synergetic for correlating chemical composition to the structural aspects of MoS 2 films provided by GIXRD. We present the usefulness of the combined diagnostic technique via exemplar MoS 2 film samples and provide a discussion regarding data extraction techniques of grazing angle series measurements.

MoS2↗

DOE OSTI2020

Dynamics of Hydroxyl Anions Promotes Lithium Ion Conduction in Antiperovskite Li 2 OHCl

Li 2 OHCl is an exemplar of the antiperovskite family of ionic conductors, for which high ionic conductivities have been reported, but in which the atomic-level mechanism of ion migration is unclear. The stable phase is both crystallographically defective and disordered, having similar to 1/3 of the Li sites vacant, while the presence of the OH(- )anion introduces the possibility of rotational disorder that may be coupled to cation migration. In this work, complementary experimental and computational methods are applied to understand the relationship between the crystal chemistry and ionic conductivity in Li 2 OHCl , which undergoes an orthorhombic to cubic phase transition near 311 K (approximate to 38 degrees C) and coincides with the more than a factor of 10 change in ionic conductivity (from 1.2 x 10 -5 mS/cm at 37 degrees C to 1.4 x 10 -3 mS/cm at 39 degrees C). X ray and neutron experiments conducted over the temperature range 20-200 degrees C, including diffraction, quasi-elastic neutron scattering (QENS), the maximum entropy method (MEM) analysis, and ab initio molecular dynamics (AIMD) simulations, together show conclusively that the high lithium ion conductivity of cubic Li2OHCl is correlated to "paddlewheel" rotation of the dynamic OH - anion. The present results suggest that in antiperovskites and derivative structures a high cation vacancy concentration combined with the presence of disordered molecular anions can lead to high cation mobility.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Accelerating Embedding Potential Optimization by Reconstructing the Pseudo-Valence Electron Density

Density functional embedding theory (DFET) enables use of electronic structure methods with higher accuracy than density functional theory in a local region, with applications thus far ranging from (photo/electro)catalysis to reactions in solution. DFET partitions a large collection of atoms into smaller groups that interact via a shared embedding (interaction) potential V emb , determined via functional optimization. The optimized effective potential (OEP) process used to optimize V emb is time-consuming and becomes a computational bottleneck due to sharp, oscillating features of V emb near nuclei. Here, similar to pseudopotential theory, by reconstructing electron densities used in the OEP process from smoother pseudo-valence-only (PVO) electron densities as proxies for total densities of the full system and subsystems, we can retain accuracy in the embedded electronic structure calculations while potentially reducing the overhead of V emb construction, within the projector augmented-wave (PAW) formalism. We explore three different chemical reactions as exemplars to test PVO–DFET, namely, H 2 dissociative adsorption on a Cu(111) surface, H 2 O adsorption on a Pt(111) surface, and aqueous [Ca 2+ –SO 4 2– ] ion-pair formation. The PVO approximation works well for all three systems with minimal loss of accuracy (∼10–70 meV error relative to the original exact-derivative (ED) approach) while accelerating V emb generation for the Cu and Pt systems respectively by 20× and 5×. Given proper numerical convergence parameters, the spatial distributions of differences between PVO- and ED-based V emb outside the core regions are small, explaining the exceptional agreement between the two approaches. Finally, we anticipate that this more efficient PVO–DFET approximation will be useful whenever computation of V emb is much more expensive than subsequent embedded high-level electron correlation calculations.

approximation↗

DOE OSTI2022

Temperature Effects in Flexible Adsorption Processes for Amorphous Microporous Polymers

A collection of atomistic molecular simulations is reported that illustrate the impact of adsorption temperature on species uptake and adsorbate-induced structural rearrangement for amorphous polymers of intrinsic microporosity. Temperature-sensitive structural rearrangement is evaluated by contrasting two methods: standard grand canonical Monte Carlo simulations using a rigid framework approximation and a combined Monte Carlo/molecular dynamics approach that fully incorporates framework flexibility. We report single-component gas phase adsorption isotherms for CH 4 , C 2 H 4 , C 2 H 6 , C 3 H 6 , C 3 H 8 , and CO 2 across a temperature range of 250–400 K for models of an archetypal polymer of intrinsic microporosity, PIM-1. A quadratic model is presented that captures two main mechanisms of temperature-dependent adsorption-induced deformation of PIM-1 up to a relative swelling of 1.15: thermal expansion and an increased propensity to swell as a function of species uptake. Two case studies are reported that highlight the critical role of operating temperature in industrial storage and separation applications. The first study focuses on methane storage and delivery applications using a pressure–temperature swing adsorption application (PTSA). We demonstrate that larger working capacities are accompanied by increased volumetric strain between adsorption–desorption steps. The second case study considers PIM-1 as an adsorbent to separate an exemplar ternary syngas mixture at operating temperatures ranging 300–550 K. Here, a temperature threshold of ~400 K is identified, beyond which adsorption-induced PIM-1 swelling is negligible and the solubility selectivity-loading curve transitions to exhibiting a nearly linear relationship.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗