Monte Carlo Iteration

After setting up the simulation, TARDIS runs the simulation using the .run_convergence() method. This runs several Monte Carlo iterations (which will be described in the links below), corresponding to one less than the number of iterations specified in the Monte Carlo Configuration. As will be described in Updating Plasma and Convergence and Spectrum Generation, these iterations will eventually be used to calculate the steady-state plasma properties. TARDIS runs the last iteration of the simulation using the .run_final() method to determine the spectrum.

The following pages provide a very basic introduction to Monte Carlo radiative transfer techniques as they are used in TARDIS. All the information listed here can also be found in various papers by L. Lucy and in the main TARDIS publication (cf. [AbbottLucy85], [MazzaliLucy93], [Lucy99a], [LongKnigge02], [Lucy02], [Lucy03], [Lucy05], [KerzendorfSim14]).

• Monte Carlo Radiative Transfer - Basic Principles
  • Random Sampling Basics
• Energy Packet Initialization
  • Theory
  • Code Example
  • Custom Packet Source
• Packet Propagation
  • Propagation in a Spherical Domain
  • Supernova Expansion
  • Reference Frames
  • Numerical and Physical Events
  • Implementation: Main Propagation Loop
• Line Interaction Treatments
  • Resonant Scattering
  • Downbranching
  • Macro Atom Scheme
  • Comparison
• Volume-Based Estimators
  • Theory
  • Implementation