Search NASA⌕ Search

SEARCH · Search NASA

Results for “data and reaction mechanisms”

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 91 records · Page 5

Leaching Kinetics of Rare Earth Elements from Fire Clay Seam Coal

Recovery of rare earth elements (REEs) from coal samples collected from the Fire Clay coal seam using diluted mineral acid solutions was investigated. The initial processing step was coal recovery using conventional froth flotation which concentrated the REEs in tailing material resulting in an upgrade to values around 700 ppm on a dry whole mass basis. Leaching experiments were performed on the flotation tailings material using a 1.2 M sulfuric acid solution adjusted to a temperature of 75 °C to study the extractability of REEs from coal material. The effect of particle size, leaching time, leaching temperature, and solid concentration on REE leaching recovery were evaluated. The kinetic data obtained from leaching over a range of temperatures suggested that the leaching process follows the shrinking core model with possibly a mixed control mechanism that may be a result of several heterogenous materials leaching simultaneously. Leaching recovery increased rapidly at the beginning of the reaction then slowed as the system reached equilibrium. The apparent activation energy determined from test data obtained over a range of temperatures using 1 M sulfuric acid was 36 kJ/mol for the first 20 min of reaction time and 27 kJ/mol for the leaching period between 20 and 120 min. The leaching of light REEs during the initial stage was determined to be driven by a chemical reaction, followed by the formation of a product layer, which required lower activation energy in the later stage of leaching. In regards to the heavy REEs, the major mechanism for leaching is desorption and the product layer formation does not affect the heavy REEs significantly.

58 GEOSCIENCES↗

Unconventional gas-phase preparation of the prototype polycyclic aromatic hydrocarbon naphthalene (C 10 H 8 ) via the reaction of benzyl (C 7 H 7 ) and propargyl (C 3 H 3 ) radicals coupled with hydrogen-atom assisted isomerization

Polycyclic aromatic hydrocarbons (PAHs) are ubiquitous in the interstellar medium and in meteorites such as Murchison and Allende and signify the missing link between resonantly stabilized free radicals and carbonaceous nanoparticles (soot particles, interstellar grains). However, the predicted lifetime of interstellar PAHs of some 108 years imply that PAHs should not exist in extraterrestrial environments suggesting that key mechanisms of their formation are elusive. Exploiting a microchemical reactor and coupling these data with computational fluid dynamics (CFD) simulations and kinetic modeling, we reveal through an isomer selective product detection that the reaction of the resonantly stabilized benzyl (C 7 H 7 ) and the propargyl (C 3 H 3 ) synthesizes the simplest representative of PAHs – the 10π Hückel aromatic naphthalene (C 10 H 8 ) molecule – via the novel Propargyl Addition–BenzAnnulation (PABA) mechanism. The gas-phase preparation of naphthalene affords a versatile concept of the reaction of combustion and astronomically abundant propargyl radicals with aromatic radicals carrying the radical center at the methylene moiety (aromatic-CH 2 ) as a previously passed over source of aromatics in high temperature environments thus bringing us closer to an understanding of the aromatic universe we live in.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Kinetic study of the azo – Hydrazone photoinduced mechanism in complexes with a -Re(CO) 3 L 0/+ core by flash photolysis

The photoprocesses of three rhenium(I) tricarbonyl complexes, fac-Re(NN)(CO) 3 L 0/+ , bearing an azo/hydrazone ligand, NN, and L = Cl, Br or CH 3 CN as axial ligand, were studied by time-resolved spectroscopic techniques. Flash irradiation of deaerated solutions of the complexes in acetonitrile photogenerated intermediates with intense optical absorptions in the UV–vis spectral region of the hydrazone ligand. The decay of the intermediates photogenerated with the Cl, Br and CH 3 CN ligands in the axial position proceeds by three different mechanisms. For the fac-Re(NN)(CO) 3 Cl, a three-step mechanism involving isomerization, tautomerization and rotation accounts for the spectroscopic transformations. In contrast to the flash photolysis result, the continuous photolysis of the complex in acetonitrile does not show spectroscopic changes associated with the generation of a stable product, i.e., the complex exhibits a photoreversible behavior. The first step observed in the flash irradiation of fac-Re(NN)(CO) 3 Cl, is missing in the decay of the photogenerated intermediates with fac-Re(NN)(CO) 3 Br and fac-[Re(NN)(CO) 3 CH 3 CN] + . A possible reason for this absence is that the process is too fast for the time-response of the flash photolysis set-up. Time-resolved spectra and associated reaction rate data for the Re(I) complexes are communicated and discussed. Finally, DFT calculations support a plausible E → Z mechanism with thermal back-isomerization mediated by tautomerization.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Simplifying the creation of thermal decomposition mechanisms for designing mixed fibre composites

