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

Time Resolved Raman Scattering of Molecules: A Quantum Mechanics Approach with Stochastic Schroedinger Equation

Raman scattering is a very powerful tool employed to characterize molecular systems. Here we propose a novel theoretical strategy to calculate the Raman cross-section in time domain, by computing the cumulative Raman signal emitted during the molecular evolution in time. Our model is based on a numerical propagation of the vibronic wave function under the effect of a light pulse of arbitrary shape. This approach can therefore tackle a variety of experimental setups. Both resonance and nonresonance Raman scattering can be retrieved, and also the time-dependent fluorescence emission is computed. The model has been applied to porphyrin considering both resonance and nonresonance conditions and varying the incident field duration. Moreover the effect of the vibrational relaxation, which should be taken into account when its time scale is similar to that of the Raman emission, has been included through the stochastic Schroedinger equation approach.

Static Electron Correlation in Anharmonic Molecular Vibrations: A Hybrid TAO-DFT Study

Hybrid thermally-assisted-occupation density functional theory is used to examine the effects of static electron correlation on the prediction of a benchmark set of experimentally observed molecular vibrational frequencies. The B3LYP and B97-1 thermally-assisted-occupation measure of static electron correlation is important for describing the vibrations of many of the molecules that make up several popular test sets of experimental data. Shifts are seen for known multireference systems and for many molecules containing atoms from the second row of the periodic table of elements. Several molecules only show significant shifts in select vibrational modes, and significant improvements are seen for the prediction of hydrogen stretching frequencies throughout the test set.

Investigation of the Failure of the MP2 Method to Describe the Out-of-Plane Bending Motions of Carbon-Carbon Double-Bonded Molecules: The Role of Atomic Orbitals

We present a detailed study on the role of atomic orbitals in the unphysical description of the out-of-plane bending (OOPB) vibrations with the MP2 method. The anharmonicities of the OOPB vibrations are found to be unphysically large for several basis sets with the MP2 method. We find that the inclusion of additional valence s and p orbitals to the basis set leads to the lowering of π* orbitals' energies along with the generation of several spurious low-energy virtual molecular orbitals (VMOs). Such erroneous VMOs create a disbalance between the σ and π correlations for the planar structure of the molecules. Since the OOPB vibrations distort the π* orbitals, the disbalance in the σ and π correlations leads to errors in the derivatives of the correlation energy with respect to the displacement along with OOPB modes. The diffuse functions in the basis set enhance the unphysical anharmonicities of the OOPB modes since these functions are of s- and p-type atomic orbitals. The polarization atomic orbitals of higher angular momentum improve the description of the π* orbitals and thereby reduce the disbalance between the σ and π correlations in the MP2 calculations. We find that the unphysical frequencies for the OOPB modes are significantly eliminated with the d and f orbitals on the C atoms in the basis set.

Understanding of the C-H stretch region of infra-red spectroscopy: an analysis of the final state wavefunctions

The nature of the wavefunctions associated with the final states in the CH stretching region of several medium sized molecules is analysed. The number of optically bright transitions is much larger than the number of CH oscillators present in the molecule, and they are spread over a range of about 300 cm^-1. Several of them are clustered together within about 5 cm^-1 with near equal intensities. The final states of all these transitions are superpositions of multiple zeroth order states. In almost all of such superpositions, no single zeroth order state has more than 50% weight. Several multiquantum states, with three to four quanta of excitation dominate the final states, with the CH chromophore contributing only a small weightage. Thus the band structure of the CH stretch region is due to several optically bright transitions whose final states are superpositions of low frequency multiquantum states with the CH chromophore contributing only a small weight to make them spectroscopically active.

Intrinsic structure of pentapeptide Leu-enkephalin: geometry optimization and validation by comparison of VSCF-PT2 calculations with cold ion spectroscopy

