Crystal and geometry‐optimized structure of an anthracene‐based Diels–Alder adduct

Computational calculations of an anthracene‐based Diels–Alder adduct, namely, 17‐ethyl‐1‐hydroxymethyl‐17‐azapentacyclo[6.6.5.0^2,7.0^9,14.0^15,19]nonadeca‐2,4,6,9,11,13‐hexaene‐16,18‐dione, C₂₁H₁₉NO₃, predicting density functional theory (DFT) optimized geometries in the gas phase are compared in terms of accuracy relative to the solid‐state crystal structure and computational cost. Crystal structure determination and Hirshfeld surface analysis of the racemic product reveal that the molecules are linked by O—H…O=C hydrogen bonds between the hydroxy and carbonyl groups, accounting for 9.5% of the intermolecular contacts, while H…H contacts represent 56.9% of the total. Boltzmann population analysis of computed relative rotamer abundances in the gas phase are based on lower‐level geometry optimization and thermochemical corrections coupled with higher‐level electronic energy calculations using the B2PLYP double‐hybrid functional. As expected, the choice of density functional has a greater effect than the basis set on accuracy for all levels of theory. For any given functional, increasing the basis set size did not always correlate with increasingly accurate structures. The hybrid functional B3LYP without dispersion correction routinely gave the closest approximations to the crystal structure where the B3LYP/aug‐cc‐pVDZ combination afforded the best structure (r.m.s. deviation = 0.1314 Å). However, the B3LYP/6‐31+G(d,p) level of theory represents the best compromise between accuracy (r.m.s. deviation = 0.1388 Å) and cost as it yielded appreciably accurate results in a fraction of the time compared to other method combinations.     

Glycine Dimers: Structure,Stability, and Medium Effects

We report a study on different ionization states and conformations of the bimolecular (Gly)₂ system by means of quantum mechanical calculations. Optimized geometries for energy minima of the glycine dimer, as well as relative energies and free energies were computed as a function of the medium: gas phase, nonpolar polarizable solvent, and aqueous solution. The polarizable continuum model was employed to account for solvation effects. Energy calculations were done using the MP2/aug‐cc‐pVTZ and B3LYP/6‐311+G(2df,2p) methods on B3LYP/6‐31+G(d,p) optimized structures (some single‐point energy calculations were also done using the B3PW91 and PBE1KCIS methods). Ionized forms of the glycine dimer (either zwitterion–zwitterion or neutral–zwitterion) are predicted to exist in all media, in contrast to amino acid monomers. In aqueous solution, dimerization is an exergonic process (?4 kcal mol^?1). Thus, according to our results, zwitterion–zwitterion Gly dimers might be abundant in supersaturated glycine aqueous solutions, a fact that has been connected with the structure of α‐glycine crystals but that remains controversial in the literature. Another noticeable result is that zwitterion–zwitterion interactions are substantially underestimated when computed using methods based on density functional theory. For comparison, some calculations for the dimer of the simplest chiral amino acid alanine were done as well and differences to the glycine dimer are discussed.     

Reactions of guanine with methyl chloride and methyl bromide: O6‐methylation versus charge transfer complex formation

