Performance of density functional theory in computing nonresonant vibrational (hyper)polarizabilities

A set of exchange‐correlation functionals, including BLYP, PBE0, B3LYP, BHandHLYP, CAM‐B3LYP, LC‐BLYP, and HSE, has been used to determine static and dynamic nonresonant (nuclear relaxation) vibrational (hyper)polarizabilities for a series of all‐trans polymethineimine (PMI) oligomers containing up to eight monomer units. These functionals are assessed against reference values obtained using the Møller–Plesset second‐order perturbation theory (MP2) and CCSD methods. For the smallest oligomer, CCSD(T) calculations confirm the choice of MP2 and CCSD as appropriate for assessing the density functionals. By and large, CAM‐B3LYP is the most successful, because it is best for the nuclear relaxation contribution to the static linear polarizability, intensity‐dependent refractive index second hyperpolarizability, static second hyperpolarizability, and is close to the best for the electro‐optical Pockels effect first hyperpolarizability. However, none of the functionals perform satisfactorily for all the vibrational (hyper)polarizabilities studied. In fact, in the case of electric field‐induced second harmonic generation all of them, as well as the Hartree–Fock approximation, yield the wrong sign. We have also found that the Pople 6–31+G(d) basis set is unreliable for computing nuclear relaxation (hyper)polarizabilities of PMI oligomers due to the spurious prediction of a nonplanar equilibrium geometry. © 2013 Wiley Periodicals, Inc.     

Conformations and rotational barriers of 2,2′‐bithiazole and 4,4′‐dimethyl‐2,2′‐bithiazole semiemperical,ab initio,and density functional theory calculations

Conformational properties of 2,2′‐bithiazole and 4,4′‐dimethyl‐2,2′‐ bithiazole have been studied by using AM1 and PM3 semiemperical methods and ab initio HF/6‐311+G* and B3LYP/6‐311+G* calculations. All methods agree that the planar s‐trans conformation is the global minimum and the perpendicular conformation is the transition state. Additional local minima were found using the Hartree–Fock (HF) and B3LYP levels for 2,2′‐bithiazole while for 4,4′‐dimethyl derivative the minima was located only at the MP2//B3LYP level. The barrier heights for rotation are 1.72, 7.69, and 7.88 kcal/mol at the PM3, HF, and B3LYP levels, respectively, and methyl substitution did not affect appreciably this value. Fourier expansion terms and bond orders were used to explain the origins of the rotational barrier in terms of π conjugation, electrostatic interaction, and steric effects, which represent the main factors in the shape of the rotational barrier. © 2000 John Wiley & Sons, Inc. Int J Quant Chem 79: 367–377, 2000     

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     

Comparative studies on the structures,infrared spectrum,and thermodynamic properties of phthalocyanine using ab initio Hartree–Fock and density functional theory methods

The molecular structure, vibrational spectrum, standard thermodynamic functions, and enthalpy of formation of free base phthalocyanine (Pc) have been studied using the density functional theory B3LYP procedure, as well as the ab initio Hartree–Fock method. Various basis sets 3‐21G, 6‐31G*, and LANL2DZ have been employed. The results obtained at various levels are discussed and compared with each other and with the available experimental data. It is shown that calculations performed at the Hartree–Fock level cannot produce a reliable geometry and related properties such as the dipole moment of Pc and similar porphyrin‐based systems. Electron correlation must be included in the calculations. The basis set has comparatively less effect on the calculated results. The results derived at the B3LYP level using the smaller 3‐21G and LANL2DZ basis sets are very close to those produced using the medium 6‐31G* basis set. The geometry of Pc obtained at the B3LYP level has D_2h symmetry and the diameter of the central macrocycle is about 4 Å. The enthalpy of formation of Pc in the gas phase has been predicted to be 1518.50 kJ/mol at the B3LYP/6‐311G(2d,2p)//B3LYP/6‐31G* level via an isodesmic reaction. © 2001 John Wiley & Sons, Inc. Int J Quantum Chem, 2001     

