In-beam superfluid-helium ultracold neutron source for the ESS

1.Introduction

Intense sources of ultracold neutrons (UCN) are a prerequisite to a variety of projects at the high-precision frontier of low-energy particle physics. Among these experimental efforts are searches for a non-vanishing permanent neutron electric dipole moment (EDM), accurate measurements of the neutron beta-decay lifetime and asymmetries, and searches for dark matter and dark energy using gravitational spectroscopy methods [6].

Most modern UCN sources are of the “superthermal” type, first proposed in 1975 by Golub and Pendlebury [11]. There, in contrast to UCN production by neutron moderation, cold neutrons are “converted” to UCN by imparting nearly their entire kinetic energy to elementary excitations – most commonly phonons – of the source medium in single scattering events. When the medium is kept at low temperature, up-scattering of UCN back to higher energies is suppressed, enabling UCN densities larger than those possible by neutron thermalization in the medium. Shortly after the 1975 proposal, superfluid ⁴He (He-II) and solid deuterium (SD₂) were identified as particularly suitable source media for neutron conversion [9,12], and have since been applied in various existing UCN sources.

More recently, solid α-oxygen [14,23] and solid α-¹⁵N [33] have been investigated as possible UCN production media. Like in He-II and SD₂, neutron conversion in these media relies on coherent excitations, magnons and phonons in the former, and only phonons in the latter case. An alternative method of UCN production would employ a gradual slowdown of neutrons via multiple incoherent inelastic scattering events that excite local low-energy modes. One such recently-considered moderation scheme proposes to employ the zero-field splitting of the magnetic triplet state of molecular oxygen to cool neutrons in a multi-step cascade of constant energy decrements [51]. Realization of this scheme might place the paramagnetic molecular oxygen as a guest within the interstices of a host solid, such as fully deuterated clathrate hydrate, a weakly-neutron-absorbing heavy-water network. Assessment of the capabilities of such materials for the production of UCN and very cold neutrons (VCN) is a topic of current research within the HighNESS project [34], and is complementary to work on He-II based UCN sources.

While the exploration of new media for UCN production, described in the previous paragraph, is still at an early stage of development, the concept of an “in-beam” UCN source based on He-II as a well-established conversion medium is rather mature nowadays. It is the purpose of the present paper to point out possibilities for a competitive in-beam UCN source at the European Spallation Source (ESS), adapted to take advantage of a high-intensity liquid-deuterium (LD₂) moderator accessible under large solid angle through a wide opening in the shielding monolith, aptly-named the Large Beam Port (LBP). Details of the relevant ESS infrastructure may be found in Ref. [50].

2.Superfluid-helium converters: “in-pile” or “in-beam”?

Unique among all possible UCN production materials, ⁴He is the only stable nuclide with vanishing neutron absorption cross section. Consequently, isotopically pure He-II-based converters enable UCN storage lifetimes τ several orders of magnitude larger than are possible in all other media. This makes it possible to gradually build up a high density of UCN in the converter, which may then be used to great advantage in a variety of experimental situations.

The upper limit to τ is set by the beta decay mean free neutron lifetime τβ≈880s. In order to achieve storage times in superfluid helium that are of this order, it is necessary to suppress strongly temperature dependent neutron up-scattering. For He-II at temperatures T⩽1K, the rate constant for up-scattering is [8,10]

Broadly speaking, there are two possible implementations of He-II-based neutron converters. The first places the He-II “in-pile” in the intense field of neutrons pre-moderated by a cold source located near a primary neutron source, such as the core of a nuclear reactor or a spallation target [20,26,40]. This is the standard location for neutron cooling devices, in which the large cold-neutron flux leads to a high UCN production rate. Unfortunately, the large heat load on the converter from neutrons and gamma radiation makes it unrealistic to maintain temperatures below 1K near a high-power primary neutron source, even when a powerful cryogenic plant is available [22,39]. Therefore, in this configuration, the converter is directly connected to an extraction guide, from which UCN stream continuously, delivering a high steady flux of UCN either to an experiment operated in a flow-through mode or to several storage experiments in a time-shared mode [43]. However, it should be noted that, if the primary neutron source is not too strong [41], temperatures well below 1K are still possible in-pile. In such cases, the UCN source can partly compensate for a lower primary neutron source power by making use of a UCN valve at the extraction guide entrance to build up, or accumulate, a high UCN density in the converter [32,54].

