Spatial non-cyclic geometric phase in neutron interferometry.

We present a split-beam neutron interferometric experiment to test the non-cyclic geometric phase tied to the spatial evolution of the system: the subjacent two-dimensional Hilbert space is spanned by the two possible paths in the interferometer and the evolution of the state is controlled by phase shifters and absorbers. A related experiment was reported previously by some of the authors to verify the cyclic spatial geometric phase. The interpretation of this experiment, namely to ascribe a geometric phase to this particular state evolution, has met severe criticism. The extension to non-cyclic evolution manifests the correctness of the interpretation of the previous experiment by means of an explicit calculation of the non-cyclic geometric phase in terms of paths on the Bloch-sphere. The theoretical treatment comprises the cyclic geometric phase as a special case, which is confirmed by experiment.

Key words: geometric phase; neutron interferometry.

1. Introduction

Since the discovery of a geometric effect by Berry [1] in the shape of an additional phase factor after an adiabatic and cyclic transport of a quantum system, Berry's phase has been intensively investigated and generalized: the extension to degenerate subspaces by Wilckzek [2], the removal of the adiabatic constraint by Aharonov and Anandan [3] and the cyclic condition by Samuel and Bhandari [4] using the early ideas of Pancharatnam [5] and the kinematic approach to geometric phases by Mukunda and Simon [6], to name a few. In all theses contexts the geometric phase is dependent only on the geometry of the subjacent Hilbert space, but not on the particular dynamics of the system under consideration. Furthermore, Manini and Pistolesi [7] proposed an off-diagonal geometric phase to exhibit the geometry of state space in situations where the usual (diagonal) geometric phase is undefined. This has been verified experimentally by some of the authors [8].

In the course of the development of quantum mechanics it has become clear that the concept of pure states is not sufficient when taking environmental influences causing decoherence effects into account. Then one has to use the concept of mixed states. Probably the first treatise of a geometric phase for mixed states is due to Uhlmann [9] in a quantum algebraic context. Another definition of a mixed state geometric phase is due to Sjoqvist et al. [10] using an interferometric approach for its definition. For these concepts one has to keep in mind that there exist points in parameter space for which the mixed state geometric phases remain undefined provoking an extension to off-diagonal mixed state geometric phases [11].

The geometric phase is associated with an evolution of a system governed by an Hamiltonian, e. g., a neutron in a magnetic field where the geometric phase arises by the spinor evolution due to the coupling with the magnetic field. Here, we observe a geometric phase as an effect of the change in the spatial degrees of freedom in an interferometry setup. A proposal to verify the spatial geometric phase is due to Sjoqvist [12] using polarized neutrons by reversing the roles of the magnetic field and the spatial degrees of freedom. Moreover, an experiment using unpolarized neutrons has been performed by Hasegawa et al. [13, 14] to test the cyclic spatial geometric phase by inducing a relative phase shift of 2[pi] between the interfering neutron beams in a perfect silicon single-crystal interferometer. The geometric interpretation of this experiment has been dismissed by Wagh [15] demanding further investigations, namely in the non-cyclic case, which is the purpose of the current article.

2. Geometric Phases

Let us briefly review the basic concepts of geometric phases: A geometric phase is a quantity which is deeply connected to the curvature of some underlying (state- or parameter-) space. A two-dimensional plane in three-dimensional real space does not have an intrinsic curvature, but when considering a sphere embedded in euclidean real space, we have to take the curvature of this manifold into account. In geometry this curvature is reflected, for example, in the angle difference of a vector transported around a loop along geodesics, i. e., great circles: If a vector is pinned onto a sphere and then transported along a meridian to the equator, for some angle [alpha] along the equator and back to the initial point without changing its length and its direction in the tangent plane to the surface of the sphere, the vector will point in a different direction with a relative angle of [alpha] as the holonomy associated with the loop. If we do the same on a two-dimensional plane the initial and the final vector will point in the same direction.

