Search NASA⌕ Search

SEARCH · Search NASA

Results for “vibrational frequencies”

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
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

The Remarkable I 2 O 3 Molecule: A New View from Theory

Atmospheric iodine chemistry has garnered increasing attention as a result of increased iodine emissions. A key subset of this chemistry involves iodine oxides (I 2 O 2–5 ), which serve as precursors to particle formation. Among these, I 2 O 3 is the simplest iodine oxide involved in particle formation, but it has remained undetected in the atmosphere. Previous theoretical studies have characterized this peculiar molecule, primarily using energies to refine geometries obtained at low levels of theory. Due to the reemerging interest in I 2 O 3 , this study presents geometries optimized at the CCSD(T)/aug-cc-pwCVTZ-PP level of theory─marking the first instance, to the best of our knowledge, where this system has been studied exclusively with CCSD(T). Harmonic vibrational frequencies were computed at the same level of theory. Final energetics were obtained using the very high level CCSDT(Q) method with basis sets up to quintuple-zeta cardinality (aug-cc-pwCV5Z-PP) and extrapolated to the CBS limit to yield CCSDT(Q)/CBS//CCSD(T)/aug-cc-pwCVTZ-PP energies. These energies include harmonic zero-point vibrational energy corrections and scalar relativistic energy corrections. Additionally, this study discovers new isomers along the I 2 O 3 potential energy surface, a novel contribution to the field. The performance of different computational methods and DFT functionals commonly used in atmospheric chemistry is also assessed relative to high-level theoretical methods.

basis sets↗
DOE OSTI2023

Helical Organic and Inorganic Polymers

Despite being a staple of synthetic plastics and biomolecules, helical polymers are scarcely studied with Gaussian-basis-set ab initio electron-correlated methods on an equal footing with molecules. This article introduces an ab initio second-order many-body Green’s-function [MBGF(2)] method with nondiagonal, frequency-dependent Dyson self-energy for infinite helical polymers using screw-axis-symmetry-adapted Gaussian-spherical-harmonics basis functions. Together with the Gaussian-basis-set density-functional theory for energies, analytical atomic forces, translational-period force, and helical-angle force, it can compute correlated energy, quasiparticle energy bands, structures, and vibrational frequencies of an infinite helical polymer, which smoothly converge at the corresponding oligomer results. These methods can handle incommensurable structures, which have an infinite translational period and are hard to characterize by any other method, just as efficiently as commensurable structures. We apply them to polyethylene (2/1 helix), polyacetylene (Peierls’ system), and polytetrafluoroethylene (13/6 helix) to establish the quantitative accuracy of MBGF(2)/cc-pVDZ in simulating their (angle-resolved) ultraviolet photoelectron spectra, and of B3LYP/cc-pVDZ or 6-31G** in reproducing their structures, infrared and Raman band positions, phonon dispersions, and (coherent and incoherent) inelastic neutron scattering spectra. We then predict the same properties for infinitely catenated chains of nitrogen or oxygen and discuss their possible metastable existence under ambient conditions. They include planar zigzag polyazene (N2)x (Peierls’ system), 11/3-helical isotactic polyazane (NH)x, 9/4-helical isotactic polyfluoroazane (NF)x, and 7/2-helical polyoxane (O)x as potential high-energy-density materials.

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

Observation of bound valence excited electronic states of deprotonated 2-hydroxytriphenylene using photoelectron, photodetachment, and resonant two-photon detachment spectroscopy of cryogenically cooled anions

