Dose calculations in proton therapy require the fast and accurate solution of a high-dimensional transport equation for a large number of (pencil) beams with different energies and directions. Deterministically solving this transport problem at a sufficient resolution can however be prohibitively expensive, especially due to highly forward peaked scattering of the protons. We propose using a model order reduction approach, the dynamical low-rank approximation (DLRA), which evolves the solution on the manifold of low-rank matrices in (pseudo-)time. For this, we compare a collided-uncollided split of the linear Boltzmann equation and its Fokker-Planck approximation. We treat the uncollided part using a ray-tracer and combine high-order phase space discretizations and a mixture model for materials with DLRA for the collided equation. Our method reproduces the results of a full-rank reference code at significantly lower rank, and thus computational cost and memory, and further makes computations feasible at much higher resolutions. At higher resolutions, we also achieve good accuracy with respect to TOPAS MC in homogeneous as well as heterogeneous materials. Finally, we demonstrate that several beam sources with different angles can be computed with little cost increase compared to individual beams. ... Read more

Deterministically solving charged particle transport problems at a sufficient spatial and angular resolution is often prohibitively expensive, especially due to their highly forward peaked scattering. We propose a model order reduction approach which evolves the solution on a low-rank manifold in time, making computations feasible at much higher resolutions and reducing the overall run-time and memory footprint. For this, we use a hybrid dynamical low-rank approach based on a collided-uncollided split, i.e., the transport equation is split through a collision source method. Uncollided particles are described using a ray tracer, facilitating the inclusion of boundary conditions and straggling, whereas collided particles are represented using a moment method combined with the dynamical low-rank approximation. Here the energy is treated as a pseudo-time and a rank adaptive integrator is chosen to dynamically adapt the rank in energy. We can reproduce the results of a full-rank reference code at a much lower rank and thus computational cost and memory usage. The solution further achieves comparable accuracy with respect to TOPAS MC as previous deterministic approaches. ... Read more