Density functional theory (DFT) at the B3LYP/6‐31+G* and B3LYP/AUG‐cc‐pVDZ levels was employed to study O6‐methylation of guanine due to its reactions with methyl chloride and methyl bromide and to obtain explanation as to why the methyl halides cause genotoxicity and possess mutagenic and carcinogenic properties. Geometries of the various isolated species involved in the reactions, reactant complexes (RCs), and product complexes (PCs) were optimized in gas phase. Transition states connecting the reactant complexes with the product complexes were also optimized in gas phase at the same levels of theory. The reactant complexes, product complexes, and transition states were solvated in aqueous media using the polarizable continuum model (PCM) of the self‐consistent reaction field theory. Zero‐point energy (ZPE) correction to total energy and the corresponding thermal energy correction to enthalpy were made in each case. The reactant complexes of the keto form of guanine with methyl chloride and methyl bromide in water are appreciably more stable than the corresponding complexes involving the enol form of guanine. The nature of binding in the product complexes was found to be of the charge transfer type (O6mG⁺ · X^?, X?Cl, Br). Binding of HCl, HBr, and H₂O molecules to the PCs obtained with the keto form of guanine did not alter the positions of the halide anions in the PCs, and the charge transfer character of the PCs was also not modified due to this binding. Further, the complexes obtained due to the binding of HCl, HBr, and H₂O molecules to the PCs had greater stability than the isolated PCs. The reaction barriers involved in the formation of PCs were found to be quite high (～50 kcal/mol). Mechanisms of genotoxicity, mutagenesis and carcinogenesis caused by the methyl halides appear to involve charge transfer‐type complex formation. Thus the mechanisms of these processes involving the methyl halides appear to be quite different from those that involve the other strongly carcinogenic methylating agents. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2007     

Computational thermochemistry of glycolaldehyde

We use a variant of the focal point analysis to refine estimates of the relative energies of the four low‐energy torsional conformers of glycolaldehyde. The most stable form is the cis‐cis structure which enjoys a degree of H‐bonding from hydroxyl H to carbonyl O; here dihedral angles τ₁ (O?C? C? O) and τ₂ (C? C? O? H) both are zero. We optimized structures in both CCSD(T)/aug‐cc‐pVDZ and aug‐cc‐pVTZ; the structures agree within 0.01 Å for bond lengths and 1.0 degrees for valence angles, but the larger basis brings the rotational constants closer to experimental values. According to our extrapolation of CCSD(T) energies evaluated in basis sets ranging to aug‐cc‐pVQZ the trans‐trans form (180°, 180°) has a relative energy of 12.6 kJ/mol. The trans‐gauche conformer (160°, ±75°) is situated at 13.9 kJ/mol and the cis‐trans form (0°, 180°) at 18.9 kJ/mol. Values are corrected for zero point vibrational energy by MP2/aug‐cc‐pVTZ frequencies. Modeling the vibrational spectra is best accomplished by MP2/aug‐cc‐pVTZ with anharmonic corrections. We compute the Watsonian parameters that define the theoretical vibrational‐rotational spectra for the four stable conformers, to assist the search for these species in the interstellar medium. Six transition states are located by G4 and CBS‐QB3 methods as well as extrapolation using energies for structures optimized in CCSD(T)/aug‐cc‐pVDZ structures. We use two isodesmic reactions with two well‐established thermochemical computational schemes G4 and CBS‐QB3 to estimate energy enthalpy and Gibbs energy of formation as well as the entropy of the gas phase system. Our extrapolated electronic energies of species appearing in the isodesmic reactions produce independent values of thermodynamic quantities consistent with G4 and CBS‐QB3. © 2013 Wiley Periodicals, Inc.     

Uncatalyzed peptide bond formation between two double amino acid molecules in the gas phase

The gas phase mechanism for peptide bond formation between two double amino acid (DAA) molecules ((NH₂)₂C(COOH)₂) is investigated in the absence of any catalysts. Two different paths, concerted and stepwise, each leading to both cis and trans DAA‐DAA dipeptide products (four mechanisms total) are examined on the basis of theoretical calculations carried out at the CCSD(T)/aug‐cc‐pVDZ//MP2/aug‐cc‐pVDZ level. The investigation indicates that the concerted mechanism leading to the trans configuration of the peptide bond in the DAA‐DAA dipeptide product is thermodynamically favored by about 5 kcal mol^?1 and requires slightly less energy than the remaining pathways considered. Moreover, the peptide bond formation process between two DAA molecules in the gas phase resembles the analogous reactions between two natural amino acids.     

Reactivities of hydroxyl and perhydroxyl radicals toward cytosine and thymine: A comparative study

