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

Journal logoJOURNAL OF
SYNCHROTRON
RADIATION
ISSN: 1600-5775

Enhancement of the spin-dependent effect of π/2-angle Compton scattering using elliptically polarized synchrotron radiation

aDepartment of Material Science, Himeji Institute of Technology, 3-2-1 Kouto, Kamigori,Ako-gun, Hyogo 678-1297, Japan
*Correspondence e-mail: n_sakai@sci.himeji-tech.ac.jp

(Received 8 November 1999; accepted 5 May 2000)

Noticeable enhancement of the spin-dependent Compton scattering intensity in right-angle scattering has been achieved by using high-energy X-rays having elliptical polarization. A promising momentum resolution of better than 0.4 atomic units for Compton-profile measurements is firmly predicted by means of an ordinary Ge solid-state detector when the scattering-angle divergence at 90° is crucially restricted. It is pointed out that the angle between a scattering vector and the direction of sample magnetization can be chosen as 90° without seriously weakening the spin-dependent effect. Comparison is made between 274 keV and 122 keV experiments.

1. Introduction

The presence of a magnetic interaction between the electromagnetic field and the electron orbital currents and intrinsic spins makes synchrotron radiation one of the most useful tools for the study of magnetism. Many theoretical and experimental studies have been carried out and valuable information has been obtained. The magnetic interaction, however, is weak in comparison with the charge interaction, except for some resonance magnetic absorption and scattering. In the case of Compton scattering on ordinary ferromagnetic samples, the magnetic effect contributes to, at most, a few percent of the total scattering intensity for circularly polarized photons. This limitation has been improved to some extent by adopting right-angle scattering geometry (RSG), in which the charge scattering is greatly suppressed, as described below. RSG has been commonly adopted in magnetic Bragg-diffraction experiments for the same purpose (Brunel et al., 1983[Brunel, M., Patrat, G., de Bergevin, F., Rousseaux, F. & Lemonnier, M. (1983). Acta Cryst. A39, 84-88.]; Ito & Hirano, 1997[Ito, M. & Hirano, H. (1997). J. Phys. Condens. Matter, 9, L613-L618.]).

An energy spectrum of Compton-scattered photons in a fixed direction (Compton profile) contains valuable information on the momentum distribution of electrons in matter. Hitherto, Compton profiles have been measured by using the backscattering geometry (BSG). This is based on the fact that, firstly, the Doppler effect on the scattered photons is large and, secondly, the scattering-angle divergence does not affect the momentum resolution of the Compton profile very much, and thus BSG gives a good counting efficiency. In addition, the spin-dependent Compton-scattering cross section is large for BSG. After utilizing high-energy and high-flux synchrotron radiation emitted from a low-emittance storage ring, like at the ESRF and SPring-8, it has been recognized that these advantages of BSG (McCarthy et al., 1997[McCarthy, J. E., Cooper, M. J., Lawson, P. K., Timms, D. N., Manninen, S. O., Hamalainen, K. & Suortii, P. (1997). J. Synchrotron Rad. 4, 102-109.]) are also fulfilled by RSG. The present paper reports a noticeable enhancement of the magnetic effect on the spin-dependent Compton profile of Fe of up to 3.4% with a reasonable momentum resolution of 0.5 atomic units (a.u.) by using an ordinary Ge solid-state detector (SSD) as an energy analyzer.

2. Theoretical analysis

The differential scattering cross section for Compton scattering for electrons at rest can be written as (Lipps & Tolhoek, 1954[Lipps, F. W. & Tolhoek, H. A. (1954). Physica, 20, 85-98; 395-404.])

[\eqalign{ & {{{{\rm d}\sigma } \over {{\rm d}\Omega }}} = r\,_0^2 \left({{{k^2 } \over {k_0 ^2 }}} \right)\left [{\Phi _0 + P_L \Phi _1 + P_C \Phi _2 } \right], \cr & \Phi _0 (\theta, E_0, E) = (1/8)\left [{(1 + \cos ^2 \theta) + {{E_0 - E} \over {mc^2 }}(1 - \cos \theta)} \right], \cr & \Phi _1 (\theta) = (1/8)\sin ^2 \theta, \cr & \Phi _2 (\theta, {\bf k}_0, {\bf k}, \boldsigma) = (1/8)(1 - \cos \theta)\boldsigma \cdot {\hbar \over {mc}}({\bf k}_0 \cos \theta + {\bf k}),}\eqno(1)]

