Abstract. We consider a simple electromagnetic gyrokinetic model for collisionless plasmas and show that it possesses a Hamiltonian structure. Subsequently, from this model we derive a two-moment gyrofluid model by means of a procedure which guarantees that the resulting gyrofluid model is also Hamiltonian. The first step in the derivation consists of imposing a generic fluid closure in the Poisson bracket of the gyrokinetic model, after expressing such bracket in terms of the gyrofluid moments. The constraint of the Jacobi identity, which every Poisson bracket has to satisfy, selects then what closures can lead to a Hamiltonian gyrofluid system. For the case at hand, it turns out that the only closures (not involving integro/differential operators or an explicit dependence on the spatial coordinates) that lead to a valid Poisson bracket are those for which the second order parallel moment, independently for each species, is proportional to the zero order moment. In particular, if one chooses an isothermal closure based on the equilibrium temperatures and derives accordingly the Hamiltonian of the system from the Hamiltonian of the parent gyrokinetic model, one recovers a known Hamiltonian gyrofluid model for collisionless reconnection. The proposed procedure, in addition to yield a gyrofluid model which automatically conserves the total energy, provides also, through the resulting Poisson bracket, a way to derive further conservation laws of the gyrofluid model, associated with the so called Casimir invariants. We show that a relation exists between Casimir invariants of the gyrofluid model and those of the gyrokinetic parent model. The application of such Hamiltonian derivation procedure to this two-moment gyrofluid model is a first step toward its application to more realistic, higher-order fluid or gyrofluid models for tokamaks. It also extends to the electromagnetic gyrokinetic case, recent applications of the same procedure to Vlasov and driftkinetic systems.
IntroductionFluid models represent a widespread and effective tool for investigating important phenomena in fusion plasmas such as instabilities, turbulence and reconnection events. Indeed, fluid models offer a considerable advantage, in terms of required computational resources, with respect to kinetic models. On the other hand, compared to these, they obviously suffer from limitations in the range of scales and frequencies of the phenomena that they can describe. The derivation of sophisticated fluid models aiming at partially remedying such limitations is an active line of research since a long time. An essential part of the problem is related to the closure adopted to truncate the fluid hierarchy of equations obtained by taking moments of a parent kinetic model. A considerable effort, directed also toward applications to space plasma turbulence, has been carried out to derive fluid and gyrofluid models (see, e.g. Refs. [1][2][3][4][5][6][7][8]) by taking moments of