In fact, such an accumulation concept, first pointed out in Ref. [12], is the basis for the second possible implementation of a He-II-based UCN source, in which the converter is irradiated “in-beam” by a beam of cold neutrons while situated far away from the hot zone of the primary source. In this configuration, the significantly lower heat load on the converter makes temperatures below 0.7K feasible without excessive cryogenic efforts. Despite the smaller solid angle under which the cold neutron beam can irradiate the converter, and hence the lower UCN production rate, UCN densities accumulated in-beam can compete with in-pile converter densities.

As will be discussed in greater detail in the next section, this competitiveness requires the suppression of neutron absorption and UCN losses due to wall collisions to the same degree as up-scattering. However, the complexity of this task is substantially reduced for a converter placed far from the primary neutron source, where its accessibility enables not only adaptation of the neutron transport system and UCN source to match the specific needs of an experiment, but makes it possible to supply UCN thereto with a relatively short guide.

In fact, a UCN guide can even be rendered unnecessary in the extreme variant of an in-beam source in which the experiment is immersed in the superfluid-helium of the converter itself. This strategy avoids UCN transport losses and, in some cases, offers additional advantages, using the helium either as a neutron detection medium, a means to increase resistivity and prevent electrical breakdown in neutron EDM searches, or a combination of the two [1,17,46]. The application of these concepts to multiple UCN storage chambers is also discussed in Ref. [5].

3.Basic principles of an in-beam UCN source

There are three main components to an in-beam UCN source: a cold neutron moderator, a neutron optical delivery system (NODS), and a converter vessel, filled in this case with He-II. The problem of optimizing the UCN production according to a particular set of desired characteristics requires a careful consideration of not only the many parameters describing these three components, but also the interplay between them. The present section is intended to serve as a point of entry into this problem, and will focus on a simplified picture of the in-beam source in order to discuss the concepts that are key to UCN production. Here, we will mainly consider the case in which one desires a high density of UCN in the source. The discussion is ordered proceeding in the upstream direction, that is, starting from the converter surface upon which cold neutrons are incident and working back toward the moderator.

Let us first consider the cold neutron beam incident on the converter. The parameters thereof relevant to UCN production are the converter surface area Aconv, which we take to be planar and perpendicular to the beam axis, and the “brilliance” of neutrons incident upon it, bconv (units of s−1cm−2sr−1A˚−1). After multiplying this latter quantity by the differentials dΩ and dλ, one obtains the number of neutrons per second and per unit area with wavelengths between λ and λ+dλ and in the solid angle dΩ about a particular direction that penetrate a unit area perpendicular to this direction. In general, the brilliance varies as a function of wavelength, direction, and position, and can be used to define the “flux spectrum” (units of s−1cm−2A˚−1) at the converter by integrating over the solid angle Ωconv of all incident neutrons, that is,

Let us consider the differential surface element dAconv at a particular position on the converter surface and restrict our attention to all neutrons of a given wavelength that are incident there from a particular direction. At the moderator surface, which we also take to be planar and perpendicular to the beam axis, there is a corresponding position and differential surface element dAmod – not necessarily equal to dAconv – from which are emitted all those neutrons that are capable of being transported by the NODS to dAconv. Let us call bmod the brilliance of such neutrons, and make the approximation that they are uniformly distributed within the solid angle Ωmod, which is small in practice. In this case, one may show that Liouville’s theorem, which dictates that the phase space density of a free neutron beam cannot increase, implies that

Although not true in general, for the special case that (dΦ/dλ)mod is uniform over the moderator surface, the flux spectra in Eq. (4) can be replaced by their area-averages. Equation (4) highlights the purely geometrical effect of the NODS in controlling the ratio Ωconv/Ωmod to potentially change the mean flux spectrum at the converter surface compared to that at the moderator. This point is important because, in comparing the performance of candidate NODS with different geometries, one should be careful to distinguish such geometric effects from other sources of loss, which can be included in the above formulation as

Before continuing along the main line of our discussion, we would like to address one point that may have occurred to the reader. In general, those portions of the flux spectrum at the physical emission surface of the moderator that are not transported to Aconv by the geometry of the NODS may nevertheless reach the plane of the converter surface. In principle then, one could capture more of the total rate of incident neutrons by increasing the size of Aconv. However, as we shall see later, increasing Aconv generally results in a decrease of the mean flux spectrum at the converter, and thus in an unwanted reduction in the UCN density of the source.

