On the vibrational linear and nonlinear optical properties of compounds involving noble gas atoms: HXeOXeH,HXeOXeF, and FXeOXeF

The vibrational (hyper)polarizabilities of some selected Xe derivatives are studied in the context of Bishop–Kirtman perturbation theory (BKPT) and numerical finite field methodology. It was found that for this set of rare gas compounds, the static vibrational properties are quite large, in comparison to the corresponding electronic ones, especially those of the second hyperpolarizability. This also holds for the dc‐Pockels β(?ω;ω,0), Kerr γ(?ω;ω,0,0) and electric field second harmonic generation γ (?2ω;ω,ω,0) effects, although the computed nuclear relaxation (nr) vibrational contributions are smaller in magnitude than the static ones. HXeOXeH was used to study the effects of electron correlation, basis set, and geometry. Geometry effects were found to lead to noticeable changes of the vibrational and electronic second hyperpolarizability. A limited study of the effect of Xe insertion to the nr vibrational properties is also reported. Assessment of the results revealed that Xe insertion has a remarkable effect on the nr (hyper)polarizabilities. In terms of the BKPT, this is associated with a remarkable increase of the electrical and mechanical anharmonicity terms. The latter is consistent with the anharmonic character of several vibrational modes reported for rare gas compounds. © 2013 Wiley Periodicals, Inc.     

Classical and quantum mechanical infrared echoes from resonantly coupled molecular vibrations

The nonlinear response function associated with the infrared vibrational echo is calculated for a quantum mechanical model of resonantly coupled, anharmonic oscillators at zero temperature. The classical mechanical response function is determined from the quantum response function by setting variant Planck's over 2pi-->0, permitting the comparison of the effects of resonant vibrational coupling among an arbitrary number of anharmonic oscillators on quantum and classical vibrational echoes. The quantum response function displays a time dependence that reflects both anharmonicity and resonant coupling, while the classical response function depends on anharmonicity only through a time-independent amplitude, and shows a time dependence controlled only by the resonant coupling. In addition, the classical response function grows without bound in time, a phenomenon associated with the nonlinearity of classical mechanics, and absent in quantum mechanics. This unbounded growth was previously identified in the response function for a system without resonant vibrational energy transfer, and is observed to persist in the presence of resonant coupling among vibrations. Quantitative agreement between classical and quantum response functions is limited to a time scale of duration inversely proportional to the anharmonicity.     

Pressure and temperature effects on the vibrational frequencies of crystalline polymethylenes. II. Lattice anharmonicity and the equation of state

The lattice anharmonicity of crystalline polymethylene is interpreted from the observed pressure and temperature dependence of Raman active interchain lattice frequencies of the n-paraffins C₂₃H₄₈ and C₄₄H₉₀. The temperature dependence of the L_c′ interchain lattice frequency is separated into quasiharmonic and self-energy shifts. The former is due to the volume dependence of the force constant of the oscillator. The latter is due to the anharmonicity of the dynamic potential, and is obtained as a function of volume and phonon population. The setting angle of the carbon skeleton is predicted to be temperature-sensitive. While the potential surface of the crystal is asymmetric along the L_c′ normal coordinate, it is essentially symmetrical along the T_b′ coordinate. The well-known Mie–Gruneisen equation of state is generalized to include anharmonicities of oscillators through the temperature dependence of their vibrational frequencies.     

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.     

Derivation of class II force fields. I. Methodology and quantum force field for the alkyl functional group and alkane molecules

