Hamiltonian magnetic reconnection with parallel electron heat flux dynamics

We analyse, both analytically and numerically, a two-dimensional six-field fluid model for collisionless magnetic reconnection, accounting for temperature and heat flux fluctuations along the direction of the magnetic guide field. We show that the model possesses a Hamiltonian structure with a non-canonical Poisson bracket. This bracket is characterized by the presence of six infinite families of Casimirs, associated with Lagrangian invariants. This reveals that the model can be reformulated as a system of advection equations, thus generalizing previous results obtained for Hamiltonian isothermal fluid models for reconnection. Numerical simulations indicate that the presence of heat flux and temperature fluctuations yields slightly larger growth rates and similar saturated island amplitudes, with respect to the isothermal models. For values of the sonic Larmor radius much smaller than the electron skin depth, heat flux fluctuations tend to be suppressed and temperature fluctuations follow density fluctuations. Increasing the sonic Larmor radius results in an increasing fraction of magnetic energy converted into heat flux, at the expense of temperature fluctuations. In particular, heat flux fluctuations tend to become relevant along the magnetic island separatrices. The qualitative structures associated with the electron field variables are also reinterpreted in terms of the rotation of the Lagrangian invariants of the system.