Modeling Singlet Oxygen-Induced Degradation Pathways Including Environmental Effects of 1,2-Dimethoxyethane in Li-O2 Batteries through Density Functional Theory.

We employ density functional theory (DFT) to examine reaction mechanisms involving singlet oxygen 1Δg (1O2) and 1,2-dimethoxyethane (DME) to probe potential parasitic reactions occurring in Li-O2 batteries. First, we investigate the attack of 1O2 on the ethylene group (-CH2-CH2-) to form H2O2 and a C-C double bond in a single step. Second, we look at hydroperoxide formation that occurs via a two-step mechanism. We employ an implicit solvent model, Li+ coordination, and external electric fields to model the complex electrolyte environment near the cathode of a Li-O2 battery. The initial barriers for these reactions are decreasing functions of the dielectric constant of the implicit solvent model as well as the strength of the electric field. These initial barriers range between 17 and 26 kcal mol-1 for large dielectric constants and in the presence of electric fields. We discuss the implications of these results on ether-based electrolytes for Li-O2 batteries.

Reaction of Singlet Oxygen with the Ethylene Group: Implications for Electrolyte Stability in Li-Ion and Li-O2 Batteries.

Recent experimental and computational evidence indicates that singlet oxygen (1O2) attacks the ethylene group (-CH2-CH2-) in ethylene carbonate (EC) leading to degradation in Li-ion batteries employing EC as the electrolyte solvent [J. Phys. Chem. A 2018, 122, 8828-8839]. Here, we employ computational quantum chemistry to explore this mechanism in detail for a large set of organic molecules. Benchmark calculations comparing density functional theory to the complete active space second-order perturbation theory and internally contracted multireference configuration interaction indicate that the M11 functional adequately captures trends in the transition-state energies for this mechanism. Based on our results, we recommend that solvents which include the ethylene group should be avoided in Li-ion and Li-O2 batteries where 1O2 is generated unless neighboring functional groups raise the reaction barrier to avoid this decomposition pathway.

Heterogeneous CPU + GPU Algorithm for Variational Two-Electron Reduced-Density Matrix-Driven Complete Active-Space Self-Consistent Field Theory.

We present a heterogeneous central processing unit (CPU) + graphical processing unit (GPU) algorithm for the direct variational optimization of the two-electron reduced-density matrix (2RDM) under two-particle N-representability conditions. This variational 2RDM (v2RDM) approach is the driver for a polynomially scaling approximation to configuration-interaction-driven complete active-space self-consistent field (CASSCF) theory. For v2RDM-based CASSCF computations involving an active space consisting of 50 electrons in 50 orbitals, we observe a speedup of a factor of 3.7 when the code is executed on a combination of an NVIDIA TITAN V GPU and an Intel Core i7-6850k CPU, relative to the case when the code is executed on the CPU alone. We use this GPU-accelerated v2RDM-CASSCF algorithm to explore the electronic structure of the 3,k-circumacene and 3,k-periacene series (k = 2-7) and compare indicators of polyradical character in the lowest-energy singlet states to those observed for oligoacene molecules. The singlet states in larger circumacene and periacene molecules display the same polyradical characteristics observed in oligoacenes, with the onset of this behavior occurring at smallest k for periacenes, followed by the circumacenes and then the oligoacenes. However, the unpaired electron density that accumulates along the zigzag edge of the circumacenes is slightly less than that which accumulates in the oligoacenes, while periacenes clearly exhibit the greatest buildup of unpaired electron density in this region.

Analytic Energy Gradients for Variational Two-Electron Reduced-Density Matrix Methods within the Density Fitting Approximation.