A new method for deriving force fields for molecular simulations has been developed. It is based on the derivation and parameterization of analytic representations of the ab initio potential energy surfaces. The general method is presented here and used to derive a quantum mechanical force field (QMFF) for alkanes. It is based on sampling the energy surfaces of 16 representative alkane species. For hydrocarbons, this force field contains 66 force constants and reference values. These were fit to 128,376 quantum mechanical energies and energy derivatives describing the energy surface. The detailed form of the analytic force field expression and the values of all resulting parameters are given. A series of computations is then performed to test the ability of this force field to reproduce the features of the ab initio energy surface in terms of energies as well as the first and second derivatives of the energies with respect to molecular deformations. The fit is shown to be good, with rms energy deviations of less than 7% for all molecules. Also, although only two atom types are employed, the force field accounts for the properties of both highly strained species, such as cyclopropane and methylcyclopropanes, as well as unstrained systems. The information contained in the quantum energy surface indicates that it is significantly anharmonic and that important intramolecular coupling interactions exist between internals. The representation of the nature of these interactions, not present in diagonal, quadratic force fields (Class I force fields), is shown to be important in accounting accurately for molecular energy surfaces. The Class II force field derived from the quantum energy surface is characterized by accounting for these important intramolecular forces. The importance of each 4.2 to 18.2%. This fourfold increase in the second derivative error dramatically demonstrates the importance of bond anharmonicity in the ab initio potential energy surface. The Class II force field derived from the quantum energy surface is characterized by accounting for these important intramolecular forces. The importance of each of the interaction terms of the potential energy function has also been assessed. Bond anharmonicity, angle anharmonicity, and bond/angle, bond/torsion, and angle/angle/ torsion cross-term interactions result in the most significant overall improvement in distorted structure energies and energy derivatives. The implications of each energy term for the development of advanced force fields is discussed. Finally, it is shown that the techniques introduced here for exploring the quantum energy surface can be used to determine the extent of transferability and range of validity of the force field. The latter is of crucial importance in meeting the objective of deriving a force field for use in molecular mechanics and dynamics calculations of a wide range of molecules often containing functional groups in novel environments. © 1994 by John Wiley & Sons, Inc.     

The vibrational energy manifold of the lowest lying bending mode of tricarbon oxide sulfide,OCCCS,as determined by relative intensity measurements

The determination of a precise vibrational energy level scheme for the two-dimensional bending mode of tricarbon oxide sulfide (3-thioxo-1,2-propadiene-1-one), OCCCS, has been carried out by relative intensity measurements of rotational transitions up to the seventh excited vibrational state of ν₇. The harmonic wavenumber ω₇ was determined to be 84.50 ± 0.63 cm^?1 while the anharmonicity constant χ₇₇ was found to be ?0.62 ± 0.11 cm^?1, respectively. A linear dependence of the expectation value of the electric dipole moment on the vibrational quantum number υ₇ was found. All results confirm that in O CCCS the potential function describing the two-dimensional oscillator of ν₇ is very harmonic without a perturbing barrier to linearity as was found in the case of OCCCO.     

Derivation of class II force fields: V. Quantum force field for amides,peptides, and related compounds