As the hydroxyl (OH) and perhydroxyl (OOH) radicals are known to play important roles in biological systems, their reactions with cytosine and thymine were studied. Addition reactions of these radicals at different sites of cytosine and thymine, and hydrogen abstraction reactions by each of the two radicals from the different sites of the two molecules were studied at the B3LYP/6‐31G(d,p), B3LYP/AUG‐cc‐pVDZ and BHandHLYP/AUG‐cc‐pVDZ levels of density functional theory. Effect of solvation in aqueous media on the reactions was studied at all these levels of theory using single point energy calculations using the polarizable continuum model. The present study shows that whereas the OH radical would abstract H atoms from the various sites of cytosine and thymine efficiently, the OOH radical would have poor reactivity in this regard. The OH radical is also predicted to be much more reactive than the OOH radical with regard to addition reactions at the C5 and C6 sites of both thymine and cytosine, though the OOH radical is also predicted to have significant reactivity in this respect. © 2012 Wiley Periodicals, Inc.     

Direct MP2 molecular dynamics studies of H atom reaction with CD4 and CH4

The mechanism and dynamics of the H + CD₄ → CD₃ + HD (I) and H + CH₄ → CH₃ + H₂ (II) reactions have been investigated by electronic structure methods. The minimum‐energy path and vibrational frequencies along the intrinsic reaction coordinate are calculated at MP2/cc‐pVDZ level. Energy distributions of the products are also obtained by the direct classical trajectory calculations at the MP2/ cc‐pVDZ level. It is found that most of the available energy appears as product translational energy, and very little of the available energy is partitioned into internal excitation of the HD (H₂) product for reaction I (II), which is in agreement with the experimental evidence. The results indicate that the experimental results could be reproduced by the direct MP2 molecular dynamics calculations. The rotational state distributions of the products show the HD (H₂) products are formed with lower rotational quantum numbers than the CD₃ (CH₃) products. © 2011 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

Hydrogen atom abstraction reactions of the sugar moiety of 2′‐deoxyguanosine with an OH radical: A quantum chemical study

Mechanisms of hydrogen atom abstraction reactions of the sugar moiety of 2′‐deoxyguanosine with an OH radical were investigated using the B3LYP and BHandHLYP functionals of density functional theory and the second order Møller–Plesset Perturbation (MP2) theory in gas phase and aqueous media. The 6‐31+G* and AUG‐cc‐pVDZ basis sets were used. Gibbs free barrier energies and rate constants of the reactions in aqueous media suggest that an OH radical would abstract the hydrogen atoms of the sugar moiety of 2′‐deoxyguanosine in the following order of preference: H5′ ≈ H5″ > H3′ > H4′ > H1′ ≈ H2′ > H2″, the rate constant for H5′ abstraction being 10³–10⁵ times greater than that for H2″ at the different levels of theory. Relative stabilities of the different deoxyribose radicals are also discussed. The most and least favored hydrogen abstraction reactions found here are in agreement with experimental observation. © 2010 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

An ab initio study of water clusters in gas phase and bulk aqueous media: (H2O)n,n=1–12

Geometries of several clusters of water molecules including single minimum energy structures of n‐mers (n=1–5), several hexamers and two structures of each of heptamer to decamer derived from hexamer cage and hexamer prism were optimized. One structural form of each of 11‐mer and 12‐mer were also studied. The geometry optimization calculations were performed at the RHF/6‐311G* level for all the cases and at the MP2/6‐311++G** level for some selected cases. The optimized cluster geometries were used to calculate total energies of the clusters in gas phase employing the B3LYP density functional method and the 6‐311G* basis set. Frequency analysis was carried out in all the cases to ensure that the optimized geometries corresponded to total energy minima. Zero‐point and thermal free energy corrections were applied for comparison of energies of certain hexamers. The optimized cluster geometries were used to solvate the clusters in bulk water using the polarized continuum model (PCM) of the self‐consistent reaction field (SCRF) theory, the 6‐311G* basis set, and the B3LYP density functional method. For the cases for which MP2/6‐311++G** geometry optimization was performed, solvation calculations in water were also carried out using the B3LYP density functional method, the 6‐311++G** basis set, and the PCM model of SCRF theory, besides the corresponding gas‐phase calculations. It is found that the cage form of water hexamer cluster is most stable in gas phase among the different hexamers, which is in agreement with the earlier theoretical and experimental results. Further, use of a newly defined relative population index (RPI) in terms of successive total energy differences per water molecule for different cluster sizes suggests that stabilities of trimers, hexamers, and nonamers in gas phase and those of hexamers and nonamers in bulk water would be favored while those of pentamer and decamer in both the phases would be relatively disfavored. © 2001 John Wiley & Sons, Inc. Int J Quant Chem 81: 90–104, 2001     