Among the most challenging aspects of simulating thermal decomposition of fibre reinforced polymers is the determination of appropriate reaction parameters. Thermogravimetric analysis is typically used to generate decomposition data, to which the reaction parameters are then fit. When designing mixed fibre materials (e.g. combinations of glass and carbon fibres), the number of TGA experiments needed to explore the entire design space may be intractable. Here, we demonstrate the creation of two candidate proxy mixed fibre mechanisms and compare them to a mechanism created by fitting parameters from TGA on the mixed fibre composite. These mechanisms are then demonstrated in a 2D axisymmetric numerical decomposition, heat transfer, and porous flow model. We find a maximum 11% uncertainty in mass and 4% in temperature difference when using a proxy mechanism.

Scott, Sarah N. [Sandia National Lab. (SNL-CA), Li↗

Understanding particle size effect on fast-charging behavior of graphite anode using ultra-thin-layer electrodes

Extreme fast charging (≤ 15 min) of lithium-ion batteries is highly desirable to accelerate mass-market adoption of electric vehicles. However, significant capacity fading, as well as safety issues due to the lithium plating caused by the fast charging rate, limit its implementation. In this study, we investigated the fast-charging capability of graphite materials with various particle sizes. To eliminate the Li + ion concentration gradient effect across the thickness of the electrode, ultra-thin-layer graphite electrodes were developed to investigate the "real" fast-charging capability of graphite at the particle level. Electrochemical assessments as well as microscopic characterizations revealed that smaller particles exhibited superior fast-charging performance, featuring enhanced capacity reversibility, faster charging rate, and less lithium plating under the same fast-charging conditions. It is shown that small-particle graphite (mean radius of 3.3 μm) could withstand a 4C charge (to 80 % state-of-charge) without plating, with minimal plating occurring at 6C. Thicker particles exhibited plating at lower C-rates. Since the experimental data could not directly explain whether intra-particle diffusion limitations or interfacial reaction limitations dominated the plating mechanism, the pseudo-2-dimensional model was used to evaluate the most likely plating mechanism. The model suggested that particle-level diffusion is the dominant mechanism contributing to plating at high rates. Finally, this work provides comprehensive insights into the particle size effects on fast-charging capability, offering a better understanding of fast-charging behavior and valuable guidance for designing optimal electrode architecture for high-rate lithium-ion batteries.

25 ENERGY STORAGE↗

Nanoconfinement of Molecular Magnesium Borohydride Captured in a Bipyridine-Functionalized Metal-Organic Framework

The lower limit of metal hydride nanoconfinement is demonstrated through the coordination of a molecular hydride species to binding sites inside the pores of a metal-organic framework (MOF). Magnesium borohydride, which has a high hydrogen capacity, is incorporated into the pores of UiO-67bpy (Zr6O4(OH)4(bpydc)6 with bpydc2- = 2,2'-bipyridine-5,5'- dicarboxylate) by solvent impregnation. The MOF retained its long-range order, and transmission electron microscopy and elemental mapping confirmed the retention of the crystal morphology and revealed a homogeneous distribution of the hydride within the MOF host. Notably, the B-, N-, and Mg-edge XAS data confirm the coordination of Mg(II) to the N atoms of the chelating bipyridine groups. In situ 11B MAS NMR studies helped elucidate the reaction mechanism and revealed that complete hydrogen release from Mg(BH4)2 occurs as low as 200 °C. Sieverts and thermogravimetric measurements indicate an increase in the rate of hydrogen release, with the onset of hydrogen desorption as low as 120 °C, which is approximately 150 °C lower than that of the bulk material. Furthermore, density functional theory calculations support the improved dehydrogenation properties and confirm the drastically lower activation energy for B-H bond dissociation.

Nanoconfinement, Metal Hydrides, Metal-Organic Fra↗

Oxidation and pyrolysis of methyl propyl ether

The ignition, oxidation, and pyrolysis chemistry of methyl propyl ether (MPE) was probed experimentally at several different conditions, and a comprehensive chemical kinetic model was constructed to help understand the observations, with many of the key parameters computed using quantum chemistry and transition state theory. Experiments were carried out in a shock tube measuring time variation of CO concentrations, in a flow tube measuring product concentrations, and in a rapid compression machine (RCM) measuring ignition delay times. The detailed reaction mechanism was constructed using the Reaction Mechanism Generator software. Sensitivity and flux analyses were used to identify key rate and thermochemical parameters, which were then computed using quantum chemistry to improve the mechanism. Validation of the final model against the 1-20 bar 600-1500 K experimental data is presented with a discussion of the kinetics. The model is in excellent agreement with most of the shock tube and RCM data. Strong non-monotonic variation in conversion and product distribution is observed in the flow-tube experiments as the temperature is increased, and unusually strong pressure dependence and significant heat release during the compression stroke is observed in the RCM experiments. These observations are largely explained by a close competition between radical decomposition and addition to O2 at different sites in MPE; this causes small shifts in conditions to lead to big shifts in the dominant reaction pathways. The validated mechanism was used to study the chemistry occurring during ignition in a diesel engine, simulated using Ignition Quality Test (IQT) conditions. At the IQT conditions, where the MPE concentration is higher, bimolecular reactions of peroxy radicals are much more important than in the RCM.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A Coverage Self-Consistent Microkinetic Model for Vapor-Phase Formic Acid Decomposition over Pd/C Catalysts

An iterative approach utilizing density functional theory (DFT, PW91-GGA)-informed mean-field microkinetic models and reaction kinetics experiments is used to determine the reaction mechanism and the active site for formic acid (HCOOH, FA) decomposition over a Pd/C catalyst. Models parametrized using DFT energetics on clean Pd(100) and Pd(111) required large corrections to the DFT energetics for capturing our experimental data. Further, both Pd(111) and Pd(100) models predicted a high coverage of adsorbed CO (CO*), inconsistent with the assumption of a clean surface at which the rate parameters for these models were calculated. To better represent the active site under reaction conditions and explicitly account for the presence of CO*, subsequent microkinetic models were formulated using DFT energetics that were calculated on partially (5/9 ML) CO*-covered Pd (111) and (100) facets. Upon parameter adjustment, the resultant 5/9 ML CO*-covered Pd(100) model, although consistent in terms of CO* coverage, was unable to capture the dehydration path measured in the experiments and was, therefore, deemed not to offer an accurate representation of the active site for FA decomposition over Pd/C. In contrast, a partially CO*-covered Pd(111) model was better at representing the catalytic active site, as in addition to being consistent in terms of CO* coverages, it required small adjustments of the DFT parameters to accurately capture the experimental data set (both dehydrogenation and dehydration). Our results suggest that the reaction occurs via the spectroscopically elusive carboxyl (COOH*) intermediate and that spectator CO*-assisted decomposition pathways play an important role under typical experimental conditions. In addition, our study highlights the importance of striving for coverage self-consistent microkinetic models and for including spectator-assisted mechanisms in order to develop an improved picture of the active site under reaction conditions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Structural Characterization of the Reaction and Substrate Specificity Mechanisms of Pathogenic Fungal Acetyl-CoA Synthetases

Acetyl CoA synthetases (ACSs) are Acyl-CoA/NRPS/Luciferase (ANL) superfamily enzymes that couple acetate with CoA to generate acetyl CoA, a key component of central carbon metabolism in eukaryotes and prokaryotes. Normal mammalian cells are not dependent on ACSs, while tumor cells, fungi, and parasites rely on acetate as a precursor for acetyl CoA. Consequently, ACSs have emerged as a potential drug target. As part of a program to develop antifungal ACS inhibitors, we characterized fungal ACSs from five diverse human fungal pathogens using biochemical and structural studies. ACSs catalyze a two-step reaction involving adenylation of acetate followed by thioesterification with CoA. Our structural studies captured each step of these two half-reactions including the acetyl-adenylate intermediate of the first half-reaction in both the adenylation conformation and the thioesterification conformation and thus provide a detailed picture of the reaction mechanism. We also used a systematic series of increasingly larger alkyl adenosine esters as chemical probes to characterize the structural basis of the exquisite ACS specificity for acetate over larger carboxylic acid substrates. Consistent with previous biochemical and genetic data for other enzymes, structures of fungal ACSs with these probes bound show that a key tryptophan residue limits the size of the alkyl binding site and forces larger alkyl chains to adopt high energy conformers, disfavoring their efficient binding. Together, our analysis provides highly detailed structural models for both the reaction mechanism and substrate specificity that should be useful in designing selective inhibitors of eukaryotic ACSs as potential anticancer, antifungal, and antiparasitic drugs.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Exploring the $A$-V-Sb landscape beyond $A$$\mathrm{V_3}$ $\mathrm{Sb_5}$: A case study on the $\mathrm{KV_6}$ $\mathrm{Sb_6}$ Kagome compound

Kagome compounds have garnered attention in the past few years for their intriguing magnetic properties arising from spin frustration dictated by the geometry of the Kagome sublattice. In this paper, we highlight the success of the unconventional hydride route for the fast and easy synthesis of the Kagome compound $\mathrm{KV_6}$ $\mathrm{Sb_6}$. Here, high-temperature in situ powder x-ray diffraction (PXRD) studies proved to be useful in hinting at the existence of $\mathrm{KV_6}$ $\mathrm{Sb_6}$, identifying its synthesis conditions, and understanding the reaction mechanism. The crystal structure for $\mathrm{KV_6}$ $\mathrm{Sb_6}$ was determined from high-resolution PXRD data. The compound has a layered structure [$R\overline{3}m, a = 5.5318$ (9) Å, $c = 34.23$ (3)Å 3 $Z = 3$] at room temperature] and features a Kagome bilayer of V atoms. $\mathrm{KV_6}$ $\mathrm{Sb_6}$ is isostructural to the previously reported $\mathrm{RbV_6}$ $\mathrm{Sb_6}$ and $\mathrm{CsV_6}$ $\mathrm{Sb_6}$ compounds. $\mathrm{KV_6}$ $\mathrm{Sb_6}$ is thermally stable in vacuum up to 1173 K, as evident from the high-temperature in situ PXRD and differential scanning calorimetric analysis. Investigation of magnetic properties for $\mathrm{KV_6}$ $\mathrm{Sb_6}$ between 2 and 300 K reveals temperature-independent paramagnetism and an absence of superconductivity, like the Rb and Cs analogs. Furthermore, we compare the magnetic properties of K V 3 Sb 5 , another ternary Kagome compound, synthesized via two different methods: the hydride route and the traditional route from elements. Low-temperature transport property measurements of $\mathrm{KV_6}$ $\mathrm{Sb_6}$ indicate metallic behavior and an intrinsically low thermal conductivity of $\mathrm{1.0 W K^-1}$ $\mathrm{m^-1}$ at 300 K. The layered structure of $\mathrm{KV_6}$ $\mathrm{Sb_6}$ makes it an attractive candidate for deintercalation and doping studies to tune both magnetic and transport properties, laying a foundation for further studies.