Polycyclic aromatic hydrocarbons (PAHs) are common atmospheric pollutants, and they are also ubiquitous in the interstellar medium. Here, we report the study of a complex O-containing PAH anion, the deprotonated 2-hydroxytriphenylene (2-OtPh–), using high-resolution photoelectron imaging and photodetachment spectroscopy of cryogenically cooled anions. Vibrationally resolved photoelectron spectra yield the electron affinity of the 2-OtPh radical as 2.629(1) eV and several vibrational frequencies for its ground electronic state. Photodetachment spectroscopy reveals bound valence excited electronic states for the 2-OtPh– anion, with unprecedentedly rich vibronic features. Evidence is presented for a low-lying triplet state (T1) and two singlet states (S1 and S2) below the detachment threshold. Single-color resonant two-photon photoelectron spectroscopy uncovers rich photophysics for the 2-OtPh– anion, including vibrational relaxation in S1, internal conversion to the ground state of 2-OtPh–, intersystem crossing from S2 to T1, and a long-lived autodetaching shape resonance about 1.3 eV above the detachment threshold. The rich electronic structure and photophysics afforded by the current study suggest that 2-OtPh– would be an interesting system for pump–probe experiments to unravel the dynamics of the excited states of this complex PAH anion.

Chemistry↗
DOE OSTI2025

Hard–Soft Acid–Base Theory Explains Photoexcited Carrier Dynamics in Porphyrin/CNT Nanohybrids: Time-Domain Atomistic Analysis

We employ the fundamental chemical concepts of hard− soft acid−base to formulate general principles governing excited-state dynamics in zinc porphyrin (ZnP)/carbon nanotube (CNT) hybrids for energy photoconversion. Atomistic quantum dynamics simulations demonstrate that electron-withdrawing and donating substituents at the ZnP β-pyrrolic position strongly influence the dynamics. ZnP photoexcitation produces subpicosecond electron transfer (ET) from ZnP to CNT, in agreement with the experiment. Substitutions of CN by H and tBu accelerate the ET. The trend is directly related to the hard−soft acid−base concept because the soft−soft interaction between the t Bu- ZnP acid and the mild CNT base enhances the donor−acceptor coupling. Longer coherence and more active vibrational modes facilitate the ET in t Bu-ZnP/CNT. Electron−hole recombination in CN-ZnP/CNT occurs on a hundred picosecond time scale, nicely corroborated by the experiment. The exciton lifetime is extended beyond a nanosecond by the substitutions. The soft−soft interaction in t Bu-ZnP/ CNT increases the splitting between the highest occupied orbitals of the two subsystems, reduces their mixing, and decreases nonadiabatic coupling between the ground and excited states. Rapid decoherence and involvement of low-frequency vibrations favor longer lifetimes. Our investigation reveals that larger p K a of the β-pyrrolic acid gives rapid ET and slow recombination and provides detailed mechanistic information, essential for future optoelectronic applications.

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

Adsorption Properties of Au−Ni Surface Alloys with a Nonstoichiometric Moiré Structure: A Density Functional Theory Study

Due to the large lattice mismatch between gold and nickel, gold–nickel surface alloys can form unique nonstoichiometric overlayer structures characterized by a moiré pattern and subsurface defects. For this work, we performed density functional theory (DFT) calculations to study the adsorption of molecular oxygen, atomic hydrogen, and atomic carbon on a gold–nickel(111) surface alloy with 0.46 monolayer gold randomly distributed in the surface layer. We observed six distinct adsorption structures for molecular oxygen characterized by intramolecular stretching frequencies of <700, 729, 795, 857, 929, and 1004 cm –1 , which describe well the experimentally observed high-resolution electron energy-loss spectra. Surface atomic hydrogen adsorption is associated with adsorbate–surface modes in the ∼1000 cm –1 range, while subsurface hydrogen can have features as low as ∼400 cm –1 . We observed a unique adsorption structure for atomic carbon inside the surface dislocation loop defect, which explains the experimentally observed low carbon-surface mode at ∼340 cm –1 . Our study sheds light on the unique adsorption properties of the gold–nickel surface alloys and helps with rationalizing vibrational frequency experimental studies for this system.

adsorption↗
DOE OSTI2025

Vibrational Signatures of Electronic Properties in Renewable-Energy Catalysis

