Variational calculation of vibrational linear and nonlinear optical properties

A variational approach for reliably calculating vibrational linear and nonlinear optical properties of molecules with large electrical and/or mechanical anharmonicity is introduced. This approach utilizes a self-consistent solution of the vibrational Schrodinger equation for the complete field-dependent potential-energy surface and, then, adds higher-level vibrational correlation corrections as desired. An initial application is made to static properties for three molecules of widely varying anharmonicity using the lowest-level vibrational correlation treatment (i.e., vibrational M?ller-Plesset perturbation theory). Our results indicate when the conventional Bishop-Kirtman perturbation method can be expected to break down and when high-level vibrational correlation methods are likely to be required. Future improvements and extensions are discussed.     

Automatic generation of force fields and property surfaces for use in variational vibrational calculations of anharmonic vibrational energies and zero-point vibrational averaged properties

An automatic and general procedure for the calculation of geometrical derivatives of the energy and general property surfaces for molecular systems is developed and implemented. General expressions for an n-mode representation are derived, where the n-mode representation includes only the couplings between n or less degrees of freedom. The general expressions are specialized to derivative force fields and property surfaces, and a scheme for calculation of the numerical derivatives is implemented. The implementation is interfaced to electronic structure programs and may be used for both ground and excited electronic states. The implementation is done in the context of a vibrational structure program and can be used in combination with vibrational self-consistent field (VSCF), vibrational configuration interaction (VCI), vibrational Moller-Plesset, and vibrational coupled cluster calculations of anharmonic wave functions and calculation of vibrational averaged properties at the VSCF and VCI levels. Sample calculations are presented for fundamental vibrational energies and vibrationally averaged dipole moments and frequency dependent polarizabilities and hyperpolarizabilities of water and formaldehyde.     

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.     

Theory of time-resolved sum-frequency generation and its applications to vibrational dynamics of water

A molecular theory of time-resolved sum-frequency generation (SFG) has been developed. The theoretical framework is constructed using the coupled-oscillator model in the adiabatic approximation. This theory can treat not only the vibrational spectroscopy but also vibrational dynamics. An application of this theory is also provided for estimation of the time constants of the intermolecular vibrational energy transfer between water molecules. This approach can be used for molecular analysis of the experimental results of Shen at al. on the SFG studies of vibrational dynamics of water.     

Moment expansion method to calculate vibrational frequencies of molecules and solids

The vibrational properties of a quantum system are determined by the density response matrix. In linear response theory this quantity is connected to the polarizability matrix, which can be expressed in terms of a double summation over one-particle energies and wave functions. In has been shown that this expression is not useful in the calculation of vibrational frequencies because of the very slow convergence of the summation in terms of the unoccupied states. In this paper, a different but equivalent expression is presented using a continued fraction. The resulting expression contains only one summation over the occupied states, solving in this way all the problems connected with the sum-over-states expression of the polarizability matrix. The elimination of all the unoccupied states via the use of the moment formula turns out to be a crucial step in the solution of the problem of the first-principles calculation of the vibrational spectra of molecules and solids.     

Relating normal vibrational modes to local vibrational modes with the help of an adiabatic connection scheme

Information on the electronic structure of a molecule and its chemical bonds is encoded in the molecular normal vibrational modes. However, normal vibrational modes result from a coupling of local vibrational modes, which means that only the latter can provide detailed insight into bonding and other structural features. In this work, it is proven that the adiabatic internal coordinate vibrational modes of Konkoli and Cremer [Int. J. Quantum Chem. 67, 29 (1998)] represent a unique set of local modes that is directly related to the normal vibrational modes. The missing link between these two sets of modes are the compliance constants of Decius, which turn out to be the reciprocals of the local mode force constants of Konkoli and Cremer. Using the compliance constants matrix, the local mode frequencies of any molecule can be converted into its normal mode frequencies with the help of an adiabatic connection scheme that defines the coupling of the local modes in terms of coupling frequencies and reveals how avoided crossings between the local modes lead to changes in the character of the normal modes.     

Rapid anharmonic vibrational corrections derived from partial Hessian analysis