Assessing the performance of density functional theory for the dynamic polarizabilities of amino acids: Treatment of correlation and role of exact exchange

The theoretical determination of electric response properties of the biological systems is a field where the application of density functional theory (DFT) appears to be quite promising. In this work, the performance of 41 density functional methods is evaluated in predicting dynamic polarizabilities of an experimental benchmark set of 20 proteinogenic amino acids. The behavior of a large number of density functionals, including various types of the local spin density approximation (LSDA), generalized gradient approximation (GGA), meta‐GGA (m‐GGA), hybrid‐GGA (h‐GGA), hybrid meta‐GGA (hm‐GGA), and range‐separated hybrid‐GGA (rsh‐GGA), has been assessed for the purpose. Analyzing the results of our DFT benchmarking, we found that these computationally economical methods show very diverse predictive capability and a careful selection of DFT functionals is very important in the polarizability calculations. Considering the role of exchange, correlation, dispersion and long‐range corrections, it turned out that in the LSDA class, SVWN3 gives better results than SPL and SVWN5 toward the reference values. Of the GGA methods, OPBE outperforms all other functionals. The M06‐L is the best method of m‐GGA class. The B3LYP and TPSSh are the best functionals of h‐GGA and hm‐GGA lineages, respectively. Finally, CAM‐B3LYP is the best method of rsh‐GGA functionals that predicts the most accurate polarizability for amino acids by a large margin with respect to others. Overall, the best performing functionals turn out to be hm‐GGAs TPSSh, TPSS1KCIS, M05, tau‐HCTHhyb, and h‐GGA B3LYP. Hopefully, the results of this investigation might provide the useful guidance to propose a new exchange‐correlation functional for calculating the optical properties of biomolecular materials. © 2013 Wiley Periodicals, Inc.     

Evaluation of DFT methods for computing the interaction energies of homomolecular and heteromolecular dimers of monosubstituted benzene

We present density functional theory (DFT) interaction energies for the sandwich and T‐shaped conformers of substituted benzene dimers. The DFT functionals studied include TPSS, HCTH407, B3LYP, and X3LYP. We also include Hartree–Fock (HF) and second‐order Møller–Plesset perturbation theory calculations (MP2), as well as calculations using a new functional, P3LYP, which includes PBE and HF exchange and LYP correlation. Although DFT methods do not explicitly account for the dispersion interactions important in the benzene–dimer interactions, we find that our new method, P3LYP, as well as HCTH407 and TPSS, match MP2 and CCSD(T) calculations much better than the hybrid methods B3LYP and X3LYP methods do. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2006     

FTIR and FT Raman, molecular geometry, vibrational assignments, ab initio and density functional theory calculations for 1,5-methylnaphthalene

The FTIR and FT Raman vibrational spectra of 1,5-methylnaphthalene (1,5-MN) have been recorded using Brunker IFS 66 V Spectrometer in the range 3600-10 cm(-1) in the solid phase. 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, harmonic frequencies, electronic polarizability, atomic charges, dipole moment, rotational constants and several thermodynamic parameters in the ground state were calculated using ab initio Hartree Fock (HF) and density functional B3LYP methods (DFT) with 6-311++ G(d) basis set. With the help of different scaling factors, the observed vibrational wavenumbers in FTIR and FT Raman spectra were analyzed and assigned to different normal modes of the molecule. Most of the modes have wavenumbers in the expected range. The results of the calculations were applied to simulated infrared and Raman spectra of the title compound which showed excellent agreement with the observed spectra.     

All‐electron LCAO calculations of the LiF crystal phonon spectrum: Influence of the basis set,the exchange‐correlation functional,and the supercell size