The intrinsic structure of an opioid peptide [Ala2, Leu5]-leucine enkephalin (ALE) has been investigated using first-principles based vibrational self-consistent field (VSCF) theory and cold ion spectroscopy. IR-UV double resonance spectroscopy revealed the presence of only one highly abundant conformer of the singly protonated ALE, isolated and cryogenically cooled in the gas phase. High-level quantum mechanical calculations of electronic structures in conjunction with a systematic conformational search allowed for finding a few low-energy candidate structures. In order to identify the observed structure, we computed vibrational spectra of the candidate structures and employed the theory at the semi-empirically scaled harmonic level and at the first-principles based anharmonic VSCF levels. The best match between the calculated "anharmonic" and the measured spectra appeared, indeed, for the most stable candidate. An average of two spectra calculated with different quantum mechanical potentials is proposed for the best match with experiment. The match thus validates the calculated intrinsic structure of ALE and demonstrates the predictive power of first-principles theory for solving structures of such large molecules.

Development and Validation of the B3LYP/N07D Computational Model for Structural Parameter and Magnetic Tensors of Large Free Radicals

Extensive calculations on a large set of free radicals containing atoms of the second and third row show that the B3LYP/N07D computational model provides remarkably accurate structural parameters and magnetic tensors at reasonable computational costs. The key of this success is the optimization of core-valence s functions for hyperfine coupling constants, while retaining (and even improving) the good performances of the parent 6-31+G(d,p) basis set for valence properties through reoptimization of polarization and diffuse p functions.

On the Calculation of Vibrational Frequencies for Molecules in Solution Beyond the Harmonic Approximation

We report some results on the calculation of vibrational spectra of molecules in condensed phase with accounting simultaneously for anharmonicity and solute-solvent interactions, the latter being described by means of the polarizable continuum model (PCM). Density functional theory force fields are employed as well as a new implementation of the PCM cavity and its derivatives. The results obtained for formaldehyde and simple peptide prototypes show that our approach is able to yield a quantitative agreement with experiments for vacuo-to-solvent harmonic and anharmonic frequency shifts.

Pushing the limits in accurate vibrational structure calculations: anharmonic frequencies of lithium fluoride clusters (LiF)n, n = 2-10

The vibrational spectra of a series of small lithium fluoride clusters, i.e. (LiF)n, n = 2-10, were studied by vibrational configuration interaction (VCI) calculations relying on potential energy surfaces including three-mode coupling terms and being obtained from explicitly correlated local coupled cluster calculations. Due to the account for anharmonicity effects, the simulated spectra allow for a direct comparison with experimental data and may thus help to identify clusters in the experiments. Even structurally closely related clusters can clearly be distinguished by infrared spectroscopy. The largest system in this study required more than 1000 basis functions in the electronic structure calculations and more than 10(7) configurations in the vibrational structure calculations and became computationally feasible only due to a combination of different approximations and highly parallelized algorithms.

Solvent induced conformational fluctuation of alanine dipeptide studied by using vibrational probes

The solvation effect on the three dimensional structure and the vibrational feature of alanine dipeptide (ALAD) was evaluated by applying the implicit solvents from polarizable continuum solvent model (PCM) through ab initio calculations, by using molecular dynamic (MD) simulations with explicit solvents, and by combining these two approaches. The implicit solvent induced potential energy fluctuations of ALAD in CHCl3, DMSO and H2O are revealed by means of ab initio calculations, and a global view of conformational and solvation environmental dependence of amide I frequencies is achieved. The results from MD simulations with explicit solvents show that ALAD trends to form PPII, αL, αR, and C5 in water, PPII and C5 in DMSO, and C5 in CHCl3, ordered by population, and the demonstration of the solvated structure, the solute-solvent interaction and hydrogen bonding is therefore enhanced. Representative ALAD-solvent clusters were sampled from MD trajectories and undergone ab initio calculations. The explicit solvents reveal the hydrogen bonding between ALAD and solvents, and the correlation between amide I frequencies and the CO bond length is built. The implicit solvents applied to the ALAD-solvent clusters further compensate the solvation effect from the bulk, and thus enlarge the degree of structural distortion and the amide I frequency red shift. The combination of explicit solvent in the first hydration shell and implicit solvent in the bulk is helpful for our understanding about the conformational fluctuation of solvated polypeptides through vibrational probes.

