Dust Scattering in Protoplanetary Disks

Using the Monte Carlo radiative transfer method, the following code, written in C++ and building on the Athena++ framework, adds a program for dust scattering simulations. We then compare the output of the code to a tested dust simulation code written in FORTRAN. For this comparison we use a problem generator which simulates a protoplanetary disc of varying density which surrounds a binary star system with the stars having different luminosities. We then discuss the issues with our code and how we might correct the code and improve on it.