We describe an experiment designed to probe the electronic damage induced in nanocrystals of Buckminsterfullerene (C60) by intense, femtosecond pulses of X-rays. The experiment found that, surprisingly, rather than being stochastic, the X-ray induced electron dynamics in C60 are highly correlated, extending over hundreds of unit cells within the crystals1.
The precise details of the interaction of intense X-ray pulses with matter are a topic of intense interest to researchers attempting to interpret the results of femtosecond X-ray free electron laser (XFEL) experiments. An increasing number of experimental observations have shown that although nuclear motion can be negligible, given a short enough incident pulse duration, electronic motion cannot be ignored. The current and widely accepted models assume that although electrons undergo dynamics driven by interaction with the pulse, their motion could largely be considered 'random'. This would then allow the supposedly incoherent contribution from the electronic motion to be treated as a continuous background signal and thus ignored. The original aim of our experiment was to precisely measure the change in intensity of individual Bragg peaks, due to X-ray induced electronic damage in a model system, crystalline C60. Contrary to this expectation, we observed that at the highest X-ray intensities, the electron dynamics in C60 were in fact highly correlated, and over sufficiently long distances that the positions of the Bragg reflections are significantly altered. This paper describes in detail the methods and protocols used for these experiments, which were conducted both at the Linac Coherent Light Source (LCLS) and the Australian Synchrotron (AS) as well as the crystallographic approaches used to analyse the data.
One of the major purposes of X-ray free electron lasers (XFELs) is to develop a high throughput, high resolution approach to molecular imaging and dynamics. Structural biology depends on atomic scale information, traditionally limited to lower resolution X-ray crystallography techniques performed at third generation synchrotrons. Long exposure times which cause significant radiation damage in crystals, greatly affect the resolution achieved using traditional techniques. The snapshot diffraction imaging scheme2,3,4 employed at XFELs, involves collecting diffraction images from short pulse X-rays hitting either fixed target samples (which are translated across the beam focus) or samples injected into the path of the beam.
The XFEL pulse-sample interaction ultimately destroys the samples, due to the onset of severe radiation damage. The diffraction images are collected before the onset of this destruction due to the sub-100 fs pulse durations. The ability to determine high resolution structures from nanocrystals is rapidly becoming well established. However, dynamical processes which occur on femtosecond timescales under experimental imaging conditions offer deeper insights into atomic physics and can have a macroscopic effect on nanocrystals and their diffraction patterns5,6,7.
Whilst catastrophic structural damage is avoided on the femtosecond timescale during which a snapshot diffraction image is recorded, the power density of an XFEL pulse may be high enough to modify the electronic properties of the sample with which the X-rays interact7,8,9. An exploration of the physics of the interaction of intense coherent X-ray pulses with matter is not only of intrinsic scientific interest, but will be critically important to the interpretation of any experiment in which light from an XFEL pulse is used to explore structure.
In X-ray imaging experiments performed on single molecules, small clusters, or nanocrystals comprised of a few unit cells, perturbative analysis indicates that one should observe both a reduction in the apparent coherence of the scattered signal8, and the growth of a structureless background signal as a result of electrodynamical processes9. This experiment sought to assess the degree to which decoherence due to electrodynamical processes, occurs in powdered nanocrystalline C60 due to interaction with short XFEL pulses.
In this article, we provide details regarding the experimental procedure in which a highly ordered transient electronic structure from C60 nanocrystals is observed due to interaction with an XFEL pulse1. The diffraction pattern produced under these conditions is significantly different from that observed when the same sample is illuminated by lower power, but otherwise identical XFEL pulses, or when a synchrotron beam at the same photon energy is used. This difference is marked by the presence of Bragg peaks that are not seen in the two diffraction profiles corresponding to the low-power and synchrotron diffraction images. We demonstrate our analysis and model-fitting approach, used to confirm the presence of a dynamic electronic distortion induced by XFEL pulse-nanocrystal interaction.
1. C60 Powder Sample Preparation
2. Preliminary Australian Synchrotron Studies
3. Request XFEL Parameters for Beamline Setup
4. Record the Darkfield
5. Record the 10% Incident XFEL Flux Run
6. Record the 100 % XFEL Flux Run
7. XFEL Data Post-processing and Peak Analysis
XFEL Powder Diffraction
The data presented for the 100% incident flux XFEL powder diffraction is the result of summing more than 1000 single-shot measurements to produce a complete powder-ring with a resolution of better than 2 Å.
Powder diffraction profiles comparison
The Bragg peaks for the diffraction rings were identified and scaled to the first (most intense) peak reflection (111). Figure 3 shows the three different diffraction line profiles. By comparing the line profiles of the three diffraction patterns, we observe that the diffraction data recorded at the Australian Synchrotron is almost identical to the Bragg profile seen in the 10% XFEL data. Some very minor differences in the relative heights of the Bragg peaks, but not their positions are observed. In stark contrast, the profile of the 100% power XFEL powder diffraction data reveals the presence of additional peaks not seen in the 10% XFEL data profile, nor in the Synchrotron data profile. The locations of these extra reflections are identified in Table 1. In order to interpret these differences, an adjustment to the model of expected diffraction from a room temperature FCC C60 crystal was constructed.
X-ray diffraction modelling of the room temperature FCC C60 structure
The intensity of powder diffraction peaks associated with Bragg reflections from a crystal is given by
(1),
where is the scattering vector, K is the scale factor, is the multiplicity factor, Lp is the Lorentz-polarization factor, W() is the peak profile function and M is the number of C60 molecules contained in the scattering volume located at positions rm. The molecular form factor (MFF), , for a C60 molecule is given by
(2),
where rj is the position of the jth carbon atom in the molecule and fc is the atomic scattering factor of the carbon atom.
The unit cell parameters of the crystal define the positions of allowed reflections for an X-ray powder diffraction pattern. Using the known room temperature FCC parameters (unit cell length, molecule positions within the unit cell) of C60, together with the experimental geometry in the X-ray diffraction experiment, the expected positions of peaks (Bragg reflections) can be calculated using the MFF for C60 and Eq. 1 and Eq. 2.
X-ray diffraction modelling of 100% XFEL data
We begin by assuming that significant distortions/transformations or displacements of the nuclei from their ideal positions do not occur during the 32 fs duration of the incident pulse as suggested in prior studies23,24. Rather, that significant change in the intensities seen in the 100% XFEL data must instead be driven by motions of the electronic structure of the C60 molecules. In the following we describe a model that reproduces the experimentally observed features of the 100% XFEL diffraction data, via a modification of the centro-symmetric distribution of the C60 molecules.
In its normal, neutral state, the crystalline structure of C60 is maintained by dipolar forces that are induced by instantaneous fluctuations in its electron density. Under the experimental conditions described here, however, the ionization of the system generates a strong internal electric field that induces electric dipole moments in the molecules by polarization. Previously the formation of dipoles in C60 has only been observed in single molecules and small clusters using optical techniques such as UV spectroscopy25. Here however, the redistribution of the electron density observed is evidently both long-range and long-lived relative to the duration of the XFEL pulse so that its effects are observed in the crystallographic X-ray diffraction pattern.
This results in the alignment of neighbouring dipoles via a Coulomb interaction, and a decoupling of the electronic structure from the underlying nuclear structure on timescales on the order of 10 fs. This charged alignment affects the resulting symmetry of the C60 molecule (see Figure 4). The loss of the spherical symmetry of the molecule leads to an additional phase contribution to the scattering amplitude, since the MFFs of C60 molecules are no longer real but complex functions.
A periodically varying MFF was used to model the occurrence of an asymmetric molecular charge distribution in which the distribution of the electron density of the mth molecule is displaced relative to its position in the crystal structure. With this modification to the C60 MFF, we were able to replicate the intensity profile seen in the 100% XFEL data.
Eq. 2 provides the basis for constructing an expression for the scattering factor, which captures the long-range electronic correlations formed from the XFEL-induced dipoles in the 100% XFEL data. From this a new MFF function, modified to account for the polarized C60 molecules, can be constructed:
(3),
where is the MFF of the ideal C60 molecules (given by Eq. 2) and defines the polarization vector of the XFEL induced dipole. In the limit , Eq. 3 approximates Eq. 2, and the room temperature 10% power diffraction data is recovered. As increases, the symmetry of the molecule is altered, and the ratios of all the possible diffraction peaks begin to vary. The actual distribution of polarized molecules in a cubic lattice affects the resulting diffraction pattern.
When , the symmetry of the C60 molecule is altered and the ratios of all the possible diffraction peaks begin to vary relative to the low-power diffraction pattern. To fit the data to this model, values of were explored, showing good agreement in the 20° ≤ 2θ ≤ 30° range of scattering angles for .
The intended purpose of this experiment was to measure the degree to which stochastic photoionisation of the K-shell in carbon atoms affects the diffracted intensities measured for FCC C60 nanocrystals. Photoionisation of the K-shell electrons in carbon atoms (electron binding energy = 284 eV) modifies the atomic scattering factors, fc, seen as a reduced scattering amplitude within the high scattering regions. K-shell holes in carbon atoms within C60 molecules arranged in a crystalline lattice causes modifications of the scattering amplitudes of the Bragg reflections.
We expected to observe a growing isotropic background, dependent on the photon flux applied to powdered nanocrystal samples according to the following fundamental assumptions: 1) that the photoionisation of the K-shell in carbon is the dominant process in the sample-XFEL interaction, 2) that photoionisation of individual carbon atoms is not correlated to any other atoms in the crystal, 3) that photoionized electrons remain delocalized for the duration of the pulse and hence contribute to the continuous background signal.
What we actually observed in the experiment was the presence of strong, forbidden reflections in room temperature, FCC nanocrystals of C60 when the sample was subjected to the 100% power XFEL pulses. Delocalized, random ionization events cannot account for the observed forbidden reflections.
Figure 3 shows the appearance of these forbidden reflections, coinciding with a substantial reduction in the intensities of the allowed FCC reflections. These changes cannot be described by any specific orientational ordering of ideal C60 molecules in the crystal lattice.
According to our analysis1, a correlated, non-centrosymmetric charge distribution on each C60 molecule (Eq. 4), has proved the only means of generating a model powder diffraction profile which matches the experimental data (seen in Figure 5). For comparison, all data and models are shown together, but offset vertically with respect to one another, on the same axis in Figure 6.
Figure 1. XFEL Powder Diffraction Sample Setup and Geometry
(a) The sample holder used for the fixed target scanning mode of C60 crystal powder. The sample frame is constructed from aluminium. Measurements indicated are in units of mm. Approximate dimensions of sample cells are 2 mm x 12 mm. (b) Photograph of C60 crystal powder applied in three of the cells (seen as darkly coloured cells) with polyimide backing applied as a support (the yellow film on top of the sample holder). (c) Schematic of the C60 experiment. The sample is raster scanned in x-y directions in the snapshot imaging scheme. K-B mirrors focus the XFEL beam to a spot size of 300 nm x 300 nm at the sample. Samples are held in vacuum to stabilize the sample conditions and minimize the possibility of X-ray interaction with scattering sources other than the sample. Incoming XFEL pulses hit the crystal powder held in the sample holder cells, and a diffraction pattern is recorded at the CSPAD detector. A resolution of 1.5 Å is achieved by setting the sample to detector distance to 79 mm. Please click here to view a larger version of this figure.
Figure 2. The CSPAD
Note that the white scale bar in a), b) and d) represents 40 mm.
(a) CSPAD darkfield. The detector is composed of 32 rectangular modules, the positions of which can be altered by moving concentrically outward to record high-angle diffraction as needed. (b) Summed raw data frames (top right-hand quadrant, over 1000 frames summed) prior to background and darkfield correction. (c) Individual diffraction snapshots demonstrating sparsity of the diffraction signal. (d) Diffraction profile showing well defined powder diffraction rings performed by summing 1500 diffraction frames with background signal subtraction applied to individual frames.
Figure 3. Powder Diffraction Data
(a) Azimuthally averaged diffraction patterns for the 10% XFEL dataset, 100% XFEL dataset and the Synchrotron dataset. Positions of FCC Bragg peaks are indicated consistent with a room temperature C60 FCC structure. (b) Inset region showing reflections present in the 100% FCC structure between scattering angles 10⁰ ≤ 2θ ≤ 13⁰ not seen in the other two profiles. (c) Inset region showing the different peak profile in the 100% XFEL data between the scattering angles 20⁰ ≤ 2θ ≤ 28⁰. The 10% XFEL data and the synchrotron data both satisfy the selection rules for FCC structures composed of electronically centrosymmetric molecules. However the presence of extra peaks (reflections) seen in the 100% XFEL data violate these selection rules. Please click here to view a larger version of this figure.
Figure 4. Transient Distortion of C60
Visualisation of the alignment of the dipoles within the FCC lattice structure during the correlated electronic transient stage. C60 molecules are represented by blue spheres and the red tips represent the direction of the ordered dipoles.
Figure 5. Powder Diffraction Model
Powder diffraction profile generated by modelling the FCC structure for C60 (using Eq 1 and 2) compared to a model of the C60 FCC structure subjected to a 100% intensity XFEL pulse (using Eq 1 and 3). Identified Bragg peaks are labelled. A region of interest (20° ≤ 2θ ≤ 30°) is highlighted by the dotted line. Although the FCC model describes the intensity of the allowed reflections well, it does not explain the presence of a number of additional peaks (see Figure 2a and b) observed for the 100% intensity XFEL data. The reason for this is that the simple translation of the molecular cluster (Figure 3) along the crystallographic axis of the cubic lattice gives us an incomplete picture of the orientational ordering of polarized C60 molecules in the cubic lattice. By contrast the 100% XFEL model, which takes into account ionisation-induced alignment of the dipoles within the FCC lattice (as shown in Figure 4), reproduces all of the additional peaks observed in the 100% intensity XFEL data. Please click here to view a larger version of this figure.
Figure 6. Powder Profile Comparison Between Model and Data
A qualitative comparison of the line profiles for the three diffraction patterns recorded under different illumination conditions experimentally. In addition, the line profiles calculated using Equations 2 and 3 using our model are shown. It is clear that the introduction of a periodically modified MFF, the 100% XFEL model line profile agrees with our 100% XFEL data.
Measured scattering angles of extra reflections (deg.) | Calculated scattering angles of extra reflections (deg.) |
21.31 | 21.25, 21.45 |
23.23 | 22.99, 23.02, 23.39 |
24.44 | 24.29, 24.43, 24.47, 24.64 |
26.6 | 26.51, 26.67 |
Table 1. Bragg Reflections Seen in XFEL Data
The set of Bragg reflections measured within the 20⁰ ≤ 2θ ≤ 30⁰ for the 100% XFEL diffraction data as well as those calculated using Eqns. 1 – 4.
Position of the molecule | Alignment |
(0,0,0) | |
(0.5,0.5,0) | |
(0.5,0,0.5) | |
(0,0.5,0.5) |
Table 2. FCC Molecular Alignment During Transient Correlated Phase
This table describes the alignment of polarized C60 molecules during the transient correlated phase of the crystal experienced during the XFEL pulse.
Calibration of diffraction data frames.
The .XTC files (which contain the data from a complete run) contain calibration parameters that define the geometrical arrangement of CSPAD modules (shown in Figure 2a) during the experiment. The correct arrangement of data recorded on individual modules is crucial to assemble the individual diffraction data images comprising data recorded in each run. At the time the experiment was performed the location of the calibration file containing the correct parameters was not automatically set up and manual computation was required by the team to correct the issue. Due to the extra time spent performing calibration of the data there was a time-lag between setting a snapshot run dataset and checking the success of the run via a darkfield and background subtracted summation of image frames in the data set.
Crystal sizes.
In some of the initial XFEL snapshot runs, strong single crystal Bragg reflections were seen in some of the image frames. This resulted from some of the C60 sample not being crushed finely enough. Observing optical reflections from crushed powder indicates that the crystal facets are too large (correspond to the wavelength of visible light ~ 400-700 nm). The powder should be checked for these reflections at the crushing stage, and if strong, single crystal Bragg reflections are seen in the data the powder needs to be further crushed.
Since the results of this experiment were not expected or planned for successful powder diffraction data collection for the C60 sample was only obtained at two extreme intensity settings (10% and 100% flux). Beam time at the facility is limited and hence any set-up, computation, or sample processing errors and issues have a large impact on an experimental plan. The two most widely separated incident intensity points were prioritized and there was insufficient beam time available to collected reliable statistics for any intermediate points. Therefore, we were not able to experimentally assess the trigger point in terms of XFEL flux at which this transient phase change occurs.
Preliminary studies.
Collecting powder diffraction data at the Australian Synchrotron, from the same C60 sample as measured at the XFEL. Synchrotrons are routinely used to screen for suitable XFEL targets26, and in the present case positively confirmed that at 10% XFEL intensity, the diffraction data was consistent with the ground state FCC structure of C60.
Sample and detector attenuation.
Calibration of the incident flux through adjustment of the silicon attenuators upstream of the sample was essential, especially since the effect being studied was intensity dependent. Construction of a suitable aluminium attenuator at the detector, matched to the incident flux was also critical.
Hitting the sample at the location of the beam focal point.
The location of the KB focal spot at the XFEL was also essential to observe the reported phenomenon, since the flux density on the sample must be sufficient to induce the formation of dipoles throughout the crystal. Measuring the size of craters created by the XFEL beam in a YAG crystal using optical microscopy, as well as performing a fine sample scan along the optical axis and looking at the diffraction intensity was used to determine the location of the focal plane.
In future implementations of this work a larger number of incident intensities as well as pulse durations will be explored. This work has potential implications for upcoming experiments analyzing the diffraction data collected from nanocrystals at XFEL sources. It also provides new insights into the fundamental interaction of XFELs with matter, highlighting that XFELs have the potential to explore new physics not accommodated within conventional crystallography.
The authors have nothing to disclose.
The authors acknowledge the support of the Australian Research Council Centre of Excellence in Advanced Molecular Imaging. Portions of this research were carried out at the LCLS, a national user facility operated by the Stanford University on behalf of the U.S. Department of Energy, Office of Basic Energy Sciences. We acknowledge the travel funding provided by the International Synchrotron Access Program managed by the AS and by the Australian government. In addition, some of this research was undertaken on the MX1 and MX2 beamlines at the AS, Victoria, Australia. Author contributions: B.A. was responsible for planning and managing all experimental aspects of the project. Experiments were designed by B.A., R.A.D., V.S., C.D., and G.J.W. B.A., H.M.Q., K.A.N., and R.A.D. wrote the original LCLS proposal. D.W., R.A.D., R.A.R., A.V.M., E.C., and S.W. carried out the simulation work. B.A., R.A.D., C.D., V.S., M.W.M.J., R.A.R., N.G., F.H., G.J.W., S.B., M.M., M.M.S., A.G.P., C.T.P., A.V.M., and K.A.N. collected the experimental data at the LCLS. S.W., V.A.S. and R.A.D collected experimental data at the Australian Synchrotron. C.T.P. and A.V.M. led the experimental data conversion and analysis. B.A., C.D., N.G., and E.B. were responsible for the sample holder design and testing. R.A.R, B.A., S.W., A.V.M and H.M.Q wrote this manuscript. The formulation of electronic damage within coherence theory is performed by H.M.Q. and K.A.N.; R.A.D. conceived the idea to apply this formalism to C60.
Macroscopic 99.5+ % pure C60 | SES RESEARCH | ||
Pestle and mortar | Sigma Aldrich | used for crushing C60 powder; | |
Aluminium sheet | used for constructing sample holder | ||
kapton polyimide film | Du Pont | http://www.dupont.com/products-and-services/membranes-films/polyimide-films/brands/kapton-polyimide-film/ | |
CXI beamline | SLAC | http://scripts.iucr.org/cgi-bin/paper?yi5003 | |
safety glasses | |||
biosafety cabinet |