Effects of different basis sets and donor‐acceptor groups on linear and second‐order nonlinear optical properties and molecular frontier orbital energies

The calculation of molecular hyperpolarizability, molecular frontier orbital energies of some donor‐acceptor oxadiazoles ( 5a – f , 8a – f , and 9a – f ) have been investigated using ab initio methods and different basis sets. Ab initio optimizations were performed at the Hartree–Fock (HF) and density functional (Beckee‐3–Lee–Yang–Parr; B3LYP) levels of theory with 6‐31G basis set. The polarizability (<α>), anisotropy of polarizability (Δα), and ground‐state dipole moment (μ), first hyperpolarizability (β), and molecular frontier orbital (HOMO, highest occupied molecular orbital and LUMO, lowest unoccupied molecular orbital) energies of 5a – f , 8a – f , and 9a – f have been calculated at the HF and B3LYP methods with 6‐31G, 6‐31G(d), 6‐31+G(d), 6‐31++G(d,p), 6‐311G, 6‐311G(d), 6‐311+G(d), and 6‐311++G(d,p) basis sets. Also, the molecular hardness (η) and electronegativity (χ) parameters have been obtained using molecular frontier orbital energies. The <α>, Δα, μ, β, HOMO, LUMO energies, η and χ parameters have been investigated as dependence on the choice of method and basis set. The variation graphics of <α>, Δα, μ, β, η, and χ parameters using HF and B3LYP methods with different basis sets are presented. We have examined the frontier molecular orbital pictures of 5a – f , 8a – f , and 9a – f using B3LYP/6‐31++G(d,p) level. The 5a – f , 8a – f , and 9a – f display significant linear, second‐order molecular nonlinearity, and molecular parameters and provide the basis for future design of efficient nonlinear optical materials having the 1,3,4‐oxadiazole core. © 2009 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

New scale factors for harmonic vibrational frequencies using the B3LYP density functional method with the triple-zeta basis set 6-311+G(d,p)

We have calculated optimal frequency scaling factors for the B3LYP/ 6-311+G(d,p) method for fundamental vibrational frequencies on the basis of a set of 125 molecules. Using the new scaling factor, the vibrational frequencies calculated with the triple-zeta basis set 6-311+G(d,p) give significantly better accuracy than those calculated with the double-zeta 6-31G(d) basis set. Scale factors were also determined for low-frequency vibrations using the molecular set of 125 molecules and for zero-point energies using a smaller set of 40 molecules. We have studied the effect on the calculated vibrational frequencies for various combinations of diffuse and polarization functions added to the triple-zeta 6-311G basis set. The 6-311+G(d,p) basis set is found to give almost converged frequencies for most molecules, and we conclude that our optimum scaling factors are valid for the basis sets 6-311G(d,p) to 6-311++G(3df,3pd). The new scale factors are 0.9679 for vibrational frequencies, 1.0100 for low-frequency vibrations, and 0.9877 for zero-point vibrational energies.     

An experimental and theoretical approach to the molecular structure of 3‐{[4‐(3‐Methyl‐3‐phenyl‐cyclobutyl)‐thiazol‐2‐yl]‐hydrazono}‐1,3‐dihydro‐ indol‐2‐one

The title molecule, 3‐{[4‐(3‐methyl‐3‐phenyl‐cyclobutyl)‐thiazol‐2‐yl]‐hydrazono}‐1,3‐dihydro‐indol‐2‐one (C₂₂H₂₀N₄O₁S₁), was prepared and characterized by ¹H NMR, ¹³C NMR, IR, UV–visible, and single‐crystal X‐ray diffraction. The compound crystallizes in the monoclinic space group P21 with a = 8.3401(5), b = 5.6976(3), c = 20.8155(14) Å, and β = 95.144(5)°. Molecular geometry from X‐ray experiment and vibrational frequencies of the title compound in the ground state has been calculated using the Hartree–Fock with 6‐31G(d, p) and density functional method (B3LYP) with 6‐31G(d, p) and 6‐311G(d, p) basis sets, and compared with the experimental data. The calculated results show that optimized geometries can well reproduce the crystal structural parameters, and the theoretical vibrational frequencies values show good agreement with experimental data. Density functional theory calculations of the title compound and thermodynamic properties were performed at B3LYP/6‐31G(d, p) level of theory. © 2011 Wiley Periodicals, Inc. Int J Quantum Chem, 2012     

复合物HNO···H2O2内N-H···O蓝移氢键的理论研究