As the field of biomolecular structure advances, there is an ever-growing need for accurate modeling of molecular energy surfaces to simulate and predict the properties of these important systems. To address this need, a second generation amide force field for use in simulations of small organics as well as proteins and peptides has been derived. The critical question of what accuracy can be expected from calculations in general, and with this class II force field in particular, is addressed for structural, dynamic, and energetic properties. The force field is derived from a recent methodology we have developed that involves the systematic use of quantum mechanical observables. Systematic ab initio calculations were carried out for numerous configurations of 17 amide and related compounds. Relative energies and first and second derivatives of the energy of 638 structures of these compounds resulted in 140,970 ab initio quantum mechanical observables. The class II peptide quantum mechanical force field (QMFF), containing 732 force constants and reference values, was parameterized against these observables. A major objective of this work is to help establish the role of anharmonicity and coupling in improving the accuracy of molecular force fields, as these terms have not yet become an agreed upon standard in the ever more extensive simulations being used to probe biomolecular properties. This has been addressed by deriving a class I harmonic diagonal force field (HDFF), which was fit to the same energy surface as the QMFF, thus providing an opportunity to quantify the effects of these coupling and anharmonic contributions. Both force field representations are assessed in terms of their ability to fit the observables. They have also been tested by calculating the properties of 11 stationary states of these amide molecules. Optimized structures, vibrational frequencies, and conformational energies obtained from the quantum calculations and from both the QMFF and the HDFF are compared. Several strained and derivatized compounds including urea, formylformamide, and butyrolactam are included in these tests to assess the range of applicability (transferability) of the force fields. It was found that the class II coupled anharmonic force field reproduced the structures, energies, and vibrational frequencies significantly more faithfully than the class I harmonic diagonal force field. An important measure, rms energy deviation, was found to be 1.06 kcal/mol with the class II force field, and 2.30 kcal/mol with the harmonic diagonal force field. These deviations represent the error in relative configurational energy differences for strained and distorted structures calculated with the force fields compared with quantum mechanics. This provides a measure of the accuracy that might be expected in applications where strain may be important such as calculating the energy of a system as it approaches a (rotational) barrier, in ligand binding to a protein, or effects of introducing substituents into a molecule that may induce strain. Similar results were found for structural properties. Protein dynamics is becoming of ever-increasing interest, and, to simulate dynamic properties accurately, the dynamic behavior of model compounds needs to be well accounted for. To this end, the ability of the class I and class II force fields to reproduce the vibrational frequencies obtained from the quantum energy surface was assessed. An rms deviation of 43 cm⁻¹ was achieved with the coupled anharmonic force field, as compared to 105 cm⁻¹ with the harmonic diagonal force field. Thus, the analysis presented here of the class II force field for the amide functional group demonstrates that the incorporation of anharmonicity and coupling terms in the force field significantly improves the accuracy and transferability with regard to the simulation of structural, energetic, and dynamic properties of amides. © 1998 John Wiley & Sons, Inc. J Comput Chem 19: 430–458, 1998     

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.     

Electrostatic calculation of linear and non-linear optical properties of ice Ih,II, IX and VIII

Macroscopic first- and third-order susceptibilities of ice Ih, ice II, ice IX and ice VIII are calculated using static and frequency-dependent electronic and static vibrational molecular (hyper)polarizabilities at the MP2 level. The molecular properties are in good agreement with experiment and with high-level ab initio calculations. Intermolecular electrostatic and polarization effects due to induced dipoles are taken into account using a rigorous local-field theory. The electric field due to permanent dipoles is used to calculate effective in-crystal (hyper)polarizabilities. The polarizability depends only weakly on the permanent field, but the dipole moment and the hyperpolarizabilities are strongly affected. The calculated linear susceptibility is in good agreement with available experimental data for ice Ih, and the third-order susceptibility for a third harmonic generation experiment is in reasonable agreement with experimental values for liquid water. The molecular vibrational contributions have a small effect on the susceptibilities. The electric properties of a water tetramer are calculated and used to estimate the effect of non-dipolar interactions on the susceptibilities of ice Ih, which are found to be small.     

Spectroscopic and theoretical optical properties of indoleninyl-substituted dibenzotetraaza[14]annulenes

Two new macrocyclic dibenzotetraaza[14]annulene (DBTAA) compounds with indolenine ( 5 ) and pyridoindolenine ( 6 ) moieties were synthesized and characterized by spectroscopy. Both DBTAAs exhibit strong UV-Vis absorption properties in the Soret band region. The theoretical second-order nonlinear optical property, electric dipole moment (μ), dispersion-free dipole polarizability (α) and first hyper-polarizability values were calculated by density functional theory and time dependent density functional theory. The ab-initio quantum mechanical calculation by time-dependent Hartree-Fock method was utilized to investigate the dynamic dipole polarizabilities, dynamic second-order, static, and dynamic third-order (γ) hyper-polarizabilities of the DBTAAs. The configuration interaction technique of all doubly occupied molecular orbitals possesses theoretically defined single-photon absorption (OPA) specifications for the examined structures. The computed maximum OPA wavelengths on both macrocyclic compounds coincide with the preceding measurement outcomes.     

