Nature of the interaction between ammonia derivatives and carbon disulfide. A theoretical investigation

The interactions between NH₃, its methylated and chlorinated derivatives and CS₂ are investigated by ab initio CCSD(T) and density functional BLYP‐D3 methods. The CCSD(T)/aug‐cc‐pVTZ calculated interaction energies of complexes characterized by the S···N chalcogen bonds range between ?1.71 and ?2.78 kcal mol^?1. The S···N bonds are studied by atoms in molecules, natural bond orbital, and noncovalent interaction methods. The lack of correlation between the interaction energies of methylated amines complexes and the electrostatic potential results from the lone pair effect in aliphatic amines. Different structures of CS₂ complexed with ammonia derivatives, stabilized by other than the S···N chalcogen bonds, are also predicted. These structures are characterized by interaction energies ranging between 1.15 and 3.46 kcal mol^?1. The results show that the complexing ability of CS₂ is not very high but this molecule is able to attack the electrophilic or nucleophilic sites of a guest molecule.     

Is the DPT tautomerization of the long A·G Watson–Crick DNA base mispair a source of the adenine and guanine mutagenic tautomers? A QM and QTAIM response to the biologically important question

Herein, we first address the question posed in the title by establishing the tautomerization trajectory via the double proton transfer of the adenine·guanine (A·G) DNA base mispair formed by the canonical tautomers of the A and G bases into the A*·G* DNA base mispair, involving mutagenic tautomers, with the use of the quantum‐mechanical calculations and quantum theory of atoms in molecules (QTAIM). It was detected that the A·G ? A*·G* tautomerization proceeds through the asynchronous concerted mechanism. It was revealed that the A·G base mispair is stabilized by the N6H···O6 (5.68) and N1H···N1 (6.51) hydrogen bonds (H‐bonds) and the N2H···HC2 dihydrogen bond (DH‐bond) (0.68 kcal·mol^?1), whereas the A*·G* base mispair—by the O6H···N6 (10.88), N1H···N1 (7.01) and C2H···N2 H‐bonds (0.42 kcal·mol^?1). The N2H···HC2 DH‐bond smoothly and without bifurcation transforms into the C2H···N2 H‐bond at the IRC = ?10.07 Bohr in the course of the A·G ? A*·G* tautomerization. Using the sweeps of the energies of the intermolecular H‐bonds, it was observed that the N6H···O6 H‐bond is anticooperative to the two others—N1H···N1 and N2H···HC2 in the A·G base mispair, while the latters are significantly cooperative, mutually strengthening each other. In opposite, all three O6H···N6, N1H···N1, and C2H···N2 H‐bonds are cooperative in the A*·G* base mispair. All in all, we established the dynamical instability of the А*·G* base mispair with a short lifetime (4.83·10^?14 s), enabling it not to be deemed feasible source of the A* and G* mutagenic tautomers of the DNA bases. The small lifetime of the А*·G* base mispair is predetermined by the negative value of the Gibbs free energy for the A*·G* → A·G transition. Moreover, all of the six low‐frequency intermolecular vibrations cannot develop during this lifetime that additionally confirms the aforementioned results. Thus, the A*·G* base mispair cannot be considered as a source of the mutagenic tautomers of the DNA bases, as the A·G base mispair dissociates during DNA replication exceptionally into the A and G monomers in the canonical tautomeric form. © 2013 Wiley Periodicals, Inc.     

Rotation about the C?N bond in 2-aza-1,3-butadiene and the N?N bond in 2,3-diaza-1,3 butadiene: A molecular orbital study

The geometry and energy of 2-aza-1,3-butadiene and 2,3-diaza-1,3-butadiene have been calculated using the 6-31G* basis set as a function of the CNCC and CNNC dihedral angles, respectively. With the 2-aza derivative potential minima are located at 0° (trans) and at about 130° for a gauche structure approximately 9.5 kJ mol^?1 less stable than the trans. Potential maxima are at about 75° giving a gauche barrier height of approximately 19 kJ mol^?1 relative to the trans structure, and at 180° (cis) giving a barrier height of approximately 14.5 kJ mol^?1 relative to the 130° gauche structure. With the 2,3-diaza derivative the gauche barrier has disappeared and there are a series of gauche structures in the region 70°–100° of almost equal energy 12.5-15 kJ mol^?1 less stable than the trans. In addition the cis barrier is much greater, nearly 70 kJ mol^?1 relative to the trans structure. Inclusion of electron correlation, accounting for about 50% of the correlation energy, produces no significant changes in the shape of the potential energy curves. There are systematic and progressive changes in almost all the geometrical parameters as the ?CH? groups in butadiene are replaced by ?N? . The outward tilt and compression within the methylene groups show adverse steric interactions to be operative in the cis structures. The values of V_nn indicate that gauche structures of both the 2-aza and the 2,3-diaza derivatives near the cis structure are more compact (as with butadiene), and gauche structures of the 2-aza derivative near the trans structure are less compact (as with butadiene). Originating in the changes in bond lengths and bond angles, rotation-independent nuclear–nuclear interactions again play an important role.     

