
This page provides interface to a program called BRL which can calculate multiple Bragg diffraction of Xrays by perfect crystals. BRL is based on the paper by S.Stepanov and A.Ulyanenkov "A New Algorithm for Computation of Multiple Bragg Diffraction", Acta Cryst. (1994) A50 579585. It can simulate up to 12wave dynamical Bragg diffraction of Xrays from a plateshaped crystal including the cases of Xray waves grazing along the plate surface and Bragg angles being close to 90^{o}.
Typically the calculations of multiple Bragg diffraction are reduced to the eigenvalue problem for a scattering matrix. If there are no grazing waves involved into the diffraction geometry so that the specular reflection effects could be neglected, then the size of the scattering matrix is 2N*2N for Nwave diffraction. Here the factor of two appears due to the two polarizations (sigma and pi) of Xrays. A description of respective theory can be found e.g. in the paper by V.G.Kohn [Phys. Stat. Sol. (a) 54, 375384].
However, if the diffraction geometry involves an Xray wave grazing along the crystal surface, the task becomes more complicated. For this case R.Colella [Acta Cryst. (1974) A30, 413423] suggested a theory where the calculations of multiple Bragg diffraction are reduced to the eigenvalue problem for 4N*4N scattering matrix. This holds even if only one of e.q. 24 waves is grazing which is obviously not very effective in terms of calculations.
In the paper by S.Stepanov and A.Ulyanenkov it was made possible to reduce the calculations of multiple Bragg diffraction to a generalized eigenvalue problem for 2(N+N_{s})*2(N+N_{s}) scattering matrix where N_{s} is the number of grazing waves. Thus, if there are no grazing waves, the matrix size is 2N*2N and if all of the waves are grazing it becomes 4N*4N. In some cases like only one of 24 waves is grazing the calculations are reduced dramatically.
BRL has been successfully used to simulate the applications of Bragg and Lauecase Renninger effect to Xray doubleplane collimation and the multiple diffraction effect in Xray surface back diffraction.
The BRL consists of two steps. At the first step one has to enter the parameters required to build incident wavevector. As one can see on Fig.1, the specification of two reflections with the reciprocal vectors h_{1} and h_{2} determines the axis OO' on which the incident wave vector k_{0} should originate. Then, there are the following three choices:
The other parameters on the form specify the crystal, the direction of external normal to the crystal surface, and the conditions to search for additional Bragg reflections that may simultaneously occur for given k_{0}. Note that in principle BRL should work for crystals that belong to any syngony. However, the program has been tested cubic crystals only. Any comments on the calculations for noncubic crystals will be greatly appreciated. The crystal code refers to respective record in the X0h database.
Once the form is submitted, Xray Server validates the input and attempts to build k_{0}. If everything is found to be correct and in particular it is ensured that k_{0} can enter the crystal (i.e. it makes an obtuse angle with the external surface normal), one proceeds to the Step2. Otherwise, an error is reported.
At Step2 one has to enter the crystal plate thickness, the incident Xray polarization, and the scan parameters. Note that when the direction of Theta_{2} is specified, the direction of Theta_{1} is given as a vector product k_{0} and Theta_{2}. If both "mirrored" k_{0} and k'_{0} were found to be able entering the crystal, the form provides a control for choosing one of the two.
The submission of Step2 form initiates the BRL. If no error is found in the input, the calculation starts and the result is returned as an HTML page containing data plot and a link to zipped data.
In the case when the number of scan points over either Theta_{1} or Theta_{2} is one, the data are written into the .DAT files (one per reflection) with the format:
x1, y1 x2, y2 .... xn, yn
Otherwise the .GRD files are produced which are the matrix files compatible with the Golden Software Surfer. These files can be used with GbuPlot too, after removing the first 3 header lines.
Click here to proceed to the BRL input form.
Please, report problems and send comments to  