Analyzing atomic liquids and solids by means of two-dimensional Raman spectra in frequency domain

A practical method to evaluate the contributions of the nonlinear polarizability and anharmonicity of potentials from the experimental and simulation data by using double Fourier transformation is presented. In a Lennard-Jones potential system, an approximated expression of the fifth-order response function using the ratio between nonlinear polarizability and anharmonicity exhibits a good agreement with the results of the molecular dynamics simulation. In a soft-core case, the fifth-order Raman signal indicates that the system consists of the delocalized and localized modes, and only the delocalized mode affects the dramatic change of the fifth-order Raman response functions between solid and liquid phases through nonlinear polarizability.     

Electronic polarization effect on low-frequency infrared and Raman spectra of aprotic solvent: molecular dynamics simulation study with charge response kernel by second order Møller-Plesset perturbation method

Low-frequency infrared (IR) and depolarized Raman scattering (DRS) spectra of acetonitrile, methylene chloride, and acetone liquids are simulated via molecular dynamics calculations with the charge response kernel (CRK) model obtained at the second order M?ller-Plesset perturbation (MP2) level. For this purpose, the analytical second derivative technique for the MP2 energy is employed to evaluate the CRK matrices. The calculated IR spectra reasonably agree with the experiments. In particular, the agreement is excellent for acetone because the present CRK model well reproduces the experimental polarizability in the gas phase. The importance of interaction induced dipole moments in characterizing the spectral shapes is stressed. The DRS spectrum of acetone is mainly discussed because the experimental spectrum is available only for this molecule. The calculated spectrum is close to the experiment. The comparison of the present results with those by the multiple random telegraph model is also made. By decomposing the polarizability anisotropy time correlation function to the contributions from the permanent, induced polarizability and their cross term, a discrepancy from the previous calculations is observed in the sign of permanent-induce cross term contribution. The origin of this discrepancy is discussed by analyzing the correlation functions for acetonitrile.     

Raman spectrum modification in the presence of solvated electrons. Observation by polarization CARS

The multiple-scattering series for depolarized light scattering from simple liquids is investigated. It is found that the leading term in the series, due to double—double scattering, is probably not large compared to double—triple (DT) scattering contributions. Existing theories, which neglect DT scattering, have concluded that the experimental data cannot be explained using the dipole—induced dipole (DID) model for the pair polarizability. It is suggested that a theory which includes DT scattering will not lead to rejection of the DID model.     

Surface enhanced Raman optical activity of molecules on orientationally averaged substrates: theory of electromagnetic effects

We present a model for electromagnetic enhancements in surface enhanced Raman optical activity (SEROA) spectroscopy. The model extends previous treatments of SEROA to substrates, such as metal nanoparticles in solution, that are orientationally averaged with respect to the laboratory frame. Our theoretical treatment combines analytical expressions for unenhanced Raman optical activity with molecular polarizability tensors that are dressed by the substrate's electromagnetic enhancements. We evaluate enhancements from model substrates to determine preliminary scaling laws and selection rules for SEROA. We find that dipolar substrates enhance Raman optical activity (ROA) scattering less than Raman scattering. Evanescent gradient contributions to orientationally averaged ROA scale to first or higher orders in the gradient of the incident plane-wave field. These evanescent gradient contributions may be large for substrates with quadrupolar responses to the plane-wave field gradient. Some substrates may also show a ROA contribution that depends only on the molecular electric dipole-electric dipole polarizability. These conclusions are illustrated via numerical calculations of surface enhanced Raman and ROA spectra from (R)-(-)-bromochlorofluoromethane on various model substrates.     

Two-dimensional Raman spectra of atomic solids and liquids

We calculate third- and fifth-order Raman spectra of simple atoms interacting through a soft-core potential by means of molecular-dynamics (MD) simulations. The total polarizability of molecules is treated by the dipole-induced dipole model. Two- and three-body correlation functions of the polarizability at various temperatures are evaluated from equilibrium MD simulations based on a stability matrix formulation. To analyze the processes involved in the spectroscopic measurements, we divide the fifth-order response functions into symmetric and antisymmetric integrated response functions; the symmetric one is written as a simple three-body correlation function, while the antisymmetric one depends on a stability matrix. This analysis leads to a better understanding of the time scales and molecular motions that govern the two-dimensional (2D) signal. The 2D Raman spectra show novel differences between the solid and liquid phases, which are associated with the decay rates of coherent motions. On the other hand, these differences are not observed in the linear Raman spectra.     

