Systematic Theoretical Study of Non-nuclear Electron Density Maxima in Some Diatomic Molecules

Prediction of Covalent Metal-Metal Bonding in Cp-M-M'-Nacnac Complexes of Group 2 and 12 Metals (Be, Mg, Ca, Zn, Cd, Hg)

Bimetallic CpMM'Nacnac molecules with group 2 and 12 metals (M=Be, Mg, Ca, Zn, Cd, Hg) that contain novel metal-metal bonding have been investigated in a theoretical study of their molecular and electronic structure, thermodynamic stability, and metal-metal bonding. In all cases the metal-metal bonds are characterized as electron-sharing covalent single bonds from natural bond orbital (NBO) and energy-decomposition analysis with natural orbitals of chemical valence (EDA-NOCV) analysis. The sum of [MM'] charges is relatively constant, with all complexes exhibiting a [MM']2+ core. Quantum theory of atoms in molecules (QTAIM) analysis indicates the presence of non-nuclear attractors (NNA) in the metal-metal bonds of the BeBe, MgMg, and CaCa complexes. There is substantial electron density (0.75-1.33 e) associated with the NNAs, which indicates that these metal-metal bonds, while classified as covalent electron-sharing bonds, retain significant metallic character that can be associated with reducing reactivity of the complex. The predicted stability of these complexes, combined with their novel covalent metal-metal bonding and potential as reducing agents, make them appealing targets for the synthesis of new metal-metal bonds.

On the sensitivity of computed partial charges toward basis set and (exchange-)correlation treatment

Partial charges are a central concept in general chemistry and chemical biology, yet dozens of different computational definitions exist. In prior work [Cho et al., ChemPhysChem 21, 688-696 (2020)], we showed that these can be reduced to at most three 'principal components of ionicity'. The present study addressed the dependence of computed partial charges q on 1-particle basis set and (for WFT methods) n -particle correlation treatment or (for DFT methods) exchange-correlation functional, for several representative partial charge definitions such as QTAIM, Hirshfeld, Hirshfeld-I, HLY (electrostatic), NPA, and GAPT. Our findings show that semi-empirical double hybrids can closely approach the CCSD(T) 'gold standard' for this property. In fact, owing to an error compensation in MP2, CCSD partial charges are further away from CCSD(T) than is MP2. The nonlocal correlation is important, especially when there is a substantial amount of nonlocal exchange. Employing range separation proves to be "mostly" not advantageous, while global hybrids perform optimally for 20%-30% Hartree-Fock exchange across all charge types. Basis set convergence analysis shows that an augmented triple-zeta heavy-aug-cc-pV(T+d)Z basis set or a partially augmented jun-cc-pV(T+d)Z basis set is sufficient for Hirshfeld, Hirshfeld-I, HLY, and GAPT charges. In contrast, QTAIM and NPA display slower basis set convergence. It is noteworthy that for both NPA and QTAIM, HF exhibits markedly slower basis set convergence than the correlation components of MP2 and CCSD. Triples corrections in CCSD(T), denoted as CCSD(T)-CCSD, exhibit even faster basis set convergence.

On the nature of the two-positron bond: evidence for a novel bond type

The nature of the newly proposed two-positron bond in (PsH)2, which is composed of two protons, four electrons and two positrons, is considered in this contribution. The study is done at the multi-component-Hartree-Fock (MC-HF) and the Diffusion Monte Carlo (DMC) levels of theory by comparing ab initio data, analyzing the spatial structure of the DMC wavefunction, and applying the multi-component quantum theory of atoms in molecules and the two-component interacting quantum atoms energy partitioning schemes to the MC-HF wavefunction. The analysis demonstrates that (PsH)2 to a good approximation may be conceived of as two slightly perturbed PsH atoms, bonded through a two-positron bond. In contrast to the usual two-electron bonds, the positron exchange phenomenon is quite marginal in the considered two-positron bond. The dominant stabilizing mechanism of bonding is a novel type of classical electrostatic interaction between the positrons, which are mainly localized between nuclei, and the surrounding electrons. To emphasize its uniqueness, this mechanism of bonding is proposed to be called gluonic which has also been previously identified as the main driving mechanism behind formation of the one-positron bond in [H-,e+,H-]. We conclude that the studied two-positron bond should not be classified as a covalent bond and it must be seen as a brand-new type of bond, foreign to the electronic bonding modes discovered so far in the purely electronic systems.

Dipolar 1,3-cycloaddition of thioformaldehyde S-methylide (CH2 SCH2 ) to ethylene and acetylene. A comparison with (valence) isoelectronic O3 , SO2 , CH2 OO and CH2 SO

Methods rooted in the density functional theory and in the coupled cluster ansatz were employed to investigate the cycloaddition reactions to ethylene and acetylene of 1,3-dipolar species including ozone and the derivatives issued from replacement of the central oxygen atom by the valence-isoelectronic sulfur atom, and/or of one or both terminal oxygen atoms by the isoelectronic CH₂ group. This gives rise to five different 1,3-dipolar compounds, namely ozone itself (O₃ ), sulfur dioxide (SO₂ ), the simplest Criegee intermediate (CH₂ OO), sulfine (CH₂ SO), and thioformaldehyde S-methylide (CH₂ SCH₂ , TSM). The experimental and accurate theoretical data available for some of those molecules were employed to assess the accuracy of two last-generation composite methods employing conventional or explicitly correlated post-Hartree-Fock contributions (jun-Cheap and SVECV-f12, respectively), which were then applied to investigate the reactivity of TSM. The energy barriers provided by both composite methods are very close (the average values for the two composite methods are 7.1 and 8.3 kcal mol^-1 for the addition to ethylene and acetylene, respectively) and comparable to those ruling the corresponding additions of ozone (4.0 and 7.7 kcal mol^-1 , respectively). These and other evidences strongly suggest that, at least in the case of cycloadditions, the reactivity of TSM is similar to that of O₃ and very different from that of SO₂ .

New venues in electron density analysis

We provide a comprehensive overview of the chemical information within the electron density: how to extract information, but also how to obtain and how to assess the quality of the...

Model systems for screening and investigation of lithium metal electrode chemistry and dendrite formation

The use of lithium metal as an electrode for electrochemical energy storage will provide a significant impact on practical energy storage technology. Unfortunately, the use of lithium metal is plagued with challenging chemical problems. Specifically, the formation of a solid electrolyte interphase layer and the nucleation and growth of lithium dendrites: both must be addressed and controlled in order to achieve a practically useable pure lithium metal electrode. Currently sophisticated experimental techniques and computationally expensive simulations are being used to probe these problems but these methods are arduous and time consuming which delays timely evaluation and insight into the rapidly changing field of advanced energy storage. Here, we report the use of DFT simulations of lithium nanoclusters to investigate and explore lithium metal chemistry with inexpensive computational methods to gain greater insight into electrochemical reductions and the nucleation and growth of dendrites. DME, LiTFSI, and LiFSI reduction energetics and structures with electrode effects from lithium metal are reported providing better physical description of the absolute reduction potential characterization. The electronic structure of the lithium nanoclusters were used to investigate the nucleation and growth of lithium dendrites from an ab initio perspective. The results demonstrate that kinetic processes have more control over non uniform growth than thermodynamic processes. Based on this information, a non ab initio model was created in Matlab that shows the initial stages of dendrite nucleation considering approximately 2000 atoms.

Core Electron Topologies in Chemical Compounds: Case Study of Carbon versus Silicon

The similarities and differences between carbon and silicon have attracted the curiosity of chemists for centuries. Similarities and analogies can be found in their saturated compounds, but carbon exhibits a cornucopia of unsaturated compounds that silicon (and most other elements) cannot replicate. While this qualitative difference is empirically well known, quantum chemistry has previously only described quantitative differences related to orbital overlap, steric effects, or orbital energies. We study C₂ and Si₂ and their hydrides X₂ H_2n (X=C, Si; n=1, 2, 3) by first-principles quantum chemical calculation, and find a qualitative difference in the topologies of the core electrons: carbon has the propensity to alter its core electron topology when forming unsaturated compounds, and silicon has not. We draw a connection between the core electron topologies and ionization energies, and identify other elements we expect to have similarly flexible core topologies as carbon.

Computational Investigation on the Role of Disilene Substituents Toward N2O Activation

The effect of substituents in disilene mediated N₂O activation was studied at the M06-2X/QZVP//ωB97xD/TZVP level of theory. The relationship between structural diversity and the corresponding reactivity of six disilenes (I_A-F^t) in the presence of four different substituents (-NMe₂, -Cl, -Me, -SiMe₃) is addressed in this investigation. We primarily propose two plausible mechanistic routes: Pathway I featuring disilene → silylene decomposition followed by N₂O coordination and Pathway II constituting the N₂O attack without Si-Si bond cleavage. Depending on the fashion of N₂O approach the latter route was further differentiated into Pathway IIa and Pathway IIb detailing the "end-on" and "side-on" attack to the disilene scaffold. Interestingly, the lone pair containing substituents (-NMe₂, -Cl,) facilitates disilene → silylene dissociation; on the contrary it reduces the electrophilicity at Si center in silylene, a feature manifested with higher activation barrier during N₂O attack. In the absence of any lone-pair influence from substituents (-Me, -SiMe₃), the decomposition of disilenes is considerably endothermic. Therefore, Pathway I appears to be the less preferred route for both types of substituents. In Pathway IIa, the N₂O moiety uniformly approaches via O-end to both the silicon centers in disilenes. However, the calculations reveal that Pathway IIa, although not operational for all disilenes, is unlikely to be a viable route due to the predominantly higher transition barrier (ca. 36 kcal/mol). The most feasible route in this current study accompanying moderately low activation barriers (∼19-26 kcal/mol) is Pathway IIb, which involves successive addition of two N₂O units proceeding via terminal N, O toward the Si centers and is applicable for all disilenes. The reactivity of substituted disilenes can be estimated in terms of the first activation barrier of N₂O attack. Surprisingly, in Pathway IIb, the initial activation barrier and hence the reactivity shows negligible correlation with Si-Si bond strength, indicating toward the versatility of the reaction route.

Atomic charge transfer-counter polarization effects determine infrared CH intensities of hydrocarbons: a quantum theory of atoms in molecules model

Atomic charge transfer-counter polarization effects determine most of the infrared fundamental CH intensities of simple hydrocarbons, methane, ethylene, ethane, propyne, cyclopropane and allene. The quantum theory of atoms in molecules/charge-charge flux-dipole flux model predicted the values of 30 CH intensities ranging from 0 to 123 km mol(-1) with a root mean square (rms) error of only 4.2 km mol(-1) without including a specific equilibrium atomic charge term. Sums of the contributions from terms involving charge flux and/or dipole flux averaged 20.3 km mol(-1), about ten times larger than the average charge contribution of 2.0 km mol(-1). The only notable exceptions are the CH stretching and bending intensities of acetylene and two of the propyne vibrations for hydrogens bound to sp hybridized carbon atoms. Calculations were carried out at four quantum levels, MP2/6-311++G(3d,3p), MP2/cc-pVTZ, QCISD/6-311++G(3d,3p) and QCISD/cc-pVTZ. The results calculated at the QCISD level are the most accurate among the four with root mean square errors of 4.7 and 5.0 km mol(-1) for the 6-311++G(3d,3p) and cc-pVTZ basis sets. These values are close to the estimated aggregate experimental error of the hydrocarbon intensities, 4.0 km mol(-1). The atomic charge transfer-counter polarization effect is much larger than the charge effect for the results of all four quantum levels. Charge transfer-counter polarization effects are expected to also be important in vibrations of more polar molecules for which equilibrium charge contributions can be large.