Simultaneous Aromatic–Beryllium Bonds and Aromatic–Anion Interactions: Naphthalene and Pyrene as Models of Fullerenes,Carbon Single‐Walled Nanotubes,and Graphene

The possibility of forming stable BeR₂:ArH:Y^? (R=H, F, Cl; ArH=naphthalene, pyrene; Y=Cl, Br) ternary complexes in which the beryllium compounds and anions are located on the opposite sides of an extended aromatic system is explored by means of MP2/aug‐cc‐pVDZ ab initio calculations. Comparison of the electron‐density distribution of these ternary complexes with the corresponding BeR₂:ArH and ArH:Y^? binary complexes reveals the existence of significant cooperativity between the two noncovalent interactions in the triads. The energetic effects of this cooperativity are quantified by evaluation of the three‐body interaction energy Δ³E in the framework of the many‐body interaction‐energy (MBIE) approach. Although an essential component of the interaction energies is electrostatic and is well reflected in the changes in the molecular electrostatic potential of the aromatic system on complexation, strong polarization effects, in particular for the BeR₂:ArH interactions, also play a significant role. The charge transfers associated with these polarization effects are responsible for significant distortion of both the BeR₂ and the aromatic moieties. The former are systematically bent in all the complexes, and the latter are curved to a degree that depends on the nature of the R substituents of the BeR₂ subunit.     

A New Approach to the Design of Neutral 10‐C‐5 Trigonal‐Bipyramidal Carbon Compounds: A “π‐Electron Cap” Effect

DFT (B3LYP, M06‐2X) and MP2 methods are applied to the design of a wide series of the potentially 10‐C‐5 neutral compounds based on 6‐azabicyclotetradecanes: XC¹(YCH₂CH₂CH₂)₃N 1 – 3 , XC¹(YC₆H₄CH₂)₃N 4 – 6 , XC¹[Y(tBuC₆H₃)CH₂]₃N 7 – 9 and carbatranophanes 10 – 25 (X=Me, F, Cl; Y=O, NH, CH₂, SiH₂; Z=O, CH₂, (CH₂)₂, (CH₂)₃). Carbatranophanes 10 – 25 are characterized by a sterical compression of their axial 3c–4e XC¹←N fragment with respect to that in the parent molecules 4 – 6 . A magnitude of the revealed effect depends on a valence surrounding of the central carbon atom C¹, the size and the nature of the side chains (Z) that link the “π‐electron cap” with a tetradecane backbone. This circumstance allowed us to obtain 10‐C‐5 structures with the configuration of the bonds around the C¹ atom, which corresponds to practically an ideal trigonal bipyramid. In these compounds, the values of the covalence ratio χ of approximately 0.6 for the coordination C¹←N contacts with a covalent contribution (atoms in molecules (AIM) and natural bond orbital (NBO)) are record in magnitude. These values lie close to a low limit of the interval of the χ_Si←D change (0.6–0.9) being characteristic of the dative and ionic‐covalent (by nature) Si←D bond (D=N, O) in the known 10‐Si‐5 silicon compounds.     

Reaction of Paramagnetic Synthon,Lithiated 4,4,5,5‐Tetramethyl‐4,5‐dihydro‐1H‐imidazol‐1‐oxyl 3‐oxide,with Cyclic Aldonitrones of the Imidazole Series

It was shown that dipole‐stabilized paramagnetic carbanion lithiated 4,4,5,5‐tetramethyl‐4,5‐dihydro‐1H‐imidazol‐1‐oxyl 3‐oxide can be attached in a nucleophilic manner to either isolated or conjugated aldonitrones of the 2,5‐dihydroimidazole 3‐oxide and 2H‐imidazole 1‐oxide series to afford adducts the subsequent oxidation of which leads to polyfunctional mono‐ and diradicals. According to XRD, at least two polymorphic modifications can be formed during crystallization of the resulting paramagnetic compounds, and for each of them, geometric parameters of the molecules are similar. An EPR spectrum of the diradical in frozen toluene has a complicated lineshape, which can be fairly well reproduced by using X‐ray diffraction structural analysis and the following set of parameters: D=14.9 mT, E=1.7 mT; tensor a(¹⁴N)=[0.260 0.260 1.625] mT, two equivalent tensors for the nitronyl nitroxide moiety a(¹⁴N)=[0.198 0.198 0.700] mT, and g≈2.007. According to our DFT and ab initio calculations, the intramolecular exchange in the diradical is very weak and most likely ferromagnetic.     