Vibrational analysis within a partial Hessian framework can successfully describe the vibrational properties of a variety of systems where the vibrational modes of interest are localized within a specific region of the system. We have developed a new approach to calculating anharmonic frequencies based on vibrational frequencies and normal modes obtained from a partial Hessian analysis using second-order vibrational perturbation theory and the transition optimized shifted Hermite method. This allows anharmonic frequencies for vibrational modes that are spatially localized to be determined at a significantly reduced computational cost. Several molecular systems are examined in order to demonstrate the effectiveness of this method including organic molecules adsorbed on the Si(100)-2×1 surface, model peptides in solution, and the C-H stretching region of polycyclic aromatic hydrocarbons. Overall, for a range of systems, anharmonic frequencies calculated using the partial Hessian approach are found to be in close agreement with the results obtained using full anharmonic calculations while providing a significant reduction in computational cost.     

Analytic calculations of vibrational hyperpolarizabilities in the atomic orbital basis

We present an analytic scheme for the calculation of pure vibrational contributions to linear and nonlinear optical properties such as the polarizability and the first and second hyperpolarizabilities. The formalism is fully expressed in terms of a perturbation- and time-dependent atomic orbital basis, using the elements of the density matrix in the atomic orbital basis as the basic variables. We calculate perturbed densities up to third order with respect to the electric field in accordance with the n + 1 rule, and the approach is therefore applicable for the calculation of pure vibrational contributions involving all vibrational coordinates in large molecular complexes. In the case of static electric fields, we therefore only need to calculate 19 response equations, independent of the size of the molecule. If we can determine the molecular energy and force field, the calculation of pure vibrational contributions to the nonlinear optical properties of the molecule is therefore a rather straightforward task. We illustrate the implementation by calculating pure vibrational contributions to the first and second hyperpolarizabilities of molecules containing up to 66 atoms using basis sets of good quality.     

Diamondoids approach to electronic,structural, and vibrational properties of GeSi superlattice nanocrystals: a first-principles study

Germanium silicide diamondoids are used to determine electronic, structural, and vibrational properties of GeSi superlattice nanocrystals and bulk as their building block limit. Density functional theory at the generalized gradient approximation level of Perdew, Burke, and Ernzerhof (PBE) with 6-31G(d) basis including polarization functions is used to investigate the electronic structure of these diamondoids. The investigated molecules and diamondoids range from GeSiH₆ to Ge₆₃Si₆₃H₉₂. The variation of the energy gap is shown from nearly 7 eV toward bulk value which is slightly higher than the average of Si and Ge energy gaps. Variations of bond lengths, tetrahedral, and dihedral angles as the number of atoms increases are shown taking into account the effect of shape fluctuations. Localized and delocalized electronic charge distribution and bonds for these molecules are discussed. Vibrational radial breathing mode (RBM) converges from its initial molecular value at 332 cm^?1 to its bulk limit at 0 cm^?1 (blue shift). Longitudinal optical-highest reduced mass mode (HRMM) converges from its initial molecular value 332 cm^?1 to experimental bulk limit at 420.7 cm^?1 (red shift). Hydrogen vibrational modes are nearly constant in their frequencies as the size of diamondoids increases in contrast with lower frequency Ge–Si vibrational modes. GeSi diamondoids can be identified from surface hydrogen vibrational modes fingerprint, while the size of these diamondoids can be identified from Ge–Si vibrational modes.     

Enhancement factor averaging and the photostability of probes in SERS vibrational pumping

The technique of temperature dependent vibrational pumping in surface enhanced Raman scattering (SERS) has been recently demonstrated as a promising new tool to estimate SERS cross-sections. In this paper we expand on the previous developments and study several details around the implementation and physics of the vibrational pumping technique in SERS. Here we concentrate on two specific aspects related to: (i) the different averaging properties (over the distribution of enhancements) of the Stokes and anti-Stokes signals in the pumping regime; and (ii) the role of the finite photostability of the probes. The fact that the anti-Stokes signal is averaged differently from the Stokes counterpart leads to some unique phenomena in Raman spectroscopy that can only be observed under the conditions of vibrational pumping in SERS.     

Towards fast computations of correlated vibrational wave functions: vibrational coupled cluster response excitation energies at the two-mode coupling level