Synthesis,Structure, and Solution Study of a Mercury(II) Complex with the Ligand [1‐(2‐Methoxyphenyl)‐3‐(4‐chlorophenyl)]triazene

The title ligand, [1‐(2‐methoxyphenyl)‐3‐(4‐chlorophenyl)]triazene, H L ( 1 ), was prepared. In a reaction with Hg(NO₃)₂ it forms the complex [Hg(C₂₆H₂₂Cl₂N₆O₂)], [Hg L ₂] ( 2 ). Both compounds were characterized by means of X‐ray crystallography, CHN analysis, FT‐IR, ¹H NMR, and ¹³C NMR spectroscopy. In the structure of compound 1 , two independent fragments are present in the unit cell. They exhibit trans arrangement about the –N=N– double bond. The dihedral angles between two benzene rings in both fragments are 4.36 and 18.79 Å, respectively. Non‐classic C–H ··· N hydrogen bonding and C–H ··· π interactions form a layer structure along the crystallographic ab plane [110]. In compound 2 , the Hg^II atom is hexacoordinated by two tridentate [1‐(2‐methoxyphenyl)‐3‐(4‐chlorophenyl)]triazenide ligands through a N₂O₂ set. In addition, in the structure of 2 , monomeric complexes are connected to each other by C–H ··· π stacking interactions, resulting in a 2D architecture. These C–H ··· π edge‐to‐face interactions are present with H ··· π distances of 3.156 and 3.027 Å. The results of studies of the stoichiometry and formation of complex 2 in methanol solution were found to support its solid state stoichiometry.     

Atomistic understanding of the C·T mismatched DNA base pair tautomerization via the DPT: QM and QTAIM computational approaches

It was established that the cytosine·thymine (C·T) mismatched DNA base pair with cis‐oriented N1H glycosidic bonds has propeller‐like structure (|N3C4C4N3| = 38.4°), which is stabilized by three specific intermolecular interactions–two antiparallel N4H…O4 (5.19 kcal mol^?1) and N3H…N3 (6.33 kcal mol^?1) H‐bonds and a van der Waals (vdW) contact O2…O2 (0.32 kcal mol^?1). The C·T base mispair is thermodynamically stable structure (ΔG_int = ?1.54 kcal mol^?1) and even slightly more stable than the A·T Watson–Crick DNA base pair (ΔG_int = ?1.43 kcal mol^?1) at the room temperature. It was shown that the C·T ? C*·T* tautomerization via the double proton transfer (DPT) is assisted by the O2…O2 vdW contact along the entire range of the intrinsic reaction coordinate (IRC). The positive value of the Grunenberg's compliance constants (31.186, 30.265, and 22.166 Å/mdyn for the C·T, C*·T*, and TS_{C·T ? C*·T*}, respectively) proves that the O2…O2 vdW contact is a stabilizing interaction. Based on the sweeps of the H‐bond energies, it was found that the N4H…O4/O4H…N4, and N3H…N3 H‐bonds in the C·T and C*·T* base pairs are anticooperative and weaken each other, whereas the middle N3H…N3 H‐bond and the O2…O2 vdW contact are cooperative and mutually reinforce each other. It was found that the tautomerization of the C·T base mispair through the DPT is concerted and asynchronous reaction that proceeds via the TS_{C·T ? C*·T*} stabilized by the loosened N4? H? O4 covalent bridge, N3H…N3 H‐bond (9.67 kcal mol^?1) and O2…O2 vdW contact (0.41 kcal mol^?1). The nine key points, describing the evolution of the C·T ? C*·T* tautomerization via the DPT, were detected and completely investigated along the IRC. The C*·T* mispair was revealed to be the dynamically unstable structure with a lifetime 2.13·× 10^?13 s. In this case, as for the A·T Watson–Crick DNA base pair, activates the mechanism of the quantum protection of the C·T DNA base mispair from its spontaneous mutagenic tautomerization through the DPT. © 2013 Wiley Periodicals, Inc.     

