Aza‐Diels–Alder reaction between cyclopentadiene and protonated N‐phenylethyliminoacetates of 8‐phenylmenthol and 8‐phenylneomenthol: a density functional theory study

The cycloaddition between glyoxylate imines possessing two chiral auxiliaries, N‐(R)‐ or N‐(S)‐1‐phenylethyl and 8‐phenylmenthyl or 8‐phenylneomenthyl, and cyclopentadiene is described. Computational calculations using density functional theory with the Becke, three‐parameter, Lee–Yang–Parr functional and the 6‐31G(d) basis set were performed to better understand the highly diastereoselective mechanism and the exo‐selectivity observed experimentally for these ionic aza‐Diels–Alder reactions. Copyright © 2011 John Wiley & Sons, Ltd.     

Photophysical properties of a series of 4‐aryl substituted 1,4‐dihydropyridines

In this article, a series of Hantzsch 1,4‐dihydropyridines with different substituted aryl groups were synthesized and its spectral data obtained by UV–Vis absorption and fluorescence emission spectroscopies in solution. The dihydropyridines present absorption located around 350 nm and fluorescence emission in the blue–green region. A higher Stokes’ shift could be observed for the derivative 3b because of an intramolecular charge transfer in the excited state from the dimethylaniline to the dihydropyridine chromophores, which was corroborated by a linear relation of the fluorescence maxima (ν_max) versus the solvent polarity function (Δf) from the Lippert–Mataga correlation. A comparison between the experimental data and time‐dependent density functional theory‐polarizable continuum model calculations of the vertical transitions was performed to help on the elucidation of the photophysics of these compounds. For these calculations, the S₀ and S₁ states were optimized using Becke, three‐parameter, Lee–Yang–Parr/6‐31 G* and Configuration Interaction Singles/6‐31 G*, respectively. The predicted absorption maxima are in good agreement with the experimental; however, the theoretical fluorescence emission maxima do not match the experimental, which means that the excited specie cannot be related to neither a locally excited state nor to an aromatized structure. Copyright © 2012 John Wiley & Sons, Ltd.     

Microsolvation of morpholine,a bidentate base – the importance of cooperativity

The structure, relative energies, and bonding in morpholine(water)_n (n = 1–4) clusters have been investigated at the Becke, three‐parameter, Lee–Yang–Parr/6–311 + G(3df,2p)//Becke, three‐parameter, Lee–Yang–Parr/6–311 + G(d,p) level of theory. Cooperative effects have been analyzed through the use of structural, energetic, and electron density indexes. Our analysis shows that these effects are crucial to trace the relative stability of the complexes formed. In all cases water molecules prefer to self‐associate forming chains in which each individual molecule behaves as a hydrogen bond (HB) donor and HB acceptor. The chain so formed behaves in turn as HB donor and HB acceptor with respect to morpholine, being the most stable arrangements those in which the NH group of morpholine behaves simultaneously as HB donor and HB acceptor. Higher in energy lie complexes in which the HB acceptor continues to be the NH group, but the HB donor is a CH group, or alternative structures in which the HB acceptor is the ether‐like oxygen of morpholine and the HB donor its NH group. Cooperativity increases with the number of solvent molecules, but there is a clear attenuation effect. Thus, whereas the additive interaction energy on going from dihydrated to trihydrated species increases by a factor of 3, this increase is about half on going from trihydrated to tetrahydrated complexes. Copyright © 2012 John Wiley & Sons, Ltd.     

Density functional theory study and vibrational analysis of FT‐IR and FT‐Raman spectra of N‐hydroxyphthalimide

