No Arabic abstract
The difference between the guiding center phase-space Lagrangians derived in [J.W. Burby, J. Squire, and H. Qin, Phys. Plasmas {bf 20}, 072105 (2013)] and [F.I. Parra, and I. Calvo, Plasma Phys. Control. Fusion {bf 53}, 045001 (2011)] is due to a different definition of the guiding center coordinates. In this brief communication the difference between the guiding center coordinates is calculated explicitly.
This paper had no abstract originally. A second-order symplectic integration algorithm for guiding center motion is presented. The algorithm is based on the Poincare (mid-point) generating function.
Upon combining Northrops picture of charged particle motion with modern liquid crystal theories, this paper provides a new description of guiding center dynamics (to lowest order). This new perspective is based on a rotation gauge field (gyrogauge) that encodes rotations around the magnetic field. In liquid crystal theory, an analogue rotation field is used to encode the rotational state of rod-like molecules. Instead of resorting to sophisticated tools (e.g. Hamiltonian perturbation theory and Lie series expansions) that still remain essential in higher-order gyrokinetics, the present approach combines the WKB method with a simple kinematical ansatz, which is then replaced into the charged particle Lagrangian. The latter is eventually averaged over the gyrophase to produce Littlejohns guiding-center equations. A crucial role is played by the vector potential for the gyrogauge field. A similar vector potential is related to liquid crystal defects and is known as `wryness tensor in Eringens micropolar theory.
Electromagnetically induced transparency (EIT) cooling has established itself as one of the most widely used cooling schemes for trapped ions during the past twenty years. Compared to its alternatives, EIT cooling possesses important advantages such as a tunable effective linewidth, a very low steady state phonon occupation, and applicability for multiple ions. However, existing analytic expression for the steady state phonon occupation of EIT cooling is limited to the zeroth order of the Lamb-Dicke parameter. Here we extend such calculations and present the explicit expression to the second order of the Lamb-Dicke parameter. We discuss several implications of our refined formula and are able to resolve certain difficulties in existing results.
We present a program for the reduction of large systems of integrals to master integrals. The algorithm was first proposed by Laporta; in this paper, we implement it in MAPLE. We also develop two new features which keep the size of intermediate expressions relatively small throughout the calculation. The program requires modest input information from the user and can be used for generic calculations in perturbation theory.
First-order accurate degenerate variational integration (DVI) was introduced in C. L. Ellison et. al, Phys. Plasmas 25, 052502 (2018) for systems with a degenerate Lagrangian, i.e. one in which the velocity-space Hessian is singular. In this paper we introducing second order accurate DVI schemes, both with and without non-uniform time stepping. We show that it is not in general possible to construct a second order scheme with a preserved two-form by composing a first order scheme with its adjoint, and discuss the conditions under which such a composition is possible. We build two classes of second order accurate DVI schemes. We test these second order schemes numerically on two systems having noncanonical variables, namely the magnetic field line and guiding center systems. Variational integration for Hamiltonian systems with nonuniform time steps, in terms of an extended phase space Hamiltonian, is generalized to noncanonical variables. It is shown that preservation of proper degeneracy leads to single-step methods without parasitic modes, i.e. to non-uniform time step DVIs. This extension applies to second order accurate as well as first order schemes, and can be applied to adapt the time stepping to an error estimate.