Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “highly deviated”

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 487 records · Page 27
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

A Geo‐Structurally Based Correction Factor for Apparent Dissolution Rates in Fractured Media

Abstract Field measurements of apparent geochemical weathering reaction rates in subsurface fractured porous media are known to deviate from laboratory measurements by multiple orders of magnitude. To date, there is no geologically based explanation for this discrepancy that can be used to predict reaction rates in field systems. Proposed correction factors are typically based on ad hoc characterizations related to geochemical kinetic models. Through a series of high‐fidelity reactive transport simulations of mineral dissolution within explicit 3D discrete fracture networks, we are able to link the geo‐structural attributes with reactive transport observations. We develop a correction factor to linear transition state theory for the prediction of the apparent dissolution rate based on measurable geological properties. The modified rate law shows excellent agreement with numerical simulations, indicating that geological structure could be a primary reason for the discrepancy between laboratory and field observations of apparent dissolution rates in fractured media.

58 GEOSCIENCES↗
DOE OSTI2021

Search for pair production of scalar leptoquarks decaying into first- or second-generation leptons and top quarks in proton–proton collisions at $\sqrt{s}$ = 13 TeV with the ATLAS detector

A search for pair production of scalar leptoquarks, each decaying into either an electron or a muon and a top quark, is presented. This is the first leptoquark search using ATLAS data to investigate top-philic cross-generational couplings that could provide explanations for recently observed anomalies in B meson decays. This analysis targets high leptoquark masses which cause the decay products of each resultant top quark to be contained within a single high-p T large-radius jet. The full Run 2 dataset is exploited, consisting of 139 fb -1 of data collected from proton–proton collisions at √s=13 TeV from 2015 to 2018 with the ATLAS detector at the CERN Large Hadron Collider. In the absence of any significant deviation from the background expectation, lower limits on the leptoquark masses are set at 1480 GeV and 1470 GeV for the electron and muon channel, respectively.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Oxygen diffusion in garnet: Experimental calibration and implications for timescales of metamorphic processes and retention of primary O isotopic signatures

Knowledge of oxygen diffusion in garnet is crucial for a correct interpretation of oxygen isotope signatures in natural samples. A series of experiments was undertaken to determine the diffusivity of oxygen in garnet, which remains poorly constrained. The first suite included high-pressure (HP), nominally dry experiments performed in piston-cylinder apparatus at: (1) T = 1050–1600 °C and P = 1.5 GPa and (2) T = 1500 °C and P = 2.5 GPa using yttrium aluminum garnet (YAG; Y3Al5O12) cubes. Second, HP H2O-saturated experiments were conducted at T = 900 °C and P = 1.0–1.5 GPa, wherein YAG crystals were packed into a YAG + Corundum powder, along with 18O-enriched H2O. Third, 1 atm experiments with YAG cubes were performed in a gas-mixing furnace at T = 1500–1600 °C under Ar flux. Finally, an experiment at T = 900 °C and P = 1.0 GPa was done using a pyrope cube embedded into pyrope powder and 18O-enriched H2O. Experiments using grossular were not successful. Profiles of 18O/(18O+16O) in the experimental charges were analyzed with three different secondary ion mass spectrometers (SIMS): sensitive high-resolution ion microprobe (SHRIMP II and SI), CAMECA IMS-1280, and NanoSIMS. Considering only the measured length of 18O diffusion profiles, similar results were obtained for YAG and pyrope annealed at 900 °C, suggesting limited effects of chemical composition on oxygen diffusivity. However, in both garnet types, several profiles deviate from the error function geometry, suggesting that the behavior of O in garnet cannot be fully described as simple concentration-independent diffusion, certainly in YAG and likely in natural pyrope as well. The experimental results are better described by invoking O diffusion via two distinct pathways with an inter-site reaction allowing O to move between these pathways. Modeling this process yields two diffusion coefficients (D values) for O, one of which is approximately two orders of magnitude higher than the other. Taken together, Arrhenius relationships are:logDm2s-1=-7.2(±1.3)+(-321(±32)kJmol-12.303RT) for the slow pathway, andlogDm2s-1=-5.4(±0.7)+(-321(±20)kJmol-12.303RT) for the fast pathway. We interpret the two pathways as representing diffusion following vacancy and inter-stitial mechanisms, respectively. Regardless, our new data suggest that the slow mechanism is prevalent in garnet with natural compositions, and thus is likely to control the retentivity of oxygen isotopic signatures in natural samples. The diffusivity of oxygen is similar to Fe-Mn diffusivity in garnet at 1000–1100 °C and Ca diffusivity at 850 °C. However, the activation energy for O diffusion is larger, leading to lower diffusivities at P-T conditions characterizing crustal metamorphism. Therefore, original O isotopic signatures can be retained in garnets showing major element zoning partially re-equilibrated by diffusion, with the uncertainty caveat of extrapolating the experimental data to lower temperature conditions.

