How important is the dispersion interaction for cyclobis(paraquat‐p‐phenylene)‐based molecular “shuttles”? A theoretical study

Inclusion complexes of cyclobis(paraquat‐p‐phenylene) and various aromatic molecules in their neutral and oxidized form were studied at the LMP2/6‐311+G**//BHandHLYP/6‐31G* level of theory, which represents the highest level theoretical study to date for these complexes. The results show that it is dispersion interaction that contributes most to the binding energy. One electron oxidation of a guest molecule leads to complete dissociation of inclusion complex generating strong repulsion potential between guest and host molecules. Electrostatic interactions also can play an important role, provided the guest molecule has a dipole moment; however, dispersion interactions always dominate in binding energy. © 2004 Wiley Periodicals, Inc. Int J Quantum Chem, 2005     

A theoretical study on three conformational structures of 2,6‐distyrylpyridine

Ab initio methods at the levels HF/cc‐pVDZ, HF/6‐31G(d,p), MP2/cc‐pVDZ, and MP2/6‐31G(d,p), as well as methods based on density functional theory (DFT) employing the hybrid functional B3LYP with the basis sets cc‐pVDZ and 6‐31G(d,p), have been applied to study the conformers of 2,6‐distyrylpyridine. Bond distances, bond angles, and dihedral angles have been calculated at the B3LYP level. The calculated values were in good agreement with those measured by X‐ray diffraction analysis of 2,6‐distyrylpyridine. The values calculated using the Hartree‐Fock method and second‐order perturbation theory (MP2) were inconsistent. The optimized lowest‐energy geometries were calculated from the reported X‐ray structural data by the B3LYP/cc‐pVDZ method. Three conformations, A, B, and C, were proposed for 2,6‐distyrylpyridine. Calculations at the three levels of theory indicated that conformation A was the most stable structure, with conformations C and B being higher in energy by 1.10 and 2.57 kcal/mol, respectively, using the same method and basis function. The same trend in the relative energies of the three possible conformations was observed at the two levels of theory and with the different basis sets employed. The reported X‐ray data were utilized to optimize total molecular energy of conformation A at the different calculation levels. The bond lengths, bond angles, and dihedral angles were then obtained from the optimized geometries by ab initio methods and by applying DFT using the two basis functions cc‐pVDZ and 6‐31G(d,p). The values were analyzed and compared. The calculated total energies, the relative energies of the molecular orbitals, the gap between them, and the dipole moment for each conformational structure proposed for 2,6‐distyrylpyridine are also reported. © 2009 Wiley Periodicals, Inc. Int J Quantum Chem, 2010     

Pressure‐Tuneable Visible‐Range Band Gap in the Ionic Spinel Tin Nitride

The application of pressure allows systematic tuning of the charge density of a material cleanly, that is, without changes to the chemical composition via dopants, and exploratory high‐pressure experiments can inform the design of bulk syntheses of materials that benefit from their properties under compression. The electronic and structural response of semiconducting tin nitride Sn₃N₄ under compression is now reported. A continuous opening of the optical band gap was observed from 1.3 eV to 3.0 eV over a range of 100 GPa, a 540 nm blue‐shift spanning the entire visible spectrum. The pressure‐mediated band gap opening is general to this material across numerous high‐density polymorphs, implicating the predominant ionic bonding in the material as the cause. The rate of decompression to ambient conditions permits access to recoverable metastable states with varying band gaps energies, opening the possibility of pressure‐tuneable electronic properties for future applications.     

Theoretical Studies on the Intermolecular Interactions of Potentially Primordial Base‐Pair Analogues

Recent experimental studies on the Watson–Crick type base pairing of triazine and aminopyrimidine derivatives suggest that acid/base properties of the constituent bases might be related to the duplex stabilities measured in solution. Herein we use high‐level quantum chemical calculations and molecular dynamics simulations to evaluate the base pairing and stacking interactions of seven selected base pairs, which are common in that they are stabilized by two N? H???O hydrogen bonds separated by one N? H???N hydrogen bond. We show that neither the base pairing nor the base stacking interaction energies correlate with the reported pK_a data of the bases and the melting points of the duplexes. This suggests that the experimentally observed correlation between the melting point data of the duplexes and the pK_a values of the constituent bases is not rooted in the intrinsic base pairing and stacking properties. The physical chemistry origin of the observed experimental correlation thus remains unexplained and requires further investigations. In addition, since our calculations are carried out with extrapolation to the complete basis set of atomic orbitals and with inclusion of higher electron correlation effects, they provide reference data for stacking and base pairing energies of non‐natural bases.     