On the Reaction of Glycerol Dehydratase with But‐3‐ene‐1,2‐diol

Intriguing inactivation : Calculations suggest that the ability of relatively high‐energy radical intermediates to inactivate glycerol dehydratase (GDH) may reflect a general and hitherto unidentified inactivation mechanism in the reaction of coenzyme B₁₂‐dependent enzymes and 3‐unsaturated 1,2‐diols (see scheme; AdoCbl: adenosylcobalamin or coenzyme B₁₂).

     

Subsystem‐Based Theoretical Spectroscopy of Biomolecules and Biomolecular Assemblies

The absorption properties of chromophores in biomolecular systems are subject to several fine‐tuning mechanisms. Specific interactions with the surrounding protein environment often lead to significant changes in the excitation energies, but bulk dielectric effects can also play an important role. Moreover, strong excitonic interactions can occur in systems with several chromophores at close distances. For interpretation purposes, it is often desirable to distinguish different types of environmental effects, such as geometrical, electrostatic, polarization, and response (or differential polarization) effects. Methods that can be applied for theoretical analyses of such effects are reviewed herein, ranging from continuum and point‐charge models to explicit quantum chemical subsystem methods for environmental effects. Connections to physical model theories are also outlined. Prototypical applications to optical spectra and excited states of fluorescent proteins, biomolecular photoreceptors, and photosynthetic protein complexes are discussed.     

Double Pancake Versus Long Chalcogen–Chalcogen Bonds in Six‐Membered C,N,S‐Heterocycles

The double “pancake” bonding in the dimers of the six‐membered heterocycles 1,3‐dithia‐2,4,6‐triazine ( 4 ) and 1,3‐dithia‐2,4‐diazine ( 16 ) were investigated by means of high‐level quantum chemical calculations (B3LYP and CCSD(T)). It was found that the S–S dimers, 20 a and 27 , are not the most stable isomers, but the dimers showing short S?N ( 21 a ) and S?C ( 25 , 28 ) bonds. An investigation of the 5‐phenyl‐1,3‐dithia‐2,4,6‐triazine ( 4 b ) yields that the syn dimer with two S?S bonds (2.57 Å) is the most stable one. In this dimer, the phenyl groups are placed on top of each other. The additional dispersion energy of the phenyl rings causes a stabilization of the syn‐S–S (C_2v‐like) isomer. As a result, two weak albeit relevant single S?S bonds (2.57 Å) are predicted. These findings contradict the recently published concept of double “pancake” bonding in the dimer 4 b ₂.     

Intermolecular Interactions in Weak Donor–Acceptor Complexes from Symmetry‐Adapted Perturbation and Coupled‐Cluster Theory: Tetracyanoethylene–Benzene and Tetracyanoethylene–p‐Xylene

The interactions in the complexes of tetracyanothylene (TCNE) with benzene and p‐xylene, often classified as weak electron donor–acceptor (EDA) complexes, are investigated by a range of quantum chemical methods including intermolecular perturbation theory at the DFT‐SAPT (symmetry‐adapted perturbation theory combined with density functional theory) level and explicitly correlated coupled‐cluster theory at the CCSD(T)‐F12 level. The DFT‐SAPT interaction energies for TCNE–benzene and TCNE–p‐xylene are estimated to be ?35.7 and ?44.9 kJ mol^?1, respectively, at the complete basis set limit. The best estimates for the CCSD(T) interaction energy are ?37.5 and ?46.0 kJ mol^?1, respectively. It is shown that the second‐order dispersion term provides the most important attractive contribution to the interaction energy, followed by the first‐order electrostatic term. The sum of second‐ and higher‐order induction and exchange–induction energies is found to provide nearly 40 % of the total interaction energy. After addition of vibrational, rigid‐rotor, and translational contributions, the computed internal energy changes on complex formation approach results from gas‐phase spectrophotometry at elevated temperatures within experimental uncertainties, while the corresponding entropy changes differ substantially.     

Interaction Modes and Absolute Affinities of α‐Amino Acids for Mn2+: A Comprehensive Picture