Geochemistry & Geophysics↗
DOE OSTI2025

New NDA Methods for Thorium Fuel Cycle Safeguards (Final Report)

This project developed portable Neutron Resonance Transmission Analysis (pNRTA) as a new non-destructive assay (NDA) method for thorium fuel cycles safeguards and other applications where multiple isotopes must be measured when present together. pNRTA leverages epithermal neutron resonances to assay multiple safeguards-relevant isotopes (e.g., 233 U and 235 U) when they are present together in a sample. Existing techniques are challenged by this task, driving the need for new active interrogation methods. With selected detectors, pNRTA works in high gamma-ray backgrounds from fission and activation products and 232 U progeny expected in thorium fuel cycle samples. This project leveraged a pNRTA system developed at Pacific Northwest National Laboratory (PNNL) and collaboration with the Massachusetts Institute of Technology (MIT). The system uses a commercially available deuterium-tritium (DT) neutron generator at short standoff (2 m). Key achievements in this project included: first-of-a-kind pNRTA quantitative measurements of 233 U oxide samples, an assessment of neutron detector technologies suitable for pNRTA in high gamma-ray background environments, experimentally demonstrating quantitative assay of samples containing 233 U and 235 U, and modeling studies showing the applicability of pNRTA to a wide range of material forms. Further, a custom algorithm was developed at MIT, which provided mean bias of 9% and relative standard deviation of 36% in assaying 233 U, 235 U, 238 U, and 232 Th content in eight measured samples. These outcomes form a solid technical basis for pNRTA as a new promising capability for international safeguards verification that is portable, non-destructive, quantitative, and isotopic specific.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

In-situ/operando study of Cu-based nanocatalysts for CO 2 electroreduction using electrochemical liquid cell TEM

The structure of a nanocatalyst during electrocatalytic reactions often deviates from its pristine structure due to intrinsic properties, or physical and chemical adsorption at the catalytic surfaces. Taking Cu-based catalysts for CO 2 electroreduction reactions (CO 2 RR) as an example, they often experience segregation, leaching, and alloying during reactions. With the recent breakthrough development of high-resolution polymer electrochemical liquid cells, in-situ electrochemical liquid cell transmission electron microscopy (EC-TEM) alongside other advanced microscopy techniques, has become a powerful platform for revealing electrocatalysts restructuring at the atomic level. Considering the complex reactions involving electrified solid-liquid interfaces and catalyst structural evolution with intermediates, systematic studies with multimodal approaches are crucial. In this article, we demonstrate a research protocol for the study of electrocatalysts structural evolution during reactions using the in-situ EC-TEM platform. Using Cu and CuAg nanowire catalysts for CO 2 RR as model systems, we describe the experimental procedures and findings. We highlight the platform’s crucial role in elucidating atomic-scale pathways of nanocatalyst restructuring and identifying catalytic active sites, as well as avoiding potential artifacts to ensure unbiased conclusions. Using the multimodal characterization toolbox, we provide the opportunity to correlate the structure of a working catalyst with its performance. Finally, we discuss advancements as well as the remaining gap in elucidating the structural-performance relationship of working catalysts. We expect this article will assist in establishing guidelines for future investigations of complex electrochemical reactions, such as CO 2 RR and other catalytic processes, using the in-situ EC-TEM platform.

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