IR Low‐Temperature Matrix and ab Initio Study on β‐Alanine Conformers

For the first time the argon‐matrix low‐temperature FTIR spectra of β‐alanine are recorded. They reveal a quite complicated spectral pattern which suggests the presence of several β‐alanine conformers in the matrix. To interpret the spectra, the eighteen β‐alanine conformers, stable in the gas phase, are estimated at the B3LYP and MP2 levels combined with the aug‐cc‐pVDZ. Ten low‐energy structures are reoptimized at the QCISD/aug‐cc‐pVDZ and B3LYP and MP2 levels by using the aug‐cc‐pVTZ basis sets. Assignment of the experimental spectra is undertaken on the basis of the calculated B3LYP/aug‐cc‐pVDZ anharmonic IR frequencies as well as careful estimation of the conformer population. The presence of at least three β‐alanine conformers is demonstrated. The detailed analysis of IR spectra points to the possible presence of five additional β‐alanine conformers.     

Interaction of angelicin with DNA‐bases (III) DFT and ab initio second‐order Moeller–Plesset study

Angelicin geometry was optimized at MP2/6‐31+G(d,p) level and compared with X‐ray experimental data. The highest π‐electron density was found to be localized on C1? C2 and on C13? C14 as confirmed by the calculated bond length and bond order values. Spectrophptometric properties of angelicin were measured and compared with the computed within the TD‐DFT. Quantum chemical methods were used to study the interaction of angelicin, as a nonlinear furocoumarin, with DNA bases and base pairs. The interactions with DNA bases and base pairs were studied to shade more light on the nature of the intercalation binding forces between angelicin and DNA. Comparing computed electronic properties of angelicin with that of linear psoralens show that the former is a weaker intercalator. The geometry of complexes of angelicin with adenine, thymine, adenine–thymine base pair, cytocine, guanine as well as cytocine–guanine base pair have been optimized in two main orientations, planar and stacked, at the levels of B3LYP/cc‐pVDZ, MP2/6‐31G(d,p) and MP2/cc‐pVDZ. Effect of vertical distance and rotational angle between the stacked molecules on the interaction energy were investigated by the aforementioned methods in gas phase and water media. It was found that ab initio methods which account for the electron correlation effects are the minimum level for studying the noncovalent interactions. © 2009 Wiley Periodicals, Inc. Int J Quantum Chem, 2010     

Photo‐Induced Hydrogen Exchange Reaction between Methanol and Glyoxal: Formation of Hydroxyketene

We study the structure and photochemistry of the glyoxal–methanol system (G–MeOH) by means of FTIR matrix isolation spectroscopy and ab initio calculations. The FTIR spectra show that the non‐hydrogen‐bonded complex, G–MeOH‐1, is present in an inert environment of solid argon. MP2/aug‐cc‐pVDZ calculations indicate that G–MeOH‐1 is the most stable complex among the five optimized structures. The interaction energy partitioned according to the symmetry‐adapted perturbation theory (SAPT) scheme demonstrates that the dispersion energy gives a larger contribution to the stabilization of a non‐hydrogen‐bonded G–MeOH‐1 complex than compared to the hydrogen‐bonded ones. The irradiation of G–MeOH‐1 with the filtered output of a mercury lamp (λ>370 nm) leads to its photo‐conversion into the hydroxyketene–methanol complex HK–MeOH‐1. The identity of HK–MeOH‐1 is confirmed by both FTIR spectroscopy and MP2/aug‐cc‐pVDZ calculations. An experiment with deuterated methanol (CH₃OD) evidences that hydroxyketene is formed in a photo‐induced hydrogen exchange reaction between glyoxal and methanol. The pathway for the photo‐conversion of G–MeOH‐1 to HK–MeOH‐1 is studied by a coupled‐cluster method [CR–CC(2,3)]. The calculations confirm our experimental findings that the reaction proceeds via hydrogen atom exchange between the OH group of methanol and CH group of glyoxal.     