A theoretical study on the blue-shifted H-bond N-H…O and red-shifted H-bond O-H…O in the complexHNO…H_2O_2 was conducted by employment of both standard and counterpoise-corrected methods to calculate thegeometric structures and vibrational frequencies at the MP2/6-31G(d),MP2/6-31 G(d,p),MP2/6-311 q G(d,p),B3LYP/6-31G(d),B3LYP/6-31 G(d,p) and B3LYP/6-311 G(d,p) levels.In the H-bond N-H…O,the calcu-lated blue shift of N-H stretching frequency is in the vicinity of 120 cm~(-1) and this is indeed the largest theoreticalestimate of a blue shift in the X-H…Y H-bond ever reported in the literature.From the natural bond orbital analy-sis,the red-shifted H-bond O-H…O can be explained on the basis of the dominant role of the hyperconjugation.For the blue-shifted H-bond N-H…O,the hyperconjugation was inhibited due to the existence of significant elec-tron density redistribution effect,and the large blue shift of the N-H stretching frequency was prominently due tothe rehybridization of sp~n N-H hybrid orbital.     

Spectroscopic (FTIR and FT Raman) analysis and vibrational study on 2,3-dimethyl naphthalene using ab-initio HF and DFT calculations

A combined experimental and theoretical study on molecular and vibrational structure of 2,3-dimethyl naphthalene (2,3-DMN) has been undertaken in the present work. The FTIR and FT Raman spectra of 2,3-DMN were recorded in the region 4000-100 cm(-1). The optimized geometries were calculated by HF and DFT (B3LYP) methods with 6-31++G (d, p), 6-311G (d, p) and 6-311++G (d, p) basis sets. The harmonic vibrational frequencies, infrared intensities and Raman activities of the 2,3-DMN were evaluated with these methods. After scaling the computational wave numbers are in very good agreement with the experimental values. A detailed interpretation of the infrared and Raman spectra of 2,3-DMN is presented. The effects of substitution of methyl group on the molecule have also been discussed.     

FT-IR and FT-Raman investigation, computed vibrational intensity analysis and computed vibrational frequency analysis on m-Xylol using ab-initio HF and DFT calculations

The FT-IR and FT-Raman spectra of m-Xylol molecule have been recorded using Bruker IFS 66V spectrometer in the range 4000-100cm(-1). The molecular geometry and vibrational frequencies in the ground state are evaluated using the Hartree-fock (HF) and B3LYP with 6-31+G (d, p), 6-31++G (d, p) and 6-311++G (d, p) basis sets. The computed frequencies are scaled using a suitable scale factors to yield good agreement with the observed values. The HF and DFT analysis agree well with experimental observations. Comparison of the fundamental vibrational frequencies with calculated results by HF and B3LYP methods indicate that B3LYP/6-311++G (d, p) is superior to HF/6-31+G (d, p) for molecular vibrational problems. The complete data of this title compound provide some useful information for the study of substituted benzenes. The influences of Methyl groups on the geometry of benzene and its normal modes of vibrations have also been discussed.     

Quantum chemical studies on prototautomerism of 1H‐imidazo[4,5‐c]pyridine

The structural features of the 1H‐imidazo[4,5‐c]pyridine (ICPY) tautomers and homodimers of the most stable tautomers have been studied by quantum chemical methods. FTIR and Raman spectra of the ICPY were recorded in the range of 4000–60 cm^?1 and 3500–5 cm^?1. The predominant tautomer among four possible isomers of ICPY were determined. The optimized geometries and vibrational frequencies of possible ICPY tautomers and dimers were computed by B3LYP/DFT method with 6‐311++G(d,p) and 6‐31G(d) basis sets. All vibrational frequencies assigned in detail with the help of total energy distribution (TED) and isotopic shifts. ICPY dimeric forms were also characterized according to their hydrogen bonding interactions, and it has been found that the most stable ICPY homodimer establishes moderate strong N ? H …N type hydrogen bond. ¹H NMR, ¹³C NMR, and ¹⁵N NMR properties have been calculated for all tautomeric forms using the gauge independent atomic orbital (GIAO) method. © 2010 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

DFT study of the gas‐phase thermal decomposition kinetics of 2‐ethoxypyridine into 2‐pyridone