Spectroscopic diagnostics of high-temperature plasma in stellar coronae using Fe XXIV-XXIV K-shell lines with XRISM

The RS CVn type binary star GT Muscae was observed during its quiescence using the Resolve X-ray microcalorimeter spectrometer onboard XRISM. The main and satellite lines of the Fe XXIV-XXIV K-shell transitions were resolved for the first time from stellar sources. We conducted line ratio analysis to investigate any deviations from collisional ionization equilibrium and Maxwell electron energy distribution with a single temperature. By using five combinations of direct excitation lines and dielectronic recombination satellite lines in three line complexes (Fe He$\alpha$, Ly$\alpha$, and He$\beta$), we found that the plasma is well characterized by two-temperature thermal plasmas with temperatures of 1.7 and 4.3 keV, which is consistent with thermal broadening of Fe xxv, and the broad-band fitting results in the 1.7–10 keV band. Other forms of deviation from a single-temperature plasma, such as different ionization and electron temperatures or the $\kappa$ distribution for the electron energy distributions, are not favored, which is reasonable for stellar coronae at quiescence. This study demonstrates the utility of the Fe K-shell line ratio diagnostics to probe plasma conditions using X-ray microcalorimeters.

X-rays: stars↗
DOE OSTI2020

High precision control and deep learning-based corn stand counting algorithms for agricultural robot

This paper presents high precision control and deep learning-based corn stand counting algorithms for a low-cost, ultra-compact 3D printed and autonomous field robot for agricultural operations. Currently, plant traits, such as emergence rate, biomass, vigor, and stand counting, are measured manually. This is highly labor-intensive and prone to errors. The robot, termed TerraSentia, is designed to automate the measurement of plant traits for efficient phenotyping as an alternative to manual measurements. In this paper, we formulate a Nonlinear Moving Horizon Estimator that identifies key terrain parameters using onboard robot sensors and a learning-based Nonlinear Model Predictive Control that ensures high precision path tracking in the presence of unknown wheel-terrain interaction. Moreover, we develop a machine vision algorithm designed to enable an ultra-compact ground robot to count corn stands by driving through the fields autonomously. The algorithm leverages a deep network to detect corn plants in images, and a visual tracking model to re-identify detected objects at different time steps. We collected data from 53 corn plots in various fields for corn plants around 14 days after emergence (stage V3 - V4). The robot predictions have agreed well with the ground truth with C robot =1.02×C human -0.86 and a correlation coefficient R=0.96. The mean relative error given by the algorithm is -3.78%, and the standard deviation is 6.76%. These results indicate a first and significant step towards autonomous robot-based real-time phenotyping using low-cost, ultra-compact ground robots for corn and potentially other crops.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Anisotropic positive linear and sub-linear magnetoresistivity in the cubic type-II Dirac metal Pd 3 In 7

We report a transport study on Pd 3 In 7 which displays multiple Dirac type-II nodes in its electronic dispersion. Pd 3 In 7 is characterized by low residual resistivities and high mobilities, which are consistent with Dirac-like quasiparticles. For an applied magnetic field (μ 0 H) having a non-zero component along the electrical current, we find a large, positive, and linear in μ 0 H longitudinal magnetoresistivity (LMR). The sign of the LMR and its linear dependence deviate from the behavior reported for the chiral-anomaly-driven LMR in Weyl semimetals. Interestingly, such anomalous LMR is consistent with predictions for the role of the anomaly in type-II Weyl semimetals. In contrast, the transverse or conventional magnetoresistivity (CMR for electric fields E⊥μ 0 H) is large and positive, increasing by 10 3 –10 4 % as a function of μ 0 H while following an anomalous, angle-dependent power law ρ xx α (μ 0 H) n with n(θ) ≤ 1. The order of magnitude of the CMR, and its anomalous power-law, is explained in terms of uncompensated electron and hole-like Fermi surfaces characterized by anisotropic carrier scattering likely due to the lack of Lorentz invariance.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Multiband effects on the upper critical field angular dependence of 122-family iron pnictide superconductors

