Submillisecond X-ray photon correlation spectroscopy from a pixel array detector with fast dual gating and no readout dead-time

Zhang, Q.; Dufresne, E.M.; Grybos, P.; Kmon, P.; Maj, P.; Narayanan, S.; Deptuch, G.W.; Szczygiel, R.; Sandy, A.

doi:10.1107/S1600577516005166

research papers

JOURNAL OF
SYNCHROTRON
RADIATION

ISSN: 1600-5775

Volume 23| Part 3| May 2016| Pages 679-684

https://doi.org/10.1107/S1600577516005166

Open

access

Submillisecond X-ray photon correlation spectroscopy from a pixel array detector with fast dual gating and no readout dead-time

Qingteng Zhang,^a ^* Eric M. Dufresne,^a ^* Pawel Grybos,^b Piotr Kmon,^b Piotr Maj,^b Suresh Narayanan,^a Grzegorz W. Deptuch,^c Robert Szczygiel ^b and Alec Sandy ^a

^aAdvanced Photon Source, Argonne National Laboratory, Lemont, IL 60439, USA, ^bAGH University of Science and Technology, av. Mickiewicza 30, Krakow 30-059, Poland, and ^cFermi National Accelerator Laboratory, Batavia, IL 60510, USA
^*Correspondence e-mail: qzhang234@aps.anl.gov, dufresne@anl.gov

(Received 16 October 2015; accepted 26 March 2016; online 19 April 2016)

Small-angle scattering X-ray photon correlation spectroscopy (XPCS) studies were performed using a novel photon-counting pixel array detector with dual counters for each pixel. Each counter can be read out independently from the other to ensure there is no readout dead-time between the neighboring frames. A maximum frame rate of 11.8 kHz was achieved. Results on test samples show good agreement with simple diffusion. The potential of extending the time resolution of XPCS beyond the limit set by the detector frame rate using dual counters is also discussed.

Keywords: pixel array detector (PAD); high frame rate; photon counting; two counters; X-ray photon correlation spectroscopy (XPCS); small-angle X-ray scattering (SAXS).

1. Introduction