The mechanism of the gas‐phase elimination kinetics of 2‐ethoxypyridine has been studied through the electronic structure calculations using density functional methods: B3LYP/6‐31G(d,p), B3LYP/6‐31++G(d,p), B3PW91/6‐31G(d,p), B3PW91/6‐31++G(d,p), MPW1PW91/6‐31G(d,p), MPW1PW91/6‐31++G(d,p), PBEPBE/6‐31G(d,p), PBEPBE/6‐31++G(d,p), PBE1PBE1/6‐31G(d,p), and PBE1PBE1/6‐31++G(d,p). The elimination reaction of 2‐ethoxypyridine occurs through a six‐centered transition state geometry involving the pyridine nitrogen, the substituted carbon of the aromatic ring, the ethoxy oxygen, two carbons of the ethoxy group, and a hydrogen atom, which migrates from the ethoxy group to the nitrogen to give 2‐pyridone and ethylene. The reaction mechanism appears to occur with the participation of π‐electrons, similar to alkyl vinyl ether elimination reaction, with simultaneous ethylene formation and hydrogen migration to the pyridine nitrogen producing 2‐pyridone. © 2011 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

FT-IR and FT-Raman spectroscopic investigation, computed vibrational frequency analysis and IR intensity and Raman activity peak resemblance analysis on 2-nitroanisole using HF and DFT (B3LYP and B3PW91) calculations

Fourier-transform Raman and infrared spectra of 2-nitroanisole are recorded (4000-100 cm(-1)) and interpreted by comparison with respective theoretical spectra calculated using HF and DFT method. The geometrical parameters with C(S) symmetry, harmonic vibrational frequencies, infrared and Raman scattering intensities are determined using HF/6-311++G (d, p), B3LYP/6-311+G (d, p), B3LYP/6-311++G (d, p) and B3PW91/6-311++G (d, p) level of theories. A detailed vibrational spectral analysis has been carried out and assignments of the observed fundamental bands have been proposed on the basis of peak positions and relative intensities. The results of the calculations have been used to simulate IR and Raman spectra for the molecule that showed good agreement with the observed spectra. The SQM method, which implies multiple scaling of the DFT force fields has been shown superior to the uniform scaling approach. The vibrational frequencies and the infrared intensities of the C-H modes involved in back-donation and conjugation are also investigated.     

FT-IR, FT-Raman vibrational spectra and molecular structure investigation of 2-chloro-4-methylaniline: A combined experimental and theoretical study

In this work, the experimental and theoretical vibrational spectra of 2-chloro-4-methylaniline (2Cl4MA, C₇H₈NCl) were studied. FT-IR and FT-Raman spectra of 2Cl4MA in the liquid phase have been recorded in the region 4000–400 cm⁻¹ and 3500–50 cm⁻¹, respectively. The structural and spectroscopic data of the molecule in the ground state have been calculated by using Hartree-Fock (HF) and density functional method (B3LYP) with the 6-31G(d), 6-31G(d,p), 6-31+G(d,p), 6-31++G(d,p) and 6-311G(d), 6-311G(d,p), 6-311+G(d,p), 6-311++G(d,p) basis sets. The vibrational frequencies have been calculated and scaled values have been compared with experimental FT-IR and FT-Raman spectra. The observed and calculated frequencies are found to be in good agreement. The complete assignments were performed on the basis of the total energy distribution (TED) of the vibrational modes, calculated with scaled quantum mechanics (SQM) method. The DFT-B3LYP/6-311++G(d,p) calculations have been found more reliable than the ab initio HF/6-311++G(d,p) calculations for the vibrational study of 2Cl4MA. The optimized geometric parameters (bond lengths and bond angles) were compared with experimental values of aniline and p-methylaniline molecules.     

Theoretical calculations on the mechanisms of the gas‐phase elimination kinetics of 2‐chloro‐1‐phenylethane, 3‐chloro‐1‐phenylpropane, 4‐chloro‐1‐phenylbutane, 5‐chloro‐1‐phenylpentane,and their corresponding chloroalkanes: The effect of the phenyl ring

The kinetics and mechanisms of the dehydrochlorination of 2‐chloro‐1‐ phenylethane, 3‐chloro‐1‐phenylpropane, 4‐chloro‐1‐phenylbutane, 5‐chloro‐1‐phenylpentane, and their corresponding chloroalkanes were examined by means of electronic structure calculation using density functional theory methods B3LYP/6–31G(d,p), B3LYP/6–31++G(d,p), MPW1PW91/6–31G(d,p), MPW1PW91/6–31++G(d,p), PBEPBE/6–31G(d,p), and PBEPBE/6–31++G(d,p). The potential energy surface was investigated for the minimum energy path. Calculated enthalpies and energies of activation are in good agreement with experimental values using the MPW1PW91 and B3LYP methods. The transition state of these reactions is a four‐centered cyclic structure. The reported experimental results proposing neighboring group participation by the phenyl group was not supported by theoretical calculations. The rate‐determining process in these reactions is the breaking of Cl? C bond. The reactions are described as concerted moderately polar and nonsynchronous. © 2011 Wiley Peiodicals, Inc. Int J Chem Kinet 43: 292–302, 2011     