For the first time the convergence of the phonon frequencies and dispersion curves in terms of the supercell size is studied in ab initio frozen phonon calculations on LiF crystal. Helmann–Feynman forces over atomic displacements are found in all‐electron calculations with the localized atomic functions (LCAO) basis using CRYSTAL06 program. The Parlinski–Li–Kawazoe method and FROPHO program are used to calculate the dynamical matrix and phonon frequencies of the supercells. For fcc lattice, it is demonstrated that use of the full supercell space group (including the supercell inner translations) enables to reduce essentially the number of the displacements under consideration. For Hartree–Fock (HF), PBE and hybrid PBE0, B3LYP, and B3PW exchange‐correlation functionals the atomic basis set optimization is performed. The supercells up to 216 atoms (3 × 3 × 3 conventional unit cells) are considered. The phonon frequencies using the supercells of different size and shape are compared. For the commensurate with supercell k ‐points the best agreement of the theoretical results with the experimental data is found for B3PW exchange‐correlation functional calculations with the optimized basis set. The phonon frequencies at the most non‐commensurate k ‐points converged for the supercell consisting of 4 × 4 × 4 primitive cells and ensures the accuracy 1–2% in the thermodynamic properties calculated (the Helmholtz free energy, entropy, and heat capacity at the room temperature). © 2009 Wiley Periodicals, Inc. J Comput Chem 2009     

Theoretical investigation of the (hyper)polarizabilities of pyrrole homologues C4H4XH (X = N, P, As, Sb, Bi). A coupled-cluster and density functional theory study

Geometries, inversion barriers, static and dynamic electronic and vibrational dipole polarizability (alpha), and first (beta) and second (gamma) hyperpolarizability of the pyrrole homologues C(4)H(4)XH (X = N, P, As, Sb, Bi) have been calculated by Hartree-Fock, M?ller-Plesset second-order perturbation theory, coupled-cluster theory accounting for singles, doubles, and noniterative triple excitations methods, as well as density functional theory using B3LYP and PBE1PBE functionals and Sadlej's Pol and 6-311G basis sets. Relativistic effects on the heavier homologues stibole and bismole have been taken into account within effective core potential approximation. The results show that the electronic (hyper)polarizabilities monotonically increase with the atomic number of the heteroatom, consistent with the decrease in the molecular hardness. Ring planarization reduces the carbon-carbon bond length alternation of the cis-butadienic unit, enhancing the electronic polarizability values (alpha(e)) by 4-12% and the (hyper)polarizability values (and gamma(e)) by 30-90%. Pure vibrational and zero-point vibrational average contributions to the (hyper)polarizabilities have been determined within the clamped nucleus approach. In the static limit, the pure vibrational hyperpolarizabilities have a major contribution. Anharmonic corrections dominate the pure vibrational hyperpolarizabilities of pyrrole, while they are less important for the heavier homologues. Static and dynamic electronic response properties of the pyrrole homologues are comparable to or larger than the corresponding properties of the furan and cyclopentadiene homologue series.     

Theoretical analysis of vibrational spectra and scaling‐factor of 2‐aryl‐1,3,4‐oxadiazole derivatives

In this study, the direct molecular structure implementations for calculating vibrational spectra and scaling factors, and infrared intensities at both the Hartree–Fock (HF) and density functional (B3LYP) levels of theory with 6‐31G(d), 6‐311G(d), 6‐31++G(d,p), and 6‐311++G(d,p) basis sets are presented. Also, vibrational frequencies have been investigated as dependence on the choice of method and basis set. The parameters of molecular geometry and vibrational frequencies values of 2‐aryl‐1,3,4‐oxadiazoles 5a–g in the ground state have been calculated. Theoretical determination of vibrational frequencies is quite useful both in understanding the relationship between the molecular structures and scaling factor. The data of 2‐aryl‐1,3,4‐oxadiazoles 5a–g display significant electronic properties provide the basis for future design of efficient materials having the oxadiazole core and theoretical IR studies. © 2008 Wiley Periodicals, Inc. Int J Quantum Chem, 2009     

Theoretical study on the molecular structures of toluene,para‐fluorotoluene,para‐chlorotoluene,and 4‐methylpyridine and their sixfold internal rotational barriers