where r0 denotes the electromagnetic radius of the electron, k0 and E0 (k and E) the wave number and energy of the incident (scattered) photons, respectively, [\boldsigma] the unit vector parallel to the direction of the electron spin, and θ the scattering angle.

Here we define a figure of merit f(θ1, θ2, θ), which is adequate to optimize the relative intensity of the spin-dependent Compton scattering with respect to the total scattering intensity. The following relation expresses the signal-to-noise ratio as a function of the Stokes parameters PL and PC, the flux density Nph of emitted photons, and the vertical width of the front-end slit represented by angles θ1 and θ2,

[\eqalignno{f(\theta _1, \theta _2, \theta)={}&\Phi _2 (\theta, {\bf k}_0, {\bf k},\boldsigma)\!\!\textstyle\int\limits_{ - \infty }^\infty \!\!{\rm d}\theta _x \textstyle\int\limits_{\theta _2 }^{\theta _1 }\!N_{\rm ph} (\theta _x, \theta _y) P_C (\theta _x, \theta _y)\,{\rm d}\theta _y \cr&\big/ \Big\{\textstyle\int\limits_{ - \infty }^\infty \!\!{\rm d}\theta _x \textstyle\int\limits_{\theta _2 }^{\theta _1 } \!N_{\rm ph} (\theta _x, \theta _y) \cr&{\times}\left[\Phi _0 (\theta, E_0, E) + P_L (\theta _x, \theta _y)\Phi _1 (\theta) \right]\,{\rm d}\theta _y \Big\}^{1/2}\!.&(2)}]

Here θx and θy are horizontal and vertical angles, respectively, of the photons as shown in Fig. 1[link]. The denominator expresses the noise, which is the statistical error of the total number of scattered photons. When the photon flux is low, as in the case for 274 keV X-rays, we are forced to use a wide slit and optimize the slit width with parameters θ1 = θH and θ2 = −θH. On the other hand, when the incident flux density is high, as in the case for 122 keV X-rays, we are able to make the front-end slit narrow (θ1θ2) to obtain the best condition. For photons of a few hundred keV, a scattering angle of 90° promises to reduce the denominator of (2)[link], and the f factor for a low photon flux is proportional to

[\eqalignno{h(\theta _H, - \theta _H, 90) ={}& \textstyle\int\limits_{ - \infty }^\infty {\rm }d\theta _x \textstyle\int\limits_0^{\theta _H } N_{\rm ph} (\theta _x, \theta _y) P\!_C (\theta _x, \theta _y)\,\,{\rm d}\theta _y\cr& \big/ \big\{\textstyle\int\limits_{ - \infty }^\infty {\rm d}\theta _x \textstyle\int\limits_0^{\theta _H } N_{\rm ph} (\theta _x, \theta _y) \left[1 + P\!_L (\theta _x, \theta _y)\right.\cr&\left. +\,\,(E_0 - E)/mc^2\right]\,\,{\rm d}\theta _y\big\}^{1/2}.&(3)}]

