Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “first principles calculation”

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 1,729 records · Page 96

Measurements of Beam Spin Asymmetries in p+p0 and p´p0 Dihadron Production at CLAS12

Semi-Inclusive Deep Inelastic Scattering (SIDIS) is a powerful experimental tool for studying the internal structure and dynamics of the proton, revealing how quarks and gluons are distributed and interact within it. SIDIS describes a process where an elec tron scatters off one of the constituent quarks within the proton, causing it to undergo hadronization, creating multiple hadrons in the final state. Through factorization, the full process can be split into probabilistic components: one which describes the internal structure of the proton using Parton Distribution Functions (PDFs), and another which describes the hadronization process using Fragmentation Functions (FFs). These functions are non-perturbative quantities of Quantum Chromodynamics (QCD), meaning they cannot be calculated directly from first principles and must instead be extracted from experimental measurements. Acommon approach for accessing PDFs and FFs using SIDIS is to measure asymmetries. In this context, asymmetries correspond to subtle differences in the angular distribution of outgoing particles that arise when the spin orientation of the incoming beam or target is reversed. Because many of these effects only appear when spin is involved, they isolate specific, nuanced properties of the proton’s spin-structure that are otherwise hidden in spin averaged measurements. In practice, they show up as specific azimuthal modulations (e.g., sin ¿R, sin(¿h ´ ¿R)), whose amplitudes isolate convolutions of PDFs and FFs at leading and subleading twist. Non-zero asymmetries of these angular distributions can be traced back to unique combinations of PDFs and FFs, offering a way to probe them directly. In this work, we measure SIDIS by analyzing high energy electron-proton scattering events using the CLAS12 detector at Jefferson Lab. This study focuses on subset of SIDIS referred to as dihadron SIDIS, where pairs of hadrons — here p+p0 and p´p0 — are observed. We analyzed these dihadrons using detector data collected during Fall 2018 and Spring 2019, where longitudinally polarized electrons from the CEBAF accelerator were incident on a liquid hydrogen target. A photon classifier using a Gradient Boosted Trees (GBTs) architecture was trained using Monte Carlo simulations to reduce the amount of iv false combinatorial background p0’s. When deployed on experimental data, the model in creases our dihadron statistics by up to five-fold compared to previous CLAS12 p0 analyses. This work reports the first measurements of beam spin asymmetries for p+p0 and p´p0 dihadron production in SIDIS. The measured asymmetries offer new insights to the spin-dependent structure and dynamics within the proton, as well as the spin-dependent properties of quark fragmentation. Non-zero twist-3 sin¿R amplitudes are observed, pro viding sensitivity to the subleading twist PDF e(x). The PDF e(x) encodes quark-gluon correlations within the proton — a property that is otherwise inaccessible at leading twist. Additionally, this work measured significant twist-2 modulations carried by sin(¿h ´ ¿R) and sin(2¿h ´2¿R), providing experimental access to the helicity dihadron fragmentation function (DiFF) GK 1 . Because there is no equivalent quark helicity-dependent FF in single pion SIDIS, the DiFF GK 1 offers a unique lens into novel spin-dependent fragmentation. For instance, the twist-2 modulations observed in this study are enhanced by vector mesons created during fragmentation — a behavior predicted by phenomenological models. This study broadens our understanding of dihadron fragmentation, revealing new details about the flavor and charge dependence of hadronization.