36 MATERIALS SCIENCE↗

An Assessment of Kinetic Models for Ammonia Flame Extinction

The growing interest in using ammonia as a carbon-neutral replacement fuel has prompted numerous recent research efforts towards understanding fundamental combustion characteristics and chemical kinetics of ammonia oxidation. Despite these recent efforts, there remains a large deficiency in extinction limit data for ammonia flames, particularly under elevated pressure and temperature conditions relevant to gas turbines. Even the literature reaction mechanisms developed specifically for ammonia have been untested under gas turbine relevant conditions, and existing kinetic models have shown large variations in laminar flame speed, extinction strain rate, and speciation predictions for NH3–H2 mixtures at ambient conditions. Recent work by Thomas and co-workers, for example, compared different ammonia kinetic models in the literature for predicting their extinction strain rate measurements for non-premixed NH3-H2 counterflow flames, and again have highlighted the discrepancies in model predictions. Here, our work has undertaken to further highlight important reaction pathways for ammonia oxidation that may be responsible for these discrepancies for extinction strain rate predictions, as well as elucidate sensitive reaction steps with large uncertainties that require further attention. This work points to a need for an improved understanding of ammonia combustion chemistry as well as new datasets at elevated pressure and temperature conditions to better constrain rate expressions for key reaction steps.