Accurate molecular structures and infrared spectra of trans-2,3-dideuterooxirane, methyloxirane, and trans-2,3-dimethyloxirane

Oxirane derivatives are the most used benchmarks for chiroptical spectroscopies in view of their small size and relative rigidity. The molecular structure, vibrational harmonic and anharmonic frequencies, and infrared intensities of the ground electronic states are analyzed in this paper. Equilibrium structure and harmonic force fields have been evaluated by means of high-level quantum-chemical calculations at the coupled-cluster level including single and double excitations together with a perturbative treatment of triples (CCSD(T)). Extrapolation to the complete basis-set limit as well as core-correlation effects have also been taken into account. Anharmonic contributions have been computed at the CCSD(T)/cc-pVTZ level for trans-2,3-dideuterooxirane. These data can serve as references to evaluate the accuracy of less expensive computational approaches rooted in the density functional theory (DFT). The latter have been used within hybrid CC/DFT approaches, which have been applied to simulate fully anharmonic infrared (IR) spectra. Finally, the best theoretical estimates of the equilibrium structures and vibrational wavenumbers are compared to the most accurate experimental data and show in all cases very good agreement, i.e., within 0.001 Å, 0.1 deg, 10 cm(-1), and 0.5 km mol(-1), for bond lengths, angles, wavenumbers, and IR intensities, respectively.

Dispersion corrected DFT approaches for anharmonic vibrational frequency calculations: nucleobases and their dimers

Computational spectroscopy techniques have become in the last few years an effective means to analyze and assign infrared (IR) spectra of molecular systems of increasing dimensions and in different environments. However, transition from compilation of harmonic data to fully anharmonic simulations of spectra is still underway. The most promising results for large systems have been obtained, in our opinion, by perturbative vibrational approaches based on potential energy surfaces computed by hybrid (especially B3LYP) density functionals and medium size (e.g. SNSD) basis sets. In this framework, we are actively developing a comprehensive and robust computational protocol aimed at quantitative reproduction of the spectra of nucleic acid base complexes and their adsorption on solid supports (organic/inorganic). In this contribution we report the essential results of the first step devoted to isolated monomers and dimers. It is well known that in order to model the vibrational spectra of weakly bound molecular complexes dispersion interactions should be taken into proper account. In this work we have chosen two popular and inexpensive approaches to model dispersion interactions, namely the semi-empirical dispersion correction (D3) and pseudopotential based (DCP) methodologies both in conjunction with the B3LYP functional. These have been used for simulating fully anharmonic IR spectra of nucleobases and their dimers through generalized second order vibrational perturbation theory (GVPT2). We have studied, in particular, isolated adenine, hypoxanthine, uracil, thymine and cytosine, the hydrogen-bonded and stacked adenine and uracil dimers, and the stacked adenine-naphthalene heterodimer. Anharmonic frequencies are compared with standard B3LYP results and experimental findings, while the computed interaction energies and structures of complexes are compared to the best available theoretical estimates.

Stereo-electronic, vibrational, and environmental contributions to polarizabilities of large molecular systems: a feasible anharmonic protocol

Reliable computations of linear and non-linear optical properties of molecular systems in condensed phases require a proper account of stereo-electronic, vibrational, and environmental effects. In the framework of density functional theory, these effects can be accurately introduced using second-order vibrational perturbation theory in conjunction with polarizable continuum models. We illustrate the combination of an anharmonic description of the ground-state potential energy surface with solvation effects treated with the polarizable continuum model (PCM) in the calculation of the electronic, zero-point, and pure vibrational polarizabilities of selected systems. The description of the solvation environment is enriched by taking into account the dynamical aspects of the solute-solvent interactions through the inclusion of both electronic and vibrational non-equilbrium effects, as well as the direct effect of the solvent on the electric field that generates the molecular response (local field effect). This treatment yields accurate results which can be directly compared with experimental findings without the need of empirical corrections.

Fully anharmonic IR and Raman spectra of medium-size molecular systems: accuracy and interpretation