Having thus outlined the cold neutron beam transportation, let us now consider UCN production in the converter vessel, and briefly describe the typical fashion in which it is used. As was previously mentioned, the basis for the implementation of the in-beam source is the accumulation of UCN therein under irradiation by the cold neutron beam. In typical configurations, the accumulated ensemble of UCN is used by releasing it as needed to an experiment through a UCN valve, enabling cyclical operation of the source. This mode is particularly suitable for UCN storage experiments, which often detect in some way those neutrons that remain after storage. Thereafter, a new ensemble of accumulated UCN is released from the replenished source, and the cycle is repeated until the desired sensitivity to the effect being investigated is reached.

The time evolution of the build-up of UCN density in the source under irradiation by cold neutrons follows from the rate equation including both the UCN production and loss, and is

Let us first consider maximization of τ. Its inverse, the “rate constant” for UCN loss, is given by

While maximization of τ concerns primarily the converter vessel itself, maximization of p also entails a consideration of the moderator and NODS. There are two mechanisms that contribute to the production rate, i.e., p=pI+pII. The dominant “single-phonon” process (pI), in which a cold neutron with initial energy in a narrow range centered on E∗≈1.04meV, corresponding to wavelength λ∗≈8.9A˚, imparts nearly all its energy to the converter medium by exciting a single phonon. The resulting neutron can be stored if its final energy is less than the “well depth” Vconv, which is the difference in the neutron optical potentials of the converter walls and the He-II. Taking the walls to be beryllium, the well depth is Vconv=VBe−VHe≈233neV, and the production of storable UCN is kinematically allowed for cold neutrons within an interval of width ΔE≈30μeV centered on E∗, corresponding to a full range width of Δλ/λ∗≈1%. Contributions from all ΔE result in a UCN production rate density of [36]

When an experiment installed at the UCN source is filled, as described above, by releasing the UCN accumulated in the converter, the UCN therein are subsequently diluted over the combined volume of the source, guides, and experiment. A larger source volume is therefore often desirable. However, increasing the source volume, in which UCN production takes place, typically reduces the mean UCN production rate density. Nevertheless, the resulting larger total UCN production rate may still increase the total number of UCN that finally arrive in the experiment. In such a situation, the high transparency of He-II for cold neutrons is very helpful, as it enables an increase of the total UCN production rate simply by elongating the converter vessel along the beam direction. Depending on the incident solid angle Ωconv and the overall length of the converter Lconv, dilution of the UCN production rate density p due to beam expansion can be minimized by lateral confinement of the incident beam in the source through the use of a neutron guide of constant cross section Aconv running down the converter length. In general, the total UCN production rate in such a system follows from integration of p over the converter volume, and is

Limiting the incident flux spectrum to the previously discussed narrow range of wavelengths centered on 8.9A˚ greatly reduces the heat load and background associated with neutrons that anyway have lower efficiency to convert to UCN [31]. Considering then only single-phonon UCN production, the volume-averaged UCN production density is

In this expression, we have replaced Λconv(λ) by Λ(λ), which is the mean free path for a neutron in He-II alone. At 8.9A˚, Λ∗=Λ(λ∗)≈16.6m [42], so that even a converter several meters in length is “short” in the sense that Lconv≪Λ∗. Thus, the attenuation factor ≈1−Lconv/(2Λconv∗) stays close to unity, enabling a large source volume, and hence large PI, at marginal cost in UCN density. One could also consider forgoing the technical complication of immersing a neutron guide within the He-II converter by increasing the converter vessel cross sectional area along the beam direction, thus allowing for beam expansion. This might, however, reduce the number of UCN in the experiment, and one should also note that the UCN extraction rate from the source is proportional to the reciprocal of the converter vessel volume. Interestingly, inclusion of the PII contribution to P does very little to improve the total rate or the UCN density; in long converters, PII falls significantly below the aforementioned 30% observed in short converters due to a strong increase in the scattering cross section at smaller wavelengths, resulting in Λ(λ)≪Λ∗ [42].

Finally, let us examine the conceptually simplest NODS, analyzed in Ref. [31], in which the converter is directly connected to the moderator by a neutron guide of constant cross section Aguide=Aconv=Amod. The glancing angle θ of neutrons of wavelength λ transported in the guide is limited by the “critical angle” θc=θc(λ), above which the probability of reflection from the guide surface quickly drops to zero. When θc is small, the dependence can be expressed to good approximation as

4.Advanced neutron delivery systems for in-beam UCN sources