An efficient implementation of vibrational coupled cluster theory with two-mode excitations and a two-mode Hamiltonian is described. The algorithm is shown to scale cubically with respect to the number of modes which is identical to the scaling of the corresponding vibrational configuration interaction algorithm. This is achieved through the use of special intermediates. The same algorithm can also be used in vibrational M?ller-Plesset calculations. To improve performance, screening techniques have been implemented as well. Test calculations on polyaromatic hydrocarbons with up to 264 coupled modes and model systems with up to 1140 modes are used to illustrate the various features of the algorithm.     

Size-extensive vibrational self-consistent field method

The vibrational self-consistent field (VSCF) method is a mean-field approach to solve the vibrational Schro?dinger equation and serves as a basis of vibrational perturbation and coupled-cluster methods. Together they account for anharmonic effects on vibrational transition frequencies and vibrationally averaged properties. This article reports the definition, programmable equations, and corresponding initial implementation of a diagrammatically size-extensive modification of VSCF, from which numerous terms with nonphysical size dependence in the original VSCF equations have been eliminated. When combined with a quartic force field (QFF), this compact and strictly size-extensive VSCF (XVSCF) method requires only quartic force constants of the ?(4)V/?Q(i)(2)?Q(j)(2) type, where V is the electronic energy and Q(i) is the ith normal coordinate. Consequently, the cost of a XVSCF calculation with a QFF increases only quadratically with the number of modes, while that of a VSCF calculation grows quartically. The effective (mean-field) potential of XVSCF felt by each mode is shown to be harmonic, making the XVSCF equations subject to a self-consistent analytical solution without matrix diagonalization or a basis-set expansion, which are necessary in VSCF. Even when the same set of force constants is used, XVSCF is nearly three orders of magnitude faster than VSCF implemented similarly. Yet, the results of XVSCF and VSCF are shown to approach each other as the molecular size is increased, implicating the inclusion of unnecessary, nonphysical terms in VSCF. The diagrams of the XVSCF energy expression and their evaluation rules are also proposed, underscoring their connected structures.     

Fast and accurate computation schemes for evaluating vibrational entropy of proteins

Standard normal mode analysis (NMA) method is able to calculate vibrational entropy of proteins, but it is computationally intensive, especially for large proteins. To evaluate vibrational entropy efficiently and accurately, we, here, propose computation schemes based on coarse-grained NMA methods. This can be achieved by rescaling coarse-grained results with a specific factor that is derived on the basis of the linear correlation of protein vibrational entropy between standard NMA and coarse-grained NMA. Our coarse-grained NMA computation schemes can repeat correctly and efficiently the results of standard NMA for large proteins.     

Calculations of intermode coupling constants and simulations of amide I,II, and III vibrational spectra of dipeptides

Amide I, II, and III vibrations of polypeptides are important marker modes whose vibrational spectra can provide critical information on structure and dynamics of proteins in solution. The extent of delocalization and vibrational properties of amide normal mode can be described by the amide local mode frequencies and intermode coupling constants between a pair of amide local modes. To determine these fundamental quantities, the previous Hessian matrix reconstruction method has been generalized here and applied to the density functional theory results for various dipeptide conformers. The calculation results are then used to simulate IR absorption, vibrational circular dichroism, and 2D IR spectra of dipeptides. The relationships between dipeptide backbone conformations and these vibrational spectra are discussed. It is believed that the present computational method and results will be of use to quantitatively simulate vibrational spectra of complicated polypeptides beyond simple dipeptides     

Vibrational amplitude profile of molecular vibrational modes for vibrational mode assignment

Ultrashort pulse lasers with 6- and 20-fs durations were utilized for phthalocyanine thin film sample to induce several vibrational modes and vibration amplitude spectra were determined by multi-wavelength measurement technique. From the spectra we could identify the electronic states, which couple to two vibrational modes with frequencies of 670 and 750 cm⁻¹. It was shown that the vibrational amplitude profile obtained by the method can be used for providing information for the assignment of the vibrational mode.     

Accuracy of gates in a quantum computer based on vibrational eigenstates

A model is developed to study the properties of a quantum computer that uses vibrational eigenstates of molecules to implement the quantum information bits and shaped laser pulses to apply the quantum logic gates. Particular emphasis of this study is on understanding how the different factors, such as properties of the molecule and of the pulse, can be used to affect the accuracy of quantum gates in such a system. Optimal control theory and numerical time-propagation of vibrational wave packets are employed to obtain the shaped pulses for the gates NOT and Hadamard transform. The effects of the anharmonicity parameter of the molecule, the target time of the pulse and of the penalty function are investigated. Influence of all these parameters on the accuracy of qubit transformations is observed and explained. It is shown that when all these parameters are carefully chosen the accuracy of quantum gates reaches 99.9%.     