Computation of full infrared (IR) and Raman spectra (including absolute intensities and transition energies) for medium- and large-sized molecular systems beyond the harmonic approximation is one of the most interesting challenges of contemporary computational chemistry. Contrary to common beliefs, low-order perturbation theory is able to deliver results of high accuracy (actually often better than those issuing from current direct dynamics approaches) provided that anharmonic resonances are properly managed. This perspective sketches the recent developments in our research group toward the development of a robust and user-friendly virtual spectrometer rooted in second-order vibrational perturbation theory (VPT2) and usable also by non-specialists essentially as a black-box procedure. Several examples are explicitly worked out in order to illustrate the features of our computational tool together with the most important ongoing developments.

A robust and effective time-independent route to the calculation of Resonance Raman spectra of large molecules in condensed phases with the inclusion of Duschinsky, Herzberg-Teller, anharmonic, and environmental effects

We present an effective time-independent implementation to model vibrational resonance Raman (RR) spectra of medium-large molecular systems with the inclusion of Franck-Condon (FC) and Herzberg-Teller (HT) effects and a full account of the possible differences between the harmonic potential energy surfaces of the ground and resonant electronic states. Thanks to a number of algorithmic improvements and very effective parallelization, the full computations of fundamentals, overtones, and combination bands can be routinely performed for large systems possibly involving more than two electronic states. In order to improve the accuracy of the results, an effective inclusion of the leading anharmonic effects is also possible, together with environmental contributions under different solvation regimes. Reduced-dimensionality approaches can further enlarge the range of applications of this new tool. Applications to imidazole, pyrene, and chlorophyll a1 in solution are reported, as well as comparisons with available experimental data.

Development of a virtual spectrometer for chiroptical spectroscopies: the case of nicotine

The impressive advances of computational spectroscopy in most recent years are providing robust and user-friendly multifrequency virtual spectrometers, which can also be used by nonspecialists to complement experimental studies. At the heart of these developments there are latest-generation models based on Density Functional Theory for the proper treatment of stereo-electronic effects, coupled to the polarizable continuum model to deal with bulk solvent effects, and low-order perturbative treatments of anharmonic effects. Continuing our efforts to increase the range of application of virtual spectrometers, we report here about chiroptical spectroscopies with special reference to optical rotation and vibrational circular dichroism. The capabilities and possible limitations of our latest tool will be analyzed for the specific case of (S)-nicotine in vacuo and in different solvents.

Infrared and Raman Spectra of and Isotopomers: A DFT-PT2 Anharmonic Study

IR and Raman spectra of selenophene and of its perdeuterated isotopomer have been obtained in gas phase through density-functional theory (DFT) computations. Vibrational wavenumbers have been calculated using harmonic and anharmonic second-order perturbation theory (PT2) procedures with the B3LYP method and the 6-311 basis set. Anharmonic overtones have been determined by means of the PT2 method. The introduction of anharmonic terms decreases the harmonic wavenumbers, giving a significantly better agreement with the experimental data. The most significant anharmonic effects occur for the C–H and C–D stretching modes, the observed H/D isotopic wavenumber redshifts being satisfactorily reproduced by the PT2 computations within 6–20 cm−1(1–3%). In the spectral region between 500 cm−1and 1500 cm−1, the IR spectra are dominated by the out-of-plane C–H (C–D) bending transition, whereas the Raman spectra are mainly characterized by a strong peak mainly attributed to the C=C + C–C bonds stretching vibration with the contribution of the in-plane C–H (C–D) bending deformation. The current results confirm that the PT2 approach combined with the B3LYP/6-311 level of calculation is a satisfactory choice for predicting vibrational spectra of cyclic molecules.

Numerical-analytic implementation of the higher-order canonical Van Vleck perturbation theory for the interpretation of medium-sized molecule vibrational spectra

