Rotation matrix of a charged symmetrical body: one-parameter family of solutions in elementary functions

Equations of motion of a charged symmetrical body in external constant and homogeneous electric and magnetic fields are deduced starting from the variational problem, where the body is considered as a system of charged point particles subject to holonomic constraints. The final equations are written for the center-of mass-coordinate, rotation matrix and angular velocity. General solution to the equations of motion is obtained for the case of a charged ball. For the case of a symmetrical charged body (solenoid), the task of obtaining the general solution is reduced to the problem of a one-dimensional cubic pseudo-oscillator. Besides, we present a one-parametric family of solutions to the problem in elementary functions.

Euler-Poisson equations of a charged spinning body in external constant and homogeneous electric and magnetic fields are deduced starting from the variational problem, where the body is considered as a system of charged point particles subject to holonomic constraints.The final equations are written for the center-of mass-coordinate, rotation matrix and angular velocity.General solution to the equations of motion is obtained for the case of a charged ball.For the case of a symmetrical charged body (solenoid), the task of obtaining the general solution to Euler equations is reduced to the problem of a one-dimensional cubic pseudo-oscillator.Besides, we present a one-parametric family of solutions to the problem in elementary functions.

I. INTRODUCTION.
Spinning bodies represent an important object of study in modern research at different scales, from problems of approaching compact stars in gravitational wave physics to levitating nano-particles and elementary particles with spin.Behavior of a spinning object in special and general relativity as well as in quantum mechanics in many cases is considered in first order by perturbation theory based on exact solutions to the classical problem [1].Therefore, the search for new integrable cases and analytical solutions in the dynamics of a spinning body is important for further progress in such problems [2][3][4].In the present work we consider a charged spinning body in external constant and homogeneous electric and magnetic fields.While the case of a ferromagnet has been discussed quite widely in the literature [5][6][7][8][9][10], much less attention has been paid to a charged dielectric [11].In this article we will try to fill this gap in the literature.
The work is organized as follows.In Sect.II we deduce equations of motion on the base of a Lagrangian action, formulated for the case under consideration.We will present a detailed derivation of the equations, since semiempirical methods applied to a spinning body in some cases lead to either inaccuracies or erroneous interpretation of the final result [12][13][14][15].In Sect.III we present general solution to obtained equations for the case of a charged ball.In Sect.IV we reduce the problem of a symmetrical charged body to the problem of a one-dimensional non-linear pseudo-oscillator.In Sect.V we present a one-parametric family of solutions in elementary functions for the motions with specially chosen initial angular velocity of the symmetrical charged body.In the Appendix, for the convenience of the reader, we summarized the motion of a point charged particle in external electromagnetic field.

