Radical Bonding: Structure and Stability of Bis(Phenalenyl) Complexes of Divalent Metals from across the Periodic Table

We examine the bonding possibilities of the bis(phenalenyl) MP2 sandwich complexes of the divalent metals M = Be, Mg, Ca, Sr, Ba, Zn, Cd, and Hg, at the B3LYP level of theory. The outcome is an extraordinarily diverse class of low symmetry bis(phenalenyl)metal complexes in which bonding preferences and binding enthalpies differ dramatically. The lowest energy group 2 metal MP2 complexes include an intriguing η1,η3 BeP2 structure, and bent η6,η6 systems for M = Ca, Sr, and Ba. The group 12 bis(phenalenyl) complexes are thermodynamically unstable η1,η1 slip-sandwich structures. To better understand changes in the structural preferences going from the (η6,η6) group 2 to the (η1,η1) group 12 complexes, we explored the bonding in the bis(phenalenyl) complexes of transition metals with stable +2 oxidations states between Ca and Zn in period 4. The computed binding enthalpies are large and negative for nearly all of the minimum energy bis(phenalenyl) complexes of the group 2 and the transition metals; they are tiny for MgP2, and are quite positive for the group 12 systems. The structural preferences and stability of the complexes is a subtle negotiation of several influences: the (un)availability of (n - 1)d and np, orbitals for bonding, the cost of the rehybridization at carbon sites in the phenalenyl rings in preparation for bonding to the metals, and the (P—P) interaction between the phenalenyl radicals

Bending Ternary Dihalides: A Single Functional Form For Linearization Energies

IntroductionAlthough the bonding in the symmetric groups 2 and 12 dihalides (MX2) has been studied extensively1,2, remarkably little work –experimental or theoretical – has been done on the mixed (ternary) dihalides, MXY. Previously, a criterion3,4 based on atomic softness (σ) was proposed for the bending of MX2 and MXY molecules. We extend this softness criterion on the slate of the mixed dihalides and the predicted separation is achieved between the bent and linear structures with almost the same cutoff, and with quasilinear species straddling the boundary. In this work, we report a complete assessment of the bonding preferences and vibrational frequencies of the mixed dihalides of the groups 2 and 12 metals MX2 and MXY

First principles predictions of van der Waals bonded inorganic crystal structures: Test case, HgCl2

We study the crystals structure and stability of four possible polymorphs of HgCl2 using first principles density functional theory. Mercury (II) halides are a unique class of materials which, depending on the halide species, form in a wide range of crystal structures, ranging from densely packed solids to layered materials and molecular solids. Predicting the groundstate structure of any member of this group from first principles, therefore, requires a general purpose functional that treats van der Waals bonding and covalent/ionic bonding adequately. Here, we demonstrate that the non-local van der Waals density functional paired with the C09 exchange functional meets this bar for HgCl2. In particular, this functional is able to predict the correct groundstate among the structures tested as well as having extremely good agreement with the experimentally known crystal structure. These results highlight the maturity of this functional and open the door to using this method for truly first principles crystal structure predictions

Tuning σ-Holes: Charge Redistribution in the Heavy (Group 14) Analogues of Simple and Mixed Halomethanes Can Impose Strong Propensities for Halogen Bonding

Halogen bonding between halide sites (in substituted organic molecules or inorganic halides) and Lewis bases is a rapidly progressing area of exploration. Investigations of this phenomenon have improved our understanding of weak intermolecular interactions and suggested new possibilities in supramolecular chemistry and crystal engineering. The capacity for halogen bonding is investigated at the MP2(full) level of theory for 100 compounds, including all 80 MH4-nXn systems (M = C, Si, Ge, Sn, and Pb; X = F, Cl, Br, and I). The charge redistribution in these molecules and the (in)stability of the σ-hole at X as a function of M and n are catalogued and examined. For the mixed MH3-mFmI compounds, we identify a complicated dependence of the relative halogen bond strengths on M and m. For m = 0, for example, the H3C-I----NH3 halogen bond is 6.6 times stronger than the H3Pb-I----NH3 bond. When m = 3, however, the F3Pb-I----NH3 bond is shorter and ∼1.6 times stronger than the F3C-I----NH3 bond. This substituent-induced reversal in the relative strengths of halogen bond energies is explained

Charge Saturation and Neutral Substitutions in Halomethanes and Their Group 14 Analogues

A computational analysis of the charge distribution in halomethanes and their heavy analogues (MH4-nXn: M= C, Si, Ge, Sn, Pb; X = F, Cl, Br, I) as a function of n uncovers a previously unidentified saturation limit for fluorides when M ≠ C. We examine the electron densities obtained at the CCSD, MP2(full), B3PW91, and HF levels of theory for 80 molecules for four different basis sets. A previously observed substituent independent charge at F in fluoromethanes is shown to be a move toward saturation that is restricted by the low polarizability of C. This limitation fades into irrelevance for the more polarizable M central atoms such that a genuine F saturation is realized in those cases. A conceptual model leads to a function of the form [qM(n) -- qM(n)] = a[χA\u27 -- χA] + b that links the electronegativities (χ) of incoming and leaving atoms (e.g., A\u27 = X and A = H for the halogenation of MH4-nXn) and the associated charge shift at M. We show that the phenomenon in which the charge at the central atom, qM, is itself independent of n (e.g., at carbon in CH4-nBrn) is best described as an “M-neutral substitution”—not saturation. Implications of the observed X saturation and M-neutral substitutions for larger organic and inorganic halogenated molecules and polymeric materials are identified