Berry [1] was the first who addressed this issue in quantum mechanics: He considered a system initially in an eigenstate |n(R(t))>[.sub.t=0] of the governing Hamiltonian H(R(t)) dependent on the parameters R(t) changing with time t. As a demonstrative example one may consider a neutron coupling to a magnetic field H(R(t)) = -[mu] * B(R(t)) due to its magnetic dipole moment [mu] = [[mu].sub.n][sigma], where [sigma] = {[[sigma].sub.x], [[sigma].sub.y], [[sigma].sub.z]} are the Pauli matrices and [[mu].sub.n] denotes the magnetic moment of a neutron. Suppose now that the neutron is initially polarized in the direction of the magnetic field. If the direction of the magnetic field is changed adiabatically, i. e., slowly enough to avoid transitions to an orthogonal state, the system will stay in the eigenstate |n(R(t))> at all times t. Furthermore, when tracing out a loop in parameter space the final state |[PSI]([tau])> at time [tau] will be the same as the initial state up to an additional phase factor:

|[PSI]([tau])> = [e.sup.i[phi]d][e.sup.i[phi]g]|n(R([tau]))> = [e.sup.-[i/h][[integral].sub.0.sup.[tau]][E.sub.n](t)dt][e.sup.-[[??].sub.C]<n(R)|[[nabla].sub.R]|n(R)>dR]|n(R([tau]))>. (1)

The first phase value [[phi].sub.d] = [1/h][[integral].sub.0.sup.[tau]][E.sub.n](R(t))dt is dependent on the time needed to traverse the loop and on the instantaneous energy [E.sub.n](t) = <[PSI](t)|H(t)|[PSI](t)> of the system, whereas the second phase [[phi].sub.g] = i[[??].sub.C] <n(R)|[[nabla].sub.R]|n(R)> dR is dependent only on the circuit integral in parameter space revealing the geometric structure. The latter is termed Berry phase or more general geometric phase in contrast to the former dynamical phase [[phi].sub.d]. [[phi].sub.g] can be rewritten as a surface integral by use of Stoke's Theorem yielding [[phi].sub.g] = -Im[[integral].sub.F]dS[V.sub.n](R), where F is the surface enclosed by the loop in parameter space with dS denoting the area element and [V.sub.n] = [nabla] X <n|[[nabla].sub.n]> in an obvious abbreviated notation. For the neutron example--or more generally for any spin-1/2 particle--[[phi].sub.g] equals half of the solid angle enclosed by the loop as seen from the degeneracy point |R| = 0 in parameter space. This can also be related to the example from geometry above where the holonomy after the transport of the vector pinned initially to the north pole of a sphere equals the solid angle as seen from the origin of the sphere.

Several restrictions have been relaxed in course of the years, e.g., extensions to nonadiabatic [3], noncyclic and nonunitary [4], and nonpure [9, 10] geometric phases have been made. Important in our case are the generalizations to the nonadiabatic regime by Aharonov and Anandan and to noncyclic paths by Samuel and Bhandari.

In this case we have to introduce the Projective Hilbert space (Ray space) R by identifying all state vectors in Hilbert space H which differ only by an overall phase factor:

|[phi]> [approximately] |[phi]'>:|[phi]'> = [e.sup.i[alpha]]|[phi]>, [alpha][member of]R. (2)

The stress is therefore shifted from the parameter space of the Hamiltonian in case of Berry's construction to state space. We are not interested in the changes of the driving parameters (as the direction and strength of the magnetic field) but in the changes of the state itself. In Berry's considerations these two spaces are identical since the state follows the changes in parameters due to the adiabaticity condition.

In the construction by Aharonov and Anandan one considers an open path in Hilbert space which is projected to a path in Ray space by use of the equivalence relation in Eq. (2), i. e., the curve C: t[member of][0, [tau]] [right arrow] |[phi](t)> [member of] H is projected to [~.C]: t [member of] [0, [tau]] [right arrow] [pi](|[phi](t)>) [equivalent to] |[phi](t)> <[phi](t)| [member of] R with |[phi](0)> [approximately] |[phi]([tau])>. For [~.C] an absolute phase factor of |[phi](t)> is immaterial, since the curve [~.C] is defined via the evolution of the projection operator |[phi](t)><[phi](t)|. The geometric phase is a property of Ray space, where [~.C] is closed due to the equivalence of the initial and the final state |[phi](0)> [approximately] |[phi]([tau])> and can be calculated via a surface integral over the area enclosed by [~.C].