II. CHARGED BODY IN CONSTANT AND HOMOGENEOUS ELECTRIC AND MAGNETIC FIELDS.
Consider a rigid body that consist of n particles y N (t) of charge e N and mass m N , N = 1, 2, . . ., n.Its Lagrangian action reads [12] The first term is kinetic energy of all particles, while the remaining terms account the presence of constraints, that guarantee that distances and angles among the particles do not change with time 1 .The constraints were added with help of Lagrangian multipliers λ AN (t).In all calculations these auxiliary variables should be treated on equal footing with y N (t).In particular, looking for the equations of motion, we take variations with respect to y N and all λ AN .The 3 × 3 -block λ AB of λ AN was chosen to be the symmetric matrix.The variations with respect to λ AN imply the constraints, which therefore arise as a part of conditions of extreme of the action functional.So the presence of λ AN allows y i N to be treated as unconstrained variables, that should be varied independently in obtaining the equations of motion.
We consider the body immersed into constant and homogeneous electric and magnetic fields with scalar potential A 0 (y N ) and vector potential A(y N ), see Appendix.Generally, its movement will produce time-dependent distribution of charges and currents leading to radiation of an electromagnetic fields, see Chapter 17 in [17].We neglect the resulting radiation and radiation damping effects.Then, summing up the potential energies (103) of body's particles, we obtain its total potential energy Adding it to the action (1), we obtain a variational problem for the body in external electric and magnetic fields.We assume that all particles of the body have the same charge to mass ratio, e N m N = µ for any N = 1, 2, . . ., n.Then our action implies the following dynamical equations Introducing the center of mass, that is center-of-mass behaves like charged point particle discussed in Appendix.In particular, rotational motion of the body does not affect its translational motion.Substituting y N = y 0 +x N into Eq.( 3) and taking into account (4), we rewrite these equations in the center-of-mass coordinate system Each solution x i N (t) to these equations is of the form , where R ij (t) is an orthogonal matrix that, by construction, obeys the universal initial data R ij (0) = δ ij .Substituting this expression into Eqs.(5), then multiplying the equation with number N by x j N (0) and taking their sum, we obtain the following second-order equations for determining the rotation matrix R ik : It was denoted Besides, g ik is the mass matrix while λ jk (t) is the symmetric matrix where all x i N are taken at the instant t = 0. Due to the identity m N x N = 0, satisfied for the center-of-mass coordinates, second term on r. h.s. of Eq. ( 6) vanishes.In the result, the electric field do not affect the motion of rotational degrees of freedom.Besides, the center-of-mass variable do not enter into this equation, so the translational motion does not affect the rotational motion of the body.
The variables λ jk (t) in equations ( 6) depend on the unknown dynamical variables λ AN (t).Fortunately, we do not need to know λ jk (λ AN (t)), because of these equations determine λ jk algebraically, as some functions of R and Ṙ.This result, obtained with use of the procedure described in [12], can be formulated as follows Affirmation.Consider the second-order system for determining the variables R ij (t) and λ kj (t) where A ik (R ab , Ṙab ) is some given matrix that does not depend on g kj and λ kj (as before, g kj is a numerical symmetric non-degenerate matrix, and λ kj (t) = λ jk (t)).The problem (10) is equivalent to the following Cauchy problem for the first-order system, written for the mutually independent variables R ij (t) and Ω i (t): where I is inertia tensor of the body with the components Note that this system is composed of SO(3) vectors and tensors, so it is covariant under the rotations.We will work with these equations assuming that the mass matrix and inertia tensor are of diagonal form.This implies [12], that at initial instant t = 0 the Laboratory basis vectors e i were taken in the directions of axes of inertia R i (t) taken as the body-fixed frame: e i = R i (0).We also recall that the body-fixed basis vectors are columns of the rotation matrix: R = (R 1 , R 2 , R 3 ).Eigenvalues of mass matrix and inertia tensor are related as follows: 2g 1 = I 2 + I 3 − I 1 , and so on.
Our equations ( 6) are of the form (10), so they are equivalent to the first-order system (11) and (12) with A written in (7).To obtain an explicit form of Eq. ( 11), we use (12) to rewrite the quantity A ia as follows: Using the latter expression for A ia in Eq. (11), after direct calculations this acquires the form The vector composed of last three terms can be written in a more compact form in terms of the mass matrix.Indeed, writing the first component of this vector in explicit form we get and similar expressions for the second and third components.Then Eq. ( 11) acquires the final form For completeness we also present equations for the vector of angular momentum m = RIΩ and for its components M = IΩ in the body-fixed frame Let us consider Eqs. ( 16) in the Laboratory system with third axis in the direction of magnetic vector B. Denote by I ′ and g ′ the non diagonal tensors of inertia and mass that will appear in this coordinates.Then Eqs. ( 16) acquire the form , where G 3 is the third row of the rotation matrix R ij .They coincide with those deduced by G. Grioli, see [1,11].
In resume, we have succeeded in obtaining equations of motion ( 16), ( 12) and ( 4) for a spinning charged body in external electric and magnetic fields.The solution R ij (t), y i 0 (t) to these equations contains complete information on the evolution of the body with respect to Laboratory frame: dynamics of the body's point y i N (t) with initial position