Investigating hydrogen bonds in liquid ethylene glycol structure by means of molecular dynamics

Models of liquid ethylene glycol are built by means of molecular dynamics at temperatures ranging between 268 and 443 K, with 1000 molecules in rectangular parallelepiped basic cells. The dependences of structures of O-H…O hydrogen bonds on modeling time and temperature are analyzed. It is found that the hydrogen bonds emerge at different sites of a model, thus forming a hydrogen bonds network that is continuously rebuilt under the action of thermal fluctuations. The number of hydrogen bonds in the models is observed to decrease when the temperature is raised. The energy of hydrogen bond formation is found to be ?20.0 ± 2.6 kJ mol^?1, the average bond lifetime is 370 ps at 268 K and 147 ps at 323 K, and the activation energy of hydrogen bond rupture at these temperatures is ～12.1 kJ mol^?1. It is concluded that the data on the breaking of H-bonds at temperatures of 323 to 443 K can be explained by the molecules moving away from each other as a result of diffusive motion, accompanied by rearrangement of the hydrogen bonds network. The concentration of dimers in the models is shown to be rather low, while the average energy of forming a dimer from two ethylene glycol molecules is ?35.4 kJ mol^?1.     

Nature of MoH···I bonds in Cp2Mo(L)H···I‐C≡C‐R Complexes (L=H,CN, PPh2, C(CH3)3; R=NO2, Cl,Br, H,OH, CH3, NH2)

The nature of the MoH···I bond in Cp₂Mo(L)H···I‐C≡C‐R (L= H, CN, PPh₂, C(CH₃)₃; R=NO₂, Cl, Br, H, OH, CH₃, NH₂) was investigated using electrostatic potential analysis, topological analysis of the electron density, energy decomposition analysis and natural bond orbital analysis. The calculated results show that MoH···I interactions in the title complexes belong to halogen‐hydride bond, which is similar to halogen bonds, not hydrogen bonds. Different to the classical halogen bonds, the directionality of MoH···I bond is low; Although electrostatic interaction is dorminant, the orbital interactions also play important roles in this kind of halogen bond, and steric interactions are weak; the strength of H···I bond can tuned by the most positive electrostatic potential of the I atom. As the electron‐withdrawing ability of the R substituent in the alkyne increases, the electrostatic potential maximum of the I atom increases, which enhances the strength of the H···I halogen bond, as well as the electron transfer.     

Van Der Waals heterogeneous layer‐layer carbon nanostructures involving π···H‐C‐C‐H···π···H‐C‐C‐H stacking based on graphene and graphane sheets

Noncovalent interactions involving aromatic rings, such as π···π stacking, CH···π are very essential for supramolecular carbon nanostructures. Graphite is a typical homogenous carbon matter based on π···π stacking of graphene sheets. Even in systems not involving aromatic groups, the stability of diamondoid dimer and layer‐layer graphane dimer originates from C − H···H − C noncovalent interaction. In this article, the structures and properties of novel heterogeneous layer‐layer carbon‐nanostructures involving π···H‐C‐C‐H···π···H‐C‐C‐H stacking based on [n ]‐graphane and [n ]‐graphene and their derivatives are theoretically investigated for n = 16–54 using dispersion corrected density functional theory B3LYP‐D3 method. Energy decomposition analysis shows that dispersion interaction is the most important for the stabilization of both double‐ and multi‐layer‐layer [n ]‐graphane@graphene. Binding energy between graphane and graphene sheets shows that there is a distinct additive nature of CH···π interaction. For comparison and simplicity, the concept of H‐H bond energy equivalent number of carbon atoms (noted as NHEQ), is used to describe the strength of these noncovalent interactions. The NHEQ of the graphene dimers, graphane dimers, and double‐layered graphane@graphene are 103, 143, and 110, indicating that the strength of C‐H···π interaction is close to that of π···π and much stronger than that of C‐H···H‐C in large size systems. Additionally, frontier molecular orbital, electron density difference and visualized noncovalent interaction regions are discussed for deeply understanding the nature of the C‐H···π stacking interaction in construction of heterogeneous layer‐layer graphane@graphene structures. We hope that the present study would be helpful for creations of new functional supramolecular materials based on graphane and graphene carbon nano‐structures. © 2017 Wiley Periodicals, Inc.     

