Gyrofluid simulations of collisionless reconnection in the presence of diamagnetic effects

The effects of the ion Larmor radius on magnetic reconnection are investigated by means of numerical simulations, with a Hamiltonian gyrofluid model. In the linear regime, it is found that ion diamagnetic effects decrease the growth rate of the dominant mode. Increasing ion temperature tends to make the magnetic islands propagate in the ion diamagnetic drift direction. In the nonlinear regime, diamagnetic effects reduce the final width of the island. Unlike the electron density, the guiding center density does not tend to distribute along separatrices and at high ion temperature, the electrostatic potential exhibits the superposition of a small scale structure, related to the electron density, and a large scale structure, related to the ion guiding-center density

Gyrofluid simulations of collisionless reconnection in the presence of diamagnetic effects

The effects of the ion Larmor radius on magnetic reconnection are investigated by means of numerical simulations, with a Hamiltonian gyrofluid model. In the linear regime, it is found that ion diamagnetic effects decrease the growth rate of the dominant mode. Increasing ion temperature tends to make the magnetic islands propagate in the ion diamagnetic drift direction. In the nonlinear regime, diamagnetic effects reduce the final width of the island. Unlike the electron density, the guiding center density does not tend to distribute along separatrices and at high ion temperature, the electrostatic potential exhibits the superposition of a small scale structure, related to the electron density, and a large scale structure, related to the ion guiding-center density

Gyrofluid simulations of collisionless reconnection in the presence of diamagnetic effects

The effects of the ion Larmor radius on magnetic reconnection are investigated by means of numerical simulations, with a Hamiltonian gyrofluid model. In the linear regime, it is found that ion diamagnetic effects decrease the growth rate of the dominant mode. Increasing ion temperature tends to make the magnetic islands propagate in the ion diamagnetic drift direction. In the nonlinear regime, diamagnetic effects reduce the final width of the island. Unlike the electron density, the guiding center density does not tend to distribute along separatrices and at high ion temperature, the electrostatic potential exhibits the superposition of a small scale structure, related to the electron density, and a large scale structure, related to the ion guiding-center density

Electronic structure and carrier transfer in B-DNA monomer polymers and dimer polymers: Stationary and time-dependent aspects of wire model vs. extended ladder model

We employ two Tight-Binding (TB) approaches to study the electronic structure and hole or electron transfer in B-DNA monomer polymers and dimer polymers made up of $N$ monomers (base pairs): (I) at the base-pair level, using the on-site energies of base pairs and the hopping integrals between successive base pairs, i.e., a wire model and (II) at the single-base level, using the on-site energies of the bases and the hopping integrals between neighboring bases, i.e., an \textit{extended} ladder model since we also include diagonal hoppings. We solve a system of $MD$ ("matrix dimension") coupled equations [(I) $MD$ = $N$, (II) $MD$ = $2N$] for the time-independent problem, and a system of $MD$ coupled $1^\text{st}$ order differential equations for the time-dependent problem. We study the HOMO and the LUMO eigenspectra, the occupation probabilities, the Density of States (DOS) and the HOMO-LUMO gap as well as the mean over time probabilities to find the carrier at each site [(I) base pair or (II) base)], the Fourier spectra, which reflect the frequency content of charge transfer (CT) and the pure mean transfer rates from a certain site to another. The two TB approaches give coherent, complementary aspects of electronic properties and charge transfer in B-DNA monomer polymers and dimer polymers.Comment: 20 pages, 23 figure

Mode signature and stability for a Hamiltonian model of electron temperature gradient turbulence

Stability properties and mode signature for equilibria of a model of electron temperature gradient (ETG) driven turbulence are investigated by Hamiltonian techniques. After deriving the infinite families of Casimir invariants, associated with the noncanonical Poisson bracket of the model, a sufficient condition for stability is obtained by means of the Energy-Casimir method. Mode signature is then investigated for linear motions about homogeneous equilibria. Depending on the sign of the equilibrium "translated" pressure gradient, stable equilibria can either be energy stable, i.e.\ possess definite linearized perturbation energy (Hamiltonian), or spectrally stable with the existence of negative energy modes (NEMs). The ETG instability is then shown to arise through a Kre\u{\i}n-type bifurcation, due to the merging of a positive and a negative energy mode, corresponding to two modified drift waves admitted by the system. The Hamiltonian of the linearized system is then explicitly transformed into normal form, which unambiguously defines mode signature. In particular, the fast mode turns out to always be a positive energy mode (PEM), whereas the energy of the slow mode can have either positive or negative sign

On the rate of convergence of the Hamiltonian particle-mesh method

The Hamiltonian Particle-Mesh (HPM) method is a particle-in-cell method for compressible fluid flow with Hamiltonian structure. We present a numer- ical short-time study of the rate of convergence of HPM in terms of its three main governing parameters. We find that the rate of convergence is much better than the best available theoretical estimates. Our results indicate that HPM performs best when the number of particles is on the order of the number of grid cells, the HPM global smoothing kernel has fast decay in Fourier space, and the HPM local interpolation kernel is a cubic spline

Noncollisional plasmoid instability based on a gyrofluid and gyrokinetic integrated approach

In this work, the development of two-dimensional current sheets with respect to tearing-modes, in collisionless plasmas with a strong guide field, is analysed. During their non-linear evolution, these thin current sheets can become unstable to the formation of plasmoids, which allows the magnetic reconnection process to reach high reconnection rates. We carry out a detailed study of the impact of a finite $\beta_e$, which also implies finite electron Larmor radius effects, on the collisionless plasmoid instability. This study is conducted through a comparison of gyrofluid and gyrokinetic simulations. The comparison shows in general a good capability of the gyrofluid models in predicting the plasmoid instability observed with gyrokinetic simulations. We show that the effects of $\beta_e$ promotes the plasmoid growth. The impact of the closure applied during the derivation of the gyrofluid model is also studied through the comparison of the energy variation

Gyrofluid simulations of collisionless reconnection in the presence of diamagnetic effects

International audienceThe effects of the ion Larmor radius on magnetic reconnection are investigated by means of numerical simulations, with a Hamiltonian gyrofluid model. In the linear regime, it is found that ion diamagnetic effects decrease the growth rate of the dominant mode. Increasing ion temperature tends to make the magnetic islands propagate in the ion diamagnetic drift direction. In the nonlinear regime, diamagnetic effects reduce the final width of the island. Unlike the electron density, the guiding center density does not tend to distribute along separatrices and at high ion temperature, the electrostatic potential exhibits the superposition of a small scale structure, related to the electron density, and a large scale structure, related to the ion guiding-center density