Observational and theoretical evidence suggests that beams of accelerated particles are produced in flaring events of all sizes in the solar atmosphere, from X-class flares to nanoflares. Current models of such particles in flaring loops assume an isolated 1D atmosphere. A more realistic environment for modelling accelerated particles can be provided by 3D radiative magnetohydrodynamics codes. Here we present a simple model for particle acceleration and propagation in the context of a 3D simulation of the quiet solar atmosphere spanning from convection zone to corona. We then examine the additional transport of energy introduced by the particle beams. Locations of particle acceleration associated with magnetic reconnection are identified by detecting changes in magnetic topology. At each location, the parameters of the accelerated particle distribution are estimated from local conditions. The particle distributions are then propagated along the magnetic field, and the energy deposition due to Coulomb collisions with the ambient plasma is computed. We find that particle beams originate in extended acceleration regions distributed across the corona. Upon reaching the transition region, they converge and produce strands of intense heating penetrating the chromosphere. Within these strands, beam heating consistently dominates conductive heating below the bottom of the transition region. This indicates that particle beams will qualitatively alter the energy transport even outside of active regions.

Context. We introduce a model for including accelerated particles in pure magnetohydrodynamics (MHD) simulations of the solar atmosphere. Aims. We show that the method is viable and produces results that enhance the realism of MHD simulations of the solar atmosphere. Methods. The acceleration of high-energy electrons in solar flares is an accepted fact, but is not included in the most advanced 3D simulations of the solar atmosphere. The effect of the acceleration is not known, and here we introduce a simple method to account for the ability of the accelerated electrons to move energy from the reconnection sites and into the dense transition zone and chromosphere. Results. The method was only run for a short time and with low reconnection energies, but this showed that the reconnection process itself changes, and that there is a clear effect on the observables at the impact sites of the accelerated electrons. Further work will investigate the effect on the reconnection sites and the impact sites in detail.