The objective of this research program is to discover and develop new approaches for ab initio computational simulations of molecular motion. Specifically, this program aims to decipher the connection between the underlying electronic structure and the accordant molecular vibrations in reactive ions and radicals for renewable-energy purposes. Recent experimental progress in ion sources and optical spectroscopies has unearthed considerable new inner-sphere detail for these complexes, but the connection between these spectral signatures and mechanistic information often remains elusive, to the continued frustration of experimentalists. For this purpose, new anharmonic vibrational frequency methods, along with a publicly deployed software package, will be developed. Working closely with committed experimental collaborators, this conceptual and computational framework will be used to explain the results of new spectroscopy experiments, focusing specifically on the inner-shell mechanisms of renewable-energy catalysis. The oxidation half of catalytic water-splitting chemistry will be a central focus, along with fundamental studies of the manner in which strong ions and radicals activate solvent as a chemical species. The resulting products of the research program will include openly available software and algorithms for the ab initio simulation of challenging vibrational spectra, as well as critical mechanistic insight into energy-focused catalytic processes that are opaque to other existing analytical techniques.

42 ENGINEERING↗
DOE OSTI2020

Reaction mechanism and kinetics for CO 2 reduction on nickel single atom catalysts from quantum mechanics

Experiments have shown that graphene-supported Ni-single atom catalysts (Ni-SACs) provide a promising strategy for the electrochemical reduction of CO 2 to CO, but the nature of the Ni sites (Ni-N 2 C 2 , Ni-N 3 C 1 , Ni-N 4 ) in Ni-SACs has not been determined experimentally. Here, we apply the recently developed grand canonical potential kinetics (GCP-K) formulation of quantum mechanics to predict the kinetics as a function of applied potential (U) to determine faradic efficiency, turn over frequency, and Tafel slope for CO and H 2 production for all three sites. We predict an onset potential (at 10 mA cm –2 ) U onset = –0.84 V (vs. RHE) for Ni-N 2 C 2 site and U onset = –0.92 V for Ni-N 3 C 1 site in agreement with experiments, and U onset = –1.03 V for Ni-N 4 . We predict that the highest current is for Ni-N 4 , leading to 700 mA cm –2 at U = –1.12 V. To help determine the actual sites in the experiments, we predict the XPS binding energy shift and CO vibrational frequency for each site.

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

Vibrational Spin-Orbit Coupling Contributions to Excited State Decay of Ligand-to-Ligand Charge Transfer States

Controlling excited state relaxation processes is important in a variety of photochemical and photophysical processes, including the generation of ground and excited state spin polarization for quantum information science applications. Here, we address how specific static distortions – based on vibrational spin-orbit active modes at C2v symmetry determined by group theory – in a series of low-symmetry ligand-to-ligand charge transfer complexes enable direct spin-orbit coupling contributions to T1 → S0 excited state decay. These results are used to address spin-vibronic coupling contributions to T1 → S0 decay in a high-symmetry (tBu2bpy)Pt(S,S) (tBu2bpy = 4,4’-di-tert-butyl-2,2’-bipyridine and S,S = benzene-1,2-dithiolate) ligand-to-ligand charge transfer complex with effective C2v symmetry, where T1 relaxation is both spin- and orbitally forbidden due to the direct spin-orbit coupling matrix element being zero by symmetry. Low-frequency vibrations that involve a pyridine-pyridine twisting motion within the bpy ligand generate large ∂/∂Qi values that will contribute significantly to T1 → S0 relaxation. The work advances ligand design strategies for the generation of tailored T1 → S0 relaxation rates, which can be utilized to optimize the generation of electron spin polarization in radical-elaborated ligand-to-ligand charge transfer complexes.

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

Convergent Concordant Mode Approach for Molecular Vibrations: CMA-2