This is a complex function of θH, because PC, PL and Nph have a distribution in space as shown in Figs. 2[link](a) and 2[link](b). The parameter Kx is proportional to the horizontal magnetic field intensity of a wiggler. This field is necessary for producing circularly polarized photons, and induces a twin-peaked flux-density distribution of the X-rays. PC is high along the z-axis, while PL is high at the peaks. Fig. 3[link] shows the evaluated h(θH, −θH, 90) as a function of θH and Kx. Numerical values of Nph, PC and PL were referred to a computer simulation performed using the program SPring-8 spectra (Tanaka & Kitamura, 1998[Tanaka, T. & Kitamura, H. (1998). http://radiant.harima.riken.go.jp/spectra/index_e.html ]). In the case of the 274 keV photon flux, as shown in Fig. 3[link], it was found that, for any Kx parameter greater than 0.3, the maximum of h is approximately the same, and the values Kx = 1.0 and θH = 0.057 mrad were adopted. We can understand from Fig. 2[link](b) that the selected value of θH covers the flux density between the two peaks.

[Figure 1]
Figure 1
Definition of θx and θy of emitted X-rays.
[Figure 2]
Figure 2
Distributions of (a) flux density Nph and (b) Stokes parameters for linear polarization PL and circular polarization PC of X-rays emitted from the elliptical multipole wiggler. Ky is fixed at 9.894.
[Figure 3]
Figure 3
Calculated h(θH, 0, 90) for 274 keV X-rays as a function of Kx and θH.

The following ratio g(θ1,θ2, θ) between the spin-dependent intensity and the total intensity is required later for comparison with experimental data,

[\eqalignno{g(\theta _1, \theta _2, \theta) ={}& \Phi _2 (\theta, \boldsigma, {\bf k}_0, {\bf k})\!\textstyle\int\limits_{ - \infty }^\infty \!\!{\rm d}\theta _x \textstyle\int\limits_{\theta _2 }^{\theta _1 } \!N_{\rm ph} (\theta _x, \theta _y) P\!_C (\theta _x, \theta _y)\,{\rm d}\theta _y\cr& \big/ \textstyle\int\limits_{ - \infty }^\infty {\rm d}\theta _x \textstyle\int\limits_{\theta _2 }^{\theta _1 } N_{\rm ph} (\theta _x, \theta _y)\cr& {\times}\left[\Phi _0 (\theta, E_0, E) + P\!_L (\theta _x, \theta _y)\Phi _1 (\theta) \right]\,\,{\rm d}\theta _y\cr&{\times}(m_{\rm spin}/n_{\rm tot}).&(4)}]

Here ntot denotes the total electron number per atom of the sample, and mspin denotes the spin component of magnetization of the sample in units of the Bohr magneton, μB; the cross section of the magnetic Compton scattering does not depend on the orbital angular momentum (Lovesey, 1993[Lovesey, S. W. (1993). Rep. Prog. Phys. 56, 257-326.]; Sakai, 1994[Sakai, N. (1994). J. Phys. Soc. Jpn, 63, 4655-4656.]).

In addition to enhancing the spin-dependent Compton scattering intensity, a good momentum resolution for Compton-profile measurements is highly desired so that details of the momentum distribution of the electrons can be clarified. The momentum component pz along the scattering vector can be experimentally evaluated by the equation

[p_z = {{E - E_0 + E_0 E(1 - \cos \theta)/mc^2}\over{\left(E\,_0^2 + E\,^2 - 2E_0 E\cos \theta\right)^{1/2}}}mc.\eqno(5)]

Here the constant mc is 137.036 a.u. According to (5)[link], a finite momentum resolution Δpz is caused by the uncertainty of E0, that of E arising from the energy resolution of a detector, and the divergence of θ,

[\Delta p_z = \left[\left({{{\partial p_z } \over {\partial E_0 }}\Delta E_0 } \right)^2 + \left({{{\partial p_z } \over {\partial E}}\Delta E} \right)^2 + \left({{{\partial p_z } \over {\partial \theta }}\Delta \theta } \right)^2\right]^{1/2},\eqno(6)]

with

[{{\partial p_z } \over {\partial E_0 }} = {{ - mcE\left [{E_0 + E - {{E(E - E_0 \cos \theta)/mc^2 }}} \right](1 - \cos \theta)} \over {(E\,_0^2 + E\,^2 - 2E_0 E\cos \theta)^{3/2} }},\eqno(7)]

[{{\partial p_z } \over {\partial E}} = {{mcE_0 \left [E_0 + E + E_0 (E_0 - E\cos \theta)/mc^2 \right](1 - \cos \theta)} \over {(E\,_0^2 + E\,^2 - 2E_0 E\cos \theta)^{3/2} }},\eqno(8)]

[\eqalignno{&{{\partial p_z} \over{\partial \theta }} =\cr&\ {{mcEE_0 \left\{E_0 \!- E + \left[E\,_0^2 + E\,^2 \!- E_0 E(1 + \cos \theta)\right]/mc^2 \right\}\sin \theta } \over {(E\,_0^2 + E\,^2 \!- 2E_0 E\cos \theta)^{3/2} }}.\cr&&(9)}]

The first and the third terms in (6)[link] induce Compton-peak shifts of Compton profiles, and deteriorate a resultant momentum resolution. The first and second terms in (6)[link] introduce the ambiguity of the scattering vector. Equations (7)[link] and (8)[link] indicate that the ambiguity of E0 and E considerably worsens the momentum resolution when scattering angles are large. Equation (9)[link] indicates that the influence of Δθ is most significant at θ = π/2. Then, Δθ should be as small as possible to achieve a good momentum resolution for experiments with RSG.

It is important to compare the amount of the peak shift ΔEθ caused by Δθ with the amount of the Doppler shift ΔED. The former can be evaluated by using (8)[link] and (9)[link] to be

[\eqalignno{ \Delta E_\theta &= {{\partial E} \over {\partial \theta }}\Delta \theta = \left[\left({{\partial p_z } \over {\partial \theta }} \right) \bigg/ \left({{\partial p_z } \over {\partial E}} \right)\right]\Delta \theta \cr & = {{E\left [\!{1 \!- {E \over {E_0 }} + {{E_0 } \over {mc^2 }}\!\left\{ {1 - {E \over {E_0 }}(1 + \cos \theta) + \left({{E \over {E_0 }}} \right)^2 }\! \right\}} \right]\sin \theta \!} \over {1 + {E \over {E_0 }} + {{E_0 } \over {mc^2 }}\left({1 - {E \over {E_0 }}\cos \theta } \right)(1 - \cos \theta)}}\Delta \theta,\cr& &(10)}]

and the latter is

[\Delta E_D = {{E_0 \left[ 1 + \left(E / E_0 \right)^2 - 2(E / E_0)\cos \theta \right]^{1/2}} \over {1 + (E_0 / mc^2 )(1 - \cos \theta)}}{{\Delta p_z } \over {mc}}.\eqno(11)]

Fig. 4[link] shows the ratios, as a function of θ, between ΔED for Δpz = 1 a.u. and an instrumental resolution ΔE that consists of an energy resolution of a detector and ΔEθ. When Δθ = 0, it can be seen that the ratio ΔED/ΔE is more enhanced at θ = π/2 by higher-energy photons. This advantage is in agreement with the enhancement of the spin-dependent Compton scattering, and beneficial to measuring magnetic Compton profiles at θ = π/2; it should be noted that this merit diminishes when Δθ becomes large, as demonstrated in Fig. 4[link].

[Figure 4]
Figure 4
Ratios in arbitrary units between the Doppler broadening and an instrumental energy uncertainty ΔE as a function of the scattering angle θ and the scattering-angle divergence Δθ.

3. Experimental details

Experiments using 274 keV and 122 keV photons have been carried out on station A and station B, respectively, at beamline BL08W of the SPring-8 facility, Japan. Fig. 5[link] shows the experimental set-up for 274 keV photons. The sample was a 5 mm-thick polycrystalline Fe plate incorporated into a C-type small electromagnet. According to (2)[link], f for θ = 90° is optimized when the magnetization, i.e. the spin direction [\boldsigma], is parallel to the scattered photon direction k. This arrangement, however, causes the serious absorption of scattered photons by the Fe sample. The figure of merit f was then re-estimated by this sample absorption, and 31.3° was evaluated as the best angle between [\boldsigma] and k. The Compton-scattered photons were detected by a Ge-SSD1 (Eurisys, EGPC 30-185-R: diameter 51.6 mm, thickness 66.6 mm) and a Ge-SSD2 (Canberra, GL0115: diameter 11.3 mm, thickness 15 mm) detector. The peak energy of the Compton profile was 178.4 keV. Ge-SSD1 has a good detection efficiency of about 70% at this energy, but its low-energy resolution gives a corresponding momentum resolution of 2.4 a.u. A relatively large angle divergence of 56.6 mrad was adopted in compliance with the detector resolution, and the Δθ-dependent momentum resolution was 2.27 a.u. On the other hand, Ge-SSD2 has a detection efficiency of 49.7% and its energy resolution gives a momentum resolution of 0.33 a.u. An angular divergence Δθ = 9.5 mrad was then adopted to realize a tolerable overall momentum resolution. The counting rate, however, was very low due to the low flux of 274 keV photons and the small scattering-angle allowance. Although the present result for Ge-SSD2 is not satisfactory in terms of the counting efficiency, the result suggests a further improvement in Δpz of better than 0.4 a.u. if a smaller Δθ could be adopted. The observed energy spectra for Iup + Idown and IupIdown are shown in Fig. 6[link], where Iup and Idown correspond to the intensities of Compton profiles with a positive and negative magnetic field, respectively. The spin-dependent effect can be evaluated by the ratio g = (IupIdown)/(Iup + Idown − BG), where BG denotes the background intensity. This g ratio is equivalent to (4)[link].

[Figure 5]
Figure 5
Experimental set-up for 274 keV X-rays. The linear polarization vector is in the scattering plane, giving PL ≃ −1.
[Figure 6]
Figure 6
Experimental total and magnetic Compton profiles measured using Ge-SSD1 and Ge-SSD2. The incident energy is 274 keV.

Similar measurements have been made using 122 keV photons and Ge-SSD2. The high photon flux density of 122 keV photons allowed a narrow width of the front-end slits to be used. The Compton-peak energy was 97.6 keV. The detector has 100% detection efficiency at this energy, and a reasonable energy resolution of the detector gives a momentum resolution of 0.42 a.u. at the Compton peak. The sample was a 2 mm-thick polycrystalline Fe plate and the magnetization direction was optimized at 31.1°. The beam spot at the sample position measured by a Polaroid film (Polapan 57) was 2 mm wide and 7 mm high. The distance between the sample and the detector was 175 cm, and a narrow slit of width 1 mm was placed in front of the detector. The resultant Δθ was 2.1 mrad. Since the flux of 122 keV photons is about 103 times higher than that of 274 keV photons, the optimized condition for θH evaluated by (3)[link] induced a tremendous amount of scattered photons, causing the Ge detector to overload: the dead time of the electronics for the energy analysis increased almost up to 100%. Then the front-end slit was crucially collimated to 3.06 µrad at θ1 = 0.041 mrad, which selects the best polarization parameters, as shown in Fig. 7[link], and the monochromated photon flux was adjusted to an adequate intensity for the detector system. The factor mmag in equation (4)[link] was experimentally evaluated from the magnetization of the sample: mmag = 2.07 for a 2 mm-thick Fe sample and mmag = 1.87 for a 5 mm-thick sample, each of which was measured by winding a 25 turn pick-up coil around the sample and connecting it to a flux meter; here the orbital magnetic moment of Fe is neglected. The reduction of mmag for the thicker sample indicates incomplete magnetization of the sample.

[Figure 7]
Figure 7
h(θH + Δ, θH, 90) and flux density distributions for 122 keV X-rays and Kx = 1.0 with a slit width Δ = 3.06 µrad.

4. Results and conclusion

The experimental results are summarized in Table 1[link]. The experimental g ratios are compared with the theoretical ones. Good agreement is found within the experimental errors, which are mainly due to the ambiguity of the background subtraction. As may be seen from Table 1[link] for Δpz(E0), which directly reflects the accuracy of E0, a single bent monochromator for 274 keV X-rays gives a good energy resolution ΔE0/E0 better than 10−3, while that for 122 keV X-rays was not so good. This was simply due to a tentative adjustment of a doubly bent monochromator for 122 keV X-rays. The highest g ratio, 3.4%, is observed for 122 keV photons. This ratio is better than the best value for BSG of 2.4% at BL08W. From the viewpoint of the statistical accuracy, this enhancement is equivalent to twice the reduction of the accumulation time. If the flux density of the 274 keV photons increased comparably with that of the 122 keV photons, the front-end slits could be narrowed, and the g ratio for the 274 keV photons would become 5%, because the ratio Φ2/(Φ0 + PLΦ1) is increased with increasing photon energy.

Table 1
Experimental results and theoretical estimation

The g ratio of experiment 2 is eliminated because of its low statistical accuracy.

Item Experiment 1 Experiment 2 Experiment 3
Accumulation time for Iup + Idwn (h) 44.3 44.3 29.4
Count rate under Compton line (s−1) 104 1.37 858
Detector Ge-SSD1 Ge-SSD2 Ge-SSD2
E0 (keV) 274 274 122
θ (°) 90 90 90
g (experiment) (%) 3.2 ± 0.5 3.38 ± 0.35
g (theory) (%) 3.52 3.52 3.49
pz/∂θ (a.u. rad−1) 40.08 40.08 20.25
Δθ (mrad) 56.6 9.5 2.1
Δpz(Δθ) (a.u.) 2.27 0.382 0.0425
Δpz(ΔE) (a.u.) 0.795 0.327 0.423
Δpz(ΔE0) (a.u.) <0.01 <0.01 0.458
Δpz (a.u.) 2.41 0.503 0.625

A shallow dip at the top of the magnetic Compton profile in Fig. 8[link] measured by using 122 keV X-rays together with Ge-SSD2 clearly demonstrates that the RSG method also guarantees a good momentum resolution comparable with the BSG method (Sakai, 1996[Sakai, N. (1996). J. Appl. Cryst. 29, 81-99.]; McCarthy et al., 1997[McCarthy, J. E., Cooper, M. J., Lawson, P. K., Timms, D. N., Manninen, S. O., Hamalainen, K. & Suortii, P. (1997). J. Synchrotron Rad. 4, 102-109.]). If the flux density of the 274 keV photons could be increased, it would be possible to improve Δpz(Δθ) by using a narrower slit, and an overall momentum resolution better than 0.4 a.u. could be achieved.

[Figure 8]
Figure 8
Experimental total and magnetic Compton profiles measured using Ge-SSD2. The incident energy is 122 keV.

According to Φ2 in (1)[link], the spin direction is desirable to be parallel to the direction of the vector k0cos θ + k. In the case of BSG, the direction of this vector can be approximated by the scattering vector K = k0k. Then the z-axis (zK) and an applied magnetic field are required to be approximately parallel. Since the direction of K is roughly approximated by k0, the shape of a sample for magnetic Compton-profile measurement is usually prepared to be long in the direction along the incident photon in consideration of the demagnetization effect. In addition, the large magnetic anisotropy sometimes makes it difficult to magnetize a single crystalline sample in the best direction. RSG experiments allow a large angle between the z-axis and the applied magnetic field, providing a degree of freedom to the experimental set-up. Suppose we use a disc as a sample, and magnetize it (the direction of [\boldsigma]) parallel to the surface. When the incident energy is 120 keV, the angle between the incident X-rays and the z-axis is 39°. If we set the sample surface at an angle of 51° off the incident direction, the direction of the z-axis can be perpendicular to the applied magnetic field.

Acknowledgements

The authors gratefully acknowledge Dr Y. Sakurai, Dr J. Mizumaki and Dr M. Itoh of the Japan Synchrotron Radiation Research Institute for their help and advice at BL08W of SPring-8. This research was performed under the SPring-8 approval numbers 1998A0122-NM-np and 1999A0141-CD-np.

References

First citationBrunel, M., Patrat, G., de Bergevin, F., Rousseaux, F. & Lemonnier, M. (1983). Acta Cryst. A39, 84–88.  CrossRef CAS IUCr Journals Google Scholar
First citationIto, M. & Hirano, H. (1997). J. Phys. Condens. Matter, 9, L613–L618.  CrossRef CAS Web of Science Google Scholar
First citationLipps, F. W. & Tolhoek, H. A. (1954). Physica, 20, 85–98; 395–404.  Google Scholar
First citationLovesey, S. W. (1993). Rep. Prog. Phys. 56, 257–326.  CrossRef CAS Web of Science Google Scholar
First citationMcCarthy, J. E., Cooper, M. J., Lawson, P. K., Timms, D. N., Manninen, S. O., Hamalainen, K. & Suortii, P. (1997). J. Synchrotron Rad. 4, 102–109.  CrossRef CAS Web of Science IUCr Journals Google Scholar
First citationSakai, N. (1994). J. Phys. Soc. Jpn, 63, 4655–4656.  CrossRef CAS Web of Science Google Scholar
First citationSakai, N. (1996). J. Appl. Cryst. 29, 81–99.  CrossRef CAS Web of Science IUCr Journals Google Scholar
First citationTanaka, T. & Kitamura, H. (1998). http://radiant.harima.riken.go.jp/spectra/index_e.html  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