One can find many different curves C', C",... in Hilbert space differing by a phase factor [e.sup.i[alpha](t)] and yielding the same curve in Ray space under the projection map [pi]. On the other hand for a given curve in Ray space, there exists one distinct curve in Hilbert space fulfilling the parallel transport conditions, namely that two neighbouring states |[phi](t)> and |[phi](t + dt)> in H have the same phase, that is to say, <[phi](t)|[phi](t + dt)> is real and positive. This implies by Taylor expansion that <[phi](t)|[d/[dt]]|[phi](t)> = 0 [18]. For this curve the dynamical phase vanishes as one can verify by inserting the Schrodinger equation H(t)|[phi](t)> = ihd/dt|[phi](t)> into the parallel transport condition.

The concept can be extended to apply to open paths in Ray space where |[phi]([tau])> [??] |[phi](0)> by closing the curve by a geodesic, i. e., a path in Ray space with the shortest distance from |[phi]([tau])> <[phi]([tau])| to |[phi](0)> <[phi](0)|. Then one obtains a well-defined surface area enclosed by the path generated by the evolution of the system plus the geodesic closure. This surface provides an expression for the geometric phase, which has been proven by Samuel and Bhandari [4].

To sum up, for a general evolution of a quantum state the state obtains a dynamical phase dependent on the energy and time as well as a geometric phase only dependent on the subjacent geometry of state space. For special Hamiltonians which fulfill the parallel transport conditions the dynamical phase vanishes, which is also the case when the state is transported along a geodesic. An example of the latter is an evolution along a great circle on a sphere for a two-level system which we will encounter in the forthcoming discussion.

3. Interferometric Setup