Detailed measurements of the in-plane resistivity were performed in a high-quality Ba(Fe 1-x Co x ) 2 As 2 (x=0.065) single crystal, in magnetic fields up to 9 T and with different orientations θ relative to the crystal c axis. A significant ρ(T) H,θ rounding is observed just above the superconducting critical temperature T c due to Cooper pairs created by superconducting fluctuations. These data are analyzed in terms of a generalization of the Aslamazov-Larkin approach, that extends its applicability to high reduced-temperatures and magnetic fields. This method allows us to carry out a criterion-independent determination of the angular dependence of the upper critical field, H c2 (θ). In spite of the relatively small anisotropy of this compound, it is found that H c2 (θ) presents a significant deviation from the single-band 3D anisotropic Ginzburg-Landau (3D-aGL) approach, particularly for large θ (typically above ~60o). These results are interpreted in terms of the multiband nature of these materials, in contrast with other proposals for similar H c2 (θ) anomalies. Our results are also consistent with an effective anisotropy factor almost temperature independent near T c , a result that differs from the ones obtained by using a single-band model.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Structure and pressure dependence of the Fermi surface of lithium

Here, we report studies of the Fermi surface (FS) of isotopically pure polycrystalline 7 Li from ambient pressure to 4.7 GPa. Shubnikov–de Haas (SdH) oscillations at 300 mK in external magnetic fields up to 35 T are measured to map the spherical parts of the FS of lithium. Our ambient pressure data show that the principal SdH frequencies consist of three distinct peaks at 41.25 ± 0.05, 41.65 ± 0.05, and 42.05 ± 0.05 kT, which are theoretically consistent with the presence of two crystal structure domains exhibiting nearly spherical FS. The size of the spherical parts of the FS is compatible with bcc and fcc crystal structures. The measured frequencies at 41.25 and 42.65 kT present direct quantitative evidence for the spherical deformation of the FS in fcc Li. Our high-pressure data show that while the FS of Li deforms under compression, it remains mostly spherical up to 4.7 GPa and the pressure dependence of the SdH frequency is consistent with the theoretically calculated pressure dependence in the fcc structure. Finally, we find that the electron effective mass does not deviate under pressure significantly from its ambient pressure value.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Local and global approaches to treat the torsional barriers of 4-methylacetophenone using microwave spectroscopy

The Fourier transform microwave spectrum of 4-methylacetophenone recorded from 8 GHz to 18 GHz under jet-cooled conditions has revealed large tunneling splittings arising from a low barrier to internal rotation of the ring methyl group and small splittings from a high torsional barrier of the acetyl methyl group. The large splittings are especially challenging to model, while the small splittings are difficult to analyze due to the resolution limit of 120 kHz. The combination of two methyl groups undergoing internal rotations caused each rotational transition to split into five torsional species, which were resolved and fitted using a modified version of the XIAM code and the newly devel- oped ntop code to a root-mean-square deviation close to measurement accuracy, providing an estimate of the V 3 potential barriers of about 22 cm -1 and 584–588 cm -1 for the ring and the acetyl methyl groups, respectively. The assignment was aided by separately fitting the five torsional species using odd-power order operators. Only one conformer in which all heavy atoms are located on a symmetry plane could be identified in the spectrum, in agreement with results from conformation analysis using quantum chemical calculations.

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

Elucidating Structural Disorder in Ultra‐Thin Bi‐Rich Bismuth Oxyhalide Photocatalysts

