Entropy versus aromaticity in the conformational dynamics of aromatic rings

Comparison of the results of Car-Parrinello molecular dynamics simulations of isolated benzene, pyrimidine and 1,2,4-triazine molecules reveals that the unusually low population of planar geometry of the benzene ring is caused by entropy effects despite its high aromaticity. The decrease in symmetry of the molecule results in smaller changes in entropy and Gibbs free energy due to out-of-plane deformations of the ring, leading to an increase in the population of planar geometry of the ring. This leads to differences in the topology of potential energy and Gibbs free energy surfaces.

A B3LYP and MP2(full) theoretical investigation into the strength of the C–NO2 bond upon the formation of the intermolecular hydrogen-bonding interaction between HF and the nitro group of nitrotriazole or its methyl derivatives

The changes of bond dissociation energy (BDE) in the C–NO₂ bond and nitro group charge upon the formation of the intermolecular hydrogen-bonding interaction between HF and the nitro group of 14 kinds of nitrotriazoles or methyl derivatives were investigated using the B3LYP and MP2(full) methods with the 6-311++G**, 6-311++G(2df,2p) and aug-cc-pVTZ basis sets. The strength of the C–NO₂ bond was enhanced and the charge of nitro group turned more negative in complex in comparison with those in isolated nitrotriazole molecule. The increment of the C–NO₂ bond dissociation energies correlated well with the intermolecular H-bonding interaction energies. Electron density shifts analyses showed that the electron density shifted toward the C–NO₂ bond upon complex formation, leading to the strengthened C–NO₂ bond and the possibly reduced explosive sensitivity.

Reactions of ketones with aromatics in acid media. The effect of trifluoromethyl groups and the acidity media. A theoretical study

The reactions of acetone, 2,2,2-trifluoroacetone and hexafluoroacetone in methanesulfonic (MSA) and triflic acids (TFSA) with benzene have been studied at M06-2X/6-311+G(d,p) level using cluster-continuum model, where the carbonyl group is explicitly solvated by acid molecules. The introduction of a trifluoromethyl group into the ketone structure reduces the activation energy of the tetrahedral intermediates formation due to an increase of the electrophilicity of the carbonyl group and raises the activation and the reaction energies of the C-O bond cleavage in formed carbinol due to the destabilization of the corresponding carbocation. The introduction of the second trifluoromethyl group inhibits the hydroxyalkylation reaction due to a very strong increase of the reaction and activation energies of the C-O bond cleavage which becomes the rate determining step. The most important catalytic effect of TFSA compared to MSA is not the protonation of the ketone carbonyl, but the reduction of the activation and reaction energies of the carbinol C-O bond cleavage due to better protosolvation properties. Even for TFSA no complete proton transfer to carbonyl oxygen has been observed for free ketones. Therefore, the protonation energies of free ketones cannot be considered as a measure of ketone reactivity in the hydroxyalkylation reaction.

Nucleus-independent chemical shift criterion for aromaticity in π-extended tetraoxa[8]circulenes

Recently synthesized π-extended symmetrical tetraoxa[8]circulenes that exhibit electroluminescent properties were calculated at the density functional theory (DFT) level using the quantum theory of atoms in molecules (QTAIM) approach to electron density distribution analysis. Nucleus-independent chemical shift (NICS) indices were used to characterize the aromaticity of the studied molecules. The tetraoxa[8]circulene molecules were found to consist of two antiaromatic perimeters (according to the Hückel “4n” antiaromaticity rule) that include 8 and 24 π-electrons. Conversely, NICS calculations demonstrated the existence of a common π-extended system (distributed like a flat ribbon) in the studied tetraoxa[8]circulene molecules. Thus, these symmetrical tetraoxa[8]circulene molecules provide examples of diatropic systems characterized by the presence of induced diatropic ring currents.

Deamination features of 5-hydroxymethylcytosine,a radical and enzymatic DNA oxidation product

