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

Fulvenallenyl cation (C7H5(+)) and its complex with an argon atom: results of high-level quantum-chemical calculations

The fulvenallenyl cation (C(7)H(5)(+)) and its complex with an argon atom have been studied by explicitly correlated coupled cluster theory at the CCSD(T)-F12x(x = a, b) level and by the double-hybrid density functional B2PLYP-D. For the free cation, an accurate equilibrium structure has been established and ground-state rotational constants of A(0) = 8116.4 MHz, B(0) = 2004.3 MHz, and C(0) = 1606.9 MHz are predicted. The equilibrium dipole moment is calculated to be μ(e) = 1.305 D, with the positive end of the dipole at the acetylenic hydrogen site. Anharmonic wavenumbers of C(7)H(5)(+) were obtained by combination of harmonic CCSD(T*)-F12a values and B2PLYP-D anharmonic contributions. The most intense vibration is the pseudoantisymmetric CC stretching vibration at 2083 cm(-1). The potential energy surface of the complex C(7)H(5)(+)·Ar is characterized by two energy minima of C(s) symmetry which are separated by a very low energy barrier. The dissociation energy of the most stable structure is predicted to be D(0) = 530 ± 30 cm(-1).

Density functional theory-based simulations of sum frequency generation spectra involving methyl stretching vibrations: effect of the molecular model on the deduced molecular orientation and comparison with an analytical approach

The knowledge of the first hyperpolarizability tensor elements of molecular groups is crucial for a quantitative interpretation of the sum frequency generation (SFG) activity of thin organic films at interfaces. Here, the SFG response of the terminal methyl group of a dodecanethiol (DDT) monolayer has been interpreted on the basis of calculations performed at the density functional theory (DFT) level of approximation. In particular, DFT calculations have been carried out on three classes of models for the aliphatic chains. The first class of models consists of aliphatic chains, containing from 3 to 12 carbon atoms, in which only one methyl group can freely vibrate, while the rest of the chain is frozen by a strong overweight of its C and H atoms. This enables us to localize the probed vibrational modes on the methyl group. In the second class, only one methyl group is frozen, while the entire remaining chain is allowed to vibrate. This enables us to analyse the influence of the aliphatic chain on the methyl stretching vibrations. Finally, the dodecanethiol (DDT) molecule is considered, for which the effects of two dielectrics, i.e. n-hexane and n-dodecane, are investigated. Moreover, DDT calculations are also carried out by using different exchange-correlation (XC) functionals in order to assess the DFT approximations. Using the DFT IR vectors and Raman tensors, the SFG spectrum of DDT has been simulated and the orientation of the methyl group has then been deduced and compared with that obtained using an analytical approach based on a bond additivity model. This analysis shows that when using DFT molecular properties, the predicted orientation of the terminal methyl group tends to converge as a function of the alkyl chain length and that the effects of the chain as well as of the dielectric environment are small. Instead, a more significant difference is observed when comparing the DFT-based results with those obtained from the analytical approach, thus indicating the importance of a quantum chemical description of the hyperpolarizability tensor elements of the methyl group.

Towards an accurate description of anharmonic infrared spectra in solution within the polarizable continuum model: reaction field, cavity field and nonequilibrium effects

We present a newly developed and implemented methodology to perturbatively evaluate anharmonic vibrational frequencies and infrared (IR) intensities of solvated systems described by means of the polarizable continuum model (PCM). The essential aspects of the theoretical model and of the implementation are described and some numerical tests are shown, with special emphasis towards the evaluation of IR intensities, for which the quality of the present method is compared to other methodologies widely used in the literature. Proper account of an incomplete solvation regime in the treatment of the molecular vibration is also considered, as well as inclusion of the coupling between the solvent and the probing field (cavity field effects). In order to assess the quality of our approach, comparison with experimental findings is reported for selected cases.

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.

Gas-phase and Ar-matrix SQM scaling factors for various DFT functionals with basis sets including polarization and diffuse functions

Scaling factors for Pulay's scaled quantum mechanical (SQM) scheme have been determined for four different widely used DFT functionals (PBE, B3LYP, B3PW91, and M06-2X) and for two basis sets (6-31++G** and aug-cc-pVTZ) by fitting computed results to 347 fundamental experimental vibrational frequencies of 33 molecules. Measurements in the gas phase and in solid argon matrices were used independently in the fitting procedure in order to provide a simple method of estimating matrix shifts. The accuracy of the new scaling factors is demonstrated on test molecules including hydrogen-bonded systems and molecules containing chlorine and sulfur atoms.

Time-dependent density functional theory investigation of the electronic spectra of hexanuclear chalcohalide rhenium(III) clusters

The UV-vis spectra of the hexanuclear chalcohalide rhenium(III) clusters of formula [Re(6)S(8)X(6)](4-) (X(-) = Cl(-), Br(-), I(-)) were investigated at the time-dependent density functional theory (TD-DFT) level employing B3LYP, PBE1PBE, and the double-hybrid B2PLYP functional in combination with the LANL2DZ basis set. We were able to reproduce the red shift experimentally observed when the halide changes from chloride to iodide. However, some discrepancies between experimental results and theory were found. First, we did not observe a remarkable dependence of the experimental ill-resolved bands on the solvent. Indeed, similar spectra were obtained taken CH(2)Cl(2) or CH(3)CN as solvents into account. Second, all calculations explained the origin of the band peaked at the low-energy region in contraposition with the one experimentally assumed by R. Long et al. (J. Am. Chem. Soc. 1996, 118, 4603), and theoretically made available by Arratia-Pérez et al. (J. Chem. Phys. 1999, 110, 2529). These authors have proposed a ligand-to-cluster charge transfer (LCCT) for all title complexes. Our findings undoubtedly allow such origin to be discarded. While the HGGA functionals lead to a cluster-to-halide ligand charge transfer (CLCT), an intracluster charge transfer (ICCT) has been considered by employing B2PLYP. This contribution showed B2PLYP in the presence of the solvent to be the best performer in studying the UV-vis spectra of large complexes of rhenium(III) containing the Re-S bond. We strongly recommended the use of the double-hybrid B2PLYP in studying UV-vis spectrum of rhenium complexes of size making the computational cost affordable. We expect that our work stimulates new experimental and theoretical investigations of the title complexes to confirm our assignment.