  The radiative transfer equation (RTE) is a cornerstone for describing the
propagation of electromagnetic radiation in a medium, with applications
spanning atmospheric science, astrophysics, remote sensing, and biomedical
optics. Despite its importance, an exact analytical solution to the RTE has
remained elusive, necessitating the use of numerical approximations such as
Monte Carlo, discrete ordinate, and spherical harmonics methods.
  In this paper, we present an exact solution to the one-dimensional
time-dependent RTE. We delve into the moments of the photon distribution,
providing a clear view of the transition to the diffusion regime. This analysis
offers a deeper understanding of light propagation in the medium.
  Furthermore, we demonstrate that the one-dimensional RTE is equivalent to the
Klein-Gordon equation with an imaginary mass term determined by the inverse
reduced scattering length. Contrary to naive expectations of superluminal
solutions, we find that our solution is strictly causal under appropriate
boundary conditions, determined by the light transport problem.
  We validate the found solution using Monte Carlo simulations and benchmark
the performance of the latter. Our analysis reveals that even for highly
forward scattering, dozens of random light scatterings are required for an
accurate estimate, underscoring the complexity of the problem.
  Moreover, we propose a method for faster convergence by adjusting the
parameters of Monte Carlo sampling. We show that a Monte Carlo method sampling
photon scatterings with input parameters $(\mu_s,g)$, where $\mu_s$ is the
inverse scattering length and $g$ is the scattering anisotropy parameter, is
equivalent to that with $(\mu_s(1-g)/2,-1)$. This equivalence leads to a
significantly faster convergence to the exact solution, offering a substantial
improvement of the Monte Carlo method for the one-dimensional RTE.