While transportation of the solid angle Ω(λ) in Eq. (14) concerns a point source inside the guide volume, in reality, the guide entrance is separated from the moderator emission surface by a gap. As the physical extent of the moderator surface is finite, increasing this gap will eventually result in the guide entrance ceasing to be “fully illuminated”, i.e., irradiated by the source from all directions in which neutrons can be transported at a given wavelength, at which point the transported solid angle becomes limited by geometry rather than by θc(λ). For a UCN source that converts primarily 8.9A˚ neutrons, the critical angle of interest is

Fig. 1.

Schematic of an elliptic nested mirror optical (NMO) system. M and M’, separated by 2f, are the common foci of a set of ellipses, one of which is indicated by a dashed line. The mirror surfaces are formed by truncating the ellipses to the common length l. See Ref. [53] for details of the mathematical construction. The two types of symmetry about the optical axis of the NMO are indicated in the upper corners; the views show the NMO surfaces projected onto a surface normal to the optical axis z. Throughout the paper, x denotes the horizontal direction transverse to z, and y the vertical direction.

The basic idea of an elliptic NMO, shown in Fig. 1, is to image neutrons by a single reflection from the moderator surface (M) to a target region (M’), i.e., the converter. The device consists of an assembly of mirrors of lengths l, whose surfaces are short sections of a set of ellipses, that are located in a plane between the common focal points M and M’, separated by twice the focal length, 2f, of the set of ellipses. A central absorber blocks the direct view onto the source along the beam, or “optical” axis, shielding the target region from fast neutrons and gamma radiation. The mirror assembly produces an “image” at M’ of the neutron beam at M, which becomes “sharper” with decreasing l/f. The extent of the image is increased (reduced) in size if the mirrors are located closer to M (M’). Thus, the mirrors shown in Fig. 1, which are halfway between M and M’, produce a non-magnified image of the moderator.

Figure 2 shows examples of two basic implementations of NMO, having either planar (Figs 2a and 2b) or toroidal (Fig. 2c) symmetry. In the latter case, the device is rotationally symmetric about the optical axis, with each mirror surface being a section of an ellipsoid of revolution. In this implementation, a neutron at M is transported to M’ by a single reflection. For NMO with planar symmetry, each mirror surface is a section of an ellipse having local translational symmetry transverse to both the optical axis and the mirror surface normal. In such a construction, two planar NMO with a relative orientation of 90∘ about the optical axis are required, as shown in Figs 2a and 2b. Such systems image the beam by two reflections, one for each transverse dimension. A first experimental demonstration of a single-planar elliptic NMO and a discussion of further applications of NMO for neutron delivery and beam focusing can be found in Ref. [16]. To be noted is also that the NMO discussed here are conceptually simpler than (and involve only half of the reflections occurring in) Wolter imaging systems. Initially developed for X-rays, they focus a parallel beam by double reflections off a combination of parabolic with hyperbolic or elliptical mirrors. Adapted to slow neutrons [29], a microscope based on three nested coaxial Wolter mirrors has already been demonstrated [24].

One major advantage of NMO over neutron guides is their ability to efficiently transport neutrons with larger angles to the beam axis, which, recalling that Ω(λ)∝m2, constitute a large portion of the beam. Furthermore, the NMO transports all neutrons with only one (toroidal) or two (double-planar) reflections, so that, for each range of angles from the optical axis that is covered by one mirror section of the NMO, the transport efficiency varies only linearly (toroidal) or quadratically (double-planar) with imperfect reflectivity, and not in higher powers as for large angles in a guide.

Fig. 2.

Examples of NMO types for neutron transport by “imaging”: a) two separated planar mirror assemblies; b) two intersecting planar assemblies (double-planar); c) toroidal assembly.

Another advantage of NMO over neutron guides is the ability of the former to tailor the beam spectrum to the desired application (see Ref. [52] for several illustrations). This is possible because reflections at a given surface occur only within a narrow band of angles for which they are kinematically possible, allowing each mirror to be optimized to reflect those neutrons by which it is illuminated. This opens up new possibilities for producing a short-wavelength cutoff to the transported beam spectrum. Setting this cutoff to a value slightly below 8.9A˚ would remove the main part of the spectrum from the beam incident on the converter and hence eliminate the heat load and background associated with neutrons that increases UCN production by only 30% at most (see the discussion of multi-phonon UCN production in Section 3). Note that achieving a common cutoff requires the m-values of the supermirrors to decrease with decreasing semi-minor axis, which is also economically interesting, as lower-m mirrors require fewer layers in the multi-layer structure.