The Fourier transform infrared (FT‐IR) spectrum of N‐hydroxyphthalimide has been recorded in the range of 4000–400 cm⁻¹, and the Fourier transform Raman (FT‐Raman) spectrum of N‐hydroxyphthalimide has been recorded in the range of 4000–50 cm⁻¹. With the hope of providing more and effective information on the fundamental vibrations, the Density Functional Theory (DFT)‐Becke3‐Lee‐Yang‐Parr (B3LYP) level with 6‐31G* basis set has been employed in quantum chemical analysis, and normal coordinate analysis has been performed on N‐hydroxyphthalimide by assuming C_s symmetry. The computational wavenumbers are in good agreement with the observed results. The theoretical spectra obtained along with intensity data agree well with the observed spectra. Copyright © 2009 John Wiley & Sons, Ltd.     

A Curtin–Hammett pentamethylene chain symmetrization process in the Bergman cyclization of an 11‐membered ring enediyne

Pentamethylene chain conformational effects for the Bergman cyclization of the 11‐membered ring enediyne, (3Z)‐3‐cycloundecene‐1,5‐diyne, 2, are examined theoretically with unrestricted Becke, three‐parameter, Lee–Yang–Parr/6‐31 G(d,p) calculations. A C₁ symmetric enediyne conformation was found to be the global minimum, where its nonsymmetric pentamethylene chain prevented π‐orbital alignment of the acetylene groups for C–C σ bond product formation. The Bergman cyclization of 2 was found to be conformationally dependent. In a Curtin–Hammett type process, the C₁ symmetric 2 inverts to one of the C_S or C₂ symmetric conformers required for the Bergman cyclization, which produced a C_S or C₂ symmetric 1,4‐diradical intermediate. The activation energy for the cyclization is slightly higher to reach the C₂ symmetry diradical compared with the C_S symmetry diradical. Copyright © 2012 John Wiley & Sons, Ltd.     

Theoretical study of keto–enol tautomerism by quantum mechanical calculations (the QM/MC/FEP method)

In this study, the tautomeric equilibrium between the keto and enol forms has been studied for five typical ketones and aldehydes: i‐butanal, acetaldehyde, acetone, acetylacetone, and dimedone. The level of theory used in the gas‐phase calculation was Becke, three‐parameter, Lee–Yang–Parr/6‐311G(d,p)//Becke, three‐parameter, Lee–Yang–Parr/6‐31G(d). The free energies of solvation were included in the calculation by using the free‐energy perturbation method based on Monte Carlo simulation, that is, the quantum mechanical/Monte Carlo/free‐energy perturbation method. Three different models, incorporating no‐water, one‐water, and two‐waters, were adopted. The results showed that in the gas phase the addition of water molecules to the reaction mechanism caused the activation barriers (ΔG^?_gas) to decrease by half relative to the water‐free mechanism, but there was no effect on the relative difference in free energy, ΔG_gas. The solvation effects (ΔG_sol), based on quantum mechanical/Monte Carlo/free‐energy perturbation calculations, were added to those of the gas‐phase results of the one‐water and two‐waters models. The two‐waters model produced values that were very consistent with the experimental data for all of the tautomers. The differences in the relative Gibbs free energy (ΔG_rxn) were less than 1.0 kcal mol^–1. In summary, the inclusion of solvent molecules in gas‐phase calculations plays a very important role in producing results consistent with experimental data. Copyright © 2012 John Wiley & Sons, Ltd.     

Theoretical study of energetic trinitromethyl‐substituted tetrazole and tetrazine derivatives

Density functional theory method was used to study the heats of formation, energetic properties, and thermal stability for a series of trinitromethyl‐substituted tetrazole and tetrazine derivatives with different substituents. It is found that the group ―NO₂, ―NHNO₂, or ―NF₂ play a very important role in increasing the heats of formation of the derivatives. The calculated detonation velocities and pressures indicate that the group ―CF₂NF₂, ―NHNO₂, ―1H‐tetrazolyl, ―2H‐tetrazolyl, or ―1,2,4,5‐tetrazinyl is an effective structural unit for enhancing their detonation performance. An analysis of the bond dissociation energies for several relatively weak bonds indicates that incorporating the group ―NHNO₂ and ―NH₂ into parent ring decreases their thermal stability. Considering the detonation performance and thermal stability, 37 compounds may be considered as the potential high‐energy compounds. Their oxygen balances are close to zero. These results provide basic information for the molecular design of novel high‐energy compounds. Copyright © 2013 John Wiley & Sons, Ltd.     

