Applications and Limits of Time-to-Energy Mapping of Protein Crystal Diffraction Using Energy-Chirped Polychromatic XFEL Pulses

: A broadband energy-chirped hard X-ray pulse has been demonstrated at the SwissFEL (free electron laser) with up to 4% bandwidth. We consider the characteristic parameters for analyzing the time dependence of stationary protein diffraction with energy-chirped pulses. Depending on crystal mosaic spread, convergence, and recordable resolution, individual reﬂections are expected to spend at least ≈ 50 attoseconds and up to ≈ 8 femtoseconds in reﬂecting condition. Using parameters for a chirped XFEL pulse obtained from simulations of 4% bandwidth conditions, ray-tracing simulations have been carried out to demonstrate the temporal streaking across individual reﬂections and resolution ranges for protein crystal diffraction. Simulations performed at a higher chirp (10%) emphasize the importance of chirp magnitude that would allow increased observation statistics for the temporal separation of individual reﬂections for merging and structure determination. Finally, we consider the fundamental limitation for obtaining time-dependent observations using chirped pulse diffraction. We consider the maximum theoretical time resolution achievable to be on the order of 50–200 as from the instantaneous bandwidth of the chirped SASE pulse. We then assess the ability to propagate ultrafast optical pulses for pump-probe cross-correlation under characteristic conditions of material dispersion; in this regard, the limiting factors for time resolution scale with crystal thickness. Crystals that are below a few microns in size will be necessary for subfemtosecond time resolution. three dimensional, unaveraged free electron laser simulation code.


