Molecular core-valence correlation effects involving the post-d elements Ga–Rn: Benchmarks and new pseudopotential-based correlation consistent basis sets

Low-lying states of MX2 (M = Ag, Au; X = Cl, Br and I) with coupled-cluster approaches: effect of the basis set, high level correlation and spin–orbit coupling

Spin–orbit coupling, electron correlation level and basis set are important in describing Renner–Teller and pseudo-Jahn–Teller effects and properties of MX₂.

Balancing Donor-Acceptor and Dispersion Effects in Heavy Main Group Element π Interactions: Effect of Substituents on the Pnictogen⋅⋅⋅π Arene Interaction

High-level ab initio calculations using the DLPNO-CCSD(T) method in conjunction with the local energy decomposition (LED) were performed to investigate the nature of the intermolecular interaction in bismuth trichloride adducts with π arene systems. Special emphasis was put on the effect of substituents in the aromatic ring. For this purpose, benzene derivatives with one or three substituents (R=NO2 , CF3 , OCHO, OH, and NH2 ) were chosen and their influence on donor-acceptor interaction as well as on the overall interaction strength was examined. Local energy decomposition was performed to gain deeper insight into the composition of the interaction. Additionally, the study was extended to the intermolecular adducts of arsenic and antimony trichloride with benzene derivatives having one substituent (R=NO2 and NH2 ) in order to rationalize trends in the periodic table. The analysis of natural charges and frontier molecular orbitals shows that donor-acceptor interactions are of π→σ* type and that their strength correlates with charge transfer and orbital energy differences. An analysis of different bonding motifs (Bi⋅⋅⋅π arene, Bi⋅⋅⋅R, and Cl⋅⋅⋅π arene) shows that if dispersion and donor-acceptor interaction coincide as the donor highest occupied molecular orbital (HOMO) of the arene is delocalized over the π system, the M⋅⋅⋅π arene motif is preferred. If the donor HOMO is localized on the substituent, R⋅⋅⋅π arene bonding motifs are preferred. The Cl⋅⋅⋅π arene bonding motif is the least favorable with the lowest overall interaction energy.

Beyond chemical accuracy in the heavy p-block: The first ionization potentials and electron affinities of Ga-Kr, In-Xe, and Tl-Rn

A relativistic coupled-cluster version of the Feller-Peterson-Dixon composite method has been used to accurately calculate the first ionization potentials (IPs) and electron affinities (EAs) of the post-d, p-block elements Ga-Rn. Complete basis set extrapolations including outer-core correlation at the CCSD(T) level of theory were combined with contributions from higher order electron correlation up to CCSDTQ, quantum electrodynamic effects (Lamb shift), and spin-orbit (SO) coupling including the Gaunt contribution. Several methods for including SO were investigated, in which all involved the four-component (4c) Dirac-Coulomb (DC) Hamiltonian. The treatment of SO coupling was the contribution that limited the final accuracy of the present results. In the cases where 4c-DC-CCSD(T) could be reliably used for the SO contributions, the final composite IPs and EAs agreed with the available experimental values to within an unsigned average error of just 0.16 and 0.20 kcal/mol, respectively. In all cases, the final IPs and EAs were within 1 kcal/mol of the available experimental values, except for the EAs of the group 13 elements (Ga, In, and Tl), where the currently accepted experimental values appear to be too large by as much as 4 kcal/mol. The values predicted in this work, which have estimated uncertainties of ±0.5 kcal/mol, are 5.25 (Ga), 7.69 (In), and 7.39 (Tl) kcal/mol. For the EAs of Po and At, which do not have experimental values, the current calculations predict values of 34.2 and 55.8 kcal/mol with estimated uncertainties of ±0.6 and ±0.3 kcal/mol, respectively.

The CUAGAU Set of Coupled-Cluster Reference Data for Small Copper, Silver, and Gold Compounds and Assessment of DFT Methods

We have obtained benchmark data for a set of small molecular systems of Cu, Ag, and Au using coupled-cluster methods. Using this collection of reference data (that we termed the CUAGAU set) for assessing DFT-type methods, we find the MN15-L nonhybrid DFT to be cost-effective for geometry optimization [mean absolute deviation (MAD) in bond length = 0.20 Å], with an accuracy that is comparable to that for the double-hybrid (DH) DFT method DSD-PBEP86 (MAD = 0.19 Å). For the computation of thermochemical properties, among "conventional" (non-MP2-based) DFT methods, the best performance is found for the global-hybrid meta-GGA functional MN15, with an MAD of 11.4 kJ mol^-1. We also find the nonhybrid method B97M-rV to have a reasonable performance (MAD = 14.4 kJ mol^-1), and it may serve as a cost-effective means for qualitative study. If we look beyond conventional functionals, we find DSD-PBEP86 (MAD = 7.3 kJ mol^-1) to be more accurate than even MN15. Nonetheless, this level of accuracy is still not sufficient for quantitative studies. In this regard, high-level wave function methods such as composite procedures that are based on coupled cluster are still indispensable for obtaining reliable reference data for transition-metal species.