Theoretical study of the low‐lying electronic states of the molecular ion KRb+

The potential energy curves of the molecular ion KRb⁺ have been investigated for the 60 lowest molecular states of symmetry ²Σ⁺, ²Π, ²Δ, and Ω = 1/2, 3/2, and 5/2. Using an ab initio method, the calculation has been done in a one active electron approach based on nonempirical pseudopotentials with core valence effects taken into account through parameterized l‐dependent polarization potentials. Using the canonicals functions approach a rovibrational study is done by calculating the eigenvalues E_v, the rotational constants B_v, the centrifugal distortion constants D_v (up to 135 vibrational levels), and the spectroscopic constants ω_e and B_e for the five electronic states (1)²Σ⁺, (3)²Σ⁺, (1)²Π, (1)Ω = 1/2, and (1)Ω = 3/2. No comparison of these values with other results is yet possible because they are given here for the first time. Extensive tables of energy values of E_v, B_v, and D_v are displayed at http://hplasim2.univ‐lyon1.fr/allouche . © 2003 Wiley Periodicals, Inc. Int J Quantum Chem, 2003     

The Low‐Lying Excited States of 2,2′‐Bithiophene: A Theoretical Analysis

The low‐energy regions of the singlet→singlet, singlet→triplet, and triplet→triplet electronic spectra of 2,2′‐bithiophene are studied using multiconfigurational second‐order perturbation theory (CASPT2) and extended atomic natural orbitals (ANO) basis sets. The computed vertical, adiabatic, and emission transition energies are in agreement with the available experimental data. The two lowest singlet excited states, 1¹B_u and 2¹B_u, are computed to be degenerate, a novel feature of the system to be borne in mind during the rationalization of its photophysics. As regards the observed high triplet quantum yield of the molecule, it is concluded that the triplet states 2³A_g and 2³B_u, separated about 0.4 eV from the two lowest singlet excited states, can be populated by intersystem crossing from nonplanar singlet states.     

A C‐linked Glycomimetic in the Gas Phase and in Solution: Synthesis and Conformation of the Disaccharide Manα(1,6)‐C‐ManαOPh

The effect of carbon is subtle but sweet : The flexible C‐linkage in the newly synthesised C‐glycosyl mimetic, Manα(1,6)‐C‐ManαOPh allows OH? π bonding, both in the gas phase and in aqueous solution. This interaction is absent in the O‐linked disaccharide (see figure).

     

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.     

Theoretical Study on Cooperativity Effects between Anion–π and Halogen‐Bonding Interactions

This article analyzes the interplay between lone pair–π (lp–π) or anion–π interactions and halogen‐bonding interactions. Interesting cooperativity effects are observed when lp/anion–π and halogen‐bonding interactions coexist in the same complex, and they are found even in systems in which the distance between the anion and halogen‐bond donor molecule is longer than 9 Å. These effects are studied theoretically in terms of energetic and geometric features of the complexes, which are computed by ab initio methods. Bader′s theory of “atoms in molecules” is used to characterize the interactions and to analyze their strengthening or weakening depending upon the variation of charge density at critical points. The physical nature of the interactions and cooperativity effects are studied by means of molecular interaction potential with polarization partition scheme. By taking advantage of all aforementioned computational methods, the present study examines how these interactions mutually influence each other. Additionally, experimental evidence for such interactions is obtained from the Cambridge Structural Database (CSD).     

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.     

Strong CH⋅⋅⋅Halide Hydrogen Bonds from 1,2,3‐Triazoles Quantified Using Pre‐Organized and Shape‐Persistent Triazolophanes