Theoretical studies on nitrogen‐rich furoxan‐based heterocyclic derivatives

The density functional theory method was used to study the heats of formation (HOFs), energetic properties, electronic structure of a series of 4,4″‐dinitro(3,3′:4′,3′′)tris([1,2,5]oxadiazole)‐2′‐oxide (3,4‐bis[4′‐nitrofurazan‐3′‐yl]furoxan) derivatives. The results show that the substitution of the nitro group is very useful for improving their HOFs and detonation performances. The HOFs of the title compounds are all positive and larger than those of 1,3,5‐trinitro‐1,3,5‐triazinane and 1,3,5,7‐tetranitro‐1,3,5,7‐tetrazocane. The analysis of oxygen balance shows that the studied compounds need the oxygen in the explosive. Compound A1 has larger detonation velocity and detonation pressure than those of 1,3,5,7‐tetranitro‐1,3,5,7‐tetrazocane and can be regarded as a potential candidate for high‐energy compounds because of the moderate heat of detonation, high density, and high N. In addition, the energy gaps between the highest occupied molecular orbital and lowest unoccupied molecular orbital of the studied compounds are further investigated.     

Density functional theory study of the bis‐3‐benzo‐crown ethers and their complexes with alkali metal cations Na+, K+, Rb+ and Cs+

The binding interactions of bis‐3‐benzo‐15‐crown‐5 ethers and bis‐3‐benzo‐18‐crown‐6 ethers (neutral hosts) with a series of alkali metal cations Na⁺, K⁺, Rb⁺ and Cs⁺ (charged guests) were investigated using quantum chemical density functional theory. Different optimized structures, binding energies and various thermodynamic parameters of free crown ethers and their metal cation complexes were obtained based on the Becke, three‐parameter, Lee–Yang–Parr functional using mixed basis set (C, H, O, Na⁺ and K⁺ using 6‐31 g, and the heavier cation Rb⁺ and Cs⁺ using effective core potentials). Natural bond orbital analysis is conducted on the optimized geometric structures. The main types of driving force host–guest interactions are investigated. The electron donating O offers a lone pair of electrons to the contacting LP* (1‐center valence antibond lone pair) orbitals of metal cations. The bis‐3‐benzocrown ethers are assumed to have sandwich‐like conformations, considering the binding energies to gauge the exact interactions with alkali cations. It is found that there are two different types of complexes: one is a tight ion pair and the other is a separated ion pair. Copyright © 2011 John Wiley & Sons, Ltd.     

Investigation of cation–anion interaction in 1‐(2‐hydroxyethyl)‐3‐methylimidazolium‐based ion pairs by density functional theory calculations and experiments

Gas‐phase structure, hydrogen bonding, and cation–anion interactions of a series of 1‐(2‐hydroxyethyl)‐3‐methylimidazolium ([HOEMIm]⁺)‐based ionic liquids (hereafter called hydroxyl ILs) with different anions (X = [NTf₂]^–, [PF₆]^–, [ClO₄]^–, [BF₄]^–, [DCA]^–, [NO₃]^–, [AC]^– and [Cl]^–), as well as 1‐ethyl‐3‐methylimizolium ([EMIm]⁺)‐based ionic liquids (hereafter called nonhydroxyl ILs), were investigated by density functional theory calculations and experiments. Electrostatic potential surfaces and optimized structures of isolated ions, and ion pairs of all ILs have been obtained through calculations at the Becke, three‐parameter, Lee–Yang–Parr/6‐31 + G(d,p) level and their hydrogen bonding behavior was further studied by the polarity and Kamlet–Taft Parameters, and ¹H‐NMR analysis. In [EMIm]⁺‐based nonhydroxyl ILs, hydrogen bonding preferred to be formed between anions and C2–H on the imidazolium ring, while in [HOEMIm]⁺‐based hydroxyl ILs, it was replaced by a much stronger one that preferably formed between anions and OH. The O–H···X hydrogen bonding is much more anion‐dependent than the C2–H···X, and it is weakened when the anion is changed from [AC]^– to [NTf₂]^–. The different interaction between [HOEMIm]⁺ and variable anion involving O–H···X hydrogen bonding resulted in significant effect on their bulk phase properties such as ¹H‐NMR shift, polarity and hydrogen‐bond donor ability (acidity, α). Copyright © 2011 John Wiley & Sons, Ltd.     