Electronic Structure Predictions of the Energetic Properties of Tellurium Fluorides

The heats of formation, bond dissociation energies (BDEs), fluoride affinities (FA), fluorocation affinities (FCA), electron affinities (EA), and ionization energies (IP) of TeF _n ( n = 1-6) have been predicted using the Feller-Peterson-Dixon (FPD) approach. To benchmark the approach, the bond dissociation energies of Te₂ and TeO, the heats of formation of Te₂, TeH₂, TeO, and TeO₂, and the electron affinity for TeO and TeO₂ were calculated as there are experimental thermodynamic data available for these tellurium compounds, which allow confirmation of the heat of formation of Te gas as Δ H_f,0K(Te) = 50.7 ± 0.6 kcal/mol. Spin-orbit corrections are required for good results and cannot be ignored. A comparison among fluoride affinities, fluorocation affinities, electron affinities, and ionization energies of TeF _n and SeF _n is reported.

Gas Phase Silver Thermochemistry from First Principles

Domain-based local pair natural orbital coupled cluster approach with single, double, and perturbative triple excitations, DLPNO-CCSD(T), has been applied within a framework of a reduced version of the reaction-based Feller-Peterson-Dixon (FPD) scheme to predict gas phase heats of formation and absolute entropies of silver inorganic and organometallic compounds. First, we evaluated all existing experimental data currently limited by thermodynamic functions of 10 silver substances (AgH, AgF, AgBr, AgI, Ag₂, Ag₂S, Ag₂Se, Ag₂Te, AgCN, AgPO₂). The mean average deviation between computed and experimental heats of formation was found to be 1.9 kcal/mol. Notably, all predicted heats of formation turned out to be within the error bounds of their experimental counterparts. Second, we predicted heats of formation and entropies for additional 90 silver species with no experimental data available, substantially enriching silver thermochemistry. Combination of gas phase heats of formation Δ H_f and entropies S° of AgNO₂, AgSCN, Ag₂SO₄, and Ag₂SeO₄ obtained in this work, with respective solid-state information, resulted in accurate sublimation thermochemistry of these compounds. Complementation of predicted Δ H_f with heats of formation of some neutrals and positive ions produced 33 silver bond strengths of high reliability. Obtained thermochemical data are promising for developing the concepts of silver chemistry. In addition, derived heats of formation and bond dissociation enthalpies, due to their high diversity, are found to be relevant for testing and training of computational chemistry methods.

Laser-cooling with an intermediate electronic state: Theoretical prediction on bismuth hydride

The possibility of laser cooling of bismuth hydride (BiH) molecules has been investigated based on high-level ab initio calculations by considering the core-valence and the spin-orbit coupling (SOC) effects. The potential energy curves of the 12 Λ-S states as well as the 25 Ω states that split from them via SOC are obtained by multireference configuration interaction plus the Davidson correction. The properties of b-X transition are investigated. Based on our calculations, we show that the transition between Ω states b0⁺-X₁0⁺ of BiH is a possible candidate for laser cooling, with consideration of the intermediate Ω state X₂1. An optical cycling scheme is proposed by utilizing four lasers at wavelengths around 471 and 601 nm with 5400 cycles for photon absorption/emission and a sub-microkelvin temperature. Our study should shed some light on searching for possible molecular candidates for laser cooling with the existence of an intermediate electronic state.

Real-time observation of X-ray-induced intramolecular and interatomic electronic decay in CH2I2

The increasing availability of X-ray free-electron lasers (XFELs) has catalyzed the development of single-object structural determination and of structural dynamics tracking in real-time. Disentangling the molecular-level reactions triggered by the interaction with an XFEL pulse is a fundamental step towards developing such applications. Here we report real-time observations of XFEL-induced electronic decay via short-lived transient electronic states in the diiodomethane molecule, using a femtosecond near-infrared probe laser. We determine the lifetimes of the transient states populated during the XFEL-induced Auger cascades and find that multiply charged iodine ions are issued from short-lived (∼20 fs) transient states, whereas the singly charged ones originate from significantly longer-lived states (∼100 fs). We identify the mechanisms behind these different time scales: contrary to the short-lived transient states which relax by molecular Auger decay, the long-lived ones decay by an interatomic Coulombic decay between two iodine atoms, during the molecular fragmentation.

