Collective molecular rotation in D2O

Liquid Water: When Hyperpolarizability Fluctuations Boost and Reshape the Second Harmonic Scattering Intensities

Second harmonic scattering (SHS) is a method of choice to investigate the molecular structure of liquids. While a clear interpretation of SHS intensity exists for diluted solutions of dyes, the scattering due to solvents remains difficult to interpret quantitatively. Here, we report a quantum mechanics/molecular mechanics (QM/MM) approach to model the polarization-resolved SHS intensity of liquid water, quantifying different contributions to the signal. We point out that the molecular hyperpolarizability fluctuations and correlations cannot be neglected. The intermolecular orientational and hyperpolarizability correlations up to the third solvation layer strongly increase the scattering intensities and modulate the polarization-resolved oscillation that is predicted here by QM/MM without fitting parameters. Our approach can be generalized to other pure liquids to provide a quantitative interpretation of SHS intensities in terms of short-range molecular ordering.

Comment on "Water-water correlations in electrolyte solutions probed by hyper-Rayleigh scattering" [J. Chem. Phys. 147, 214505 (2017)]

The work by Shelton [J. Chem. Phys. 147, 214505 (2017)] discussed and interpreted differences with a previous study by Chen et al. [Sci. Adv. 2, e1501891 (2016)] regarding the influence of electrolytes on the structure of water. It is argued by Shelton [J. Chem. Phys. 147, 214505 (2017)] that impurities and hyper-Raman scattering contributions are the reasons for differences in the measured second harmonic intensity between the above two studies. Here, we show that these proposed effects are not relevant and discuss the influence of pulse parameters, focusing on pulse duration, since these two sets of experiments are performed with substantially different pulse durations, 100 ns and 190 fs, respectively. We show that inelastic higher-order effects play a role in the experiment with 100 ns laser pulses (the probed structure is that of the electrolyte solution that is modified by a laser pulse), while in the experiment with 190 fs laser pulses, only the elastic second-order response is measured (probing the unperturbed water structure).

Second-Harmonic Scattering as a Probe of Structural Correlations in Liquids

Second-harmonic scattering experiments of water and other bulk molecular liquids have long been assumed to be insensitive to interactions between the molecules. The measured intensity is generally thought to arise from incoherent scattering due to individual molecules. We introduce a method to compute the second-harmonic scattering pattern of molecular liquids directly from atomistic computer simulations, which takes into account the coherent terms. We apply this approach to large-scale molecular dynamics simulations of liquid water, where we show that nanosecond second-harmonic scattering experiments contain a coherent contribution arising from radial and angular correlations on a length scale of ≲1 nm, much shorter than had been recently hypothesized ( Shelton , D. P. J. Chem. Phys. 2014 , 141 ). By combining structural correlations from simulations with experimental data ( Shelton , D. P. J. Chem. Phys. 2014 , 141 ), we can also extract an effective molecular hyperpolarizability in the liquid phase. This work demonstrates that second-harmonic scattering experiments and atomistic simulations can be used in synergy to investigate the structure of complex liquids, solutions, and biomembranes, including the intrinsic intermolecular correlations.

Formation of Ferroelectric Domains Observed in Simulation of Droplets of Dipolar Particles

In this work it is shown that domains of ordered dipoles are formed in large droplets made from dipolar particles provided that the dipole-dipole interaction between nearest neighbors is larger than the thermal energy. The size of the domains grows almost linearly with the size of the droplets for droplets containing 1000-30 000 particles. The largest domains occupy around 25-35% of the droplet volume. The total dipole moment of a domain is of the order of 3-10% of the maximum dipole moment possible if all dipoles in the domain were parallel. The finding offers an explanation to the observation that different boundary conditions yield different long-range order for dipolar liquids and challenges the present view of a short-range dipolar order in polar solvents.

Are dipolar liquids ferroelectric? Simulation studies

The observation of a very sharp low frequency spike in the hyper-Rayleigh spectrum (HRS) of strongly dipolar fluids, such as acetonitrile and water, has been interpreted as reflecting a very slowly relaxing component in the transverse dipole density. This suggestion is at variance with the expectation of the dielectric theory for an isotropic fluid and has led to the speculation that the slow relaxation is associated with the reorganization of ferroelectric domains. Very large-scale molecular-dynamics simulation (approximately 28,000 molecules) have been carried out using a three-site potential model of acetonitrile. The simulated fluid shows no suggestion of strong dipole correlations and domain structure. The dipole density correlations behave as predicted by the normal dielectric theory and their spectra do not show the low-frequency feature seen in the HRS. In order to examine the characteristics of the spectra which would be seen in a ferroelectric domain, the acetontrile model was transmuted to more closely resemble a Stockmayer-type fluid with the same dipole density and a ferroelectric phase was observed. In this phase the dielectric spectra show (i) a high-frequency spectral feature due to librational motion of the molecules within a domain, and (ii) slowly-relaxing longitudinal and transverse polar modes, again at variance from the experimental HRS characteristics.

Ferroelectric domains in nitrobenzene-nitromethane solutions measured by hyper-Rayleigh scattering

Hyper-Rayleigh scattering (HRS) spectra were measured for liquid solutions of C6H5NO2 and CH3NO2 at T=300 K. The depolarized HRS spectra at small frequency shift are dominated by two components due to reorientation of the nitrobenzene molecules. One is a Lorentzian with spectral width nu1=0.16-0.45 cm(-1) and corresponding orientation relaxation time tau=33-12 ps. The second component is a narrow spike with spectral width <2 MHz and corresponding relaxation time tau>80 ns, attributed to HRS from slowly relaxing ferroelectric domains. The dipole order parameter g0=0.053+/-0.005, saturation parameter p=0.9+/-0.1, and volume V=20+/-6 nm3 for these domains in nitromethane were determined from measurements of the nitrobenzene-concentration dependence of the intensity ratio for these two spectral components. Orientation of the 230 nitromethane molecules within each domain is inhomogenous but highly ordered.

Hyper-Rayleigh scattering spectrum of liquid nitromethane

The vertical vertical (VV), horizontal vertical (HV), and vertical horizontal (VH) hyper-Rayleigh scattering (HRS) spectra were measured for liquid CH3NO2 at T= 300 K. The main HRS spectral component has a width upsilon1=1.28 +/- 0.04 cm(-1), which gives an orientation relaxation time pi=4.1 +/- 0.1 ps in good agreement with other experiments. However, the VH spectrum also contains a previously unobserved strong narrow peak at zero-frequency shift, absent from the VV and HV spectra, which is due to a slowly relaxing longitudinal orientation mode. The upper bound on the width of this peak is 5 MHz, which corresponds to a relaxation time pi > 30 ns.

Are dipolar liquids ferroelectric?

VH and HV depolarized hyper-Rayleigh scattering spectra were measured for liquid solutions of dipolar CH3CN in nondipolar C2Cl4 at T=300 K. The VH spectrum contains a strong narrow peak due to a slowly relaxing longitudinal orientation mode. This peak is absent in the HV spectrum, and it disappears from the VH spectrum when the CH3CN concentration is reduced to 8%. This observation is consistent with a ferroelectric phase transition predicted to occur when rho mu0(2)=9epsilon0kT=49 D2 M.