Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Solvation”

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 883 records · Page 49

Toward Quantitatively Accurate Calculation of the Redox-Associated Acid–Base and Ligand Binding Equilibria of Aquacobalamin

In redox processes in complex transition metal-containing species are often intimately associated with changes in ligand protonation states and metal coordination number. Moreover, a major challenge is therefore to develop consistent computational approaches for computing pH-dependent redox and ligand dissociation properties of organometallic species. Reduction of the Co center in the vitamin B12 derivative aquacobalamin can be accompanied by ligand dissociation, protonation, or both, making these properties difficult to compute accurately. We examine this challenge here by using density functional theory and continuum solvation to compute Co ligand binding equilibrium constants (Kon/off), pKas and reduction potentials for models of aquacobalamin in aqueous solution. We consider two models for cobalamin ligand coordination: the first follows the hexa, penta, tetra coordination scheme for Co III , Co II , and Co I species, respectively, and the second model features saturation of each vacant axial coordination site on Co II and Co I species with a single, explicit water molecule to maintain six directly interacting ligands or water molecules in each oxidation state. Comparing these two coordination schemes in combination with five dispersion-corrected density functionals, we find that the accuracy of the computed properties is largely independent of the scheme used, but including only a continuum representation of the solvent yields marginally better results than saturating the first solvation shell around Co throughout. PBE performs best, displaying balanced accuracy and superior performance overall, with RMS errors of 80 mV for seven reduction potentials, 2.0 log units for five pK a s and 2.3 log units for two log K on/off values for the aquacobalamin system. Furthermore, we find that the BP86 functional commonly used in corrinoid studies suffers from erratic behavior and inaccurate descriptions of Co axial ligand binding, leading to substantial errors in predicted pK a s and K on/off values. Finally, these findings demonstrate the effectiveness of the present approach for computing electrochemical and thermodynamic properties of a complex transition metal-containing cofactor.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The Influence of Water on the Alkyl Region Structure in Variable Chain Length Imidazolium-Based Ionic Liquid/Water Mixtures

In this work, solutions of room temperature ionic liquids (RTILs) and water were studied by observing the reorientational dynamics of the fluorescent probe perylene. Perylene is solvated in the alkyl regions of the RTILs. Its D 2h symmetry made it possible to extract dynamical information on both in-plane and out-of-plane reorientation from time resolved fluorescence anisotropy measurements. Perylene reorientation reports on its interactions with the alkyl chains. The RTILs were a series of 1-alkyl-3-methylimidazolium tetrafluoroborates (C n mimBF 4 , where n is the number of carbons in the alkyl chain), and the effects on perylene’s dynamics were observed when varying the alkyl chain length of the cation (n = 4, 6, 8, and 10; butyl, hexyl, octyl, decyl) and varying the water content from pure RTIL to roughly three water molecules per RTIL ion pair. Time correlated single photon counting was used to measure the fluorescence anisotropy decays to determine the orientational dynamics. The friction coefficients for both the in-plane and out-of-plane reorientation were determined to eliminate the influence of changes in viscosity caused by both the addition of water and the different alkyl chain lengths. The friction coefficients provided information on the interactions of the perylene with it alkyl environment and how these interactions changed with chain length and water content. As chain length increased, the addition of water had less of an effect on the local alkyl environment surrounding the perylene. The friction coefficients generally increased with higher water contents; the inplane orientational motion was hindered significantly more than the out-of-plane motion. The restructuring of the alkyl regions is likely a consequence of a rearrangement of the ionic imidazolium head groups to accommodate partial solvation by water, which results in a change in the arrangement of the alkyl chains. At very high water content, BmimBF4 broke this general trend, with both in-plane and out-of-plane rotational friction decreasing above a water content of one water per ion pair. This decrease indicates a major reorganization of the overall liquid structure in high water content mixtures. In contrast to BmimBF 4 , the longer chain length RTILs are not infinitely miscible with water, and do not show evidence of a major reorganization before reaching saturation and phase-separating. The results suggest that phase separation in longer chain length BF 4 RTILs is a consequence of their inability to undergo the reorganization of the alkyl regions necessary to accommodate high water concentrations.