Abstract Advancing the field of photocatalysis requires the elucidation of structural properties that underpin the photocatalytic properties of promising materials. The focus of the present study is layered, Bi‐rich bismuth oxyhalides, which are widely studied for photocatalytic applications yet poorly structurally understood, due to high levels of disorder, nano‐sized domains, and the large number of structurally similar compounds. By connecting insights from multiple scattering techniques, utilizing electron‐, X‐ray‐ and neutron probes, the crystal phase of the synthesized materials is allocated as layered Bi 24 O 31 X 10 (X = Cl, Br), albeit with significant deviation from the reported 3D crystalline model. The materials comprise anisotropic platelet‐shaped crystalline domains, exhibiting significant in‐plane ordering in two dimensions but disorder and an ultra‐thin morphology in the layer stacking direction. Increased synthesis pH tailored larger, more ordered crystalline domains, leading to longer excited state lifetimes determined via femtosecond transient absorption spectroscopy (fs‐TAS). Although this likely contributes to improved photocatalytic properties, assessed via the photooxidation of benzylamine, increasing the overall surface area facilitated the most significant improvement in photocatalytic performance. This study, therefore, enabled both phase allocation and a nuanced discussion of the structure‐property relationship for complicated, ultra‐thin photocatalysts.

Chemistry↗
DOE OSTI2020

Interlaboratory Study of Digital Volume Correlation Error Due to X-Ray Computed Tomography Equipment and Scan Parameters: an Update from the DVC Challenge

Background: The quality of Digital Volume Correlation (DVC) full-field displacement measurements depends directly on the characteristics of the X-ray Computed Tomography (XCT) equipment, and scan procedures used to acquire the tomographic images. Objective: In this work, we seek to experimentally study the effects of XCT equipment and tomographic scan procedures on the quality of these images for DVC analysis, and to survey the level of DVC error that may be achieved using standard XCT operating procedures. Methods: Six participants in an interlaboratory study acquired high-quality XCT scans of a syntactic foam before and after rigid body motion. The resulting images were correlated using commercial DVC software to quantify error sources due to random image noise, reconstruction artifacts, as well as systematic spatial or temporal distortion. Results: In the absence of rigid body motion, the standard deviation of the displacement measurements ranged from 0.012 to 0.043 voxels using a moderate subvolume size, indicating that subvoxel measurement resolution could readily be achieved with a variety of XCT equipment and scan recipes. Comparison of consecutive scans without rigid body motion showed transient dilatational displacement gradients due to self-heating of the X-ray source and/or thermal expansion of the foam. Evaluation of the scans after rigid body motion showed significant, machine-specific spatial distortion in the displacement fields of up to 0.5 voxels; new approaches to remove this error need to be developed. Conclusions: Analysis of the scan protocols used in the interlaboratory study, as well as a complementary parametric sensitivity study, showed that the DVC error was strongly influenced by the XCT equipment, but could be mitigated by adjusting the total scan duration.

47 OTHER INSTRUMENTATION↗
DOE OSTI2020

K-shell internal conversion coefficient for M4 decay of the 30.8 keV isomer in 93 Nb

The 30.8 keV (1/2) – isomeric state in $^{93}_{41}$Nb decays by an M4 transition to the (9/2) + ground state. Here, we have measured the K-shell internal conversion coefficient α K for this transition in order to test the validity of current methods of calculation for a case with relatively low atomic number, low transition energy, and a high value of α K . Taking the fluorescence yield of niobium to be 0.751(4), we obtain α K = 2.56(9) × 10 4 , a result that agrees with the Dirac-Fock calculations that take into account the presence of a K vacancy in the final state, and disagrees by 1.7 standard deviations with calculations that ignore the vacancy. We also determine the energy of the isomeric state to be 30.760(5) keV, a fourfold improvement in precision over the currently accepted value.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Comparing pion production in transport simulations of heavy-ion collisions at 270⁢𝐴 MeV under controlled conditions