Molecular dynamics study of polarizability anisotropy relaxation in aromatic liquids and its connection with local structure

The collective polarizability anisotropy dynamics in a set of three aromatic liquids, benzene (Bz), hexafluorobenzene (HFB), and 1,3,5-trifluorobenzene (TFB), has been studied by molecular dynamics simulation. These liquids have very similar shapes, but different electrostatic interactions due to opposite polarities of C-H and C-F bonds, giving rise to different local intermolecular structures in the liquid phase. We have investigated how these structural arrangements affect polarizability anisotropy dynamics observed in optical Kerr-effect (OKE) spectroscopy. We have modeled the interaction-induced polarizability with the first-order dipole-induced dipole approximation, with the molecular polarizability distributed over the carbon sites. Local contributions to the librational OKE spectrum were computed separately for molecules participating in parallel or perpendicular relative orientations within the first coordination shell. We found that the relative locations of parallel and perpendicular librational bands of the OKE spectra are closely related to the corresponding pair energy distributions of the closest four neighbors of a given molecule, corresponding to a model of a harmonic oscillator in a cage of nearest neighbors. This model predicts higher librational frequencies for more attractive intermolecular interactions, which in all three liquids correspond to parallel local arrangements. On the diffusive orientational time scale, all three liquids exhibit slower relaxation of molecules in parallel arrangements, although the difference in relaxation rates is substantial only in TFB, which has the strongest tendency toward parallel stacking. The analysis of the collective polarizability relaxation was performed using two different approaches, the projection scheme (J. Chem. Phys. 1980, 72, 2801) and the theory developed by Steele (Mol. Phys. 1987, 61, 1031) for the second time derivatives applied to collective time correlations. Both approaches allow the decomposition of the OKE response into contributions from orientational relaxation and other dynamical processes. We find that they lead to different predictions on how the response depends on collective reorientation and processes arising from fluctuations in the interaction-induced polarizability. We discuss the reasons for these differences and the advantages and disadvantages of the two analysis schemes.     

Gaussian induced dipole polarization model

A new induced dipole polarization model based on interacting Gaussian charge densities is presented. In contrast to the original induced point dipole model, the Gaussian polarization model is capable of finite interactions at short distances. Aspects of convergence related to the Gaussian model will be explored. The Gaussian polarization model is compared with the damped Thole-induced dipole model and the point dipole model. It will be shown that the Gaussian polarization model performs slightly better than the Thole model in terms of fitting to molecular polarizability tensors. An advantage of the model based on Gaussian charge distribution is that it can be easily generalized to other multipole moments and provide effective damping for both permanent electrostatic and polarization models. Finally, a method of parameterizing polarizabilities is presented. This method is based on probing a molecule with point charges and fitting polarizabilities to electrostatic potential. In contrast to the generic atom type polarizabilities fit to molecular polarizability tensors, probed polarizabilities are significantly more accurate in terms of reproducing molecular polarizability tensors and electrostatic potential, while retaining conformational transferability.     

Geometrical derivatives of dipole moments and polarizabilities

Molecular dipole moments and polarizabilities, as well as their geometrical derivatives, are given analytical expressions for multiconfiguration self-consistent-field and configuration interaction wavefunctions. By considering the response of the electronic wavefunction induced by electric field and geometrical displacement terms in the Hamiltonian, the response of the total electronic energy to these terms is analyzed. The dipole moment and polarizability are then identified through the factors in the energy which are linear and quadratic in the electric field, respectively. Derivatives with respect to molecular deformation are obtained by identifying factors in these moments which are linear, quadratic, etc., in the distortion parameter. The analytical derivative expressions obtained here are compared to those which arise through finite-difference calculations, and it is shown how previous configuration-interaction-based finite difference dipole moment and polarizability derivatives are wrong. The proper means of treating such derivatives are detailed.     