Spectroscopic and computational studies on self‐assembly complexes of bis(pyrrol‐2‐ ylmethyleneamine) ligands linked by alkyl spacers with Cu(II)

Bis(pyrrol‐2‐ylmethyleneamine) ligands and their mononuclear monomeric and dinuclear dimeric self‐assembly complexes with Cu(II) were investigated by means of IR and Raman spectroscopies and density functional theory. The ground‐state geometries were calculated by using the Becke Lee Yang Parr composite exchange‐correlation functional (B3LYP) and a combined basis set (LanL2DZ for Cu; 6–31G(d) for C, H, N), and they were compared with the single‐crystal X‐ray diffraction (XRD) structures. The DFT‐calculated Cu N bond lengths are generally higher by 0.001–0.040 Å than those determined through XRD. The vibrational spectra were also calculated at the same level of theory for the optimized geometries. The calculated wavenumbers were scaled by a uniform scaling factor and compared with the experimental fundamentals. The predicted spectra are in good agreement with the experimental ones with the deviations generally less than 30 cm⁻¹. In comparison with the spectra of the ligands, the coordination effect shifts the υ(CN) wavenumber by about 50 cm⁻¹ toward a lower value. Because of the weak intermolecular C H···Cu hydrogen bond, the Cu N stretching mode is shifted toward a lower wavenumber. Copyright © 2006 John Wiley & Sons, Ltd.     

Comparative theoretical studies of high energetic cyclic nitramines

Density functional theory studies on cyclic nitramines were performed at B3LYP/6‐311G(d,p) level. The crystal structures were obtained by molecular packing calculations. Heats of formation (HOFs) were predicted through designed isodesmic reactions. Results indicate that the value of HOF relates to the number of =N–NO₂ group and aza nitrogen atom and increases with the augment of the number of =N–NO₂ group and aza nitrogen atom for cyclic nitramines. Detonation performance was evaluated by using the Kamlet–Jacobs equations based on the calculated densities and HOFs. All the cyclic nitramines exhibit better detonation performance than 1,3,5‐trinitro‐1,3,5‐triazacyclohexane and 1,3,5,7‐tetranitro‐1,3,5,7‐tetraazacyclooctane. The stability of cyclic nitramines was investigated by the bond dissociation energies. The result shows that the increase of =N?NO₂ group or aza nitrogen atom reduces the stability of the title compounds. These results provide basic information for molecular design of novel high energetic density materials. Copyright © 2013 John Wiley & Sons, Ltd.     

Electronic and structural effects on the nonlinear optical behavior in π‐conjugated structure bis(4‐nitrophenyl) carbonate: a vibrational spectroscopic approach

Fourier transform (FT)‐Raman and infrared (IR) spectra of the nonlinear optical (NLO) material bis(4‐nitrophenyl) carbonate were recorded and analyzed. The geometry, first hyperpolarizability and harmonic vibrational wavenumbers were calculated with the help of Becke3–Lee–Yang–Parr density functional theory method. The detailed interpretation of the vibrational spectra was carried out with the aid of normal coordinate analysis following the scaled quantum mechanical force field methodology. The second‐order NLO properties of the molecule were studied by the Kurtz and Perry powder reflection technique. Stability of the molecule arising from hyperconjugative interactions leading to its NLO activity and charge delocalization were analyzed using natural bond orbital analysis. Copyright © 2010 John Wiley & Sons, Ltd.     