Strong N−X⋅⋅⋅O−N Halogen Bonds: A Comprehensive Study on N‐Halosaccharin Pyridine N‐Oxide Complexes

A study of the strong N?X???^?O?N⁺ (X=I, Br) halogen bonding interactions reports 2×27 donor×acceptor complexes of N‐halosaccharins and pyridine N‐oxides (PyNO). DFT calculations were used to investigate the X???O halogen bond (XB) interaction energies in 54 complexes. A simplified computationally fast electrostatic model was developed for predicting the X???O XBs. The XB interaction energies vary from ?47.5 to ?120.3 kJ mol^?1; the strongest N?I???^?O?N⁺ XBs approaching those of 3‐center‐4‐electron [N?I?N]⁺ halogen‐bonded systems (ca. 160 kJ mol^?1). ¹H NMR association constants (K_XB) determined in CDCl₃ and [D₆]acetone vary from 2.0×10⁰ to >10⁸ m ^?1 and correlate well with the calculated donor×acceptor complexation enthalpies found between ?38.4 and ?77.5 kJ mol^?1. In X‐ray crystal structures, the N‐iodosaccharin‐PyNO complexes manifest short interaction ratios (R_XB) between 0.65–0.67 for the N?I???^?O?N⁺ halogen bond.     

Structural and Solution Studies of two Mercury(II) Complexes with [1,3‐Bis(2‐ethoxy)benzene]triazene

The reaction of [1,3‐bis(2‐ethoxy)benzene]triazene, [ HL ], with Hg(SCN)₂ and Hg(CH₃COO)₂, resulted in the formation of the complexes [Hg L (SCN)] ( 1 ) and [Hg L ₂] · CH₃OH ( 2 ). They were characterized by means of X‐ray crystallography, CHN analysis, FT‐IR, ¹H NMR, and ¹³C NMR spectroscopy. The structure of compound 1 consists of two independent complexes in which the Hg^II atoms are stacked along the crystallographic a axis to form infinite chains. Each Hg^II atom is chelated by one L ligand and one SCN ligand, whereas in compound 2 , the Hg^II atom is surrounded by two L ligands. In addition, 1D chains formed by metal–π interactions are connected to each other by C–H ··· π stacking interactions in the structure of 1 , which results in a 2D architecture. An interesting feature of compound 2 is the presence of C–H ··· π edge‐to‐face interactions.     

A theoretical study of 1:1 and 1:2 complexes of acetylene with nitrosyl hydride

Ab initio calculations at MP2 computational level using aug-cc-pVTZ basis set were used to analyze the interactions between 1:1 and 1:2 complexes of acetylene and nitrosyl hydride. The structures obtained have been analyzed with the atoms in molecules and the density functional theory–symmetry adapted perturbation theory methodologies. Four minima were located on the potential energy surface of the 1:1 complex. Twenty-four different structures have been obtained for the 1:2 complexes. Five types of interactions are observed, CH···O, CH···N, NH···π hydrogen bonds and orthogonal interactions between the π clouds of triple bond, or the lone pair of oxygen with the electron-deficient region of the nitrogen atom. Stabilization energies of the 1:1 and 1:2 clusters including basis set superposition error and ZPE are in the range 3–8 and 6–17 kJ mol⁻¹ at MP2/aug-cc-pVTZ computational level, respectively. Blue shift of NH bond upon complex formation in the ranges between 18–30 and 20–96 cm⁻¹ is predicted for 1:1 and 1:2 clusters, respectively. The total nonadditive energy in the 1:2 cluster, calculated as the sum of the supermolecular nonadditive MP2 energy and the three-body dispersion energy, presents values between −1.48 and 1.20 kJ mol⁻¹.     

Unexpected A·T(WC)↔A·T(rWC)/A·T(rH) and A·T(H)↔A·T(rH)/A·T(rWC) conformational transitions between the classical A·T DNA base pairs: A QM/QTAIM comprehensive study

