The development of Laue techniques for single-pulse diffraction of chemical complexes: time-resolved Laue diffraction on a binuclear rhodium metal-organic complex

doi:10.1107/S0108767311011883

Acta Crystallogr A. 2011 July 1; 67(Pt 4): 319–326.

Published online 2011 May 10. doi: 10.1107/S0108767311011883

PMCID: PMC3121236

The development of Laue techniques for single-pulse diffraction of chemical complexes: time-resolved Laue diffraction on a binuclear rhodium metal-organic complex

Anna Makal,^a Elzbieta Trzop,^a Jesse Sokolow,^a Jaroslaw Kalinowski,^a Jason Benedict,^a and Philip Coppens^a^*

^aChemistry Department, University at Buffalo, State University of New York, Buffalo, NY 14260-3000, USA

Correspondence e-mail: coppens/at/buffalo.edu

Received March 7, 2011; Accepted March 30, 2011.

This article has been cited by other articles in PMC.

Abstract

A modified Laue method is shown to produce excited-state structures at atomic resolution of a quality competitive with those from monochromatic experiments. The much faster data collection allows the use of only one or a few X-ray pulses per data frame, which minimizes crystal damage caused by laser exposure of the samples and optimizes the attainable time resolution. The method has been applied to crystals of the α-modification of Rh₂(μ-PNP)₂(PNP)₂ (BPh₄)₂ [PNP = CH₃N(P(OCH₃)₂)₂, Ph = phenyl]. The experimental results show a shortening of the Rh—Rh distance in the organometallic complex of 0.136 (8) Å on excitation and are quantitatively supported by quantum-mechanical (QM)/molecular-mechanics (MM) theoretical calculations which take into account the confining effect of the crystal environment, but not by theoretical results on the isolated complex, demonstrating the defining effect of the crystal matrix.

Keywords: Laue techniques, single-pulse diffraction, quantum-mechanical/molecular-mechanics calculations, QM/MM calculations, time-resolved X-ray crystallography

1. Introduction

Time-resolved (TR) X-ray crystallography at advanced light sources is a frontier area bound to undergo dramatic new developments with the advent of ultrabright X-ray free-electron lasers (XFELs). Earlier TR studies of chemical systems in which atomic resolution was reached have used monochromatic radiation for which methods are well developed and a reasonable accuracy can be achieved (Coppens, Benedict et al., 2010

; Cailleau et al., 2010

). The use of monochromatic data avoids the wavelength dependence of both the scattering process and the detector response and does not suffer from the broadening of the reflection maxima inherent in the Laue technique. However, with a narrow bandwidth of monochromatic radiation only a very small fraction of the photons in the synchrotron beam are productively used, even at undulator-equipped beamlines. Thus a longer exposure time and therefore more pulses are required, which limits the time resolution that can be achieved. Furthermore, the longer exposure time implies more extensive laser exposure, thereby enhancing the temperature increase of the sample due to light absorption.