The prevalence of high intensity, partially coherent synchrotron X-ray sources worldwide has led to great advances in the development of coherent X-ray scattering techniques, including transmission and Bragg geometry X-ray ptychography (Rodenburg et al., 2007 ; Hruszkewycz et al., 2012 ), coherent diffraction imaging (Robinson & Harder, 2009 ) and X-ray photon correlation spectroscopy (XPCS) (Dierker et al., 1995 ; Mochrie et al., 1997 ). XPCS measures fluctuations at the mesoscale via changes in the coherent scattering patterns from samples. XPCS is derived from the well established dynamical light scattering (DLS) technique, where a laser beam scatters from fluctuations in the dielectric constant of the sample and generates a scattering pattern with sharp angular intensity variations (a so-called `speckle' pattern) due to interference of the scattered radiation. Dynamics in the system are measured from the temporal decorrelation of the speckle patterns. As an extension of DLS into the X-ray regime, XPCS benefits from both the ability of X-ray beams to (weakly) scatter from materials that are opaque to visible light (Sutton et al., 1991 ) and to probe length scales as small as 1 Å (Leitner et al., 2009 ; Ruta et al., 2014 ). Access to a wide dynamic range in delay times is also critically important and has enabled investigations of novel length-scale-dependent dynamics in soft condensed matter (Falus et al., 2006 ; Lu et al., 2008 ). Studies of dynamics in various systems all benefit from rapid, high fidelity recording of fluctuating speckle patterns and this need will become increasingly important with proposed multibend achromat upgrades of the European Synchrotron Radiation Facility (ESRF) and Advanced Photon Source (APS) set to increase the time-averaged coherent flux by up to two orders of magnitude.

For multi-speckle XPCS, the time resolution is most typically limited by the detector frame rate, as the frame readout time (the reciprocal of frame rate) is the minimum time separation between the detector frames (Fig. 1a). Currently, the fastest `few' X-ray-sensitive frame rates are provided by the AGIPD detector (Becker et al., 2011 , 2013 ; Schwandt et al., 2013 ) and the Keck-PAD (Philipp et al., 2016 ). Both of these devices are integrating detectors with the former being developed for the European X-ray Free-Electron Laser (XFEL) and the latter being developed for both synchrotron and FEL applications. The AGIPD is capable of operating at a 4.7 MHz frame rate whereas the Keck-PAD can operate at up to 10 MHz. In both cases, however, only a limited number of frames can be collected before a longer delay for readout and storage which results in too limited a dynamic range for typical XPCS measurements at a synchrotron source. To date, the fastest continuous frame rate tested for multi-speckle XPCS is 22 kHz from the research version of the Eiger detector (Johnson et al., 2012 ). Though very fast, the Eiger detector has one counter per pixel and 5 µs of readout dead-time between frames. The dead-time limitation can be removed and also more novel acquisition strategies can be implemented when two (or more) counters are available for each pixel. For instance, both the newly developed Pixirad (Bellazzini et al., 2013 , 2015 ) and Lambda (Pennicard et al., 2012 , 2013 ) detectors have dual counters assigned to each pixel that can be assigned to provide either `color' vision when different thresholds are applied to each pixel or to provide nearly readout dead-time-free operation when the counters are alternately acquiring or digitizing data (Fig. 1b). In the readout dead-time-free operation mode, the throughput for the Pixirad is 40 Mpix s⁻¹, whereas for the Lambda it is 1.5 Gpix s⁻¹. Another novel capability provided by dual counters is, with sufficiently fast gating and sensor response dead-time, the ability to resolve two closely spaced signals in time, i.e. recording the sample scattering immediately before and after a laser impulse. This capability is provided by a prototype detector developed by Voxtel (Ross et al., 2016 ). Very recently, the applicability of this detector for short delay time XPCS has also been explored (Dufresne et al., 2016 ).

Figure 1
Time infrastructure of (a) continuous acquisition with one counter per pixel, (b) dual counter acquisition discussed in this study with no readout dead-time between the frames and (c) future upgrade of dual counter acquisition where the separation between the frames is smaller than the digitization time associated with each counter.

Here we present XPCS results obtained using the new Ultra Fast X-ray Camera 32k (UFXC32k) which is a single-photon-counting PAD with dual counters that can operate at a frame rate as high as 11.8 kHz (0.4 Gpix s⁻¹). The dual counters can be read out alternately so there is essentially no readout dead-time between neighboring frames. Though the total pixel throughput is not yet as high as that for the Lambda, this detector provides the shortest delay times of any dual counter detector and with anticipated near-term upgrades (see below) will provide a total pixel throughput comparable with the Lambda and a minimum delay time approximately a factor of two better than was achieved with the research version of the Eiger. To test the applicability of this detector for fast time-resolved coherent scattering, XPCS measurements were performed on latex nanoparticles suspended in glycerol/water and the results were in good agreement with simple diffusion. Our work demonstrates the capacity of the UFXC32k to perform XPCS experiments with very high frame rates and no readout dead-time.

2. Instrumentation

The detector has two parts: a reverse-biased pixelated silicon PN diode as the sensor (128 × 256 square-shaped pixels, pixel size of 75 µm and thickness of 320 µm) and a readout integrated circuit with built-in pixel architecture (see Fig. 2a). Each pixel on the PN diode is bump-bonded to its own dedicated electronics. X-ray photons arrive at the depleted region of the PN diode and generate electron–hole pairs. Fig. 3 is the schematic of the readout circuit of a single pixel. The current signal from the collected electron–hole pairs under the reverse bias voltage is integrated into a charge-sensitive preamplifier (CSA). The preamplifier outputs a voltage pulse with its amplitude proportional to the total charge generated in the pixel. The voltage pulse is then fed to the main amplifier (SHAPER), which reshapes the pulse based on the specifications of noise filtering and timing. If the amplitude of the reshaped pulse is above the preset threshold (TH_SET_L and TH_SET_H), a photon count is recorded in the counter. With respect to charge sharing, like other PADs, during typical operation the threshold is set to half the mean pulse amplitude of a pixel-isolated event. In this way, all photon hits are recorded except for those near the vertex of four pixels (the slope = 1 on the left-hand side of Fig. 4 indicates a very small number of lost events). With the use of trim digital-to-analog converters in each pixel, the effective threshold spread at the discriminator input is only 8.5 e⁻ r.m.s. For the nominal gain of about 50.3 µV e⁻ the measured gain spread is only 1.9% (standard deviation/mean). As a result, no flat-field correction is required, which is not typically the case for PAD detectors (Maj et al., 2014 ; Grybos et al., 2016 ). The measured noise for the nominal settings is 123 e⁻ r.m.s.

Figure 2
(a) Illustration of the pixelated architecture of the Si PN diode sensor and the readout chip of UFXC32k. (b) Photo of UFXC32k.

Figure 3
Schematic of the readout circuit of a single pixel.

Figure 4
Dead-time measurement for about 1200 pixels of UFXC32k.

The dead-time of a detector pixel is defined as the time window after a registered photon event within which all arrival photon events will be ignored (Müller, 1991 ). The pixel dead-time of UFXC32k is calibrated using a high-power copper-target rotating-anode X-ray generator. Fig. 4 shows the relation between the count rate measured by the detector and the rate of the X-ray photons impinging on the detector. The ratio between the measured and the incoming X-ray flux is nearly 1 for all pixels up to 10⁶ counts per pixel per second, which is 10³ times stronger than even the strongest scattering intensity in this study. The dead-time is evaluated using the paralyzable dead-time model (Müller, 1991; Walko et al., 2010 ):

$[{R}_{\rm meas} = {R}_{\rm in} \exp\left(-R_{\rm in}\tau_{\rm P}\right).\eqno(1)]$

Here, R_meas is the measured count rate, R_in is the averaged incoming X-ray pulse rate and $[{\tau _{\rm P}}]$ is the pixel dead-time. The $[{\tau_{\rm P}}]$ of the UFXC32k detector has a mean of 85.1 ns and a standard deviation of 7.08 ns. So far the shortest pixel dead-time reported from a PAD is 67 ns measured on the integrated circuit of the PILATUS3 with instant retrigger capability and a pixel size of 172 µm × 172 µm (Loeliger et al., 2012 ). In the UFXC32k detector the dead-time can be as low as 85 ns and, given that the pixel area is 5.2 times smaller (75 µm × 75 µm), the count rate per detector unit area is significantly higher than that of the PILATUS3.

During the data readout phase, the counters in each column of the pixel array form a shift register. The data from the register is loaded bit by bit into the peripheral fast 128-bit registers and shifted out of the UFXC32k chip via eight LVDS parallel lines. Each pixel has two discriminators (DISCR_L and DISCR_H with thresholds) and two 14-bit ripple counters (COUNTER_L and COUNTER_H). The counters in a pixel can operate in two different modes. The first one is the dual-counter mode, in which each counter is connected to a different discriminator. The level of thresholds (TH_SET_L and TH_SET_H) can be set to different values to allow for energy discrimination or the same value to perform fast acquisition with no readout dead-time. The second one is the long counter mode, in which the two counters are linked serially to form one 28-bit long counter for measurements that require very large dynamic range. The readout dead-time of a single frame t_ro can be expressed in the following equation:

$[{t}_{\rm ro} = {{{N}_{\rm pix}\,{N}_{\rm bit}}\over{{N}_{\rm out}\,{f}_{\rm clk}}}+{t}_{\rm ctrl}.\eqno (2)]$

Here, N_pix = 128 × 256 = 32768 is the total number of pixels, N_bit is the bit depth of the pixel (can be switched between 2, 4, 8 or 14 bits), N_out = 8 is the number of LVDS parallel lines, f_clk = 100 MHz is the clock frequency of the circuit and t_ctrl = 2.56 µs is the additional time required for starting the readout sequence at 100 MHz clock frequency. During the experiment, N_bit is set to 2 (the dynamic range of the detector is from zero to 3) to achieve the minimum t_ro of 84.48 µs and an ultrafast frame rate of 11.8 kHz. The probability of a pixel receiving four photons in the region of interest with the strongest scattering intensity is measured to be 5.4 ×10^{- 7} in this study, which is equivalent to 0.017 pixels per frame. We therefore do not expect significant truncation of the scattering intensities. The acquisition time t₀ for the detector frame is set to be equal to t_ro to maximize the counting statistics. As long as t₀ is greater than 2 µs for both counters, the switching time between counters is shorter than the bunch separation time at the APS (153 ns), leading to essentially no readout dead-time between the neighboring frames (Fig. 1b). A more detailed list of detector parameters is given by Grybos et al. (2016).

3. Experiment and discussion

Small-angle X-ray scattering (SAXS) measurements were performed at the experimental station 8-ID-I of the APS. Two Undulator As were set to produce 7.4 keV X-rays. A water-cooled double-crystal Ge(111) monochromator selected a monochromatic beam with a bandpass of about 0.03% and a 20 µm × 150 µm (H×V) slit was used to select the coherent part of the beam, which had a flux of approximately 2 × 10¹⁰ photons s⁻¹. The beam was then vertically focused at the sample with a full width at half-maximum (FWHM) of about 3 µm using 10 one-dimensional parabolic compound refractive X-ray lenses made of beryllium with a radius of curvature of 0.2 mm (RX Optics). Guard slits between the lens and the sample were used to remove background scattering from both the lenses and the blades on the main collimating slit. A tungsten beamstop immediately upstream from the detector (4 m downstream of the sample) blocked the direct beam.

The sample is 4% (volume fraction) 67 nm-radius latex nanoparticles dispersed in glycerol (Lurio et al., 2000 ). The solvent may contain traces of water. The sample was sealed in 1 mm-diameter glass capillaries, integrated into the vacuum of the experiment setup and measured at room temperature. During each acquisition sequence, a series of 100000 frames of coherent scattering patterns are collected from the sample in about 8.5 s. The pixels were converted into coordinates in reciprocal space based on the sample-to-detector distance, the photon energy and the transverse position of the detector pixels relative to the direct beam. A total of nine acquisition sequences were performed to improve the correlation function statistics. Fig. 5(a) shows the time-averaged two-dimensional small-angle X-ray scattering pattern from the sample. The number of dead pixels is only 10 out of a total of 32768 pixels. Fig. 5(b) is the azimuthal average of Fig. 5(a) plotted as a function of q. The average radius of the particles calculated from the spacing between the fringes is 66.7 nm, which is very close to the expected value of 67 nm.

Figure 5
(a) Time-averaged scattering from the latex nanoparticle suspension. The scattering intensity is indicated by the logarithmic color bar. (b) Azimuthal average of Fig. 5

(a).

The dynamics of the latex nanoparticles are determined from the correlation coefficient g₂(τ) between speckle patterns separated by a time τ, which is defined as:

$[{g}_{2}(q,\tau) = {{{\left\langle I(q,t)\,I(q,t+\tau)\right\rangle }_{t}}\over{{\left\langle I(q,t)\right\rangle }_{t}}}. \eqno(3)]$

Here, I(q,t) is the photon count of a pixel at wavevector q in the frame collected at time t and g₂ is evaluated within the multi-τ framework (Cipelletti & Weitz, 1999 ) for improved efficiency of the calculation. The time average $[\langle \cdots\rangle _{\rm t}]$ goes over all pairs of frames with a temporal separation of τ. The statistics of g₂ are improved by first azimuthally averaging over pixels with similar q to calculate g₂ for each data set and then averaging the g₂ values from the nine data sets. The correlation time τ₀(q) can be extracted from fitting the following equation to the measured correlation decays:

$[{g}_{2}(q,\tau) = 1+\beta \exp\left\{2\left[-{{\tau }\over{{\tau }_{0}(q)}}\right]\right\}. \eqno(4)]$

Here, β is the instrumental beamline contrast that depends on the pixel size, speckle size and beamline factors, and is measured to be about 10% for this experiment. It is clear from equation (4) that g₂ decays to the baseline of 1 at long delay times. Fig. 6(a) shows g₂(τ) calculated using equation (3) at a few different values of q. The baseline (at 1) has been subtracted and the resulting correlation decays have been normalized by the short-time contrast. It can be seen from Fig. 6(a) that the sample decorrelates faster at smaller length scales (larger q values).

Figure 6
(a) Dynamics of latex nanoparticles indicated by g₂(τ) at different q. (b) Decorrelation time τ(q) versus q. The red line shows the inverse-square decay of the correlation time.

For particles undergoing Brownian motion, the correlation time and the wavevector q at which it is measured satisfy the inverse-square relation:

$[{\tau }_{0}(q) = {{1}\over{D{q}^{2}}},\eqno(5)]$

where D is the diffusion coefficient. The q dependence in Fig. 6(b) is fitted using an inverse-square function. The good agreement between experimental results and the prediction for simple diffusion demonstrates the capacity of the UFXC32k detector to capture fast dynamics with high fidelity. Further comparison between the calculated and the measured diffusion constant is not possible (nor the point of this paper) owing to the ambiguity of the solvent viscosity from traces of water mixed in the glycerol.

4. Conclusion

The UFXC32k is a photon-counting dual-counter PAD with a very fast frame rate of 11.8 kHz. The dual-counter setup can be used for energy discrimination, expansion of the dynamic range or frame acquisition with zero readout dead-time (Fig. 1b). In addition, both the length and position of the gating signal for dual counters are fully adjustable, meaning that the minimum temporal separation between acquisitions is not restricted by the digitization time of one of the counters (Fig. 1c). In this case, XPCS studied using PADs with time resolution nearing 1 µs is possible given enough coherent flux (Dufresne et al., 2016). Future upgrades will include: (i) increasing the frame rate to 47.2 kHz by using a 400 MHz clock signal; (ii) reducing the data transfer time and increasing the largest number of frames allowed in each acquisition through sparsification (saving only pixels with non-zero readouts). With anticipated increases in coherent flux from synchrotron X-ray sources becoming a global trend, i.e. the recent start-ups of Petra-III and NSLS-II and the expected future upgrades of the ESRF and APS, development and implementation of fast frame rate, dual counter detectors such as the UFXC32k will surely help expand the capacity of XPCS and open doors to considerable scientific opportunities.

Acknowledgements

This research was performed on beamline 8-ID-I of the Advanced Photon Source, a US Department of Energy (DOE) Office of Science User Facility operated for the DOE Office of Science by Argonne National Laboratory under contract No. DE-AC02-06CH11357. AGH University of Science and Technology was supported by the National Center for Research and Development, Poland, PBS1/A3/12/2012, in the years 2012–2015. Fermilab is operated by Fermi Research Alliance, LLC, under contract No. DE-AC02-07CH11359 with the US Department of Energy. We wish to thank the ASIC Development Group at Fermilab for technical support. QZ acknowledges help with the sample preparation from Professor Laurence Lurio of Northern Illinois University.

References

Becker, J. et al. (2013). J. Instrum. 8, C01042. Google Scholar
Becker, J., Gutt, C. & Graafsma, H. (2011). J. Instrum. 6, P11005. Google Scholar
Bellazzini, R., Brez, A., Spandre, G., Minuti, M., Pinchera, M., Delogu, P., de Ruvo, P. L. & Vincenzi, A. (2015). J. Instrum. 10, C01032. CrossRef Google Scholar
Bellazzini, R., Spandre, G., Brez, A., Minuti, M., Pinchera, M. & Mozzo, P. (2013). J. Instrum. 8, C02028. CrossRef Google Scholar
Cipelletti, L. & Weitz, D. A. (1999). Rev. Sci. Instrum. 70, 3214–3219. CrossRef CAS Google Scholar
Dierker, S., Pindak, R., Fleming, R. M., Robinson, I. K. & Berman, L. E. (1995). Phys. Rev. Lett. 75, 449–452. CrossRef PubMed CAS Google Scholar
Dufresne, E. M., Narayanan, S., Sandy, A. R., Kline, D. M., Zhang, Q., Landahl, E. C. & Ross, S. (2016). Opt. Express, 24, 355–364. CrossRef PubMed Google Scholar
Falus, P., Borthwick, M. A., Narayanan, S., Sandy, A. & Mochrie, S. (2006). Phys. Rev. Lett. 97, 066102. CrossRef PubMed Google Scholar
Grybos, P., Kmon, P., Maj, P. & Szczygiel, R. (2016). IEEE Trans. Nucl. Sci. 63, DOI:10.1109/TNS.2016.2523260. Google Scholar
Hruszkewycz, S. O., Holt, M. V., Murray, C. E., Bruley, J., Holt, J., Tripathi, A., Shpyrko, O. G., McNulty, I., Highland, M. J. & Fuoss, P. H. (2012). Nano Lett. 12, 5148–5154. Web of Science CrossRef CAS PubMed Google Scholar
Johnson, I., Bergamaschi, A., Buitenhuis, J., Dinapoli, R., Greiffenberg, D., Henrich, B., Ikonen, T., Meier, G., Menzel, A., Mozzanica, A., Radicci, V., Satapathy, D. K., Schmitt, B. & Shi, X. (2012). J. Synchrotron Rad. 19, 1001–1005. Web of Science CrossRef CAS IUCr Journals Google Scholar
Leitner, M., Sepiol, B., Stadler, L.-M., Pfau, B. & Vogl, G. (2009). Nat. Mater. 8, 717–720. Web of Science CrossRef PubMed CAS Google Scholar
Loeliger, T., Broennimann, C., Donath, T., Schneebeli, M., Schnyder, R. & Trueb, P. (2012). Proceedings of the 2012 IEEE Nuclear Science Symposium and Medical Imaging Conference, 27 October 2012–3 November 2012, Anaheim, California, USA, pp. 610–615. Google Scholar
Lu, X., Mochrie, S. G. J., Narayanan, S., Sandy, A. R. & Sprung, M. (2008). Phys. Rev. Lett. 100, 045701. Web of Science CrossRef PubMed Google Scholar
Lurio, L. B., Lumma, D., Sandy, A. R., Borthwick, M. A., Falus, P., Mochrie, S. G. J., Pelletier, J. F., Sutton, M., Regan, L., Malik, A. & Stephenson, G. B. (2000). Phys. Rev. Lett. 84, 785–788. Web of Science CrossRef PubMed CAS Google Scholar
Maj, P., Grybos, P., Kmon, P. & Szczygiel, R. (2014). Proceedings of the 40th European Solid-State Circuit Conference (ESSCIRC), pp. 147–150. Google Scholar
Mochrie, S. G. J., Mayes, A. M., Sandy, A. R., Sutton, M., Brauer, S., Stephenson, G. B., Abernathy, D. L. & Grübel, G. (1997). Phys. Rev. Lett. 78, 1275–1278. CrossRef CAS Web of Science Google Scholar
Müller, J. W. (1991). Nucl. Instrum. Methods Phys. Res. A, 301, 543–551. Google Scholar
Pennicard, D., Lange, S., Smoljanin, S., Hirsemann, H. & Graafsma, H. (2012). J. Instrum. 7, C11009. Web of Science CrossRef Google Scholar
Pennicard, D., Lange, S., Smoljanin, S., Hirsemann, H., Graafsma, H., Epple, M., Zuvic, M., Lampert, M. O., Fritzsch, T. & Rothermund, M. (2013). J. Phys. Conf. Ser. 425, 062010. CrossRef Google Scholar
Philipp, H. T., Tate, M. W., Purohit, P., Shanks, K. S., Weiss, J. T. & Gruner, S. M. (2016). J. Synchrotron Rad. 23, 395–403. CrossRef CAS IUCr Journals Google Scholar
Robinson, I. & Harder, R. (2009). Nat. Mater. 8, 291–298. Web of Science CrossRef PubMed CAS Google Scholar
Rodenburg, J., Hurst, A., Cullis, A., Dobson, B., Pfeiffer, F., Bunk, O., David, C., Jefimovs, K. & Johnson, I. (2007). Phys. Rev. Lett. 98, 034801. Web of Science CrossRef PubMed Google Scholar
Ross, S., Haji-Sheikh, M., Huntington, A., Kline, D., Lee, A., Li, Y., Rhee, J., Tarpley, M., Walko, D. A., Westberg, G., Williams, G., Zou, H. & Landahl, E. (2016). J. Synchrotron Rad. 23, 196–205. CrossRef CAS IUCr Journals Google Scholar
Ruta, B., Baldi, G., Chushkin, Y., Rufflé, B., Cristofolini, L., Fontana, A., Zanatta, M. & Nazzani, F. (2014). Nat. Commun. 5, 3939. CrossRef PubMed Google Scholar
Schwandt, J., Fretwurst, E., Klanner, R., Poehlsen, T. & Zhang, J. (2013). Nucl. Instrum. Methods Phys. Res. A, 731, 252–254. CrossRef CAS Google Scholar
Sutton, M., Mochrie, S. G. J., Greytak, T., Nagler, S. E., Berman, L. E., Held, G. A. & Stephenson, G. B. (1991). Nature (London), 352, 608–610. CrossRef Web of Science Google Scholar
Walko, D. A., Arms, D. A., Dufresne, E. M., Landahl, E. C., Garrett, R., Gentle, I., Nugent, K. & Wilkins, S. (2010). AIP Conf. Proc. 1234, 856–859. CrossRef Google Scholar

This is an open-access article distributed under the terms of the Creative Commons Attribution (CC-BY) Licence, which permits unrestricted use, distribution, and reproduction in any medium, provided the original authors and source are cited.