Axisymmetric global gravitational equilibrium for magnetized, rotating hot plasma

> We present analytic solutions for three-dimensional magnetized axisymmetric equilibria confining rotating hot plasma in a gravitational field. Our up-down symmetric solution to the full Grad-Shafranov equation can exhibit equatorial plane localization of the plasma density and current, resulting in disk equilibria for the plasma density. For very weak magnetic fields and high plasma pressure, we find strongly rotating thin plasma disk gravitational equilibria that satisfy strict Keplerian motion provided the gravitational energy is much larger than the plasma pressure, which must be large compared to the magnetic energy of the poloidal magnetic field. When the rotational energy exceeds the gravitational energy and it is larger than the plasma pressure, diffuse disk equilibrium solutions continue to exist provided the poloidal magnetic energy remains small. For stronger magnetic fields and lower plasma pressure and rotation, we can also find gravitational equilibria with strong localization to the equatorial plane. However, a toroidal magnetic field is almost always necessary to numerically verify these equilibria are valid solutions in the presence of gravity for the cases considered in Catto & Krasheninnikov (J. Plasma Phys., vol. 81, 2015, 105810301). In all cases both analytic and numerical results are presented.