Accurate Vibrational Spectra of Large Molecules by Density Functional Computations beyond the Harmonic Approximation: The Case of Azabenzenes

Reliable structural, thermodynamic, and spectroscopic properties of organic molecules adsorbed on silicon surfaces from computational modeling: the case of glycine@Si(100)

Chemisorption of glycine on Si(100) has been studied by an integrated computational strategy based on perturbative anharmonic computations employing geometries and harmonic force fields evaluated by hybrid density functionals coupled to purposely tailored basis sets. It is shown that such a strategy allows the prediction of spectroscopic properties of isolated and chemisorbed molecules with comparable accuracy, paving the route toward a detailed analysis of surface-induced changes of glycine vibrational spectra.

Environmental effects in computational spectroscopy: accuracy and interpretation

Spectroscopic techniques are valuable tools for understanding the structure and dynamics of complex systems, such as biomolecules or nanomaterials. Most of the current research is devoted to the development of new experimental techniques for improving the intrinsic resolution of different spectra. However, the subtle interplay of several different effects acting at different length and time scales still makes the interpretation and analysis of such spectra a very difficult task. In this respect, computational spectroscopy is becoming a needful and versatile tool for the assignment and interpretation of experimental spectra. It is in fact possible nowadays to model with relatively high accuracy the physical-chemical properties of complex molecules in different environments, and to link spectroscopic evidence directly to the structural and dynamical properties of optically or magnetically active solvated probes. In this Review, significant steps toward the simulation of entire spectra in condensed phases are presented together with some basic aspects of computational spectroscopy, which highlight how intramolecular and intermolecular degrees of freedom influence several spectroscopic parameters.

Harmonic and Anharmonic Vibrational Frequency Calculations with the Double-Hybrid B2PLYP Method: Analytic Second Derivatives and Benchmark Studies

This work aims to provide reliable benchmark data on the accuracy of harmonic and anharmonic vibrational frequencies computed with the B2PLYP double-hybrid density functional method. The exchange-correlation contributions required for the B2PLYP analytical second derivatives are presented here, which allow for the effective calculation of harmonic frequency as well as cubic and semidiagonal quartic force fields. The latter, in turn, are necessary to compute the anharmonic vibrational frequencies with the perturbative approach (VPT2). The quality of harmonic vibrational frequencies computed in conjunction with basis sets of double- to quadruple-ζ quality has been checked against reference data from the F38 benchmark set. Then, for an additional set of small closed- and open-shell systems, both harmonic frequencies and anharmonic contributions computed at the B2PLYP/N07D and the B2PLYP/aug-cc-pVTZ levels have been compared to their CCSD(T) counterparts. Moreover, for selected medium-size molecules (furan, pyrrole, thiophene, uracil, anisole, phenol, and pyridine), anharmonic frequencies have been compared to well established experimental results. Such benchmark studies have shown that the B2PLYP/N07D model provides good quality harmonic frequencies and describes correctly anharmonic contributions, the latter being of similar accuracy to their B3LYP/N07D counterparts, but obtained at significantly larger computational cost. Additionally, increased accuracy can be obtained by adopting hybrid models where the B2PLYP/N07D anharmonic contributions are combined with harmonic frequencies computed with more accurate quantum mechanical (QM) approaches or by B2PLYP with larger basis sets. This work confirmed also that most of the recently developed density functionals are significantly less suited for vibrational computations, while the B2PLYP method can be recommended for spectroscopic studies where a good accuracy of vibrational properties is required.

The anharmonic force field of 1,3-cyclopentadienes

The cubic and quartic semi-diagonal anharmonic force fields of 1,3-cyclopentadienes, C(5)H(6) and C(5)D(6), of the 1,3-cyclopentadiene-5-d and 1,3-cyclopentadiene-1,2,3,4,5-d(5) have been obtained using density functional theory, DFT, with the B97-1 functional and a triple-zeta plus double polarization (TZ2P) basis set. The 108 fundamental frequencies computed by second order vibrational perturbation theory are in very good agreement with the experimental data, with a root mean square deviation (RMSD) equal to 7 cm(-1). Some fundamental frequencies have been reassigned and two have been assigned for the first time on the basis of the present results. The RMSD value, identical to that obtained for naphthalene (E. Cané, A. Miani and A. Trombetti, J. Phys. Chem. A, 2007, 111, 8218), indicates that computations based on this method can be confidently extended to larger polycyclic aromatic hydrocarbons.