Theoretical Study on Gas—phase Pyrolytic Reactions of N—Ethyl,Isopropyl and N—t—Butyl Substituted 2—Aminopyrazine

Density functional theory (DFT) and ab initio methods were used to study gas‐phase pyrolytic reaction mechanisms of iV‐ethyl, N‐isopropyl and N‐t‐butyl substituted 2‐aminopyrazine at B3LYP/6–31G* and MP2/6–31G*, respectively. Single‐point energies of all optimized molecular geometries were calculated at B3LYP/6–311 + G(2d,p) level. Results show that the pyrolytic reactions were carried out through a unimolecular first‐order mechanism which were caused by the migration of atom H(17) via a six‐member ring transition state. The activation energies which were verified by vibrational analysis and correlated with zero‐point energies along the reaction channel at B3LYP/6–311 + G(2d,p) level were 252.02 kJ. mo^?1 (N‐ethyl substituted), 235.92 kJ‐mol^?1 (N‐t‐isopropyl substituted) and 234.27 kJ‐mol^?1 (N‐t‐butyl substituted), respectively. The results were in good agreement with available experimental data.     

Model Study of the Photochemical Rearrangement Pathways of 1,2,4‐Oxadiazole

The mechanisms of photochemical isomerization reactions are investigated theoretically by using a model system of 1,2,4‐ oxadiazole with the CAS(14,9)/6‐311G(d) and MP2‐CAS‐(14,9)/ 6‐311++G(3df,3pd)//CAS(14,9)/6‐311G(d) methods. Three reaction pathways are examined, including 1) the direct mechanism, 2) the ring contraction–ring expansion mechanism, and 3) the internal cyclization–isomerization mechanism, which lead to two types of photoisomers. The theoretical findings suggest that conical intersections play a crucial role in the photorearrangement of 1,2,4‐oxadiazoles. These model investigations also indicate that the preferred reaction route for 1,2,4‐oxadiazole, which leads to phototransposition products, is as follows: reactant → Franck‐Condon region → conical intersection → photoproduct. In other words, the direct mechanism is a one‐step process that has no barrier. These theoretical results agree with the available experimental observations.     

Relaxed energetic maps of κ‐carrabiose: A DFT study

The B3LYP density function was used with the 6‐31G(d) basis set to perform relaxed energetic contour maps of the charged form of κ‐carrabiose in the gas phase and for the neutral form first in the gas phase and then by simulating the presence of water as solvent using the Onsager model. Only one starting conformation has been considered to perform all the calculations. Rigid energetic maps have been then constructed either by addition of diffuse or polarization functions to the basis set obtaining in that way 6‐31+G(d)//6‐31G(d), 6‐31+G(d,p)//6‐31G(d), and 6‐311++G(d,p)//6‐31G(d) energetic maps that have been carefully examined. The obtained structures corresponding to the lower energy conformers have been then fully optimized using different basis sets with the B3LYP method, a reversion in term of energy has been observed for the two first minima in the case of the charged disaccharide in the gas phase, this was attributed to the large grid of 30° that could lead to the exclusion of an intermediate value corresponding to the real minimum of energy. We thus suggest that after establishing potential energy maps it is essential to proceed to full optimizations of the lower energy conformers. Calculations using the more accurate correlated method MP2 with the 6‐31G(d) basis set have also been performed for conformers of the two disaccharides in the gas phase. © 2010 Wiley Periodicals, Inc. J Comput Chem 2010     

Scaling factors for vibrational frequencies and zero-point vibrational energies of some recently developed exchange-correlation functionals

The scaling factors for the vibrational frequencies and zero-point vibrational energies evaluated at various combinations of recently developed exchange-correlation functionals and various basis sets are reported. The exchange-correlation functionals considered are B972, B98, HCTH, OLYP, O3LYP, G96LYP, PBE0 and VSXC functionals; the basis sets employed are 3-21G, 6-31G*, 6-31G**, 6-31+G, 6-311G*, 6-311G**, 6-311G(df,p), 6-311+G(df,p), cc-pVDZ and aug-cc-pVDZ. The experimental harmonic frequencies of 122 small molecules and the zero-point vibrational energies of 39 small molecules are used to determine the scaling factors through the least-square fitting procedure. It was found that the scaling factors do not depend significantly on the basis sets considered. The vibrational frequency scaling factors evaluated by using the B98 and PBE0 functionals are recommended on the basis of smallest root mean square error. The zero-point vibrational energy scaling factors evaluated from the B972 functional with Pople's double-zeta basis set and the HCTH functional with Pople's triple-zeta basis set are recommended on the basis of smallest root mean square error.     