Matousek, Gregory [Duke Univ., Durham, NC (United ↗

Doping Effects on the Ductility of a Lightweight Refractory High-Entropy Alloy: Grain Boundary and Bulk Lattice Aspects

Doping elements in small amounts often segregate to grain boundaries (GBs) in alloys and can significantly impact mechanical properties and performance. Refractory high-entropy alloys (RHEAs) are known for their poor ductility, especially at low temperatures. Promoting GB cohesion through segregation can be an effective approach to mitigate embrittlement. Here, in this study, first-principles density functional theory (DFT) calculations were performed to examine the effects of important interstitial dopants (O, B, C, and N) and substitutional dopants (Cr, Y, La, and Ce) on the Σ5(310) [001] tilt GB of a lightweight RHEA Nb 32.5 Ti 27.5 Mo 22.5 Ta 12.5 Hf 2.5 Zr 2.5 . The DFT calculations reveal that certain dopants, such as B, C, Cr, N, and Y, exhibit favorable GB strengthening effects by improving bonding interactions with the bulk alloy. The impact of doping on the ductility parameter of the bulk lattice, defined as the ratio of surface energy to unstable stacking fault energy for the {110} <111> slip system, was also studied; and the results show that doping reduces the intrinsic ductility of the alloy, decreasing the D-parameter from 2.90 to 2.55, depending on the specific dopant. The present findings provide a foundational understanding at atomic level of the effect of representative dopants on mechanical properties of RHEAs and can be used to guide future alloy design for improved mechanical properties.

deformation charge density↗

First-principles study on the mechanical and thermodynamic properties of MoNbTaTiW

Refractory high-entropy alloys (RHEAs) are emerging as new materials for high temperature structural applications because of their stable mechanical and thermal properties at temperatures higher than 2273 K. In this study, the mechanical properties of MoNbTaTiW REDEA are examined by applying calculations based on first-principles density functional theory (DFT) and using a large unit cell with 100 randomized atoms. The phase calculation of MoNbTaTiW with CALPHAD method shows the existence of a stable body-centered cubic structure at a high temperature and a hexagonal closely packed phase at a low temperature. The predicted phase, shear modulus, Young’s modulus, Poisson’s ratio, and hardness values are consistent with available experimental results. The linear thermal expansion coefficient, vibrational entropy, and vibrational heat capacity of MoNbTaTiW RHEA are investigated in accordance with Debye-Grüneisen theory. These results may provide a basis for future research related to the application of RHEAs.

36 MATERIALS SCIENCE↗

Deducing subnanometer cluster size and shape distributions of heterogeneous supported catalysts

Abstract Infrared (IR) spectra of adsorbate vibrational modes are sensitive to adsorbate/metal interactions, accurate, and easily obtainable in-situ or operando. While they are the gold standards for characterizing single-crystals and large nanoparticles, analogous spectra for highly dispersed heterogeneous catalysts consisting of single-atoms and ultra-small clusters are lacking. Here, we combine data-based approaches with physics-driven surrogate models to generate synthetic IR spectra from first-principles. We bypass the vast combinatorial space of clusters by determining viable, low-energy structures using machine-learned Hamiltonians, genetic algorithm optimization, and grand canonical Monte Carlo calculations. We obtain first-principles vibrations on this tractable ensemble and generate single-cluster primary spectra analogous to pure component gas-phase IR spectra. With such spectra as standards, we predict cluster size distributions from computational and experimental data, demonstrated in the case of CO adsorption on Pd/CeO 2 (111) catalysts, and quantify uncertainty using Bayesian Inference. We discuss extensions for characterizing complex materials towards closing the materials gap.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Exploiting Machine Learning to Efficiently Predict Multidimensional Optical Spectra in Complex Environments

The excited-state dynamics of chromophores in complex environments determine a range of vital biological and energy capture processes. Time-resolved, multidimensional optical spectroscopies provide a key tool to investigate these processes. Although theory has the potential to decode these spectra in terms of the electronic and atomistic dynamics, the need for large numbers of excited-state electronic structure calculations severely limits first-principles predictions of multidimensional optical spectra for chromophores in the condensed phase. In this work, we leverage the locality of chromophore excitations to develop machine learning models to predict the excited-state energy gap of chromophores in complex environments for efficiently constructing linear and multidimensional optical spectra. By analyzing the performance of these models, which span a hierarchy of physical approximations, across a range of chromophore–environment interaction strengths, we provide strategies for the construction of machine learning models that greatly accelerate the calculation of multidimensional optical spectra from first principles.

36 MATERIALS SCIENCE↗

Strengthening model development and effects of low diffusing solutes to coarsening resistance in aluminum alloys

A modified Orowan strengthening model is proposed to account for finite rod-shaped precipitates with hemispherical caps in aluminum alloy systems. A combined computational and experimental approach is used to study the influences of anisotropic Orowan looping and solute-dislocation interaction on temperature-dependent yield strength. Here, the strengthening model is validated with a dataset containing 297 experimental precipitate geometries, chemistries, temperatures, and strength measurements, and achieves a strong predictive correlation of 0.8713 with experimentally measured yield strengths. Under conditions that are applicable to coarsening, constant particle volume fraction and aspect ratio, the model predicts that short rod precipitates provide far superior strengthening effects compared to plate precipitates. A cast Al-Si-Mg-Cu alloy with rod-shaped Q-phase (Al 3 Cu 2 Mg 9 Si 7 ) precipitates was developed with a novel chemistry exploring the use of low-diffusivity elements (Mn, Ni, V, Zr) to limit precipitate coarsening. The thermodynamic behavior of Mn, Ni, V, and Zr across the Q-phase interface is examined using transmission electron microscopy (TEM), first-principles density-functional theory (DFT) calculations. and atom-probe tomography (APT). DFT calculations utilizing TEM identified Q-phase/Al-matrix interfaces show that Mn, Ni, V, and Zr preferentially segregate to the Q-phase precipitate boundaries which suggests inhibition of precipitate coarsening and higher strengths after temperature exposure. Atom-probe tomography confirms solute atom partitioning/segregation at the Q-phase/Al-matrix interface, found in the modified commercial AS7GU alloy (A356 +0.5%Cu), which supports these observations.

36 MATERIALS SCIENCE↗

Electron mobility of SnO 2 from first principles

The transparent conducting oxide SnO 2 is a wide bandgap semiconductor that is easily n-type doped and widely used in various electronic and optoelectronic applications. Experimental reports of the electron mobility of this material vary widely depending on the growth conditions and doping concentrations. In this work, we calculate the electron mobility of SnO 2 from first principles to examine the temperature and doping concentration dependence and to elucidate the scattering mechanisms that limit transport. We include both electron–phonon scattering and electron-ionized impurity scattering to accurately model scattering in a doped semiconductor. Furthermore, we find a strongly anisotropic mobility that favors transport in the direction parallel to the c-axis. At room temperature and intrinsic carrier concentrations, the low-energy polar-optical phonon modes dominate scattering, while ionized-impurity scattering dominates above 10 18 cm −3 .

36 MATERIALS SCIENCE↗

The Role of Surface Oxygen Vacancies in BiVO 4

Bismuth vanadate (BiVO 4 ) is a widely studied oxide in solar water splitting, known for its ease of synthesis, high charge extraction yields, and advantageous band alignment with water. We present a combined first-principles and experimental study of the electronic structure of the (010) surface of BiVO 4 aimed at disentangling the impact of the surface and bulk oxygen vacancies on the electronic structure and transport properties. We found that oxygen vacancies are deep donors at the surface as they are in the bulk; our calculations on defect and polaron formation energies suggest that, while polarons formed from oxygen vacancies in the bulk can contribute to conductivity, those at the surface likely do not. Our results also show that out-of-plane structural relaxations at the surface contribute to the relatively immobile nature of electron polarons derived from surface oxygen vacancies. The structural model derived from first-principles calculation was validated by comparing computed results with experimental measurements of single-crystal and epitaxially grown single-crystalline BiVO 4 samples. Finally, we also found a reasonably good agreement between our calculated and measured work functions for BiVO 4 samples with and without oxygen vacancies.

36 MATERIALS SCIENCE↗

Spin dynamics in the itinerant antiferromagnet SrCr 2 ⁢As 2

SrCr 2 ⁢As 2 is an itinerant antiferromagnet in the same structural family as the SrFe 2 ⁢As 2 high-temperature superconductors. Here, we report our calculations of exchange-coupling parameters 𝐽 𝑖⁢𝑗 for SrCr 2 ⁢As 2 using a static linear-response method based on first-principles electronic-structure calculations. We find that the dominant nearest-neighbor exchange coupling 𝐽 1 >0 is antiferromagnetic whereas the next-nearest-neighbor exchange coupling 𝐽 2 <0 is ferromagnetic with 𝐽 2 /𝐽 1 = −0.68, reinforcing the checkerboard in-plane magnetic structure. Thus, unlike other transition-metal arsenides based on Mn, Fe, or Co, we find no competing magnetic interactions in SrCr 2 ⁢As 2 , which aligns with experimental findings. Moreover, the orbital resolution of exchange interactions shows that 𝐽 1 and 𝐽 2 are dominated by direct exchange mediated by the Cr 𝑑 orbitals. To validate the calculations we conduct inelastic neutron-scattering measurements on powder samples that show steeply dispersive magnetic excitations arising from the magnetic Γ points and persisting up to energies of at least 175 meV. The spin-wave spectra are then modeled using the Heisenberg Hamiltonian with the theoretically calculated exchange couplings. In conclusion, the calculated neutron-scattering spectra are in good agreement with the experimental data.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

First-principles study of the T center in silicon

The T center in silicon is a well-known carbon-based color center that has been recently considered for quantum technology applications. Using first-principles computations, we show that the excited state is formed by a defect-bound exciton made of a localized defect state occupied by an electron to which a hole is bound. The localized state is of strong carbon p character and reminiscent of the localization of the unpaired electron in the ethyl radical molecule. The radiative lifetime for the defect-bound exciton is calculated to be on the order of microseconds, much longer than for other quantum defects such as the nitrogen vacancy center in diamond and in agreement with experiments. The longer lifetime is associated with the small transition dipole moment as a result of the very different nature of the localized and delocalized states forming the defect-bound exciton. Finally, we use first-principles calculations to assess the stability of the T center. Here, we find the T center to be stable against decomposition into simpler defects when keeping the stoichiometry fixed. However, we identify that the T center is easily prone to (de)hydrogenation and so requires very precise annealing conditions (temperature and atmosphere) to be efficiently formed.

36 MATERIALS SCIENCE↗

Reexamination of hydrodynamic phonon transport in thin graphite

The recent hydrodynamic phonon transport theory for graphitic materials has been supported by the measurements of the second sound at temperatures up to about 100 K. When boundary scattering becomes comparable to momentum-conserving normal phonon scattering processes that are responsible for phonon hydrodynamics, Poiseuille phonon flow phenomena can emerge to give rise to unique size-dependent thermal conductivity in thin graphite. Here, we examine the thickness range for the Poiseuille phonon flow to become observable in thin graphite with the use of both deviational Monte Carlo simulation of the Peierls-Boltzmann transport equation and four-probe thermal transport measurements. As the basal-plane thermal conductivity calculated by prior first-principles theories saturates to the graphite value when the thickness is increased to five graphene layers, the phonon dispersion of graphite is used in the current calculations of thin graphite of micrometer thickness and a 23-layer thick ultrathin graphite (UTG) sample. The calculations show that diffuse surface scattering by surface defects can lead to Poiseuille phonon flow at 50 K in thin graphite with the thickness close to several micrometers but not in the 65 μm thin graphite and 23-layer UTG, where phonon scattering by the top and bottom surfaces become, respectively, much less and more frequent than the normal processes. In addition, the calculation results with the bulk graphite dispersion and diffuse surface scattering show decreased basal-plane thermal conductivity with decreasing thickness, opposite to recent thermocouple measurements of thin graphite samples. In comparison, the calculation results reveal that partially diffuse surface defect scattering can yield the four-probe measurement results of UTG samples, which are prepared here with an improved process to minimize surface contamination.

36 MATERIALS SCIENCE↗

Two-dimensional Dirac semimetal based on the alkaline earth metal CaP 3

Using an evolutionary algorithm in combination with first-principles density-functional theory calculations, we identify a two-dimensional (2D) CaP 3 monolayer as a new Dirac semimetal due to inversion and nonsymmorphic spatial symmetries of the structure. This new topological material, composed of light elements, exhibits high structural stability (higher than the phase known in the literature), which is confirmed by thermodynamic and kinetic stability analysis. Moreover, it satisfies the electron filling criteria, so that its Dirac state is located near the Fermi level. The existence of the Dirac state predicted by the theoretical symmetry analysis is also confirmed by first-principles electronic band structure calculations. We find that the energy position of the Dirac state can be tuned by strain, while the Dirac state is unstable against an external electric field since it breaks the spatial inversion symmetry. In conclusion, our findings should be instrumental in the development of 2D Dirac fermions based on light elements for their application in nanoelectronic devices and topological electronics.

2-dimensional systems↗

First-Principles Modeling Of Electromagnetic Scattering By Discrete and Discretely Heterogeneous Random Media

A discrete random medium is an object in the form of a finite volume of a vacuum or a homogeneous material medium filled with quasi-randomly and quasi-uniformly distributed discrete macroscopic impurities called small particles. Such objects are ubiquitous in natural and artificial environments. They are often characterized by analyzing theoretically the results of laboratory, in situ, or remote-sensing measurements of the scattering of light and other electromagnetic radiation. Electromagnetic scattering and absorption by particles can also affect the energy budget of a discrete random medium and hence various ambient physical and chemical processes. In either case electromagnetic scattering must be modeled in terms of appropriate optical observables, i.e., quadratic or bilinear forms in the field that quantify the reading of a relevant optical instrument or the electromagnetic energy budget. It is generally believed that time-harmonic Maxwell's equations can accurately describe elastic electromagnetic scattering by macroscopic particulate media that change in time much more slowly than the incident electromagnetic field. However, direct solutions of these equations for discrete random media had been impracticable until quite recently. This has led to a widespread use of various phenomenological approaches in situations when their very applicability can be questioned. Recently, however, a new branch of physical optics has emerged wherein electromagnetic scattering by discrete and discretely heterogeneous random media is modeled directly by using analytical or numerically exact computer solutions of the Maxwell equations. Therefore, the main objective of this Report is to formulate the general theoretical framework of electromagnetic scattering by discrete random media rooted in the Maxwell- Lorentz electromagnetics and discuss its immediate analytical and numerical consequences. Starting from the microscopic Maxwell-Lorentz equations, we trace the development of the first principles formalism enabling accurate calculations of monochromatic and quasi-monochromatic scattering by static and randomly varying multiparticle groups. We illustrate how this general framework can be coupled with state-of-the-art computer solvers of the Maxwell equations and applied to direct modeling of electromagnetic scattering by representative random multi-particle groups with arbitrary packing densities. This first-principles modeling yields general physical insights unavailable with phenomenological approaches. We discuss how the first-order-scattering approximation, the radiative transfer theory, and the theory of weak localization of electromagnetic waves can be derived as immediate corollaries of the Maxwell equations for very specific and well-defined kinds of particulate medium. These recent developments confirm the mesoscopic origin of the radiative transfer, weak localization, and effective-medium regimes and help evaluate the numerical accuracy of widely used approximate modeling methodologies.

Michael I Mishchenko↗

First-principles characterization of the magnetic properties of Cu 2 (OH) 3 Br

Low-dimensional spin-1/2 transition metal oxides and oxyhalides continue to be at the forefront of research investigating nonclassical phases such as quantum spin liquids. In this work, we examine the magnetic properties of the oxyhalide Cu 2 (OH) 3 Br in the botallackite structure using first-principles density functional theory, linear spin-wave theory, and exact diagonalization calculations. This quasi-two-dimensional system consists of Cu 2+ S = 1/2 moments arranged on a distorted triangular lattice. Our exact diagonalization calculations, which rely on a first-principles-based magnetic model, generate spectral functions consistent with inelastic neutron scattering data. By performing computational experiments to disentangle the chemical and steric effects of the halide ions, we find that the dominant effect of the halogen ions is steric in the Cu 2 (OH) 3 X series of compounds.

36 MATERIALS SCIENCE↗

Black hole induced false vacuum decay from first principles

We provide a method to calculate the rate of false vacuum decay induced by a black hole. The method uses complex tunneling solutions and consistently takes into account the structure of different quantum vacua in the black hole metric via boundary conditions. The latter are connected to the asymptotic behavior of the time-ordered Green’s function in the corresponding vacua. We illustrate the technique on a two-dimensional toy model of a scalar field with inverted Liouville potential in an external background of a dilaton black hole. We analytically derive the exponential suppression of tunneling from the Boulware, Hartle-Hawking and Unruh vacua and show that they are parametrically different. The Unruh vacuum decay rate is exponentially smaller than the decay rate of the Hartle-Hawking state, though both rates become unsuppressed at high enough black hole temperature. We interpret the vanishing suppression of the Unruh vacuum decay at high temperature as an artifact of the two-dimensional model and discuss why this result can be modified in the realistic case of black holes in four dimensions.

2D Gravity↗

Electronic and optical properties of transition metal impurities in GaN and AlN

The optical properties of Cr(III) impurities are well studied and exploited in 𝛼−Al 2 ⁢O 3 , which enables ruby lasers. Photoluminescence spectra resembling the famous 1.78 eV emission of Cr in 𝛼−Al 2 ⁢O 3 have been observed in other wide-band-gap host materials, in particular AlN. However, it is unclear whether Cr would produce the same emission in the nitrides as in 𝛼−Al 2 ⁢O 3 , especially given the different local environment (tetrahedral vs octahedral coordination). We investigate the electronic and optical properties of transition metal substitutional impurities (Cr, Co, Ni, and Cu) in wurtzite GaN and AlN from first principles, utilizing hybrid density functional theory. We calculate the optical transition energies using a constrained occupation approach that accounts for the multideterminant nature of the excited states. Our results indicate that substitutional Cr(III) is not responsible for the emission near 1.78 eV observed in the nitrides. Among the examined impurities, only substitutional Co(II) produces emission near 1.78 eV. Based on the favorable properties of Co(II) (low formation energy, small Huang-Rhys factor, and an 𝑆 = 3/2 ground state) we propose it as a promising single-photon emitter or spin qubit. We additionally assign the PL peak observed near 1.05 eV in GaN and AlN to the 2 𝐸→ 4 𝐴 2 transition of the Ni(III) substitution.

36 MATERIALS SCIENCE↗

Thermal Neutron Scattering Law for Beryllium Hydride and Critical Mass Calculations

The thermal neutron scattering law (TSL) for crystalline beryllium hydride (BeH 2 ) is developed from first-principles ab initio lattice dynamics calculations and the impact of neutron thermalization in this material on critical mass is estimated. BeH 2 has a body-centered orthorhombic crystal structure with 12 molecules per unit cell and a theoretical density of 0.755 g/cm 3 . The vibrational (phonon) densities of states for H and Be bound in BeH 2 are determined using VASP density functional theory and PHONON lattice dynamics calculations. The TSLs for H bound in BeH 2 , H(BeH 2 ), and Be bound in BeH 2 , Be(BeH 2 ), are then evaluated in the incoherent approximation from the calculated H and Be partial phonon density of states using FLASSH. Finally, critical mass as a function of 235 U loading density for bare and reflected BeH 2 moderated spheres is predicted from MC21 Monte Carlo neutron transport calculations using ENDF/B-VIII.0 cross sections and the H(BeH 2 ) and Be(BeH 2 ) TSL evaluations. Comparisons are made to water (H 2 O), polyethylene (CH 2 ), and beryllium oxide (BeO) as moderators. These critical mass predictions are a refinement upon the prior work by Rao and Srinivasan that neglected thermal neutron scattering effects. The minimum critical mass of a BeH 2 moderated assembly is estimated to be 0.207 kg 235 U for a 0.20 m thick BeO reflected sphere and 0.178 kg 235 U for a 0.40 m thick BeO reflected sphere.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Study of Tritium Diffusivity in Pure and Sn-Defective Zr: A First-Principles Density Functional Theory Approach

Zirconium alloys (e.g., zircaloy-4) are used as tritium ( 3 H) getter materials in tritium-producing burnable absorber rods (TPBARs) owing to their ability to capture 3 H and chemically convert them into metal hydrides. Understanding of 3 H diffusion mechanisms in zircaloy is crucial for the optimal design of material performance in nuclear technology. Here, for this work, we perform first-principles density functional theory calculations to study the 3 H diffusion mechanism in pure and impure Zr with a low concentration of tin (Sn) atoms to determine the impact of the presence of Sn on the movement of 3 H atoms through the material. First, we calculated the diffusion barriers for 3 H in pure Zr by taking different migration pathways. We then introduced a low concentration of Sn impurity and systematically explored the impurity effect on the diffusion barriers for 3 H. Using our calculated diffusion energy barriers, we further obtained the diffusion coefficients and analyzed the results by comparing them with the experimental and previously calculated values. A diffusion coefficient of the order of 10 –8 m 2 /s is predicted. We also found that the presence of a Sn impurity could reduce the diffusivity up to 4 orders of magnitude. In conclusion, our results could serve as guidelines for further experimental investigations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