Understanding strong X-ray induced phenomena is important for applications of X-ray free-electron laser imaging. Here, the authors show time-resolved measurements of X-ray free-electron laser induced electronic decay of CH₂I₂ molecule probed with NIR pulses and identify mechanisms behind different transient states lifetimes.

Ab initio interaction potentials of the Ba, Ba+ complexes with Ar, Kr, and Xe in the lowest excited states

The complexes of the Ba atom and Ba⁺ cation with the rare gas atoms Ar, Kr, and Xe in the states associated with the 6s → 5d, 6p excitations are investigated by means of the multireference configuration interaction techniques. Scalar relativistic potentials are obtained by the complete basis limit extrapolation through the sequence of aug-cc-pwCVnZ basis sets with the cardinal numbers n = Q, T, 5, combined with the suitable effective core potentials and benchmarked against the coupled cluster with singles, doubles, and non-iterative triples calculations and the literature data available for selected electronic states. Spin-orbit coupling is taken into account by means of the state-interacting multireference configuration interaction calculations performed for the Breit-Pauli spin-orbit Hamiltonian. The results show weak spin-orbit coupling between the states belonging to distinct atomic multiplets. General trends in the interaction strength and long-range anisotropy along the rare gas series are discussed. Vibronic spectra of the Ba and Ba⁺ complexes in the vicinity of the ¹S → ¹P° and ²S → ²P° atomic transitions and diffusion cross sections of the Ba(¹S₀, ³D_J) atom in high-temperature rare gases are calculated. Comparison with available experimental data shows that multireference calculations tend to underestimate the interaction strength for excited complexes.

Theoretical investigation on the low-lying electronic states of beryllium antimonide

The Λ-S electronic states with respect to the lowest four dissociation limits of BeSb are investigated theoretically on the icMRCI + Q level employing basis set of quintuple-ζ quality. The geometrical parameters, potential energy curves, vibrational energy levels, spectroscopic constants for the twelve Λ - S states are obtained, analyzed and compared with those of the Beryllium-VA group diatomic family species where data are available. The permanent dipole moments, transition dipole moments, Einstein emission coefficients, radiative lifetimes and Franck-Condon factors for interested Λ - S states are also derived. Further assessments of the spin-orbit coupling effect are performed for states associated with the first two dissociation asymptotes of BeSb. Four Λ - S states split into seven Ω states, and some of the PECs are distorted significantly through the spin-orbit coupling effect, which is similar to its isovalent diatomics BeAs. In consideration of potential risks of manipulating beryllium-containing species directly, the information associated with molecular structures, spectroscopic parameters as well as transition properties that provide in this paper is anticipated to serve as guidelines for further researches of BeSb.

Detailed benchmark ab initio mapping of the potential energy surfaces of the X + C2H6 [X = F, Cl, Br, I] reactions

We provide benchmark relative energies for the stationary points of three different channels of the halogen atom + ethane reactions.

Dispersion interactions between neighboring Bi atoms in (BiH3 )2 and Te(BiR2 )2

Triggered by the observation of a short Bi⋯Bi distance and a BiTeBi bond angle of only 86.6° in the crystal structure of bis(diethylbismuthanyl)tellurane quantum chemical computations on interactions between neighboring Bi atoms in Te(BiR₂ )₂ molecules (R = H, Me, Et) and in (BiH₃ )₂ were undertaken. Bi⋯Bi distances atoms were found to significantly shorten upon inclusion of the d shells of the heavy metal atoms into the electron correlation treatment, and it was confirmed that interaction energies from spin component-scaled second-order Møller-Plesset theory (SCS-MP2) agree well with coupled-cluster singles and doubles theory including perturbative triples (CCSD(T)). Density functional theory-based symmetry-adapted perturbation theory (DFT-SAPT) was used to study the anisotropy of the interplay of dispersion attraction and steric repulsion between the Bi atoms. Finally, geometries and relative stabilities of syn-syn and syn-anti conformers of Te(BiR₂ )₂ (R = H, Me, Et) and interconversion barriers between them were computed. © 2018 Wiley Periodicals, Inc.

Interaction potentials and transport properties of Ba, Ba+, and Ba2+ in rare gases from He to Xe