Introduction
With serial femtosecond crystallography (SFX) at X-ray Free Electron Lasers (XFELs), pump-probe data of protein structural dynamics with femtosecond time resolution can be obtained. In an SFX experiment, protein structures are determined by merging thousands of stationary diffraction patterns from single femtosecond XFEL pulse exposures. Several time-resolved experiments have been reported with SFX [1][2][3][4][5][6][7] and have allowed, for example, measurement of difference electron density maps identifying conformational changes as early as approximately 140 fs for the photoactive yellow protein (PYP) [1], 49-406 fs in bacteriorhodopsin [7], and 1 ps in the rsEGFP2 fluorescent protein [4]. Pump-probe SFX can therefore provide a description of the mechanistic connection between structural dynamics at different timescales (from ultrafast to millisecond) in photoactivated systems.
The bandwidth for an XFEL pulse, when operated in self-amplified spontaneous emission (SASE) mode, is on the order of ∆E/E ≈ 0.2% [8][9][10]. The SASE process relies on the emission of spontaneous radiation by electrons in the first section of an undulator that is then itself the X-ray seed for the FEL [11][12][13][14]. As spontaneous radiation emission from the electrons is stochastic, the SASE radiation arises from microbunching with spectral and energy fluctuations [13][14][15][16]. With a typical SASE bandwidth in the range of 0.2%, the majority of the reflections in protein SFX are partials [17], meaning that the full intensities of Bragg reflections are not recorded from a single crystal. A Monte Carlo integration technique is commonly used in order to determine the average values of structure factor amplitudes [18,19]. Based on the postrefinement procedure for oscillation data [20,21], methods for partiality assignment have been developed [22][23][24][25][26], in addition to integration, that can significantly improve merging statistics of SFX data [22,27]. These define partialities by modeling the overlap between the reciprocal space lattice point volume and the set of Ewald spheres that constitutes the beam [24]. Initial estimates of the full intensities are obtained from partiality values and these estimates are then improved in an iterative manner through agreement across the entire dataset [22]. An alternative approach would be to use polychromatic X-ray beams that fully integrate many reflections, as in the Laue technique. In Laue crystallography, the bandwidth is large enough that many reflections sit entirely within the Ewald sphere volume in reciprocal space and are therefore fully recorded. A serial polychromatic crystallography approach has been demonstrated with "pink"-beam undulator radiation and a bandwidth of 2.5% at synchrotrons [28,29]. At XFELs, the use of increased bandwidth for serial crystallography could improve the partiality determination and postrefinement results, though structure factor amplitude noise from the spectral fluctuation intrinsic to SASE operation remains present. A source of increased uncertainty is also the reduced flux under conditions of reflection, as the pulse energy is essentially divided over a larger frequency range.
The new compact XFEL source at the Paul Scherrer Institute, SwissFEL, is able to operate in a broad-bandpass mode: simulated ultrafast hard X-ray SASE pulses have a bandwidth as large as ∆E/E ≈ 4% [30,31]. Such pulses at SwissFEL display an energy chirp, meaning that the central frequency varies with time. The X-ray energy chirp is obtained through an energy chirp of the electron beam: in contrast to standard operation, the energy chirp of the electron bunch is purposefully increased through overcompression in the last bunch compressor [30]. This scheme can maximize the chirp of the electron beam, and thus that of the X-rays. Currently, 4% is the maximum achievable chirp expected for hard X-rays at SwissFEL, though the possibility exists to increase the chirp with a dechirper device [32,33].
The possibility of exploiting an energy chirp to obtain time resolution within the probing period of the applied pulse in macromolecular crystallography was first proposed by Moffat in 2002 Faraday Discussions [34], shortly before the construction of the first XFEL. The energy chirp leads to streaking of time into reciprocal space in the diffraction pattern-because different wavelengths diffract to different locations on the detector, diffraction from the earliest energies is physically separated from diffraction caused by the energies arriving later in time. Through Laue indexing, an individual reflection can be assigned to the wavelength that stimulated it, and thus to its specific time within the pulse if there is an energy chirp. With sufficient bandwidth, a deconvolution of multiplicity would additionally be needed, such that the time-to-energy relation could still be extracted. Time-dependent information of the pump-probe experiment is therefore intrinsically present in individual images collected with chirped broadband pulses. The time interval across which structure amplitude differences are recorded is set by the chirp magnitude of the X-ray probe pulse. Moffat additionally proposed that if the pump pulse occurs within the probe, each image will also include unpumped reference structure factor amplitudes which will aid the scaling of time-resolved differences [34].
Here, we extend this method to XFEL pulses in order to obtain ultrafast time resolution with the availability of femtosecond energy-chirped hard X-ray pulses and consider the time resolution that could be obtainable. Several numerical simulation codes for SASE FELs have been developed to explore the effects of different beamline components on performance and to aid in the interpretation of experimental results [35][36][37][38]: they provide three-dimensional representations of the pulse by solving the FEL equations [39], which describe both the propagation of the electrons and the evolution of the radiation field. Calculation methods are constantly being improved. A recent version of GENESIS, one of the most widely used codes, is able to simulate the propagation and SASE of the entire electron beam through the undulator [40]. Pulse profiles computed in this way are presented here in order to simulate and assess the time dependence of experimental diffraction.
The aim of this work is thus to evaluate the time resolution, ∆t res , that would theoretically be achievable in a polychromatic XFEL experiment from considerations of physical parameters for protein crystallography. The X-ray pulse length, spectral bandwidth over time, and chirp magnitude determine the resolution and are obtained through simulations of the XFEL provided by SwissFEL. A ray-tracing simulation demonstrates the chirped diffraction patterns generated by the polychromatic XFEL pulse.
It has been shown that a 4% "extra pink" beam mode at SwissFEL will probe a large enough volume of reciprocal space to fully integrate multiple reflections [41]. As the polychromatic beam is chirped, it should be possible to identify which spots correspond to different energies and therefore to different points in time. Here, we consider a negative energy chirp, such that diffraction spots that appear at the leading edge of the pulse (t in ) are stimulated by the highest energies. When evaluating full reflection integration, it is noted that, while the total spectral bandwidth that generates the polychromatic image is equal to the magnitude of the chirp (e.g. 4%), we analyze the "instantaneous" frequency (IF) and "instantaneous" bandwidth (BW in ) associated with every instance of the pulse that is obtained from Wigner-Ville distribution analysis of the simulated field [42,43]. The determination of IF in a pulse becomes useful when its spectral characteristics vary over time. The Wigner-Ville distribution (WVD) is a two-dimensional distribution of the signal in the time-frequency domain. Projection of the distribution on one axis results in either the temporal or the spectral structure of the pulse. The IF corresponds to the first moment of the WVD with respect to frequency [44,45], while BW in is defined as the standard deviation, or spread of frequencies, present at a given time [46].
The shortest resolvable time interval (∆t res ) should be dictated by the "instantaneous duration" of BW in . This is given by the time-bandwidth product (∆ν∆t = 0.44 for a Gaussian pulse): the time duration obtained from such time-bandwidth product corresponds to the minimum value at which a pulse with specific FWHM can theoretically be fully compressed. The bandwidth of the SASE pulse is increased relative to that of the contributions of individual spikes in a way that is analogous to how the experimental bandwidth of optical pulses can be expanded by adding together multiple optical sources each with a different carrier frequency [47]. The spectral width of a single isolated SASE spike at 9 keV has been recorded to be ≈ 10 eV [48] at LCLS, which, through the time-bandwidth product, corresponds to a duration on the order of 200 as. The chirped pulse at SwissFEL consists of the superposition of "single spikes", varying in time and central frequency. The instantaneous bandwidth is determined by the pulse length and chirp of an individual spike. Therefore, the maximum time resolution that is possible for analyzing diffraction using chirped XFEL pulses is limited by the duration and chirp of the individual spike, which is analyzed from simulations here. Experimentally, based on the SwissFEL simulations presented, intervals shorter than 200 as will not be readily distinguishable with a SASE pulse, as the autocorrelation shows the intensity distribution within the 200-as duration. A coherence time shorter than 100 as has been discussed for harder X-ray (24.8 KeV) and low bunch charge conditions at the European XFEL [49].
The time duration that each reflection is stimulated for will depend on the chirp and on the bandwidth required to fully integrate it (BW req ). The spot bandwidth is given as in [50]: where η + γ is the effective mosaic spread (the sum of mosaicity and beam convergence) and θ is the Bragg angle. There is a dependence on resolution because the Ewald spheres diverge from the reciprocal space origin ( Figure 1). Figure 1. Effect of beam convergence (γ) and mosaicity (η) on the reflection bandwidth (BW req ). BW req depends on the sum of γ and η and on the distance of the reciprocal lattice point from the reciprocal space origin (d*). The minimum (λ min ) and maximum (λ max ) wavelengths of the spot bandwidth are shown together with the energy gradient across the reflection.

