Noncovalent Interactions of Silatranes and Germatranes with the Surface of Stoichiometric and Hydroxylated 2D Silica.

Analysis of noncovalent interactions formed by the surface of a 2D silica bilayer and atrane molecules, as precursors of functional layers immobilized on a surface of silicatene, was performed. For this purpose, the systems of substituted silatranes and germatranes adsorbed on silicatene surfaces with different amounts of hydroxyl groups (0, 2, 4, and 30 per cell) were simulated by using quantum chemical modeling with periodic boundary conditions and full-electron basis sets. The observation results for interaction energy changes in the systems "atrane molecule-silicatene surface" with increasing silanol number of the silicatene surface can be used to predict the optimal degree of silicatene hydroxylation in order to control the effective progress of atrane deprotection on activated 2D silica materials. In addition to the typical hydrogen bonds, the ability of atranes to form noncovalent O···O and O···Hal interactions was discovered. In these bonds, the oxygen or halogen atoms of atranes act as electron-donor centers in relation to the silicatene oxygen atoms. The observed weakening of the Ge-O covalent bonds in germatranes, on which further deprotection reaction depends, is more manifested than for the Si-O bonds in adsorbed silatranes.

Electron delocalization in defect-containing graphene and its influence on tetrel bond formation.

Using the advanced analyses of electron density and fermionic potential, we show how electron delocalization influences the ability of defect-containing graphene to form tetrel bonds. The Cg atoms of a vacancy defect can produce one nonpolar interaction, alongside a peculiar polar Cg⋯Cg bond. The latter stems from the presence of a localized electron pair on a vacancy defect Cg atom and the local depletion of electron localization on another Cg atom. This interaction is an example of intralayer tetrel bond. In the presence of an absorbed molecule of bisphenol A diglycidyl ether (DGEBA), graphene is able to form incipient tetrel Cg⋯O bonds with an ether group oxygen. In contrast to an epoxy group oxygen, the disposition of the ether oxygen often causes the orientation of electron-rich π-domains of graphene carbon on the weakly expressed electrophilic region of the oxygen. In the case of graphene with a point Si defect, the Si atom can form quite strong Si⋯C interactions with the DGEBA aryl carbons. In contrast to other noncovalent bonds, this interaction significantly alters the electron (de)localization on the Si atom and in the aryl ring. The reliability of the obtained results is enhanced by the use of multiple 2D periodic models with defects located at different positions along the DGEBA skeleton.

Molecular Dynamics Study of Cured ED-20 Epoxy Resin for Predicting the Glass Transition Temperature and Relationship with Structure Features.

A structural model of ED-20 epoxy resin cured by triethylenetetramine was constructed using a molecular dynamics (MD) simulation method. In order to model the epoxy resin hardening, we modified a typical stepwise protocol introducing cross-links between the amino groups and epoxide carbon atoms that allowed us to reproduce experimentally observed glass transition temperature at the value of 360-364 K. The set of MD trajectories of the final molecular-mechanical model can be useful for analysis of alterations in structural and physical properties of the epoxy resin depending on scaleable parameters. The relationships among qualitative alterations in macromolecular structure, the glass transition temperature, and the number of imposed cross-links were elucidated. Analysis of intermolecular interactions between the largest macromolecules and other molecules of a system made it possible to observe the features of the transition from glassy to the viscoelastic state of the studied polymer during the temperature rise.

Can We Merge the Weak and Strong Tetrel Bonds? Electronic Features of Tetrahedral Molecules Interacted with Halide Anions.

Using the orbital-free quantum crystallography approach, we have disclosed the quantitative trends in electronic features for bonds of different strengths formed by tetrel (Tt) atoms in stable molecular complexes consisting of electrically neutral tetrahedral molecules and halide anions. We have revealed the role of the electrostatic and exchange-correlation components of the total one-electron static potential that are determined by the equilibrium atomic structure and by kinetic Pauli potential, which reflects the spin-dependent electron motion features of the weak and strong bonds. The gap between the extreme positions in the electrostatic and total static potentials along the line linking the Tt atom and halide anion is wide for weak bonds and narrow for strong ones. It is in very good agreement with the number of minima in the Pauli potential between the bounded atoms. This gap exponentially correlates with the exchange-correlation potential in various series with a fixed nucleophilic fragment. A criterion for categorizing the noncovalent tetrel bonds (TtB) based on the potential features is suggested.