Manganese is involved as a cofactor in the activation of numerous enzymes as well as the oxygen‐evolving complex of photosystem II. Full understanding of the role played by the Mn²⁺ ion requires detailed knowledge of the interaction modes and energies of manganese with its various environments, a knowledge that is far from complete. To bring detailed insight into the local interactions of Mn in metallopeptides and proteins, theoretical studies employing first‐principles quantum mechanical calculations are carried out on [Mn‐amino acid]²⁺ complexes involving all 20 natural α‐amino acids (AAs). Detailed investigation of [Mn‐serine]²⁺_, [Mn‐cysteine]²⁺, [Mn‐phenylalanine]²⁺, [Mn‐tyrosine]²⁺, and [Mn‐tryptophan]²⁺ indicates that with an electron‐rich side chain, the most stable species involves interaction of Mn²⁺ with carbonyl oxygen, amino nitrogen, and an electron‐rich section of the side chain of the AA in its canonical form. This is in sharp contrast with aliphatic side chains for which a salt bridge is formed. For aromatic AAs, complexation to manganese leads to partial oxidation as well as aromaticity reduction. Despite multisite binding, AAs do not generate strong enough ligand fields to switch the metal to a low‐ or even intermediate‐spin ground state. The affinities of Mn²⁺ for all AAs are reported at the B3LYP and CCSD(T) levels of theory, thereby providing the first complete series of affinities for a divalent metal ion. The trends are compared with those of other cations for which affinities of all AAs have been previously obtained.     

Static and Frequency‐Dependent Dipole–Dipole Polarizabilities of All Closed‐Shell Atoms up to Radium: A Four‐Component Relativistic DFT Study

We test the performance of four‐component relativistic density functional theory by calculating the static and frequency‐dependent electric dipole–dipole polarizabilities of all (ground‐state) closed‐shell atoms up to Ra. We consider 12 nonrelativistic functionals, including three asymptotically shape‐corrected functionals, by using two smooth interpolation schemes introduced by the Baerends group: the gradient‐regulated asymptotic connection (GRAC) procedure and the statistical averaging of (model) orbital potentials (SAOP). Basis sets of doubly augmented triple‐zeta quality are used. The results are compared to experimental data or to accurate ab initio results. The reference static electric dipole polarizability of palladium has been obtained by finite‐field calculations using the coupled‐cluster singles, doubles, and perturbative triples method within this work. The best overall performance is obtained using hybrid functionals and their GRAC shape‐corrected versions. The performance of SAOP is among the best for nonhybrid functionals for Group 18 atoms but its precision degrades when considering the full set of atoms. In general, we find that conclusions based on results obtained for the rare‐gas atoms are not necessarily representative of the complete set of atoms. GRAC cannot be used with effective core potentials since the asymptotic correction is switched on in the core region.     

The Stability of α‐Hydroperoxyalkyl Radicals

High‐level ab initio and Born–Oppenheimer molecular dynamic calculations have been carried out on a series of hydroperoxyalkyl (α‐QOOH) radicals with the aim of investigating the stability and unimolecular decomposition mechanism into QO+OH of these species. Dissociation was shown to take place through rotation of the C?O(OH) bond rather than through elongation of the CO?OH bond. Through the C?O(OH) rotation, the unpaired electron of the radical overlaps with the electron density on the O?OH bond, and from this overlap the C=O π bond forms and the O?OH bond breaks spontaneously. The CH₂OOH, CH(CH₃)OOH, CH(OH)OOH, and α‐hydroperoxycycloheptadienyl radical were found to decompose spontaneously, but the CH(CHO)OOH has a decomposition energy barrier of 5.95 kcal mol^?1 owing to its steric and electronic features. The systems studied in this work provide the first insights into how structural and electronic effects govern the stabilizing influence on elusive α‐QOOH radicals.     

Structural Characteristics of Graphane‐Type C and BN Nanostructures by Periodic Local MP2 Approach

The structural characteristics of fully‐hydrogenated carbon and boron nitride mono‐ and multilayer slabs, together with nanotubes derived from the slabs, are investigated mainly by means of periodic local second‐order Møller–Plesset perturbation (LMP2) calculations and the results are compared with Hartree–Fock (HF), density functional theory (DFT), and dispersion function‐augmented DFT (DFT‐D) obtained ones. The investigated systems are structurally analogous to (111) and (110) slabs of diamond, where the hydrogenated (111) slab of diamond corresponds to the experimentally known graphane. Multilayering of monolayers and nanotubes is energetically favorable at the LMP2 level for both C and BN, while HF and DFT are not able to reproduce this behavior for CH systems. The work highlights the importance of utilizing methods capable of properly describing weak interactions in the investigation of dispersively‐bound systems such as the multilayered graphanes and the corresponding nanotubes.     

Realistic Energy Surfaces for Real‐World Systems: An IMOMO CCSD(T):DFT Scheme for Rhodium‐Catalyzed Hydroformylation with the 6‐DPPon Ligand