Two kinds of sixfold internal rotational configurations of toluene, para‐fluorotoluene, para‐chlorotoluene, and 4‐methylpyridine were calculated using Hartree–Fock (HF), second‐order Møller–Plesset (MP2), and Beck's three parameter hybrid functional using the LYP correlation functional (B3LYP) theory methods with various high‐level basis sets. Structures and energies were compared for different configurations. Calculations indicate that the orthogonal configuration has a local minimum while the planar configuration is a transition structure. Furthermore, geometries of the orthogonal and the planar configurations are quite similar, except for a methyl CH bond. Sixfold internal rotational barriers were calculated from the energy difference of two different configurations. For the calculated results, HF methods underestimated the rotational barriers, but MP2 calculations overestimated them. However, the density functional theory (DFT) method is a reliable method since the calculated internal rotational barriers are similar to the experimental ones. © 2000 John Wiley & Sons, Inc. Int J Quant Chem 77: 772–778, 2000     

Compatibility of correlation‐consistent basis sets with a hybrid Hartree–Fock/density functional method

The present study examines the feasibility of combining the correlation‐consistent basis sets developed by Dunning and coworkers with the hybrid Hartree–Fock/density functional method B3LYP. Furthermore, extrapolation to the complete basis set (CBS) limit minimizes errors due to the presence of an incomplete basis set and can act as a rigorous test of the limitations of the B3LYP method. Equilibrium geometries, energies, and harmonic vibrational frequencies were determined for a series of well‐studied, yet computationally challenging, small inorganics and their respective ions. The results were then extrapolated to the CBS limit, where applicable, and compared to experiment. It was found that a union between the hybrid Hartree–Fock/density functional B3LYP method and Dunning's augmented correlation‐consistent basis sets gave results that were comparable to molecular orbital methods that explicitly account for electron correlation. Furthermore, the minimum basis set necessary to attain reasonable results for the systems studied was aug‐cc‐pVTZ. Upgrading to the aug‐cc‐pVQZ level and subsequent extrapolation to the CBS limit further improved the overall agreement with the experiment. ©1999 John Wiley & Sons, Inc. J Comput Chem 20: 207–216, 1999     

Substitution Effect in NICS Values in Tropylium Ion Derivatives: An Ab Initio Study

Hartree–Fock (HF) and hybrid density functional theory (B3LYP) calculations were performed on tropylium ion and 19 of its mono- and diheteroatomic derivatives. The aromaticity in this class of compounds is evaluated based on the nucleus independent chemical shift (NICS) values. The NICS values are calculated at the center of the rings NICS (0) and at 1Å above the molecular plane NICS (1). The geometry optimization and NICS calculations were carried out at the HF/6–311+G?? and at the B3LYP/6–311+G (2d, p) density functional level, respectively. These calculations in the effects of heteroatoms such as N, B, P, and Si are considered on aromaticity, molecular properties, NICS values, and structural parameters.     

Theoretical study of the vibrational frequencies of carbon disulfide

A benchmark comparison for different computational methods and basis sets has been presented. In this study, five computational methods (Hartree–Fock (HF), MP2, B3LYP, MPW1MP91, and PBE1PBE) along with 18 basis sets have been applied to optimize the geometry of carbon disulfide (CS₂), and further calculate the vibrational frequencies of the optimized geometries. The differences between the calculated frequencies and corresponding experimental data are used to evaluate the efficiency of each combination of computational method and basis set. The comparison of frequency difference indicates that B3LYP generally gives the best prediction of frequencies for CS₂, whereas the other two density functional theory (DFT) methods, i.e., MPW1PW91 and PBE1PBE, often give parallel results. Although MP2 predicts the frequencies with accuracy almost as good as those from DFT methods, in a particular case, HF calculation outperforms MP2 as well as MPW1PW91 and PBE1PBE for prediction of the frequency of asymmetrical stretching for CS₂. © 2013 Wiley Periodicals, Inc.     