Density functional theory and time‐dependent density functional theory study on a series of iridium complexes with tetraphenylimidodiphosphinate ligand

A series of heteroleptic cyclometalated Ir(III) complexes for organic light‐emitting diodes (OLEDs) application have been investigated theoretically to explore their electronic structures and spectroscopic properties. The geometries, the electronic structures, the lowest‐lying singlet absorptions and triplet emissions of Ir(dfppy)₂(tpip), Ir(tfmppy)₂(tpip), and theoretically designed models of Ir(ppy)₂(tpip) were investigated with the density functional theory (DFT)‐based approaches, where ppy = 2‐phenylpyridine, dfppy = 4,6‐difluorophenylpyridine, tfmppy = 4‐trifluoromethylphenylpyridine, and tpip = tetraphenylimidodiphosphinate. Their structures in the ground and their excited states have been optimized at the DFT/Becke 3‐parameter Lee Yang Parr (B3LYP)/Los Alamos National Laboratory 2‐double‐z (LANL2DZ) and time‐dependent DFT/B3LYP/LANL2DZ levels, and the lowest absorptions and emissions were evaluated at B3LYP and M062X level of theory, respectively. Furthermore, the energy transfer mechanism together with the advantage of low efficiency roll‐off for these complexes also can be analyzed here. Copyright © 2013 John Wiley & Sons, Ltd.     

DFT studies on a high-energy density cage compound 1, 3, 5, 7, 9, 11-hexo(N(CH3)NO2)-2, 4, 6, 8, 10, 12-hexaazatetracyclo[5, 5, 0, 0, 0] dodecane

The infrared and Raman spectra, heat of formation (HOF) and thermodynamic properties were investigated by B3LYP/6-31G^** method for a new designed polynitro cage compound 1,3,5,7,9,11-hexo(N(CH₃)NO₂)-2,4,6,8,10,12-hexaazatetracyclo[5,5,0,0,0]dodecane. The detonation velocity (D) and pressure (P) were predicted by the Kamlet–Jacobs equations based on the theoretical density and condensed HOF. The bond dissociation energies and bond orders for the weakest bonds were analysed to investigate the thermal stability of the title compound. The computational result shows that the detonation velocity and pressure of the title compound are superior to those of hexahydro-1,3,5-trinitro-1,3,5-triazine (RDX), but inferior to those of 1,3,5,7-tetranitro-1,3,5,7-tetraazacyclooctane (HMX) and hexanitrohexaazaisowurtzitane (HNIW). And the analysis of thermal stability shows that the first step of pyrolysis is the rupture of the N₇–NO₂ bond. The crystal structure obtained by molecular mechanics belongs to the P2₁ space group, with the lattice parameters Z = 2, a = 11.8246 Å, b = 10.4632 Å, c = 15.9713 Å, ρ = 1.98 g cm^?3.     

A theoretical study of hemiacetal formation from the reaction of methanol with derivatives of CX3CHO (X = H,F, Cl,Br and I)

A theoretical study of the hemiacetal formation reaction between methanol and CX₃CHO (X = H, F, Cl, Br, and I) has been carried out using density functional theory and Becke, three‐parameter, Lee–Yang–Parr/6‐311++G(d,p) computational methods. The stationary points of the reaction between the isolated molecules and the reaction catalyzed by an additional methanol molecule have been characterized. Because the final products present a stereogenic center, the potential autocatalysis of the reaction has been examined and also the possibility of spontaneous generation of chirality when the hemiacetal molecules are involved in the transition state structure. High barriers are found in the reaction between the isolated molecules that are reduced by the assistance of an additional molecule (methanol or hemiacetal product). The reactions catalyzed by the hemiacetal products show higher barriers than the one catalyzed by methanol. Copyright © 2012 John Wiley & Sons, Ltd.     

Experimental and theoretical studies of Raman spectroscopy on 4‐mercaptopyridine aqueous solution and 4‐mercaptopyridine/Ag complex system