Mechanisms and kinetics for preparing carbohydrazide by reacting dimethyl carbonate with hydrazine: A theoretical study

The mechanism and kinetic modeling for preparing carbohydrazide from dimethyl carbonate and hydrazine has been declared. The geometries of all the stationary points (reactants, intermediates, transition states, and products) are optimized by using the B3LYP method with the cc‐pVDZ basis set, and the harmonic vibrational frequencies as well as infrared intensities are predicted with the same method. The minimum‐energy paths are obtained by using the intrinsic reaction coordinate (IRC) theory at the B3LYP/cc‐pVDZ level of theory with the step length 0.02 (amu)^1/2·bohr. The rate constants are evaluated by using the TST, TST/Eckart, and RRKM (T)/Eckart methods. © 2008 Wiley Periodicals, Inc. Int J Quantum Chem, 2008     

硫的取代对基态硫代乙酸单分子解离反应影响的理论研究

应用B3LYP方法,结合6-31G**、cc-pVDZ、aug-cc-pVDZ和cc-pVTZ基组对硫代乙酸的两种异构体CH3C(O)SH和CH3C(S)OH在基态势能面上的9个单分子反应进行了研究。本文计算预测硫代乙酸主要以CH3C(O)SH的形式存在,两种异构体均以顺式构象为优势构象。通过对比CH3C(O)SH、CH3C(S)OH和 CH3C(O)OH的反应性差异,我们可以得出结论：CH3C(O)OH中-OH基团的O被S取代后,只有当-SH作为一个整体参加反应时才对分子解离过程有较大影响;而C=O或C=S对反应性影响较小。     

An ab initio study of excited states of guanine in the gas phase and aqueous media: Electronic transitions and mechanism of spectral oscillations

Ground state geometries of the four tautomeric forms keto‐N9H, keto‐N7H, enol‐N9H, and enol‐N7H of guanine were optimized in the gas phase at the RHF level using a mixed basis set consisting of the 4‐31G basis set for all the atoms except the nitrogen atom of the amino group for which the 6‐311+G* basis set was used. These calculations were also extended to hydrogen‐bonded complexes of three water molecules with each of the keto‐N9H (G9‐3W) and keto‐N7H (G7‐3W) forms of guanine. Relative stabilities of the four above‐mentioned tautomers of guanine as well as those of G9‐3W and G7‐3W complexes in the ground state in the gas phase were studied employing the MP2 correlation correction. In aqueous solution, relative stabilities of these systems were studied using the MP2 correlation correction and polarized continuum model (PCM) or the isodensity surface polarized continuum model (IPCM) of the self‐consistent reaction field (SCRF) theory. Geometry optimization in the gas phase at the RHF level using the 6‐31+G* basis set for all atoms and the solvation calculations in water at the MP2 level using the same basis set were also carried out for the nonplanar keto‐N9H and keto‐N7H forms of guanine. Thus, it is shown that among the different tautomers of guanine, the keto‐N7H form is most stable in the gas phase, while the keto‐N9H form is most stable in aqueous solution. It appears that both the keto‐N9H and keto‐N7H forms of guanine would be present in the ground state, particularly near the aqueous solution–air interface. Vertical excitation and excited state geometry optimization calculations were performed using configuration interaction involving single electron excitation (CIS). It is found that the absorption spectrum of guanine would arise mainly due to its keto‐N9H form but the keto‐N7H form of the same would also make some contribution to it. The enol‐N9H and enol‐N7H forms of the molecule are not expected to occur in appreciable abundance in the gas phase or aqueous media. The normal fluorescence spectrum of guanine in aqueous solution with a peak near 332 nm seems to originate from the lowest singlet excited state of the keto‐N7H form of the molecule while the fluorescence of oxygen‐rich aqueous solutions of guanine with a peak near 450 nm appears to originate from the lowest singlet excited state of the keto‐N9H form of the molecule. The origin of the slow damped spectral oscillation observed in the absorption spectrum of guanine has been explained. © 2000 John Wiley & Sons, Inc. J Comput Chem 21: 826–846, 2000     

