research papers\(\def\hfill{\hskip 5em}\def\hfil{\hskip 3em}\def\eqno#1{\hfil {#1}}\)

Journal logoJOURNAL OF
SYNCHROTRON
RADIATION
ISSN: 1600-5775

Optimizing the signal-to-noise ratio for X-ray photon correlation spectroscopy

CROSSMARK_Color_square_no_text.svg

aInstitute Laue Langevin, Grenoble F-38042, France, bDepartment of Physics, Northern Illinois University, DeKalb, IL 60115, USA, and cDepartments of Physics and Applied Physics, Yale University, New Haven, CT 06520, USA
*Correspondence e-mail: falus@ill.fr

(Received 8 June 2005; accepted 23 February 2006)

An analysis is presented of how to optimize the experimental beamline configuration for achieving the best possible signal-to-noise ratio (SNR) in X-ray photon correlation spectroscopy experiments using area detectors. It is shown that there exists an optimum detector distance; namely, the highest SNR is achieved by matching the angular pixel size with the angular source size. Binning several pixels together can increase the SNR by permitting to match the shape of a detector pixel to the shape of the source. It is also shown that collimating slits several times wider than the effective transverse coherence length are optimal; further, it is demonstrated that the energy dependence of the SNR is dictated by the energy dependence of detector efficiency and source brilliance. Ultimately the effects of focusing and low longitudinal coherence are discussed.

1. Introduction

X-ray photon correlation spectroscopy is an emerging new method that can now be used at high-brilliance synchrotron beamlines for characterizing the dynamics of slowly moving condensed matter (Falus et al., 2005[Falus, P., Borthwick, M. A. & Mochrie, S. G. J. (2005). Phys. Rev. Lett. 94, 016105.]; Sutton et al., 1991[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.]; Brauer et al., 1995[Brauer, S., Stephenson, G. B., Sutton, M., Brüning, R., Dufresne, E., Mochrie, S. G. J., Grübel, G., Als-Nielsen, J. & Abernathy, D. L. (1995). Phys. Rev. Lett. 74, 2010.]; Dierker et al., 1995[Dierker, S. B., Pindak, R., Fleming, R. M., Robinson, I. K. & Berman, L. E. (1995). Phys. Rev. Lett. 75, 449-452.]; Thurn-Albrecht et al., 1996[Thurn-Albrecht, T., Steffen, W., Patkowski, A., Meier, G., Fisher, E. W., Grübel, G. & Abernathy, D. L. (1996). Phys. Rev. Lett. 77, 5437-5440.]; Mochrie et al., 1997[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.]; Tsui & Mochrie, 1998[Tsui, O. K. C. & Mochrie, S. G. J. (1998). Phys. Rev. E, 57, 2030.]; Thurn-Albrecht et al., 1999[Thurn-Albrecht, T., Meier, G., Müller-Buschbaum, P., Patkowski, A., Steffen, W., Grübel, G., Abernathy, D. L., Diat, O., Winter, M., Koch, M. G. & Reetz, M. T. (1999). Phys. Rev. E, 59, 642-649.]; Price et al., 1999[Price, A. C., Sorensen, L. B., Kevan, S. D., Toner, J., Poniewierski, A. & Holyst, R. (1999). Phys. Rev. Lett. 82, 755.]; Lurio et al., 2000[Lurio, L. B., Lumma, D., Falus, P., Borthwick, M. A., Mochrie, S. G. J., Pelletier, J.-F., Sutton, M., Malik, A. & Stephenson, G. B. (2000). Phys. Rev. Lett. 84, 785.]; Lumma, Lurio, Sandy et al., 2000[Lumma, D., Lurio, L. B., Sandy, A. R., Borthwick, M., Falus, P., Mochrie, S. G. J., Sutton, M. & Regan, L. (2000). Phys. Rev. E, 62, 8258.]; Lumma et al., 2001[Lumma, D., Borthwick, M. A., Falus, P., Lurio, L. B. & Mochrie, S. G. J. (2001). Phys. Rev. Lett. 86, 2042.]; Seydel et al., 2001[Seydel, T., Madsen, A., Tolan, M., Grubel, G. & Press, W. (2001). Phys. Rev. B, 63, 073409.]; Livet et al., 2001[Livet, F., Bley, F., Caudron, R., Geissler, E., Abernathy, D., Detlefs, C., Grubel, G. & Sutton, M. (2001). Phys. Rev. E, 63, 036108.]; Dufresne et al., 2002[Dufresne, E. M., Nurushev, T., Clarke, R. & Dierker, S. B. (2002). Phys. Rev. E, 65, 061507.]; Sikharulidze et al., 2002[Sikharulidze, I., Dolbnya, I. P., Fera, A., Madsen, A., Ostrovskii, B. I. & de Jeu, W. H. (2002). Phys. Rev. Lett. 88, 115503.]; Kim et al., 2003[Kim, H. J., Ruehm, A., Lurio, L. B., Basu, J., Lal, J., Lumma, D., Mochrie, S. G. J. & Sinha, S. K. (2003). Phys. Rev. Lett. 90, 068302.]; Madsen et al., 2003[Madsen, A., Als-Nielsen, J. & Grübel, G. (2003). Phys. Rev. Lett. 90, 085701.]). Typically, to carry out X-ray photon correlation spectroscopy (XPCS) experiments, the X-ray beam from a high-brilliance X-ray source is monochromated and apertured to a size approximately equal to the transverse coherence length of the source. This creates a partially coherent X-ray beam, which illuminates a sample of interest. Density fluctuations within the sample cause scattering, giving rise to a time-varying X-ray speckle pattern, which is captured moment-by-moment by a detector. The two-time correlation of the intensities in each pixel yields the intensity auto-correlation function, g2, which shows the characteristic relaxation times of the sample. Specifically, g2 is related to the sample's intermediate scattering function (ISF), S(Q, t), via

[g_{2}=1+A[S(Q,t)]^2.\eqno(1)]

where A is the optical contrast, which depends on the parameters of the experimental set-up. At present, XPCS is a signal-limited technique and therefore it is important to optimize as far as possible the signal-to-noise ratio (SNR). Minimizing sample damage also calls for the most efficient use of X-rays. Moreover, there are presently several major proposals to construct new very high brilliance X-ray sources: (NSLS2 Brookhaven, ERL Cornell, MIT-Bates XFEL, LCLS Stanford). High-brilliance experiments, of which XPCS is one important example, provide the scientific rationale for these machines. It is essential to understand the ingredients that allow one to optimize the XPCS SNR in order to make informed decisions concerning the layout and design of beamlines that might be destined for XPCS experiments. In this paper we elaborate calculations (Lumma, Lurio, Mochrie & Sutton, 2000[Lumma, D., Lurio, L. B., Mochrie, S. G. J. & Sutton, M. (2000). Rev. Sci. Instrum. 71, 3274.]) of the anticipated SNR in XPCS experiments in order to answer the following questions: what is the optimum aperture before the sample for a given source size? Is there an optimum sample-to-detector distance for given CCD pixel dimensions and given number of pixels? Is it possible to improve the SNR by binning together the intensity in spatially separated pixels, and if so how many? How does the SNR change with X-ray energy?

2. Signal-to-noise ratio

We consider a synchrotron beamline that consists of a partially coherent source, a sample, adjustable collimating slits close to the sample, and an area detector (CCD). Throughout this paper we will use the numerical values relevant to beamline 8-ID at the Advanced Photon Source (APS), where our XPCS experiments have been carried out, and also to a hypothetical XPCS beamline at the proposed NSLS-II synchrotron at Brookhaven National Laboratory. In the case of 8-ID the beam-defining slits are located 65 m from the source and the sample is 1.3 m further downstream of the slits. As long as the sample is in the near field of the slits (in this case if the slits are bigger than 15 µm) there is a negligible difference between the slits size and the beam width at the sample. We denote the horizontal and vertical source size as σx and σz, respectively. At 8-ID we take σx = 220 µm and σz = 15 µm, while at NSLS-II we use σx = 54.5 µm and σz = 3.9 µm. These are the standard deviations of an assumed Gaussian profile for the source intensity, so that the full width at half-maximum of the source intensity is (8ln2)1/2σ ≃ 2.35σ. At 8-ID the source-to-sample distance is R′ = 65 m. For a hypothetical beamline at NSLS-II, we assume R′ = 32 m, simply scaling the source-to-sample distance by the size of the ring. The angular source sizes are defined as θz = (8ln2)1/2σz/R′, θx = (8ln2)1/2σx/R′. The transverse coherence length (Sandy et al., 1999[Sandy, A. R., Lurio, L. B., Mochrie, S. G. J., Malik, A., Stephenson, G. B., Pelletier, J.-F. & Sutton, M. (1999). J. Synchrotron Rad. 6, 1174.]) is defined as

[\xi_{x}={{\lambda{R'}}\over{2\pi\sigma_{x}}}={{(8\ln2)^{1/2}\,}\over{2\pi}}{{\lambda}\over{\theta_{x}}},\eqno(2)]

where λ = 0.185 nm is the X-ray wavelength. The CCD detector consists of pixels of linear size Ux by Uz, and we define the angular size of pixels as ωx = Ux/Rdet and ωz = Uz/Rdet, where Rdet is the slit-to-detector distance. Using these parameters it may be shown (Lumma, Lurio, Mochrie & Sutton, 2000[Lumma, D., Lurio, L. B., Mochrie, S. G. J. & Sutton, M. (2000). Rev. Sci. Instrum. 71, 3274.]; Sandy et al., 1999[Sandy, A. R., Lurio, L. B., Mochrie, S. G. J., Malik, A., Stephenson, G. B., Pelletier, J.-F. & Sutton, M. (1999). J. Synchrotron Rad. 6, 1174.]) that we can introduce the effective transverse coherence length in the horizontal direction as

[\eqalignno{\Xi_{x}&={}{{\xi_{x}}\over{\left[1+\left(U_{x}^{2}R'^{\,2}\right)/\left(8\ln2\,\sigma_{x}^{2}R_{\rm{det}}^{2}\right)\right]^{1/2}}}\cr&={}{{(8\ln2)^{1/2}\,\lambda}\over{2\pi}}{{1}\over{\left(\omega_{x}^{2}+\theta_{x}^{2}\right)^{1/2}}}.&(3)}]

Analogously defined is Ξz, which is the effective coherence length in the z-direction. The effective coherence lengths fold the non-zero detector resolution into the source size (Abernathy et al., 1998[Abernathy, D. L., Grübel, G., Brauer, S., McNulty, I., Stephenson, G. B., Mochrie, S. G. J., Sandy, A. R., Mulders, N. & Sutton, M. (1998). J. Synchrotron Rad. 5, 37-47.]). The longitudinal coherence length is given by

[\Lambda\simeq\lambda\left(E/\Delta{E}\right).\eqno(4)]

For the germanium [111] monochromator currently used at 8-ID, the relative bandwidth is almost energy independent with E/ΔE ≃ 3 × 103. In a small-angle-scattering geometry, the longitudinal coherence length will only effect the speckle contrast when it becomes comparable with the maximum pathlength difference in the sample, [\delta], which is given by

[\eqalignno{\delta={}&2\sin\psi\cos\psi\cos\varphi\,L_x+2\sin\psi\cos\psi\sin\varphi\,L_z\cr&+2\sin^2\psi\,W,&(5)}]

where Lx, Lz and W are the dimensions of the scattering volume in the horizontal, vertical and along-the-beam directions, respectively, ψ is the scattering angle and φ is the azimuthal angle from the horizontal scattering plane. The contrast is reduced by a finite longitudinal coherence length by a factor of approximately

[\exp\left(-\left|\delta\right|/\Lambda\right).\eqno(6)]

Typically, at 8-ID, we have Lx = Lz = 50 µm, while the sample thickness, W, is comparable with the absorption length which depends on the sample material and energy. For a typical hydrocarbon studied using 6.7 keV X-rays the attenuation length (and consequently a typical sample thickness) is ∼1.6 mm. Choosing Q = 0.2 nm−1 gives δ/Λ ≃ 0.1, resulting in a less than 10% correction to the contrast for the typical case of Q < 0.2 nm−1. The effects of a finite longitudinal coherence length will be discussed in more detail in §6[link].

Further, following the calculation of Lumma, Lurio, Mochrie & Sutton (2000[Lumma, D., Lurio, L. B., Mochrie, S. G. J. & Sutton, M. (2000). Rev. Sci. Instrum. 71, 3274.]) and ignoring, for now, the effects of longitudinal coherence we can write the optical contrast (A) of the speckle patterns as

[A=F\left({{L_{x}}/{\Xi_{x}}}\right)F\left({{L_{z}}/{\Xi_{z}}}\right),\eqno(7)]

where Lx and Lz denote the horizontal and vertical slit widths, respectively, and F is the contrast function, given by

[F(x)=\left(1/x^2\right)\left[x\,\pi^{1/2}\,{\rm{erf}}\left(x\right)+\exp\left(-x^2\right)-1\right],\eqno(8)]

where erf(…) is the error function. It will be useful to note that, for x [\gg] 1, F(x) ≃ π1/2/x, and that, for x [\ll] 1, F(x) ≃ 1 − x2/6.

In considering the SNR, the low-count-rate limit is the appropriate limit for XPCS. Jakeman has shown that, in this limit, the variance of g2 is given by

[{\rm{var}}\,g_{2}=g_{2}/\left(M{\bar{n}}^{2}\right),\eqno(9)]

where [\bar{n}] is the mean number of counts detected per accumulation time and M is the total number of correlated pairs averaged to yield g2 (Jakeman, 1973[Jakeman, E. (1973). Photon Correlation and Light Beating Spectroscopy, edited by H. Z. Cummins and E. R. Pike, pp. 75-149. New York: Plenum.]). This may be re-written in terms of the count rate ([\bar{I}]) per CCD pixel, the accumulation time (τ), the number of pixels (nx × nz) and the total experimental duration (T) as

[{\rm{var}}\,g_{2}=g_{2}/\left(n_xn_zT\tau{\bar{I}}^{\,2}\right).\eqno(10)]

A sensible measure of the SNR is

[R_{\rm{sn}}=\left(g_2-1\right)/\left({\rm{var}}\,g_2\right)^{1/2}=\left(T{\tau}\,n_xn_z/g_2\right)^{1/2}\,\,{\bar{I}}\left(g_{2}-1\right).\eqno(11)]

A simplified expression, suitable for comparing different experimental parameters, is obtained by combining (10)[link] and (11)[link], substituting g2 with its zero time limit 1 + A, and assuming that g21/2 ≃ 1, i.e. that the contrast is small,

[R_{\rm{sn}}=A{\bar{I}}\left(T{\tau}\,n_xn_z\right)^{1/2}.\eqno(12)]

In fact, we define (12)[link] to be the SNR in XPCS experiments and seek to maximize it. Equation (12)[link] shows that for XPCS the SNR is linearly proportional to the scattering intensity, in contrast to ordinary intensity measurements where the SNR is proportional to the square root of the intensity. On the other hand, the XPCS SNR is proportional to the square root of the experimental duration (T), just as for scattering experiments.

In an XPCS experiment the goal is to measure g2(Q, t) (and thus the ISF) as a function of time (t) and wavevector (Q). With an area detector it is possible to acquire data at many different wavevectors simultaneously. However, data with insufficient SNR at many different values of Q are worthless. It follows that in assessing the consequences of (12)[link] we should focus on the SNR for a particular value of Q within some desired Q-resolution, not the whole area of the detector. In the context of an experiment, this corresponds to the scattered X-rays that fall within a definite solid-angular range of particular scattering angles. Thus, in (12)[link], the values of nx and ny correspond to this solid angular range. It is therefore convenient to introduce the angle subtended by the Q of interest in the x-direction (Ωx) and in the z-direction (Ωz). In terms of the angular extent of each pixel, ωx and ωz in the x- and z-directions, respectively, and Ωx and Ωz, we may express the number of pixels appearing in (12)[link] as nx = Ωx/ωx and nz = Ωz/ωz. In this paper we will usually assume that nx and nz do not exceed the number of rows and columns available on the detector (or detectors).

The count rate per pixel ([\bar{I}]) depends on the incident flux density (Φ), the solid angle subtended by each pixel, the scattering cross section per unit volume of the sample (Σ), the detector efficiency (K), the illuminated sample area (Lx × Lz) and the sample thickness (W) and attenuation length (Λ) via

[{\bar{I}}=K\Phi{L}_{x}L_{z}\omega_{x}\omega_{z}\Sigma{W}\exp\left(-W/\Lambda\right).\eqno(13)]

[Equation (13)[link] assumes that the flux is constant over the illuminated sample area.] The flux density, in turn, can be expressed in terms of the source brilliance B~, the relative energy bandwidth ΔE/E and the angular source sizes, yielding

[\Phi=\tilde{B}\left(\Delta{E}/E\right)\left[2\pi/\left(8\ln2\right)\right]\,\theta_x\theta_z.\eqno(14)]

Using (2)[link], (7)[link], (12)[link], (13)[link] and (14)[link] we obtain

[R_{\rm{sn}}=K\left(T\tau\,\Omega_x\Omega_z\right)^{1/2}\Sigma{W}\exp\left(-W/\Lambda\right)\tilde{B}(\Delta{E}/E)r_{{\rm{sn}}x}r_{{\rm{sn}}z},\eqno(15)]

where rsnx depends solely on the horizontal pixel size, the horizontal source size and the horizontal slit size. Similarly, rsnz depends solely on the vertical parameters. Specifically,

[r_{{\rm{sn}}x}=F\left(L_x/\Xi_x\right)\,\omega_x^{1/2}\,\theta_{x}L_{x},\eqno(16)]

with an analogous expression for rsnz. The factorization of (15)[link] permits us to optimize the SNR separately for the vertical and the horizontal directions. Equation (15)[link] further informs us that the SNR is linearly proportional to the source brilliance, but proportional to the square root of the accumulation time (τ). This scaling is especially interesting in the context of proposals for more-brilliant X-ray sources. It implies that, for a sample with a given scattering cross section, if one increases the source brilliance by a factor f then one may decrease the accumulation time by a factor of f1/2 to achieve the same SNR. For example, with a 100-fold increase in [\tilde{B}], one may keep the total measurement time constant, and reduce the shortest correlation time τ by (100)2 to achieve an unchanged SNR. For the present, (15)[link] also reflects the important facts that the SNR is linearly proportional to the detector efficiency and the sample scattering strength. The SNR is also proportional to Wexp(−W/Λ), which is maximized by choosing the sample thickness to be equal to the attenuation length.

3. Optimizing the angular pixel size

In this section we explore whether there is an optimal angular pixel size and whether we can improve the SNR by binning several pixels together or by changing the detector-to-sample distance. In these calculations we will consider the source sizes and the slit sizes to be constant. Our single degree of freedom will be the angular pixel size. As we noted in the previous section, we can optimize the SNR separately for the horizontal and vertical direction. Thus, we seek to optimize (16)[link].

Fig. 1[link] plots (16)[link] as a function of ω for several slit sizes. It is apparent from this graph that for large slits (L [\gg] ξ) rsn is independent of L. We explore this case first.

[Figure 1]
Figure 1
SNR versus the ratio of the angular pixel size (ω) and the angular source size (θ) at different slit sizes (L).

3.1. Large slit sizes

For large values of x, F(x) ≃ π1/2/x. This allows us to simplify (16)[link] to become

[r_{\rm{sn}}=\left({{2\ln2\,\lambda^2}\over{\pi}}\right)^{1/2}\,\left[{{1}\over{\theta\left(\theta/\omega+\omega/\theta\right)}}\right]^{1/2}.\eqno(17)]

This informs us first that for large slit sizes the SNR is independent of slit size. The second notable feature is that rsn exhibits a peak at ω = θ, so that there indeed exists an optimum detector distance for a given pixel size, or vice versa. Usually CCD pixel sizes are between 10 and 30 µm; it follows that it may be difficult to achieve the optimum sample-to-detector distance. For example, the optimum sample-to-detector distance for a CCD camera at 8-ID with 14 µm pixels is 24.8 m, dictated by the excellent vertical coherence at the APS. On the other hand the maximum possible sample-to-detector distance in the 8-ID endstation hutch is 7 m. Clearly the SNR of XPCS can only be improved if adequate higher-resolution detectors are developed. Although the most pressing problem of CCD area detectors is speed (Falus et al., 2004[Falus, P., Borthwick, M. A. & Mochrie, S. G. J. (2004). Rev. Sci. Instrum. 75, 4383-4400.]), the SNR can only be optimized if the size of the detector pixels is much smaller then the present sizes. Hopefully, parallel to developing new X-ray sources, new high-efficiency detectors are also being developed.

We should emphasize that Fig. 1[link] is not valid for point detectors. When we have a single pixel, the detector size ΩzΩx is not constant [unlike in (15)[link]]; it scales with the pixel size ωxωz. Then Fig. 1[link] is modified such that for large pixel sizes the SNR is constant; for small ω the SNR will decay linearly with ω. Thus we will have a lower limit on pixel size, but no higher limit. The SNR is optimal for point detectors as long as the detector size is bigger than the speckle size.

It is also interesting to note that to achieve the maximum SNR requires pixels with an aspect ratio that is the same as the aspect ratio of the source. For the usual case of cameras with square pixels, this can be achieved by binning or adding the signal from several pixels in the horizontal direction. These considerations suggest that the optimum experimental set-up at 8-ID would be achieved by moving the camera as far back as possible in the hutch, so as to approach as closely as possible the condition ωz = θz. At this distance, ωx will be smaller than θx but this can be compensated by binning several pixels together until ωx = θx.

To elucidate how important is it to achieve these optimum conditions, it is useful to employ the approximation that rsn[r_{{\rm{sn}}_{\rm{max}}}(2\omega/\theta)^{1/2}], which is approximately valid for θ/ω > 3. Then we may examine the following three 8-ID-based examples:

(i) Using a camera with 14 µm × 14 µm pixels at 1.7 m from the sample, the horizontal size is optimal, the vertical size is larger than optimal by a factor of 14 causing the SNR to be 71/2 = 2.6 smaller than optimal. This is a significant loss of SNR from that optimally possible.

(ii) Using the same camera with 14 µm × 14 µm pixels at 7 m, the vertical pixel size is a factor of 3.5 larger than optimal, while the horizontal pixel size is a factor of four smaller than optimal. It follows that the SNR is [(3.5/2)(4/2)]1/2 ≃ 1.9 less than optimal. This too represents an important loss of SNR.

(iii) Using the same camera with 14 µm × 14 µm pixels at 7 m but now binning four pixels horizontally, the vertical pixel size is a factor of three larger than optimal, while the effective horizontal pixel size (56 µm) is optimal. In this case the SNR is (3.5/2)1/2 = 1.3 times smaller than optimal, which is a considerably better situation than in (i) and (ii).

These considerations are modified in the case where the detector area is not large enough to accept all of the scattering at the wavevector of interest. This is not a fundamental issue, since in principle we may tile together a number of smaller detectors to make an arbitrarily large detector, which would be able to accept all of the scattering at a given wavevector. Nevertheless, for a single detector the finite detector acceptance is of practical relevance. In particular, for isotropic systems, which give rise to circularly symmetric rings of scattering, we can average g2 from all the pixels within rings that correspond to the same wavenumber. If the detector acceptance is not large enough to detect the entire ring, the detected solid angle is reduced by a factor proportional to R−1. Consequently, for high wavenumbers the number of pixels in the ring varies as Rdet, not as Rdet2, resulting in a Rdet-1/2 correction to our calculation or the SNR.

Again taking into account the scaling form of rsnRdet1/2 for ω [\gg] θ, and rsnRdet-1/2 for ω [\ll] θ (large Rdet), noting that the azimuthal correction varies as Rdet-1/2, and noting that in case of optimum binning rsn is constant for ω < θ, we distinguish four cases which are demonstrated in Fig. 2[link]:

[Figure 2]
Figure 2
SNR versus sample-to-detector distance, calculated with parameters appropriate to 8-ID. Specifically, the detector pixel size is 14 µm, the source size is 220 µm × 15 µm, the slit size is 1000 µm and the photon energy is 6.7 keV. For reference, with these parameters θx = ωx occurs at a distance of 1.7 m, θz = ωz at 25 m. The wiggles on the binned curves are the result of taking the maximum of several peaks corresponding to all possible binnings.

(i) The detector intercepts all of the scattering, but we do not use binning. The SNR has a wide peak between the ωx = θx and ωz = θz points. With 8-ID parameters the peak is between Rdet = Rx = 1.7 m and Rdet = Rz = 25 m. This case is shown by curve `Q = 0' in Fig. 2[link].

(ii) The detector intercepts all of the scattering and we use optimal binning. As shown in curve `Q = 0 bin' in Fig. 2[link], the SNR will increase with increasing sample-to-detector distance until we have moved the detector as far as Rz. Beyond Rz, the SNR remains approximately constant (the discernible dips in the theoretical SNR from twofold to threefold binning and from threefold to fourfold binning etc.).

(iii) The detector does not intercept all of the ring of scattering of interest and we do not use binning. Here the SNR peaks at a detector distance of Rx and decreases sharply for larger sample-to-detector distances, as shown by curve `Q = 0.2 nm−1'.

(iv) The detector does not intercept all of the ring of scattering and we use optimal binning. In this case the SNR shows a broad maximum for sample-to-detector distances between Rx and Rz, and we have a wide range within which to set the sample-to-detector distance (see curve `Q = 0.2 nm−1 bin' in Fig. 2[link]).

3.2. Small slit sizes

It is clear from Fig. 1[link] that, for slit sizes smaller than about three coherence lengths, the SNR displays a significant L dependence. In fact, the maximum SNR decreases with decreasing L, while the optimum pixel acceptance ω increases with decreasing L. This is highlighted in Fig. 3[link], which plots versus L the values of ω at which the SNR is a maximum for a particular L.

[Figure 3]
Figure 3
Ratio of the angular pixel size ω and the angular source size θ at maximum signal-to-noise ratio versus slit size. The optimum values are deduced by numerically finding the maximum of curves similar to the ones in Fig. 1[link].

An analytically tractable expression for the L-dependence of the maximum SNR follows from the approximation that F(x) ≃ 1 − x2/6, which is valid for L/Ξ [\ll] 1. The location of the maximum SNR in this limit is

[\omega_{\rm{opt}}=\left(12\ln2\right)^{1/2}\lambda/\left(5^{1/2}\,\pi{L}\right)\simeq0.41\lambda/L,\eqno(18)]

which varies inversely with the slit width and, interestingly, is independent of the angular source size. Careful inspection of the curve in Fig. 3[link] for L/ξ [\ll] 1 reveals an ωopt ≃ 0.78λ/L dependence. Thus, the first-order analytical approximation is not precise enough to find the exact maximum; the scaling is preserved though. This limit is especially important in the case of X-ray free-electron lasers, where the fully coherent beam corresponds to θ = 0 and the optimum pixel size will always be determined by the beam size (speckle size).

4. Determining the best slit size

As already noted, for large slits the SNR is independent of slit size, while for small slits the SNR decreases linearly with the slit size. In this section we will explore the dependence of the SNR on the slit size for constant pixel size (ω) and constant source size (θ). Since we know that the maximum SNR is achieved in the large-slit limit [equation (17)[link]] for any ω, it is instructive to plot the ratio

[{{r_{\rm{sn}}(L)}\over{r_{\rm{sn}}(\infty)}}=F\left({{L}\over{\Xi}}\right){{1}\over{\pi^{1/2}}}\,{{L}\over{\Xi}}.\eqno(19)]

This function depends only on L/Ξ. It is apparent from Fig. 4[link] that the SNR increases monotonically with L, reaching 0.5 for L/Ξ = 1, 0.8 for L/Ξ = 3, and 0.9 for L/Ξ = 6. This confirms that, while widening the slits is generally beneficial, opening them too much past a few times the effective coherence length will not significantly increase the SNR. It is not surprising either that, for L/Ξ [\ll] 1, rsn is proportional to L, since in this case closing the slits does not increase the contrast but merely reduces the intensity.

[Figure 4]
Figure 4
Relative SNR versus slit size while keeping θ and ω fixed

If L/Ξ is very large, the optical contrast is correspondingly very small, and it may be preferable to instead maintain a certain level of contrast in order to insure against any systematic errors that might create small baseline shifts. Some examples are given in Table 1[link].

Table 1
SNR and contrast for several slit settings at 8-ID

Slit width Slit height Contrast SNR relative to maximum
2Ξx 2Ξz 0.40 0.52
4Ξx 4Ξz 0.14 0.74
6Ξx 6Ξz 0.07 0.82
11Ξx 3Ξz 0.07 0.77

As a compromise between the best SNR and non-zero contrast, it is sensible to set the slit width to be several times (three to six times) wider than the coherence length, which will yield contrast values between 0.07 and 0.14. It is also clear from the last two lines that, while setting the slit width symmetrically produces the best SNR, asymmetric set-ups yield a SNR that is only slightly lower.

5. Optimizing energy

Since the efficiency of direct-detection CCD detectors is invariably X-ray energy dependent, and since at 8-ID and other beamlines one generally has considerable freedom in choosing at what energy to carry out an XPCS experiment, it is valuable to investigate the effects of varying the X-ray energy on the SNR. We first consider how to optimize the energy in the limit of large slits. We furthermore suppose that when we tune the energy we keep fixed the range of wavevectors contributing to the particular Q of interest. It follows that ΩxΩz must be varied to maintain this range of wavevectors, while ω and θ are fixed (and optimally equal to each other). In addition, we suppose that we adjust the sample thickness (W) in order to keep the sample attenuation constant, i.e. we set W = Λ(E). Equation (17)[link] then informs us that rsnxrsnz varies as E−2 because of the factor λ2. Other energy-dependent factors in (15)[link] are (ΩxΩz)1/2, which is proportional to E−1, and Wexp(−W/Λ), which is material-dependent but is often proportional to E3. Thus, to leading order, these factors cancel each other out. The remaining factors are the source brilliance and the detector efficiency. In addition, at 8-ID we should recognize that there is an energy-dependent attenuation owing to Be windows in the beamline. For the APS, the brilliance varies weakly with energy with a peak at about 7 keV. Reducing the X-ray energy can significantly increase the CCD detector's quantum efficiency (at least within the range 4–10 keV) and can reduce the spreading of photons (Hashimotodani, 1998[Hashimotodani, K. (1998). Rev. Sci. Instrum. 69, 3746-3750.]). In its simplest form the CCD detector efficiency varies as η = 1 − exp(−ddep/ΛSi) where ddep is the depletion layer thickness, and the silicon attenuation length ΛSiE3, strongly favoring low energies. Both of these effects tend to improve the SNR at lower energies. By contrast, the attenuation of Be windows increases sharply for energies below about 5 keV. Thus, for 8-ID or other beamlines in which there are Be windows, we expect that the optimum X-ray energy for XPCS is about 6 keV (see Fig. 5[link]). Since Be windows are far from perfect optics, it would be desirable in any case to find an improved X-ray window material. With a different window material the calculus leading to the prediction that the SNR will be optimal at 6 keV will be changed somewhat. In the case of a sample with isotropic scattering, so that at high wavenumbers the scattering ring does not completely fit on the CCD, an extra factor in the SNR proportional to E1/2 becomes necessary. However, this moves the optimum energy only slightly. While in Fig. 5[link] the longitudinal and azimuthal correction was calculated for all wavenumbers, neither of the corrections have an effect on the solid line, and the dashed line is only affected by the azimuthal correction. For anisotropic systems at wavenumbers where the longitudinal correction is negligible, the solid line should be used.

[Figure 5]
Figure 5
Energy dependence of the SNR for parameters relevant to 8-ID. The brilliance of the APS undulator A, the attenuation of four 0.25 mm Be windows, and the measured efficiency and detector size of the SMD 1M60 detector (Falus et al., 2004[Falus, P., Borthwick, M. A. & Mochrie, S. G. J. (2004). Rev. Sci. Instrum. 75, 4383-4400.]) are taken into account. The SNR is calculated taking into account the effect of longitudinal coherence and the correction for azimuthal averaging. The solid line is unaffected both by the longitudinal and azimuthal correction.

Our conclusions are different from those of Thurn-Albrecht et al. (2003[Thurn-Albrecht, T., Grübel, G., Müller-Buschbaum, W. S. P. & Patkowski, A. (2003). Phys. Rev. E, 68, 031407.]), who predict that the SNR in an XPCS experiment has an overall E2 dependence. The chief points of difference between our analysis and that of Thurn-Albrecht et al. involve the factors dealing with the efficiency of a CCD camera, the angular factor of (ΩxΩz)1/2, and the assumption, made by those authors, that the total measurement time is limited by the rate of sample damage. The first two differences result from the choice of detector: an area detector in the present case and a point detector by Thurn-Albrecht et al. In this context we note that a current state-of-the-art area detector yields a much superior SNR than a point detector for correlation times greater than 2 ms (Falus et al., 2004[Falus, P., Borthwick, M. A. & Mochrie, S. G. J. (2004). Rev. Sci. Instrum. 75, 4383-4400.]). With regard to the sample damage issue, we assume that the time required to damage the sample is at least a few times longer than the relevant correlation times. Consequently, using an area detector, the sample can be periodically moved to illuminate a fresh spot, so that there is no limitation on overall measurement time (Lumma, Lurio, Mochrie & Sutton, 2000[Lumma, D., Lurio, L. B., Mochrie, S. G. J. & Sutton, M. (2000). Rev. Sci. Instrum. 71, 3274.]).

6. Longitudinal coherence

The preceding discussion of energy optimization has been made under the assumption that the longitudinal coherence length (Λ) is large in comparison with the maximum optical pathlength difference (δ) within the sample. A detailed calculation of the effect of longitudinal coherence on the contrast requires an integral of the relevant sample volume and energy distributions (Sandy et al., 1999[Sandy, A. R., Lurio, L. B., Mochrie, S. G. J., Malik, A., Stephenson, G. B., Pelletier, J.-F. & Sutton, M. (1999). J. Synchrotron Rad. 6, 1174.]) and is beyond the scope of this paper. In brief, however, with the approximation that the X-ray power spectrum is a Lorentzian, and that the dimensions of the sample transverse to the beam are sufficiently small, then the X-ray contrast will vary as exp(−|δ|/Λ). If the sample thickness is adjusted to match the attenuation length, then δ will depend on both the energy and wavevector of the measurement. Specifically, in the range of wavevectors where the sample thickness is much bigger than the beam size, equation (5)[link] may be approximated as

[\delta=W\lambda^2Q^2/8\pi^2.\eqno(20)]

The energy dependence of (20)[link] results from the factor of λ2 and the E3 dependence of the attenuation length. Consequently, since Λ varies as E−1, the ratio δ/Λ will vary as Q2E2. Based on this observation, we note that for a scattering wavevector of 0.5 nm−1, for example, it is only sensible to increase the sample thickness to match the attenuation length of the sample (for typical organic-based materials) out to about 11 keV. At 1 nm−1 the effect is so strong that it influences the optimum energy calculation (see Fig. 5[link]).

7. Focusing

So far we have not considered the possibility of focusing in this paper. Since the vertical coherence lengths at synchrotron sources are much greater than the horizontal coherence lengths, it is first sensible to examine the effect of vertical focusing on the coherence properties of the beam. To this end we consider a focusing optic a distance R1 from the source that focuses to an image at an additional distance R2 downstream. The image, which is of size [\sigma_z^\prime] = R2σz/R1, may be viewed as a new source, so that the coherence length at the sample position, which we take to be a further distance R3 downstream, is ζF = λR3/2π[\sigma_z^\prime] = λR3R1/2πσzR2. This may be compared with the coherence length at the sample position in the absence of the focusing optic, which is ζ = λ(R1 + R2 + R3)/2πσz, since R1 + R2 + R3 is the source-to-sample position. Although the value of R1 may be constrained by the location of the beamline's first optics enclosure, there remains considerable freedom in the choice of the ratio ζ/ζF = (R1 + R2 + R3)R2/(R3R1). At 8-ID-I, in particular, it might be convenient to pick ζ/ζF = 5, so that the optimum sample-to-detector distance would become 5 m, which would then lie within the 8-ID-I hutch. This could be arranged, for example, with R1 = 35 m, R2 = 22 m and R3 = 8 m, or with R1 = 55 m, R2 = 8 m and R3 = 2 m. A similar focusing at NSLS-II would reduce the optimum sample-to-detector distance there from 50 m to a more manageable 10 m.

8. Conclusion

In this paper we have explored how to optimize the signal-to-noise ratio in an X-ray photon correlation spectroscopy experiment by revisiting and analyzing the calculation of the optical contrast at a synchrotron source. We demonstrated that for fixed slit size the signal-to-noise ratio for data at a particular wavevector exhibits a maximum versus the angular pixel size. For slits wider than the transverse coherence length, we showed that the optimum signal-to-noise ratio is achieved when the angular source size and the angular pixel size are equal to each other. We quantified the well known fact that opening up the collimating slits increases the signal-to-noise ratio but that, once the slit size is larger than several times the effective coherence length, the signal-to-noise ratio does not improve significantly by further increasing the slit size. Finally, we showed that the signal-to-noise ratio depends on the energy mainly through the detector efficiency and window attenuation. Thus, it is sensible to select the X-ray energy accordingly to take advantage of higher detector efficiencies.

In summary, the best signal-to-noise ratio can be achieved by moving the detector sufficiently far back from the sample so that the vertical angular source size is the same as the angular detector pixel size, or, at 8-ID, as far back as possible within the hutch. Then, since the horizontal source size is usually much larger than the vertical source size, we have shown that one should bin together the counts in pixels in the horizontal direction so as to match the effective horizontal pixel size with the horizontal source size.

Acknowledgements

We thank Matthew Borthwick, Dirk Lumma, Mark Sutton and Brian Stephenson for valuable discussions. This work was supported by the NSF via Grant DMR 0453856.

References

First citationAbernathy, D. L., Grübel, G., Brauer, S., McNulty, I., Stephenson, G. B., Mochrie, S. G. J., Sandy, A. R., Mulders, N. & Sutton, M. (1998). J. Synchrotron Rad. 5, 37–47.  Web of Science CrossRef CAS IUCr Journals Google Scholar
First citationBrauer, S., Stephenson, G. B., Sutton, M., Brüning, R., Dufresne, E., Mochrie, S. G. J., Grübel, G., Als-Nielsen, J. & Abernathy, D. L. (1995). Phys. Rev. Lett. 74, 2010.  CrossRef PubMed Google Scholar
First citationDierker, S. B., Pindak, R., Fleming, R. M., Robinson, I. K. & Berman, L. E. (1995). Phys. Rev. Lett. 75, 449–452.  CrossRef PubMed CAS Web of Science Google Scholar
First citationDufresne, E. M., Nurushev, T., Clarke, R. & Dierker, S. B. (2002). Phys. Rev. E, 65, 061507.  CrossRef Google Scholar
First citationFalus, P., Borthwick, M. A. & Mochrie, S. G. J. (2004). Rev. Sci. Instrum. 75, 4383–4400.  Web of Science CrossRef CAS Google Scholar
First citationFalus, P., Borthwick, M. A. & Mochrie, S. G. J. (2005). Phys. Rev. Lett. 94, 016105.  Web of Science CrossRef PubMed Google Scholar
First citationHashimotodani, K. (1998). Rev. Sci. Instrum. 69, 3746–3750.  Web of Science CrossRef CAS Google Scholar
First citationJakeman, E. (1973). Photon Correlation and Light Beating Spectroscopy, edited by H. Z. Cummins and E. R. Pike, pp. 75–149. New York: Plenum.  Google Scholar
First citationKim, H. J., Ruehm, A., Lurio, L. B., Basu, J., Lal, J., Lumma, D., Mochrie, S. G. J. & Sinha, S. K. (2003). Phys. Rev. Lett. 90, 068302.  CrossRef PubMed Google Scholar
First citationLivet, F., Bley, F., Caudron, R., Geissler, E., Abernathy, D., Detlefs, C., Grubel, G. & Sutton, M. (2001). Phys. Rev. E, 63, 036108.  Web of Science CrossRef Google Scholar
First citationLumma, D., Borthwick, M. A., Falus, P., Lurio, L. B. & Mochrie, S. G. J. (2001). Phys. Rev. Lett. 86, 2042.  CrossRef PubMed Google Scholar
First citationLumma, D., Lurio, L. B., Mochrie, S. G. J. & Sutton, M. (2000). Rev. Sci. Instrum. 71, 3274.  CrossRef Google Scholar
First citationLumma, D., Lurio, L. B., Sandy, A. R., Borthwick, M., Falus, P., Mochrie, S. G. J., Sutton, M. & Regan, L. (2000). Phys. Rev. E, 62, 8258.  CrossRef Google Scholar
First citationLurio, L. B., Lumma, D., Falus, P., Borthwick, M. A., Mochrie, S. G. J., Pelletier, J.-F., Sutton, M., Malik, A. & Stephenson, G. B. (2000). Phys. Rev. Lett. 84, 785.  CrossRef PubMed Google Scholar
First citationMadsen, A., Als-Nielsen, J. & Grübel, G. (2003). Phys. Rev. Lett. 90, 085701.  Web of Science CrossRef PubMed Google Scholar
First citationMochrie, 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.  CrossRef Web of Science Google Scholar
First citationPrice, A. C., Sorensen, L. B., Kevan, S. D., Toner, J., Poniewierski, A. & Holyst, R. (1999). Phys. Rev. Lett. 82, 755.  CrossRef Google Scholar
First citationSandy, A. R., Lurio, L. B., Mochrie, S. G. J., Malik, A., Stephenson, G. B., Pelletier, J.-F. & Sutton, M. (1999). J. Synchrotron Rad. 6, 1174.  Web of Science CrossRef IUCr Journals Google Scholar
First citationSeydel, T., Madsen, A., Tolan, M., Grubel, G. & Press, W. (2001). Phys. Rev. B, 63, 073409.  Web of Science CrossRef Google Scholar
First citationSikharulidze, I., Dolbnya, I. P., Fera, A., Madsen, A., Ostrovskii, B. I. & de Jeu, W. H. (2002). Phys. Rev. Lett. 88, 115503.  Web of Science CrossRef PubMed Google Scholar
First citationSutton, M., Mochrie, S. G. J., Greytak, T., Nagler, S. E., Berman, L. E., Held, G. A. & Stephenson, G. B. (1991). Nature (London), 352, 608.  CrossRef Google Scholar
First citationThurn-Albrecht, T., Grübel, G., Müller-Buschbaum, W. S. P. & Patkowski, A. (2003). Phys. Rev. E, 68, 031407.  Google Scholar
First citationThurn-Albrecht, T., Meier, G., Müller-Buschbaum, P., Patkowski, A., Steffen, W., Grübel, G., Abernathy, D. L., Diat, O., Winter, M., Koch, M. G. & Reetz, M. T. (1999). Phys. Rev. E, 59, 642–649.  Web of Science CrossRef CAS Google Scholar
First citationThurn-Albrecht, T., Steffen, W., Patkowski, A., Meier, G., Fisher, E. W., Grübel, G. & Abernathy, D. L. (1996). Phys. Rev. Lett. 77, 5437–5440.  PubMed CAS Google Scholar
First citationTsui, O. K. C. & Mochrie, S. G. J. (1998). Phys. Rev. E, 57, 2030.  CrossRef Google Scholar

© International Union of Crystallography. Prior permission is not required to reproduce short quotations, tables and figures from this article, provided the original authors and source are cited. For more information, click here.

Journal logoJOURNAL OF
SYNCHROTRON
RADIATION
ISSN: 1600-5775
Follow J. Synchrotron Rad.
Sign up for e-alerts
Follow J. Synchrotron Rad. on Twitter
Follow us on facebook
Sign up for RSS feeds