Theoretical investigation on the structure and thermodynamic properties of the 2,4-dinitroimidazole complex with methanol

The structure and thermodynamic properties of the 2, 4-dinitroimidazole complex with methanol were investigated using the B3LYP and MP2(full) methods with the 6-31++G(2d,p) and 6-311++G(3df,2p) basis sets. Four types of hydrogen bonds [N–H?O, C–H?O, O–H?O (nitro oxygen) and O–H?π] were found. The hydrogen-bonded complex having the highest binding energy had a N–H?O hydrogen bond. Analyses of natural bond orbital (NBO) and atoms-in-molecules (AIM) revealed the nature of the intermolecular hydrogen-binding interaction. The changes in thermodynamic properties from monomers to complexes with temperatures ranging from 200.0 to 800.0 K were investigated using the statistical thermodynamic method. Hydrogen-bonded complexes of 2,4-dinitroimidazole with methanol are fostered by low temperatures.

Synthetic and quantum chemical study on the regioselective addition of amines to methyl maleamate

Synthetic and theoretical studies were performed to gain insight into the regioselectivity in the mechanism of aspartyl-isoaspartyl formation, modeled by additions of ammonia and primary amines to methyl maleamate. Reactions between maleamate and aliphatic, araliphatic amines or O-methyl acetimidate lead to the formation of N-substituted isoasparaginates. The size of the amine and the activating effect of the amide and ester group on the double bond are the determining factors of the site of addition. The formation of both isomers was observed only in the case of ammonia addition. The regioselectivity was predicted on the basis of the charge distribution for low-energy methyl maleamate conformers, calculated at the B3LYP/6-311++G(2df,2pd)//B3LYP/6-31+G(d) level, both in gas phase and in methanol. The methyl isoasparaginate over methyl asparaginate product ratio was computed based on the free energy Boltzmann distribution of their conformers. The calculated 2 : 1 ratio is in agreement with the experimental regioselectivity of the addition of nitrogen nucleophiles.

Structural characterization of the (MeSH)4 potential energy surface

A random walk on the PES for (MeSH)₄ clusters produced 50 structural isomers held together by hydrogen-bonding networks according to calculations performed at the B3LYP/6–311++G** and MP2/6–311++G** levels. The geometric motifs observed are somewhat similar to those encountered for the methanol tetramer, but the interactions responsible for cluster stabilization are quite different in origin. Cluster stabilization is not related to the number of hydrogen bonds. Two distinct, well-defined types of hydrogen bonds scattered over a wide range of distances are predicted.

A B3LYP and MP2(full) theoretical investigation into the strength of the C-NO2 bond upon the formation of the molecule-cation interaction between Na+ 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 molecule-cation interaction between Na⁺ 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 in comparison with that in the isolated nitrotriazole molecule upon the formation of molecule-cation interaction. The increment of the C–NO₂ bond dissociation energy (ΔBDE) correlated well with the molecule-cation interaction energy. Electron density shifts analysis 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.

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.

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.     

Density functional theory investigation of electrophilic addition reaction of bromine to tricyclo[4.2.2.22,5]dodeca-1,5-diene1

The geometry and the electronic structure of tricyclo[4.2.2.2^2,5]dodeca-1,5-diene (TCDD) molecule were investigated by DFT/B3LYP and /B3PW91 methods using the 6-311G(d,p) and 6-311++G(d,p) basis sets. The double bonds of TCDD molecule are syn-pyramidalized. The structure of π-orbitals and their mutual interactions for TCDD molecule were investigated. Potential energy surface (PES) of the TCDD-Br₂ system was studied by B3LYP/6-311++G(d,p) method and the configurations [molecular charge-transfer (CT) complex, transition states (TS1 and TS2), intermediate (INT) and product (P)] corresponding to the stationary points (minima or saddle points) were determined. Initially, a molecular CT-complex forms between Br₂ and TCDD. With a barrier of 22.336 kcal mol^-1 the CT-complex can be activated to an intermediate (INT) with energy 15.154 kcal mol^-1 higher than that of the CT-complex. The intermediate (INT) then transforms easily (barrier 5.442 kcal mol^-1) into the final, N-type product. The total bromination is slightly exothermic. Accompanying the breaking of Br-Br bond, C1-Br, C5-Br and C2-C6 bonds are formed, and C1 = C2 and C5 = C6 double bonds transform into single bonds. The direction of the reaction is determined by the direction of intramolecular skeletal rearrangement that is realized by the formation of C2-C6 bond.

Relationship between structure and entropy contributions in an anthraquinone mercapto derivative

