A theoretical investigation on the conformation and the interaction of CHF2OCF2CHF2 (desflurane II) with one water molecule

The conformation and the interaction of CHF₂OCF₂CHF₂ (desflurane II) with one water molecule is investigated theoretically using the ab initio MP2/aug-cc-pvdz and DFT-based M062X/6-311++G(d,p) methods. The calculations include the optimized geometries, the harmonic frequencies of relevant vibrational modes along with a natural bond orbital (NBO) analysis including the NBO charges, the hybridization of the C atom and the intra- and intermolecular hyperconjugation energies. In the two most stable conformers, the CH bond of the F₂HCO- group occupies the gauche position. The hyperconjugation energies are about the same for both conformers and the conformational preference depends on the interaction between the non-bonded F and H atoms. The deprotonation enthalpies of the CH bonds are about the same for both conformers, the proton affinity of the less stable conformer being 3 kcal mol^?1 higher. Both conformers of desflurane II interact with water forming cyclic complexes characterized by CH…O and OH…F hydrogen bonds. The binding energies are moderate, ranging from ?2.4 to ?3.2 kcal mol^?1 at the MP2 level. The origin of the blue shifts of the ν(CH) vibrations is analyzed. In three of the complexes, the water molecule acts as an electron donor. Interestingly, in these cases a charge transfer is also directed to the non bonded OH group of the water molecule. This effect seems to be a property of polyfluorinated ethers.     

A new era of the prototropic tautomerism of the quercetin molecule: A QM/QTAIM computational advances

Abstract

In this study for the first time we have revealed and investigated in details 123 different prototropic tautomers of the most stable conformer of the quercetin molecule using quantum-mechanical calculations at the MP2/6-311++G(2df,pd)//B3LYP/6-311++G(d,p) level of QM theory. We have found that in the most energetically favorable prototropic tautomer mobile hydrogen atoms are localized at the О3, О3′, О4′, О5, and О7 exocyclic oxygen atoms. Molecular tautomers are in the range of the Gibbs free energies from 0.0 to 69.8?kcal·mol^?1, while zwitterionic ones – from 30.1 до 172.8?kcal·mol^?1 at normal conditions. It was also reliably established that the weakest point causing the decyclization of the molecule is its C ring – this reaction is launched by the transition of the proton from the C8H group to the endocyclic O1 oxygen atom. All prototropic tautomers, except two cases, are joined by the intramolecular cooperative specific interactions (from 1 to 5) – H-bonds and attractive van der Waals contacts, which have been revealed and characterized by QTAIM analysis.

Communicated by Ramaswamy H. Sarma     

Zwitterionic conformers of pyrrolysine and their interactions with metal ions—a theoretical study

A total of 16 pyrrolysine conformers in their zwitterionic forms are studied in gas and simulated aqueous phase using a polarizable continuum model (PCM). These conformers are selected on the basis of our study on the intrinsic conformational properties of non-ionic pyrrolysine molecule in gas phase [Das and Mandal (2013) J Mol Model 19:1695?1704]. In aqueous phase, the stable zwitterionic pyrrolysine conformers are characterized by full geometry optimization and vibrational frequency calculations using B3LYP/6-311++G(d,p) level of theory. Single point calculations are also carried out at MP2/6-311++G(d,p) level. Characteristic intramolecular hydrogen bonds present in each conformer, their relative energies, theoretically predicted vibrational spectra, rotational constants and dipole moments are systematically reported. The calculated relative energy range of the conformers at B3LYP/6-311++G(d,p) level is 5.19 kcal mol^?1 whereas the same obtained by single point calculations at MP2/6-311++G(d,p) level is 4.58 kcal mol^?1. A thorough analysis reveals that four types of intramolecular H-bonds are present in the conformers; all of which play key roles in determining the energetics and in imparting the observed conformations to the conformers. The vibrational frequencies are found to shift invariably toward the lower side of frequency scale corresponding to the presence of the H-bonds. This study also points out that conformers with diverse structural motifs may differ in their thermodynamical stability by a narrow range of relative energy. The effects of metal coordination on the relative stability order and structural features of the conformers are examined by complexing five zwitterionic conformers of pyrrolysine with Cu⁺² through their carboxylate groups. The interaction enthalpies and Gibbs energies, rotational constants, vibrational frequencies and dipole moments of the metal complexes calculated at B3LYP level are also reported. The zwitterionic conformers of pyrrolysine are not stable in gas phase; after geometry optimization they are converted to the non-ionic forms.     