It follows that for picosecond-timescale TR diffraction at synchrotron sources the use of polychromatic radiation is imperative. To achieve this goal we have developed a number of methods to improve the accuracy and interpretation of Laue measurements. They include the RATIO technique in which the measured I _ON/I _OFF ratios are used in combination with a set of monochromatic data collected at the same temperature (Coppens et al., 2009

), photo-Wilson plots to estimate the temperature increase due to the laser exposure of the sample (Schmøkel et al., 2010

), the definition of R factors specific for dynamic structure crystallography (Coppens, Kamiński & Schmøkel, 2010

), and a scaling technique for relative scaling within multi-crystal data sets collected on the same substance. The scaling is required for calculation of the photodifference maps using all available reflections. The technique is based on the moduli of the absolute fractional intensity change and is described in §2.3. We describe a single-pulse Laue diffraction experiment of a binuclear Rh^I complex in which the techniques summarized above have been applied and a comparison of the results with theoretical calculations. A preliminary account of the results of this study has been published (Benedict et al., 2011

2. The experiment

2.1. Crystal structure and photochemical properties

The crystal structures of the three polymorphs of Rh₂(μ-PNP)₂(PNP)₂(BPh₄)₂, (1), where PNP = CH₃N(P(OCH₃)₂)₂ and Ph = phenyl (Fig. 1

), and related salts were reported by Mague and coworkers (Blakley et al., 1989

; Mague, 1995

). On excitation of the α form with 337 nm light at 10 K a strong emission band was observed at 730 nm, but at ~ 50 K a second band develops at 780 nm, while the band at 730 nm weakens and disappears at 80 K. The change is attributed to a thermal transition to a second triplet state of somewhat lower energy, separated from the higher state by an activation barrier (Blakley et al., 1989

). The large Stokes shift observed indicates a significant structural rearrangement on excitation. Owing to a broadening of the spot shape in the Laue patterns on cooling attributed to strain induced in the specimen crystals, the current experiment was performed at 225 K by cooling in a nitrogen stream. Unit-cell information from the monochromatic experiment and luminescence properties are summarized in Table 1

. Full crystallographic information on the ground-state structure is given in the supplementary material.¹ The 225 K temperature of the experiment reported here implies that the excited state with the 780 nm emission was studied. The resulting shortening of the emission lifetime from the value of 46.5 (3) µs reported at 10 K (Blakley et al., 1989

) does not affect the current experiment as the delay between the 35 ps pump and 100 ps probe pulses was only 100 ps.

Figure 1

Molecular diagram of (1) and labeling of atoms.

Table 1

Cell dimensions and emission properties

2.2. Data collection

Time-resolved intensity data were collected at the 14-ID BioCARS station at the Advanced Photon Source. Two different undulator settings were used with peak energies at 12 and 15 keV, respectively. The corresponding spectral distributions are illustrated in Fig. 2

. Use of the polychromatic beam (ΔE/E = ~ 8%) allows collection of a frame of data with a single 100 ps synchrotron pulse (Fig. 3

). The samples were pumped with 35 ps pulses of a Ti:sapphire laser tuned to a wavelength of 390 nm. The laser power varied between 0.45 and 0.6 mJ mm⁻². It is noteworthy that laser powers beyond 0.6 mJ mm⁻² led to a decrease in the ON/OFF intensity changes, an effect attributed to secondary excitation of the photo-induced excited state possibly creating on a non-radiative relaxation pathway. To maximize the number of weak reflections observed, the pump–probe cycle for each frame was repeated three times before detector read-out. The pump–probe delay time was set at 100 ps. Much longer delay times caused spot extension in the ON frames, as observed in earlier experiments (Coppens, Benedict et al., 2010

). Laser-OFF and Laser-ON frames were collected in immediate succession to minimize the effect of long-range fluctuations in the beam parameters. The OFF/ON pump–probe cycle was repeated ten times for each frame to allow subsequent statistical filtering of the data. This redundancy is made possible by the extremely fast data collection, limited essentially by the ~ 4 s read-out time used with the MAR detector. Six sufficiently large data sets covering at least half of the reciprocal space within the resolution limit, totaling three each at 12 and 15 keV, were collected each at successive values of the ϕ angle of the diffractometer with Δϕ = 5° and were used in the subsequent analysis. Other runs had to be aborted because of premature laser damage of the samples as evident from pronounced weakening of the diffraction pattern and smearing of the diffraction spots. Details of the six data sets are summarized in Table 2

and Table S1 of the supplementary material.

Figure 2

Experimental spectral distribution at BIOCARS beamline 14-ID prior to the experiments.

Figure 3

Single-pulse Laue image collected at the APS ID-14, BIOCARS beamline.

Table 2

Description of the data sets

The mutual consistence of the data sets was examined by the construction of correlation plots (Coppens et al., 2005

), two of which are shown in Fig. 4

Figure 4

Examples of plots showing the correlation between ON/OFF ratios of individual hkls in different data sets. Red lines correspond to ±1σ(R). The green line represents a trend line for an ideal correlation ( An external file that holds a picture, illustration, etc.
Object name is a-67-00319-efi3.jpg

An external file that holds a picture, illustration, etc.
Object name is a-67-00319-efi3.jpg

). The data in the upper figure (more ...)

2.3. Data processing

The reflections were indexed using the program PRECOGNITION (Ren, 2010

), as well as using a new program LAUEUTIL currently being developed in our laboratory (to be published). Spot integration was performed by the program LAUEGUI (Messerschmidt & Tschentscher, 2008

; Peters, 2003

). The program is based on the seed–skewness method (Bolotovsky & Coppens, 1997

), which does not require profile fitting. Only reflections with I > 10σ(I), where σ(I) is the statistical standard deviation of the measured intensity I, were considered in the calculation of the ON/OFF ratios R and related response ratios η, respectively defined as

A mathematical equation, expression, or formula.
Object name is a-67-00319-efd1.jpg

and

A mathematical equation, expression, or formula.
Object name is a-67-00319-efd2.jpg

Several filtering steps were applied. For each hkl on a frame the ten repeated ratio measurements were rejected whenever the standard deviation for the sample was larger than 0.5. Reflections with less than ten repeated measurements were also rejected. The resulting ratios were averaged over repeated and symmetry-equivalent measurements with the program SORTAV (Blessing, 1997

). Outliers were down-weighted by the use of the robust/resistant Tukey biweighting criterion based on the deviation of individual measurements from the median (Blessing, 1997

). Details are given in Table S2. The ten individual intensities of two reflections with I _ON > I _OFF are shown in Fig. S1.

2.4. Preliminary scaling of the data sets collected on different crystals

To obtain a photodifference map for visualization of the light-induced changes and for extraction of initial atomic positions for the excited-state least-squares refinement, it is necessary to combine the data sets, as reliable Fourier maps essentially require a complete set of terms in the Fourier summation. In the RATIO method the ON data are obtained by multiplication of the monochromatic intensities by the synchrotron-determined ratios R,

A mathematical equation, expression, or formula.
Object name is a-67-00319-efd3.jpg

For multi-sample data sets with different fractional excitation the ratios must be reduced to a common scale before the Fourier summation is performed. To accomplish this the η values of individual reflections are scaled in SORTAV (Blessing, 1997

) by

|η|

_set/

|η|

_all, in which left angle bracket

|η|

_all and

|η|

_set are the averages over all measured reflections and over the reflections in the specific set respectively. Corrected values of R are obtained from

A mathematical equation, expression, or formula.
Object name is a-67-00319-efd4.jpg

The scaled and merged ratios are then used in equation (2) to generate uniform I _ON values for the calculation of photodifference maps.

2.5. Analysis of the sample temperature during the ON measurements

An estimate of the temperature increase caused by the laser pulse can be obtained from a photo-Wilson plot (Schmøkel et al., 2010

; Vorontsov & Coppens, 2005

). The plot gives an estimate of ΔB, which with the value of left angle bracket

from the monochromatic data measured at the same temperature can be converted to a temperature scale factor representing the relative temperature increase with the equation

A mathematical equation, expression, or formula.
Object name is a-67-00319-efd5.jpg

As may be expected, the values of k_B for each set are correlated with left angle bracket

|η|

_set as shown in Fig. 5

. However, they become almost equal after the scaling according to equations (3c) and (2) is implemented.

Figure 5

Relation between the temperature scale factor k _B and 1/k(η), the inverse of the response ratio scale factor for the 12 and the 15 keV data sets.

The value of k _B for the set of combined data equals 1.07.

2.6. Least-squares refinement

The data were refined with the program LASER, which is based on the refinement of the response ratios (Vorontsov et al., 2010

). The program assumes a random spatial distribution of the excited-state species with a scattering formalism

A mathematical equation, expression, or formula.
Object name is a-67-00319-efd6.jpg

and

A mathematical equation, expression, or formula.
Object name is a-67-00319-efd7.jpg

in which F _ON and F _OFF are the structure factors for the light-ON and light-OFF measurements, GS and ES are the ground and excited state, respectively, and P is the percentage conversion to the light-induced state. An alternative algorithm assuming clustering of the light-induced species leads to quite different calculated values of the ON structure factors and has not been successful in any of our previous studies (Vorontsov & Coppens, 2005

). Differences between An external file that holds a picture, illustration, etc.
Object name is a-67-00319-efi1.jpg

and

An external file that holds a picture, illustration, etc.
Object name is a-67-00319-efi2.jpg

, which would imply a motion or deformation of the GS molecules in the light-exposed crystal are expected to be minimal in this work because of the low conversion percentages achieved (Table S3). The LASER program allows for the refinement of up to six different data sets with individual values of the conversion percentage P and the temperature scale factor k_B and calculates R factors specifically designed for photocrystallographic experiments (Vorontsov et al., 2010

; Coppens, Kamiński & Schmøkel, 2010

Refinement was performed on the six data sets listed in Table 2

. Only data with |η|/σ(η) > 1 were used in the refinement giving a total of 2256 independent reflections for the refinement of 42 excited-state structure parameters. They included the six positional parameters of the two Rh atoms, temperature scale factors k_B, and the excited-state occupancies for each of six experimental runs plus three rotational and translational parameters of the four rigid PNP groups. The ground-state Rh positions were taken as starting points in the excited-state refinement of the Rh atoms. The ligands pivoted around anchor atoms, the positions of which were allowed to translate. For the two terminal ligands the initial positions of the anchor atoms were the locations of the excited-state Rh atoms to which the ligands were attached, obtained in a preliminary refinement of the Rh positions only. For the bridging ligands the midpoint of the Rh atoms was used as the starting position in the subsequent refinement. Structural results are presented in Tables 3

and 4

and S2–S4. Final agreement factors, defined elsewhere (Coppens, Kamiński & Schmøkel, 2010

), are R(η) = 0.342 and R(R) = 0.034. Further numerical details on the refinement are given in Tables S5–S7.

Table 3

Shifts of the Rh and P atoms on excitation and standard deviations (Å)

Table 4

Experimental excited-state (ES) and ground-state (GS) bond lengths and their differences (Å)

3. Discussion of results

3.1. Photodifference maps

The photodifference map based on all independent reflections, calculated at isosurfaces of ± 0.25 e Å⁻³, is shown in Fig. 6

(a). The map indicates a sideward displacement of Rh1 and a displacement of Rh2 towards Rh1, thus suggesting both a rotation and a shortening of the Rh—Rh vector. In addition it shows a depletion of the electron density in the regions of P atoms. The compensating positive electron density at these atoms becomes visible when contours are drawn at the lower 0.175 e Å⁻³. One of these maps is shown in Fig. 6

(b).

Figure 6

Photodifference maps with isosurfaces (red positive, blue negative) of ±0.25 e Å⁻³ (a) and ±0.175 e Å⁻³ (b) (calculated with XDGRAPH of the XD program set; Volkov et al., (more ...)

3.2. The least-squares results

The results of the refinement on the excited-state Rh atom positions and the ligand translations and rotations are listed in Table 3

. In agreement with the photodifference maps a significant change in the positions of the Rh atoms is observed. As Rh2 moves towards Rh1, but the latter is displaced in a perpendicular direction, the resulting bond shortening of 0.136 (8) Å (Table 4

) is essentially a result of the Rh2 shift. The displacements of the P atoms are generally significant, but hardly affect the Rh—P distances, indicating that the ligands follow the metal atoms. The rigid-body rotations do not exceed 2.5°, but are more pronounced for the terminal ligand attached to Rh1 than for that linked to Rh2 as may be expected from the shift perpendicular to the Rh—Rh bond of the former atom.

A diagram of the molecule illustrating the atomic displacements on excitation is shown in Fig. 7

Figure 7

Molecular diagram indicating atomic shifts on laser excitation. Blue lines: ground state; green lines: excited state. The excited-state atoms are marked with ‘e’ after the label. Ellipsoids are drawn at the 5% probability level for the (more ...)

3.3. Comparison with theory

Theoretical DFT calculations were performed both with GAUSSIAN09 (Frisch et al., 2009

) and with the Amsterdam Density Functional program (Velde et al., 2001

). Two different functionals BP86 and B3LYP were used with a number of different basis functions. It is noticeable that even the isolated molecule Rh—Rh distance, determined as 3.183 Å, is calculated considerably longer by all theoretical isolated molecule treatments as shown in Tables 5

and 6

. The results are slightly better with the BP86 functional applied with both programs (Table 6

), but the bond length and its shortening of 0.136 (8) Å on transition to the triplet state are not well reproduced with either functional, independent of the basis set used. The large discrepancies are at least in part attributed to the effect of the crystal environment on the geometry of both the ground and excited states.

Table 5

Isolated-molecule theoretical bond lengths obtained with GAUSSIAN09 B3LYP functional (Å)

Table 6

Isolated-molecule theoretical bond lengths obtained with the BP86 functional

In an earlier study we applied the mixed quantum-mechanical (QM)/molecular-mechanics (MM) technique to explore the effect of the crystalline environment on the molecular dimensions. In our application the crystal environment is kept constant at the experimental values, while the complex in the crystalline cavity is optimized subject to the constraint imposed by the neighboring atoms (Kaminski et al., 2010

). Results of the QM/MM optimization for the cation of (1) performed with the ADF program using the BP86 functional with different basis sets are tabulated in Table 7

Table 7

Theoretical results (BP86 functional, ADF program) for (1) confined in a crystal cavity as calculated with the QM/MM method (Å)

The results confirm that the molecule in the crystal is severely constrained by its crystalline environment, even in the ground state. The ground-state Rh—Rh bond length is best approximated by the smaller basis set DZ (double zeta) calculation, but this result is likely to be fortuitous as the HOMO and LUMO diagrams obtained with this basis set significantly deviate from those from the more advanced calculations. The contraction calculated with the larger TZP basis is within three standard deviations from the experimental result. The HOMO and LUMO orbital diagrams shown in Fig. 8

indicate an antibonding-to-bonding transition. They are remarkably similar to the diagrams observed in the time-resolved study of the [Rh^II ₂(1,8-diisocyano-p-menthane)₄]²⁺ ion (2). However, in the latter case a much larger contraction of 0.85 (5) Å was observed in crystals of the PF₆ ⁻ salt. We note that the ground-state Rh—Rh bond length of (2) at 4.496 (1) Å is much longer than that of (1), leaving room for a much larger contraction.

Figure 8

The HOMO (a) and the LUMO (b) of the title complex (TZP basis, BP86 functional, QM/MM results, isolated-molecule orbitals are essentially identical). Isovalue surfaces at ±0.0275 a.u. Blue/blue–grey positive, red/orange negative. (more ...)

The experimental changes in the Rh—P distance as determined by the rigid-body refinement of the PNP ligands are listed in Table 8

. The changes are of the same magnitude as the standard deviations derived from those of the rotation and translation parameters of the refinement. Nevertheless, except for the Rh1—P8 distance there seems to be a clear correlation with the changes determined theoretically with the QM/MM method.

Table 8

Calculated (QM/MM) and observed changes in the Rh—P distances Δ on excitation (Å)

4. Conclusions

We conclude that the improved Laue method is capable of producing excited-state structures at atomic resolution of a quality at least comparable with those of monochromatic experiments. The much faster data collection allows the use of one of a few X-ray pulses per frame of data, and thus minimizes crystal damage caused by laser exposure of the samples. Equally important, the shorter exposure time optimizes the attainable time resolution. The experimental results are quantitatively supported by the QM/MM calculations which take into account the confining effect of the crystal environment, but not by theoretical results on the isolated complex. This effect has to be taken into account in the interpretation of photochemical properties of molecules in the solid state.

Supplementary Material

Crystal structure: contains datablock(s) I. DOI: 10.1107/S0108767311011883/sh5129sup1.cif

Click here to view.^{(34K, cif)}

Structure factors: contains datablock(s) I. DOI: 10.1107/S0108767311011883/sh5129Isup2.hkl

Click here to view.^{(848K, hkl)}

Supplementary material file. DOI: 10.1107/S0108767311011883/sh5129sup3.pdf

Click here to view.^{(135K, pdf)}

Acknowledgments

We would like to thank Robert Henning and Tim Graber of BioCARS for assistance and advice with the data collection. This work was funded by the Division of Chemical Sciences, Geosciences, and Biosciences, Office of Basic Energy Sciences of the US Department of Energy, through grant DEFG02-ER15372. Use of the BioCARS Sector 14 was supported by the National Institutes of Health, National Center for Research Resources, under grant number RR007707. The time-resolved facility at Sector 14 was funded in part through a collaboration with Philip Anfinrud (NIH/NIDDK). The Advanced Photon Source is supported by the US Department of Energy, Office of Basic Energy Sciences, under Contract No. W-31-109-ENG-38.

Footnotes

¹Crystallographic information on the ground-state structure discussed in this paper is available from the IUCr electronic archives (Reference: SH5129). Services for accessing these data are described at the back of the journal.

References

Benedict, J. B., Makal, A., Sokolow, J. D., Trzop, E., Scheins, S., Henning, R., Graber, T. & Coppens, P. (2011). Chem. Commun. 47, 1704–1706. [PMC free article] [PubMed]
Blakley, R. L., Yin, Y., Lloyd, C., Mague, J. T. & McPherson, G. L. (1989). Chem. Phys. Lett. 157, 398–402.
Blessing, R. H. (1997). J. Appl. Cryst. 30, 421–426.
Bolotovsky, R. & Coppens, P. (1997). J. Appl. Cryst. 30, 244–253.
Cailleau, H., Lorenc, M., Guérin, L., Servol, M., Collet, E. & Buron-Le Cointe, M. (2010). Acta Cryst. A66, 189–197. [PubMed]
Coppens, P., Benedict, J., Messerschmidt, M., Novozhilova, I., Graber, T., Chen, Y.-S., Vorontsov, I., Scheins, S. & Zheng, S.-L. (2010). Acta Cryst. A66, 179–188. [PMC free article] [PubMed]
Coppens, P., Kamiński, R. & Schmøkel, M. S. (2010). Acta Cryst. A66, 626–628. [PubMed]
Coppens, P., Pitak, M., Gembicky, M., Messerschmidt, M., Scheins, S., Benedict, J., Adachi, S., Sato, T., Nozawa, S., Ichiyanagi, K., Chollet, M. & Koshihara, S. (2009). J. Synchrotron Rad. 16, 226–230. [PMC free article] [PubMed]
Coppens, P., Vorontsov, I. I., Graber, T., Gembicky, M. & Kovalevsky, A. Y. (2005). Acta Cryst. A61, 162–172. [PubMed]
Frisch, M. J. et al. (2009). GAUSSIAN09, Revision A.1. Gaussian Inc., Pittsburgh, Pennsylvania, USA.
Kaminski, R., Schmøkel, M. S. & Coppens, P. (2010). J. Phys. Chem. Lett. 1, 2349–2353.
Mague, J. T. (1995). Inorg. Chim. Acta, 229, 17–25.
Messerschmidt, M. & Tschentscher, T. (2008). Acta Cryst. A64, C611.
Peters, J. (2003). J. Appl. Cryst. 36, 1475–1479.
Ren, Z. (2010). Precognition User Guide RenZ Research, Inc, Illinois, USA.
Schmøkel, M. S., Kamiński, R., Benedict, J. B. & Coppens, P. (2010). Acta Cryst. A66, 632–636. [PubMed]
Velde, G. T., Bickelhaupt, F. M., Baerends, E. J., Guerra, C. F., Gisbergen, S. J. A. V., Snijders, J. G. & Ziegler, T. (2001). J. Comput. Chem. 22, 931–967.
Volkov, A., Macchi, P., Farrugia, L. J., Gatti, C., Mallinson, P. R., Richter, T. & Koritsanszky, T. (2006). XD2006. Middle Tennessee State University, USA, Università di Milano and CNR-ISTM Milano, Italy, University of Glasgow, Scotland, State University of New York at Buffalo, USA, and Freie Universität Berlin, Germany.
Vorontsov, I. I. & Coppens, P. (2005). J. Synchrotron Rad. 12, 488–493. [PubMed]
Vorontsov, I., Pillet, S., Kamiński, R., Schmøkel, M. S. & Coppens, P. (2010). J. Appl. Cryst. 43, 1129–1130.

Articles from Acta Crystallographica Section A: Foundations of Crystallography are provided here courtesy of
International Union of Crystallography