Geometric consequences of electron delocalization were studied for all possible adenine tautomers in aqueous solution by means of ab initio methods {PCM(water)//DFT(B3LYP)/6-311+G(d,p)} and compared to those in the gas phase {DFT(B3LYP)/6-311+G(d,p)}. To measure the consequences of any type of resonance conjugation (π-π, n-π, and σ-π), the geometry-based harmonic oscillator model of electron delocalization (HOMED) index, recently extended to the isolated (DFT) and hydrated (PCM//DFT) molecules, was applied to the molecular fragments (imidazole, pyrimidine, 4-aminopyrimidine, and purine) and also to the whole tautomeric system. For individual tautomers, the resonance conjugations and consequently the bond lengths strongly depend on the position of the labile protons. The HOMED indices are larger for tautomers (or their fragments) possessing the labile proton(s) at the N rather than C atom. Solvent interactions with adenine tautomers slightly increase the resonance conjugations. Consequently, they slightly shorten the single bonds and lengthen the double bonds. When going from the gas phase to water solution, the HOMED indices increase (by less than 0.15 units). There is a good relation between the HOMED indices estimated in water solution and those in the gas phase for the neutral and ionized forms of adenine. Subtle effects, being a consequence of intramolecular interactions between the neighboring groups, are so strongly reduced by solvent that the relation between the HOMED indices and the relative energies for the neutral adenine tautomers seems to be better in water solution than in the gas phase.

Correlation of macrocycle distortion with oxidation potentials of iron(III) porphyrins: molecule structure of the sterically crowded chloro-iron(III) 7,8,17,18-tetrabromo-5,10,15,20-tetraphenylporphyrin

The chloro-iron(III) complex of the tetraphenylporphyrin tetrabrominated at the antipodal β-pyrrole positions [(7,8,17,18-tetrabromo-5,10,15,20-tetraphenyl)porphyrin] has been synthesized and characterized by spectroscopy and X-ray crystallography. The iron atom is bonded to the chloride ion and the four pyrrole nitrogens. The Fe-Cl bond distance is 2.209(4) Å, and the mean value of the two opposite Fe-N_p lengths at the brominated pyrrole rings is 2.079(8) Å, whereas the mean value of the two opposite Fe-N_p′ bond distances at the non-brominated pyrrole rings is 2.041(8) Å. The X-ray structure determination and the analysis of the UV-Vis spectra obtained in solution and on thin films indicate that |FeCl(tpp-Br₄)| (1) is principally saddle-shaped in the solid state and in solution. Variable-temperature (195–325 K) ¹H NMR spectroscopy confirms the high-spin state (S=5/2) of the iron(III) center and indicates that the saddle-shaped conformation of 1 is maintained in solution. EPR spectra obtained in frozen CH₂Cl₂ solution and in the solid state show a rhombic symmetry with g values of 6.25, 5.70 and 1.99. Kadish et al. have shown that the one-electron oxidation potential of 1 increases only by 0.06?V relative to that of the non-brominated complex |FeCl(tpp)|. The present study indicates that the increase of the first oxidation potential of 1 is related to the non-planar distortion of the porphyrin. Relative to the unbrominated derivative |FeCl(tpp)|, this distortion destabilizes the π system of the macrocycle and thus compensates for the effects of the four electron-withdrawing bromine substituents.     

Quantitative analysis of molecular surfaces: areas, volumes, electrostatic potentials and average local ionization energies

We describe a procedure for performing quantitative analyses of fields f(r) on molecular surfaces, including statistical quantities and locating and evaluating their local extrema. Our approach avoids the need for explicit mathematical representation of the surface and can be implemented easily in existing graphical software, as it is based on the very popular representation of a surface as collection of polygons. We discuss applications involving the volumes, surface areas and molecular surface electrostatic potentials, and local ionization energies of a group of 11 molecules.

A new interaction mechanism of LiNH2 with MgH2: magnesium bond

Quantum chemical calculations were performed for LiNH₂–HMgX (X?=?H, F, Cl, Br, CH₃, OH, and NH₂) complexes to propose a new interaction mechanism between them. This theoretical survey showed that the complexes are stabilized through the combinative interaction of magnesium and lithium bonds. The binding energies are in the range of 63.2–66.5 kcal mol^?1, i.e., much larger than that of the lithium bond. Upon complexation, both Mg–H and Li–N bonds are lengthened. Substituents increase Mg-H bond elongation and at the same time decrease Li-N bond elongation. These cyclic complexes were characterized with the presence of a ring critical point and natural population analysis charges.

A DFT study on the initial stage of thermal degradation of Poly(methyl methacrylate)/carbon nanotube system

DFT calculations, with VWN exchange correlation functional and double numeric basis set, were used to evaluate the energies required for the scission reactions taking place in the initial stage of the thermal degradation of Poly(methyl methacrylate) (PMMA) in the presence of a carbon nanotube (CNT). Side group and main chain scissions were investigated. The results averaged from five configurations of pure PMMA (DP?=?5) were used as references and compared to the results obtained for the five same configurations of PMMA grafted on three carbon nanotubes of similar diameter (1.49 nm). The bond dissociation energies (BDE) of main chain scission evaluated for grafted PMMA was 4 % less endothermic than for pure PMMA. These results seemed independent of the tested chirality (11,11); (12,10) and (16,5) of the carbon nanotubes. Comparisons with the BDE of the weakest bonds due to intrinsic defaults (head to head and unsaturated end chain) were performed.

Stereoselectivity of chalcone isomerase with chalcone derivatives: a computational study

Chalcone isomerase (CHI) catalyzes the intramolecular cyclization of chalcones into flavonoids. The activity of CHI is essential for the biosynthesis of flavonoids precursors of floral pigments and phenylpropanoid plant defense compounds. In the present study, we explored the detailed binding structures and binding free energies for two different active site conformations of CHI with s-cis/s-trans conformers of three chalcone compounds by performing molecular dynamics (MD) simulations and binding free energy calculations. The computational results indicate that s-cis/s-trans conformers of chalcone compounds are orientated in the similar binding position in the active site of CHI and stabilized by the different first hydrogen bond network and the same second hydrogen bond network. The first hydrogen bond network results in much lower binding affinity of s-trans conformer of chalcone compound with CHI than that of s-cis conformer. The conformational change of the active site residue T48 from indirectly interacting with the substrate via the second hydrogen bond network to directly forming the hydrogen bond with the substrates cannot affect the binding mode of both conformers of chalcone compounds, but remarkably improves the binding affinity. These results show that CHI has a strong stereoselectivity. The calculated binding free energies for three chalcone compounds with CHI are consistent with the experimental activity data. In addition, several valuable insights are suggested for future rational design and discovery of high-efficiency mutants of CHI.

Modulating weak intramolecular interactions through the formation of beryllium bonds: complexes between squaric acid and BeH2

The electronic structure of the two most stable isomers of squaric acid and their complexes with BeH₂ were investigated at the B3LYP/6-311?+?G(3df,2p)// B3LYP/6-31?+?G(d,p) level of theory. Squaric acid forms rather strong beryllium bonds with BeH₂, with binding energies of the order of 60 kJ?mol^?1. The preferential sites for BeH₂ attachment are the carbonyl oxygen atoms, but the global minima of the potential energy surfaces of both EZ and ZZ isomers are extra-stabilized through the formation of a BeH···HO dihydrogen bond. More importantly, analysis of the electron density of these complexes shows the existence of significant cooperative effects between the beryllium bond and the dihydrogen bond, with both becoming significantly reinforced. The charge transfer involved in the formation of the beryllium bond induces a significant electron density redistribution within the squaric acid subunit, affecting not only the carbonyl group interacting with the BeH₂ moiety but significantly increasing the electron delocalization within the four membered ring. Accordingly the intrinsic properties of squaric acid become perturbed, as reflected in its ability to self-associate.

Theoretical study on cooperative effects between X⋯N and X⋯Carbene halogen bonds (X = F,Cl,Br and I)

Quantum chemical calculations are performed to study the interplay between halogen?nitrogen and halogen?carbene interactions in NCX?NCX?CH₂ complexes, where X?=?F, Cl, Br and I. Molecular geometries and interaction energies of dyads and triads are investigated at the MP2/aug-cc-pVTZ level of theory. It is found that the X?N and X?C_carbene interaction energies in the triads are larger than those in the dyads, indicating that both the halogen bonding interactions are enhanced. The estimated values of cooperative energy E _coop are all negative with much larger E _coop in absolute value for the systems including iodine. The nature of halogen bond interactions of the complexes is analyzed using parameters derived from the quantum theory atoms in molecules methodology and energy decomposition analysis.

A molecular dynamics study on sI hydrogen hydrate

A molecular dynamics simulation is carried out to explore the possibility of using sI clathrate hydrate as hydrogen storage material. Metastable hydrogen hydrate structures are generated using the LAMMPS software. Different binding energies and radial distribution functions provide important insights into the behavior of the various types of hydrogen and oxygen atoms present in the system. Clathrate hydrate cages become more stable in the presence of guest molecules like hydrogen.

Competition between halogen, dihalogen and hydrogen bonds in bromo- and iodomethanol dimers

O-H…X and O-H…O H-bonds as well as C-X…X dihalogen and C-X…O halogen bonds have been investigated in halomethanol dimers (bromomethanol dimer, iodomethanol dimer, difluorobromomethanol…bromomethanol complex and difluoroiodomethanol…iodomethanol complex). Structures of all complexes were optimized at the counterpoise-corrected MP2/cc-pVTZ level and single-point energies were calculated at the CCSD(T)/aug-cc-pVTZ level. Energy decomposition for the bromomethanol dimer complex was performed using the DFT-SAPT method based on the aug-cc-pVTZ basis set. OH…O and OH…X H-bonds are systematically the strongest in all complexes investigated, with the former being the strongest bond. Halogen and dihalogen bonds, being of comparable strength, are weaker than both H-bonds but are still significant. The strongest bonds were found in the difluoroiodomethanol…iodomethanol complex, where the O-H…O H-bond exceeds 7 kcal mol^-1, and the halogen and dihalogen bonds exceed 2.5 and 2.3 kcal mol^-1, respectively. Electrostatic energy is dominant for H-bonded structures, in halogen bonded structures electrostatic and dispersion energies are comparable, and, finally, for dihalogen structures the dispersion energy is clearly dominant.

Electronic structures of bisnoradamantenyl and bisnoradamantanyl dications and related species

The highly pyramidalized molecule bisnoradamantene is extremely reactive toward nucleophiles and dienes. In this work, we studied the electronic structure of bisnoaradamantene, as well as those of its cation and dication, which are previously unreported carbonium ions. According to QTAIM and MO analysis, there is a 3c-2e bonding system in the bisnoradamantenyl cation and a 4c-2e bonding system in the bisnoradamantenyl dication. A topological study indicated that, on going from bisnoradamantene to its dication, π-bond interaction with the bridgehead carbon atom increases. Additional study of the bisnoradamantanyl dication also indicated that it has two multicenter bonding systems. Comparison of the D3BIA and NICS aromaticity indices of these molecules and other derivatives indicates that these indices are well correlated, and analysis of these indices shows that the cationic and dicationic bisnoradamantenyl species are homoaromatic.

Dissociation quenching using exceptional points

We examine a short way to reach an exceptional point that corresponds to a coalescence of two resonance energies. The application concerns the photodissociation of the Na₂ molecule exposed to a laser field. In this case, the resonances can be correlated with the field-free vibrational states of the diatomic species. The resonances are due to the field-induced coupling with the continuum of a repulsive potential. We also draw attention to a new kind of exceptional point involving a resonance originating from a vibrational state coalescing with a shape-type resonance of the repulsive potential. A laser control scenario, aiming at the adiabatic transport from this field-free decaying state to a stable field-free vibrational state, is discussed in terms of field-induced dissociation quenching.

Low-energy conformers of pamidronate and their intramolecular hydrogen bonds: a DFT and QTAIM study

Extensive DFT and ab initio calculations were performed to characterize the conformational space of pamidronate, a typical pharmaceutical for bone diseases. Mono-, di- and tri-protic states of molecule, relevant for physiological pH range, were investigated for both canonical and zwitterionic tautomers. Semiempirical PM6 method were used for prescreening of the single bond rotamers followed by geometry optimizations at the B3LYP/6-31++G(d,p) and B3LYP/6-311++G(d,p) levels. For numerous identified low energy conformers the final electronic energies were determined at the MP2/6-311++G(2df,2p) level and corrected for thermal effects at B3LYP level. Solvation effects were also considered via the COSMO and C-PCM implicit models. Reasonable agreement was found between bond lengths and angle values in comparison with X-ray crystal structures. Relative equilibrium populations of different conformers were determined from molecular partition functions and the role of electronic, vibrational and rotational degrees of freedom on the stability of conformers were analyzed. For no level of theory is a zwitterionic structure stable in the gas-phase while solvation makes them available depending on the protonation state. Geometrically identified intramolecular hydrogen bonds were analyzed by QTAIM approach. All conformers exhibit strong inter-phosphonate hydrogen bonds and in most of them the alkyl-amine side chain is folded on the P-C-P backbone for further hydrogen bond formation.

Non-covalent interactions – QTAIM and NBO analysis

MP2(full)/6-311++G(3df,3pd) calculations were carried out on complexes linked through various non-covalent Lewis acid – Lewis base interactions. These are: hydrogen bond, dihydrogen bond, hydride bond and halogen bond. The quantum theory of ´atoms in molecules´ (QTAIM) as well as the natural bond orbitals (NBO) method were applied to analyze properties of these interactions. It was found that for the A-H…B hydrogen bond as well as for the A-X…B halogen bond (X designates halogen) the complex formation leads to the increase of s-character in the A-atom hybrid orbital aimed toward the H or X atom. In opposite, for the A…H-B hydride bond, where the H-atom possesses negative charge, the decrease of s-character in the B-atom orbital is observed. All these changes connected with the redistribution of the electron charge being the effect of the complex formation are in line with Bent´s rule. The numerous correlations between energetic, geometrical, NBO and QTAIM parameters were also found.