Conformational transitions of the quercetin molecule via the rotations of its rings: a comprehensive theoretical study

Abstract

Quercetin is an important flavonoid compound, usually extracted from plants, vegetables and fruits such as blueberries, apples, green tea, wine, onions and possessing broad range of pharmacological properties, in particular, powerful antioxidant, antitoxic, antiinflammation and antimicrobial effects due to its various reactive sites. The structure of this phenolic compound consists of three (A?+?C) and B rings, bearing five hydroxyl groups. Primarily, the chemical structure of quercetin determines its physico-chemical properties. Earlier, it was established that isolated quercetin molecule can acquire 48 stable conformations (24 planar and 24 non-planar) due to the mobility of its hydroxyl groups and (A?+?C) and B rings with relative Gibbs free energies in the range of 0.0–25.3?kcal·mol^?1 under normal conditions (Brovarets’ et al., 2019c Brovarets’, O. O., & Hovorun, D. M. (2019c). Conformational diversity of the quercetin molecule: A quantum-chemical view. Journal of Biomolecular Structure and Dynamics. doi: 10.1080/07391102.2019.1645734[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]). In this work by quantum-mechanical calculations at the MP2/6-311++G(2df,pd)//B3LYP/6-311++G(d,p) level of theory and Bader’s ‘Quantum Theory of Atoms in Molecules’, we have theoretically modeled the interconversions in the 24 pairs of the conformers of the quercetin molecule, occuring via the rotation of its non-deformable (A+С) and B rings around the С2-С1' bond through the quasi-orthogonal transition state with low values of the imaginary frequencies (28–33/29–36?cm^?1) and Gibbs free energies of activation in the range of 2.17–5.68/1.86–4.90?kcal·mol^?1 in the continuum with dielectric permittivity ε?=?1/ε?=?4 under normal conditions. Also, we studied the changes of the number of physico-chemical characteristics of all intramolecular-specific contacts – hydrogen bonds and attractive van der Waals contacts during these conformational rearrangements.

Communicated by Ramaswamy H. Sarma     

Quantum-chemical and empirical calculations on phospholipids. III. Hydration of the dimethylphosphate anion

The binding of water to the dimethylphosphate anion (DMP^?) was calculated using the PCILO method. We found binding energies of 25.95 kcal·mol^?1 in the O1-P-O3 plane bridging the anionic oxygen atoms and 19.3 kcal·mol^?1 for the one-site association of a water molecule to an anionic oxygen atom of DMP^?. In this range one water molecule added to DMP^? in the O1 … O3 bridged configuration has a significantly higher binding energy to DMP^? than water molecules added to other binding sites. The total binding energy of 5 water molecules to DMP^? is 92 kcal·mol^?1, a quantity which is about 10% less than the sum of the binding energies of the corresponding monohydrates.     

Anesthetic-protein interaction Random versus helix polylysine monolayers and interaction with 1-alkanols