Computation of expectation values from vibrational coupled-cluster at the two-mode coupling level

In this work we show how the vibrational coupled-cluster method at the two-mode coupling level can be used to calculate zero-point vibrational averages of properties. A technique is presented, where any expectation value can be calculated using a single set of Lagrangian multipliers computed solving iteratively a single linear set of equations. Sample calculations are presented which show that the resulting algorithm scales only with the third power of the number of modes, therefore making large systems accessible. Moreover, we present applications to water, pyrrole, and para-nitroaniline.     

Nonperturbative vibrational energy relaxation effects on vibrational line shapes

A general formulation of nonperturbative quantum dynamics of solutes in a condensed phase is proposed to calculate linear and nonlinear vibrational line shapes. In the weak solute-solvent interaction limit, the temporal absorption profile can be approximately factorized into the population relaxation profile from the off-diagonal coupling and the pure-dephasing profile from the diagonal coupling. The strength of dissipation and the anharmonicity-induced dephasing rate are derived in Appendix A. The vibrational energy relaxation (VER) rate is negligible for slow solvent fluctuations, yet it does not justify the Markovian treatment of off-diagonal contributions to vibrational line shapes. Non-Markovian VER effects are manifested as asymmetric envelops in the temporal absorption profile, or equivalently as side bands in the frequency domain absorption spectrum. The side bands are solvent-induced multiple-photon effects which are absent in the Markovian VER treatment. Exact path integral calculations yield non-Lorentzian central peaks in absorption spectrum resulting from couplings between population relaxations of different vibrational states. These predictions cannot be reproduced by the perturbative or the Markovian approximations. For anharmonic potentials, the absorption spectrum shows asymmetric central peaks and the asymmetry increases with anharmonicity. At large anharmonicities, all the approximation schemes break down and a full nonperturbative path integral calculation that explicitly accounts for the exact VER effects is needed. A numerical analysis of the O-H stretch of HOD in D(2)O solvent reveals that the non-Markovian VER effects generate a small recurrence of the echo peak shift around 200 fs, which cannot be reproduced with a Markovian VER rate. In general, the nonperturbative and non-Markovian VER contributions have a stronger effect on nonlinear vibrational line shapes than on linear absorption.     

Diabatic ordering of vibrational normal modes in reaction valley studies

Diabatic ordering of the normal model of a reaction complex along the reaction path has several advantages with regard to adiabatic ordering. The method is based on rotations of the vibrational normal modes at one point, s, of the reaction path to maximize overlap with the vibrational modes at a neighboring point. Global rotations precede the rotations of degenerate modes so that changes in the direction of the reaction path and changes in the force constant matrix, which represent the two major effects for changes in mode ordering, can be separated. Overlap criteria identify resolved and unresolved avoided crossings of normal modes of the same symmetry. Diabatic mode ordering (DMO) can be used to resolve the latter by reducing the step size, thus guaranteeing correct ordering of normal modes in dependence of s. DMO is generally applicable to properties of the reaction complex that depend on s such as normal mode frequencies, orbital energies, the energy of excited states, etc. Additional applications are possible using a generalized reaction path vector, which may describe the change in atom masses, geometrical parameters, and/or the force constant matrix. In this way, the vibrational spectra of isotopomers can be investigated or the vibrational frequencies of different molecules correlated. © 1997 John Wiley & Sons, Inc. J Comput Chem 18: 1282–1294, 1997     

Amide I vibrational circular dichroism of dipeptide: Conformation dependence and fragment analysis

The amide I vibrational circular dichroic response of alanine dipeptide analog (ADA) was theoretically investigated and the density functional theory calculation and fragment analysis results are presented. A variety of vibrational spectroscopic properties, local and normal mode frequencies, coupling constant, dipole, and rotational strengths, are calculated by varying two dihedral angles determining the three-dimensional ADA conformation. Considering two monopeptide fragments separately, we show that the amide I vibrational circular dichroism of the ADA can be quantitatively predicted. For several representative conformations of the model ADA, vibrational circular dichroism spectra are calculated by using both the density functional theory calculation and fragment analysis methods.