The concordant mode approach (CMA) is a promising new scheme for dramatically increasing the system size and level of theory achievable in quantum chemical computations of molecular vibrational frequencies. Here, we achieve advances in the CMA hierarchy by computations targeting CCSD(T)/cc-pVTZ (coupled cluster singles and doubles with perturbative triples using a correlation-consistent polarized-valence triple-ζ basis set) benchmarks within the G2 molecular test set, executing a statistical analysis for 1501 frequencies from 111 compounds and then separately solving the refractory case of pyridine. First, MP2/cc-pVTZ (second-order Møller–Plesset perturbation theory with the same basis set) proves to be an excellent and preferred choice for generating the underlying (Level B) normal modes of the CMA scheme. Utilizing this Level B within the CMA-0A method reproduces the 1501 benchmark frequencies with a mean absolute error (MAE) of only 0.11 cm –1 and an attendant standard deviation of 0.49 cm –1 . Second, a convergent CMA-2 method is constituted that allows efficient computation of higher level (Level A) frequencies to any reasonable accuracy threshold by using only Hartree–Fock (HF) and MP2 or density functional theory (DFT) data to generate ξ parameters, which select the sparse off-diagonal force field elements for explicit evaluation at Level A. When Level B = MP2/cc-pVTZ, a cutoff of ξ = 0.02 provides an average maximum absolute error per molecule of only 0.17 cm –1 by incurring merely a 33% increase in average cost over CMA-0A. This CMA-2 method also eradicates the 4 problematic CMA-0A outliers of pyridine with even less effort (ξ = 0.04, 22% increase). Finally, the newly developed CMA procedures are shown to be highly successful when applied to 1-(1H-pyrrol-3-yl)ethanol, a new test molecule with diverse types of vibration.

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

High-resolution infrared spectroscopy of supersonically cooled singlet carbenes: Bromomethylene (HCBr) in the CH stretch region

First high-resolution spectra of cold (~35 K) singlet bromomethylene HCBr in the CH stretching (v 1 ) region from 2770 to 2850 cm-1 are reported using near quantum shot-noise limited laser absorption methods in a slit jet supersonic discharge expansion source. Three rovibrational bands are identified at high S/N (20:1 - 40:1) and rotationally assigned to i) the CH stretch fundamental (v 1 ) band $\tilde{X}$(1,0,0) ← $\tilde{X}$ (0,0,0) and ii) vibrational hot bands ($\tilde{X}$(1,1,0) ← $\tilde{X}$(0,1,0) and $\tilde{X}$(1,0,1) ← $\tilde{X}$(0,0,1)) arising from vibrationally excited HCBr populated in the discharge with single quanta in either the H–C–Br bend (v 2 ) or C–Br stretch (v 3 ) modes. Precision rotational constants are reported for a total of six states, with an experimentally determined CH stretch vibrational frequency (2799.38 cm -1 ) in good agreement with previous low-resolution fluorescence studies. Detailed analysis of the fundamental v 1 band highlights the presence of perturbations in the $\tilde{X}$(1,0,0) level, which we tentatively attribute to arise from the nearby triplet state $\tilde{a}$(0,0,1) through spin-orbit interaction or the multiple quanta $\tilde{X}$(0,2,1) singlet state via c-type Coriolis coupling. Reduced-Doppler resolution (60 MHz) in the slit-jet IR spectrometer permits clear observation of nuclear spin hyperfine structure, with experimental line shapes well reproduced by nuclear quadrupole/spin-rotation coupling constants from microwave studies. Lastly, the a-type to b-type transition intensity ratio for the fundamental CH stretch band is notably larger than predicted by a bond-dipole model, which from high level ab initio quantum calculations (CCSD(T)/PVQZ) can be attributed to vibrationally induced “charge-sloshing” of electron density along the polar C–Br bond.

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

Terahertz Circular Dichroism Spectroscopy of Molecular Assemblies and Nanostructures

