Arxiv:1611.03384V1

Second Harmonic Scattering as a Probe of Structural Correlations in Liquids Gabriele Tocci,1, 2, ∗ Chungwen Liang,1, 2 David M. Wilkins,1, 2 Sylvie Roke,1 and Michele Ceriotti2 1Laboratory for fundamental BioPhotonics, Institutes of Bioengineering and Materials Science and Engineering, School of Engineering, and Lausanne Centre for Ultrafast Science, Ecole´ Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland 2Laboratory of Computational Science and Modeling, Institute of Materials, Ecole´ Polytechnique Fédérale de Lausanne, 1015 Lausanne, Switzerland 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 [1]. By combining structural correlations from simulations with experimental data [1] 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 inter-molecular correlations. Nonlinear light scattering has been widely used to order response arising from collective polar modes [21], investigate aqueous interfaces, including suspensions of collective rotation of molecules [22], the presence of fer- metallic or semiconducting nanoparticles, water droplets roelectric domains [23], coupling of rotations and trans- and biological membranes (for example, Refs. [2{15]). In lations in acoustic phonons extending up to 2000 nm [1], particular, second harmonic scattering (SHS), also com- or orientational correlations found in ideal isotropic ho- monly referred to as Hyper-Rayleigh scattering (HRS; mogeneous random vector fields [24]. Clearly SHS mea- the terms SHS and HRS are used interchangeably in this surements contain information about the intermolecular work), is especially sensitive to the molecular orientation structure of liquids although it is evident that no consen- of liquids at interfaces: there is a greater contribution to sus has been reached on their origin (see e.g. Ref. [25]) radiated second harmonic light arising from molecules at and their length scale. So far, experiments in bulk wa- the interface, where inversion symmetry is broken, com- ter have been interpreted by treating intermolecular cor- pared to those in the bulk, which is on average centrosym- relations as fitting parameters [1, 24] and a framework metric. In bulk liquids the coherent term is usually as- to explicitly account for the correlations that give rise sumed to be negligible and the measured signal is mostly to coherent scattering in SHS experiments has not yet interpreted in terms of fully incoherent scattering from been developed. Recent elastic femtosecond SHS mea- individual molecules [16, 17]. Under this assumption, the surements from ionic solutions have shown that several SHS intensity of a solution can be expressed as a linear different electrolytes induce long-range orientational cor- combination of the incoherent contribution of the solute relations in water, starting at concentrations as low as 10 and of that of the solvent. Using this approach the aver- µM [26]. The strong dependence on isotopic composition age hyperpolarizability of a molecule in solution can be also suggests a link between these observations and the estimated [18]. H-bond network of water. This assumption stems from early work by Terhune The liquid-phase hyperpolarizability plays a crucial et al., who neglected the coherent contribution to the role in the interpretation of SHS experiments. Mea- arXiv:1611.03384v1 [cond-mat.soft] 10 Nov 2016 intensity because the wavelength of the laser is much surements have often been analyzed using a molecular larger than the range of molecular correlations expected non-linear response tensor of water obtained from quan- in a liquid [16]. Although Bersohn and Maker later de- tum chemical calculations [27]. However, the model that veloped a theory to include coherent scattering due to underlies the evaluation of of the hyperpolarizability in structural correlations [19, 20], these additional terms Ref. [27] described the surroundings of a water molecule have been largely ignored in the interpretation of exper- by three point charges, which does not take into account iments. More recently, in a series of articles by Shel- the many complex molecular environments that can be ton et al., a sign of coherence has been reported from found in liquid water [28]. nanosecond HRS measurements on bulk liquids includ- In this work we present a method to compute the SHS ing water and acetonitrile [1, 21{24]. A diverse range of intensity directly from atomistic computer simulations, possible origins was suggested for this coherence: a third including both the incoherent and the coherent contri- 2 butions. This framework is used to calculate the angular to zero if the instantaneous relative orientation of dif- scattering pattern from large-scale force-field molecular ferent molecules were completely random, and is often dynamics (MD) simulations of liquid water. We show assumed to be negligible compared to the incoherent that quantitative agreement with nanosecond HRS ex- term [2, 16, 18, 20] In order to test this assumption, we periments [1] cannot be achieved by using the values developed a framework to explicitly evaluate Eq. 1 for of molecular hyperpolarizability often used to interpret MD simulations. We applied this framework to MD simu- experiments, [27] obtained with Møller-Plesset pertur- lations of liquid water involving about 260,000 molecules, bation theory expanded up to the fourth order (MP4). using a cubic simulation box with a side of 20 nm. Simu- We find instead that it is possible to obtain satisfactory lations were performed with the GROMACS code v.5.0.4 agreement by combining the correlations obtained from [29], using the rigid TIP4P/2005 water model [30]. The simulations with an effective liquid phase molecular hy- equations of motion were integrated in the NVT ensem- perpolarizability, used as a fitting parameter. The values ble using the velocity Verlet algorithm for 20 ns with a obtained for this hyperpolarizability can be used to pro- 2 fs timestep. Temperature control was achieved using vide a more quantitative interpretation of other SHS ex- the stochastic velocity rescaling thermostat [31], with a periments. By computing explicitly the correlations that target temperature of 300 K. Full details of the computa- contribute to the SHS signal, we can also gather insight tional setup can be found in the Supporting Information. into the length scales that are most relevant for these We performed extensive tests on the setup to ensure that experiments. it was insensitive to finite size effects and the type of We begin by introducing the general expression for the force-field used. SHS intensity I(q; u; v) for an ensemble of N molecules Fig. 1 shows the ratios of the SHS intensities of bulk at a scattered wave-vector q and with the polarization water reported by Shelton in Ref. [1] using nanosecond directions of the incoming and outgoing beams defined experiments at different polarization combinations as a by the vectors u and v: function of the scattering angle (θ) and as computed in this work from the analysis of our MD simulations. The * 2+ X convention used for the polarization combinations follows I(q; u; v) / N −1 β(α; u; v)eiq·rα ; (1) Ref. [26] and references therein. Light polarized parallel α to the scattering plane is denoted by P and light polarized where the brackets indicate a time average over uncor- perpendicular to the scattering plane is denoted by S. A related configurations and β(α; u; v) is the component polarization combination is specified as XYY, with X (= of the hyperpolarizability of a scattering unit α in the P or S) the polarization of the outgoing beam and Y (= laboratory reference frame (L) projected onto the polar- P or S) that of the incoming beam. Fig. 1(a) illustrates ization direction of the incoming (u) and outcoming (v) the intensity ratios obtained from our simulations un- P L beams, i.e. β(α; u; v) = LMN βLMN (α)vLuM uN , and der the assumption of purely incoherent scattering. We the proportionality constant can be found in Ref. [19]. A use the constant hyperpolarizability from Ref. [27], com- schematic of the geometric setup is shown in the Support- puted using MP4 and modelling a liquid-like environment ing Information (SI) in Fig. S1. Eq. 1 can be decomposed by placing point charges in a geometry resembling the into the sum of an incoherent term and of a coherent term first solvation shell. Despite its simplicity, this model respectively [19]: has been used in several studies of the hyperpolarizability of water [32, 33]. It can be seen that a model for the * + X intensity assuming only incoherent scattering completely I(q; u; v) / N −1 β(α; u; v)2 + fails to account for the qualitative features of the exper- α * + (2) imental curves, which show an asymmetry with respect −1 X X 0 [iq·(r −r 0 )] to 90 degrees and a varying I =I ratio. N β(α; u; v)β(α ; u; v)e α α : SPP PSS α α06=α We then proceed to compute the intensity with the M full Eq. 1, using the βabc elements from Ref. [27]. As L The hyperpolarizability in the laboratory frame βLMN shown in Fig.

Load more