Theoretical study of spectroscopic and molecular properties of several low‐lying electronic states of CO molecule

The potential energy curves (PECs) of eight low‐lying electronic states (X¹Σ⁺, a³Π, a′³Σ⁺, d³Δ, e³Σ^?, A¹Π, I¹Σ^?, and D¹Δ) of the carbon monoxide molecule have been studied by an ab initio quantum chemical method. The calculations have been performed using the complete active space self‐consistent field method, which is followed by the valence internally contracted multireference configuration interaction (MRCI) approach in combination with the correlation‐consistent aug‐cc‐pV5Z basis set. The effects on the PECs by the core‐valence correlation and relativistic corrections are included. The way to consider the relativistic corrections is to use the third‐order Douglas–Kroll Hamiltonian approximation at the level of a cc‐pV5Z basis set. Core‐valence correlation corrections are performed using the cc‐pCVQZ basis set. To obtain more reliable results, the PECs determined by the MRCI calculations are corrected for size‐extensivity errors by means of the Davidson modification (MRCI+Q). The spectroscopic parameters (D_e, T_e, R_e, ω_e, ω_ex_e, ω_ey_e, B_e, α_e, and γ_e) of these electronic states are calculated using these PECs. The spectroscopic parameters are compared with those reported in the literature. Using the Breit–Pauli operator, the spin–orbit coupling effect on the spectroscopic parameters is discussed for the a³Π electronic state. With the PECs obtained by the MRCI+Q/aug‐cc‐pV5Z+CV+DK calculations, the complete vibrational states of each electronic state have been determined. The vibrational manifolds have been calculated for each vibrational state of each electronic state. The vibrational level G(ν), inertial rotation constant B_ν, and centrifugal distortion constant D_ν of the first 20 vibrational states when the rotational quantum number J equals zero are reported and compared with the experimental data. Comparison with the measurements demonstrates that the present spectroscopic parameters and molecular constants determined by the MRCI+Q/aug‐cc‐pV5Z+CV+DK calculations are both reliable and accurate. © 2012 Wiley Periodicals, Inc.     

New parallel algorithm for MP2 energy gradient calculations

A new parallel algorithm has been developed for calculating the analytic energy derivatives of full accuracy second order Møller‐Plesset perturbation theory (MP2). Its main projected application is the optimization of geometries of large molecules, in which noncovalent interactions play a significant role. The algorithm is based on the two‐step MP2 energy calculation algorithm developed recently and implemented into the quantum chemistry program, GAMESS. Timings are presented for test calculations on taxol (C₄₇H₅₁NO₁₄) with the 6‐31G and 6‐31G(d) basis sets (660 and 1032 basis functions, 328 correlated electrons) and luciferin (C₁₁H₈N₂O₃S₂) with aug‐cc‐pVDZ and aug‐cc‐pVTZ (530 and 1198 basis functions, 92 correlated electrons). The taxol 6‐31G(d) calculations are also performed with up to 80 CPU cores. The results demonstrate the high parallel efficiency of the program. © 2007 Wiley Periodicals, Inc. J Comput Chem, 2007     

A theoretical analysis of substituted aromatic compounds