Penetration of 1-alkanols into monolayers of hydrophobic polypeptides, poly(ε-benzyloxycarbonyl-l-lysine) and poly(ε-benzyloxycarbonyl-dl-lysine), was compared with their adsorption on the air/water interface in the absence of monolayers. The polypeptide prepared from l-lysine is generally considered to be in the α-helical form whereas dl-copolymer polypeptide contains random-coiled portions due to the structural incompatibility between the two isomers. The free energy of adsorption of 1-alkanols on the air/water interface at dilute concentrations was ?0.68 kcal·mol^?1 per methylene group and 0.15 kcal·mol^?1 for the hydroxyl group at 25°C. In the close-packed state, the surface area occupied by each molecule of 1-alkanols of varying carbon chain-lengths showed nearly a constant value of about 27.2 Ų, indicating perpendicular orientation of the alkanol molecules at the interface. About 75% of the water surface was covered by 1-butanol in this close-packed state. The mode of adsorption of 1-alkanols on the vacant air/water interface followed the Gibbs surface excess while the mode on the polypeptide membranes followed the Langmuir adsorption isotherm, indicating that the latter is characterized by the presence of a finite number of binding sites. The free energies of adsorption of 1-alkanols on the l-polymer monolayers were more negative than those on the vacant air/water interface and less negative than those on the dl-copolymer monolayers. Thus, the affinity of 1-alkanols to the interface was in the order of vacant air/water interface <l-polymer <dl-copolymer. The difference between the air/water interface and l-polymer was about 0.54 kcal·mol^?1 and that between l-polymer and dl-copolymer was 0.17 kcal·mol^?1 at 25°C: the adsorption of 1-alkanols to the dl-copolymer was favored compared to the l-polymer. The polar moieties of the backbone of the dl-copolymer may be exposed to the aqueous phase at the disordered portion. Dipole interaction between this portion and 1-alkanol molecules may account for the enhanced adsorption of the alkanols to the dl-copolymer.     

DFT studies on the intrinsic conformational properties of non-ionic pyrrolysine in gas phase

B3LYP/6-31G(d,p) level of theory is used to carry out a detailed gas phase conformational analysis of non-ionized (neutral) pyrrolysine molecule about its nine internal back-bone torsional angles. A total of 13 minima are detected from potential energy surface exploration corresponding to the nine internal back-bone torsional angles. These minima are then subjected to full geometry optimization and vibrational frequency calculations at B3LYP/6-31++G(d,p) level. Characteristic intramolecular hydrogen bonds present in each conformer, their relative energies, theoretically predicted vibrational spectra, rotational constants and dipole moments are systematically reported. Single point calculations are carried out at B3LYP/6-311++G(d,p) and MP2/6-31++G(d,p) levels. Six types of intramolecular H-bonds, viz. O…H–O, N…H-O, O…H–N, N…H–N, O…H–C and N…H–C, are found to exist in the pyrrolysine conformers; all of which contribute to the stability of the conformers. The vibrational frequencies are found to shift invariably toward the lower side of frequency scale corresponding to the presence of intramolecular H-bond interactions in the conformers.     

DPT tautomerization of the long A∙A* Watson-Crick base pair formed by the amino and imino tautomers of adenine: combined QM and QTAIM investigation

Combining quantum-mechanical (QM) calculations with quantum theory of atoms in molecules (QTAIM) and using the methodology of sweeps of the energetic, electron-topological, geometric and polar parameters, which describe the course of the tautomerization along the intrinsic reaction coordinate (IRC), we showed for the first time that the biologically important A?A* base pair (C_s symmetry) formed by the amino and imino tautomers of adenine (A) tautomerizes via asynchronous concerted double proton transfer (DPT) through a transition state (TS), which is the A⁺?A^? zwitterion with the separated charge, with C_s symmetry. The nine key points, which can be considered as electron-topological “fingerprints” of the asynchronous concerted A?A*?A*?A tautomerization process via the DPT, were detected and completely investigated along the IRC of the A?A*?A*?A tautomerization. Based on the sweeps of the H-bond energies, it was found that intermolecular antiparallel N6Н?N6 (7.01 kcal mol^?1) and N1H?N1 (6.88 kcal mol^?1) H-bonds are significantly cooperative and mutually reinforce each other. It was shown for the first time that the A?A*?A*?A tautomerization is assisted by the third C2H?HC2 dihydrogen bond (DHB), which, in contrast to the two others N6H?N6 and N1H?N1 H-bonds, exists within the IRC range from ?2.92 to 2.92 Å. The DHB cooperatively strengthens, reaching its maximum energy 0.42 kcal mol^?1 at IRC?=??0.52 Å and minimum energy 0.25 kcal mol^?1 at IRC?=??2.92 Å, and is accompanied by strengthening of the two other aforementioned classical H-bonds. We established that the C2H?HC2 DHB completely satisfies the electron-topological criteria for H-bonding, in particular Bader’s and all eight “two-molecule” Koch and Popelier’s criteria. The positive value of the Grunenberg’s compliance constant (5.203 Å/mdyn) at the TS_A?A*?A*?A proves that the C2H?HC2 DHB is a stabilizing interaction. NBO analysis predicts transfer of charge from σ(C2–H) bonding orbital to σ*(H–C2) anti-bonding orbital; at this point, the stabilization energy E⁽²⁾ is equal to 0.19 kcal mol^?1 at the TS_A?A*?A*?A.     