14 SOLAR ENERGY↗

Restrained Proton Indicator in Combined Quantum-Mechanics/Molecular-Mechanics Dynamics Simulations of Proton Transfer through a Carbon Nanotube

Recently, a collective variable “proton indicator” was purposed for tracking an excess proton solvated in bulk water in molecular dynamics simulations. In this work, we demonstrate the feasibility of utilizing the position of this proton indicator as a reaction coordinate to model an excess proton migrating through a hydrophobic carbon nanotube in combined quantum-mechanics/molecular-mechanics simulations. Our results indicate that applying a harmonic restraint to the proton indicator in the bulk solvent near the nanotube pore entrance leads to the recruitment of water molecules into the pore. This is consistent with an earlier study that employed a multistate empirical valence bond potential and a different representation (center of excess charge) of the proton. We attribute this water recruitment to the delocalized nature of the solvated proton, which prefers to be in high-dielectric bulk solvent. While water recruitment into the pore is considered an artifact in the present simulations (because of the artificially imposed restraint on the proton), if the proton were naturally restrained, it could assist in building water wires prior to proton transfer through the pore. The potential of mean force for a proton translocation through the water-filled pore was computed by umbrella sampling, where the bias potentials were applied to the proton indicator. The free energy curve and barrier heights agree reasonably with those in the literature. Finally, the results suggest that the proton indicator can be used as a reaction coordinate in simulations of proton transport in confined environments.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Comprehensive Experimental and Computational Spectroscopic Study of Hexacyanoferrate Complexes in Water: From Infrared to X-ray Wavelengths

We present a joint experimental and computational study of the hexacyanoferrate aqueous complexes at equilibrium in the 250 meV to 7.15 keV regime. The experiments include the vibrational spectroscopy of the cyanide ligands, the valence electronic absorption spectra and Fe 1s core hole spectra using X-ray absorption and emission techniques. These experiments are accompanied by density functional theory (DFT) based quantum mechanics/molecular mechanics (QM/MM) molecular dynamics (MD) simulations and spectroscopy calculations of the experiments spanning the IR to X-ray wavelengths. Our simulations, including explicit solvent effects, across this large energy window allows for a systematic comparison of vibrational, valence electronic and core-level spectra. The analysis reveals how the local solvation environment influences the geometric and electronic structure and impacts the line shape and intensity of the observed spectral features in the infrared, optical and X-ray wavelengths.

Ross, Matthew↗

Ab Initio Molecular Dynamics Reveal Spectroscopic Siblings and Ion Pairing as New Challenges for Elucidating Prenucleation Aluminum Speciation

We report the characterization of prenucleation species is essential to understand crystallization mechanisms across many chemical systems and often involves the use of vibrational spectroscopy. Nowhere is this more evident than in the development of “green” aluminum processing technologies, where detailed understanding of the speciation of aluminum and its polynuclear analogues in highly alkaline, low water solutions is elusive. The aluminate anion Al(OH) 4 – predominates in alkaline conditions, yet equilibrium with dimeric species, either μ-oxo Al 2 O(OH) 6 2– or di-μ-hydroxo Al 2 (OH) 8 2– , can be assumed. Using ab initio molecular dynamics with full solvation and the presence of counterions, this work reconciles previous contradictory studies that had concluded only a single species under relevant solution conditions. We reveal that the two dimers are energetically separated by 2 kcal/mol in pure water but that the stability of each can be reversed by ion pairing expected in saturated salt solutions. Simulated Raman and IR spectra for each species (accounting for anharmonicity and the fluctuating solvating environment) provide the first proof that the considered species are “spectroscopic siblings”, whose multiple overlapping bands prevent definitive assertions in terms of speciation when compared to the experimental spectra. Lastly, these observations are likely to hold in higher order aluminate oligomers and as such present a massive challenge toward understanding the crystallization mechanisms relevant to aluminum processing.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Rate versus Free Energy Change for Attaching Highly Mobile Electrons to Molecules in Nonpolar Liquids