Analytic energy gradients are presented for a variational two-electron reduced-density-matrix-driven complete active space self-consistent field (v2RDM-CASSCF) procedure that employs the density fitting (DF) approximation to the two-electron repulsion integrals. The DF approximation significantly reduces the computational cost of v2RDM-CASSCF gradient evaluation, in terms of both the number of floating-point operations and memory requirements, enabling geometry optimizations on much larger chemical systems than could previously be considered at this level of theory [ Maradzike et al., J. Chem. Theory Comput. , 2017 , 13 , 4113 - 4122 ]. The efficacy of v2RDM-CASSCF for computing equilibrium geometries and harmonic vibrational frequencies is assessed using a set of 25 small closed- and open-shell molecules. Equilibrium bond lengths from v2RDM-CASSCF differ from those obtained from configuration-interaction-driven CASSCF (CI-CASSCF) by 0.62 and 0.05 pm, depending on whether the optimal reduced-density matrices from v2RDM-CASSCF satisfy two-particle N-representability conditions (PQG) or PQG plus partial three-particle conditions (PQG+T2), respectively. Harmonic vibrational frequencies, which are obtained by finite differences of v2RDM-CASSCF analytic energy gradients, similarly demonstrate that quantitative agreement between v2RDM- and CI-CASSCF requires the consideration of partial three-particle N-representability conditions. Lastly, optimized geometries are obtained for the lowest-energy singlet and triplet states of the linear polyacene series up to dodecacene (C50H28), in which case the active space is comprised of 50 electrons in 50 orbitals. The v2RDM-CASSCF singlet-triplet energy gap extrapolated to an infinitely long linear acene molecule is found to be 7.8 kcal mol-1.

Radicals derived from acetaldehyde and vinyl alcohol.

Vinyl alcohol and acetaldehyde are isoelectronic products of incomplete butanol combustion. Along with the radicals resulting from the removal of atomic hydrogen or the hydroxyl radical, these species are studied here using ab initio methods as complete as coupled cluster theory with single, double, triple, and perturbative quadruple excitations [CCSDT(Q)], with basis sets as large as cc-pV5Z. The relative energies provided herein are further refined by including corrections for relativistic effects, the frozen core approximation, and the Born-Oppenheimer approximation. The effects of anharmonic zero-point vibrational energies are also treated. The syn conformer of vinyl alcohol is predicted to be lower in energy than the anti conformer by 1.1 kcal mol-1. The alcoholic hydrogen of syn-vinyl alcohol is found to be the easiest to remove, requiring 84.4 kcal mol-1. Five other radicals are also carefully considered, with four conformers investigated for the 1-hydroxyvinyl radical. Beyond energetics, we have conducted an overhaul of the spectroscopic literature for these species. Our results also provide predictions for fundamental modes yet to be reported experimentally. To our knowledge, the ν3 (3076 cm-1) and ν4 (2999 cm-1) C-H stretches for syn-vinyl alcohol and all but one of the vibrational modes for anti-vinyl alcohol (ν1-ν14) are yet to be observed experimentally. For the acetyl radical, ν6 (1035 cm-1), ν11 (944 cm-1), ν12 (97 cm-1), and accounting for our changes to the assignment of the 1419.9 cm-1 experimental mode, ν10 (1441 cm-1), are yet to be observed. We have predicted these unobserved fundamentals and reassigned the experimental 1419.9 cm-1 frequency in the acetyl radical to ν4 rather than to ν10. Our work also strongly supports reassignment of the ν10 and ν11 fundamentals of the vinoxy radical. We suggest that the bands assigned to the overtones of these fundamentals were in fact combination bands. Our findings may be useful in constructing improved combustion models of butanol and in spectroscopically characterizing these molecules further.

Intermolecular interactions and proton transfer in the hydrogen halide-superoxide anion complexes.

The superoxide radical anion O2(-) is involved in many important chemical processes spanning different scientific disciplines (e.g., environmental and biological sciences). Characterizing its interaction with various substrates to help elucidate its rich chemistry may have far reaching implications. Herein, we investigate the interaction between O2(-) (X[combining tilde] (2)Πg) and the hydrogen halides (X[combining tilde] (1)Σ) with coupled-cluster theory. In contrast to the short (1.324 Å) hydrogen bond formed between the HF and O2(-) monomers, a barrierless proton transfer occurs for the heavier hydrogen halides with the resulting complexes characterized as long (>1.89 Å) hydrogen bonds between halide anions and the HO2 radical. The dissociation energy with harmonic zero-point vibrational energy (ZPVE) for FHO2(-) (X[combining tilde] (2)A'') → HF (X[combining tilde] (1)Σ) + O2(-) (X[combining tilde] (2)Πg) is 31.2 kcal mol(-1). The other dissociation energies with ZPVE for X(-)HO2 (X[combining tilde] (2)A'') → X(-) (X[combining tilde] (1)Σ) + HO2 (X[combining tilde] (2)A'') are 25.7 kcal mol(-1) for X = Cl, 21.9 kcal mol(-1) for X = Br, and 17.9 kcal mol(-1) for X = I. Additionally, the heavier hydrogen halides can form weak halogen bonds H-XO2(-) (X[combining tilde] (2)A'') with interaction energies including ZPVE of -2.3 kcal mol(-1) for HCl, -8.3 kcal mol(-1) for HBr, and -16.7 kcal mol(-1) for HI.