Going a step further, one could even prepare monochromatic beams. Concerning UCN production, NMO equipped with bandpass supermirrors [7,25,38] that have high reflectivity for only a narrow range of neutron wavelengths centered on 8.9A˚ are of particular interest. First, in addition to the common short-wavelength cutoff discussed before, such mirrors would also remove most of the longer wavelengths. Second, bandpass supermirrors have higher reflectivity at the target wave-vector transfer than broadband supermirrors, corresponding to better transport efficiency for 8.9A˚ neutrons. Moreover, they consist of far fewer bilayers, making them cheaper and easier to manufacture. The use of an NMO equipped with such supermirrors would be significantly more efficient than, for example, a velocity selector at the end of a guide to “clean” the spectrum incident on the converter [31].

The advantages of NMO do not end there however. Using NMO with polarizing broadband supermirrors, which would not introduce additional efficiency losses as might a dedicated polarization device, would enable polarization of the cold neutron beam incident on the converter. As has been experimentally demonstrated [46], due to the absence of nuclear spin or atomic magnetic fields in ⁴He, He-II-based UCN production proceeds without “spin-flipping”, allowing cold neutrons to maintain their polarization as they are converted to UCN. Thus, a beam of polarized cold neutrons produces an ensemble of polarized UCN, necessary for experiments like the multi-cell approach to the neutron EDM search described in Ref. [5]. Additionally, for those experiments that permit it, wall losses may be significantly reduced by surrounding the storage container with a multipole magnet [55].

Assuming no drastic decrease in the mean flux spectrum at the converter results from increasing the area of the beam Abeam at the converter entrance, such a modification has several benefits for UCN production. First, a larger Abeam results in a larger total rate of neutrons incident on the converter, and thus a larger total UCN production rate. Since the mean free path of UCN in the converter vessel scales with Abeam, a larger beam also reduces the rate constant of UCN wall collisions, which, according to Eq. (8) with Eq. (9), increases the saturated UCN density. A larger Abeam also facilitates implementation of a multipole magnet around the converter, as there is more space for either the placement of a larger number of poles, which increases the phase space of storable “low-field seeking” polarized UCN, or the placement of stronger poles at lower multipole order. For an elongated source, a larger Abeam also reduces the number of reflections during beam transport through the converter, and hence reduces beam attenuation due to imperfect reflectivity at the converter walls. Furthermore, a larger Abeam increases the source volume, and reduces the severity of dilution losses for UCN extraction to an external volume.

Fig. 3.

Partially integrated brilliance transfer B, defined in Ref. [16], for a single-planar elliptic NMO scaled in proportion to the separation 2f, where f is the focal length (see the text). Solid and dashed curves correspond to simulations with and without gravity, respectively.

The size of the beam extracted from the emission surface of the moderator, which is at most as large as the physical area of the moderator surface, influences the transport efficiency of an NMO and determines its optimal dimensions. To illustrate this, we show in Fig. 3 the results of Monte-Carlo simulations performed for a single-planar elliptic NMO with its mirrors oriented to reflect neutrons in the vertical plane. There, the “partially integrated brilliance transfer” B (see Ref. [16] for a precise definition) is plotted on the vertical axis as a function of focal length. B gives a quantitative measure of the device’s ability to image neutrons from a horizontal slit of height h centered on the focal line M (not a focal point as in the 2-D picture of Fig. 1) onto a slit of the same height h at the focal line M’, and separated from M by twice the focal length, 2f. For example, B=1 would corresponds to 100% neutron transport from slit to slit. The three pairs of curves in Fig. 3 correspond to simulations performed for three different values of h, in which the dimension of the NMO in the direction transverse to the optical axis is scaled in proportion to the focal length f, i.e., constant eccentricity.

The solid curves in Fig. 3 all follow the same general trend: for a given slit height h, the efficiency B monotonically increases with increasing f up to a maximum, whereafter it monotonically decreases. This maximum occurs at increasing f for increasing h. Much of this behavior can be explained by several considerations. First, for purely geometrical reasons, some neutrons starting with an offset from the focal line either pass the NMO unreflected or undergo two reflections, and hence do not arrive within the target slit. The proportion of such neutrons scales roughly with h/bn, where bn is the semi-minor axis, i.e., in the direction transverse to the optical axis, of the nth ellipse. Thus, the effect is initially greatest for the largest slit (green curve), and decreases as the size of the NMO increases, i.e., with increasing f. Another effect favoring large NMO is the decreased importance of the “divergence hole” for small angles; see Ref. [16] for an explanation. Second, large losses are present for large separations because, there, neutron trajectories are significantly curved by gravity. “Switching off” gravity, see the dashed curves in Fig. 3, results in a restoration of linear neutron trajectories, and hence in increased efficiency with increasing NMO size.

