This study proposes a new approximating formula for calculating gamma-ray buildup factors in multilayer shields. The formula combines the buildup factors of single-layer shields with products and quotients. The feasibility of the formula for reproducing the buildup factors was tested by using point isotropic buildup factors calculated with the SN1D discrete ordinates code as reference data. The dose buildup factors of single-, double-, and multilayer shields composed of water, aluminum, iron, and lead were calculated for a spherical geometry in the energy range between 10 MeV and 40 keV and for total thicknesses of up to 30 mean free paths. The calculation of the buildup factors takes into account the bound electron effect of Compton scattering (incoherent scattering), the coherent scattering, the pair production, and the secondary sources of bremsstrahlung and fluorescence.

The tests have shown that the approximating formula reproduces the reference data of double-layer shields very well for most cases. With the same parameters and with a new physical consideration that takes into account in a global way the degradation of the gamma-ray energy spectrum, the buildup factors of three- and five-layer shields were also very well reproduced.