The substituents ? CH₃, ? F, ? NO₂, ? OCH₃, and ? CH₂?CH₂ were placed at the ortho, meta, and para positions on the aromatic molecules aniline, benzaldehdye, nitrobenzene, and phenol. MMFF94, AM1, B3LYP, M06, M06‐2X, ωB97X, ωB97X‐d, and RI‐MP2 using cc‐pVDZ and cc‐pVTZ and CCSD(T) with cc‐pVDZ basis sets were used to calculate the geometries and energies of all regiomers of the molecules. Relative energies of the ortho and meta regiomers relative to the para regiomers were calculated and compared to the CCSD(T) values. A good basis set correlation between cc‐pVDZ and cc‐pVTZ was observed in RI‐MP2. Overall, RI‐MP2 gave the best correlation with the CCSD(T) results. All of the hybrid functionals showed similar accuracy and could effectively describe the intramolecular hydrogen‐bonding interactions of these compounds. The methoxy group at the para position in methoxyaniline, methoxyphenol, methoxynitrobenzene, and methoxybenzaldehyde was rotated around the phenyl‐O bond. HF, along with the cc‐pVDZ basis with the other methods, generated inaccurate energy profiles for p‐methoxyphenol. For the density functional theory methods, it was necessary to use improved grids to get smooth curves. © 2013 Wiley Periodicals, Inc.     

Electronic properties of the low‐lying spin states of dimethylnitrosamine coordinated to Fe(III) heme models: An ab initio study

The unusual O‐coordination mode of nitrosamines to Fe(III) heme models has been observed in the bis(dimethylnitrosamine)(meso‐tetraphenylporphyrinate)iron(III) cation. For the first time, this latter as well as the simpler bis(dimethylnitrosamine)(porphinate)iron(III) heme model cations have been studied through ab initio methods. The sextet, quartet, and doublet spin states of both cations have been studied through single‐point calculations based on the experimental (X‐ray) geometry. Their energies, charges, and spin densities have been analyzed. The obtained results (at the UHF/cc‐pVDZ and ROHF/cc‐pVDZ levels) indicate that the peripheral benzene rings are of secondary importance for the coordination of dimethylnitrosamine to the Fe(III) porphyrin core. The obtained energy ordering is sextet < quartet < doublet, at all computational levels. The UHF, ROHF, and UMP2 results indicate an excess of alpha spin density around the Fe atom, a low covalency for the Fe? O bond and a substantial charge transfer to the Fe atom. Our best estimates [obtained at ROMP2 level with the mixed cc‐pVDZ/cc‐pVTZ‐DK(Fe) basis set] for the energy differences (in eV) between the three spin states considered are 0.929 for the sextet‐quartet gap and 0.812 for the quartet‐doublet gap, which indicate that the spin crossover (at room temperature) is very unlikely. These results represent the substantial decrease in the uncorrelated values. The implications of spin contaminations at the UHF and UMP2 levels for subsequent geometry optimizations to be performed in the smaller cation have also been discussed. © 2013 Wiley Periodicals, Inc.     

A New Glimpse into the CO2‐Philicity of Carbonyl Compounds

We report a theoretical study on non‐conventional structures of 1:1 complexes between carbon dioxide and carbonyl compounds. These structures have never been reported before but are relevant for understanding the solubility of carbonyl compounds in supercritical CO₂. The work is based on the results of ab initio calculations at the MP2 and CCSD(T) levels using aug‐cc‐pVDZ and aug‐cc‐pVTZ basis sets. Investigated systems include aldehydes, ketones and esters, together with some fluorinated derivatives. The results are interpreted in terms of natural bond orbital analyses. Harmonic vibrational frequency calculations have also been done in order to compare them with available experimental data. We show for the first time that complexes where CO₂ behaves globally as a Lewis base are stable in the case of ketones and esters, but not in the case of aldehydes, and their stability is similar to that of traditional complexes in which CO₂ behaves as a Lewis acid. This finding considerably modifies the concept of CO₂‐philicity and may have important ramifications in the development of green reactions in supercritical CO₂.