Anharmonic Force Fields of Naphthalene-h8and Naphthalene-d8

The cubic and the quartic semidiagonal anharmonic force fields of naphthalene-h8 and -d8 are obtained using density functional theory (DFT) with the B9-71 functional and a triple-zeta plus double polarization (TZ2P) basis set. The fundamental frequencies computed by second-order vibrational perturbation theory are in very good agreement with the experimental data, with a mean absolute deviation (MAD) of 4 cm(-1) for C(10)H(8) and 6 cm(-1) for C(10)D(8). Some of the fundamental frequencies have been reassigned on the basis of the present results. Only CH stretchings seem to be significantly affected by Fermi resonances, with two shifts larger than 10 cm(-1). Calculated infrared harmonic intensities reproduce the experimental data within 15%, with the exception of CH stretchings affected by a larger error. Scale factors from C(10)H(8) have been tested by deriving the fundamental frequencies of C(10)D(8) from the theoretical harmonic ones. These fundamentals are in nice agreement with those obtained from the C(10)D(8) anharmonic force field. These results support the use of scale factors to calculate the vibration spectra of larger polycyclic aromatic hydrocarbons of great astrophysical interest.

Anharmonicity of amide modes

The principal contributions to the anharmonic coupling of amide vibrations are explored with the objective of comparing recent experiments with density functional theory and evaluating simple models of mode coupling. Experimental information obtained by means of two-dimensional infrared spectroscopy (2D IR) is reasonably well predicted by the computed one- and two-quantum anharmonic modes of amide-A, -I, and -II types in mono-, di- and tripeptides. The expansion of the vibrational energy up to the cubic and quartic coupling of harmonic modes suggested criteria to assess how localized are the forces determining the anharmonicity. The off-diagonal anharmonicity between an amide-A and one other amide mode was shown to be mainly determined by forces involving only these two modes, whereas the off-diagonal anharmonicity of two amide-I modes in peptides depended significantly on forces due to motions other than those of the amide-I type. Both the diagonal and off-diagonal anharmonicities exhibit sensitivity to peptide structures. These results should prove useful in linking 2D IR experimental results to secondary structure. Further, the results are used to evaluate the vibrational exciton model for the mixed-mode anharmonicities of the amide-I transitions.

Anharmonic analysis of the vibrational spectrum of ketene by density functional theory using second-order perturbative approach

The paper reports main results of a comprehensive study of the vibrational spectrum of ketene computed using second-order perturbation theory treatment based on quartic, cubic and semidiagonal quartic force constants. Two different models--a homogeneous model using the same density functionals and basis functions for the harmonic calculations and anharmonic corrections, and a hybrid model in which the two parts of the calculation are conducted using different density functionals and basis sets--have been employed in the present calculations. Different DFT and CCSD methods and DZ and TZ extended basis sets involving diffuse and polarization functions have been used to calculate optimized and vibrationally averaged geometrical parameters, the harmonic and anharmonic vibrational frequencies and the spectroscopic constants such as anharmonicity constants, rotational constants, rotation-vibration coupling constants, Nielsen's centrifugal distortion constants and Coriolis coupling constants. Homogeneous model is found to be superior to the hybrid model in several respects. Difficulties in the hybrid model may arise due to one of the following reasons: (a) the basic requirement that the geometry optimization and frequency calculations must be done at the same level of theory to have valid frequencies is not met in the hybrid model; (b) the molecular structure gets reoptimized at the low level for anharmonic corrections; (c) in addition, the perturbation could also diverge for the above reasons, particularly for the very low, very anharmonic terms where the harmonic approximation is not close enough to make the perturbation work.

High Resolution Rovibrational Spectroscopy of Chiral and Aromatic Compounds

The analysis of selected rovibrationally resolved infrared spectra of some relatively heavy and large polyatomic molecules is reviewed. A short historical summary of the development of high resolution interferometric Fourier transform infrared (FTIR) spectrometers is given and the possibilities of the currently most highly resolving FTIR spectrometer, which is commercially available in the Bruker IFS 125 series, are discussed. The computational tools necessary to analyse FTIR spectra are described briefly. As examples of rovibrational analysis the spectra of three selected molecules CHCl(2)F, CDBrClF, and pyridine (C(5)H(5)N) are discussed. The spectrum of CHCl(2)F, a fluorochlorohydrocarbon, is of interest for a better understanding of the chemistry of the Earth's atmosphere. CDBrClF is a chiral molecule and therefore the analysis of its rovibrational spectra provides the basis for carrying out further experiments towards the detection of molecular parity violation. The analysis of the pyridine FTIR spectra illustrates the potential of the new generation of FTIR spectrometers in the study of spectra and rovibrational dynamics of aromatic systems and molecules of potential biological interest.

