We present a novel Dynamic Monte Carlo code (DynaMo code) which self-consistently simulates the Compton scattering driven dynamic evolution of a plasma. We use the DynaMo code to investigate the time-dependent expansion and acceleration of dissipationless GRB fireballs by varying their initial opacities and baryonic content. We study the opacity and energy density evolution of an initially optically thick, radiation-dominated fireball across its entire phase space - in particular during the $R_{\rm ph} < R_{\rm sat}$ regime. Our results reveal new phases of fireball evolution: a transition phase with a radial extent of several orders of magnitude - the fireball transitions from $\Gamma \propto R$ to $\Gamma \propto R^0$, a post-photopsheric acceleration phase - where fireballs accelerate beyond the photosphere, and a Thomson-dominated acceleration phase - characterized by slow acceleration of optically thick, matter-dominated fireballs due to Thomson scattering. We quantify the new phases by providing analytical expressions of Lorentz factor evolution, which will be useful for deriving jet parameters.

It has been proposed that blazar jets are structured, with a fast spine surrounded by a slower sheath or layer. This structured jet model explains some properties of their emission and morphology. Because of their relative motion, the radiation produced by one component is seen amplified by the other, thus enhancing the inverse Compton emission of both. Radiation is emitted anisotropically in the comoving frames, and causes the emitting plasma to recoil. As seen in the observer frame, this corresponds to a deceleration of the fastest component (the spine) and an acceleration of the slower one (the layer). While the deceleration of the spine has already been investigated, here we study for the first time the acceleration of the sheath and find self-consistent velocity profile solutions for both the spine and the sheath while accounting for radiative cooling. We find that the sheath can be accelerated to the velocities required by the observations if its leptons remain energetic in the acceleration region, assumed to be of the order of 100 Schwarzschild radii, demanding continuous injection of energetic particles in that region.

One of the most important unresolved issues in gamma-ray burst physics is the origin of the prompt gamma-ray spectrum. Its general non-thermal character and the softness in the X-ray band remain unexplained. We tackle these issues by performing Monte Carlo simulations of radiation-matter interactions in a scattering dominated photon-lepton plasma. The plasma -- initially in equilibrium -- is driven to non-equilibrium conditions by a sudden energy injection in the lepton population, mimicking the effect of a shock wave or the dissipation of magnetic energy. Equilibrium restoration occurs due to energy exchange between the photons and leptons. While the initial and final equilibrium spectra are thermal, the transitional photon spectra are characterized by non-thermal features such as power-law tails, high energy bumps, and multiple components. Such non-thermal features are observed at infinity if the dissipation occurs at small to moderate optical depths, and the spectrum is released before thermalization is complete. We model the synthetic spectra with a Band function and show that the resulting spectral parameters are similar to observations for a frequency range of 2-3 orders of magnitude around the peak. In addition, our model predicts correlations between the low-frequency photon index and the peak frequency as well as between the low- and high-frequency indices. We explore baryon and pair dominated fireballs and reach the conclusion that baryonic fireballs are a better model for explaining the observed features of gamma-ray burst spectra.