Electron paramagnetic resonance line shifts and line shape changes due to heisenberg spin exchange and dipole-dipole interactions of nitroxide free radicals in liquids 8. Further experimental and theoretical efforts to separate the effects of the two interactions

The work in part 6 of this series (J. Phys. Chem. A 2009, 113, 4930), addressing the task of separating the effects of Heisenberg spin exchange (HSE) and dipole-dipole interactions (DD) on electron paramagnetic resonance (EPR) spectra of nitroxide spin probes in solution, is extended experimentally and theoretically. Comprehensive measurements of perdeuterated 2,2,6,6-tetramethyl-4-oxopiperidine-1-oxyl (pDT) in squalane, a viscous alkane, paying special attention to lower temperatures and lower concentrations, were carried out in an attempt to focus on DD, the lesser understood of the two interactions. Theoretically, the analysis has been extended to include the recent comprehensive treatment by Salikhov (Appl. Magn. Reson. 2010, 38, 237). In dilute solutions, both interactions (1) introduce a dispersion component, (2) broaden the lines, and (3) shift the lines. DD introduces a dispersion component proportional to the concentration and of opposite sign to that of HSE. Equations relating the EPR spectral parameters to the rate constants due to HSE and DD have been derived. By employing nonlinear least-squares fitting of theoretical spectra to a simple analytical function and the proposed equations, the contributions of the two interactions to items 1-3 may be quantified and compared with the same parameters obtained by fitting experimental spectra. This comparison supports the theory in its broad predictions; however, at low temperatures, the DD contribution to the experimental dispersion amplitude does not increase linearly with concentration. We are unable to deduce whether this discrepancy is due to inadequate analysis of the experimental data or an incomplete theory. A new key aspect of the more comprehensive theory is that there is enough information in the experimental spectra to find items 1-3 due to both interactions; however, in principle, appeal must be made to a model of molecular diffusion to separate the two. The permanent diffusion model is used to illustrate the separation in this work. In practice, because the effects of DD are dominated by HSE, negligible error is incurred by using the model-independent extreme DD limit of the spectral density functions, which means that DD and HSE may be separated without appealing to a particular model.     

A density functional theory-based chemical potential equalisation approach to molecular polarizability

The electron density changes in molecular systems in the presence of external electric fields are modeled for simplicity in terms of the induced charges and dipole moments at the individual atomic sites. A chemical potential equalisation scheme is proposed for the calculation of these quantities and hence the dipole polarizability within the framework of density functional theory based linear response theory. The resulting polarizability is expressed in terms of the contributions from individual atoms in the molecule. A few illustrative numerical calculations are shown to predict the molecular polarizabilities in good agreement with available results. The usefulness of the approach to the calculation of intermolecular interaction needed for computer simulation is highlighted.     

Infrared absorption by collisional CH4+X pairs, with X=He, H2, or N2

Existing measurements of the collision-induced rototranslational absorption spectra of gaseous mixtures of methane with helium, hydrogen, or nitrogen are compared to theoretical calculations, based on refined multipole-induced and dispersion force-induced dipole moments of the interacting molecular pairs CH4-He, CH4-H2, and CH4-N2. In each case the measured absorption exceeds the calculations substantially at most frequencies. We present the excess absorption spectra, that is the difference of the measured and the calculated profiles, of these supramolecular CH4-X systems at various gas temperatures. The excess absorption spectra of CH4-X pairs differ significantly for each choice of the collision partner X, but show common features (spectral intensities and shape) at frequencies from roughly 200 to 500 cm(-1). These excess spectra seem to defy modeling in terms of ad hoc exchange force-induced dipole components attempted earlier. We suggest that besides the dipole components induced by polarization in the electric molecular multipole fields and their gradients, and by exchange and dispersion forces, other dipole induction mechanisms exist in CH4-X complexes that presumably are related to collisional distortion of the CH4 molecular frame.     

Dynamic polarizabilities and raman intensities of CO,N2, HCl and Cl2