Within the TMEP, we present a detailed study of the performance of different transport models in Sn+Sn collisions at 270A MeV, and put particular emphasis on the production of pions and Delta resonances, which have been used as probes of the nuclear symmetry energy. We prescribe a common and rather simple physics model, and follow in detail the results of 4 BUU models and 6 QMD models. The nucleonic evolution of the collision and the nucleonic observables in these codes do not completely converge, but the differences among the codes can be understood as being due to several reasons: the basic differences between BUU and QMD models in the representation of the phase-space distributions, computational differences in the mean-field evaluation, and differences in the adopted strategies for the Pauli blocking in the collision integrals. For pionic observables, we find that a higher maximum density leads to an enhanced pion yield and a reduced pi- pi+ yield ratio, while a more effective Pauli blocking generally leads to a slightly suppressed pion yield and an enhanced pi- pi+ yield ratio. We specifically investigate the effect of the Coulomb force, and find that it increases the total pi- pi+ yield ratio but reduces the ratio at high pion energies, although differences in its implementations do not have a dominating role in the differences among the codes. Taking into account only the results of codes that strictly follow the homework specifications, we find a convergence of the codes in the final charged pion yield ratio to a 1 sigma deviation of about 5%. Furthermore, the uncertainty is expected to be reduced to about 1.6% if the same or similar strategies and ingredients, i.e., an improved Pauli blocking and calculation of the non-linear term in the mean-field potential, are similarly used in all codes.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Dual-Band Polarization Control with Pairwise Positioning of Polarization Singularities in Metasurfaces

Emerging applications of metasurfaces in classical and quantum optics are driving the need for precise polarization control of nearly-degenerate, high quality (𝑄)-factor modes. However, current approaches to creating specifically polarized pairs of modes force a trade-off between maintaining high 𝑄 factors and robustness. Here, we solve this challenge by employing pairwise generation, annihilation, and positioning of polarization singularities, derived from symmetry-guaranteed pairs of symmetry-protected bound states in the continuum. We experimentally demonstrate this design paradigm in silicon metasurfaces with mode splittings of ≈ 20 nm, mode splitting deviations as low as 1 nm, and 𝑄 factors up to 200. Furthermore, this approach opens new avenues for enhancing metasurface performance across a diverse range of applications, including sensing, modulating, nonlinear mixing, and generating quantum light.

Metasurfaces↗
DOE OSTI2020

Ultrarelativistic quark-nucleus scattering in a light-front Hamiltonian approach

In this work, we investigate the scattering of a quark on a heavy nucleus at high energies using the time-dependent basis light-front quantization (tBLFQ) formalism, which is the first application of the tBLFQ formalism in QCD. We present the real-time evolution of the quark wave function in a strong classical color field of the relativistic nucleus, described as the color glass condensate. The quark and the nucleus color field are simulated in the QCD SU(3) color space. We calculate the total and the differential cross sections, and the quark distribution in coordinate and color spaces using the tBLFQ approach. We recover the eikonal cross sections in the eikonal limit. We find that the differential cross section from the tBLFQ simulation is in agreement with a perturbative calculation at large $p_⊥$, and it deviates from the perturbative calculation at small $p_⊥$ due to higher-order contributions. In particular, we relax the eikonal limit by letting the quark carry realistic finite longitudinal momenta. We study the sub-eikonal effect on the quark through the transverse coordinate distribution of the quark with different longitudinal momentum, and we find the sub-eikonal effect to be sizable. Our results can significantly reduce the theoretical uncertainties in small $p_⊥$ region which has important implications to the phenomenology of the hadron-nucleus and deep inelastic scattering at high energies.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Simulated production of OH, HO 2 , CH 2 O, and CO 2 during dilute fuel oxidation can predict 1st-stage ignition delays

Chemical kinetics simulations are used to explore whether detailed measurements of relevant chemical species during the oxidation of very dilute fuels (less than 1 Torr partial pressure) in a high-pressure plug flow reactor (PFR) can predict autoignition propensity. We find that for many fuels the timescale for the onset of spontaneous oxidation in dilute fuel/air mixtures in a simple PFR is similar to the 1st-stage ignition delay time (IDT) at stoichiometric engine-relevant conditions. For those fuels that deviate from this simple trend, the deviation is closely related to the peak rate of production of OH, HO 2 , CH 2 O, and CO 2 formed during oxidation. We use these insights to show that an accurate correlation between simulated profiles of these species in a PFR and 1st-stage IDT can be developed using convolutional neural networks. Our simulations suggest that the accuracy of such a correlation is 10–50%, which is appropriate for rapid fuel screening and may be sufficient for predictive fuel performance modeling.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