Anharmonic vibrational states of semirigid polyatomic molecules are often studied using the second-order vibrational perturbation theory (VPT2). For efficient higher-order analysis, an approach based on the canonical Van Vleck perturbation theory (CVPT), the Watson Hamiltonian and operators of creation and annihilation of vibrational quanta is employed. This method allows analysis of the convergence of perturbation theory and solves a number of theoretical problems of VPT2, e.g., yields anharmonic constants y(ijk), z(ijkl), and allows the reliable evaluation of vibrational IR and Raman anharmonic intensities in the presence of resonances. Darling-Dennison and higher-order resonance coupling coefficients can be reliably evaluated as well. The method is illustrated on classic molecules: water and formaldehyde. A number of theoretical conclusions results, including the necessity of using sextic force field in the fourth order (CVPT4) and the nearly vanishing CVPT4 contributions for bending and wagging modes. The coefficients of perturbative Dunham-type Hamiltonians in high-orders of CVPT are found to conform to the rules of equality at different orders as earlier proven analytically for diatomic molecules. The method can serve as a good substitution of the more traditional VPT2.

Vibrational averaging of the chemical shift in crystalline α-glycine

Averaging of the chemical shift over the molecular motion improves the simulated data and provides additional information about the temperature dependence and system dynamics. However, crystal modeling is difficult due to the limited precision of the plane-wave density functional theory (DFT) methods and approximate vibrational schemes. On the glycine example, we investigate how the averaging can be achieved within the periodic boundary conditions at the DFT level. The nuclear motion is modeled with the vibrational configuration interaction, with other simplified quantum anharmonic schemes, and the classical Born-Oppenheimer molecular dynamics (BOMD). The results confirm a large vibrational contribution to the isotropic shielding values. Both the first and second derivatives of the shielding were found important for the quantum averaging. The first derivatives influence the shielding mostly due to the anharmonic character of the CH and NH stretching modes, whereas second derivatives produce most vibrational corrections associated with the lower-frequency vibrational modes. Temperature excitations of the lowest-frequency vibrational states and the expansion of the crystal cell both determine the temperature dependence of nuclear magnetic resonance parameters. The vibrational quantum approach as well as classical BOMD schemes provided temperature dependencies of the chemical shifts that are consistent with the previous experimental data.

Structures, energetics and vibrational frequency shifts of hydrated pyrimidine

More than 70 unique micro-hydrated structures of pyrimidine, ranging in size from 1 to 7 water molecules, have been characterized with the B3LYP density functional and the 6-311++G(2df,2pd) triple-ζ split-valence basis set. Explicitly correlated MP2-F12 single-point computations were performed on each structure with a correlation consistent triple-ζ basis set to estimate the relative and dissociation energies at the MP2 complete basis set (CBS) limit. Many of these new structures have significantly lower energies than those previously reported (by as much as 12.66 kcal mol(-1)). For clusters with 1 and 2 water molecules, the MP2-F12 relative and dissociation energies are virtually identical to the corresponding CCSD(T)-F12 values. As the number of hydrating waters increases, the structures in which the water molecules are clustered together at one of the N atoms have lower energies than those where the water molecules are more distributed around the pyrimidine ring. Micro-hydrated structures that effectively extend the low-energy hydrogen-bonding motifs to both sides of the ring, as would be expected in the bulk phase, reproduce the experimentally observed vibrational frequency shifts of ν(1) and ν(8b) in very dilute aqueous pyrimidine solutions to within 1 cm(-1) . Micro-hydrated structures of pyrimidine in which water molecules are clustered together have lower energies than structures in which the water molecules are more evenly spread around the pyrimidine ring.

Structure and vibrational spectroscopy investigation of 2-(4-chlorophenyliminomethyl)-8-hydroxyquinoline

2-(4-Chlorophenyliminomethyl)-8-hydroxyquinoline was synthesized and crystal data was obtained in the orthorhombic space group P-1, with Z=4. Unit cell parameters a=4.744(7) Å, b=9.981(15) Å, c=27.27(4) Å and V=12915(3) Å(3). In this paper the structural properties and vibrational frequencies of 8-hydroxyquinoline derivative, 2CP8HQ, are studied with the B3LYP and HF methods. Two stable conformers are obtained. The calculated frequencies are in good agreement with the experiment results. It is indicated that both of theoretical calculations were suitable for molecular vibrational frequencies study and the scaled B3LYP method was superior to the scaled HF methods.