We have performed first-order (identical to coupled Hartree-Fock) and second-order polarization propagator calculations of the dynamic dipole polarizability tensor for the CO, N₂, HCl and Cl₂ molecules. The derivatives of the polarizability tensor with respect to the internuclear distance at the equilibrium internuclear separation are compared with related data obtained from non-resonance Raman spectra. In most cases the correlation contributions beyond the coupled Hartree-Fock approximation have a more pronounced effect on the derivatives of the polarizability tensor than on the polarizability tensor itself. At both levels of approximation we found that derivatives of the dynamic polarizability tensor with respect to the internuclear separation increase with 10–15% for variation of the frequency from 0 (static polarizability) to about 28500 cm^?1. The depolarization ratio calculated from the polarizability derivatives shows no variation with frequency in the same frequency range.     

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.     

Nonequilibrium molecular dynamics simulations with a backward-forward trajectories sampling for multidimensional infrared spectroscopy of molecular vibrational modes

A full molecular dynamics (MD) simulation approach to calculate multidimensional third-order infrared (IR) signals of molecular vibrational modes is proposed. Third-order IR spectroscopy involves three-time intervals between three excitation and one probe pulses. The nonequilibrium MD (NEMD) simulation allows us to calculate molecular dipoles from nonequilibrium MD trajectories for different pulse configurations and sequences. While the conventional NEMD approach utilizes MD trajectories started from the initial equilibrium state, our approach does from the intermediate state of the third-order optical process, which leads to the doorway-window decomposition of nonlinear response functions. The decomposition is made before the second pump excitation for a two-dimensional case of IR photon echo measurement, while it is made after the second pump excitation for a three-dimensional case of three-pulse IR photon echo measurement. We show that the three-dimensional IR signals are efficiently calculated by using the MD trajectories backward and forward in time for the doorway and window functions, respectively. We examined the capability of the present approach by evaluating the signals of two- and three-dimensional IR vibrational spectroscopies for liquid hydrogen fluoride. The calculated signals might be explained by anharmonic Brownian model with the linear-linear and square-linear system-bath couplings which was used to discuss the inhomogeneous broadening and dephasing mechanism of vibrational motions. The predicted intermolecular librational spectra clearly reveal the unusually narrow inhomogeneous linewidth due to the one-dimensional character of HF molecule and the strong hydrogen bond network.     

External electric field effects on absorption, fluorescence, and phosphorescence spectra of diphenylpolyynes in a polymer film

External electric field effects on absorption, fluorescence, and phosphorescence spectra of a series of unsubstituted diphenylpolyynes have been examined in a PMMA film. The analysis of the electroabsorption spectra indicates that the shorter diphenylpolyynes exhibit only the change in molecular polarizability, whereas the longer ones exhibit the change both in dipole moment and in molecular polarizability following absorption. The finding of the change in dipole moment following absorption of centrosymmetric diphenylpolyynes is interpreted in terms of the symmetry distortion upon doping a polymer film. When the external electric field is applied, the fluorescence yield is reduced and enhanced, respectively, in diphenylacetylene and diphenyloctatetrayne, indicating that the rate of the nonradiative process from the fluorescence state is accelerated in diphenylacetylene and decelerated in diphenyloctatetrayne by an external electric field. All of the diphenylpolyynes used in the present study exhibit the change in molecular polarizability following the phosphorescence process.     

Effective molecular polarizabilities and crystal refractive indices estimated from x-ray diffraction data

Although it was proposed some time ago that (hyper)polarizabilities might be estimated from the results of x-ray charge density refinements, early results were unconvincing. In this work we show that the one particle density obtained from the usual multipole refinement model does not contain sufficient information to determine these response properties and instead pursue the "constrained wave function" approach of fitting to x-ray structure factors. Simplified sum-over-states expressions are derived for determining the dipole polarizability from these wave functions, and these clearly show that the earlier work ignored important two-electron expectation values for the dipole polarizability, and two- and three-electron terms for Beta, etc. Correction factors for the simplified sum-over-states polarizability tensors from the constrained wave function are obtained by calibration against coupled Hartree-Fock ab initio results to yield in-crystal effective polarizability tensors. Results obtained for benzene, urea, and 2-methyl-4-nitroaniline demonstrate that the effective molecular polarizabilities clearly include the effects of intermolecular interactions and electron correlation, especially for urea where the effects on the polarizability are known to be quite large. We also carefully consider the way in which the linear bulk susceptibility, chi((1)), and refractive indices are determined from the x-ray fitted polarizabilities, employing three models based on a rigorous treatment of the local field. Incorrect results are obtained for the sort of molecules that are of interest in nonlinear optical applications if the molecules are approximated by single point dipoles. In contrast, the use of Lorentz-factor tensors averaged over several sites yields excellent results, with refractive indices obtained using this model in remarkably good agreement with optical measurements extrapolated to zero frequency.     