From these simulations, based on NMO using m=6 broadband supermirrors with 67% edge reflectivity, it is apparent that transport of a beam with extent as great as 20cm to a position farther than 30m from the moderator surface is possible with efficiency in excess of 75%. Such a separation places the NMO completely outside the largely-inaccessible region within the shielding bunker at the ESS, and still-larger distances would increase the efficiency to its maximum at 2f≈100m. However, such additional gains come at the price of an increased technical complexity and monetary cost associated with constructing a proportionally larger NMO.

5.In-beam UCN source implementation at the ESS

Having discussed in the previous two sections many of the concepts that are key to optimizing in-beam UCN production, we shall now apply them in order to estimate the UCN densities that may be possible at the ESS using NMO to transport cold neutrons to a He-II converter vessel. As discussed previously, the ideal location for such a UCN source would make use of the LBP, through which moderators in the vicinity of the spallation target wheel can be viewed with large solid angle (see “option 5” in Ref. [50]). During initial operation of the ESS spallation source, only the compact “bi-spectral” moderator system above the target wheel will be implemented [49]. This system includes the flat para-hydrogen moderator, which has an emission surface area and an expected surface averaged brilliance at λ∗=8.9A˚ of

Figure 4 shows a simulated intensity map of neutrons with wavelengths near 8.9A˚ emitted by the surface of this larger moderator. Although its mean brightness is smaller by a factor of 1.5 in comparison to the para-hydrogen moderator, its surface area is larger by a factor of 40. Thus, efficient transport to the converter surface of cold neutrons leaving the LD₂ moderator would result in a much larger total UCN production rate in comparison to the para-hydrogen moderator. For this reason, the LD₂ moderator would be preferred for the subsequent operational stages at the ESS. While the small size of the para-hydrogen moderator does not preclude its use for an in-beam UCN source, it would restrict the total UCN production rate, and thereby the possible scope of such a source.

Fig. 4.

Simulated intensity map at the ESS LD₂ moderator surface of neutrons with wavelengths near 8.9A˚. The data shown, which are normalized to the maximum, correspond to neutrons detected at 4m from the emission surface after having passed through a “pinhole” halfway between the detector and moderator. Note that an image produced by such a “camera obscura” is inverted, and that the map shown has been mirrored up-down and left-right to represent the moderator surface as it would be seen looking through the LBP. The maximum intensity occurs closer to the upper end, above which the target wheel is situated.

Extraction of the divergent beam at 8.9A˚ from the LD₂ moderator by an elliptic NMO, discussed in Section 4, requires a large solid angle. Although this is possible through the LBP, the rather “anisotropic” access geometry there makes this task less than straightforward. We note, however, that no horizontally reflecting mirrors close to the moderator are necessary, keeping the view of the moderator surface at neighboring beamlines unobstructed and thus available for other purposes. Let us first consider the horizontal direction. The maximal angle θx at which an NMO can view a surface element within a central region of the moderator depends on the width wmod of this region. Taking wmod=22cm, for example, a choice that is motivated further below, would correspond to θx=5.6∘, and would require an m=6.4 coating on the outermost NMO mirrors.

In the vertical direction, access is limited from the bottom by a 2m long horizontal shielding plate, which, unlike that leading away from the top of the moderator, is flat rather than tapered; see Fig. 5 in Ref. [50]. Therefore, in order to minimize losses of neutrons emitted downward, the bottom plate should be covered by a mirror starting close to the moderator. If one then places a second mirror, parallel to this, one obtains a vertical, “one-dimensional” guide of height hg and length lg (see Fig. 5). A planar elliptic NMO with horizontal mirrors can then extract neutrons from the “virtual source” defined by the end of this guide. To minimize the fraction of neutrons transported by multiple reflections in the guide, hg should be maximized, which is also in line with the aim of maximizing the viewable portion of the moderator surface. Choosing hg=22cm, which is slightly less than the height of the moderator emission surface (24cm), would permit full vertical illumination of the guide when there is a gap dg of several centimeters, likely necessary for technical reasons, between its entrance and the moderator surface.

Fig. 5.