Structure and Dynamics of the Aniline−Argon Complex as Derived from its Potential Energy Surface

The structure and dynamics of the van der Waals (vdW) complex of aniline (An) with argon (Ar) are studied using ab initio methods. The inversion potential of the aniline-argon (AnAr) complex perturbed by the weak vdW interaction is calculated taking into account subtle corrections from the zero-point energy of the vdW modes and from the frequency shifts of the An normal modes modified by the complexation. The intermolecular potential energy surface (PES) of the AnAr complex is determined by performing a large-scale computation of the interaction energy and the fitting of the analytical many-body expansion to the set of single-point interaction energies. The PES determined shows two deep local minima corresponding to the anti and syn AnAr conformers. The difference in the energies of these two minima is only 15 cm-1, but it is sufficient to localize the inversion wave functions and to form the two conformers. In the conformer anti (syn) of lower (higher) energy, Ar is bound to the An ring opposite (adjacent) the amino-hydrogens. In the additional local minima higher in energy, Ar approaches the aniline ring between the C-H bonds near its plane. An additional local minimum is located opposite of nitrogen between the two N-H bonds. The high-energy minima are, however, too flat to form stable conformers. The perturbation of the interaction of Ar with the phenyl ring by the NH2 group is described by the vdW hole, which is responsible for unusually strong intermode mixing in the excited intermolecular vibrational states. The analysis of these states calculated for the ground (S0) as well as the first excited electronic state (S1) resolves difficulties faced earlier with the assignment of the observed vibronic bands of AnAr.

Specific heat of molecular crystals from atomic mean square displacements with the Einstein, Debye, and Nernst-Lindemann models

Analysis of atomic displacement parameters (ADPs) from multitemperature diffraction data provides mean-field molecular translation and libration frequencies. These quantities have been combined with molecular deformation frequencies calculated ab initio, e.g. by DFT methods, to calculate the specific heat Cv of molecular crystals of naphthalene, anthracene, and hexamethylenetetramine. If the difference Cp - Cv is represented by the Nernst-Lindemann relation, Cp curves from diffraction experiments and ab initio calculations agree well with those based on calorimetry. Agreement is better if the Debye rather than the Einstein model is chosen to represent the contribution of the translational vibrations. Compressibilities estimated from the differences Cp - Cv are 2-5 times higher than those obtained from compressibility measurements at 298 K and Grüneisen constants derived from the temperature dependence of ADPs.

An Experimental Raman and Theoretical DFT Study on the Self-Association of Acrylonitrile

The liquid structure of acrylonitrile (propenenitrile) has been investigated using Raman spectroscopy and density functional theory (DFT) ab initio calculations with the 6-311++G** basis set. Two different and complementary experimental approaches were undertaken: FT-Raman spectra of 13 acrylonitrile solutions in carbon tetrachloride (concentration range=0.25-12.0 mol.L-1) were studied in detail including principal component analysis (PCA) of the CN stretching band. Furthermore, dispersive Raman spectra of neat acrylonitrile were obtained at eight different temperatures from 238 up to 343 K. The complex and asymmetric acrylonitrile Raman CN stretching band can be decomposed into two components attributed to monomeric and self-associated forms. Ab initio results fully support this assignment and suggest that the self-associated complex is a nonplanar trimer held together by dipole-dipole interactions. At ambient temperature, the composition of acrylonitrile can be expressed as a mixture of 25% monomers and 75% trimers. Close to the boiling point, trimers still represent 65% of the liquid composition. The corresponding enthalpy of association was estimated to be -22+/-2 kJ.mol-1.

Anharmonic vibrational properties by a fully automated second-order perturbative approach

This paper describes the implementation of a fully automated code for the building of anharmonic force constants and their use in a second-order perturbative evaluation of vibrorotational parameters. Next, a number of test applications are discussed, which show the strengths and limits of various computational levels.

Comparison of the numerical stability of methods for anharmonic calculations of vibrational molecular energies