Chemical, physical, biological and materials engineering disciplines use a variety of chiroptical spectroscopies to probe geometrical and optical asymmetry in molecules and particles. Electronic (ECD) and vibrational (VCD) circular dichroism are the most common of these techniques and collectively enable the studies of electronic and vibronic transitions with energies between 0.1 and 5.0 eV. The vibrational states with characteristic energies in the range of 0.001–0.01 eV carry valuable information about concerted intermolecular motions in molecules and crystals involving multiple atoms. These vibronic transitions located in the terahertz (THz) part of the spectrum become increasingly more important for the chemistry, physics, and biology of complex molecules and materials However, the methodology and hardware of THz circular dichroism (TCD) are much less developed than the chiroptical spectroscopies for ultraviolet, visible, near- and mid infrared photons. In this report we provide theoretical foundations, practical implementations, comparative assessments, and exemplary applications of TCD spectroscopy. We show that the sign, intensity, and position of TCD peaks are highly sensitive to the three-dimensional structure and long-range organization of molecular crystals, which offer unique capabilities to study (bio) molecules, their crystals, and nanoscale assemblies and apply the novel data processing methodologies. TCD also offers a convenient toolbox to identify new physical phenomena, such as chiral phonons and their propagation in nanostructured matter. We also discuss the major challenges, emerging opportunities and promising research directions, including broad investigation of chiral phonons observed in chiral (nano) crystals and emerging machine learning methodologies for TCD in biological and nanoscale structures. Ubiquity of low-frequency vibrations with rotational components in biomolecular structures, combined with sharpness of peaks in TCD spectra, enables a variety of technological translations.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Soft Lattice and Defect Covalency Rationalize Tolerance of β‐CsPbI 3 Perovskite Solar Cells to Native Defects

Abstract Although all‐inorganic metal halide perovskites (MHPs) have shown tremendous improvement, they are still inferior to the hybrid organic–inorganic MHPs in efficiency. Recently, a conceptually new β‐CsPbI 3 perovskite reached 18.4 % efficiency combined with good thermodynamic stability at ambient conditions. We use ab initio non‐adiabatic molecular dynamics to show that native point defects in β‐CsPbI 3 are generally benign for nonradiative charge recombination, regardless of whether they introduce shallow or deep trap states. These results indicate that MHPs do not follow the simple models used to explain defect‐mediated charge recombination in the conventional semiconductors. The strong tolerance is due to the softness of the perovskite lattice, which permits separation of electrons and holes upon defect formation, and only allows carriers to couple to the low‐frequency vibrations. Both factors decrease notably the non‐adiabatic coupling and slow down the dissipation of energy to heat.

Chu, Weibin↗
DOE OSTI2021

Rhodium Single–Atom Catalysts on Titania for Reverse Water Gas Shift Reaction Explored by First Principles Mechanistic Analysis and Compared to Nanoclusters

The thermocatalytic reduction of CO2 by H2 often proceeds via two competing reaction mechanisms—the reverse water gas shift reaction (rWGSR, CO2 + H2 ⇌ CO + H2O) and methanation (CO2 + 4H2 ⇌ CH4 + 2H2O). Atomically dispersed Rh1 catalysts on TiO2 show high selectivity toward the rWGSR compared with larger Rh nanoclusters, but the origin of this size-dependent selectivity has not been fully explained. Here we report density functional theory (DFT) calculations and microkinetic simulations that clarify the Rh1 active sites and rWGSR pathway on anatase TiO2(101), as well as the high rWGSR selectivity of Rh1 compared with supported Rhx (x = 2–8 atoms) nanoclusters. DFT-computed formation energies, vibrational frequency analysis, and microkinetic modeling suggest three plausible active sites: Rh1 on titania (Rh1/TiO2(101)), Rh1 with a nearby hydroxyl group (Rh1OH/TiO2(101)), and Rh1 near an oxygen vacancy at a three-fold coordinated site (Rh1 near O3cvac). Predicted turnover frequencies and apparent activation barriers for Rh1 indicate a faster reaction involving CO2 dissociation assisted by a support oxygen vacancy via Rh1 near O3cvac, as well as slower reactions involving Rh1OH/TiO2(101) or Rh1/TiO2(101) through a COOH intermediate. Furthermore, these Rh1 sites are selective toward CO rather than CH4 because of the weak adsorption of CO, large barrier for C-O bond dissociation, and the lack of nearby metal sites for H2 dissociation, in contrast to Rhx nanoclusters, including Rh2 dimers.

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

Insights into the Mechanism of Neptunium Oxidation to the Heptavalent State