The structural and thermodynamic properties of an anthraquinone derivative were studied by means of quantum-chemical calculations. Conformational analysis using ab initio and density functional theory methods revealed 14 low-energy conformers. In order to discuss similarities and differences in entropy of the conformers, the rotational and vibrational contributions to entropy were correlated with changes in conformer structure. The component of the moment of inertia perpendicular to the molecular plane gives significant input to ΔS _rot , whereas the largest contributions to the ΔS _vib have vibrations associated with the τ _S1C20 coordinate.

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.

Predicting the preferred conformations of luteolin-4′-O-β-D-glucoside in gas phase: a comparison of two computational approaches

A tree-step computational approach has been applied to determine the lowest-energy conformers of luteolin-4′-O-β-D-glucoside (L4′G). Fifty-seven starting structures of the L4′G have been built, and then by performing with density functional theory (DFT) optimizations and second-order Møller-Plesset (MP2) calculations, the preferred conformations of L4′G are predicted. In order to test the accuracy of the computational approach, a hybrid Monte-Carlo multiple minimum (MCMM)/quantum mechanical (QM) approach is applied to determine the favorable conformers of L4′G. The alternative classification is employed to put similar conformations into the same catalogue according to the dihedral angles among the luteolin rings, glycosidic dihedral angles, and the orientations of hydroxyl and hydroxymethyl groups. The low-energy conformations are located after the optimizations at the HF/6-31G(d) and B3LYP/6-311+G(d) levels. Compared with the hybrid MCMM/QM approach, the tree-step computational approach not only remains accurate but also saves a lot of computing resources.

New insights into the stability of alkenes and alkynes, fluoro-substituted or not: a DFT, G4, QTAIM and GVB study

Many undergraduate organic chemistry books do not agree with the order of relative stability of alkenes towards hydrogenation reactions. Although they ascribe the stability of alkenes to the number and spatial position of the alkyl groups attached to the vinyl carbon atoms, results from the quantum theory of atoms in molecules indicate that the influence of an alkyl substituent on the stability of unsaturated hydrocarbons arises from the slight removal of electron density of the π bond, not from donation of their charge density to unsaturated carbon atoms as stated in many text books. There is an inverse relation between delocalization index—the number of shared electrons between two atoms, or Wiberg bond index of C=C bond—and the number of methyl groups attached to the vinyl carbon atoms. Electron withdrawing groups (EWGs) attached to unsaturated carbon atoms of alkenes and alkynes have two different behaviors: slight EWGs (alkyl groups) stabilize unsaturated carbon atoms while the strong EWG destabilizes the unsaturated carbon atoms. Generalized valence bond theory was also used to study the ambiguous behavior of fluorine substituents bonded to vinyl carbon atoms.

DFT studies of the adsorption and dissociation of H2O on the Al13 cluster: origins of this reactivity and the mechanism for H2 release

A theoretical study of the chemisorption and dissociation pathways of water on the Al₁₃ cluster was performed using the hybrid density functional B3LYP method with the 6-311+G(d, p) basis set. The activation energies, reaction enthalpies, and Gibbs free energy of activation for the reaction were determined. Calculations revealed that the H₂O molecule is easily adsorbed onto the Al₁₃ surface, forming adlayers. The dissociation of the first H₂O molecule from the bimolecular H₂O structure via the Grotthuss mechanism is the most kinetically favorable among the five potential pathways for O–H bond breaking. The elimination of H₂ in the reaction of an H₂O molecule with a hydrogen atom on the Al cluster via the Eley–Rideal mechanism has a lower activation barrier than the elimination of H₂ in the reaction of two adsorbed H atoms or the reaction of OH and H. Following the adsorption and dissociation of H₂O, the structure of Al₁₃ is distorted to varying degrees.

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.

Penta- and heteropentadienyl ligands coordinated to beryllium

In this work we have performed a systematic study of new organometallic complexes containing penta- and heteropentadienyl (CH₂CHCHCHX, X?=?CH₂, O, NH, S) ligands coordinated to beryllium. Calculated complexes were studied using the density functional theory (PBE) in combination with the 6-311++G(3d,2p) basis set. The coordination number on the beryllium atom varies according to the type of ligand. Pentadienyl ligand shows hapticities η¹ and η⁵, while heteropentadienyl ligands display η¹ and η² hapticities. A Wiberg bond indices study was performed in order to get information about their bond orders.

Exploration of various electronic properties along the reaction coordinate for hydration of Pt(II) and Ru(II) complexes; the CCSD, MPx, and DFT computational study