MCSCF/MCLR Studies of potential energy surfaces,spectra, and properties of the X1A1 and a3B2 states of ozone

Potential energy surfaces, properties, and spectra of singlet (X¹A₁) and triplet (a³B₂) ozone are investigated by means of MCSCF /MCLR analytical response theory calculations. MCSCF analytical gradients and Hessians are used to locate equilibrium and transition-state structures and to obtain associated vibrational and rotational constants, infrared intensities, and dipole moments. By means of MC linear response functions, electronic excitation energies, and oscillator strengths, static and dynamic polarizabilities as well as dispersion (C₆) coefficients are obtained. Good agreement is achieved, in some cases within experimental error margins, for properties where experimental data are known. A very low IR intensity for triplet ozone is predicted.     

Vibronic absorption spectra of acenaphthene in solution and its excited state electric dipole moments and polarizabilities

The vibronic spectra of acenaphthene in solution have been studied in detail in the region 27778–50000 cm⁻¹. A vibronic analysis of the two longest-wavelength absorption bands was made to reveal the vibrational modes that contribute to the enhancement of the intensities of these bands. The oscillator strengths of the various electronic transitions and the electric dipole moments and polarizabilities of several excited states were determined, the latter two by the solvent spectral frequency shift method.     

Second-order Brillouin-Wigner perturbation theory: size-extensivity correction

A convergence study of intensities of transitions from vibrational ground state to the lower lying states is done with respect to the rank of the dipole moment surface (DMS) in the Taylor series expansion of the DMS. The relative roles of the mechanical and electrical anharmonicity are analyzed in the calculation of the intensities of vibrational transitions from ground state. We find that at least a quadratic expansion of the dipole moment functional is necessary to predict the intensities of vibrational transitions. The mechanical anharmonicity becomes important when the resonances between the vibrational states are significant.     

Static electric dipole polarizability of small sodium aggregates

We present pseudopotential local-spin-density calculations of the static electric polarizability of sodium dimers and trimers and their respective cations. The electronic polarizabilities are obtained from self-consistent calculations in the presence of an external electric field, which is kept sufficiently small to avoid non-linear effects. The calculated polarizability tensor has a strong anisotropy directly related to the geometric and electronic structures of the molecules, the anisotropy being larger for the neutral clusters. The polarizabilities are averaged over the vibrational motion and rotations of the aggregates in order to be compared with the experimental measurements. The obtained values show an improvement in the agreement with experiment with respect to the values calculated in the spherical approximation.     

Study of the Influence of the Nature and Position of the Halogen Atom in the N-Halogenophenyl-2,5-dimethylpyrrole Upon the Vibrational (Hyper)polarizability: Theoretical Study

The two contributions, vibrational and electronic, to the electrical properties polarizability and first hyperpolarizability of the N-[(2, 3, or 4)-fluorophenyl]-2,5-dimethylpyrrole are evaluated theoretically at the HF/6-31G level within the double harmonic oscillator approximation. The calculations demonstrate that, with the exception of the second harmonic generation, the vibrational contribution to the first hyperpolarizability is important. However, the vibrational polarizability, contributes at most, 10% to the total electric polarizability. The analysis upon the sum-over-states expressions shows that there are only few modes which contribute strongly and, generally, have small vibrational energies. The effect of the fluorine substitution by an other substituent is also addressed.     

Tests of semiclassical treatments of vibrational-translational energy transfer collinear collisions of helium with hydrogen molecules