The hydroformylation of terminal alkenes is one of the most important homogeneously catalyzed processes in industry, and the atomistic understanding of this reaction has attracted enormous interest in the past. Herein, the whole catalytic cycle for rhodium‐catalyzed hydroformylation with the 6‐diphenylphosphinopyridine‐(2H)‐1‐one (6‐DPPon) ligand 1 was studied. This catalytic transformation is challenging to describe computationally, since two requirements must be met: 1) changes in the hydrogen‐bond network must be modeled accurately and 2) bond‐formation/bond‐breaking processes in the coordination sphere of the rhodium center must be calculated accurately. Depending on the functionals used (BP86, B3LYP), the results were found to differ strongly. Therefore, the complete cycle was calculated by using highly accurate CCSD(T) computations for a PH₃ model ligand. By applying an integrated molecular orbital plus molecular orbital (IMOMO) method consisting of CCSD(T) as high level and DFT as low‐level method, excellent agreement between the two functionals was achieved. To further test the reliability of the calculations, the energetic‐span model was used to compare experimentally derived and computed activation barriers. The accuracy of the new IMOMO method apparently makes it possible to predict the catalytic potential of real‐world systems.     

Model studies of the optical rotation,and theoretical determination of its sign for β‐pinene and trans‐pinane

We have carried out extensive studies on the basis set dependence of the calculated specific optical rotation (OR) in molecules at the level of the time–dependent Hartree–Fock and density functional approximations. To reach the limits of the basis set saturation, we have devised an artificial model, the asymmetrically deformed (chiral) methane (CM) molecule. This small system permits to use basis sets which are prohibitively large for real chiral molecules and yet shows all the important features of the basis set dependence of the OR values. The convergence of the OR has been studied with n‐aug‐cc‐pVXZ basis sets of Dunning up to the 6–ζ. In a parallel series of calculations, we have used the recently developed large polarized (LPolX) basis sets. The relatively small LPolX sets have been shown to be competitive to very large n‐aug‐cc‐pVXZ basis sets. The conclusions reached in calculations of OR in CM concerning the usefulness of LPolX basis sets have been further tested on (S)‐methyloxirane and (S)‐fluoro‐oxirane. The smallest set of the LPolX family (LPol–ds) has been found to yield OR values of similar quality as those obtained with much larger Dunning's aug‐cc‐pVQZ basis set. These results have encouraged us to carry out the OR calculations with LPol–ds basis sets for systems as large as β‐pinene and trans‐pinane. In both cases, our calculations have lead to the correct sign of the OR value in these molecules. This makes the relatively small LPol–ds basis sets likely to be useful in OR calculations for large molecules. © 2009 Wiley Periodicals, Inc. J Comput Chem, 2010     

Anion Recognition with Hydrogen‐Bonding Cyclodiphosphazanes

Modular cyclodiphosph(V)azanes are synthesised and their affinity for chloride and actetate anions were compared to those of a bisaryl urea derivative ( 1 ). The diamidocyclodiphosph(V)azanes cis‐[{ArNHP(O)(μ‐tBu)}₂] [Ar=Ph ( 2 ) and Ar=m‐(CF₃)₂Ph ( 3 )] were synthesised by reaction of [{ClP(μ‐NtBu)}₂] ( 4 ) with the respective anilines and subsequent oxidation with H₂O₂. Phosphazanes 2 and 3 were obtained as the cis isomers and were characterised by multinuclear NMR spectroscopy, FTIR spectroscopy, HRMS and single‐crystal X‐ray diffraction. The cyclodiphosphazanes 2 and 3 readily co‐crystallise with donor solvents such as MeOH, EtOH and DMSO through bidentate hydrogen bonding, as shown in the X‐ray analyses. Cyclodiphosphazane 3 showed a remarkably high affinity (log[K]=5.42) for chloride compared with the bisaryl urea derivative 1 (log[K]=4.25). The affinities for acetate (AcO^?) are in the same range ( 3 : log[K]=6.72, 1 : log[K]=6.91). Cyclodiphosphazane 2 , which does not contain CF₃ groups, exhibits weaker binding to chloride (log[K]=3.95) and acetate (log[K]=4.49). DFT computations and X‐ray analyses indicate that a squaramide‐like hydrogen‐bond directionality and C_α?H interactions account for the efficiency of 3 as an anion receptor. The C_α?H groups stabilise the Z,Z‐ 3 conformation, which is necessary for bidentate hydrogen bonding, as well as coordinating with the anion.