The Distance between Minima of Electron Density and Electrostatic Potential as a Measure of Halogen Bond Strength.

In this study, we present results of a detailed topological analysis of electron density (ED) of 145 halogen-bonded complexes formed by various fluorine-, chlorine-, bromine-, and iodine-containing compounds with trimethylphosphine oxide, Me3PO. To characterize the halogen bond (XB) strength, we used the complexation enthalpy, the interatomic distance between oxygen and halogen, as well as the typical set of electron density properties at the bond critical points calculated at B3LYP/jorge-ATZP level of theory. We show for the first time that it is possible to predict the XB strength based on the distance between the minima of ED and molecular electrostatic potential (ESP) along the XB path. The gap between ED and ESP minima exponentially depends on local electronic kinetic energy density at the bond critical point and tends to be a common limiting value for the strongest halogen bond.

Acylative kinetic resolution of racemic methyl-substituted cyclic alkylamines with 2,5-dioxopyrrolidin-1-yl (R)-2-phenoxypropanoate.

The diastereoselective acylation of a number of racemic methyl-substituted cyclic alkylamines with active esters of 2-phenoxypropanoic acid was studied in detail. The ester of (R)-2-phenoxypropanoic acid and N-hydroxysuccinimide was found to be the most selective agent. The highest stereoselectivity was observed in the kinetic resolution of racemic 2-methylpiperidine in toluene at -40 °C (selectivity factor s = 73) with the predominant formation of (R,R)-amide (93.7% de). To explain the observed stereoselectivity, DFT modelling of the transition states in the reactions of the title acylating agent with 2-methylpiperidine and 2-methylpyrrolidine was performed. The calculated values were in good agreement with experimental data. It has been demonstrated that the acylation proceeds via a concerted mechanism, in which the addition of an amine occurs simultaneously with the elimination of the hydroxysuccinimide fragment. The high stereoselectivity of the (R,R)-amide formation is largely ensured by the lower steric hindrances in the transition states as compared to the formation of (R,S)-amide.

Pseudocyclic Form of 4-Hydroxypyrrolidine-2-carboxanilide Podands with Trioxyethylene Chain: Modeling, Conformational Search, and NMR Analysis.

The 4-hydroxypyrrolidine-2-carboxanilide podand salt demonstrates catalytic activity in asymmetric Biginelli reaction. The systematic search for prevalent conformational state of the cation was carried out by computer simulations in combination with one- and two-dimensional NMR experiments. For that purpose, we proposed a novel algorithm for the generation and selection of conformers based on molecular dynamics and clustering in the space of principal components. The search had found an important trend of the podand to form a pseudocyclic structure with a horseshoe-shaped conformation of the oligooxyethylene fragment. This conformation is stabilized by different types of intramolecular hydrogen bonds between the acidic and basic centers of the two 4-hydroxypyrrolidine-2-carboxanilide residuals (branches). The proposed approach had made it possible to identify the major structural factors, providing a correlation between the calculated and experimental chemical shifts of hydrogen atoms in the 1H NMR spectra of the protonated podand.

Variations of quantum electronic pressure under the external compression in crystals with halogen bonds assembled in Cl3-, Br3-, I3-synthons.

The inner-crystal quantum electronic pressure was estimated for unstrained C6Cl6, C6Br6, and C6I6 crystals and for those under external compression simulated from 1 to 20 GPa. The changes in its distribution were analyzed for the main structural elements in considered crystals: for triangles of the typical halogen bonds assembled in Hal3-synthons, where Hal = Cl, Br, I; for Hal...Hal stacking interactions, as well as for covalent bonds. Under simulated external compression, the quantum electronic pressure in the intermolecular space reduces as the electron density increases, indicating spatial areas of relatively less crystal resistance to external compression. The most compliant C6Cl6 crystal shows the largest changes of quantum electronic pressure in the centre of Cl3-synthon while the deformation of rigid I3-synthon under external compression depends only on the features of I...I halogen bonds.