III. GENERAL SOLUTION TO THE EQUATIONS OF A CHARGED BALL.
Consider a totally symmetric charged body: This could be a charged ball.Then its center moves according to Eq. ( 101).From kinematic relations m between angular momentum m, angular velocity ω and its components Ω i in the body-fixed frame, together with Eq. ( 19), we get m = I 0 ω = I 0 RΩ.The first equality means that for all t the instantaneous rotation axis ω remains parallel with the vector of angular momentum.Substituting (19) into equations of previous section, we get The last equation implies that angular momentum precesses around B with Larmor's frequency α = µ|B| 2c .Besides, the length of angular momentum and its projection on B -axis are first integrals Components Ω i of angular velocity ω with respect to body-fixed frame precess with the same frequency around the components (R T B) i of magnetic field B in body-fixed frame.
Substituting Ω = 1 I0 R T m into (21) and taking into account (20), we arrive at the equation ( 22).So the system (20), ( 21 Further, the first (non linear on R) equation of latter system can be replaced on the linear equation Ṙij = 1 I0 ϵ iab m a R bj .The resulting system an the initial one have the same solutions in the set of orthogonal matrices.In the result, instead of Eqs. ( 20) and ( 21), a totally symmetric body can be described by the equations where all quantities are defined with respect to the Laboratory system.The initial conditions are R ij (0) = δ ij .
According to these equations, m(t) precesses around constant vector B, while the vectors R j (t) instantaneously precess around m(t).
Let us obtain the general solution to the system (24), (25).For a totally symmetric body we can choose the directions of Laboratory axes as convenient, this will not violate the diagonal form of inertia tensor.Using this freedom, we choose the Laboratory system so that at t = 0 the vectors B and m(0) lie in the plane of e 2 , and e 3 , and B is directed along e 3 , see Figure 1.Then B = (0, 0, |B|) T , m(0) = (0, m 2 , m 3 ) T , and is a solution to Eq. ( 25), where the precession frequency is the Larmor's frequency To solve Eq. ( 24) with this m(t), we look for its solution in the form To fix k and γ, we substitute the ansatz (28) into (24) an then take t = 0 in the resulting expressions.They determine k as follows: k2 = − m2 γI0 , k3 = − αI0+m3 γI0 . Then k2 The obtained equalities allow us to represent k and γ through m and α as follows: By direct calculations, it can be verified that the expression (28) with these k and γ satisfies the equations (24).
In resume, we obtained analytical solution for a charged ball launched with initial angular velocity ω = 1 I0 m = 1 I0 (0, m 2 , m 3 ) T in constant and homogeneous electric E and magnetic B fields.It is given by the double-frequency rotation matrix (28), (29).The total motion can be thought as a superposition of two rotations: the first around unit vector k = (0, k2 , k3 ) T with the frequency γ, and the second around the axis of magnetic field B = (0, 0, |B|) T with the frequency α = µ|B| 2c .The angular momentum vector m(t) precesses around the vector B with the Larmor's frequency α = µ|B| 2c .Let us consider the ball launched with initial vector of angular velocity parallel to the vector of magnetic field B. That is the initial conditions are m(0) = (0, 0, m 3 ) T .Then m(t) = (0, 0, m 3 ) T , k = (0, 0, −1) T and γ = α + m3 I0 .With these values, the rotation matrix (28) reduces to As it should be expected, the ball experiences a stationary rotation around the vector of magnetic field B with the frequency m3 I0 .
IV. SYMMETRICAL CHARGED BODY AND ONE-DIMENSIONAL NON LINEAR PSEUDO-OSCILLATOR.
In this section we start to study the symmetrical charged body.We show that for any solution R ij (t), Ω i (t) to Euler-Poisson equations, the function Ω 3 (t) obeys to the equation of a one-dimensional cubic pseudo-oscillator, see Eq. (47) below.Besides, when Ω 3 (t) is known, the functions Ω 1 (t) and Ω 2 (t) can be found by quadratures, see Eqs. (44) and (45).
Consider Eqs. ( 12) and ( 16) for the symmetrical body2 : I = diagonal(I 2 , I 2 , I 3 ).This implies the following mass matrix: g = 1 2 diagonal(I 3 , I 3 , 2I 2 − I 3 ).We consider the positively charged body, then the charge-mass ratio is a positive number, µ > 0. We assume that at t = 0 the third inertia axis of the body is vertical.Then, without spoiling the diagonal form of inertia tensor in our equations, the Laboratory system can be chosen as shown on Figure 2. The basis vector e 3 is directed along the third inertia axis R 3 (0), the vectors e 2 and e 3 lie on the plane of paper sheet together with the vector of constant magnetic field B = (0, B 2 > 0, B 3 > 0) T .The initial instantaneous angular velocity of the body is It is convenient to introduce the following notation: Contracting the Poisson equations ( 12) with B ′ i , we get K = −[Ω, K].This equation together with (16) give us the auxiliary system of 3 + 3 closed equations for determining the variables Ω(t) and K(t) By construction, the initial conditions for K(t) are K(0) = B ′ .Any solution R ij (t), Ω i (t) to Euler-Poisson equations obeys to this system.So we can use the latter to look for the angular velocity Ω i (t).This system admits four integrals of motion.Two of them are To obtain two more integrals, we write our system in components The equations with Ω3 and K3 imply the third integral Combining the equations with Ω1 , Ω2 , K1 and K2 we get one more integral of motion We written them through the integrations constants c 3 and c 5 , as well as through the initial data ω i and B ′ i of the problem.
Using (40), (41) and the equation with Ω3 of the system (38), we represent the variables K i through Ω i as follows Substituting them into equations for Ω1 and Ω2 from (38), we get where ϕ(t) and f (t) turn out to be the following functions of Ω 3 (t): If Ω 3 (t) is known, the equations (43) can be immediately integrated as follows where Φ(t) is indefinite integral of ϕ(t), while f 0 and ϕ 0 are the integration constants.So it remains to find the third component Ω 3 (t).To this aim we compute the time derivative of the last equation from (38), and use other equations of the system (38), (39) in the resulting expression, presenting it as follows Using the integrals of motion (36), (37) and (41), we obtain closed equation for determining Ω 3 (t), that can be called the equation of cubic pseudo-oscillator where the numeric coefficients a i are functions of initial data of original problem It is not difficult to obtain a two-parametric family of simple solutions to the equation (47).Note that Ω 3 (t) = ω 3 will be (constant) solution to (47) if the third component ω 3 if initial angular velocity is a root of the third degree polynomial on the right side of (47).Substituting Ω 3 (t) = ω 3 into Eq.(47), we obtain the condition on initial data ω 1 , ω 2 , ω 3 under which ω 3 satisfies this equation.To obtain this condition, it is convenient to represent Eq. ( 47) in terms of the initial data, keeping the combinations like ω 3 − Ω 3 (t) as follows: Substituting Ω 3 (t) = ω 3 , we get that (47) will be satisfied only for the initial data ω i obeying the following equation: or, equivalently This is a surface of second order.Since the point with ω i = 0 obeys this equation, the surface always pass through the origin of coordinate system.Resolving (50) with respect to ω 3 , we get the following two-parametric family of constant solutions to the equation (47) of cubic pseudo-oscillator: Let us find out which quadric is defined by the equation ( 51), by writing it in the canonical form.Following the standard procedure [16], we arrive at the new coordinates ω ′ i : where β is the angle between the vectors e 2 and I −1 B ′ , see Figure 3.The new coordinates are obtained from ω i by shifting the origin of coordinate system to the point −B ′ , and subsequent rotation counter-clockwise by the angle In these coordinates Eq. (51) acquires the form Depending on the relationship between the inertia moments I 2 and I 3 , it describes different surfaces.
1. Let 1 < I ≡ I3 I2 < 2. This body could be a charged sufficiently short cylindrical surface.If it rotates around its coaxial axis, it will produce a magnetic field corresponding to a short solenoid.The equation (54) turn into Hence the surface of initial data is a hyperboloid of two sheets.Its upper sheet is shown in Figure 3 (a).In the limiting case I ≡ I3 I2 = 2 we have a plane body, that could be charged circular loop.In this case the sheets of the hyperboloid are tangent to the horizontal planes ω 3 = 0 and ω 3 = −2B ′ 3 .

For the totally symmetric body I ≡ I3
I2 = 1, the equation (54) turn into the cone with semi axes a 2 and a 3 written in equation (54).
3. Let 0 < I ≡ I3 I2 < 1.This body could be a charged long cylindrical surface.If it rotates around its coaxial axis, it will produce a magnetic field corresponding to a long solenoid.The equation (54) turn into Hence the surface of initial data is a hyperboloid of one sheet shown in Figure 3 (b).
In resume, we have shown that for any solution to the Euler-Poisson equations ( 12) and ( 16) of a symmetrical charged body, the function Ω 3 (t) obeys the equation of cubic pseudo-oscillator (47).We obtained a two-parameter family of constant solutions (52) to this equation.Not all of them generate solutions to the original problem.In the next section they will help us to obtain a one-parameter family of solutions to the original Euler-Poisson equations in elementary functions.

V. ROTATION MATRIX: ONE-PARAMETRIC FAMILY OF SOLUTIONS IN ELEMENTARY
FUNCTIONS.
As we saw in previous section, our problem ( 12), ( 16) probably admits solutions with constant Ω 3 (t).So let us search for solutions of the auxiliary task (34), (35) of the form Substituting this ansatz into the equations ( 34) and ( 35), they turn into These equations admit three integrals of motion: Using the equations that lie on the plane ω 1 = 0. Geometrically, these are hyperbolas that appear as a result of the intersection of the hyperboloids in Figure 3 with this plane.
For the circular loop or short solenoid they are For the long solenoid they are At last, for a totally symmetric body they degenerate into the straight lines Resolving Eq. ( 67) with respect to ω 3 we get With this ω 3 , the two systems (63) and (64) depend on the same frequency and imply the following solution This means that vector of angular velocity in the body-fixed frame R i precesses around the third axis R 3 with the frequency ϕ.The next step is to solve the Poisson equations (12).We consider them in the form with Ω(t) specified by Eqs. ( 71)-( 73).Here γ is any one of rows of the rotation matrix.
In components this reads This system admits the integral of motion where B2 ≡ B 2 /|B| and B3 ≡ B 3 /|B| are components of unit vector in the direction of magnetic vector B. Preservation in time of the quantity (76) can be verified by direct computation of its time-derivative, with use the identities We need to find the general solution to the equations (75).Then, according to [12], the rows of the rotation matrix R ij can be obtained taking the following three particular solutions.The row R 11 , R 12 , R 13 is γ(t) with the initial data γ(0) = (1, 0, 0) and with c = Ω 1 (0) = 0.The row R 21 , R 22 , R 23 is γ(t) with the initial data γ(0) = (0, 1, 0) and with c = Ω 2 (0) = ω 2 .At last, the row R 31 , R 32 , R 33 is γ(t) with the initial data γ(0) = (0, 0, 1) and with c = B3 ω 2 / B2 .
First we solve algebraically the equations , representing γ 1 and γ 2 as follows: Substituting them into the equation for γ1 , we get closed equation of second order for γ 3 (t) This is equation of harmonic oscillator with constant frequency k, under the action of an external constant force.Its general solution with the integration constants b and k 0 is Substituting this result into the expressions (78), we obtain the remaining variables At t = 0 we get Solving the equations (82) with the data described below Eq. (77), we get, in each case Substituting these values into Eqs.( 80) and (81) we get the rotation matrix of a symmetrical charged body, immersed into the magnetic field B = (0, B 2 , B Two frequencies in the problem are: ϕ written in Eq. ( 72), and k = ω 2 / B2 .The dependence of the rotation matrix on the inertia moments I 2 , I 3 as well as on the charge-mass ratio µ is hidden in the frequency ϕ.By direct substitution of obtained functions (73) and (84) into the equations ( 12) and ( 16), I verified that they are satisfied.
The rotation matrix can be decomposed as follows: Then the position x(t) of any point of the body at the instant t is: It is obtained by rotating the initial position vector x(0) first around the laboratory axis OZ by the angle ϕt and then around the B -axis by the angle kt.
It can be said that the motion is the composition of a proper rotation around third inertia axis with precession of this axis around the vector of magnetic field B. The final answer (84) admits the limit of totally symmetric body I 2 = I 3 , this implies ϕ = 0.The resulting motion is the precession around the magnetic vector B without a proper rotation.
Combining Eqs. ( 71), ( 72) and (79) we get the relation between two frequencies of the motion (84) We recall that the most general motion of a free symmetrical body is the precession without nutation [12].Observe that the rotation matrix (84) coincides with Eq. ( 132) of this work if we replace B, k B and ϕ B on m, k and ϕ.The physical meaning of this coincidence can be formulated as follows.
Affirmation.If a symmetrical charged body in the magnetic field B = (0, B 2 , B 3 ) moves according (84) with the precession frequency k B around B and the proper rotation frequency ϕ B , then in the absence of a magnetic field its precession with the same frequency k B around the unit vector B will happen with the proper rotation frequency Indeed, consider the motion (84) with initial angular velocity ω(0) = (0, ω 2 , ω 3 (ω 2 )) T .Let it then was launched in the absence of a magnetic field with initial angular velocity Ω(0) = (0, ω 2 , I2ω2B3 I3B2 ) T .According to [12], it will precess around the vector of conserved angular momentum m = IΩ(0) = (0, Components of angular momentum m(t) for our solutions in elementary functions are not conserved quantities.But using the integrals of motion (36), ( 40) and (41) with Ω 3 That is the angular momentum always lies in the plane orthogonal to the constant vector of magnetic field.

VI. CONCLUSION.
In this work we deduced equations of motion of a charged symmetrical body in external constant and homogeneous electric and magnetic fields starting from the variational problem (1) and (2), where the body is considered as a system of charged point particles subject to holonomic constraints.The final equations are written in terms of centerof mass-coordinate, rotation matrix and angular velocity.They are (4), ( 12) and ( 16).According to them, rotational motion of the body does not perturb its translational motion and vice-versa.In particular, the center of mass obeys to Eq. ( 4) and behaves as a point charged particle in the electromagnetic field.Besides, the electric field does not affect the rotational motion of the body.
For the case of a totally symmetrical body (charged ball) we found general solution (28), (29) to the equations of motion.The resulting motion can be thought as a superposition of two rotations: the first around unit vector k written in Eq. ( 29), with the frequency γ determined by initial values of angular velocity and Larmor's frequency; and the second around the axis of magnetic field B = (0, 0, |B|) T with the frequency α = µ|B| 2c .The angular momentum vector m(t) precesses around the vector B with the Larmor's frequency α = µ|B| 2c .Analysing the equations ( 12) and ( 16) for the case of a symmetrical charged top, we demonstrated that the task to find the components of angular velocity Ω i (t) can be reduced to solving the equation of a one-dimensional cubic pseudo-oscillator (47).We found a two-parametric family of solutions (52) to this equation.This helped us later find a one-parametric family of solutions (84) for the rotation matrix of a symmetrical charged body, immersed into the magnetic field B = (0, B 2 , B 3 ) T , and launched with initial angular velocity (73).The resulting motions turn out to be the composition of a proper rotation around third inertia axis with precession of this axis around the vector of magnetic field B.
Consider a particle with mass m and charge e (e < 0 for the electron), moving subject to constant and homogeneous electric E and magnetic B fields.Its dynamics is governed by the Lorentz-force equation The universal constant c is the speed of light in vacuum, and charge to mass ratio we denoted by µ = e m .This equation can be solved in elementary functions [18].One integration can be performed directly, giving a first-order equation with three integration constants being components of the vector When B = 0, particle moves with the acceleration directed along E: y(t) = y 0 + Ct + 1 2 µEt 2 .When E = 0, particle moves along a helical line located on a cylinder, the axis of which is directed along the vector B.
To find an explicit form of the solution when both E and B are presented, we choose the Laboratory system with axis z in the direction of B, so that the vector E lies in z, y plane, see Fig. 4. Consider the problem in the variables y(t) = y ∥ (t) + y ⊥ (t), where y ∥ (t) is the projection of our particle on axis B while y ⊥ (t) is the projection on the plane orthogonal to B. Similarly, we separate C = C ∥ + C ⊥ and E = E ∥ + E ⊥ .Substituting these decompositions into Eq.(90), we get separate equations for two projections In obtaining of second term on r. h.s. we used that [E ⊥ , B] = |E ⊥ ||B|e 1 , see Fig. 4. By choosing the time reference point, we can give any desired value to the constant ω 0 , so we put ω 0 = 0. Choosing C ⊥ = c[E ⊥ , B]/B 2 , the first two terms in Eq. ( 98) vanish.This corresponds to the choice of initial position y(0) with y ⊥ (0) = ae 2 .That is at t = 0 the particle lie in the plane z, y.
These are parametric equations of a plane curve called trochoid, see [18].In resume, trajectory of a charged particle in constant and homogeneous electric E and magnetic B fields is where the point y ∥ (t) moves along the axis B according Eq. (93) while the point y ⊥ (t) moves along the trochoid (99) on the plane orthogonal to B.
Let us write a variational problem for the equations (89).To this aim, we introduce the scalar A 0 and vector A potentials at each spatial point y as follows:

FIG. 1 :
FIG.1: Choice of Laboratory system ei for analysis of a charged ball subject to electric and magnetic fields.

3 FIG. 2 :
FIG.2: Choice of Laboratory system for analysis of a charged symmetrical body.

FIG. 3 :
FIG.3:The third components of initial data vectors ω, lying on drawn hyperboloids, turn out to be solutions to the equation of cubic pseudo-oscillator (47).Figure (a) shows the upper sheet of the hyperboloid of a short solenoid.Figure (b) shows the hyperboloid of a long solenoid.All drawn vectors and line segments lie in the plane of paper sheet.The basis vector e1 is orthogonal to the plane of paper sheet and is not shown in the figure.

2 µE ∥ t 2 . 2 ,FIG. 4 :e 1
FIG.4: Choice of Laboratory system for analysis of charged particle subject to electric and magnetic fields.