Three kinds of semiclassical theory are tested against quantum mechanical results for vibrational transition probabilities and average vibrational energy transfers in collinear collisions of atoms with harmonic and Morse vibrators for the He-H₂ mass combination. The interaction potential is assumed to be a repulsive exponential function with an exponential parameter which is realistic for He-H₂ collisions. The energy range studied is total energies of 2–8 in units of ?ω_e. The uniform semiclassical approximations of classical S matrix theory are tested only for classically allowed transitions, i.e., for transition probabilities greater than about 0.2. They are accurate quantitatively for both harmonic and Morse vibrators. The integral expressions of classical S matrix theory are found to be quantitatively accurate for classically allowed and weakly classically forbidden transitions, i.e., for transition probabilities greater than about 0.01–0.05, and to be unreliable for strongly classically forbidden transitions. Quasiclassical trajectory methods yield qualitatively accurate results only for classically allowed transitions but the phase-averaged energy transfer in quasiclassical collisions may be accurate even when classically forbidden transition probabilities are important for the calculation of the average energy transfer. Forced quantum oscillator methods using a classical path whose initial velocity is the average of the initial and final velocities corresponding to the transition of interest are accurate for transition probabilities as small as 4 × 10^?8 for harmonic vibrators but do not seem to accurately account for the effect of anharmonicity.     

Rotational-vibrational rainbows in impulsive electron-diatomic molecule collisions

Rotational and vibrational rainbow effects in electron-diatomic molecule scattering at intermediate impact energies (≈10² eV) are discussed in a simple quantum mechanical spectator model within the rigid rotor/harmonic oscillator approximation. The total vibrational (summed over all final rotational quantum numbers) and rotational (vibrationally summed) transition probabilities show vibrational or rotational rainbow patterns, characteristic steps, and rainbow singularities, which are analyzed and interpreted in terms of classical cross sections.     

The effects of static quartic anharmonicity on the quantum dynamics of a linear oscillator with time-dependent harmonic frequency: Perturbative analysis and numerical calculations

The effects of quartic anharmonicity on the quantum dynamics of a linear oscillator with time-dependent force constant (K) or harmonic frequency (ω) are studied both perturbatively and numerically by the time-dependent Fourier grid Hamiltonian method. In the absence of anharmonicity, the ground-state population decreases and the population of an accessible excited state (k = 2, 4, 6…) increases with time. However, when anharmonicity is introduced, both the ground- and excited-state populations show typical oscillations. For weak coupling, the population of an accessible excited state at a certain instant of time (short) turns out to be a parabolic function of the anharmonic coupling constant (λ), when all other parameters of the system are kept fixed. This parabolic nature of the excited-state population vs. the λ profile is independent of the specific form of the time dependence of the force constant, K_t. However, it depends upon the rate at which K_t relaxes. For small anharmonic coupling strength and short time scales, the numerical results corroborate expectations based on the first-order time-dependent perturbative analysis, using a suitably repartitioned Hamiltonian that makes H₀ time-independent. Some of the possible experimental implications of our observations are analyzed, especially in relation to intensity oscillations observed in some charge-transfer spectra in systems in which the dephasing rates are comparable with the time scale of the electron transfer. © 1995 John Wiley & Sons, Inc.     

Application of the effective harmonic oscillator method to computing the thermal averages for a system of coupled anharmonic oscillators

By treating the Hamiltonian for coupled oscillators with polynomial anharmonicity by the Gibbs-Bogoliubov inequality, the effective harmonic oscillator (EHO) method is developed and applied to computing the thermal averages for polyatomic molecules. Practical utility is demonstrated with calculations of electron diffraction quantities, namely the distance r_a and amplitude l, and of the vibrational partition functions for CO₂, CS₂, SO₂ and H₂O from spectroscopic data on the force fields. The results are compared with those in the literature obtained by more accurate techniques. A comparison of r_a and l was also made with the results of electron diffraction measurements.     

Quantum rigid dipole in a permanent electric field. II. Model of librational motion in liquid water and ice Ih: Preliminary results

A theoretical one-body model of librational motion in liquid water and ice Ih is proposed within the quantum mechanical treatment of a rigid and pointlike dipole in a permanent electric field and the concept of the internal electric field justified for V-structure of liquid water.