The inverted region of the Marcus theory, usually absent for bimolecular electron transfer reactions, is clearly observed for electron attachment reactions to molecules in nonpolar fluids. Application of pressure increased the energies of the solvated electron reactants letting us continuously adjust the free energy change. Inverted behavior is enabled by the very high mobilities of the solvated electrons which raise the diffusion-controlled encounter rates so high that they do not limit the reaction rates. The nonpolar media used in these experiments reduce reorganization energies, enhancing inverted behavior. Still, for every case showing an inverted region, the presence of low-lying excited states in the product radical anions led to regions of increasing rate constants that began at the energies of excited states of –0.54 to –1.2 eV. Here, while continuum models predict no solvent reorganization energy in nonpolar liquids, fits to the data found solvent reorganization energies of 0.05–0.4 eV supporting ideas advanced in theories of Matyushov.

14 SOLAR ENERGY↗

Electron Trapping and Ion Leaching at the Li-Modified Quartz–Water Interface

Electrons and ions at metal oxide-water interfaces have a critical role in many phenomena and applications, thus making their properties of considerable interest. We here investigate the behavior of a lithium impurity at the water interface with Li-doped SiO 2 using hybrid density functional based ab initio molecular dynamics simulations. We find that the excess electron donated by the lithium dopant localizes on a surface silicon atom and is partially solvated at the aqueous interface. While the position of the excess electron does not change during our ~60 ps simulation, the lithium cation diffuses from its initial position in the SiO 2 subsurface toward the interface and eventually leaches out of the oxide surface and becomes solvated by four water molecules, forming an aqueous Li+- electron complex at the interface. The degree of interaction between the localized electron and the lithium ion controls not only the energy level of the excess electron but also the extent of structural distortion and hydrogen-bonding at the interface. The results of our study can be relevant for understanding the initial stages of glass corrosion, where alkali ions are leached from the surface region by interactions with water.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Surface-Mediated Solvent Decomposition in Li–Air Batteries: Impact of Peroxide and Superoxide Surface Terminations

A viable Li/O 2 battery will require the development of stable electrolytes that do not continuously decompose during cell operation. In some recent experiments it is suggested that reactions occurring at the interface between the liquid electrolyte and the solid lithium peroxide (Li 2 O 2 ) discharge phase are a major contributor to these instabilities. To clarify the mechanisms associated with these reactions, a variety of atomistic simulation techniques, classical Monte Carlo, van der Waals-augmented density functional theory, ab initio molecular dynamics, and various solvation models, are used to study the initial decomposition of the common electrolyte solvent, dimethoxyethane (DME), on surfaces of Li 2 O 2 . Comparisons are made between the two predominant Li 2 O 2 surface charge states by calculating decomposition pathways on peroxide-terminated (O 2 2– ) and superoxide-terminated (O 2 1– ) facets. For both terminations, DME decomposition proceeds exothermically via a two-step process comprised of hydrogen abstraction (H-abstraction) followed by nucleophilic attack. In the first step, abstracted H dissociates a surface O 2 dimer, and combines with a dissociated oxygen to form a hydroxide ion (OH – ). In the remaining surface oxygen then attacks the DME, resulting in a DME fragment that is strongly bound to the Li 2 O 2 surface. DME decomposition is predicted to be more exothermic on the peroxide facet; nevertheless, the rate of DME decomposition is faster on the superoxide termination. The impact of solvation (explicit vs implicit) and an applied electric field on the reaction energetics are investigated. Finally, our calculations suggest that surface-mediated electrolyte decomposition should out-pace liquid-phase processes such as solvent auto-oxidation by dissolved O 2 .

25 ENERGY STORAGE↗

Effects of High and Low Salt Concentration in Electrolytes at Lithium–Metal Anode Surfaces