ammonia, flame, combustion, kinetic, kinetics↗

Impact of Two-Body Currents on Semi-Exclusive Lepton-Nucleus Reactions

We generalize the spectral-function formalism to describe two-nucleon knockout processes in exclusive kinematics. Significant improvements are introduced both in the treatment of the current operators entering the $Δ$-current contribution and in the modeling of correlations between the two struck nucleons, including a consistent treatment of isospin dependence and the explicit incorporation of angular correlations. The framework is validated through comparisons with relativistic Fermi-gas calculations and with semi-exclusive electron-nucleus scattering data. Our results demonstrate that an accurate description of nuclear dynamics plays a crucial role in modeling this reaction mechanism. We further present a study of selected electroweak observables relevant to neutrino-scattering experiments.

Rocco, N. [Fermilab; Valencia U., IFIC] (ORCID:000↗

Impact of Two-Body Currents on Semi-Exclusive Lepton-Nucleus Reactions

We generalize the spectral-function formalism to describe two-nucleon knockout processes in exclusive kinematics. Significant improvements are introduced both in the treatment of the current operators entering the $Δ$-current contribution and in the modeling of correlations between the two struck nucleons, including a consistent treatment of isospin dependence and the explicit incorporation of angular correlations. The framework is validated through comparisons with relativistic Fermi-gas calculations and with semi-exclusive electron-nucleus scattering data. Our results demonstrate that an accurate description of nuclear dynamics plays a crucial role in modeling this reaction mechanism. We further present a study of selected electroweak observables relevant to neutrino-scattering experiments.

Rocco, N. [Fermilab; Valencia U., IFIC] (ORCID:000↗

Impact of Two-Body Currents on Semi-Exclusive Lepton-Nucleus Reactions

We generalize the spectral-function formalism to describe two-nucleon knockout processes in exclusive kinematics. Significant improvements are introduced both in the treatment of the current operators entering the $Δ$-current contribution and in the modeling of correlations between the two struck nucleons, including a consistent treatment of isospin dependence and the explicit incorporation of angular correlations. The framework is validated through comparisons with relativistic Fermi-gas calculations and with semi-exclusive electron-nucleus scattering data. Our results demonstrate that an accurate description of nuclear dynamics plays a crucial role in modeling this reaction mechanism. We further present a study of selected electroweak observables relevant to neutrino-scattering experiments.

Rocco, N. [Fermilab; Valencia U., IFIC] (ORCID:000↗

Impact of Two-Body Currents on Semi-Exclusive Lepton-Nucleus Reactions

We generalize the spectral-function formalism to describe two-nucleon knockout processes in exclusive kinematics. Significant improvements are introduced both in the treatment of the current operators entering the $Δ$-current contribution and in the modeling of correlations between the two struck nucleons, including a consistent treatment of isospin dependence and the explicit incorporation of angular correlations. The framework is validated through comparisons with relativistic Fermi-gas calculations and with semi-exclusive electron-nucleus scattering data. Our results demonstrate that an accurate description of nuclear dynamics plays a crucial role in modeling this reaction mechanism. We further present a study of selected electroweak observables relevant to neutrino-scattering experiments.

Rocco, N. [Fermilab; Valencia U., IFIC] (ORCID:000↗

Mechanism of Catechol Oxidation by the Coupled Bicupric Active Site of Tyrosinase: Completion of the Oxygenase/Oxidase Reaction Cycle

Tyrosinase contains a coupled binuclear copper (CBC) active site, which in its bicuprous form (deoxy-Ty) binds O 2 to form a side-on peroxide [Cu(II) 2 O 2 ] 2+ intermediate (oxy-Ty) that performs the regioselective monooxygenation of monophenols to catechols and their subsequent 2e − oxidation to quinones. Previously, we used spectroscopic, kinetic, and computational methods to elucidate the mechanism of the initial steps in the monooxygenation reaction. Herein, we investigated the final step in catecholate formation and provide experimental and computational results elucidating the catechol oxidation reaction, formally a proton-coupled 2e − process. We employed single-turnover stoppedflow absorption to observe the elusive catecholate-Cu(II) 2 OH intermediate via the anaerobic reaction of resting Cu(II) 2 OH Ty (met-Ty) with the catecholate substrate. This intermediate was cryo-trapped and characterized by electron paramagnetic resonance and X-ray absorption spectroscopies. These experimental results were correlated to quantum mechanics/molecular mechanics (QM/MM) and QM calculations to describe the complete catalytic cycle of tyrosinase, revealing that the 2e − oxidation of catechol to quinone involves two steps: (i) coupled 1H + /1e − transfer from the bound monoanionic catecholate to the bridged hydroxide of the met-Ty active site, followed by (ii) the second 1e − transfer. The latter is calculated to be the rate-limiting step in catechol oxidation, confirmed by experimental solvent kinetic isotope effect studies. Our computational data suggest that quinone release from the protein provides the necessary driving force for this reaction. These results establish a detailed description of this oxidase cycle of the tyrosinase mechanism and more broadly provide molecular-level insights into the diverse reactivity of CBC sites in biology.

QM/MM modeling↗

Neutron Insights into Sorption Enhanced Methanol Catalysis

Sorption enhanced methanol production makes use of the equilibrium shift of the CO 2 hydrogenation reaction towards the desired products. However, the increased complexity of the catalyst system leads to additional reactions and thus side products such as dimethyl ether, and complicates the analysis of the reaction mechanism. On the other hand, the unusually high concentration of intermediates and products in the sorbent facilitates the use of inelastic neutron scattering (INS) spectroscopy. Despite being a post-mortem method, the INS data revealed the change of the reaction path during sorption catalysis. Concretely, the experiments indicate that the varying water partial pressure due to the adsorption saturation of the zeolite sorbent influences the progress of the reaction steps in which water is involved. Experiments with model catalysts support the INS findings.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Effects of stereoisomeric structure and bond location on the ignition and reaction pathways of hexenes

The current work presents new experimental autoignition and speciation data on the two cis-hexene isomers: cis-2-hexene and cis-3-hexene. The new data provide insights on the effects of carbon-carbon double bond location and stereoisomeric structures on ignition delay times and reaction pathways for linear hexene isomers. Experiments were performed using the University of Michigan rapid compression facility to determine ignition delay times from pressure-time histories. Stoichiometric (Φ = 1.0) mixtures at dilution levels of inert gas to O 2 = 7.5:1 (mole basis) were investigated at an average pressure of 11 atm and temperatures from 809 to 1052 K. Speciation experiments were conducted at T = 900 K for the two cis-hexene isomers, where fast-gas sampling and gas chromatography were used to identify and quantify the two cis-hexene isomers and stable intermediate species. The ignition delay time data showed negligible sensitivity to the location of the carbon-carbon double bond and the stereoisomeric structure (cis-trans), and the species data showed no correlation with the stereoisomeric structure, but there was a strong correlation of some of the measured species with the location of the double bond in the hexene isomer. In particular, 2-hexene showed strong selectivity to propene, acetaldehyde, and 1,3-butadiene, and 3-hexene showed selectivity to propanal. Model predictions of ignition delay times were in excellent agreement with the experimental data. Furthermore, there was generally good agreement for the model predictions of the species data for 2-hexene; however, the mechanism overpredicted some of the small aldehyde (C 2 -C 4 ) species for 3-hexene. Reaction pathway analysis indicates the hexenes are almost exclusively consumed by H-atom abstraction reactions at the conditions studied (P = 11 atm, T > 900 K), and not by C 3 -C 4 scission as observed in high-temperature (>1300 K) hexene ignition studies. Improved estimates for 3-hexene + OH reactions may improve model predictions for the species measured in this work.

42 ENGINEERING↗