A highly accurate, consistent set of ab initio interaction potentials is obtained for the title systems at the coupled cluster with singles, doubles, and non-iterative triples level of theory with extrapolation to the complete basis set limit. These potentials are shown to be more reliable than the previous potentials based on their long-range behavior, equilibrium properties, collision cross sections, and transport properties.

Pair natural orbital and canonical coupled cluster reaction enthalpies involving light to heavy alkali and alkaline earth metals: the importance of sub-valence correlation

In this work, we tested canonical and domain based pair natural orbital coupled cluster methods (CCSD(T) and DLPNO-CCSD(T), respectively) for a set of 32 ligand exchange and association/dissociation reaction enthalpies involving ionic complexes of Li, Be, Na, Mg, Ca, Sr, Ba and Pb(ii). Two strategies were investigated: in the former, only valence electrons were included in the correlation treatment, giving rise to the computationally very efficient FC (frozen core) approach; in the latter, all non-ECP electrons were included in the correlation treatment, giving rise to the AE (all electron) approach. Apart from reactions involving Li and Be, the FC approach resulted in non-homogeneous performance. The FC approach leads to very small errors (<2 kcal mol^-1) for some reactions of Na, Mg, Ca, Sr, Ba and Pb, while for a few reactions of Ca and Ba deviations up to 40 kcal mol^-1 have been obtained. Large errors are both due to artificial mixing of the core (sub-valence) orbitals of metals and the valence orbitals of oxygen and halogens in the molecular orbitals treated as core, and due to neglecting core-core and core-valence correlation effects. These large errors are reduced to a few kcal mol^-1 if the AE approach is used or the sub-valence orbitals of metals are included in the correlation treatment. On the technical side, the CCSD(T) and DLPNO-CCSD(T) results differ by a fraction of kcal mol^-1, indicating the latter method as the perfect choice when the CPU efficiency is essential. For completely black-box applications, as requested in catalysis or thermochemical calculations, we recommend the DLPNO-CCSD(T) method with all electrons that are not covered by effective core potentials included in the correlation treatment and correlation-consistent polarized core valence basis sets of cc-pwCVQZ(-PP) quality.

The bismuth tetramer Bi4: the ν3 key to experimental observation

The spectroscopic identification of Bi₄ has been very elusive. Two constitutional Bi₄ isomers of T_d and C_2v symmetry are investigated and each is found to be a local energetic minimum.

Grotthuss Transport of Iodide in EMIM/I3 Ionic Crystal

Highly ionic environments can mediate unusual chemical reactions that would otherwise be considered impossible based on chemical intuition. For example, the formation of a chemical bond between two iodide anions to form a divalent polyiodide anion is seemingly prohibited due to Coulombic repulsion. Using ab initio molecular dynamics simulations, we show that in the 1-ethyl-3-methylimidazolium (EMIM)/I₃ ionic crystal, the reactive formation of divalent and even trivalent polyiodide anions occurs with extremely small energetic barriers, due to the electrostatic field of the ionic lattice. A practical consequence of this anomalous reactivity is that iodide anions are efficiently transported within the crystal through a "Grotthuss-exchange" mechanism involving bond-breaking and forming events. We characterize two distinct transport pathways, involving both I₄^2- and I₇^3- intermediates, with fast transport of iodide resulting from the release of an I^- anion on the opposite side of the intermediate species from the initial bond formation. The ordered cation arrangement in the crystal provides the necessary electrostatic screening for close approach of anions, suggesting a new counterintuitive approach to obtain high ionic conductivity. This new design principle could be used to develop better solid-state electrolytes for batteries, fuel cells, and supercapacitors.

Ground-State Gas-Phase Structures of Inorganic Molecules Predicted by Density Functional Theory Methods