The use of high concentration salts in electrolyte solutions of lithium-sulfur (Li-S) batteries has been shown beneficial for mitigating some effects such as polysulfide shuttle and dendrite growth at the Li metal anode. Such complex solutions have structural, dynamical, and reactivity associated issues that need to be analyzed for a better understanding of the reasons behind such beneficial effects. A passivation interfacial layer known as solid-electrolyte interphase (SEI) is generated during battery cycling as a result of electron transfer from the metal anode causing electrolyte decomposition. Here in this work, we investigate using density functional theory and ab initio molecular dynamics simulations the salt decomposition, solvation effects, interactions among intermediate products and other species, and potential components of the SEI layer as a function of chemical nature and concentration of the salt, for lithium bis(trifluoromethanesulfonyl)imide (LiTFSI) and lithium bis(fluorosulfonyl)imide (LiFSI) at 1M and 4M concentrations in dimethoxyethane. It is found that LiTFSI undergoes a less complete reduction and facilitates charge transfer from the anode, whereas LiFSI shows a more complete decomposition forming LiF as one of the main SEI products. In addition, the specific decomposition mechanisms of each salt clearly point to the initial SEI components and the potential main products derived from them. Finally, very complex networks are found among the salt and solvent molecules in their attempt to maximize Li ion solvation that is quantified through the determination of coordination numbers.

25 ENERGY STORAGE↗

First-Principles Investigation of Lithium Polysulfide Structure and Behavior in Solution

We present the Lithium-Sulfur battery is a promising next generation energy storage technology that could meet the demands of modern society with a theoretical specific energy near 2500 W h kg -1 . However, this battery chemistry faces unique problems such as the parasitic polysulfide shuttle reaction which hinders battery performance severely. This shuttle phenomenon is caused by solubilities of intermediate reaction products in the electrolyte during the reduction chemistry of the battery. With molecular simulation and computational chemistry tools, we studied the thermodynamics, solvation structure, and dynamics of the long-chain lithium polysulfide species Li 2 S 6 and Li 2 S 8 in dimethoxyethane and 1,3-dioxolane to gain a deeper fundamental understanding of this process. We determined the structure of the 1st solvation shell for Li + as well as those of Li 2 S 6 , Li 2 S 8 closed and Li 2 S 8 linear in pure solvents and solvents with extra Li + added. Finally, the lithium polysulfide species were found not to favor dissociation and would most likely exist as fully lithiated species in solution.

25 ENERGY STORAGE↗

Atomic-Scale Quantification of Interfacial Binding between Peptides and Inorganic Crystals: The Case of Calcium Carbonate Binding Peptide on Aragonite

Using a specific explicitly solvated interface model between a calcium carbonate binding peptide and crystalline aragonite, we investigate the electronic structure, atomic bonding, solvation effect, and the role of hydrogen bonding on the cohesion, stability, and functionality of this complex hybrid system using density functional calculation. The large interface model is strategically constructed using a stepwise procedure followed by ab initio molecular dynamics to obtain the optimal conformation. The calculated data on the electronic structure and bonding are analyzed in terms of three structural parts: aragonite, peptide, and water. Next, we focus on the binding between aragonite (001) surface and the peptide mediated by water. Finally, specific interatomic bonding between the amino acids in peptide and the (001) surface of aragonite is quantified. A single quantum mechanical metric, the total bond order density (TBOD), infers the dynamic interplay of different competing interactions. Four amino acids HIS1, ARG6, MET7, and TRP11 in the peptide sequence have strong interfacial Ca–O bonding and O···H hydrogen bonding between aragonite and peptide. The calculated Young’s modulus 33.37 GPa is in line with the measured value for nacre. Our approach for interfacial study between aragonite and a calcium carbonate binding peptide offers a broad perspective for probing complex interactions between the biomimetic interfaces. Finally, TBOD can be used as an effective parameter in ranking the efficacy of peptide–surface interactions and in providing a programmable design for bio-inspired material interfaces based on computational means.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Effect of Solvent and Substrate on the Surface Binding Mode of Carboxylate-Functionalized Aromatic Molecules