Schematic of the proposed Neutron Optical Delivery System for the Large Beam Port at the ESS. The values for the parameters quoted below are those used in the Monte-Carlo simulations shown in Fig. 6. The distance between the moderator, at left, and the converter entrance, at right, is 35m. Neutrons are emitted at the LD₂ moderator surface from an area wmod×hmod, where the physical height hmod=24cm. Note that, as the horizontal edge of the viewed area on the moderator is not sharply defined, the simulations used wmod=24cm, slightly larger than stated in the text. Neutrons are initially guided in the vertical direction by two parallel, flat mirrors (lg=2.85m, hg=22cm) separated from the moderator by a gap of dg=10cm. The beam is imaged onto Aconv by two planar elliptic NMO that consist of mirrors of lengths l=0.5m. The semi-minor axis of the outermost mirror of both devices is b0,x=b0,y=1.85m. The vertically focussing NMO has a focal length of fy=16.025m and is placed halfway between Aconv and the end of the flat mirrors, whereas the horizontally focussing NMO (fx=17.5m) is centred between Aconv and the moderator. All mirrors of the NODS are coated with the same broadband m=6.4 supermirror, having 64% edge reflectivity. The mirrors of the NMO are deposited on 0.5mm thick silicon substrates, and the simulation accounts for neutron refraction and absorption by these substrates. A beam stop of area wb2=0.25m2 placed in front of the first NMO blocks the direct view of moderator at the converter. For the purpose of illustration, the eccentricity of the NMO and mirror curvatures have been exaggerated.

Although such a situation does not correspond to a true optimization of the complete NODS, let us consider the case in which lg is selected so that the vertical range of angles under which the NMO can view the guide exit is the same as the horizontal range of angles under which it can view the width wmod on the moderator; i.e., θy=θx=5.6∘. This will be satisfied when lg is slightly less than 3m. Supposing that this guide uses the same broadband m=6.4 supermirror as in the NMO, the majority of cold neutrons incident on the guide surfaces at the critical angle θc∗=5.6∘ would be transported by only a single reflection.

For this geometry, a double-planar elliptic NMO would be best adapted to image the beam, whose horizontal and vertical extents are defined at the moderator surface and guide end, respectively, onto the He-II converter. As shown in Fig. 5, two separate NMO with different focal lengths fx and fy are required. He-II converter vessels are most often cylindrical in shape, which is advantageous for the implementation of a magnetic multipole reflector. This motivates the choice of the viewed width of the moderator surface to be wmod=hg=22cm. However, the linear extensions wmod and hg must then satisfy wmod⩾hg⩾dconv=2Aconv/π in order to avoid under-illumination of the circular converter surface. Extraction of a rectangular beam of nearly the full moderator size would also be possible by a double planar NMO. However, as the angular range of view onto the whole surface would be reduced by the biological shielding near the moderator, this option would profit from also guiding the neutrons inside the shielding in the horizontal direction.

Fig. 6.

Simulated intensity distributions at the locations as labelled in Fig. 5, for R0=64% supermirror edge reflectivity. I1, I2, and I3 are the fractional intensities encompassed by the red circles of diameters d1=22cm, d2=15cm, and d3=7.5cm, respectively, normalized to the total intensity Itot, which is taken as unity at the moderator surface. (a) Intensity of neutrons emitted from the area wmod×hmod at the moderator for uniform brilliance within angular and wavelength ranges of −5.6∘⩽θx, θy⩽5.6∘ and 8.7A˚⩽λ⩽9.1A˚, respectively. (b) Intensity at the end of the two flat horizontal mirrors. The total intensity Itot here, and also at (c), is normalized to the total intensity at (a). The sharp vertical boundaries are due to the mirrors, while the horizontal spread in intensity is due to beam expansion. (c) Intensity at the converter entrance. The red circles have the same diameters as in (a). Within the largest circle, the intensity ratio is I1(c)/I1(a)=47% (for R0=64%), while for the smallest circle it is slightly better: I3(c)/I3(a)=53%. If I1(c)/I1(a) is further reduced, conservatively, by 15% on account of imperfections of the NMO and absorption and scattering at beam windows, we can estimate the transport efficiency as in Eq. (6), with the flux spectra replaced by their area averages, to be ϵtr=40% for a converter entrance window of diameter dconv=22cm. The lower fractional intensities refer to a second simulation preformed for ideal reflectivity of both NMO, R0=100%.

The implementation of a vertical guide close to the ESS spallation target is complicated by the high flux of fast neutrons there. Even for supermirrors on metallic substrates, developed for in-pile guide sections, sufficient radiation hardness still needs to be experimentally demonstrated. If such a guide is not feasible, beam extraction would have to be accomplished solely by NMO, in which case, the aforementioned geometrical restrictions and partial obstruction of the moderator view would reduce the possible UCN yield.