SwissFEL Chirped Broadband Simulations
Simulations of the SwissFEL pulse at three chirps (1%, 2%, and 4%) are presented in this paper ( Figure 2). For the simulations, the results of measured beam parameters from the machine were used, when it was configured for large bandwidth operation [31]. This characterization of the electron beam was done at the injection point into the undulator, using a transverse deflector in order to obtain time-resolved beam parameters. We measured the current profile, the slice emittances, and the longitudinal phase space distribution. To the latter, we fitted a polynomial to determine the mean energy of the electron beam slices. With these parameters, we generated the input decks for the GENESIS simulation. Therefore, the results of the 2% chirp case actually reflect the existing machine settings, and the modeled pulse energy and spectral width agree well with the observed quantities. For the other two cases, the polynomial fit to the longitudinal phase space distribution was scaled by the factors 2.0 and 0.5 respectively, while keeping the other input parameters the same. This was done without demonstrating the feasibility of a different chirp at the machine. While it is straightforward to reduce the chirp, a larger chirp requires additional hardware, such as a dechirper, or a significant change in the compression configuration. Both approaches will be explored in upcoming machine studies at SwissFEL.
The simulation output contains the radiation power and phase as a function of the time within the time-window of the simulation. This is the basis for the calculation of the WVD, which then is used to calculate the time-resolved central frequency and the associated width of the chirped FEL pulse. From the WVD, the flux, IF, and instantaneous FWHM (FWHM in ) can be extracted for the simulated chirps ( Figure 2). We evaluate two conditions. Firstly, we consider the time intervals given by BW in at different times in the pulse. Secondly, we assess the ratio between BW in and BW req to determine the expected integration time across resolutions. Among the available simulations, the largest chirp and bandwidth of 4% is expected to result in the best time resolution and was thus chosen for the following analysis. The values for IF and FWHM in vary significantly over the course of the pulse ( Figure 2), which means that the two points above need to be evaluated separately for different time bins. For calculations, only the central 20-fs of the WVD (between 52 and 72 fs) were included, as this is where the flux is concentrated (Figure 3). The early and late times contribute only very weakly to the total pulse intensity. This slicing reduces the available chirp from 4% to approximately 1.5%, which is the realistically usable chirp. A short-range "single-shot" autocorrelation of this part of the SwissFEL simulation ( Figure 3d) shows that the width of a single spike at SwissFEL is expected to be on the same range as at LCLS (200 as).