Here, the efficiency of dye-sensitized solar cells (DSSCs) is strongly influenced by dye molecule orientation and interactions with the substrate. Understanding the factors controlling the surface orientation of sensitizing organic molecules will aid in the improvement of both traditional DSSCs and other devices that integrate molecular linkers at interfaces. Here, we describe a general approach to understand relative dye–substrate orientation and provide analytical expressions predicting orientation. We consider the effects of substrate, solvent, and protonation state on dye molecule orientation. In the absence of solvent, our model predicts that most carboxylic acid-functionalized molecules prefer to lie flat (parallel) on the surface, due to van der Waals interactions, as opposed to a tilted orientation with respect to the surface that is favored by covalent bonding of the carboxylic acid group to the substrate. When solvation effects are considered, however, the molecules are predicted to orient perpendicular to the surface. We extend this approach to help understand and guide the orientation of metal–organic framework (MOF) thin-film growth on various metal–oxide substrates. A two-part analytical model is developed on the basis of the results of DFT calculations and ab initio MD simulations that predicts the binding energy of a molecule by chemical and dispersion forces on rutile and anatase TiO 2 surfaces, and quantifies the dye solvation energy for two solvents. The model is in good agreement with the DFT calculations and enables rapid prediction of dye molecule and MOF linker binding preference on the basis of the size of the adsorbing molecule, identity of the surface, and the solvent environment. We establish the threshold molecular size, governing dye molecule orientation, for each condition.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Theoretical Approaches to Describing the Oxygen Reduction Reaction Activity of Single-Atom Catalysts

Single-atom catalysts have recently emerged as promising low-cost alternatives to Pt for the oxygen reduction reaction (ORR). Given the unique properties that distinguish these systems from traditional transition-metal electrocatalysts, it is essential to benchmark and establish appropriate computational approaches to study these novel materials. Herein, we employ multiple levels of theory, including wave function methods, density functional theory (DFT), and classical simulations, to investigate Cu-modified covalent triazine framework catalysts (Cu/CTF). We consider three major aspects of treating this system computationally. First, we present a step-wise approach to predict the ORR mechanism and adsorbate coverages on Cu/CTF. We then benchmark various DFT methods to coupled-cluster theory with the domain-based local pair natural orbital approximation, which indicates that HSE06 and PBE0 hybrid functionals most accurately describe the adsorption energies of ORR adsorbates on Cu/CTF. We finally employ thermodynamic integration and other techniques to consider solvation effects, which play significant roles in predicting the energies of reaction intermediates and the overall ORR pathway. Lastly, our findings indicate that accurate descriptions of both the electronic structure and solvation are necessary to understand the ORR activity of Cu/CTF.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Solvent–Adsorbate Interactions and Adsorbate-Specific Solvent Structure in Carbon Dioxide Reduction on a Stepped Cu Surface

In this work, the structure of water and its interactions with various carbon dioxide reduction intermediates adsorbed on a Cu(211) surface is investigated using density functional theory. We find that the presence of adsorbates has a significant and adsorbate-specific effect on the local water structure and that solvation can stabilize adsorbate conformations different from those found in vacuum. We describe relationships between the hydrogen bonding capability of an adsorbate, the dipole moment of the adsorbate, the energetic strength of water–adsorbate interactions, and the change induced in the local water orientation by the adsorbate. Mechanistic implications are discussed. Here, we investigate and quantify the error associated with using arbitrary locally optimized solvent structures in calculations of relevant physical quantities, such as solvated binding energies and work functions. Possible effects of thermal motion on calculations of the work function are investigated using ab initio molecular dynamics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Molecular Investigation of Oxidized Graphene: Anatomy of the Double-Layer Structure and Ion Dynamics

We investigated the influence of surface oxidization of planar graphene electrodes on charge storage and ion dynamics of supercapacitors. Our approach compared two distinct ionic liquid (IL) electrolytes: tetraethylammonium tetrafluoroborate solvated in acetonitrile (TEA-BF4/ACN) and 1-ethyl-3-methylimidazolium bis(trifluoromethylsulfonyl)imide solvated in acetonitrile (EMIm-TFSI/ACN). Both experimental electrochemical tests and molecular dynamics (MD) simulations showed positive, electrolyte-specific influences of hydroxyl-free electrode interfaces on capacitance. In the EMIm-TFSI/ACN system, the hydroxylated surface, because of its strong interaction with anions, impeded surface charge storage. On the other hand, in the case of TEA-BF4/ACN, the distribution and orientation of ACN across the system exerted vital influence on the capacitance, especially on the positive hydroxyl-free electrode. Furthermore, MD simulations of ion mobility with respect to the electrode surface in the lateral and perpendicular directions revealed significantly slower diffusion performance on the oxidized surface. Our efforts enhanced the level of fundamental understanding of the effects of hydroxyl groups on electrode–electrolyte interfaces and resulting supercapacitor performance.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Efficient Intermolecular Energy Exchange and Soft-ionization of Water at Nanoplatelet Interfaces