In this study, using QM/QTAIM calculations in the continuum with ε = 1 under normal conditions, we have revealed for the first time the nondissociative A·T(WC)↔A·T(rWC)/A·T(rH) and A·T(H)↔A·T(rH)/A·T(rWC) conformational transitions. It was established that they proceed via the essentially nonplanar transition states (С₁ symmetry) through the intermediates, which are wobbled conformers (С₁ symmetry) theoretically predicted in our previous work (Brovarets’ et al., Frontiers in Chemistry, 2018, 6:8, 10.3389/fchem.2018.00008) of the classical А·Т DNA base pairs—Watson–Crick А·Т(WC), reverse Watson–Crick А·Т(rWC), Hoogsteen А·Т(Н) and reverse Hoogsteen А·Т(rН). At this, the A·T(H)↔A·T(rWC) and A·T(WC)↔A·T(rH) conformational transformations are controlled by the transition states (TSs) stabilized by the participation of the intermolecular (T)N3H···N6(A) H‐bond (∼3.70 kcal·mol⁻¹) between the imino group N3H of T and pyramidilized amino group N6H₂ of A. Gibbs free energies of activation for these processes consist 12.22 and 11.11 kcal·mol⁻¹, accordingly, under normal conditions. TSs, which control the A·T(WC)↔A·T(rWC) and A·T(H)↔A·T(rH) conformational transitions are stabilized by the participation of the intermolecular (T)N3H···N6(A) H‐bond (5.82 kcal·mol⁻¹) and bifurcating intermolecular (T)N3H···N6(A) (5.00) and (T)N3H···N7(A) (0.61 kcal·mol⁻¹) H‐bonds, accordingly. Notably, in these two TSs amino group N6H₂ of A is significantly pyramidilized; Gibbs free energies of activation for these reactions are 19.07 and 19.71 kcal·mol⁻¹, accordingly.     

Determination of the acid dissociation constants of p-benzohydroquinone by the INDO method

The stepwise acid dissociation constants for p-benzohydroquinone (QH₂) in aqueous media have been explicitly calculated for the first time, with the INDO parametrized SCF –MO method. We have optimized the geometries of QH₂, QH^?, and Q^2? and of the QH₂ · 6H₂O, QH^? · (H₃O⁺) · 5H₂O, and Q^2? · (H₃O⁺)₂ · 4H₂O systems that model the solvated species. The presence of the associated water molecules (and hydronium ions) account for the stabilization due to hydrogen bonding as well as for a part of the effect of interaction of these molecules with the respective reaction fields in an aqueous medium. To simulate the first solvation shell in a more complete manner, four more water molecules have been considered to be placed above and below the quinonoid ring and the optimized geometries of the resulting hydrated species, QH₂ · 10H₂O, QH^? · (H₃O⁺) · 9H₂O, and QH^? · (H₃O⁺) · 8H₂O, have been determined. The standard free-energy changes calculated for the dissociation of QH₂ into QH^? and H⁺ is 0.0251 Hartree (65.9 kJ mol^?1) and that of QH^? into Q^2? and H⁺ is 0.0285 Hartree (74.8 kJ mol^?1). Experimentally observed dissociation constants for these two steps correspond to free-energy changes of 0.0214 Hartree (56.2 kJ mol^?1) and 0.0248 Hartree (65.1 kJ mol^?1), respectively. © 1995 John Wiley & Sons, Inc.     

Synthesis and Characterization of Iron(III) Complexes of 5‐(8‐Carboxy‐1‐naphthyl)‐10, 15, 20‐tritolyl Porphyrin

5‐(8‐Carboxy‐1‐naphthyl)‐10, 5, 20‐tritolyl porphyrin (H₃CNTTP) and its iron(III) complexes, [Fe(CNTTP)]₂ and [Fe(CNTTP)(N‐MeIm)₂], were synthesized and characterized. X‐ray crystallography revealed that the carboxylate group is “hanging” over the porphyrin plane. The rigid framework makes the distance between the carboxylate oxygen and iron in the same porphyrin too long to form a coordination bond. On the other hand, the carboxylate group is not bulky enough to block the axial binding site. In the presence of OH^–, the carboxylate oxygen is coordinated to iron in the symmetry‐related unit, which led to the dimeric structure, [Fe(CNTTP)]₂. In the presence of excess N‐methylimidazole, a six‐coordinate species, [Fe(CNTTP)(N‐MeIm)₂], was obtained. In such a structure, CH ··· O interactions between the carboxylate group and imidazole probably play an important role to determine the orientation of imidazole plane. Two imidazole planes have relative parallel orientation. For [Fe(CNTTP)(N‐MeIm)₂], ¹H NMR shows pyrrole protons at the region –10 to –25 ppm. EPR shows rhombic spectrum. Those suggest [Fe(CNTTP)(N‐MeIm)₂] is a type II low‐spin iron(III) porphyrinate.     

