research papers
Anisotropy enhanced X-ray scattering from solvated transition metal complexes
aDepartment of Physics, Technical University of Denmark, Fysikvej 307, DK-2800 Kongens Lyngby, Denmark, bPULSE Institute, SLAC National Accelerator Laboratory, Menlo Park, CA 94025, USA, cDepartment of Chemistry, Technical University of Denmark, Kemitorvet 207, DK-2800 Kongens Lyngby, Denmark, dFaculty of Physical Sciences, University of Iceland, Reykjavik, Iceland, eDepartment of Chemical Physics, Lund University, Box 118, S-22100 Lund, Sweden, and fLinac Coherent Light Source, SLAC National Accelerator Laboratory, Menlo Park, CA 94025, USA
*Correspondence e-mail: elbia@fysik.dtu.dk
Time-resolved X-ray scattering patterns from photoexcited molecules in solution are in many cases anisotropic at the ultrafast time scales accessible at X-ray free-electron lasers (XFELs). This anisotropy arises from the interaction of a linearly polarized UV–Vis pump laser pulse with the sample, which induces anisotropic structural changes that can be captured by femtosecond X-ray pulses. In this work, a method for quantitative analysis of the anisotropic scattering signal arising from an ensemble of molecules is described, and it is demonstrated how its use can enhance the structural sensitivity of the time-resolved X-ray scattering experiment. This method is applied on time-resolved X-ray scattering patterns measured upon i.e. the change in Pt—Pt bond length, and yields more information on the excitation channels than the analysis of the isotropic scattering only. Finally, it is discussed how the anisotropic transient response of the solvent can enable the determination of key experimental parameters such as the instrument response function.
of a solvated di-platinum complex at an XFEL, and the key parameters involved are explored. It is shown that a combined analysis of the anisotropic and isotropic difference scattering signals in this experiment allows a more precise determination of the main photoinduced structural change in the solute,Keywords: time-resolved; anisotropic scattering; orientational selection; XFEL; ultrafast; molecular structure.
1. Introduction
Time-resolved X-ray diffuse scattering (XDS) experiments give insight into the photoinduced structural dynamics of solvated molecules. In these experiments, a laser pulse initiates the dynamic process, which is subsequently probed by an X-ray probe pulse arriving at specific time delays after the pump event. If the laser pulse is ultra-short, the ensuing structural dynamics are coherently initiated in the molecular ensemble (Zewail, 2000), and the scattering signal can be used to retrieve structural changes occurring in the molecule after the electronic excitation (Borfecchia et al., 2013; Biasin et al., 2016; Haldrup et al., 2012; Kim et al., 2015a; van Driel et al., 2016; Chergui & Collet, 2017; Kong et al., 2008). Such experiments are of fundamental importance for understanding the structure–function relationship of, for instance, transition metal complexes whose photochemical and photophysical properties can be applied in technologies such as solar energy conversion and (Takeda et al., 2017; Esswein & Nocera, 2007; White, 1982).
In typical time-resolved XDS experiments, laser pump and X-ray probe pulses are focused onto a thin liquid jet, which is produced by pumping the liquid sample through a nozzle. The diffuse X-ray scattering is collected on a two-dimensional (2D) detector placed after the sample on the plane perpendicular to the propagation direction of the X-ray beam (Fig. 1). The measured scattering signal contains information about all the inter-nuclear distances of the sample at a specific time delay and it is usually dominated by the scattering from the solvent. The sensitivity to the solute is enhanced in the difference scattering signal, which is constructed by subtracting the signal collected without photoexciting the sample from the signal collected after In this way the unchanging background contributions cancel out, and the difference scattering signal arises from the changes in the inter-atomic distances in the probed sample volume (Borfecchia et al., 2013; Ihee et al., 2010). The established procedure for analysing scattering data from liquid samples consists of an azimuthal integration of the 2D difference scattering patterns, since the signal is usually assumed to arise from an isotropically distributed ensemble of molecules (Ihee et al., 2010; Haldrup et al., 2010). If the sample is isotropic, the azimuthal integration allows for an improvement of the signal-to-noise (S/N) without loss of information. This procedure was first established with synchrotron data and it is justified when the time resolution is longer than the molecular rotational correlation time in solution, which ranges in the 10–100 ps time scale for transition metal complexes (Lakowicz, 2006; Kim et al., 2015b).
At X-ray free-electron lasers (XFELs), that can deliver femtosecond X-ray pulses, the photoinduced structural changes in the sample can be captured at the time scale of atomic motions. This enables the observation of vibrational and rotational et al., 2016; Kim et al., 2015b; Lemke et al., 2017; Chergui & Collet, 2017). If the photoselection process (illustrated in Fig. 2 and further described below) creates an aligned excited-state ensemble of molecules, anisotropic 2D scattering patterns can be observed on time scales shorter than the rotational correlation time of the molecules (Kim et al., 2011; Kim et al., 2015b; Yang et al., 2016a,b; Glownia et al., 2016). Anisotropic scattering patterns contain information on the molecular and thus can add sensitivity to spatial compared with isotropic patterns (Hensley et al., 2012; Burger et al., 2010; Küpper et al., 2014). The theoretical foundation for interpreting anisotropic scattering contributions from aligned ensembles of molecules in time-resolved experiments has been laid out in some detail over the past ten years (Baskin & Zewail, 2006; Lorenz et al., 2010; Brinkmann & Hub, 2015; Penfold et al., 2012). However, quantitative structural analysis of anisotropic scattering data from molecular ensembles where a photoselected sub-population has been promoted to an electronic has been demonstrated only in a very few cases in the gas phase (Yang et al., 2016a,b; Glownia et al., 2016), and a robust methodology for such cannot yet be said to have been fully established as the debate surrounding the work of Glownia et al. shows (Bennett et al., 2017; Glownia et al., 2017). Specifically with respect to solution-state molecular systems, only qualitative analysis of time-resolved anisotropic scattering data have been reported, and the anisotropic contributions to the observed scattering patterns have been mostly used only to obtain the rotational correlation time of the molecules in solution (Kim et al., 2011; Kim et al., 2015b).
in, for instance, solvated transition metal complexes (BiasinIn this article, we revisit the formalism derived by Lorenz et al. (2010) for analysing scattering from a photo-aligned ensemble of molecules. We use this formalism to extract the isotropic and anisotropic contribution from 2D scattering patterns measured at an XFEL upon of a transition metal complex in solution and we explore a few of the key parameters involved in this procedure. We show that, in the presence of anisotropic scattering, the separation of the two contributions can, and should, replace the standard azimuthal integration in the data reduction procedures. As a key point, we demonstrate that the anisotropic contribution can be quantitatively analysed in order to extract structural information about the photoexcited ensemble of molecules. These results are discussed in the framework of how the information extracted from a combined analysis of the isotropic and anisotropic contributions to the measured difference scattering signals can potentially help disentangle the many inter- and intra-nuclear involved in photoinduced structural dynamics of molecules in solution.
The method is exemplified on scattering data collected at the X-ray Pump Probe (XPP) instrument (Chollet et al., 2015) at the Linac Coherent Light Source upon of the tetrakis-μ-pyrophosphitodiplatinate(II) ion [Pt2(P2O5H2)4]4-, here abbreviated as PtPOP. Fig. 3 shows the molecular structure of PtPOP: a Pt–Pt dimer is held together by four pyrophosphito ligands. The compound belongs to the C4h with approximately fourfold symmetry along the Pt–Pt axis (Zipp, 1988; Gray et al., 2017). Thus, the Pt–Pt axis has higher rotational symmetry than the other two axes, and the molecule is a symmetric top (i.e. two have the same value and the third has a unique value). This compound was first synthesized in 1977 and has been the subject of a vast amount of studies, due to its characteristic photophysical properties and activity as a photocatalyst (Roundhill et al., 1989; Stiegman et al., 1987). In the last decades, PtPOP has become a model compound for time-resolved X-ray studies, due to its high scattering power, its long singlet excited-state lifetime, and its high symmetry (Christensen et al., 2009; van der Veen et al., 2009). For this experiment, an 80 mM aqueous solution of PtPOP was circulated through a nozzle producing a 50 µm round liquid jet. Each photocycle was initiated by a 3 µJ laser pulse at 395 nm and with 50 fs pulse width (FWHM), focused onto a 150 µm-diameter spot. The scattering from the 9.5 keV X-ray probe pulses was detected by the Cornell–SLAC Pixel Array Detector (CSPAD) (Philipp et al., 2011) covering scattering vectors up to 4.5 Å−1. The time delay t between the laser and the X-ray pulses was determined for every pump–probe event with ∼10 fs (FWHM) resolution using the XPP timing-tool (Minitti et al., 2015). The detector signal was corrected according to the procedure described by van Driel et al. (2015), including corrections for the detector geometry and the horizontal polarization of the X-ray beam. The corrected scattering signal was scaled to the liquid reflecting the stoichiometry of the sample (Haldrup et al., 2010), yielding the acquired signal in electron units per solute molecule (e.u. molec−1). Individual 2D difference scattering patterns were then time-sorted and averaged. Movie S1 of the supporting information shows averaged measured difference scattering patterns upon of PtPOP as a function of increasing time delay: a strongly anisotropic difference scattering signal is visible at an early time delay and decays on a ∼100 ps time scale. The signal is expected to arise predominantly from the shortening of the Pt–Pt distance since the 395 nm transition, which has polarization along the Pt—Pt bond, involves the promotion of an electron from the antibonding to the bonding pσ orbital (Stiegman et al., 1987). Christensen et al. investigated the structure of the of PtPOP in water with time-resolved XDS and obtained a Pt—Pt contraction of 0.24 ± 0.06 Å with respect to the ground state of the molecule (Christensen et al., 2009). van der Veen et al. concluded a 0.31 ± 0.05 Å Pt—Pt contraction upon excitation of the complex in ethanol, with X-ray absorption spectroscopy (van der Veen et al., 2009). The well defined transition and the strong scattering signal arising from the contraction of the Pt—Pt bond makes PtPOP an ideal model system to benchmark the formalism required in order to interpret anisotropic difference scattering patterns.
2. Methods
In time-resolved XDS experiments, the molecules in solution are randomly oriented before the arrival of the laser pump. If the linearly polarized ultrashort laser pulse interacts with a single transition μ of the molecule, the angular distribution D of this transition in molecules being excited displays cylindrical symmetry with respect to the laser polarization E. Therefore, this distribution can be expressed through an expansion in Legendre polynomials of only the angle α between the laser polarization axis and the transition (Baskin & Zewail, 2006),
where α and η are the polar and azimuthal angles, respectively, in a spherical coordinates system with polar axis E, and cn are expansion coefficients. In the case of one-photon absorption of linearly polarized radiation by an isotropic molecular ensemble in thermal equilibrium, the distribution of the transition in molecules being excited is a cosine-squared distribution with respect to the laser polarization axis (Møller & Henriksen, 2012; Baskin & Zewail, 2006). Cast into the generic form of equation (1), the distribution of excited-state molecules immediately after the excitation event can be expressed as a sum of two terms,
where P0(x) = 1 and P2(x) = (3x2 - 1)/2, with c2 = 2c0. Properly normalizing (Kleef & Powis, 1999) the angular distribution in equation (1) on unit sphere gives
where the last equality is a consequence of the orthogonality of the Legendre polynomials.
Assuming ultrashort laser pump and X-ray probe pulses, the differential scattering t, can be expressed by weighting the squared molecular form factor with the instantaneous distribution ρ of nuclear geometries (Lorenz et al., 2010),
from an ensemble of excited molecules, at a specific pump–probe delay timewhere has dimensions 3N and describes the nuclear coordinates in the laboratory frame, is the Thomson scattering and is the scattering vector. After assuming that rotational and vibrational of the molecule are uncoupled (Lorenz et al., 2010), we express the rotational part of the distribution through in equation (1) and the vibrational distribution as a δ-distribution. Under these assumptions, the scattering from an ensemble of molecules with their transition dipole moments aligned with respect to the laser polarization axis can be calculated through
where is the vector of inter-nuclear distances. Furthermore, this expression [equation (5)] can be simplified by utilizing the fact that the molecule is a symmetric top with the transition parallel to the unique axis of symmetry. Under this assumption, as detailed by Baskin & Zewail (2005), the angular distribution of each intra-molecular distance with respect to the laser polarization axis can be expressed through the angle between and the main axis of symmetry of the molecule. These parameters are illustrated in Fig. 3: PtPOP is a symmetric top molecule with the main axis of symmetry parallel to the Pt—Pt bond, which is also the direction of the transition photoselected in this experiment and therefore meets the above assumptions. As described by Lorenz et al. (2010), under the described assumptions, the integral in equation (5) has the analytical form
where θQ is the angle between the laser polarization axis and Q (see Fig. 3) and
where jn are the spherical Bessel functions and fi is the atomic form factor of atom i, which is used to express the molecular form factor within the independent atom model (Møller & Henriksen, 2012). Assuming one-photon absorption [and thus the angular distribution in equation (2)], only the n = 0 and n = 2 terms contribute in equation (6) (Lorenz et al., 2010) and the scattering signal can be written as
Since j0(x) = (sinx)/x, S0 is recognized as the Debye formula for isotropic ensembles,
while S2 contains information about the orientation of the single bond with respect to the transition of the molecule,
with j2(x) = [(3/x2) - 1][(sinx)/x] - 3[(cosx)/x2]. From here, S0 and S2 will be referred to as the isotropic and anisotropic part of the scattering signal, respectively. c2 decays from its initial value of 2 to 0 according to the rotational correlation time of the molecules in solution.
The excitation will leave a hole in the ground-state distribution of the same rotational anisotropy as the excited-state ensemble (Baskin & Zewail, 2006; Jonas et al., 1995). Starting from equation (8), the difference scattering signal including both excited-state and ground-state contributions can, therefore, be decomposed as (Lorenz et al., 2010)
where and are, respectively, the isotropic and anisotropic difference scattering signals arising from the photoinduced structural changes in the sample at a specific pump–probe time delay t.
With respect to equation (11), S0 and S2 depend only on internal coordinates of the molecule, while depends only on the geometry of the experiment. Specifically, as detailed by Baskin & Zewail (2006), can be expressed as a function of the angles δ, θ and φ, as defined in Fig. 1,
In the case of most standard XDS experiments, and for the PtPOP experiment presented here, the laser beam is nearly collinear to the X-ray beam. This is one of the possible configurations that yields a 90° angle between the laser polarization axis and the direction of propagation of the X-ray beam (δ in Fig. 1), and equation (12) simplifies to
Since φ is the azimuthal angle on the detector surface, maps P2 anisotropically onto the detector surface. The second-order Legendre polynomial, and hence the anisotropic contribution to the scattering, vanishes when in equation (12) is equal to the magic angle. However, as extensively described by Baskin & Zewail (2006), it is not possible to simultaneously remove anisotropic contributions from the entire scattering pattern for a given choice of δ. The linear relation between scattering intensity and P2 in equation (11) can be used to retrieve, for a specific Q value, the isotropic and anisotropic difference scattering signal, as described in the following.
Fig. 4 illustrates the extraction of the isotropic and anisotropic difference scattering signals from a difference scattering pattern measured 4.5 ps after the of PtPOP in water. The scattering pattern is shown in Fig. 4(b) and it is an average of 50 difference scattering patterns collected in a ∼40 fs-wide time bin centred at 4.5 ps after the arrival of the laser pump pulse. As introduced above, such a signal arises mainly from the photoinduced contraction of the Pt atoms along their connecting vector. In the excited-state ensemble, this vector has a cosine-squared distribution with respect to the laser polarization axis, since the transition is parallel to the Pt–Pt axis. Fig. 4(a) shows P2 [P2 = and as in equation (13)] mapped onto the CSPAD, with φ = 0 along the direction of the projection of the laser polarization axis on the detector surface, which was found to be inclined by 20° with respect to the vertical for the experiment described here.
For the analysis, the scattering pattern in Fig. 4(b) was divided into 500 radial bins and 45 azimuthal bins, and the signal in each bin was calculated as the average value of the pixels. The difference scattering signal corresponding to the radial bin centred at Q = 2 Å−1 (black circle) is plotted in Fig. 4(d) as a function of at the same Q value. The red line is a least-squares fit of a straight line to the data points. This fit, according to equation (11), yields (Q = 2 Å−1) as the slope and (Q = 2 Å−1) as the intercept with the P2 = 0 axis. Repeating the procedure for all the radial bins yields the one-dimensional isotropic and anisotropic difference scattering curves for the full Q range at this specific time delay, as Fig. 5 shows.
3. Results and discussion
Fig. 5 shows the isotropic and anisotropic difference scattering signals extracted from the scattering pattern measured at 4.5 ps as a function of Q. The uncertainty estimates for and at each Q point are calculated from the covariance matrix of the coefficients of the straight-line fit. Since the difference scattering signal is usually of the order of 0.1–1% of the total scattering signal [see comparison between Fig. S1 of the supporting information and Fig. 4(b) for this experiment], averaging several difference scattering patterns acquired at the same nominal time delay allows the improvement of the S/N ratio before the extraction of the isotropic and anisotropic contributions. Experimental parameters such as the magnitude of the differential scattering and the X-ray as well as the number of radial bins considered for each 2D pattern, determine the optimum number of images to be averaged for a specific experiment. The number of radial bins should be higher than the number of independent data points in the Q range [∼tens of points (Haldrup et al., 2010)] and can be chosen arbitrarily high since re-binning can be done at a later stage of the analysis. For the present experiment, we set the number of radial bins to 500 and, for each time bin, we choose the number of images to be averaged by inspection of the decreasing uncertainties σ on the coefficients of the straight-line fit [ and in equation (11)] as a function of number of averaged images. This is exemplified in Fig. 6(a) for the signal extracted at 4.5 ps: up to 50 averaged images, decreases as the inverse square root of the number of averaged images, as expected for Gaussian noise; after 50 images, converges to a constant value (∼4 e.u. molec−1) since the noise is then dominated by systematic errors due to contributions from a non-constant background and the non-linear response of the detector rather than counting statistics (van Driel et al., 2015). Therefore averaging more than 50 images will not further decrease the uncertainties on the measured signal.
Fig. 6(b) shows the number of unique P2 values that are sampled for a specific Q value as a function of the number of azimuthal bins into which the 2D diffraction patterns is divided. The plot shows that choosing an odd number of azimuthal bins maximizes the number of unique points available for the fitting of versus P2. This behaviour arises from the dependency of P2 on the azimuthal angle φ, as shown in Fig. 4(c), and from choosing equi-angular azimuthal bins. Moreover, since P2 is non-linear as a function of φ, we introduce an error by calculating P2 as the central value rather than the mean value in each azimuthal bin. This error can cause a significant decrease of the magnitude of the anisotropic difference scattering signal when decreasing the number of azimuthal bins, as shown in Fig. 7(a). In order to calculate the mean value of P2, we should use the following expression,
where is half the size of the azimuthal bins. Fig. 7(b) shows the difference between the values of P2 (red dots) and of the magnitude of (blue dots) obtained when calculating P2 either as the central or the mean value of the azimuthal bins. From this observation, if the 2D scattering pattern is divided into more than five azimuthal bins, such a difference is within the uncertainty of the measurements (, the dashed green line) and can therefore be neglected, thus simplifying the calculations.
Since the isotropic difference scattering signal arises from the changes in the radial distribution functions (RDFs) of the sample (Dohn et al., 2015), in Fig. 5 comprises contributions arising from the changes in the structure of the solute and of the solvation shells, as well as from the changes in temperature and density of the bulk solvent (Haldrup et al., 2010; Kjaer et al., 2013; Ihee et al., 2010; Cammarata et al., 2006). The anisotropic scattering signal in Fig. 5 arises from the changes in the structure of the solute and potentially anisotropic changes in the solvation shell structure. The anisotropic response of the bulk solvent is discussed at the end of the section and, in the case of water, its contributions are negligible at time scales longer than a few hundred femtoseconds after excitation. In this experiment, the signal arising from the changes in the structure of the solute () is strongly dominated by the contraction of the Pt—Pt bond due to the very electron-rich Pt atoms. As Fig. S2 shows, in the limited Q-range available for this experiment, possible contributions to arising from intramolecular changes other than the Pt—Pt bond contraction are found within the uncertainties of the measured signal. Based on these considerations, we construct a modelling framework where the photoinduced structural changes in the structure of the solute are parametrized through a single structural parameter, , which describes the changes in the Pt–Pt distance from the ground to the of the molecule. Specifically, is simulated from a set of structures derived by varying only the Pt–Pt distance of the density functional theory (DFT) optimized ground-state geometry of PtPOP. Details of the DFT calculations are provided in the supporting information. The full model used to fit the difference scattering signal measured at 4.5 ps after the of PtPOP is
where the first term on the right-hand side is , in the case of the isotropic signal, and , in the case of the anisotropic signal. and are calculated using equations (9) and (10), respectively; and represent the fraction of excited-state molecules contributing to the isotropic and anisotropic scattering, respectively. describes the changes in the structure of the solvation shells and β is a scaling factor. For the isotropic case, was calculated from the RDFs of the solute–solvent atom pairs (Dohn et al., 2015), as further described in the supporting information, while it was not necessary to include this term in the analysis of the anisotropic signal, given the good quality of the fit without the inclusion. Finally, describes the changes arising from the heating of the bulk water and its calculation is detailed in the supporting information. We note that, since the anisotropic signal is insensitive to the isotropic changes of the bulk solvent, fewer are used in the description of the anisotropic data compared with the isotropic contribution. Equation (15) was fit to the measured difference scattering signal within a standard minimization framework (Jun et al., 2010). From a simultaneous fit of the isotropic and the anisotropic scattering signal, we find = 0.24 ± 0.04 Å, = 2.6 ± 0.2%, = 1.5 ± 0.2%. The black line in Fig. 5 shows the model constructed from the best-fit results, which well describe the data ( = 1.9). is found in agreement with previous studies (Christensen et al., 2009); the 1σ confidence interval on this parameter (0.04 Å) is found to be smaller than that obtained from fitting the isotropic and the anisotropic signals separately (which gives 1σ confidence intervals of 0.05 Å and 0.09 Å, respectively). Since changes in population fraction are not expected on these single-ps time scales (Gray et al., 2017) and taking into account the decay of the anisotropic signal, as described below, we find that following the excitation event ∼2.6% of the molecules have been photoexcited and contribute to the isotropic signal, and approximately 1.8 % contribute to the anisotropic scattering. We interpret this difference as an indication of multi-photon excitation of PtPOP, where the multi-photon excitation takes place through transitions where the is not parallel to the Pt–Pt axis (Stiegman et al., 1987). This combined analysis may further allow the disentanglement of the contributions to the difference scattering signal arising from structural changes parallel or perpendicular with respect to the transition Specifically to PtPOP, the structural changes contributing to depend directly on their displacement with respect to the Pt–Pt axis ( in Fig. 3): equation (10) is maximum (and positive) when = 0° (i.e. when the structural change is displaced parallel to the Pt–Pt axis) and minimum (and negative) when = 90° (i.e. when the structural changes occur in the plane perpendicular to the Pt–Pt axis). However, as detailed above, the difference scattering signal is mainly dominated by the Pt–Pt contraction, and this hinders the investigation of additional structural parameters in the analysis described here.
The magenta line in Fig. 5 shows the signal obtained from an azimuthal integration of the anisotropic contribution re-projected onto the detector surface []. According to equation (11), this trace corresponds to the difference between the difference scattering signal obtained from an azimuthal integration of the 2D scattering pattern and the isotropic signal (〈 − 〉φ). For the PtPOP data, such a difference is found to be ∼20% of the magnitude of the isotropic contribution at time delays immediately after the excitation event, where the anisotropy is most pronounced. These observations show that the azimuthally integrated signal, at time scales shorter than the rotational correlation time of the solute, may not be a good approximation of the scattering arising from an isotropic ensemble of photoexcited molecules. This implies that at such time scales the use of azimuthally integrated scattering signals in subsequent analysis should be justified by comparing S0 and azimuthally integrated signals.
Following the procedure described above, the anisotropic difference scattering signals were extracted from measured difference scattering patterns up to 1 µs after Movie S1). Fig. 8 shows the decay of the anisotropy as a function of increasing pump–probe time delay. We find that 20% of the anisotropic contribution to the signal decays with a time constant of 3 ps ± 2 ps; while 80% decays with a time constant of 60 ps ± 10 ps. The latter is interpreted as the rotational correlation time of the PtPOP molecule in water. Using the Stokes–Einstein–Debye model (Horng et al., 1997), the rotational correlation time of a sphere can be described as
of PtPOP (seewhere η is the viscosity of water, kB is the T is the temperature and r is the radius of the sphere. Approximating the PtPOP molecule to a sphere, and taking r as half the longest inter-atomic distance in the PtPOP molecule (r = 4 Å), the rotational correlation time at room temperature can be estimated as 50 ps. This value is in agreement with the ∼60 ps found in this analysis.
Finally, we note that the neat solvent may exhibit anisotropic scattering independently of the photoexcited solute molecules: the laser pulse weakly perturbs the equilibrium structure of the solvent through non-resonant excitations predominantly aligned with the polarization of the laser, and the X-ray pulse probes the induced structural dynamics. Optically, the sample becomes birefringent, with different indices of refraction for light polarized parallel or perpendicular to the laser polarization axis. From an X-ray point of view, the scattering patterns are observed to be anisotropic and the analysis described above can be applied to extract the anisotropic scattering as a function of Q and time delay. The time evolution of the anisotropic X-ray scattering signal can be directly compared with the impulsive nuclear-coordinate response measured through techniques such as Raman-induced Kerr effect spectroscopy (Palese et al., 1994; Castner et al., 1995). Since the time scales of this nuclear response are known from optical Kerr effect studies, the anisotropic solvent scattering signal can be used to estimate the time-zero (the arrival time of the laser pulse at the sample) and the instrument response function (IRF) of the experiment (Biasin et al., 2016). The separate determination of these parameters significantly improves the analysis of the isotropic part of the scattering signal, by reducing the number of free parameters in the model used to fit the data (Biasin et al., 2016). This method of extracting the IRF and time-zero can be used as a diagnostic for any time-resolved XDS experiment on solvated molecules where the anisotropic solvent response can be clearly identified.
4. Conclusions
In summary, the formalism introduced by Baskin & Zewail (2006) and expanded by Lorenz et al. (2010) has been applied here to separate isotropic and anisotropic contributions from 2D difference scattering patterns measured upon of the solvated PtPOP molecule at an XFEL. The presented formalism has been directly implemented in quantitative structural analysis and we find that a combined analysis of the isotropic and anisotropic difference scattering signals helps the disentanglement of the many involved in the structural response of the sample to and enhances the structural sensitivity. This analysis approach is generally applicable for all molecular systems with a well defined transition and asymmetric structural response to provided a time-resolution shorter than the rotational dephasing time. Furthermore, we have discussed how the quantitative analysis of the anisotropic scattering signal can provide access to anisotropic solvent dynamics induced by the linearly polarized pump pulse and how these measurements can lead to independent determination of time-zero and the IRF of the experiment. In conclusion, we have demonstrated a method that allows quantitative interpretation of anisotropic scattering signals measured at XFELs from aligned ensembles of molecules. The information delivered by this method can benefit the overall interpretation of high-content scattering data from solution-state molecular systems.
5. Related literature
The following references, not cited in the main body of the paper, have been cited in the supporting information: Bahn & Jacobsen (2002); Becke (1988); Che et al. (1983); Dohn et al. (2017); Enkovaara et al. (2010); Himmetoglu et al. (2012); Larsen et al. (2009, 2017); Lee et al. (1988); Maurer & Reuter (2011); Mortensen et al. (2005); Ozawa et al. (2003); Pinto et al. (1980). Sørensen & Kjaer (2013); Yasuda et al. (2004); Ziegler et al. (1977).
Supporting information
X-ray difference scattering patterns measured at LCLS upon https://doi.org/10.1107/S1600577517016964/gb5063sup1.mov
of PtPOP in water solution as a function of time delay between laser pump and X-ray probe. DOI:Supporting information containing the image of the full scattering pattern from the sample, considerations on the structural https://doi.org/10.1107/S1600577517016964/gb5063sup2.pdf
involved in the photoinduced response of PtPOP, details on the calculation of the cage and solvent terms, and description of the DFT and QM/MM calculations that support the analysis presented in the main article. DOI:Acknowledgements
The authors would like to acknowledge the beamline scientists of the XPP Instrument at LCLS. Use of the Linac Coherent Light Source (LCLS), SLAC National Accelerator Laboratory, is supported by the US Department of Energy, Office of Science, Office of Basic Energy Sciences under Contract No. DE-AC02-76SF00515.
Funding information
Funding for this research was provided by: DANSCATT; Danish Council For Independent Research (grant No. DFF 4002-00272B to MMN, KBM, EB, MGL and AOD).
References
Bahn, S. R. & Jacobsen, K. W. (2002). Comput. Sci. Eng. 4, 56–66. Web of Science CrossRef CAS Google Scholar
Baskin, J. S. & Zewail, A. H. (2005). ChemPhysChem, 6, 2261–2276. Web of Science CrossRef CAS Google Scholar
Baskin, J. S. & Zewail, A. H. (2006). ChemPhysChem, 7, 1562–1574. Web of Science CrossRef CAS Google Scholar
Becke, A. D. (1988). Phys. Rev. A, 38, 3098–3100. CrossRef CAS PubMed Web of Science Google Scholar
Bennett, K., Kowalewski, M. & Mukamel, S. (2017). Phys. Rev. Lett. 119, 069301. CrossRef Google Scholar
Biasin, E., van Driel, T. B., Kjaer, K. S., Dohn, A. O., Christensen, M., Harlang, T., Chabera, P., Liu, Y., Uhlig, J., Pápai, M., Németh, Z., Hartsock, R., Liang, W., Zhang, J., Alonso-Mori, R., Chollet, M., Glownia, J. M., Nelson, S., Sokaras, D., Assefa, T. A., Britz, A., Galler, A., Gawelda, W., Bressler, C., Gaffney, K. J., Lemke, H. T., Møller, K. B., Nielsen, M. M., Sundström, V., Vankó, G., Wärnmark, K., Canton, S. E. & Haldrup, K. (2016). Phys. Rev. Lett. 117, 013002. Web of Science CrossRef Google Scholar
Borfecchia, E., Garino, C., Salassa, L. & Lamberti, C. (2013). Philos. Trans. R. Soc. London A, 371, 20120132. Web of Science CrossRef Google Scholar
Brinkmann, L. U. L. & Hub, J. S. (2015). J. Chem. Phys. 143, 104108. Web of Science CrossRef Google Scholar
Burger, C., Hsiao, B. S. & Chu, B. (2010). Polym. Rev. 50, 91–111. Web of Science CrossRef CAS Google Scholar
Cammarata, M., Lorenc, M., Kim, T., Lee, J., Kong, Q., Pontecorvo, E., Lo Russo, M., Schiró, G., Cupane, A., Wulff, M. & Ihee, H. (2006). J. Chem. Phys. 124, 124504. Web of Science CrossRef Google Scholar
Castner, E. W., Chang, Y. J., Chu, Y. C. & Walrafen, G. E. (1995). J. Chem. Phys. 102, 653–659. CrossRef CAS Web of Science Google Scholar
Che, C. M., Herbstein, F. H., Schaefer, W. P., Marsh, R. E. & Gray, H. B. (1983). J. Am. Chem. Soc. 105, 4604–4607. CrossRef CAS Web of Science Google Scholar
Chergui, M. & Collet, E. (2017). Chem. Rev. 117, 11025–11065. Web of Science CrossRef CAS Google Scholar
Chollet, M., Alonso-Mori, R., Cammarata, M., Damiani, D., Defever, J., Delor, J. T., Feng, Y., Glownia, J. M., Langton, J. B., Nelson, S., Ramsey, K., Robert, A., Sikorski, M., Song, S., Stefanescu, D., Srinivasan, V., Zhu, D., Lemke, H. T. & Fritz, D. M. (2015). J. Synchrotron Rad. 22, 503–507. Web of Science CrossRef CAS IUCr Journals Google Scholar
Christensen, M., Haldrup, K., Bechgaard, K., Feidenhans'l, R., Kong, Q., Cammarata, M., Russo, M. L., Wulff, M., Harrit, N. & Nielsen, M. M. (2009). J. Am. Chem. Soc. 131, 502–508. Web of Science CrossRef PubMed CAS Google Scholar
Dohn, A. O., Biasin, E., Haldrup, K., Nielsen, M. M., Henriksen, N. E. & Møller, K. B. (2015). J. Phys. B, 48, 244010. Web of Science CrossRef Google Scholar
Dohn, A. O., Jonsson, E. O., Levi, G., Mortensen, J. J., Lopez-Acevedo, O., Thygesen, K. S., Jacobsen, K. W., Ulstrup, J., Henriksen, N. E., Mller, K. B. & Jonsson, H. (2017). J. Chem. Theory Comput. Accepted. Google Scholar
Driel, T. B. van, Kjaer, K. S., Biasin, E., Haldrup, K., Lemke, H. T. & Nielsen, M. M. (2015). Faraday Discuss. 177, 443–465. PubMed Google Scholar
Driel, T. B. van, Kjaer, K. S., Hartsock, R. W., Dohn, A. O., Harlang, T., Chollet, M., Christensen, M., Gawelda, W., Henriksen, N. E., Haldrup, K., Kim, K. H., Ihee, H., Kim, J., Lemke, H. T., Sun, Z., Sundström, V., Zhang, W., Zhu, D., Møller, K. B., Nielsen, M. M. & Gaffney, K. J. (2016). Nat. Commun. 7, 13678. Google Scholar
Enkovaara, J., Rostgaard, C., Mortensen, J. J., Chen, J., Dułak, M., Ferrighi, L., Gavnholt, J., Glinsvad, C., Haikola, V., Hansen, H. A., Kristoffersen, H. H., Kuisma, M., Larsen, A. H., Lehtovaara, L., Ljungberg, M., Lopez-Acevedo, O., Moses, P. G., Ojanen, J., Olsen, T., Petzold, V., Romero, N. A., Stausholm-Møller, J., Strange, M., Tritsaris, G. A., Vanin, M., Walter, M., Hammer, B., Häkkinen, H., Madsen, G. K. H., Nieminen, R. M., Nørskov, J. K., Puska, M., Rantala, T. T., Schiøtz, J., Thygesen, K. S. & Jacobsen, K. W. (2010). J. Phys. Condens. Matter, 22, 253202. Web of Science CrossRef Google Scholar
Esswein, A. J. & Nocera, D. G. (2007). Chem. Rev. 107, 4022–4047. Web of Science CrossRef PubMed CAS Google Scholar
Glownia, J. M., Natan, A., Cryan, J. P., Hartsock, R., Kozina, M., Minitti, M. P., Nelson, S., Robinson, J., Sato, T., van Driel, T., Welch, G., Weninger, C., Zhu, D. & Bucksbaum, P. H. (2016). Phys. Rev. Lett. 117, 153003. Web of Science CrossRef Google Scholar
Glownia, J. M., Natan, A., Cryan, J. P., Hartsock, R., Kozina, M., Minitti, M. P., Nelson, S., Robinson, J., Sato, T., van Driel, T., Welch, G., Weninger, C., Zhu, D. & Bucksbaum, P. H. (2017). Phys. Rev. Lett. 119, 069302. Web of Science CrossRef Google Scholar
Gray, H. B., Zli, S. & Vlek, A. (2017). Coord. Chem. Rev. 345(Suppl. C), 297–317. Google Scholar
Haldrup, K., Christensen, M. & Meedom Nielsen, M. (2010). Acta Cryst. A66, 261–269. Web of Science CrossRef IUCr Journals Google Scholar
Haldrup, K., Vankó, G., Gawelda, W., Galler, A., Doumy, G., March, A. M., Kanter, E. P., Bordage, A., Dohn, A., van Driel, T. B., Kjaer, K. S., Lemke, H. T., Canton, S. E., Uhlig, J., Sundström, V., Young, L., Southworth, S. H., Nielsen, M. M. & Bressler, C. (2012). J. Phys. Chem. A, 116, 9878–9887. Web of Science CrossRef CAS PubMed Google Scholar
Hensley, C. J., Yang, J. & Centurion, M. (2012). Phys. Rev. Lett. 109, 133202. Web of Science CrossRef PubMed Google Scholar
Himmetoglu, B., Marchenko, A., Dabo, I. & Cococcioni, M. (2012). J. Chem. Phys. 137, 154309. Web of Science CrossRef Google Scholar
Horng, M., Gardecki, J. A. & Maroncelli, M. (1997). J. Phys. Chem. A, 101, 1030–1047. CrossRef CAS Web of Science Google Scholar
Ihee, H., Wulff, M., Kim, J. & Adachi, S. (2010). Intl Rev. Phys. Chem. 29, 453–520. Web of Science CrossRef CAS Google Scholar
Jonas, D. M., Bradforth, S. E., Passino, S. A. & Fleming, G. R. (1995). J. Phys. Chem. 99, 2594–2608. CrossRef CAS Web of Science Google Scholar
Jun, S., Lee, J. H., Kim, J., Kim, J., Kim, K. H., Kong, Q., Kim, T. K., Lo Russo, M., Wulff, M. & Ihee, H. (2010). Phys. Chem. Chem. Phys. 12, 11536–11547. Web of Science CrossRef CAS Google Scholar
Kim, J., Kim, K. H., Kim, J. G., Kim, T. W., Kim, Y. & Ihee, H. (2011). J. Phys. Chem. Lett. 2, 350–356. Web of Science CrossRef CAS Google Scholar
Kim, J. G., Kim, K. H., Oang, K. Y., Kim, T. W., Ki, H., Jo, J., Kim, J., Sato, T., Nozawa, S., Adachi, S. & Ihee, H. (2015b). J. Phys. B, 48, 244005. Web of Science CrossRef Google Scholar
Kim, K. H., Kim, J. G., Nozawa, S., Sato, T., Oang, K. Y., Kim, T. W., Ki, H., Jo, J., Park, S., Song, C., Sato, T., Ogawa, K., Togashi, T., Tono, K., Yabashi, M., Ishikawa, T., Kim, J., Ryoo, R., Kim, J., Ihee, H. & Adachi, S. (2015a). Nature (London), 518, 385–389. Web of Science CrossRef CAS Google Scholar
Kjaer, K. S., van Driel, T. B., Kehres, J., Haldrup, K., Khakhulin, D., Bechgaard, K., Cammarata, M., Wulff, M., Sørensen, T. J. & Nielsen, M. M. (2013). Phys. Chem. Chem. Phys. 15, 15003–15016. Web of Science CAS PubMed Google Scholar
Kong, Q., Lee, J., Plech, A., Wulff, M., Ihee, H. & Koch, M. (2008). Angew. Chem. Int. Ed. 47, 5550–5553. Web of Science CrossRef CAS Google Scholar
Küpper, J., Stern, S., Holmegaard, L., Filsinger, F., Rouzée, A., Rudenko, A., Johnsson, P., Martin, A. V., Adolph, M., Aquila, A., Bajt, S., Barty, A., Bostedt, C., Bozek, J., Caleman, C., Coffee, R., Coppola, N., Delmas, T., Epp, S., Erk, B., Foucar, L., Gorkhover, T., Gumprecht, L., Hartmann, A., Hartmann, R., Hauser, G., Holl, P., Hömke, A., Kimmel, N., Krasniqi, F., Kühnel, K., Maurer, J., Messerschmidt, M., Moshammer, R., Reich, C., Rudek, B., Santra, R., Schlichting, I., Schmidt, C., Schorb, S., Schulz, J., Soltau, H., Spence, J. C. H., Starodub, D., Strüder, L., Thøgersen, J., Vrakking, M. J. J., Weidenspointner, G., White, T. A., Wunderer, C., Meijer, G., Ullrich, J., Stapelfeldt, H., Rolles, D. & Chapman, H. N. (2014). Phys. Rev. Lett. 112, 083002. Google Scholar
Lakowicz, J. R. (2006). Principle of Fluorescence Spectroscopy. Berlin: Springer. Google Scholar
Larsen, A. H., Mortensen, J. J., Blomqvist, J., Castelli, I. E., Christensen, R., Duak, M., Friis, J., Groves, M. N., Hammer, B., Hargus, C., Hermes, E. D., Jennings, P. C., Jensen, P. B., Kermode, J., Kitchin, J. R., Kolsbjerg, E. L., Kubal, J., Kaasbjerg, K., Lysgaard, S., Maronsson, J. B., Maxson, T., Olsen, T., Pastewka, L., Peterson, A., Rostgaard, C., Schiøtz, J., Schütt, O., Strange, M., Thygesen, K. S., Vegge, T., Vilhelmsen, L., Walter, M., Zeng, Z. & Jacobsen, K. W. (2017). J. Phys. Condens. Matter, 29, 273002. Google Scholar
Larsen, A. H., Vanin, M., Mortensen, J. J., Thygesen, K. S. & Jacobsen, K. W. (2009). Phys. Rev. B, 80, 195112. Web of Science CrossRef Google Scholar
Lee, C., Yang, W. & Parr, R. G. (1988). Phys. Rev. B, 37, 785–789. CrossRef CAS Web of Science Google Scholar
Lemke, H. T., Kjaer, K. S., Hartsock, R., van Driel, T. B., Chollet, M., Glownia, J. M., Song, S., Zhu, D., Pace, E., Nielsen, M. M., Benfatto, M., Gaffney, K. J., Collet, E. & Cammarata, M. (2017). Nat. Commun. 8, 15342. Web of Science CrossRef Google Scholar
Lorenz, U., Møller, K. B. & Henriksen, N. E. (2010). New J. Phys. 12, 113022. Web of Science CrossRef Google Scholar
Maurer, R. J. & Reuter, K. (2011). J. Chem. Phys. 135, 224303. Web of Science CrossRef Google Scholar
Minitti, M. P., Robinson, J. S., Coffee, R. N., Edstrom, S., Gilevich, S., Glownia, J. M., Granados, E., Hering, P., Hoffmann, M. C., Miahnahri, A., Milathianaki, D., Polzin, W., Ratner, D., Tavella, F., Vetter, S., Welch, M., White, W. E. & Fry, A. R. (2015). J. Synchrotron Rad. 22, 526–531. Web of Science CrossRef CAS IUCr Journals Google Scholar
Møller, K. B. & Henriksen, N. E. (2012). Struct. Bond. 142, 185–212. Google Scholar
Mortensen, J., Hansen, L. & Jacobsen, K. W. (2005). Phys. Rev. B, 71, 035109. Web of Science CrossRef Google Scholar
Ozawa, Y., Terashima, M., Mitsumi, M., Toriumi, K., Yasuda, N., Uekusa, H. & Ohashi, Y. (2003). Chem. Lett. 32, 62–63. Web of Science CSD CrossRef CAS Google Scholar
Palese, S., Schilling, L., Miller, R. J. D., Staver, P. R. & Lotshaw, W. T. (1994). J. Phys. Chem. 98, 6308–6316. CrossRef CAS Web of Science Google Scholar
Penfold, T. J., Tavernelli, I., Abela, R., Chergui, M. & Rothlisberger, U. (2012). New J. Phys. 14, 113002. Web of Science CrossRef Google Scholar
Philipp, H. T., Hromalik, M., Tate, M., Koerner, L. & Gruner, S. M. (2011). Nucl. Instrum. Methods Phys. Res. A, 649, 67–69. Web of Science CrossRef CAS Google Scholar
Pinto, M. A. F. D. R., Sadler, P. J., Neidle, S., Sanderson, M. R., Subbiah, A. & Kuroda, R. (1980). J. Chem. Soc. Chem. Commun. pp. 13–15. CrossRef Google Scholar
Roundhill, D. M., Gray, H. B. & Che, C. M. (1989). Acc. Chem. Res. 22, 55–61. CrossRef CAS Web of Science Google Scholar
Sørensen, T. J. & Kjaer, K. S., (2013). Solvent response in TRWAXS, https://sites.google.com/site/trwaxs/. Google Scholar
Stiegman, A. E., Rice, S. F., Gray, H. B. & Miskowski, V. M. (1987). Inorg. Chem. 26, 1112–1116. CrossRef CAS Web of Science Google Scholar
Takeda, H., Cometto, C., Ishitani, O. & Robert, M. (2017). ACS Catal. 7, 70–88. Web of Science CrossRef CAS Google Scholar
Van Kleef, E. H. & Powis, I. (1999). Mol. Phys. 96, 757–774. Web of Science CrossRef CAS Google Scholar
Veen, R. van der, Milne, C., El Nahhas, A., Lima, F., Pham, V., Best, J., Weinstein, J., Borca, C., Abela, R., Bressler, C. & Chergui, M. (2009). Angew. Chem. Int. Ed. 48, 2711–2714. Google Scholar
White, C. (1982). Organometallic Chemistry, Vol. 10, edited by E. W. Abel & F. G. A. Stone, pp. 315–350. The Royal Society of Chemistry. Google Scholar
Yang, J., Guehr, M., Shen, X., Li, R., Vecchione, T., Coffee, R., Corbett, J., Fry, A., Hartmann, N., Hast, C., Hegazy, K., Jobe, K., Makasyuk, I., Robinson, J., Robinson, M. S., Vetter, S., Weathersby, S., Yoneda, C., Wang, X. & Centurion, M. (2016b). Phys. Rev. Lett. 117, 153002. Web of Science CrossRef Google Scholar
Yang, J., Guehr, M., Vecchione, T., Robinson, M. S., Li, R., Hartmann, N., Shen, X., Coffee, R., Corbett, J., Fry, A., Gaffney, K., Gorkhover, T., Hast, C., Jobe, K., Makasyuk, I., Reid, A., Robinson, J., Vetter, S., Wang, F., Weathersby, S., Yoneda, C., Centurion, M. & Wang, X. (2016a). Nat. Commun. 7, 11232. Web of Science CrossRef Google Scholar
Yasuda, N., Uekusa, H. & Ohashi, Y. (2004). Bull. Chem. Soc. Jpn, 77, 933–944. Web of Science CSD CrossRef CAS Google Scholar
Zewail, A. H. (2000). Pure Appl. Chem. 72, 2219–2231. Web of Science CrossRef CAS Google Scholar
Ziegler, T., Rauk, A. & Baerends, E. J. (1977). Theor. Chim. Acta, 43, 261–271. CrossRef CAS Web of Science Google Scholar
Zipp, A. P. (1988). Coord. Chem. Rev. 84(Suppl. C), 47–83. Google Scholar
© International Union of Crystallography. Prior permission is not required to reproduce short quotations, tables and figures from this article, provided the original authors and source are cited. For more information, click here.