Theoretical studies of the interaction between enflurane and water

Increase of the atmospheric concentration of halogenated organic compounds is partially responsible for a change of the global climate. In this work we have investigated the interaction between halogenated ether and water, which is one of the most important constituent of the atmosphere. The structures of the complexes formed by the two most stable conformers of enflurane (a volatile anaesthetic) with one and two water molecules were calculated by means of the counterpoise CP-corrected gradient optimization at the MP2/6–311++G(d,p) level. In these complexes the CH…O_w hydrogen bonds are formed, with the H…O_w distances varying between 2.23 and 2.32 Å. A small contraction of the CH bonds and the blue shifts of the ν(CH) stretching vibrations are predicted. There is also a weak interaction between one of the F atoms and the H atom of water, with the H_w…F distances between 2.41 and 2.87 Å. The CCSD(T)/CBS calculated stabilization energies in these complexes are between ?5.89 and ?4.66 kcal?mol^?1, while the enthalpies of formation are between ?4.35 and ?3.22 kcal?mol^?1. The Cl halogen bonding between enflurane and water has been found in two complexes. The intermolecular (Cl···O) distance is smaller than the sum of the corresponding van der Waals radii. The CCSD(T)/CBS stabilization energies for these complexes are about ?2 kcal?mol^?1.

Molecular dynamics simulations of void defects in the energetic material HMX

A molecular dynamics (MD) simulation was carried out to characterize the dynamic evolution of void defects in crystalline octahydro-1, 3, 5, 7-tetranitro-1, 3, 5, 7-tetrazocine (HMX). Different models were constructed with the same concentration of vacancies (10 %) to discuss the size effects of void. Energetic ground state properties were determined by annealing simulations. The void formation energy per molecule removed was found to be 55–63 kcal/mol^?1, and the average binding energy per molecule was between 32 and 34 kcal/mol^?1 according to the change in void size. Voids with larger size had lower formation energy. Local binding energies for molecules directly on the void surface decreased greatly compared to those in defect-free lattice, and then gradually increased until the distance away from the void surface was around 10 Å. Analysis of 1 ns MD simulations revealed that the larger the void size, the easier is void collapse. Mean square displacements (MSDs) showed that HMX molecules that had collapsed into void present liquid structure characteristics. Four unique low-energy conformers were found for HMX molecules in void: two whose conformational geometries corresponded closely to those found in HMX polymorphs and two, additional, lower energy conformers that were not seen in the crystalline phases. The ratio of different conformers changed with the simulated temperature, in that the ratio of α conformer increased with the increase in temperature.     

Base Pairing at the Stem-loop Junction in the SL1 Kissing Complex of HIV-1 RNA: A Thermodynamic Study Probed by Molecular Dynamics Simulation

Abstract

The thermodynamics of the opening/closure process of a GC base pair located at the stem-loop junction of the SL1 sequence from HIV-1_Lai genomic RNA was investigated in the context of a loop-loop homodimer (or kissing complex) using molecular dynamics simulation. The free energy, enthalpy and entropy changes for the closing reaction are 0 kcal·mol^?1, ?11 kcal·mol^?1and ?0.037 kcal·mol^?1-K^?1 at 300° K respectively. Furthermore it is found that the free energy change is the same for the formation of a 11 nucleotide loop closed with UG and for the formation of a 9 nucleotide loop closed with GC. Our study evidences the high flexibility of the nucleotides at the stem-loop junction explaining the weak stability of this structure.     

