Vlasov tokamak equilibria with sheared toroidal flow and anisotropic pressure

By choosing appropriate deformed Maxwellian ion and electron distribution functions depending on the two particle constants of motion, i.e., the energy and toroidal angular momentum, we reduce the Vlasov axisymmetric equilibrium problem for quasineutral plasmas to a transcendental Grad-Shafranov-like equation. This equation is then solved numerically under the Dirichlet boundary condition for an analytically prescribed boundary possessing a lower X-point to construct tokamak equilibria with toroidal sheared ion flow and anisotropic pressure. Depending on the deformation of the distribution functions, these steady states can have toroidal current densities either peaked on the magnetic axis or hollow. These two kinds of equilibria may be regarded as a bifurcation in connection with symmetry properties of the distribution functions on the magnetic axis