Coupling particle-in-cell and magnetohydrodynamics methods for realistic solar flare models

Context. Simulating solar flares requires capturing both large-scale magnetohydrodynamic (MHD) evolution and small-scale kinetic processes near reconnection sites. Bridging these scales has been a significant computational challenge. Aims. This study introduces a Particle-In-Cell (PIC) solver integrated within the DISPATCH framework, facilitating seamless embedding within MHD simulations. This development aims to enable self-consistent multi-scale solar flare simulations. Methods. Our PIC solver, inspired by the PhotonPlasma code, addresses the Vlasov–Maxwell equations in collisionless plasma. We validate its accuracy through fundamental plasma tests – including plasma oscillations, two-stream instability, and current sheet reconnection. To make kinetic simulations computationally feasible, we employ physical adjustment of constants (PAC), modifying the speed of light, elementary charge, and electron mass to shift plasma scales. Additionally, we implement and validate a coupling strategy that enables smooth transitions between kinetic and fluid regimes. Results. The PIC solver successfully recovers expected plasma dynamics and electromagnetic field behaviour. Our analysis highlights the effects of PAC on reconnection dynamics, underscoring the importance of transparent and well-documented scaling choices. Test cases involving propagating waves across PIC-MHD interfaces confirm the robustness of our coupling approach. Conclusions. The integration of the PIC solver into the DISPATCH framework makes it possible to run self-consistent, multi-scale solar flare simulations. Our approach provides a computationally efficient foundation for investigating reconnection physics in large-scale astrophysical plasmas.