Due to Feynman [16] the description of any two-level quantum system is equivalent to the description of a spin-1/2 particle. Exploiting this equivalence there is in principle no difference between manipulations in the spin space of neutrons with the orthogonal basis {|[up arrow]>, |[down arrow]>} as eigenstates of [[sigma].sub.z], and momentum space with {|k>, |k'>} as orthogonal basis vectors corresponding to two directions of the neutron beam in an interferometer. In both cases one can assign a geometric phase to the particular evolution of the initial state. An even more appropriate description for the interferometric case for the forthcoming discussion is in terms of "which-way" basis states {|p>, |[p.sup.[perpendicular to]]>}, namely, if the neutron is found in the upper beam path after a beam-splitting plate it is said to be in the state |p>, or in the state |[p.sup.[perpendicular to]]>, if found in the lower beam path. In case of a 50 : 50 beam-splitting of the incident (neutron) beam into a transmitted beam and a reflected beam, the associated wave vector after the beamsplitter can be written as an equally weighted coherent superposition of the two paths |q([delta])> [equivalent to] 1/[square root of 2](|p> + [e.sup.i[delta]]|[p.sup.[perpendicular to]]>) with the relative phase [delta] [member of] R depending on the particular physical realization of the beamsplitter.

For testing the spatial geometric phase we use a double-loop interferometer (Fig. 1), where the incident unpolarized neutron beam |[psi]> is split up into a diffracted reference beam |[[psi].sub.ref]> and a transmitted beam |[[psi].sub.t]>. The transmitted beam is subjected to further evolution in the second loop of the interferometer by use of beam-splitters (BS1 and BS2), an absorber (A) with transmission coeficient T and the phase shifter PS2 generating a phase shift of [e.sup.i[phi]1] on the upper and [e.sup.i[phi]2] on the lower beam path, respectively, yielding the final state |[[psi].sub.f]> = U|[[psi].sub.t]> = U |p>. Here |[[psi].sub.t]> = |p> since before the beam splitter BS1 the beam is clearly localized as seen from the second loop of the interferometer. The unitary matrix U = U(T, [[phi].sub.1], [[phi].sub.2]) comprises all the manipulations in the second loop:

[FIGURE 1 OMITTED]

|[[psi].sub.t]>[BS.[right arrow]] [1/[square root of 2]](|p> + |[p.sup.[perpendicular to]]>)[A.[right arrow]] [1/[square root of 2]](|p> + [square root of T]|[p.sup.[perpendicular to]]>; [PS2.[right arrow]] [1/[square root of 2]]([e.sup.i[phi]1]|p> + [square root of T][e.sup.i[phi]2]|[p.sup.[perpendicular to]]>) [equivalent to] U|[[psi].sub.t]> = |[[psi].sub.f]>. (3)

The geometric phase can then be extracted from the argument of the complex valued scalar product between the initial and the final state arg <[[psi].sub.t]|[[psi].sub.f]> (when removing dynamical contributions as will be discussed later). This is where the reference beam comes into play: |[[psi].sub.ref]> is not subjected to any further evolution, but is stationary apart from adding a phase factor [e.sup.i[eta]] by use of the phaseshifter PS1. |[[psi].sub.ref]> propagates towards the beamsplitter BS2 from the upper path, thus we can assert it to be in the state [e.sup.i[eta]]|p>. Then by the variable phase shift [e.sup.i[eta]] one can measure the shift of the interference fringes reflecting the phase difference between |[[psi].sub.ref]> and |[[psi].sub.f]>.

This preparation of the states is followed by the recombination of the two beams |[[psi].sub.f]> and |[[psi].sub.ref]> at the beamsplitter BS2 and the detection at the detector [D.sub.0] in the forward beam. This step can be described by the application of the projection operator |q> <q| = 1/2(|p> + |[p.sup.[perpendicular to]]>) (<p| + <p'|) (with [delta] = 0, which can always be achieved by an appropriate choice of the phase of the basis states) to |[[psi].sub.f]> as well as to |[[psi].sub.ref]>:

|[[psi]'.sub.f]> = |q><q|[[psi].sub.f]> = K([e.sup.i[phi]1] + [square root of T][e.sup.i[phi]2])|q> |[[psi]'.sub.ref]> = |q><q|[[psi].sub.ref]> = K|q>, (4)

where K is some scaling constant.

The intensity I measured in the detector [D.sub.0] is proportional to the absolute square of the superposition |[[psi]'.sub.f]> + [e.sup.i[eta]]|[[psi]'.sub.ref]>:

I [proportional] |([e.sup.i[eta]] + [e.sup.i[phi]1] + [square root of T][e.sup.i[phi]2])|q>|[.sup.2] = <[[psi]'.sub.ref]|[[psi]'.sub.ref]> + <[[psi]'.sub.f]|[[psi]'.sub.f]>+2|<[[psi]'.sub.ref]|[[psi]'.sub.f]>|cos([eta] - arg<[[psi]'.sub.ref]|[[psi]'.sub.f]>). (5)

We notice a phase shift of the interference pattern by arg <[[psi]'.sub.ref]|[[psi]'.sub.f]>. This phase shifts corresponds to the Pancharatnam connection [5] between the state |[[psi]'.sub.f]> and the state |[[psi]'.sub.ref]> = |q><q|[[psi].sub.t]> = |q> from which we can extract the geometric phase. Explicitly we obtain

[phi] = arg<[[psi]'.sub.ref]|[[psi]'.sub.f]> = [[[[phi].sub.1] + [[phi].sub.2]]/2] - arctan [tan([[DELTA][phi]]/2)([1 - [square root of T]]/[1 + [square root of T]])]. (6)

where [DELTA][phi] [equivalent to] [[phi].sub.2]-[[phi].sub.1]. The geometric phase is defined as [6]

[[phi].sub.g] [equivalent to] arg<[[psi]'.sub.ref]|[[psi]'.sub.f]> - [[phi].sub.d], (7)

where [[phi].sub.d] denotes the dynamical part. From Refs. [13] and [15] we know that the dynamical part stemming from the phase shifter PS2 is given by a weighted sum of the phase shifts [[phi].sub.1] and [[phi].sub.2] with the weights depending on the transmission coeficient T. In particular we have

[[phi].sub.d] = [[[phi].sub.1] + T[[phi].sub.2]]/[1 + T], (8)

which vanishes by an appropriate choice of phase shifts and transmission, i. e., [[phi].sub.d] = 0 for [[phi].sub.1]/[[phi].sub.2] = -T.

By varying the relative phase [DELTA][phi] from 0 to 2[pi] and setting [[phi].sub.d] = 0 the geometric phase [[phi].sub.g] can be plotted over [DELTA][phi] (Fig. 2).

[FIGURE 2 OMITTED]

4. Bloch-Sphere Description

For every two-level system we can use the Bloch-sphere for depicting the state vectors and evolutions thereof as points and curves on a sphere. Then the results obtained above, i. e., the shift of the interference pattern in Eq. (5) without dynamical contributions, should be equal to the (oriented) surface area enclosed by the paths of the state vectors on the Bloch-sphere, or, equivalently, to the solid angle traced out by the state vectors as seem from the origin of the sphere.

As we can observe in Fig. 3, the north pole of the sphere can be identified with a state with well known path, i. e., an eigenstate of the observable |p><p|. After the beam splitter BS1 the state |[[psi].sub.t]> evolves to an equal superposition of upper path and lower path, therefore the evolution on the Bloch sphere is given by a geodesic from the north pole to the equatorial line (the particular point on the equator is arbitrary due to the arbitrary choice of the phases of the basis vectors).

[FIGURE 3 OMITTED]

The absorber changes the weights of the superposed basis states, in particular for the extremal values of T parameterized by the angle [theta] with T = [tan.sup.2] [theta]/2, we end up either again with an equally weighted superposition for no absorption (T = 1 or [theta] = [pi]/2) or the state is now on the north pole for total absorption (T = 0 or [theta] = 0), since in the latter scenario we know the particle has taken the upper path when detecting a neutron in [D.sub.0]. For T [member of] (0, 1) the state is encoded as a point on the geodesic from the north pole to the equatorial line.

Due to the phase shifter PS2 we obtain a relative phase shift between the superposing states of [DELTA][phi] = [[phi].sub.2] - [[phi].sub.1]:

[1/[square root of 2]](|p> + |[p.sup.[perpendicular to]]>) [right arrow] [1/[square root of 2]](|p> + [e.sup.i[DELTA][phi]]|[p.sup.[perpendicular to]]>). (9)

This can be depicted as an evolution along a circle of latitude on the Bloch sphere with periodicity of 2[pi].

The recombination at BS2 followed by the detection of the forward beam in [D.sub.0] is represented as a projection to the starting point on the equatorial line, i. e., we have to close the curve associated with the evolution of the state by a geodesic to the point |q> <q| on the sphere as discussed for non-cyclic paths in Sec. 2. As for the reference state |[[psi].sub.ref]> we note that the phase shift of [eta] has no impact on the position of the state on the Bloch sphere, it stays at the north pole. Due to the recombination at BS2 and the detection the state is also projected to |q> <q| contributing to the forward beam incident to the detector [D.sub.0].

The paths are depicted in Fig. 4 in detail for cyclic Fig. 4(a) as well as non-cyclic evolution, Fig. 4(b). For a relative phase difference greater than [pi]/2 we have to take the direction of the loops into account. In Fig. 4(b) the first loop is traversed clockwise, whereas the second loop is traversed counterclockwise yielding a positive or negative contribution to the geometric phase, respectively.

The numerical calculation of the surface area F enclosed by the path traversed by the neutron is straight-forward by evaluating the solid angle [OMEGA] = [[integral].sub.F] sin [theta]d([DELTA][phi])d[theta] via a surface integral and using [[phi].sub.g] = -[OMEGA]/2. For the cyclic case this integral can be solved easily by calculating the segment on the sphere according to Fig. 4(a) to obtain [[phi].sub.g] = -[OMEGA]/2 = [pi](cos[theta]-1). For the non-cyclic case no analytic expression has been found to compare the results with the phase shift of the interference fringes appearing in Eq. (5). However, the numerical results are equivalent and agree with Fig. 2.

[FIGURE 4 OMITTED]

5. Experimental Results

For an experimental test of the spatial geometric beam we have used the double-loop perfect-crystal-interferometer installed at the S18-beamline at the high-flux reactor at the Institut Laue-Langevin, Grenoble [17]. A schematic view of the setup is shown in Fig. 1. Before falling onto the skew-symmetric interferometer the incident neutron beam was collimated and monochromatized by the 220-Bragg reflection of a Si perfect crystal monochromator placed in the thermal neutron guide H25. The wavelength was tuned to give a mean value of [[lambda].sub.0] = 0.2715 nm. To eliminate the higher harmonics we have used prism-shaped silicon wedges. The beam cross-section was confined to (5 X 5) [mm.sup.2] and an isothermal box enclosed the interferometer to achieve reasonable thermal environmental isolation. For the phase shifters parallel sided aluminium plates have been used: 4 mm inserted in the first loop as PS1 and 4 mm and 0.5 mm, respectively, in the second loop for PS2 yielding a ratio of 1/8 for [[phi].sub.1]/[[phi].sub.2]. To avoid dynamical phase contributions a gadolinium solution with T = 0.118 [19] has been used as an absorber. For a comparison with the theoretically predicted values one has to keep in mind that the contrast reflecting the coherence properties is different between each of the beams in the interferometer. Accounting for this experimental fact in the theoretical derivation of the geometric phase we notice a slightly flattened curve in Fig. 5 compared to Fig. 2. Nevertheless, one can recognize the increase in geometric phase for [DELTA][phi] [member of] [0, [pi]/2] due to the positively oriented surface followed by a decrease due to appearance of a counterclockwise traversed loop on the sphere yielding a negative phase contribution.

[FIGURE 5 OMITTED]

6. Conclusions

In summary we have shown that one can ascribe a geometric phase not only to spin evolutions of neutrons, but also to evolutions in the spatial degrees of freedom of neutrons in an interferometric setup. This equivalence is evident from the description of both cases via state vectors in a two dimensional Hilbert space. However, there have been arguments contra to the experimental verification in [13] which we believe can be settled in favour of a geometric phase appearing in the setup described above. The twofold calculations of the geometric either in terms of a shift in the interference fringes or via surface integrals in an abstract state space allows for a geometric interpretation of the obtained phase shift.

Acknowledgments

The authors acknowledge support from the Austrian Science Foundation, Projects No. F1513 and F1514.

7. References

[1] M. V. Berry, Proc. Roy. Soc. London Ser. A 392, 45 (1984).

[2] F. Wilckzek and A. Zee, Phys. Rev. Lett. 52, 2111 (1984).

[3] Y. Aharonov and J. Anandan, Phys. Rev. Lett. 58, 1593 (1987).

[4] J. Samuel and R. Bhandari, Phys. Rev. Lett. 60, 2339 (1988).

[5] S. Pancharatnam, Proc. Indian Acad. Sci. A 44, 247 (1956).

[6] N. Mukunda and R. Simon, Ann. Phys. 228, 205 (1993).

[7] N. Manini and F. Pistolesi, Phys. Rev. Lett. 85, 3067 (2000).

[8] Y. Hasegawa, R. Loidl, M. Baron, G. Badurek, and H. Rauch, Phys. Rev. Lett. 87, 070401 (2001); Y. Hasegawa, R. Loidl, G. Badurek, M. Baron, N. Manini, F. Pistolesi, and H. Rauch, Phys. Rev. A 65, 052111 (2002).

[9] A. Uhlmann, Rep. Math. Phys. 24, 229 (1986).

[10] E. Sjoqvist, A. K. Pati, A. Ekert, J. S. Anandan, M. Ericsson, D. K. L. Oi, and V. Vedral, Phys. Rev. Lett. 85, 2845 (2000).

[11] S. Filipp and E. Sjoqvist, Phys. Rev. Lett. 90, 050403 (2003); S. Filipp and E. Sjoqvist, Phys. Rev. A 68, 042112 (2003); S. Filipp and E. Sjoqvist, e-print quant-ph/0403204 (2004).

[12] E. Sjoqvist, Phys. Rev. A 63, 035602 (2001).

[13] Y. Hasegawa, M. Zawisky, H. Rauch, and A. I. Ioffe. Phys. Rev. A 53, 2486 (1996).

[14] H. Rauch and S. A. Werner, Neutron Interferometry, Clarendon Press, Oxford (2000).

[15] A. G. Wagh. Phys. Rev. A 59, 1715 (1999).

[16] R. P. Feynman, F. L. Vernon, and R. W. Hellwarth. J. Appl. Phys. 28, 49 (1957).

[17] M. Zawisky, M. Baron, R. Loidl, and H. Rauch, Nucl. Instr. Meth. A 481, 406 (2002).

[18] These are natural parallel transport conditions analogue to the claim for constant length and direction of a parallel transported vector on a sphere in the example from geometry above.

[19] The deviation from the ratio 1/8 is due to the different absorptions of the 4 mm and 0.5 mm Al-phaseshifters.

Stefan Filipp

Atominstitut der Osterreichischen Universitaten, Stadionallee 2, A-1020 Vienna, Austria

and

Institut Laue Langevin, Boite Postale 156, F-38042 Grenoble Cedex 9, France

Yuji Hasegawa

Atominstitut der Osterreichischen Universitaten, Stadionallee 2, A-1020 Vienna, Austria

Rudolf Loidl

Atominstitut der Osterreichischen Universitaten, Stadionallee 2, A-1020 Vienna, Austria

and

Institut Laue Langevin, Boite Postale 156, F-38042 Grenoble Cedex 9, France

and

Helmut Rauch

Atominstitut der Osterreichischen Universitaten, Stadionallee 2, A-1020 Vienna, Austria

Accepted: August 11, 2004

Available online: http://www.nist.gov/jres