Abstract Neptunium can exist in multiple oxidation states, including the rare and poorly understood heptavalent form. In this work, we monitored the formation of heptavalent neptunium [Np(VII)O 4 (OH) 2 ] 3− during ozonolysis of aqueous M OH ( M =Li, Na, K) solutions using a combined experimental and theoretical approach. All experimental reactions were closely monitored via absorption and vibrational spectroscopy to follow both the oxidation state and the speciation of neptunium guided by the calculated vibrational frequencies for various neptunium species. The mechanism of the reaction partly involves oxidative dissolution of transient Np(VI) oxide/hydroxide solid phases, the identity of which are dependent on the co‐precipitating counter‐cation Li + /Na + /K + . Additional calculations suggest that the most favorable energetic pathway occurs through the reaction of a [Np(V)O 2 (OH) 4 ] 3− with the hydroxide radical to form [Np(VI)O 2 (OH) 4 ] 2− , followed by an additional oxidation with HO⋅ to create [Np(VII)O 4 (OH) 2 ] 3− .

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

Disentangling Acid‐Base Chemistry through Blue Shifting Hydrogen Bond Contributions

The blue shifting of vibrational frequencies in hydrogen bonded molecules, as observed in aqueous environments, has been attributed to local partial charge transfer from solvation. Here, in this study, we extrapolate the blue shift model to the stronger ionic interactions between hydrogen bond acceptors associated with protonation through augmented pH levels and competitive interactions with counter ion pairing. The chemical model we utilize in this work is the aqueous pyridine‐pyridinium equilibrium to characterize the blue shifts observed in the pyridinium chloride ionic system. The observed agreement between observed experimental and calculated spectral shifts shows that the blue shifting model can be extrapolated to stronger interactions and accurately describe the nature of the hydrogen bond.

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

Butterfly, Vinylidene-Like, Monobridged and Trans Structures of Si 2 H 2 + : Comparison to the Well-Characterized Neutral Si 2 H 2

This paper investigates four of the constitutional isomers of Si 2 H 2 + , namely the butterfly, vinylidene-like, monobridged, and trans structures. These isomer geometries were all studied using the CCSD(T) method with basis sets as large as cc-pV5Z. Higher level methods CCSDT and CCSDT(Q) were used for final energetics. It is found that the butterfly isomer has the lowest energy, followed by vinylidene-like and monobridged structures; the trans structure has the highest energy of the four. All structures were compared with their neutral counterparts. Partial charges are computed and it is found that all isomers share the positive charge between the silicon atoms. NBO analysis shows that the cation Si-Si bond order is reduced by 0.5 relative to the neutral versions of the same isomer in the case of the butterfly, vinylidene-like, and monobridged; and reduced by 0.6 for the trans isomer. Vibrational frequencies, infrared intensities, and dipole moments are predicted to encourage the spectroscopic identification of Si 2 H 2 + .

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

Ab initio exploration of low–lying electronic states of linear and bent MNX + (M = Ca, Sr, Ba, Ra; X = O, S, Se, Te, Po) and their origins

High-level multireference and coupled cluster quantum calculations were employed to analyze low-lying electronic states of linear-MNX + and side-bonded-M[NX] + (M = Ca, Sr, Ba, Ra; X = O, S, Se, Te, Po) species. Their full potential energy curves (PECs), dissociation energies (D e s), geometric parameters, excitation energies (T e s), and harmonic vibrational frequencies (ω e s) are reported. The first three chemically bound electronic states of MNX + and M[NX] + are 3 Σ – , 1 Δ, 1 Σ + and 3 A", 1 A', 1 A", respectively. The 3 Σ – , 1 Δ, 1 Σ + of MNX + originate from the M + ( 2 D) + NX( 2 Π) fragments, whereas the 3 A", 1 A', 1A" states of M[NX] + dissociate to M + ( 2 S) + NX( 2 Π) as a result of avoided crossings. The MNX + and M[NX] + are real minima on the potential energy surface and their interconversions are possible. The M 2+ NX – /M 2+ [NX] – ionic structure is an accurate representation for their low-lying electronic states. The D e s of MNX + species were found to depend on the dipole moment (μ) of the corresponding NX ligands and a linear relationship between these two parameters was observed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