Determination of ∆t res
The central 20-fs of the WVD (between 52 and 72 fs) were split into 20 bins of 1-fs duration, each with its mean IF and FWHM in . For each bin, BW req was calculated at resolutions between 0.9 and 10 Å and for effective mosaic spread values between 0.05-0.5 • using Equation 1 (Figure 4a). As mentioned above, the effective mosaic spread is a sum of the beam convergence and mosaicity contributions: at XFEL facilities, beam convergence at around 9 keV ranges from a few microradians (on the order of 10 −4 degrees) [51][52][53] for the unfocused beam, to tens of microradians if KB mirrors (e.g., ≈ 47 µrad ≈ 0.002 • [54]) or compound refractive lenses (e.g., ≈ 91 µrad ≈ 0.005 • [55]) are used for focusing. As this is considerably smaller than the mosaic spread of protein crystals, the latter is the determining factor. The range of mosaic spread used here considers very low mosaicity values, such as those from crystals of the photoactive yellow protein (PYP, η ≈ 0.04 • [56]), as well as values up to ten times larger. From the BW req of a reflection, it is also possible to calculate the amount of time it is stimulated for (Figure 4b). The ∆t IF /∆t for the chirp is 5.6 eV/fs, the required ∆t IF is given by BW req , and ∆t is solved for. At mosaicity = 0.5 • , low-resolution reflections are expected to spend up to 8 fs in reflecting condition. On the other hand, high-resolution reflections at low mosaicity values will be fully "instantaneously" integrated by BW in (Figure 4b,c). Figure 5 shows the predicted integration for a lowand a high-resolution reflection during the same interval in the pulse along with the linear chirp in IF. The 1.5 Å reflection is fully integrated in 500 as, while the 5 Å one is calculated to remain in reflecting condition for 2 fs. These time durations are dictated by BW req for the individual reflections and by the mean BW in in the time bin. Time durations were calculated from BW in values through time-bandwidth products. The instantaneous durations given by the mean FWHM in of each time bin specify the formal limit to the time resolution and vary between 50 and 150 as (Figure 6a). These values were calculated using a time-bandwidth product of 0.44 and therefore underestimate the pulse duration as the chirp of the single spike is not included in their calculation: in the presence of spectral phase, the spike duration is correspondingly longer [48]. Figure 4c shows the ratios between BW req and the available BW in given by the WVD. As expected, ratios decrease as a function of resolution, with values all below 0.4 at resolutions lower than 4 Å. From these ratios, the number of energy points per reflection (that is, how many BW in are needed to integrate it) was derived (Figure 6b). Realistically, many of the energy contributions will overlap on the diffraction pattern, so the number of points per reflection was only calculated for bins with BW in /BW req ≥ 0.2, setting a limit of 5 energy/time points per reflection. Bins with BW in /BW req < 0.2 were automatically assigned a value of 10 energy/time points per reflection to indicate that there should be more than 5 time points per reflection but these are not expected to be resolved. Figure 6b shows values that range from 1 to 10 in the 0-5 Å range. All bins at resolutions lower than 5 Å were assigned a value of 10 and are not shown in the plot. A ray-tracing simulation was carried out to evaluate the images that would be generated by the pulse and the observable gradients across reflections.

Ray-Tracing Simulation and Time Mapping
The program SIM-MX by Kay Diederichs [57] was used to simulate still diffraction images from a single PYP crystal (2PHY.pdb) with experimentally determined parameters such as unit cell dimensions and mosaic spread [58]. The full list of input values for the program can be found in Table 1 and the definition of each is described in the original publication [57]. As stills are not implemented in SIM-MX, these were simulated by using a very small oscillation range (0.01 • ). One "still" image was generated for each of the 20 time bins by tracing 10,000 rays per frame: the simulation wavelength and its standard deviation were given by the mean IF and FWHM in of each bin. These stills were then summed to obtain a diffraction pattern for the polychromatic, chirped pulse. The reflections within this image can be tagged according to their energy, and thus according to their time of arrival at the crystal. In Figures 7 and 8, intensities on the diffraction pattern are color-coded based on the energy that stimulated them, which allows one to relate them to their time within the pulse. Diffraction images were first simulated using the 1.5% chirped central slice of the SwissFEL pulse. A second simulation was then run with a hypothetical 10% chirp to emphasize the achievable effect. The diffraction patterns show new reflections appearing during the duration of the pulse as well as low-resolution reflections that contain multiple energy contributions (Videos S1-S3). The observable intensity gradients are significantly enhanced by the higher chirp ( Figure 8).