The structures associated with halide (F^?, Cl^?, Br^?) complexation inside CH hydrogen‐bonding macrocyclic receptors, called triazolophanes, are characterized using density functional theory (DFT). The associated binding energies in the gas and solution phases are evaluated. The ruffles in the empty triazolophane become smoothed‐out upon Cl^?‐ and Br^?‐ion binding directly into the middle of the cavity. The largely pre‐organized cavity morphs into an elliptical shape to facilitate shorter hydrogen bonds in the north and south regions and longer ones west and east. The smaller F^? ion sits in, and flattens‐out, only the north (or south) region. The 1,2,3‐triazoles show shorter CH???Cl^? contacts than for the phenylenes. Both Cl^? and Br^? show the same binding geometries but Cl^? has a larger binding energy consistent with its stronger Lewis basicity. Model triads were used to decompose the overall binding energy into those of its components. In the course of this triad analysis, anion polarization was identified and its contribution to the triad???Cl^? binding energy estimated. Consequently, the binding energies for the individual aryl units within the comparatively non‐polarized triazolophanes were estimated. The 1,2,3‐triazoles are twice as strong as the phenylenes thus contributing most of the interaction energy to Cl^?‐ion binding. Therefore, the 1,2,3‐triazoles appear to approach the hydrogen bond strengths of the NH donors of pyrrole units.     

Are water–aromatic complexes always stabilized due to π–H interactions? LMP2 study

Eight complexes of various aromatic molecules with water have been studied theoretically at the local Møller–Plesset 2nd order theory (LMP2)/aug‐cc‐pVTZ(‐f)//LMP2/6‐31+G* level of theory. Two types of complexes can be formed, depending on the electronic structure of aromatic molecules. Donor hydrocarbons form A‐type complexes, while aromatics bearing electron‐withdrawing substituents form B‐type complexes. A‐type complexes are stabilized due to π–H interactions with the OH bond pointing to the aromatic molecule plane, while B‐type complexes have geometry with the oxygen atom pointing to the aromatic molecule plane stabilized by the interaction of highest occupied molecular orbital (HOMO) of water molecule with π* orbitals of the aromatics. It has been found that a (? HOMO–lowest unoccupied molecular orbital (LUMO)/2 value of aromatic molecule, which can be called “molecular electronegativity,” is useful to predict the type of complex formed by aromatic molecule and water. Aromatic hydrocarbons with “molecular electronegativity” of <0.15 tend to form A‐type complexes, while aromatic molecules with “molecular electronegativity” of <0.15 a.u. form B‐type complexes. The binding energy of water–aromatic complexes undergoes a minimum in the area of switching from A‐type to B type complexes, which can be rationalize in terms of frontier orbital interactions. © 2005 Wiley Periodicals, Inc. Int J Quantum Chem, 2005     

Peptide models XXXI. Conformational properties of hydrophobic residues shaping the core of proteins. An ab initio study of N‐formyl‐L‐valinamide and N‐formyl‐L‐phenylalaninamide

Employing introductory (3‐21G RHF) and medium‐size (6‐311++G** B3LYP) ab initio calculations, complete conformational libraries, containing as many as 27 conformers, have been determined for diamide model systems incorporating the amino acids valine (Val) and phenylalanine (Phe). Conformational and energetic properties of these libraries were analyzed. For example, significant correlation was found between relative energies from 6‐311++G** B3LYP and single‐point B3LYP/6‐311++G**//RHF/3‐21G calculations. Comparison of populations of molecular conformations of hydrophobic aromatic and nonaromatic residues, based on their ab initiorelative energies, with their natural abundance indicates that, at least for the hydrophobic core of proteins, the conformations of Val (Ile, Leu) and Phe (Tyr, Trp) are controlled by the local energetic preferences of the respective amino acids. © 2001 John Wiley & Sons, Inc. J Comput Chem 22: 732–751, 2001     

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.     

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.     

Intermolecular interactions in group 14 hydrides: Beyond CH···HC contacts

In line with previous work in which we established the factors that enhance attractive C?H···H?C dihydrogen interactions in alkanes, an extended theoretical analysis of noncovalent intermolecular interactions in group 14 hydrides is presented here. Remarkably, these weak interactions may play a major role in determining the crystal structures adopted by several families of molecules. A combined structural and computational analysis at the MP2 level allowed us to identify and characterize different interactions of the type E?H···H?E and E···H?E (E = Si, Ge, Sn, and Pb), and to find also the most suitable scenario for the establishment of each particular type. The nature of the interactions has been analyzed in terms of natural charges of the atoms involved and a topological analysis of the electron density of several dimers confirms the existence of H···H and H···E attractive contacts. We have observed that the interaction strength increases when descending down the periodic group and that silicon has a marked tendency to establish Si···H?Si interactions. A size‐dependent backbone effect that reinforces H···H dihydrogen interactions in polyhedral systems has also been found.     

