Quasi-Monte Carlo Radiative Transfer

S. G. Shulman

2024-06-22

Abstract:We consider an alternative to the Monte Carlo method for dust continuous radiative transfer simulations: the Quasi-Monte Carlo method. We briefly discuss what it is, its history, and possible implementations. We compare the Monte Carlo method with four pseudo-random number generators and five Quasi-Monte Carlo implementations using different low-discrepancy sequences and the Hammersley set. For the comparison, we study different test matter geometries and problems. We present comparison results for single scatterings of radiation from a point source, multiple scatterings of radiation from a point source, and single scatterings of radiation from a spherical star. In all cases, Quasi-Monte Carlo shows better convergence than Monte Carlo. In several test cases, the gain in computation time to achieve a fixed error value reached 40 times. We obtained ten times speed up in many of the considered tests.

Solar and Stellar Astrophysics,Astrophysics of Galaxies,High Energy Astrophysical Phenomena,Instrumentation and Methods for Astrophysics

What problem does this paper attempt to address?

### Problems Addressed by the Paper

This paper aims to explore the advantages of the **Quasi-Monte Carlo (QMC)** method over the traditional **Monte Carlo (MC)** method in the simulation of dust continuous radiation transfer. Specifically, the authors compare the performance of the two methods in different test scenarios and demonstrate the significant advantage of the QMC method in terms of convergence speed.

#### Main Research Content:
1. **Method Comparison**: The authors conduct a detailed comparison of single scattering, multiple scattering, and scattering from a spherical light source using various pseudorandom number generators and quasi-random number generators.
2. **Error Estimation**: From a mathematical perspective, the paper analyzes the error convergence behavior of the Monte Carlo method, the Quasi-Monte Carlo method, and the Riemann summation method in different dimensions.
3. **Experimental Results**: Through multiple test cases, the paper verifies the significant improvement in computational efficiency of the QMC method, particularly in some cases where the computation time is reduced by more than 40 times.

#### Research Background:
- **Complexity of Radiation Transfer Problems**: Radiation transfer problems typically involve complex 3D integrals, especially when considering multiple scattering.
- **Limitations of the Monte Carlo Method**: Although the Monte Carlo method is widely used in practical applications, its convergence speed is relatively slow in high-dimensional spaces.
- **Advantages of the Quasi-Monte Carlo Method**: By using low-discrepancy sequences (such as Sobol sequences), the QMC method can converge to the exact solution more quickly.

#### Conclusions and Outlook:
- **QMC Method Superior to MC Method**: In most test cases, the QMC method exhibits better convergence speed.
- **Multiple Scattering Problems**: For multiple scattering problems, although the QMC method remains effective, its advantage diminishes as the number of scatterings increases.
- **Recommendation to Use QMC Method**: For single scattering and radiation transfer problems under specific geometric structures, the QMC method is a better choice.