On model examples, we compare the performance of the vibrational self-consistent field, variational, and four perturbational schemes used for computations of vibrational energies of semi-rigid molecules, with emphasis on the numerical stability. Although the accuracy of the energies is primarily dependent on the quality of the potential energy surface, approximate approaches to the anharmonic vibrational problem often do not converge to the same results due to the approximations involved. For furan, the sensitivity to variations of the anharmonic potential was systematically investigated by adding random noise to the cubic and quartic constants. The self-consistent field methods proved to be the most resistant to the potential variations. The second order perturbational techniques are sensitive to random degeneracies and provided the least stable results. However, their stability could be significantly improved by a simple generalization of the perturbational formula. The variational configuration interaction is practically limited by the size of the matrix that can be diagonalized for larger molecules; however, relatively fewer states need to be involved than for smaller ones, in favor of the computing.

Vibrational structure theory: new vibrational wave function methods for calculation of anharmonic vibrational energies and vibrational contributions to molecular properties

A number of recently developed theoretical methods for the calculation of vibrational energies and wave functions are reviewed. Methods for constructing the appropriate quantum mechanical Hamilton operator are briefly described before reviewing a particular branch of theoretical methods for solving the nuclear Schrödinger equation. The main focus is on wave function methods using the vibrational self-consistent field (VSCF) as starting point, and includes vibrational configuration interaction (VCI), vibrational Møller-Plesset (VMP) theory, and vibrational coupled cluster (VCC) theory. The convergence of the different methods towards the full vibrational configuration interaction (FVCI) result is discussed. Finally, newly developed vibrational response methods for calculation of vibrational contributions to properties, energies, and transition probabilities are discussed.

2-pyridone: The role of out-of-plane vibrations on the S1↔S0 spectra and S1 state reactivity