We tested a battery of density functional theory (DFT) methods ranging from generalized gradient approximation (GGA) via meta-GGA to hybrid meta-GGA schemes as well as Møller-Plesset perturbation theory of the second order and a single and double excitation coupled-cluster (CCSD) theory for their ability to reproduce accurate gas-phase structures of di- and triatomic molecules derived from microwave spectroscopy. We obtained the most accurate molecular structures using the hybrid and hybrid meta-GGA approximations with B3PW91, APF, TPSSh, mPW1PW91, PBE0, mPW1PBE, B972, and B98 functionals, resulting in lowest errors. We recommend using these methods to predict accurate three-dimensional structures of inorganic molecules when intramolecular dispersion interactions play an insignificant role. The structures that the CCSD method predicts are of similar quality although at considerably larger computational cost. The structures that GGA and meta-GGA schemes predict are less accurate with the largest absolute errors detected with BLYP and M11-L, suggesting that these methods should not be used if accurate three-dimensional molecular structures are required. Because of numerical problems related to the integration of the exchange-correlation part of the functional and large scattering of errors, most of the Minnesota models tested, particularly MN12-L, M11, M06-L, SOGGA11, and VSXC, are also not recommended for geometry optimization. When maintaining a low computational budget is essential, the nonseparable gradient functional N12 might work within an acceptable range of error. As expected, the DFT-D3 dispersion correction had a negligible effect on the internuclear distances when combined with the functionals tested on nonweakly bonded di- and triatomic inorganic molecules. By contrast, the dispersion correction for the APF-D functional has been found to shorten the bonds significantly, up to 0.064 Å (AgI), in Ag halides, BaO, BaS, BaF, BaCl, Cu halides, and Li and Na halides and hydrides. These results do not agree well with very accurate structures derived from microwave spectroscopy; we therefore believe that the dispersion correction in the APF-D method should be reconsidered. Finally, we found that inaccurate structures can easily lead to errors of few kcal/mol in single-point energies.

Accurate Vibrational-Rotational Parameters and Infrared Intensities of 1-Bromo-1-fluoroethene: A Joint Experimental Analysis and Ab Initio Study

The medium-resolution gas-phase infrared (IR) spectra of 1-bromo-1-fluoroethene (BrFC═CH₂, 1,1-C₂H₂BrF) were investigated in the range 300-6500 cm^-1, and the vibrational analysis led to the assignment of all fundamentals as well as many overtone and combination bands up to three quanta, thus giving an accurate description of its vibrational structure. Integrated band intensity data were determined with high precision from the measurements of their corresponding absorption cross sections. The vibrational analysis was supported by high-level ab initio investigations. CCSD(T) computations accounting for extrapolation to the complete basis set and core correlation effects were employed to accurately determine the molecular structure and harmonic force field. The latter was then coupled to B2PLYP and MP2 computations in order to account for mechanical and electrical anharmonicities. Second-order perturbative vibrational theory was then applied to the thus obtained hybrid force fields to support the experimental assignment of the IR spectra.

Reactivity, vibrational spectroscopy, internal rotation and thermochemical aspects of methylarsine

The aim of this investigation was to perform a characterization of the spectroscopic and thermodynamic properties of methylarsine (CH₃AsH₂). Post-Hartree-Fock, 29 DFT methods and eight different composite methodologies were employed in these analyses. A comparison between harmonic and anharmonic frequency accuracies in reproducing the observable frequencies was performed here. In addition, the CH₃AsH₂→CH₂AsH₃ isomerization barrier energy was estimated in 100kcalmol^-1, whereas the H₂-release routes barrier heights were in the 45-107kcalmol^-1 range. A rate constant of 10^-66s^-1 was predicted regarding the isomerization route, while the CH₂AsH₃ hydrogen elimination mechanism is faster than the methylarsine one. The transition state structure of the CH₃AsH₂ internal rotational barrier energy varied between 1.0 and 1.4kcalmol^-1. For the CH₂AsH₃ internal rotation the estimated barrier heights varied 0.6-2.5kcalmol^-1. The adiabatic ionization energy and the heat of formation each structure was also calculated here.

A look at the density functional theory zoo with the advanced GMTKN55 database for general main group thermochemistry, kinetics and noncovalent interactions

We present the updated and extended GMTKN55 benchmark database for more accurate and extensive energetic evaluation of density functionals and other electronic structure methods with detailed guidelines for method users.

Combining the spin-separated exact two-component relativistic Hamiltonian with the equation-of-motion coupled-cluster method for the treatment of spin–orbit splittings of light and heavy elements

An accurate and efficient treatment of spin–orbit splittings has been achieved by combining the sf-X2C+soc-DKH1 Hamiltonian with the equation-of-motion coupled-cluster method.

Experimental observation and quantum chemical investigation of thallium(i) (Z)-methanediazotate: synthesis of a long sought and highly reactive species

Thallium(i) (Z)-methanediazotate has been synthesised through TlOPr induced cleavage of ¹⁵N-labelled N-methyl-N-nitrosourea, and confirmed by NMR analysis. CCSD(T)/CBS calculations reveal a kinetic product control.

Role(s) of TMA in polymerization