Nonlinear optical properties of the M @ C60 endohedrals (M = Cs,Li, and Na)

The Hartree–Fock (HF) approximation and the Single‐excitation configuration interaction (single‐CI) method are used, with sum‐over‐state (SOS) approach, in order to determine the static (ω = 0) and frequency dependent linear and non linear optical (NLO) properties of the M @ C₆₀ endohedrals (M = Cs, Li, and Na). We discuss the effects of displacement, movement direction, and type of atom on the (hyper) polarizability of the M at C₆₀ endohedrals. A HF‐single‐CI model yields the (hyper) polarizability magnitudes and spectra, which are in agreement with experiment. Our results indicate that the movement of the atom (M) effect on the NLO spectra of M @ C₆₀ endohedrals is dramatic. These results may provide new means to design some new types of NLO materials. © 2009 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

Calculations of electric dipole polarizabilities of polyatomic molecules

Coupled and uncoupled Hartree–Fock theories are used to compute the electric dipole polarizability of water, ammonia, and methane with three different GTO basis sets. Bounds for the geometric approximation to uncoupled polarizabilities are also computed to examine the accuracy of calculated values. The results are compared with those obtained by a variational-perturbation method proposed by Rebane. The numerical tests provide some information on the correlation terms affecting total computed polarizabilities. The computed values are shown to be in fair or good agreement with experimental data for the largest basis sets. The reliability of Rebane's method with respect to coupled HF procedure is also discussed.     

An experimental and theoretical study of molecular structure and vibrational spectra of 2-chloronicotinic acid by density functional theory and ab initio Hartree–Fock calculations

In this work, the Fourier transform Raman and Fourier transform infrared spectra of 2-chloronicotinic acid (2-CNA) are recorded in the solid phase. The molecular geometry, vibrational frequencies, infrared intensities and Raman scattering activities of 2-CNA in ground state have been calculated by using ab initio Hartree–Fock (HF) and density functional (B3LYP and B3PW91) methods with 6-31G(d) and 6-311G(d) basis sets level. On the basis of the comparison between calculated and experimental results and the comparison with related molecule, assignments of fundamental vibrational modes are examined. The optimized geometric parameters (bond lengths and bond angles) obtained by using HF show the best agreement with the experimental values of 2-CNA. Comparison of the observed fundamental vibrational frequencies of 2-CNA and calculated results by density functional (B3LYP and B3PW91) and Hartree–Fock methods indicates that B3LYP is superior to the scaled Hartree–Fock and B3PW91 approach for molecular vibrational problems.     

Comparative study on formamide–water complex

The hydrogen bonding of 1:1 complexes formed between formamide and water molecule have been investigated systematically using Hartree–Fock (HF), hybrid density functional theory (B3LYP), and post‐Hartree–Fock (MP2 and CCSD(T)) methods with range of basis sets 6‐31G(d), cc‐pVXZ (X = D, T, Q) and aug‐cc‐pVYZ (Y = D, T). Three stable structures are considered on the potential energy surface of formamide and water system. The optimized geometric parameters and interaction energies for various isomers at different levels are estimated. The IR frequencies, intensities, and frequency shifts are reported. This study shows that B3LYP/aug‐cc‐pVDZ method gives better performance for formamide‐water complexes. © 2009 Wiley Periodicals, Inc. Int J Quantum Chem, 2010.     

Vibrational spectra and potential energy distributions for 3-cyclopropenecarboxaldehyde by density functional and normal coordinate calculations

The structural stability and internal rotation in 3-cylopropenecarboxaldehyde were investigated by ab initio calculations with 6-311++G^** basis set. The calculations were carried out at the restricted Hartree–Fock (HF) and the Density Functional B3LYP levels. The vibrational frequencies were computed at HF and DFT-B3LYP levels. Normal coordinate calculations were carried out and potential energy distributions were calculated for the cis and the trans conformers of the molecule.