The S(1)<-->S(0) vibronic spectra of supersonic jet-cooled 2-pyridone [pyridin-2-one (2PY)] and its N-H deuterated isotopomer (d-2PY) have been recorded by two-color resonant two-photon ionization, laser-induced fluorescence and emission, and fluorescence depletion spectroscopies. By combining these methods, the B origin of 2PY at 0(0) (0)+98 cm(-1) and the bands at +218 and +252 cm(-1) are identified as overtones of the S(1) state out-of-plane vibrations nu(1) (') and nu(2) ('), as are the analogous bands of d-2PY. Anharmonic double-minimum potentials are derived for the respective out-of-plane coordinates that predict further nu(1) (') and nu(2) (') overtones and combinations, reproducing approximately 80% of the vibronic bands up to 600 cm(-1) above the 0(0) (0) band. The fluorescence spectra excited at the electronic origins and the nu(1) (') and nu(2) (') out-of-plane overtone levels confirm these assignments. The S(1) nonplanar minima and S(1)<--S(0) out-of-plane progressions are in agreement with the determination of nonplanar vibrationally averaged geometries for the 0(0) (0) and 0(0) (0)+98 cm(-1) upper states by Held et al. [J. Chem. Phys. 95, 8732 (1991)]. The fluorescence lifetimes of the S(1) state vibrations show strong mode dependence: Those of the out-of-plane levels decrease rapidly above 200 cm(-1) excess vibrational energy, while the in-plane vibrations nu(5) ('), nu(8) ('), and nu(9) (') have longer lifetimes, although they are above or interspersed with the "dark" out-of-plane states. This is interpreted in terms of an S(1) (') state reaction with a low barrier towards a conical intersection with a prefulvenic geometry. Out-of-plane vibrational states can directly surmount this barrier, whereas in-plane vibrations are much less efficient in this respect. Analysis of the fluorescence spectra allows to identify nine in-plane S(0) (') state fundamentals, overtones of the S(0) state nu(1) (") and nu(2) (") out-of-plane vibrations, and >30 other overtones and combination bands. The B3LYP6-311++G(d,p) calculated anharmonic wave numbers are in very good agreement with the observed fundamentals, overtones, and combinations, with a deviation Delta(rms)=1.3%.

Combined Raman Spectroscopic and Theoretical Investigation of Fundamental Vibrational Bands of Furfuryl Alcohol (2-furanmethanol)

Furfuryl alcohol (FA) is a promising reactive precursor for new materials. FA reaction mechanisms, that is, self-reactions or cross reactions with other substances, can be studied by vibrational spectroscopy. We present a necessary prerequisite for such studies by a Raman spectroscopic and theoretical study of FA in weakly interacting environments. It is the first study of FA vibrational properties based on density functional theory (DFT/B3LYP), and a recently proposed hybrid approach to the calculation of fundamental frequencies, which also includes an anharmonic contribution. FA occupies five different conformational states, each with more than 5% probability, and two of these dominate at T = 298 K. Excluding one frequency, the remaining ones are predicted as a weighted average over the two dominant conformers to a best RMS error of 8 cm(-1) and are qualitatively assigned. The excluded CH stretching mode is underestimated by 65 cm(-1). This may be due to a combination of an insufficient level of theory and the neglect of Fermi interactions for properly describing this type of mode.

Prediction of accurate anharmonic experimental vibrational frequencies for water clusters, (H2O)n, n=2-5

Accurate anharmonic experimental vibrational frequencies for water clusters consisting of 2-5 water molecules have been predicted on the basis of comparing different methods with MP2/aug-cc-pVTZ calculated and experimental anharmonic frequencies. The combination of using HF/6-31G* scaled frequencies for intramolecular modes and anharmonic frequencies for intermolecular modes gives excellent agreement with experiment for the water dimer and trimer and are as good as the expensive anharmonic MP2 calculations. The water trimer, the cyclic Ci and S4 tetramers, and the cyclic pentamer all have unique peaks in the infrared spectrum between 500 and 800 cm-1 and between 3400 and 3700 cm-1. Under the right experimental conditions these different clusters can be uniquely identified using high-resolution IR spectroscopy.

Toward an integrated computational approach to CW-ESR spectra of free radicals

Interpretation of structural properties and dynamic behaviour of molecules in solution is of fundamental importance to understand their stability, chemical reactivity and catalytic action. Information can be gained, in principle, by a variety of spectroscopic techniques, magnetic as well as optical. In particular, continuous wave electron spin resonance (CW-ESR) measurements are highly informative. However, the wealth of structural and dynamic information which can be extracted from ESR spectroscopy is, at present, limited by the necessity of employing computationally efficient models, which are increasingly complex as they need to take into account diverse relaxation processes affecting the spectrum. In this paper, we address the basic theoretical tools needed to predict, essentially ab initio, CW-ESR spectra observables according to the stochastic Liouville equation (SLE) approach, combined with quantum mechanical and hybrid methods for the accurate and efficient computation of structural, spectroscopic and magnetic properties of molecular systems. We shall discuss, on one hand, the quantum mechanical calculation of magnetic observables, via density functional theory (DFT), time-dependent DFT (TD-DFT) and application of the polarizable continuum model (PCM) for the description of environmental effects, including anisotropic environments and systems where different regions are characterized by different dielectric constants. One the other hand, the explicit evaluation of dynamical effects will be discussed based on the numerically exact treatment of the SLE in the presence of several relaxation processes, which has been proven to be a challenging task.

Harmonic and anharmonic contributions to parity-violating vibrational frequency difference between enantiomers of chiral molecules

A general perturbative procedure for the computation of harmonic and anharmonic contributions to parity-violating vibrational shifts is introduced and applied to PHBrF and AsHBrF. The results point out the importance of both diagonal and off-diagonal anharmonic contributions and indicate that some parity-violating shift of AsHBrF approaches the resolution forecasted for next generation experiments. The proposed approach is sufficiently general and computationally effective to allow studies of similar and larger molecular systems.

Obtaining the structure and bond rotational potential of a substituted ethane by NMR spectroscopy of solutions in nematic liquid-crystalline solvents

Partially averaged dipolar couplings (also referred to as residual dipolar couplings) D(ij) can be obtained from the analysis of the NMR spectra of molecules dissolved in liquid-crystalline solvents. Their values for a nonrigid molecule depend upon the bond lengths and angles, the rotational potentials, and the orientational order of the molecules. The molecule studied, 1-chloro-2-bromoethane, is one of the simplest example of a substituted alkane in which the rotational potential has three minimum-energy positions, trans and gauche+/-conformations, and the present investigation explores the problems inherent in deriving the form of the potential and the molecular geometry from the set of partially averaged couplings between the protons, and between protons and (13)C nuclei. The geometrical parameters and the rotational potential obtained are compared with the results from a density-functional theory method.

Conformational Analysis of 2,2‘-Bithiophene: A 1H Liquid Crystal NMR Study Using the 13C Satellite Spectra

We have obtained a very large data set of spectral parameters from the analysis of (1)H NMR and (13)C satellite spectra of 2,2'-bithiophene dissolved in anisotropic, partially orienting mesophases. In particular, this parameter set includes 33 dipolar couplings, which are directly related to the interatomic distances, the dihedral angle phi between the two thiophenic rings, and the anisotropic solute-solvent interaction potential. This allows an exhaustive investigation of the conformational equilibrium of 2,2'-bithiophene in a liquidlike phase. Comparison with the predictions of high-level theoretical calculations for the isolated molecule provides evidence of a strong flattening as well as the sharpening effect of the medium on the conformer population. The approximations needed to apply vibrational corrections to flexible molecules are discussed in detail and some general conclusions concerning their effect on structure and conformational equilibria are proposed.

The Structure and Conformations of 2-Thiophenecarboxaldehyde Obtained from Partially Averaged Dipolar Couplings

The proton NMR spectra of samples of 2-thiophenecarboxaldehyde dissolved in a nematic liquid crystalline solvent, including those from all five singly labelled 13C isotopomers, have been obtained. These have been analysed to yield sets of partially averaged dipolar couplings which have been used to determine the structure and the relative amounts of the cis and trans forms, which are the two minimum-energy structures generated by rotation about the ring-aldehyde bond. A procedure for applying vibrational corrections to the dipolar couplings in the presence of large amplitude motions is discussed.

Vibrational Spectroscopy of the G· · ·C Base Pair: Experiment, Harmonic and Anharmonic Calculations, and the Nature of the Anharmonic Couplings

The results of harmonic and anharmonic frequency calculations on a guanine-cytosine complex with an enolic structure (a tautomeric form with cytosine in the enol form and with a hydrogen at the 7-position on guanine) are presented and compared to gas-phase IR-UV double resonance spectral data. Harmonic frequencies were obtained at the RI-MP2/cc-pVDZ, RI-MP2/TZVPP, and semiempirical PM3 levels of electronic structure theory. Anharmonic frequencies were obtained by the CC-VSCF method with improved PM3 potential surfaces; the improved PM3 potential surfaces are obtained from standard PM3 theory by coordinate scaling such that the improved PM3 harmonic frequencies are the same as those computed at the RI-MP2/cc-pVDZ level. Comparison of the data with experimental results indicates that the average absolute percentage deviation for the methods is 2.6% for harmonic RI-MP2/cc-pVDZ (3.0% with the inclusion of a 0.956 scaling factor that compensates for anharmonicity), 2.5% for harmonic RI-MP2/TZVPP (2.9% with a 0.956 anharmonicity factor included), and 2.3% for adapted PM3 CC-VSCF; the empirical scaling factor for the ab initio harmonic calculations improves the stretching frequencies but decreases the accuracy of the other mode frequencies. The agreement with experiment supports the adequacy of the improved PM3 potentials for describing the anharmonic force field of the G...C base pair in the spectroscopically probed region. These results may be useful for the prediction of the pathways of vibrational energy flow upon excitation of this system. The anharmonic calculations indicate that anharmonicity along single mode coordinates can be significant for simple stretching modes. For several other cases, coupling between different vibrational modes provides the main contribution to anharmonicity. Examples of strongly anharmonically coupled modes are the symmetric stretch and group torsion of the hydrogen-bonded NH2 group on guanine, the OH stretch and torsion of the enol group on cytosine, and the NH stretch and NH out-of-plane bend of the non-hydrogen-bonded NH group on guanine.

Theoretical Study of Pyridine and 4,4‘-Bipyridine Adsorption on the Lewis Acid Sites of Alumina Surfaces Based on Ab Initio and Density Functional Cluster Calculations

The interactions of pyridine and 4,4'-bipyridine with the Lewis acid sites of alumina surfaces are investigated using ab initio and density functional calculations. Four cluster models of different sizes and shapes are chosen to represent the Lewis acid sites: three hydrogenated clusters Al(OH)(3), Al(4)O(9)H(6), and Al(10)O(21)H(12) and one non-hydrogenated cluster Al(4)O(6). The Hartree-Fock (HF) and B3LYP approaches with two basis sets 6-31G and 6-31+G are used to calculate the geometries, the electronic structures, the vibrational frequencies, and the adsorption energies of the complexes formed upon interaction of pyridine or 4,4'-bipyridine ligands on the cluster surfaces. Electronic structures are determined by the electrostatic potential (ESP) analysis of charges. Adsorption energies are calculated with corrections made for zero-point energies (ZPE) and basis set superposition error (BSSE). The ESP analysis of atomic charges reveals that the charge-transfer effects are more important in Lewis complexes formed with Al(4)O(6) cluster than in those formed with hydrogenated clusters Al(OH)(3), Al(4)O(9)H(6), and Al(10)O(21)H(12). The significantly larger charge transferred from pyridine or 4,4'-bipyridine ligand to Al(4)O(6) cluster should increase the adsorption energy of these complexes. Consequently, at all levels of calculation, the adsorption energies of pyridine and 4,4'-bipyridine complexed to Al(4)O(6) cluster ( approximately 46 kcal/mol), which compare very well to experiment, are strongly larger than those obtained for both pyridine and 4,4'-bipyridine ligands complexed to Al(OH)(3) (32 kcal/mol), Al(4)O(9)H(6) (24 kcal/mol) and Al(10)O(21)H(12) (25 kcal/mol) clusters. The corrected adsorption energy is found to be insensitive to basis set and electron correlation effects. It essentially depends on the ionic character of the cluster model rather than on its size. For 4,4'-bipyridine complexes, similar results to those obtained for pyridine are found, and the geometry and the amount of charge of the unbound pyridyl ring are unchanged upon complexation. The calculated vibrational frequencies and frequency shifts are little sensitive to the size and shape of the cluster model. The two ring stretching modes 8a and 19b of pyridine and 4,4'-bipyridine observed in the 1400-1600 cm(-1) region are the most affected modes upon adsorption, in good agreement with the available infrared and Raman data.

Accurate ab initio predictions of ionization energies of hydrocarbon radicals: CH2, CH3, C2H, C2H3, C2H5, C3H3, and C3H5

The ionization energies for methylene (CH2), methyl (CH3), ethynyl (C2H), vinyl (C2H3), ethyl (C2H5), propargyl (C3H3), and allyl (C3H5) radicals have been calculated by the wave-function-based ab initio CCSD(T)/CBS approach, which involves the approximation to the complete basis set (CBS) limit at the coupled-cluster level with single and double excitations plus a quasiperturbative triple excitation [CCSD(T)]. When it is appropriate, the zero-point vibrational energy correction, the core-valence electronic correction, the scalar relativistic effect correction, the diagonal Born-Oppenheimer correction, and the high-order correlation correction have also been made in these calculations. The comparison between the computed ionization energy (IE) values and the highly precise experimental IE values determined in previous pulsed field ionization-photoelectron (PFI-PE) studies indicates that the CCSD(T)/CBS method is capable of providing accurate IE predictions for these hydrocarbon radicals achieving error limits well within +/-10 meV. The benchmarking of the CCSD(T)/CBS IE predictions by the PFI-PE experimental results also lends strong support for the conclusion that the CCSD(T)/CBS approach with high-level energy corrections can serve as a valuable alternative for reliable IE determination of radicals, particularly for those radicals with very unfavorable Franck-Condon factors for photoionization transitions near their ionization thresholds.

Accurate vibrational spectra and magnetic properties of organic free radicals: The case of H2CN

We present the structural, magnetic, and vibrational properties of H2CN computed using a second-order perturbative approach in which equilibrium values and harmonic frequencies evaluated at the coupled-cluster level are combined with anharmonic and vibrational averaging contributions obtained by hybrid Hartree-Fock/Kohn-Sham methods. Our computations lead to remarkably accurate results and suggest some revision of the experimental vibrational assignments.

Vibrational computations beyond the harmonic approximation: Performances of the B3LYP density functional for semirigid molecules

The performances of the B3LYP density functional in the computation of harmonic and anharmonic frequencies were tested using 14 standard basis sets of double and triple zeta quality for a set of semirigid molecules containing from 4 to 12 atoms. The quality of the results is assessed by comparison with the most reliable computations available in the literature. The study reveals that the relatively cheap 6-31+G(d,p) basis set performs a very good job for harmonic frequency calculations and that B3LYP anharmonicities are in close agreement with the reference values irrespective of the basis set used. On these grounds "hybrid force fields" are proposed to achieve the best compromise between computer time and quality of the results.