Quercetin and daidzein β-apo-14’-carotenoic acid esters as membrane antioxidants

Abstract

Esterification by β-apo-14’-carotenoic acid was found to have opposite effects on antioxidant activity of quercetin (at B4’, B3’ hydroxyl) as of daidzein (at A7 hydroxyl) in phosphatidylcholine liposomes. The daidzein ester had increased activity, while quercetin had a significant decreased activity. Quantum mechanical calculations using density function theory (DFT) indicate a modest decrease in bond dissociation enthalpy, BDE, for (weakest) hydrogen–oxygen phenolic bond in daidzein from 368.4 kJ·mol^{? 1} to 367.7 kJ·mol^{? 1} compared to a significant increase in quercetin from 329.5 kJ·mol^{? 1} to 356.6 kJ·mol^{? 1} upon derivatization. These opposite changes in tendency for hydrogen atom transfer from phenolic groups to lipid radicals combined with an increase in A-to-B dihedral angle from 0.0° to 36.4° and in dipole moment from 0.40 D to 6.01 D for quercetin upon derivatization, while less significant for daidzein (36.4°–36.7° and 3.26 D–7.87 D, respectively), together provide a rationale for the opposite effect of esterification on antioxidation.     

A semi-empirical conformational analysis of the interaction of n-alkanols with dipalmitoylphosphatidylcholine

The interaction of a homologous series of saturated aliphatic n-alkanols (containing 1–13 carbon atoms) with dipalmitoylphosphatidylcholine was studied by a semi-empirical conformational analysis. The minimal conformational energy of the isolated molecule at the hydrocarbon-water interface was calculated as the sum of the contributions resulting from the Van der Waals, torsional, electrostatic and transfer energies. From the conformers of minimal energies were calculated the hydrophilic-hydrophobic balance, the distance between hydrophilic and hydrophobic centres and the energies of interaction between homologous alkanols and between alkanols and lipids. Using these parameters, different modes of conformation, orientation and interaction of n-alkanols and dipalmitoylphosphatidylcholine were described. For methanol, ethanol and n-propanol, the gauche conformers were the most probable interacting only with the lipid polar heads. Only all-trans conformers were obtained for alkanols with longer acyl chains. n-Butanol to n-octanol form clusters in the lipid matrix. Longer n-alkanols are randomly distributed in the lipid layer. However, due to the increase in hydrocarbon chain-length, n-nonanol and higher alkanols have an interaction energy identical or superior to that found in a pure lipid monolayer, leading to a more ordered alkanol-lipid organization.     

Adsorption and decomposition mechanism of hexogen (RDX) on Al(111) surface by periodic DFT calculations

The adsorption of hexogen (RDX) molecule on the Al(111) surface was investigated by the generalized gradient approximation (GGA) of density functional theory (DFT). The calculations employ a supercell (4×4×3) slab model and three-dimensional periodic boundary conditions. The strong attractive forces between RDX molecule and aluminum atoms induce the N?O and N?N bond breaking of the RDX. Subsequently, the dissociated oxygen atoms, NO₂ group and radical fragment of RDX oxidize the Al surface. The largest adsorption energy is ?835.7 kJ mol^–1. We also investigated the adsorption and decomposition mechanism of RDX molecule on the Al(111) surface. The activation energy for the dissociation steps of V4 configuration is as large as 353.1 kJ mol^–1, while activation energies of other configurations are much smaller, in the range of 70.5–202.9 kJ mol^–1. The N?O is even easier than the N?NO₂ bond to decompose on the Al(111) surface.     

Theoretical study of the geometrical and electronic structures and thermochemistry of spherophanes