Theoretical Study on the Complexes of Benzene with Isoelectronic Nitrogen‐Containing Heterocycles

The π–π interactions between benzene and the aromatic nitrogen heterocycles pyridine, pyrimidine, 1,3,5‐triazine, 1,2,3‐triazine, 1,2,4,5‐tetrazine, and 1,2,3,4,5‐pentazine are systematically investigated. The T‐shaped structures of all complexes studied exhibit a contraction of the C? H bond accompanied by a rather large blue shift (40–52 cm^?1) of its stretching frequency, and they are almost isoenergetic with the corresponding displaced‐parallel structures at reliable levels of theory. With increasing number of nitrogen atoms in the heterocycle, the geometries, frequencies, energies, percentage of s character at C, and the electron density in the C? H σ antibonding orbital of the complexes all increase or decrease systematically. Decomposition analysis of the total binding energy showed that for all the complexes, the dispersion energy is the dominant attractive contribution, and a rather large attraction originating from electrostatic contribution is compensated by its exchange counterpart.     

Quantum‐Transport Characteristics of a p–n Junction on Single‐Layer TiS3

By using density functional theory and non‐equilibrium Green′s function‐based methods, we investigated the electronic and transport properties of a TiS₃ monolayer p–n junction. We constructed a lateral p–n junction on a TiS₃ monolayer using Li and F adatoms. An applied bias voltage caused significant variability in the electronic and transport properties of the TiS₃ p–n junction. In addition, the spin‐dependent current–voltage characteristics of the constructed TiS₃ p–n junction were analyzed. Important device characteristics were found, such as negative differential resistance and rectifying diode behaviors for spin‐polarized currents in the TiS₃ p–n junction. These prominent conduction properties of the TiS₃ p–n junction offer remarkable opportunities for the design of nanoelectronic devices based on a recently synthesized single‐layered material.     

Spectroscopic and theoretical study of vibrational spectra of hydrogen‐bonded tropolone

Theoretical simulation of the bandshape and fine structure of the ν_s stretching band is presented for tropolone‐H and tropolone‐D taking into account an adiabatic coupling between the high‐frequency O–H(D) stretching and the low‐frequency intra‐ and intermolecular OO stretching modes, and linear and quadratic distortions of the potential energies for the low‐frequency vibrations in the excited state of the O–H(D) stretching vibration. In order to determine the low‐frequency vibrations, the experimental spectra of the polycrystalline tropolone in the far‐infrared and the low‐frequency Raman range have been recorded for the first time. The experimental frequencies in the low‐frequency region are compared with the results of the HF/6‐31G** and Becke3LYP/6‐31G** calculations carried out for the tropolone dimer. ©1999 John Wiley & Sons, Inc. Int J Quant Chem 73: 275–282, 1999     

Extending the Domain‐Averaged Exchange‐Correlation Energies Within the Context of the MC‐QTAIM: Tracing Subtle Variations Induced by Isotope Substitution

Recently, it has been demonstrated that the domain‐averaged exchange‐correlation energies, V^xc, are capable of tracing the covalent character of atom–atom interactions unequivocally and thus pave the way for detailed bonding analysis within the context of the quantum theory of atoms in molecules (QTAIM) [M. García‐Revilla, E. Francisco, P. L. Popelier, A. Martín Pendás, ChemPhysChem 2013 , 14, 1211–1218]. Herein, the concept of V^xc is extended within the context of the newly developed multicomponent QTAIM (MC‐QTAIM). The extended version, , is capable of analyzing nonadiabatic wavefunctions and thus is sensitive to the mass of nuclei and can trace “locally” the subtle electronic variations induced by isotope substitution. To demonstrate this capability in practice, ab initio nonadiabatic wavefunctions for three isotopically substituted hydrogen cyanide molecules, in which the hydrogen nucleus was assumed to be a proton, deuterium, or tritium, were derived. The resulting wavefunctions were then used to compute and it emerged that for the hydrogen–carbon bond, the was distinct for each isotopic composition and varied in line with chemical expectations. Indeed, the introduction of paves the way for the investigation of vast numbers of structural and kinetic isotope effects within the context of the MC‐QTAIM.