Influence of the C—H· · ·N intramolecular interaction on the spatial structures and 1H and 13C NMR parameters of heteroaryl vinyl ethers and sulfides

A complete analysis of the ¹H and ¹³C spectra of the representative series of heteroaryl vinyl ethers and sulfides and heteroaryl styryl sulfides was carried out. The electronic and spatial structures of these compounds are discussed. It was shown that the C—H· · ·N intramolecular interactions in the investigated molecules influence significantly the spectral parameters and the conformational equilibrium. Copyright © 2003 John Wiley & Sons, Ltd.     

A computational study of 1:1 and 1:2 complexes of nitryl halides (O2NX) with HCN and HNC

Quantum calculations at the MP2/cc-pVTZ have been used to examine 1:1 and 1:2 complexes between O₂NX (X = Cl, Br) with HCN and HNC moieties. The interaction of the lone pair of the HCN(HNC) with the σ-hole and π-hole of O₂NX molecules and hydrogen bonding between lone pairs of X and O of O₂NX with H of HCN and HNC have been considered in 1:1 complexes. The 1:1 complexes can easily be differentiated using the stretching frequency of the N–X bond. Thus, those complexes with σ-hole and H···O₂NX interactions show a blue shift of the N-X bond stretching while a red shift is observed in the complexes along the π-hole and H···XNO₂ interactions. In the 1:2 complexes, the cooperative and diminutive energetic effects have been analyzed using the many-body interaction energies. The nature of the interactions has been characterized with the Atoms in Molecules (AIM) and Natural Bond Orbital (NBO) methodologies. Stabilization energies of 1:1 and 1:2 complexes including the variation of the zero-point vibrational energy (ΔZPVE) are in the range 3–9 kJ mol^?1 and 21–40 kJ mol^?1, respectively.     

New Fluoromanganate(III) Hydrates: Mn3F8 · 12H2O and AgMnF4 · 4H2O

Single crystals of fluoride hydrates Mn₃F₈ · 12 H₂O and AgMnF₄ · 4 H₂O have been prepared and characterized by X-ray methods. Mn₃F₈ · 12 H₂O crystallizes in the space group P1 (a = 623.0(3), b = 896.7(4), c = 931.8(4) pm, α = 110.07(2)°, β = 103.18(2)°, γ = 107.54(2)°, Z = 1); AgMnF₄ · 4 H₂O crystallizes in the space group P2₁/m (a = 700.9(2), b = 726.1(1), c = 749.4(3) pm, β = 107.17(3)°, Z = 2). Both structures contain Jahn-Teller-distorted [Mn(H₂O)₂F₄]^? anions as well as crystal water molecules and exhibit a complex hydrogen bond network between anions and cations, i. e. [Mn(H₂O)₆]²⁺ for the first and a polymeric [Ag(H₂O)₂]^? cation for the second compound.     

Different catalysis role of in‐loop and out‐of‐loop waters in assisting HNS/HSN proton transfer isomerizations: Bridging vs. surrounding effect