Can density cumulant functional theory describe static correlation effects?

We evaluate the performance of density cumulant functional theory (DCT) for capturing static correlation effects. In particular, we examine systems with significant multideterminant character of the electronic wave function, such as the beryllium dimer, diatomic carbon, m-benzyne, 2,6-pyridyne, twisted ethylene, as well as the barrier for double-bond migration in cyclobutadiene. We compute molecular properties of these systems using the ODC-12 and DC-12 variants of DCT and compare these results to multireference configuration interaction and multireference coupled-cluster theories, as well as single-reference coupled-cluster theory with single, double (CCSD), and perturbative triple excitations [CCSD(T)]. For all systems the DCT methods show intermediate performance between that of CCSD and CCSD(T), with significant improvement over the former method. In particular, for the beryllium dimer, m-benzyne, and 2,6-pyridyne, the ODC-12 method along with CCSD(T) correctly predict the global minimum structures, while CCSD predictions fail qualitatively, underestimating the multireference effects. Our results suggest that the DC-12 and ODC-12 methods are capable of describing emerging static correlation effects but should be used cautiously when highly accurate results are required. Conveniently, the appearance of multireference effects in DCT can be diagnosed by analyzing the DCT natural orbital occupations, which are readily available at the end of the energy computation.

Conical intersections and low-lying electronic states of tetrafluoroethylene.

The low-lying electronic states of tetrafluoroethylene (C2 F4 ) are characterized theoretically for the first time using equation-of-motion coupled cluster theory (EOM-CCSD), and complete active space self-consistent field (CASSCF) and second-order perturbation theory (CASPT2). Computations are performed for vertical excitation energies, equilibrium geometries, minimum-energy conical intersections, and potential energy curves along three geometric coordinates: 1) twisting of the FCCF dihedral angle, 2) pyramidalization of the CF2 group, and 3) migration of a fluorine atom resulting in an ethylidene-like (CF3 CF) structure. The results suggest two relaxation pathways from the Rydberg-3s excited electronic state to the ground state. These relaxation pathways are discussed in conjunction with the femtosecond photoionization spectroscopy results of Trushin et al. [ChemPhysChem- 2004, 5, 1389].

Tetragermacyclobutadiene: energetically disfavored with respect to its structural isomers.

Germanium has been a central feature in the renaissance of main-group inorganic chemistry. Herein, we present the stationary-point geometries of tetragermacyclobutadiene and its related isomers on the singlet potential energy surface at the CCSD(T)/cc-pVTZ level of theory. Three of these 12 structures are reported for the first time and one of them is predicted to lie only 0.4 kcal mol(-1) above the previously reported global minimum. Focal-point analyses has provided electronic energies at the CCSD(T) level of theory, which are extrapolated to the complete basis-set limit and demonstrate the convergence behavior of the electronic energies with improving levels of theory and increasing basis-set size. The lowest-energy structure is the bicyclic structure, which lies 35 kcal mol(-1) below the "all-Ge" cyclobutadiene structure. The reaction energies for the association of known Ge hydrides (e.g., digermene) to form Ge4H4 indicate that Ge4H4 could be observed experimentally. We investigate the bonding patterns by examining the frontier molecular orbitals. Our results demonstrate that: 1) the cyclic isomers of (GeH)4 distort to maximize the mixing of the p orbitals that are involved in the π system of tetragermacyclobutadiene and 2) the lowest-energy isomers exhibit unusual bonding arrangements (e.g., bridging H bonds) that maximize the nonbonding electron density at the Ge centers.