Discussion
The analysis presented here concludes that, given the polychromatic chirped XFEL pulse currently available at SwissFEL, the instantaneous bandwidths throughout the pulse correspond to time intervals between ≈ 50-150 as. In principle, this matches the shortest possible time that high-resolution reflections can spend in stimulated conditions ( Figure 5). In order to analyze diffraction to the maximum possible time resolution determined by the instantaneous bandwidth, the precise position in the chirped pulse needs to be determined by fully resolving the X-ray energy for each reflection, with additional deconvolution of rare multiples. This would rely on stationary analysis with full knowledge of the chirped pulse spectral parameters. The intensity profile would need to be known as wavelength normalization from polychromatic data requires multiple symmetry-related reflections to be measured which are stimulated by different energies to derive a λ-curve. For spontaneous X-ray sources, they are typically fitted to Chebyshev polynomials [59][60][61][62][63]. Time-resolving structure factor amplitudes within the chirped pulse therefore prevents λ-curve fitting. The chirp magnitude may additionally contribute to the experimental time resolution, as between 2 fs and 8 fs elapse for the instantaneous frequency to be modified by the magnitude of the instantaneous bandwidth ( Figure 3). The practical limit to the time resolution, purely from the point of X-ray diffraction observation, will therefore be affected by the ability to determine the X-ray energy for each reflection, which reports on the progression of time within the chirped XFEL pulse. Practically, a stationary analysis would find a wavelength assignment with high accuracy. While the chirp magnitude from the simulation remains relatively constant at approximately ≈ 5 eV/fs in the central part of the pulse that contains most intensity (Figure 3c), the instantaneous bandwidth varies between values of 40 eV and 10 eV (Figure 3a). Therefore, a resolution of 200 as would require experimental determination of a change of ≈ 1 eV in the instantaneous frequency from the integrated intensities. Considering the spectral structure of XFEL pulses, there may be limits arising from signal-to-noise to fully exploit the time-resolved analysis of chirped-pulse crystallography.
From the simulations performed, gradients across individual low-resolution reflections are observable as well as different high-resolution reflections integrated by different energies over the course of the single pulse. It is clear that an increase in chirp magnitude allows a better distinction of specific energies: maximizing this parameter will become crucial when addressing the signal-to-noise ratio in an experimental setting. While this result is based solely on the properties of the XFEL pulse and reciprocal space, there are still several limitations that currently lower the experimentally achievable time resolution from its theoretical limit. The main limiting factors to ∆t res will be the duration of the pump pulse, the group velocity mismatch, and the pump-probe timing at the beamline. Optical pulses used for protein excitation at XFELs up to now have been on the order of 100 fs. In addition, the absolute pump-probe time at the interaction point in the experiment, determined through a specific timing tool in SFX experiments, carries an uncertainty of 30-50 fs [64][65][66]. Development of hollow core fibers for ultrafast laser pulse generation [67][68][69] has enabled production of single-cycle pulses, with a duration of only a few femtoseconds, or attosecond half-cycle pulses [47,[70][71][72] that have also been used for ultrafast pump-probe spectroscopy experiments [73,74]. The implementation of such pulses in pumping sample at SFX experiments would improve the earliest observable time point after excitation. Furthermore, Moffat's original discussion proposed that the pump pulse occur within the X-ray probe, so that each image would contain an internal dark control as well as multiple time points. Laser pulses with attosecond or few-femtosecond durations will be a requirement to have optical excitation occur within the 100-fs XFEL pulse.
In the context of using shorter pump pulses and with the idea of pursuing femtosecond time resolution as well as Moffat's original proposal [34] of pumping within the probe pulse, chromatic dispersion of the laser light has to be considered: there is a physical limit dictated by the material dispersion of the optical pump pulse when propagating through the protein crystal and the delivery matrix or solution. Chromatic dispersion is larger in pulses with large bandwidths, and so will become significant if few-cycles ultrashort lasers are used for pumping; the fact that the protein sample is excited by the pump is also important, as group velocity dispersion (GVD) behavior differs at absorption [75]. Though crystal birefringence results in an orientation dependence of the GVD value [76], its contribution is not further considered here. Figure 9 shows final pulse durations given a range of extremely short femtosecond pulses and different group delay dispersion (GDD) values. The output pulse is calculated with the following standard relationship [77]: where GDD is the GVD multiplied by sample thickness and t in is the initial pulse duration. As there are (to our knowledge) no published GVD values for visible light absorbing protein crystals, we used the concentration dependence of the GVD for Rhodamine 6G dye solution near its 600-nm absorption resonance [78] to estimate the GDD for a protein microcrystal: for a 5-µm thick crystal with 60-mM chromophore concentration, the GVD is scaled to 2258 fs 2 /mm and the GDD is calculated to be 11.29 fs 2 . This yields an output pulse of 10.85 fs, which is more than three times the 3-fs input pulse. The method of sample delivery will also add to chromatic dispersion: the stretching introduced by a 5-µm GDVN nozzle jet (GDD = 0.20 fs 2 if the GVD of water is taken to be 39 fs 2 /mm at 600 nm [79]) is on the order of a few attoseconds, so not large for our purposes. Additional stretching, however, becomes significant with delivery methods that introduce more liquid in the beam path, such as the drop-on-demand injector [80]: this produces droplets approximately 170-µm in diameter that are expected to double the duration of a 3-fs pulse (Figure 9-table). The stretching of the optical pump is clearly relevant for visible light absorbing protein crystals and will therefore have to be evaluated depending on the experimental conditions (such as the distribution of crystal sizes and the injection method available). The same considerations apply to third-order dispersion (TOD), which can contribute to stretching as much as GVD at such short pulse durations.

Final pulse duration (fs)
Stretching for a 3 fs pulse  Additionally, the group velocity of the optical pump pulse will differ from that of the X-ray probe in a protein crystal: the X-rays will propagate with a refractive index of 1 but a larger refractive index is associated with the optical pulse (≈ 1.2-1.4 for proteins at crystal growth concentrations [81,82]), resulting in group velocity mismatch (GVM) between the two. Larger crystals will increase this effect, as they introduce a larger path length for the pulses. Figure 10 plots the expected added delay of the optical pump for a distribution of crystal thicknesses at different possible protein crystal refractive indexes: for a 2-µm crystal and n = 1.3, for example, the achievable time resolution cannot be higher than 2 fs, as this is the amount of time by which the pump will be slowed down compared to the probe pulse. The added delay is calculated by multiplying the GVM with the crystal thickness. The GVM is given by the pump and probe group velocities (v pump and v probe ) in the crystal and therefore by their refractive indexes: A final point to consider is the self-focusing that occurs if the pulse power exceeds a critical power, P cr . The nonlinearity effects that result from self-focusing can lead to a total collapse of the beam. While for CW lasers there is a mathematical relationship to calculate P cr from the nonlinear refractive index of a medium, this cannot be used accurately for ultrashort pulses [83] so we cannot predict a specific threshold for self-focusing. This effect is, however, observed regularly with femtosecond pulsed lasers [84][85][86][87] and is expected at the gigawatt peak powers reached by lasers generated through hollow core fiber compression. In a successful experiment, these incredibly high peak powers will have to be contained to avoid self-focusing, beam filamentation, as well as nonlinear absorption effects in the sample.

Materials and Methods
Simulations for the SwissFEL pulse were done using GENESIS [35,40]. Ray-tracing simulations were done using the SIM-MX software [57].

Conclusions
The operation of a polychromatic XFEL source with new chirping capabilities introduces the possibility of performing SFX experiments with sub-femtosecond time resolution via a time-to-energy mapping of the protein diffraction patterns. This could further our understanding of protein dynamics by allowing, for example, visualization of conical intersections, or be applied to investigate radiation damage in protein crystals if the X-rays are used as both the probe and as a damage-inducing pump. A number of technical obstacles lie in the way of this theoretical result, yet ongoing improvements in laser and beamline technologies suggest that the experiment presented here has high chances of being possible in the near future. The ability to maximize the X-ray frequency chirp and to accurately time the pump and probe pulses will be especially relevant in this regard. New algorithms that accept large bandwidths in processing of serial crystal data (such as the PinkIndexer library implemented in CrystFEL [88]) could then be applied for the indexing and wavelength assignment of the diffraction images generated.