Using the neutron transport efficiency of the NODS shown in Fig. 5, estimated from the simulations shown in Fig. 6 to be ϵtr=40%, we can finally give estimates for UCN production induced by the single-phonon process in a He-II converter installed at a distance of 35m from the LD₂ moderator. The production rate density in a short converter with an entrance window of diameter dconv=22cm follows from using Eq. (5) with bmod=bLD2∗ from Eq. (17) and Ωmod=4×(5.6∘)2=0.038sr (Eq. (14) with γ=4), Eq. (6) with ϵtr=0.4 and Ωconv=Ωmod, and Eq. (10), resulting in

Maximal UCN production in the converter would be obtained by matching its entrance shape to that of the rectangular beam provided by the NODS, which is shown in Fig. 6 with a circular converter entrance. Placing an m=6.4 supermirror neutron guide at the converter walls would then transport the provided neutron beam through an elongated source, within which the volume-averaged UCN production density would be decreased by imperfect reflectivity and by transmission losses, Eq. (12). In a source with a cylindrical supermirror guide, the transported neutron flux would be reduced by a factor ≈π/4, cf. Eq. (14). For such a situation, with a converter having diameter 22cm and length 3m, and assuming a loss of UCN production on the order of 20% due to imperfect mirror reflectivity, the total UCN production rate PI and total saturated UCN number Nsat in the converter would be

6.Conclusions and outlook

The parameters of the UCN in-beam type source described in this paper, quoted in Eqs (18)–(20), make it a very attractive option for the community of UCN users. As the only current in-beam type, general-purpose UCN source project beyond ILL’s SuperSUN [3], it offers a saturated UCN density at the top of the range of other current projects, which all rely on UCN production in the vicinity of a strong primary cold neutron source. The high density is the result of UCN accumulation in He-II kept at a temperature so low as to render neutron up-scattering negligible, which is possible at a position far away from the primary neutron source, where only less than 1W cooling power is needed. This kind of UCN source is thus complementary to in-pile type projects, which generally possess a much higher UCN production rate. Quoting here, only for comparison, what is perhaps the most ambitious of these projects proposed for LANSCE [22], the “inverse-geometry spallation-driven UCN source” is designed to cool away 100W from 40L He-II at 1.6K, resulting in a UCN production rate density of p=5×104s−1cm−3 and a UCN flux of 5×108s−1 through a 5m long, 18cm diameter guide which would enable UCN densities up to 104cm−3 in a large external trap.

Not to be underestimated in any source comparison are also the practical advantages of the converter being placed far away from strong radiation fields of the primary source. There, the required cooling power is much reduced, the UCN source is accessible to trouble-shooting, UCN need to be transported only over short distances, and, finally yet importantly, nuclear licensing procedures will not directly hamper the project. As a very valuable property for experimenters, in-beam sources based on NMO have great flexibility to adapt the implementation of the optics and the He-II converter to the physics project, including the in-situ strategies presented in Ref. [5].

Particularly appealing in this respect is the possibility to trade maximal total UCN production against highest possible UCN density. A class of experiments that could take advantage of this possibility is the search for dark matter and dark energy via investigation of the gravitational quantum states of the neutron above a horizontal mirror. Current experiments by the qBounce collaboration [18] are performed using a beam of UCN and have already led to constraints on new physics [4,19]. A route towards further increase of sensitivity in such investigations consists in UCN trapping to prolong the time of controlled exposure of the quantum states to the field [30]. To this end, a high-density source from which one would extract only UCN in the desired quantum state [37] would be of great interest.

Focusing of the cold neutron beam incident on the He-II converter would enable a large gain in flux density and hence UCN production rate density, for which the NMO based magnifying-optics concepts presented in Ref. [16] should be further investigated. As the focused beam illuminates the converter under a solid angle Ωconv>Ωmod that may exceed by far the solid angle transportable by a guide equipped with even the most advanced supermirror coating, the converter has to be reasonably short to avoid dilution due to beam expansion behind the focal point of the optics. The correspondingly reduced beam size reduces the size of the converter and hence the mean free path of the UCN between wall collisions (which scales as Aconv). Highly reflective walls and/or magnetic trapping may keep losses due to the accordingly increased rate of UCN wall collisions well under control. Very helpful in this respect is the large emission surface of the LD₂ moderator foreseen at the ESS. Besides the aforementioned UCN trapping in quantum states, beam focusing by NMO could result in large additional gains with respect to the performance numbers given in Eqs (18) and (19), which might also open the door to other qualitatively new experiments.