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-results-19 Crystal lattice structure diagram, atomic arrangement, [111] direction, 0.7 Å spacing.](/files/ftp_upload/56296/56296fig4.jpg)
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) | ![figure-results-22 Crystallographic direction symbol δr‖[111], concept of lattice orientation in materials science.](/files/ftp_upload/56296/56296table2_1.jpg) |
| (0.5,0.5,0) | ![figure-results-23 Crystallographic direction symbol δr‖[111], concept of lattice orientation in materials science.](/files/ftp_upload/56296/56296table2_1.jpg) |
| (0.5,0,0.5) | ![figure-results-24 Vector notation in crystal physics with δr parallel to [111], featuring a directional vector symbol.](/files/ftp_upload/56296/56296table2_2.jpg) |
| (0,0.5,0.5) | ![figure-results-25 Vector notation in crystal physics with δr parallel to [111], featuring a directional vector symbol.](/files/ftp_upload/56296/56296table2_2.jpg) |
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.