Quenched-flow data for propene polymerization with rac-Me2Si(2-Me-4-Ph-1-indenyl)2ZrCl2/MAO support a picture where removal of MAO qualitatively changes the kinetic profile from a mainly enthalpic to a mainly entropic barrier. DFT studies suggest that a not previously recognized singly-bridged end-on coordination mode of Me6Al2 to catalytically active centers may be kinetically relevant as a resting state. In contrast, the more traditional doubly-bridged complex of Me3Al is proposed to be more relevant to chain transfer to cocatalyst.

Equation-of-Motion Coupled-Cluster Theory for Excitation Energies of Closed-Shell Systems with Spin-Orbit Coupling

Excitation energies of closed-shell systems based on the equation-of-motion (EOM) coupled-cluster theory at the singles and doubles (CCSD) level with spin-orbit coupling (SOC) included in the post-Hartree-Fock treatment are implemented in the present work. SOC can be included in both the CC and EOM steps (EOM-SOC-CCSD) or only in the EOM part (SOC-EOM-CCSD). The latter approach is an economical way to account for SOC effects, but excitation energies with this approach are not size-intensive. When the unlinked term in the latter approach is neglected (cSOC-EOM-CCSD), size-intensive excitation energies can be obtained. Time-reversal symmetry and spatial symmetry are exploited to reduce the computational effort. Imposing time-reversal symmetry results in a real matrix representation for the similarity-transformed Hamiltonian, which facilitates the requirement of time-reversal symmetry for new trial vectors in Davidson's algorithm. Results on some closed-shell atoms and molecules containing heavy elements show that EOM-SOC-CCSD can provide excitation energies and spin-orbit splittings with reasonable accuracy. On the other hand, the SOC-EOM-CCSD approach is able to afford accurate estimates of SOC effects for valence electrons of systems containing elements up to the fifth row, while cSOC-EOM-CCSD is less accurate for spin-orbit splittings of transitions involving p1/2 spinors, even for Kr.

Rotating Ring-Disk Electrode and Quantum-Chemical Study of the Electrochemical Reduction of Monoiodoacetic Acid and Iodoform

This study examined the electrochemical (EC) reduction of monoiodoacetic acid (MIAA) and iodoform (CHI3), which are typical iodine-containing disinfection byproducts (I-DBPs). Experiments carried out using the method of a rotating ring-disk electrode (RRDE) with a gold working electrode showed that the reduction of CHI3 and MIAA is diffusion-controlled. The MIAA diffusion coefficient was determined to be (1.86 ± 0.24)·10(-5) cm(2) s(-1). The yield of the iodide ion formed as a result of MIAA or CHI3 reduction was affected by the presence of dissolved organic matter (DOM) and resorcinol. Increasing concentrations of DOM or resorcinol did not affect the EC reduction of the examined I-DBPs, but the formation of iodide was suppressed. This indicated that free iodine, ·I, was formed as a result of the first step in the EC reduction of MIAA and CHI3. This also indicated that the pathway of the EC reduction of MIAA and CHI3 was different from that typical for the reduction of Br- and Cl-containing DBPs, in which case Br(-) or Cl(-) tend to be formed as a result of the electron transfer. Quantum-chemical (QC) calculations confirmed the thermodynamic likelihood of and possible preference to the formation of free iodine species as a result of the EC reduction of MIAA, CHI3, and other I-DBPs.

A Surprisingly Simple Electrostatic Model Explains Bent Versus Linear Structures in M(+)-RG2 Species (M = Group 1 Metal, Li-Fr; RG = Rare Gas, He-Rn)

It is found that a simple electrostatic model involving competition between the attractive dispersive interaction and induced-dipole repulsion between the two RG atoms performs extremely well in rationalizing the M(+)-RG2 geometries, where M = group 1 metal and RG = rare gas. The Li(+)-RG2 and Na(+)-RG2 complexes have previously been found to exhibit quasilinear or linear minimum-energy geometries, with the Na(+)-RG2 complexes having an additional bent local minimum [A. Andrejeva, A. M. Gardner, J. B. Graneek, R. J. Plowright, W. H. Breckenridge, T. G. Wright, J. Phys. Chem. A, 2013, 117, 13578]. In the present work, the geometries for M = K-Fr are found to be bent. A simple electrostatic model explains these conclusions and is able to account almost quantitatively for the binding energy of the second RG atom, as well as the form of the angular potential, for all 36 titular species. Additionally, results of population analyses are presented together with orbital contour plots; combined with the success of the electrostatic model, the expectation that these complexes are all physically bound is confirmed.