First-principles determinations and investigations of the electronic absorption and third-order polarizability spectra of electron donor-acceptor chromophores tetraalkylammonium halide/carbon tetrabromide

Calculations on donor-acceptor molecular pairs of tetraalkylammonium halide/carbon tetrabromide complexes are provided to investigate structure/property-related linear and nonlinear optical properties by using the time-dependent density functional theory technique coupled with the sum-over-states method. The calculated energies of the first allowed electronic transition decrease, and the nonresonant third-order polarizabilities at the THG, EFISHG, and DFWM optical processes increase progressively from [DBU-H+Br-.CBr(4) to [NPr(4)Br.CBr(4)] to [NMe(4)Br.CBr(4)]. The obtained electronic absorption spectra show a progressive red shift with increasing donor strength from Cl to I for [NR(4)h.CBr(4)] (h = Cl, Br, and I). The charge transfers from the halogen donor to the carbon tetrabromide acceptor make significant contributions to the electronic absorption spectra in the low-energy zone and the third-order polarizabilities in the nonresonant frequency region. The counterion indirectly affects the electronic absorption and third-order polarizability spectra through the interactions between the donor and acceptor.     

Excited state dipole moments and polarizabilities of centrosymmetric and dimeric molecules. II. Polyenes, polyynes and cumulenes

Electro-optical absorption spectra are measured for a series of polyenes, polyynes and cumulenes with centrosymmetric π-chromophores in cyclohexane solution at 298 K. For all molecules the long-axis component of the polarizability tensor is considerably larger in the first dipole-allowed singlet state compared to the ground state. The transition moments are found to be parallel to the long molecular axis. All polyenes and one cumulene show a linear Stark component indicating a long-axis excited state dipole moment. Both the dipole moments and the polarizabilities are corrected within the extended Onsager model for solvent cavity and reaction field effects. It is suggested that symmetry lowering solvent perturbations are the reason for the apparent excited state dipole moments.     

Polarizable contributions to the surface tension of liquid water

Surface tension, gamma, strongly affects interfacial properties in fluids. The degree to which polarizability affects gamma in water is thus far not well established. To address this situation, we carry out molecular dynamics simulations to study the interfacial forces acting on a slab of liquid water surrounded by vacuum using the Gaussian charge polarizable (GCP) model at 298.15 K. The GCP model incorporates both a fixed dipole due to Gaussian distributed charges and a polarizable dipole. We find a well-defined bulklike region forms with a width of approximately 31 A. The average density of the bulklike region agrees with the experimental value of 0.997 g/cm3. However, we find that the orientation of the molecules in the bulklike region is strongly influenced by the interfaces, even at a distance five molecular diameters from the interface. Specifically, the orientations of both the permanent and induced dipoles show a preferred orientation parallel to the interface. Near the interface, the preferred orientation of the dipoles becomes more pronounced and the average magnitude of the induced dipoles decreases monotonically. To quantify the degree to which molecular orientation affects gamma, we calculate the contributions to gamma from permanent dipolar interactions, induced dipolar interactions, and dispersion forces. We find that the induced dipole interactions and the permanent dipole interactions, as well as the cross interactions, have positive contributions to gamma, and therefore contribute stability to the interface. The repulsive core interactions result in a negative contribution to gamma, which nearly cancels the positive contributions from the dipoles. The large negative core contributions to gamma are the result of small oxygen-oxygen separation between molecules. These small separations occur due to the strong attractions between hydrogen and oxygen atoms. The final predicted value for gamma (68.65 m/Nm) shows a deviation of approximately 4% of the experimental value of 71.972 m/Nm. The inclusion of polarization is critical for this model to produce an accurate value.