|
1.INTRODUCTIONNeutron optics started since the very early work of Fermi and Zinn in 19461. For comprehensive accounts, see Byrne2 and Sears3. Neutrons have the dual wave-particle nature and, therefore, the associated de Broglie wavelength which is, for thermal neutrons: λB=1.8 Å. The neutron is electrically neutral and unstable, its average life time being the order of 15 min. The particular phenomenon of neutron confinement, studied here, takes place within time intervals smaller than the average life time. In neutron optics, the wave function φ (x, t) for a slow neutron propagating in a material medium (or vacuum) satisfies the stationary (time-independent) Schrödinger equation: In Eq. (1) m is the neutron mass (m = 939.57 MeV/c2), V(x) is the effective scalar potential due to the medium, depending on the spatial coordinates, ħ is Planck’s constant and E is the total kinetic energy. For the current study we consider only thermal neutrons with E=0.025 eV. Eq. (1) is an equation, having analogies with scalar wave equations used for classical electromagnetic fields deduced from Maxwell’s equations, which can be exploited for characterizing the behaviour of a slow neutron beam. Neutron phenomena associated with total internal reflection in guides with macroscopic diameter (i.e.: centimeters) and neutron diffraction at shorter scales are well characterized provided that a constant potential be considered2-3. In Eq. (2), b is the coherent neutron-nucleus scattering amplitude (depending only on the material and isotope) and ρ is the density of nuclei per unit volume. Those neutron phenomena (having in mind similarities with the refractive index for the corresponding optical phenomenon of total internal reflection) can also be described in terms of an effective refractive index: Where: E0 is the energy of the neutron in the medium. The refractive index for neutrons plays a role in discussing confinement of neutron waves, under strict total internal reflection conditions. We remind here that for confined propagation in waveguides, the beam wavelength is close to characteristic length scales of the system through which the neutron wave propagates (say, transverse sizes of thin films or fibers). Eqs. (2)-(3) yield practical approximate descriptions provided that Bragg reflections do not occur and other scattering and absorption processes are negligible. There are some antecedents on slow neutron beam confinement in waveguides with small transverse dimensions. In 1973 Wames and Sinha4 proposed the coupling of neutrons in planar thin film waveguides. In 1984 and 1986, Álvarez-Estrada and Calvo studied theoretically the confinement of slow neutrons in cylindrical waveguides5-6. In 1994, Feng et al. demonstrated, experimentally, confined propagation of modes for slow neutrons in Ti films7. Important studies by Kumakhov et al.8 and Chen et al.9, both in 1992, established experimentally the focusing of slow neutrons by using polycapillary glass fibers. In all these studies, there are key parameters in the design of the guides, namely, dimensions, geometry, wavelength and energy of the neutron beam, materials to be used, values of the critical angle, in order to assure an effective confinement and low loss of energy. Since experiments with neutron beams require some special technologies, arranging neutron sources and transport to the required point for the study, it is quite necessary, previously, to develop theoretical models to implement numerical simulations. These theoretical models are not yet well developed, because of the necessity to deal with non-trivial wave phenomena in the framework of the Schrödinger equation. It has been discussed previously that the numerical solution of Schrödinger equation, for propagation in a medium with some specific potential and geometry, requires a carefully choice of the method. For example, Matlab is very popular these days but it requires time-consuming programs and has to be optimized. Moreover, the inclusion of boundary conditions is an added requirement since the waveguide has finite dimensions. These aspects are treated in the next section. We have designed10 neutron devices, behaving as couplers, which could possibly be applied in Boron Neutron Capture Therapy (BNCT) and for other purposes. 2.THEORETICAL MODELWe introduce a waveguide model that we have considered for the purpose of the numerical simulation of the thermal neutron beam confinement. Figure 1 displays the schematic representation for the neutron waveguide geometry considered in our study. An incoming neutron plane wave propagates, with V=0 in the air, at Z<0 and strikes a semi-infinite two-dimensional planar waveguide. The neutron beam propagates confined along the waveguide. The waveguide has aperture dimension x0, and confined modes are formed and propagate along the Z-axis (Z>0). The number of modes will depend on the particular geometry, transverse dimensions of the planar waveguide (say, aperture value) and incidence conditions. For suitably large incoming neutron energy, a high number of modes are formed (the neutron waveguide being a multimodal one). For simplicity, we have considered that the scalar repulsive potential (in the so-called clad) has an infinite value: V = +∞. Although, this is far from a real case, however, this approach is quite convenient in order to obtain and to interpret approximate results for the modes formation. An integral representation is obtained for the neutron wave function by using the formalism of the two-dimensional Green’s function. That integral representation includes four functions, namely, μ1 and μ2, defined at the walls of the entrance aperture of the guide, and μ3 and μ4, defined at the internal upper and lower walls of the waveguide, respectively. We omit here for brevity the mathematical equations derived from the theoretical model. We state that the four μ’s functions are uniquely determined by imposing the so-called Dirichlet boundary conditions (namely, the vanishing of the neutron wave function) at the four walls of the waveguide indicated above. Based upon them, approximate methods are applied to arrive to analytical solutions for odd and even propagation modes, respectively: Eq. (4) denotes the representation for the odd modes. For even modes: Eqs. (4)-(5) imply a semi-quantitative description of the wavenumber of those incoming neutrons which give rise to propagation modes. Both, odd and even modes are allowed into the waveguide. Figure 2 shows the results of the maximum number of modes allowed in these waveguides as a function of the neutron incoming energy for different values of the aperture of the guide: In the present study, since we are supposing an ideal waveguide with infinite repulsive potential in the clad (Dirichlet boundary conditions), one finds that the critical angle is the order of π/2. Thus, the waveguide allows almost any angle of incidence. 3.NUMERICAL SIMULATIONWe have developed an algorithm based upon the mentioned Green’s functions formalism and Dirichlet boundary conditions. Approximate analytical results are obtained by applying the analogies with classical optical diffraction and the Sommerfeld-Rayleigh integral as well as its Fourier transform. Moreover, we have implemented an approximate iterative method for obtaining μ3 and μ4, which are more important to account for the confined beam within the waveguide. Results are displayed in Fig. 3 (a) compared with FDTD results in Fig. 3 (b). REFERENCESFermi, E. and Zinn, W. H., Phys. Rev., 70 103
(1946). Google Scholar
Byrne, J.,
“Neutrons, Nuclei and Matter: An Exploration of the Physics of Slow Neutrons,”
IoP Publishing, Bristol
(1995). Google Scholar
Sears,V.F.,
“Neutron Optics,”
Oxford University Press, Oxford
(1989). Google Scholar
De Wames, R. E. and Sinha, S. K.,
“Possibility of guided-neutron-wave propagation in thin films,”
Phys. Rev. B, 917
–921
(1973). https://doi.org/10.1103/PhysRevB.7.917 Google Scholar
Alvarez-Estrada, R. F. and Calvo, M. L.,
“Neutron waveguides: a possible application of neutron optics,”
J. Phys. D: Appl. Phys., 17 475
–502
(1984). https://doi.org/10.1088/0022-3727/17/3/007 Google Scholar
Alvarez-Estrada, R. F. and Calvo, M. L.,
“Confined propagation of neutrons in solid waveguides (fibres),”
J. Physique, 45 C3-243
–C4-248
(1984). Google Scholar
Calvo, M. L. and Alvarez-Estrada, R. F.,
“Neutron fibre II: some improving alternatives and analysis of bending losses,”
J. Phys. D: Appl. Phys., 19 957
–973
(1986). https://doi.org/10.1088/0022-3727/19/6/012 Google Scholar
Feng, Y. P.,
“Direct observation of neutron-guided waves in a thin-film waveguide,”
Phys. Rev. B, 49
(15), 10814
–10817
(1994). https://doi.org/10.1103/PhysRevB.49.10814 Google Scholar
Kumakhov, M. A. and Sharov, V. A.,
“A neutron lens,”
Nature, 357 390
–391
(1992). https://doi.org/10.1038/357390a0 Google Scholar
Chen, H.,
“Guiding and focusing neutron beams using capillary optics,”
Nature, 357 391
–393
(1992). https://doi.org/10.1038/357391a0 Google Scholar
Alvarez-Estrada, R.F., Molina de la Peña, I. and Calvo, M.L,
“Focalizing slow neuron beams at and below micron scales: discussion on BNCT,”
Phosphorus, Sulfur and Silicon and the related elements, 193 64
–73
(2018). https://doi.org/10.1080/10426507.2017.1417300 Google Scholar
|