In this work, a density function theory (DFT) study is presented for the HNS/HSN isomerization assisted by 1–4 water molecules on the singlet state potential energy surface (PES). Two modes are considered to model the catalytic effect of these water molecules: (i) water molecule(s) participate directly in forming a proton transfer loop with HNS/HSN species, and (ii) water molecules are out of loop (referred to as out‐of‐loop waters) to assist the proton transfer. In the first mode, for the monohydration mechanism, the heat of reaction is 21.55 kcal · mol^?1 at the B3LYP/6‐311++G** level. The corresponding forward/backward barrier lowerings are obtained as 24.41/24.32 kcal · mol^?1 compared with the no‐water‐assisting isomerization barrier T (65.52/43.87 kcal · mol^?1). But when adding one water molecule on the HNS, there is another special proton‐transfer isomerization pathway with a transition state 10T′ in which the water is out of the proton transfer loop. The corresponding forward/backward barriers are 65.89/65.89 kcal · mol^?1. Clearly, this process is more difficult to follow than the R–T–P process. For the two‐water‐assisting mechanism, the heat of reaction is 19.61 kcal · mol^?1, and the forward/backward barriers are 32.27/12.66 kcal · mol^?1, decreased by 33.25/31.21 kcal · mol^?1 compared with T. For trihydration and tetrahydration, the forward/backward barriers decrease as 32.00/12.60 (30T) and 37.38/17.26 (40T) kcal · mol^?1, and the heat of reaction decreases by 19.39 and 19.23 kcal · mol^?1, compared with T, respectively. But, when four water molecules are involved in the reactant loop, the corresponding energy aspects increase compared with those of the trihydration. The forward/backward barriers are increased by 5.38 and 4.66 kcal · mol^?1 than the trihydration situation. In the second mode, the outer‐sphere water effect from the other water molecules directly H‐bonded to the loop is considered. When one to three water molecules attach to the looped water in one‐water in‐loop‐assisting proton transfer isomerization, their effects on the three energies are small, and the deviations are not more than 3 kcal · mol^?1 compared with the original monohydration‐assisting case. When adding one or two water molecules on the dihydration‐assisting mechanism, and increasing one water molecule on the trihydration, the corresponding energies also are not obviously changed. The results indicate that the forward/backward barriers for the three in‐loop water‐assisting case are the lowest, and the surrounding water molecules (out‐of‐loop) yield only a small effect. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2006     

A DFT study of Lp···π/halogen bond competition in complexes of perhalogenated alkenes with oxygen/nitrogen containing simple molecules

The possible noncovalent lone pair‐π/halogen bond (lp···π/HaB) complexes of perhalogenated unsaturated C₂Cl_nF_4?n (n = 0–4) molecules with four simple molecules containing oxygen or nitrogen as electron donor, formaldehyde (H₂CO), dimethyl ether (DME), NH₃, and trimethylamine (TMA), have been systematically examined at the M062X/aug‐cc‐pVTZ level. Natural bond orbital (NBO) analysis at the same level is used for understanding the electron density distributions of these complexes. The progressive introduction of Cl atom on C₂Cl_nF_4?n influences more on the lp···π complexes over the corresponding HaB ones. Within the scope of this study, gem‐C₂Cl₂F₂ is the best partner molecule for lp···π interaction with the simple molecules, coupled with the greatest interaction energy (IE) and second‐order orbital interaction [E(2) value], whereas C₂F₄ is the poorest one. The C₂Cl₃F·H₂CO and C₂Cl₄·H₂CO complexes exhibit reverse lp···π bonding, while the Z/E‐C₂Cl₂F₂·NH₃, C₂Cl₃F·NH₃ and C₂Cl₄·NH₃ complexes perform half‐lp···π bonding according to the NBO analysis. The lp···π interaction involving the oxygen/nitrogen and the π‐hole of C₂Cl_nF_4?n overwhelms the HaB involving the oxygen/nitrogen and the σ‐hole of the Cl atom. The electron‐donating methyl groups contribute significantly to the two competitive interactions, therefore, DME and TMA engage stronger in the partner molecules than H₂CO and NH₃. Our theoretical study would be useful for future experimental investigation on noncovalent complexes. © 2016 Wiley Periodicals, Inc.     

Global and local interactions in the structure of crystalline 7-(diethylamino)-2-(2-oxo-2<Emphasis Type="Italic">H</Emphasis>-chromen-3-yl)chromenium perchlorate

Computational methods were used to calculate the crystal lattice energy reflecting global interactions, predominantly long-range electrostatic interactions between ions, as well as the energy of selected specific local C–H···O, C–H···π and π···π interactions found in synthesized 7-(diethylamino)-2-(2-oxo-2H-chromen-3-yl)chromenium perchlorate, the structure of which was determined by X-ray crystallography. Local interactions occurring between specific sites of molecules, amounting to a few tens of kJ mol^?1, most likely account for the mutual arrangement of molecular ions, whereas global ones, exceeding half-a-thousand kJ mol^?1, are responsible for the thermodynamic stability of the compound investigated in the crystalline solid phase, whose potential applications are briefly outlined.