In [20]:

Gamma emissions under grazing incidence and resonance-enhanced neutron standing waves in a thin film

X-ray or neutron standing waves can be generated above a mirror surface under conditions of total external reflection by incident and reflected beams. These beams have equal amplitude and superpose to form standing waves inside the sample. With the proper sample morphology, i.e. layer thickness close to standing waves period, the amplitude of the standing waves can be many-times amplified, causing emission of X-ray fluorescense in the case of X-rays or gammas in the case of neutrons by a layer of "heavier" material or stronger absorber. In this notebook we demonstrate how ESCAPE can be used for the modelling of Gamma emissions together with neutrons specular reflectivity.

As a basis of our model and obtained results we use the following publication:

Grazing incidence prompt gamma emissions and resonance-enhanced neutron standing waves in a thin film

Huai Zhang, P. D. Gallagher, S. K. Satija, R. M. Lindstrom, R. L. Paul, T. P. Russell, P. Lambooy, and E. J. Kramer Phys. Rev. Lett. 72, 3044

Theoretical background

The calculation of the intensity of Gamma emissions starts with the the Schrödinger equation which solution is the wave function in each layer given by


where zm corresponds to the the upper interface of the m-th layer, Tm and Rm are the trasmission and reflection Fresnel coefficients.

Continuity of the wave functions and their derivatives at each interface allows to find Fresnel coefficients for each layer using so-called matrix formalism.

The analysis of conservation of probability in the case of absorption, i.e. non-zero imaginary part of the scattering potential leads to the relation for the neutron capture rate of a single isotope element given by


where kz is the transverse wave vector component, ψ(z) - neutron wave function, ρim - imaginary part of scattering length density, c - normalization constant. The inegration is performed over the layer Z - boundaries.

Further details can be found in the original publication or in the third edition of Quantum Mechahics by Leonard I. Schiff.

Sample description

According to the publication the sample consists of the following layers starting from Si substrate: Si/Ni(600 nm)/PVCH(72nm)/Gd(5nm)/PVCH(72nm)/dPS(25nm).

First we should create materials instances. For the polymer layers we use the mass density as fit parameter. The imaginary part of the scattering legth density is very small for these layers and it's influence on the reflectivity is negligible compared to Gd layer. Plus, we reduce the number of fit parameters.


In [22]:


Initial thickness and roughness values for each layer are equal to the values provided in the publication or estimated from the provided SLD profile. In our model roughness from interfaces is included as in the standard matrix method statistically averaged intensity over the interfaces roughness. In the publication interface roughness was modelled by smearing the abrupt change in the scattering potential at each interface with a hyperbolic tangent function. It make sense if a very large interface roughness is expected, when a standard statistical averaging of thickness changes perpendicular to the interface plane

In [23]:

Specular reflectivity

Gamma-ray intensities and neutron reflectivity data were measured simulataneously keeping the illuminated area of the sample constant at  35×40mm2 with the beam resolution Δkz/kz0.058. The latter is used for averaging of specular reflectivity over the resolution. The FWHM of the resolution function is a function of kz in this case.

In [24]:
In [ ]:

Gamma emissions

In ESCAPE we introduced reftrans package which implements matrix formalism and allows to calculate Fresnel coefficients inside particular layer of a multilayer sample as well as the wave function. In the next cell we demonstrate how to use them for the calculation of Gamma emissions. The neutron capture rate should be finally multiplied by a constant C, representing the counting efficiency. According to the publication, C was one of the fit parameters together with structural parameters.

In [25]:

The imaginary part of the scattering potential uim(m) is constant inside the layer and the solution of the integral Zm|ψm|2dz is straightforward:

In [26]:
In [29]:
In [30]:
Best cost: 0.0001368 Status: maximum number of iterations reached