Raman scattering and surface‐enhanced Raman scattering (SERS) have been used to study the behavior of 4‐mercaptopyridine (4‐Mpy) dissolved in water and adsorbed on silver mirrors. In order to gain the actual structure and the theoretical modes of the 4‐Mpy dissolved in water and adsorbed on the surface of silver mirror, ab initio calculation at the Hartree–Fock (HF) level and density functional theory (DFT) at Beck's three‐parameter Lee‐Yang‐Parr (B3LYP) level were performed to calculate the vibrational modes and wavenumbers. 4‐Mpy/2H₂O and 4‐Mpy/Ag complex systems were optimized, and then the corresponding Raman spectra were calculated and analyzed. Compared with the experimental results, the calculated results of 4‐Mpy and 4‐Mpy/2H₂O complex systems obtained from DFT method were more accurate. Among the results calculated with HF method, the one with three Ag atoms was economical, which took less computer time but gave equivalent results to those with more noumber of Ag atoms. Copyright © 2007 John Wiley & Sons, Ltd.     

Organic cyclic difluoramino‐nitramines: infrared and Raman spectroscopy of 3,3,7,7‐tetrakis(difluoramino)octahydro 1,5‐dinitro‐1,5‐diazocine (HNFX)

We present the first vibrational structure investigation of 3,3,7,7‐tetrakis(difluoramino)octahydro‐1,5‐dinitro‐ 1,5‐diazocine (HNFX)—and, more generally, of a member of the new class of gem‐bis(difluoramino)‐substituted heterocyclic nitramine energetic materials—using combined theoretical and experimental approaches. Optimized molecular structure and vibrational spectra of the Ci… symmetry conformer constituting the HNFX crystal were computed using density functional theory methods. Fourier transform infrared and Raman spectra of HNFX crystalline samples were also collected at ambient temperature and pressure. The average deviation of calculated structural parameters from X‐ray diffraction data is ∼1% at the B3LYP/6‐311 + + G(d,p) level of theory, suggesting the absence of significant molecular distortion induced by the crystal field. Very good agreement was found between simulated and measured spectra, allowing reliable assignment of the fundamental normal modes of vibration of the HNFX crystal. Detailed analysis of the normal modes of the C–(NF₂)₂ and N–NO₂ moieties was performed due to their critical importance in the initial steps of the molecular homolytic fragmentation process. Copyright © 2009 John Wiley & Sons, Ltd.     

Theoretical investigations on 4,4′,5,5′‐tetranitro‐2,2′‐1H,1′H‐2,2′‐biimidazole derivatives as potential nitrogen‐rich high energy materials

The ―NH₂, ―NO₂, ―NHNO₂, ―C(NO₂)₃ and ―CF(NO₂)₂ substitution derivatives of 4,4′,5,5′‐tetranitro‐2,2′‐1H,1′H‐2,2′‐biimidazole were studied at B3LYP/aug‐cc‐pVDZ level of density functional theory. The crystal structures were obtained by molecular mechanics (MM) methods. Detonation properties were evaluated using Kamlet–Jacobs equations based on the calculated density and heat of formation. The thermal stability of the title compounds was investigated via the energy gaps (?E_{LUMO ? HOMO}) predicted. Results show that molecules T5 (D = 10.85 km·s^?1, P = 57.94 GPa) and T6 (D = 9.22 km·s^?1, P = 39.21 GPa) with zero or positive oxygen balance are excellent candidates for high energy density oxidizers (HEDOs). All of them appear to be potential explosives compared with the famous ones, octahydro‐1,3,5,7‐tetranitro‐1,3,5,7‐tetraazocane (HMX, D = 8.96 km·s^?1, P = 35.96 GPa) and hexanitrohexaazaisowurtzitane (CL‐20, D = 9.38 km·s^?1, P = 42.00 GPa). In addition, bond dissociation energy calculation indicates that T5 and T6 are also the most thermally stable ones among the title compounds. Copyright © 2014 John Wiley & Sons, Ltd.