Blue‐shifted and red‐shifted hydrogen bonds: Theoretical study of the CH3CHO· · ·HNO complexes

The blue‐shifted and red‐shifted H‐bonds have been studied in complexes CH₃CHO…HNO. At the MP2/6‐31G(d), MP2/6‐31+G(d,p) MP2/6‐311++G(d,p), B3LYP/6‐31G(d), B3LYP/6‐31+G(d,p) and B3LYP/6‐311++G(d,p) levels, the geometric structures and vibrational frequencies of complexes CH₃CHO…HNO are calculated by both standard and CP‐corrected methods, respectively. Complex A exhibits simultaneously red‐shifted C? H…O and blue‐shifted N? H…O H‐bonds. Complex B possesses simultaneously two blue‐shifted H‐bonds: C? H…O and N? H…O. From NBO analysis, it becomes evident that the red‐shifted C? H…O H‐bond can be explained on the basis of the two opposite effects: hyperconjugation and rehybridization. The blue‐shifted C? H…O H‐bond is a result of conjunct C? H bond strengthening effects of the hyperconjugation and the rehybridization due to existence of the significant electron density redistribution effect. For the blue‐shifted N? H…O H‐bonds, the hyperconjugation is inhibited due to existence of the electron density redistribution effect. The large blue shift of the N? H stretching frequency is observed because the rehybridization dominates the hyperconjugation. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2006     

The use of the hybrid density functional theory (DFT) approach for calculation of vibrational spectra: nitromethane

The molecular geometry of nitromethane was optimized and its force field and vibrational spectrum were calculated by the BECKE3LYP method. The accuracy of optimization of the geometry of MeNO₂ obtained by this method using the 6–311G(d,p) and 6–311++G(d,p) basis sets is not poorer than that obtained at the second-order Møller-Plesset level of perturbation theory (MP2). The vibrational frequencies of nitromethane and its d₁, d₂, and d₃ isotopomers obtained by the BECKE3LYP method are in much better agreement with the experimental data than those calculated at the MP2 level using the same basis set. The average absolute error of calculations performed without the use of any scaling factors is ～2% for frequencies; the maximum deviation is ～4%.     

Experimental and theoretical investigation of the molecular and electronic structure of N′‐benzylidene‐N‐[4‐(3‐methyl‐3‐phenyl‐cyclobutyl)‐thiazol‐2‐yl]‐chloro‐acetic acid hydrazide

The title compound, N′‐benzylidene‐N‐[4‐(3‐methyl‐3‐phenyl‐cyclobutyl)‐thiazol‐2‐yl]‐chloro‐acetic acid hydrazide, has been synthesized and characterized by elemental analysis, IR, ¹H and ¹³C NMR, and X‐ray single crystal diffraction. The compound crystallizes in the orthorhombic space group P 21 21 21 with a = 5.8671 (3) Å, b = 17.7182 (9) Å, and c = 20.6373 (8) Å. Moreover, the molecular geometry from X‐ray experiment, the molecular geometry, vibrational frequencies, and gauge‐including atomic orbital ¹H and ¹³C chemical shift values of the title compound in the ground state have been calculated by using the Hartree–Fock and density functional methods (B3LYP) with 6‐31G(d) and 6‐31G(d,p) basis sets. The results of the optimized molecular structure are exhibited and compared with the experimental X‐ray diffraction. Besides, molecular electrostatic potential, Frontier molecular orbitals, and thermodynamic properties of the title compound were determined at B3LYP/6‐31G(d) levels of theory. © 2011 Wiley Periodicals, Inc. Int J Quantum Chem, 2012     

FT-IR and FT-Raman, vibrational assignments, molecular geometry, ab initio (HF) and DFT (B3LYP) calculations for 1,3-dichlorobenzene

The FT-IR and FT-Raman vibrational spectra of 1,3-dichlorobenzene (1,3-DCB) have been recorded using Bruker IFS 66 V Spectrometer in the range 4000-100 cm(-1). A detailed vibrational spectral analysis has been carried out and assignments of the observed fundamental bands have been proposed on the basis of peak positions and relative intensities. The optimized molecular geometry, vibrational frequencies, atomic charges, dipole moment, rotational constants and several thermodynamic parameters in the ground state were calculated using ab initio Hartree-Fock (HF) and DFT (B3LYP) methods with 6-31++G (d, p) and 6-311++G (d, p) basis sets. With the help of different scaling factors, the observed vibrational wave numbers in FT-IR and FT-Raman spectra were analyzed and assigned to different normal modes of the molecule. Most of the modes have wave numbers in the expected range. The inductive effect of Chlorine atoms in the benzene molecule has also been investigated.