Aromatic Ouroboroi: Heterocycles Involving a Sigma-Donor-Acceptor Bond and 4n+2 Pi-Electrons

The aromaticity and dynamics of a set of recently proposed neutral 5- and 6-membered heterocycles that are closed by dative (donor–acceptor) or multi-center s bonds, and have resonance forms with a Hu¨ckel number of p-electrons, are examined. The donors and acceptors in the rings include N, O, and F, and B, Be, and Mg, respectively. The planar geometry of the rings, coupled with evidence from different measures of aromaticity, namely the NICSzz, and NICSpzz components of the conventional nucleus independent chemical shifts (NICS), and ring current strengths (RCS), indicate non-trivial degrees of aromaticity in certain cases, including the cyclic C3B2OH6 and C3BOH5 isomers, both with three bonds to the O site in the ring. The former is lower in energy by at least 17.6 kcal mol1 relative to linear alternatives obtained from molecular dynamics simulations in this work. Some of the other systems examined are best described as non-aromatic. Ring opening, closing, and isomerization are observed in molecular dynamics simulations for some of the systems studied. In a few cases, the ring indeed persists

Plane and simple: planar tetracoordinate carbon centers in small moleculesw

A class of neutral 18-electron molecules with planar tetracoordinate carbon (ptC) centers is introduced. We show computationally that when n = 3 the neutral singlet molecule C(BeH)n(BH2)4-n and other isoelectronic (18-valence electron) molecules of main group elements collapse from locally tetrahedral arrangements at the C-center to (near) planar tetracoordinate structures. For C(BeH)3BH2 and C(CH3)(BH2)Li2, for example, the tetrahedral type conformation is not even a minimum on the potential energy surface at the B3PW91, MP2(full), or CCSD levels of theory. The Mg analogue C(MgH)3BH2 of the Be system also features a completely flat global minimum (with even higher energy planar minima in both cases as well). Other neutral compounds that may prefer planar geometries are apparent, and new openings for experimental investigations and theoretical analyses of planar tetracoordinate main group systems are identified. The planar conformation persists at one center in the C(BeH)3BH2 dimer, and may be identifiable in higher order clusters of ptC molecules as well

Theoretical design of stable small aluminium-magnesium binary clusters

We explore in detail the potential energy surfaces of the AlxMgy (x, y = 1–4) systems as case studies to test the utility and limitations of simple rules based on electron counts and the phenomenological shell model (PSM) for bimetallic clusters. We find that it is feasible to design stable structures that are members of this set of small Al–Mg binary clusters, using simple electron count rules, including the classical 4n + 2 Hückel model, and the most recently proposed PSM. The thermodynamic stability of the title compounds has been evaluated using several different descriptors, including the fragmentation energies and the electronic structure of the systems. Three stable systems emerge from the analysis: the Al4Mg, Al2Mg2 and Al4Mg4 clusters. The relative stability of Al4Mg is explained by the stability of the Al42- subunit to which the Mg atom donates its electrons. Here the Mg2+ sits above the aromatic 10 π-electron Al42- planar ring. The Al2Mg2 and Al4Mg4 clusters present more complicated 3D structures, and their stabilities are rationalized as a consequence of their closed shell nature in the PSM, with 10 and 20 itinerant electrons, respectively

Stabilizing carbon-lithium stars

We have explored in silico the potential energy surfaces of the C5Linn-6 (n = 5, 6, and 7) clusters using the Gradient Embedded Genetic Algorithm (GEGA) and other computational strategies. The most stable forms of C5Li5-- and C5Li6 are two carbon chains linked by two lithium atoms in a persistent seven membered ring capped by two Li atoms. The other Li atoms are arrayed on the edge of the seven membered ring. In contrast, the global minimum structure for C5Li7+ is a bicapped star of D5h symmetry. The molecular orbital analysis and computed magnetic field data suggest that electron delocalization, as well as the saturation of the apical positions of the five-membered carbon ring with lithium atoms in C5Li7+ plays a key role in the stabilization of the carbon-lithium star. In fact, the planar star sub-structure for the carbon ring are unstable without the apical caps. This is also what has been found for the Si analogues. The split of the Bindz in its σ- and π-contribution indicates that C5Li7+ is a p-aromatic and σ-nonaromatic system

Shorter Still: Compressing C-C Single Bonds

How short can a C-C single bond get? The bonding in a set of molecules that are related structurally to previously synthesized or theoretically examined systems with short C-C bonds is investigated. According to calculations, a single C-C bond could be compressed to 1.313 A! To the best of our knowledge, this is the shortest single C-C bond reported to date. This shortening is a consequence of a change in the C-C-C bond angle, θ, to minimize strain in the cages and an effort to offset the tension in the surrounding bridges