X-ray, energetic photon and electron irradiation can ionize and electronically excite target atoms and molecules. These excitations undergo complicated relaxation and energy transfer processes that ultimately determine the manifold of system responses to the deposited excess energy. In weakly bound gas- and solution-phase samples, intermolecular Coulomb decay (ICD) and electron-transfer mediated decay (ETMD) can occur with neighboring atoms or molecules leading to efficient transfer of the excess energy to the surroundings. In ionic solids such as metal-oxides, intra- and inter-atomic Auger decay produces localized final states that lead to lattice damage and typically removal of cations from the substrate. The relative importance of Auger stimulated damage (ASD) versus ICD and ETMD in micro-solvated nanoparticle interfaces is not known. Though ASD is generally expected, essentially no lattice damage and very efficient ICD and ETMD resulting from ionization and electronic excitation of micro-solvated boehmite (AlOOH) nanoplatelets has been observed. Efficient energy transfer and soft ionization of interfacial molecules is likely a general phenomenon at gas-oxyhydroxide nanoparticle interfaces where the density of states of the ionized chemisorbed species significantly overlap the core hole states of the solid.

Jones, Brant M.↗

Direct Correlation of the Salt-Reduced Diffusivities of Organic Solvents with the Solvent’s Mole Fraction

We report that lithium bis(trifluoromethanesulfonyl)imide (LiTFSI) in organic solvents (especially propylene carbonate) has demonstrated extraordinary pseudocapacitive performance as an electrolyte in the supercapacitor configuration. However, the influence of the solvated ions on the diffusivity of the solvent molecules is yet to be understood. We examine the impact of LiTFSI on the diffusivity in five organic solvents: acetonitrile (ACN), tetrahydrofuran (THF), methanol (MeOH), dimethyl sulfoxide (DMSO), and propylene carbonate (PC) using a combination of neutron scattering, conductivity measurements, and molecular dynamics simulations. The extent of the diffusivity reduction in the concentration regime of ≤1 M directly correlates with the solvent mole fraction at which the solvation shells around Li + ions are of similar size in all the solvents, resulting in a universal ~50% reduction in the solvent diffusivity. These results provide guidance for formulation of the new electrolytes to enhance the performance of energy storage devices.

conductivity↗

Near-Quantitative Predictions of the First-Shell Coordination Structure of Hydrated First-Row Transition Metal Ions Using K-Edge X-ray Absorption Near-Edge Spectroscopy

Understanding the solvation structure of transition metal ions is important due to its broad impact on the kinetics, stability, and reactivity in a variety of applications in geochemistry, biochemistry, energy storage, and environmental chemistry. Using water as the common ligand, we study the X-ray absorption pre-edge and near-edge spectra at the K-edge of a near-complete series of hydrated first-row transition metal ions with d-orbital occupancy ranging from d 2 to d 10 . Starting with optimized structures that were derived from an explicit solvation treatment at the density functional theory (DFT) level, we then compute the pre-edge X-ray absorption spectra at the metal ion K-edge at the time dependent density functional theory (TDDFT) and restricted active-space second-order perturbation theory (RASPT2) levels of theory. TDDFT calculations provide accurate results for spectra that are dominated by single excitations, while significant improvements were obtained with RASPT2 calculations that correctly distinguish between singly and doubly excited states, where relevant, with quantitative accuracy compared with experiment. We analyze and assign the dierent pre-edge features for each metal ion in order to reveal the impact of the variations in the d orbital occupancy on the rst-shell coordination environment. For completeness, we also report the lowest energy ligand-field d-d transitions for all the transition metal aqua ions considered in this study using complete active space second order perturbation theory (CASPT2).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