A set of supramolecular cage-structures—spherophanes—was studied at the density functional B3LYP level. Full geometrical structure optimisations were made with 6–31G and 6–31G(d) basis sets followed by frequency calculations, and electronic energies were evaluated at B3LYP/6–31++G(d,p). Three different symmetries were considered: C1, Ci, and Oh. It was found that the bonds between the benzene rings are very long to allow π-electron delocalisation between them. These spherophanes show portal openings of 2.596 Å in Spher1, 4.000 Å in Meth2, 3.659 Å in Oxa3, and 4.412 Å in Thia4. From the point of view of potential host–guest interaction studies, it should also be noted that the atoms nearest to the centre of the cavities are carbons bonded to X groups. These supramolecules seem to exhibit relatively large gap HOMO?LUMO: 2.89 eV(Spher1), 5.26 eV(Meth2), 5.73 eV(Oxa3), and 4.82 eV(Thia4). The calculated ΔH°_f (298.15 K) values at B3LYP/6–31G(d) are (in kcal mol^?1) 750.98, 229.78, ?10.97, and 482.49 for Spher1, Meth2, Oxa3, and Thia4, respectively. Using homodesmotic reactions, relative to Spher1, the spherophanes Meth2, Oxa3, and Thia4 are less strained by ?399.13 kcal mol^?1, ?390.40 kcal mol^?1, and ?411.38 kcal mol^?1, respectively. Their infrared and ¹³C NMR calculated spectra are reported.     

A mechanistic investigation on the formation and rearrangement of silaspiropentane: A theoretical study

The formation of silaspiropentane from addition of singlet silacyclopropylidene 1 and silacyclopropylidenoid 8 to ethylene has been investigated separately at the B3LYP, X3LYP, WB97XD, and M05–2X theories using the 6–31+G(d,p) basis set. The silacycloproylidenoid addition follows a stepwise route. In contrast, a concerted mechanism occurs for silacyclopropylidene addition. Moreover, the intramolecular rearrangements of silaspiropentane 9 to methylenesilacyclobutane 11 and 2-silaallene?+?ethylene 12 have been studied extensively. The required energy barrier for the isomerization of 9 to 10 was determined to be 44.0 kcal mol^?1 at the B3LYP/6–31+G(d,p) level. After formation of 10, the rearrangement to methylenesilacyclobutane 12 is highly exergonic by ?15.9 kcal mol^?1, which makes this reaction promising. However, the conversion of 9 to 11 is calculated to be quite endergonic, by 26.5 kcal mol^?1.     

Structural and functional studies on a variant of cystatin purified from brain of Capra hircus

Cystatins, known for their ubiquitous presence in mammalian system are thiol protease inhibitors serving important physiological functions. Here, we present a variant of cystatin isolated from brain of Capra hircus (goat) which is glycosylated but lacks disulphide bonds. Caprine brain cystatin (CBC) was isolated using alkaline treatment, ammonium sulphate fractionation (40–60%) and gel filtration chromatography on Sephacryl S-100HR column with an overall yield of 26.29% and 322-fold purification. The inhibitor gave a molecular mass of ~44 kDa as determined by SDS-PAGE and gel filtration behaviour. The Stokes radius and diffusion coefficient of CBC were 27.14 Å and 8.18 × 10^?7 cm² s^?1, respectively. Kinetic data revealed that CBC inhibited thiol proteases reversibly and competitively, with the highest inhibition towards papain (K_i = 4.10 nM) followed by ficin and bromelain. CBC possessed 34.7% α-helical content as observed by CD spectroscopy. UV, fluorescence, CD and FTIR spectroscopy revealed significant conformational change upon CBC-papain complex formation. Isothermal titration calorimetry (ITC) was used to measure the thermodynamic parameters – ΔH, ΔS, ΔG along with N (binding stoichiometry) for CBC-papain complex formation. Binding stoichiometry (N = .97 ± .07 sites) for the CBC-papain complex indicates that cystatin is surrounded by nearly one papain molecule. Negative ΔH (?5.78 kcal mol^?1) and positive ΔS (11.01 cal mol^?1 deg^?1) values suggest that the interaction between CBC and papain is enthalpically as well as entropically favoured process. The overall negative ΔG (?9.19 kcal mol^?1) value implies a spontaneous CBC-papain interaction.     

Oxidation reactivity of zinc–cysteine clusters in metallothionein