Identification of the Tetrel Bonds between Halide Anions and Carbon Atom of Methyl Groups Using Electronic Criterion.

The consideration of the disposition of minima of electron density and electrostatic potential along the line between non-covalently bound atoms in systems with Hal-···CH3⁻Y (Hal- = Cl, Br; Y = N, O) fragments allowed to prove that the carbon atom in methyl group serves as an electrophilic site provider. These interactions between halide anion and carbon in methyl group can be categorized as the typical tetrel bonds. Statistics of geometrical parameters for such tetrel bonds in CSD is analyzed. It is established that the binding energy in molecular complexes with tetrel bonds correlate with the potential acting on an electron in molecule (PAEM). The PAEM barriers for tetrel bonds show a similar behavior for both sets of complexes with Br- and Cl- electron donors.

Electronic criterion for categorizing the chalcogen and halogen bonds: sulfur-iodine interactions in crystals.

Diversity of mutual orientations of Y-S and I-X and covalent bonds in molecular crystals complicate categorizing noncovalent chalcogen and halogen bonds. Here, the different types of S...I interactions with short interatomic distances are analysed. The selection of S...I interactions for the categorization of the chalcogen and halogen bonds has been made using angles that determine the mutual orientation of electron lone pairs and σ-holes interacted S and I atoms. In complicated cases of noncovalent interactions with `hole-to-hole' of S and I orientations, distinguishing the chalcogen and halogen bonds is only possible if the atom is uniquely determined, which also provides the electrophilic site. The electronic criterion for chalcogen/halogen bonds categorizing that is based on analysis of dispositions of electron density and electrostatic potential minima along the interatomic lines has been suggested and its effectiveness has been demonstrated.

Why are reactions of 2- and 8-thioquinoline derivatives with iodine different?

The crystal structures of 1,2-dihydro-1,1'-bi[thiazolo[3,2-a]quinoline]-10a,10a'-diium diiodide hemihydrate, C22H16N2S22+·2I-·0.5H2O, and 1,2-dihydro-1,1'-bi[thiazolo[3,2-a]quinoline]-10a,10a'-diium iodide triiodide, C22H16N2S22+·I-·I3-, obtained during the reaction of 1,4-bis(quinolin-2-ylsulfanyl)but-2-yne (2TQB) with iodine, have been determined at 120 K. The crystalline products contain the dication as a result of the reaction proceeding along the iodocyclization pathway. This is fundamentally different from the previously observed reaction of 1,4-bis(quinolin-8-ylsulfanyl)but-2-yne (8TQB) with iodine under similar conditions. A comparative analysis of the possible conformational states indicates differences in the relative stabilities and free rotation for the 2- and 8-thioquinoline derivatives which lead to a disparity in the convergence of the potential reaction centres for 2TQB and 8TQB.

A comparative view on the potential acting on an electron in a molecule and the electrostatic potential through the typical halogen bonds.

This work considers the features of the electrostatic potential (ESP), and the potential acting on an electron in a molecule (PAEM) for the series of isolated dihalide molecules and for their molecular complexes. The joint analysis of these functions enriches the vision of atomic predispositions to the halogen bond formation and reveals details for their characterization. The account for the exchange-correlation interaction in PAEM retains the specific anisotropy of the ESP, which is commonly used for the halogen bonding explanation within σ-hole concept. Along the halogen bonds, the curvatures of PAEM and ESP functions are opposite. Being jointly mapped on the closed isosurfaces of the reduced density gradient, placed between bound atoms, they are significantly differed from the side facing the halogen atom and from the side looking at the electron donor atom. © 2017 Wiley Periodicals, Inc.

Testing the tools for revealing and characterizing the iodine-iodine halogen bond in crystals.

To understand what tools are really suitable to identify and classify the iodine-iodine non-covalent interactions in solid organic polyiodides, we have examined the anisotropy of the electron density within the iodine atomic basin along and across the iodine-iodine halogen bond using the Laplacian of electron density, one-electron potential and electron localization function produced by Kohn-Sham calculations with periodic boundary conditions. The Laplacian of electron density exhibits the smallest anisotropy and yields a vague picture of the outermost electronic shells. The one-electron potential does not show such a deficiency and reveals that the valence electron shell for the halogen-bond acceptor iodine is always wider than that for the halogen-bond donor iodine along its σ-hole direction. We have concluded that the one-electron potential is the most suitable for classification of the iodine-iodine bonds and interactions in complicated cases, while the electron localization function allows to distinguish the diiodine molecule bonded with the monoiodide anion from the typical triiodide anion.

An anatomy of intramolecular atomic interactions in halogen-substituted trinitromethanes.

The intramolecular interactions in substituted trinitromethanes, XC(NO2)3 (X = F, Cl, I, H) are studied and clarified by using a combination of the Quantum Theory of Atoms in Molecules (QTAIM), the non-covalent interaction analysis and the Interacting Quantum Atoms (IQA) methods. The stretching vibration modes are formed by the concerted displacements of atoms involved in the covalent bonds showing the significant multiatomic influence in substituted trinitromethanes. In agreement with that, the arrangement of the local reduced density gradient minima indicates that the electron density favors the non-covalent intramolecular interactions X···O and N···O. However, the corresponding QTAIM bond paths are not formed; instead, contacts, which we call uncompleted links in this context, are accompanied by "quasi-bonding channels" corresponding to the λ2() ≤ 0 regions on the sign[λ2(r)]ρ(r) contour maps. The intramolecular IQA energy contributions signal the appreciable electron exchange between the pairs of atoms associated with potential atomic interactions or the bond-path-free non-covalent links. The IQA analysis shows that the electrostatic term destabilizes FC(NO2)3 and distinctly stabilizes IC(NO2)3, whereas it is close to neutral in ClC(NO2)3. The exchange energy between the X atom and the NO2 groups, in contrast, stabilizes all the molecules.

Towards crystal structure prediction of complex organic compounds--a report on the fifth blind test.

Following on from the success of the previous crystal structure prediction blind tests (CSP1999, CSP2001, CSP2004 and CSP2007), a fifth such collaborative project (CSP2010) was organized at the Cambridge Crystallographic Data Centre. A range of methodologies was used by the participating groups in order to evaluate the ability of the current computational methods to predict the crystal structures of the six organic molecules chosen as targets for this blind test. The first four targets, two rigid molecules, one semi-flexible molecule and a 1:1 salt, matched the criteria for the targets from CSP2007, while the last two targets belonged to two new challenging categories - a larger, much more flexible molecule and a hydrate with more than one polymorph. Each group submitted three predictions for each target it attempted. There was at least one successful prediction for each target, and two groups were able to successfully predict the structure of the large flexible molecule as their first place submission. The results show that while not as many groups successfully predicted the structures of the three smallest molecules as in CSP2007, there is now evidence that methodologies such as dispersion-corrected density functional theory (DFT-D) are able to reliably do so. The results also highlight the many challenges posed by more complex systems and show that there are still issues to be overcome.

Determination of covalent bond orders and atomic valence indices using topological features of the experimental electron density.

We present an approach for the determination of covalent bond orders from the experimental electron density and its derivatives at the bond critical points. An application of this method to a series of organic compounds has shown that it provides a bonding quantification that is in reasonable agreement with that obtained by orbital theory. The 'experimental' atomic valence indices are also defined and their significance for the characterization of chemical problems is discussed.

Atoms-in-molecules study of intra- and intermolecular bonding in the pentaerythritol tetranitrate crystal.

Chemical bonding in the pentaerythritol tetranitrate crystal based on the experimental electron density obtained from X-ray diffraction data at 100 K and theoretical calculations at the experimental molecular geometry have been analyzed in terms of the Quantum Theory of Atoms in Molecules. Features of the intra- and intermolecular bond critical points and the oxygen atom lone-pair locations are discussed. Numerous intermolecular bonding interactions, including O...H and O...O, have been found and characterized. Atomic charges and atomic energies were integrated and compared with those for similar compounds. The N-O topological bond orders have been calculated for the first time, and the PETN atomic valences have been estimated.