In the study behavior of molecular electrostatic potential, averaged local ionization energy, and reaction electronic flux along the reaction coordinate of hydration process of three representative Ru(II) and Pt(II) complexes were explored using both post-HF and DFT quantum chemical approximations. Previously determined reaction mechanisms were explored by more detailed insight into changes of electronic properties using ωB97XD functional and MP2 method with 6–311++G(2df,2pd) basis set and CCSD/6–31(+)G(d,p) approach. The dependences of all examined properties on reaction coordinate give more detailed understanding of the hydration process.

Bis-dibenzo[a.i]fluorenylidene, does it exist as stable 1,2-diradical?

The geometries, energies, and electronic properties of the two possible configurations of bis-[dibenzo[a.i]fluorenylidene] were investigated theoretically by density functional theory DFT B3LYP at the UB3LYP/6-311?+?G(2d,p) // UB3LYP/6-31?+?G(d,p) level of theory. According to the performed calculations, it was found that the singlet is 3.4?kcal?mol^-1 lower in energy compared to triplet state at room temperature. This gap is compared with those of other alkenes like ethylene, (61.9?kcal?mol^-1) tetra-tert-butyethylene, (6.4?kcal?mol^-1) and bis-fluorenylidene (19.5?kcal?mol^-1). These results confirm the experimental findings of the paramagnetic properties determined by Franzen and Joschek. The low singlet-triplet gap in the case of bis-[dibenzo[a.i]fluorenylidene] is the result of a steric destabilization of the singlet due to strain and stabilization of the triplet electronic state by delocalization of each free electron within each aromatic moiety. This correlates with the special electronic structure of the triplet state of this compound, where facial interaction of two hydrogen atoms lying close to the lobes of each p-orbital occupied with a single electron at the distorted double bond in the triplet electronic state.

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.

A B3LYP and MP2(full) theoretical investigation into the cooperativity effect between dihydrogen-bonding and H–M∙∙∙π (M = Li, Na, K) interactions among HF, MH with the π-electron donor C2H2, C2H4 or C6H6

The DFT-B3LYP/6-311++G(3df,2p) and MP2(full)/6-311++G(3df,2p) calculations were carried out on the binary complex formed by HM (M?=?Li, Na, K) and HF or the π-electron donor (C₂H₂, C₂H₄, C₆H₆), as well as the ternary system FH???HM???C₂H₂/C₂H₄/C₆H₆. The cooperativity effect between the dihydrogen-bonding and H–M???π interactions was investigated. The result shows that the equilibrium distances R _H???H and R _M???π in the ternary complex decrease and both the H???H and H–M???π interactions are strengthened when compared to the corresponding binary complex. The cooperativity effect of the dihydrogen bond on the H–M???π interaction is more pronounced than that of the M???π bond on the H???H interaction. Furthermore, the values of cooperativity effect follow the order of FH???HNa???π?>?FH???HLi???π?>?FH???HK???π and FH???HM???C₆H₆?>?FH???HM???C₂H₄?>?FH???HM???C₂H₂. The nature of the cooperativity effect was revealed by the analyses of the charge of the hydrogen atoms in H???H moiety, atom in molecule (AIM) and electron density shifts methods.

Spectroscopic investigations and hydrogen bond interactions of 8-aza analogues of xanthine, theophylline and caffeine: a theoretical study

The structure, spectral properties and the hydrogen bond interactions of 8-aza analogues of xanthine, theophylline and caffeine have been studied by using quantum chemical methods. The time-dependent density functional theory (TD-DFT) and the singly excited configuration interaction (CIS) methods are employed to optimize the excited state geometries of isolated 8-azaxanthine, 8-azatheophylline tautomers and 8-azacaffeine in both the gas and solvent phases. The solvent phase calculations are performed using the polarizable continuum model (PCM). The absorption and emission spectra are calculated using the time-dependent density functional theory (TD-DFT) method. The results from the TD-DFT calculations reveal that the excitation spectra are red shifted relative to absorption in aqueous medium. These changes in the transition energies are qualitatively comparable to the experimental data. The examination of molecular orbital reveals that the molecules with a small H→L energy gap possess maximum absorption and emission wavelength. The relative stability and hydrogen bonded interactions of mono and heptahydrated 8-azaxanthine, 8-azatheophylline tautomers and 8-azacaffeine have been studied using the density functional theory (DFT) and Møller Plesset perturbation theory (MP2) implementing the 6-311++G(d,p) basis set. The formation of strong N-H…O bond has resulted in the highest interaction energy among the monohydrates. Hydration does not show any significant impact on the stability of heptahydrated complexes. The atoms in molecule (AIM) and natural bonding orbital (NBO) analyses have been performed to elucidate the nature of the hydrogen bond interactions in these complexes.