Evaluating the reactivity of the metal–thiolate clusters in metallothionein (MT) is a key step in understanding the biological functions of this protein. The effects of the metal clustering and protein environment on the thiolate reactivity with hydrogen peroxide (H₂O₂) were investigated by performing quantum theory calculations with chemical accuracy at two levels of complexity. At the first level, the reactivity with H₂O₂ of a model system ([(Zn)₃(MeS)₉]^3?, MeS is methanethiolate) of the β domain cluster of MT was evaluated using density functional theory (DFT) with the mPW1PW91 functional. At the second level of complexity, the protein environment was included in the reactant system and the calculations were performed with the hybrid ONIOM method combining the DFT–mPW1PW91 and the semiempirical PM6 levels of theory. In these conditions, the energy barrier for the oxidation of the most reactive terminal thiolate was 21.5 kcal mol^?1. This is 3 kcal mol^?1 higher than that calculated for the terminal thiolate in the model system [(Zn)₃(MeS)₉]^3? and about 7 kcal mol^?1 higher than that obtained for the free thiolate. In spite of this rise of the energy barrier induced by the protein environment, the thiolate oxidation by H₂O₂ is confirmed as a possible way for metal release from MT. On the other hand, the results suggest that the antioxidant role of MT in the living cell cannot be as important as that of glutathione (which bears a free thiol).     

How do halogen bonds (S–O⋯I,N–O⋯I and C–O⋯I) and halogen–halogen contacts (C–I⋯I–C,C–F⋯F–C) subsist in crystal structures? A quantum chemical insight

Thirteen X-ray crystal structures containing various non-covalent interactions such as halogen bonds, halogen–halogen contacts and hydrogen bonds (I?N, I?F, I?I, F?F, I?H and F?H) were considered and investigated using the DFT-D3 method (B97D/def2-QZVP). The interaction energies were calculated at MO62X/def2-QZVP and MP2/aug-cc-pvDZ level of theories. The higher interaction and dispersion energies (2nd crystal) of ?9.58 kcal mol^?1 and ?7.10 kcal mol^?1 observed for 1,4-di-iodotetrafluorobenzene bis [bis (2-phenylethyl) sulfoxide] structure indicates the most stable geometrical arrangement in the crystal packing. The electrostatic potential values calculated for all crystal structures have a positive σ-hole, which aids understanding of the nature of σ-hole bonds. The significance of the existence of halogen bonds in crystal packing environments was authenticated by replacing iodine atoms by bromine and chlorine atoms. Nucleus independent chemical shift analysis reported on the resonance contribution to the interaction energies of halogen bonds and halogen–halogen contacts. Hirshfeld surface analysis and topological analysis (atoms in molecules) were carried out to analyze the occurrence and strength of all non-covalent interactions. These analyses revealed that halogen bond interactions were more dominant than hydrogen bonding interactions in these crystal structures.

Open image in new window

Graphical Abstract Molecluar structure of 1,4-Di-iodotetrafluorobenzene bis(thianthrene 5-oxide) moelcule and its corresponding molecular electrostatic potential map for the view of σ-hole.

     

A DFT study of adsorption and decomposition of hexahydro-1,3,5-trinitro-1,3,5-triazine on Mg(0001) surface

The adsorption and decomposition of hexogen (RDX) molecule on the Mg(0001) surface were investigated by the generalized gradient approximation (GGA) of density functional theory (DFT). The calculations employed a supercell (4?×?4?×?4) slab model and three-dimensional periodic boundary conditions. The strong attractive forces between RDX molecule and magnesium atoms induce the RDX’s N???O bond breaking. Subsequently, the dissociated oxygen atoms and radical fragment of RDX oxidize the Mg surface. The largest adsorption energy is ?2104.0 kJ mol^-1. We also investigated the decomposition mechanism of RDX molecule on the Mg(0001) surface. The activation energy for the dissociation step of configuration V4 is as small as 2.5 kJ mol^-1, while activation energies of other configurations are much larger, in the range of 964.9–1375.1 kJ mol^-1. Mg powder is more active than Al powder, and Mg powder performs better in increasing the combustion exothermicity of RDX as well.