Nonneutral Weibel model plasma in the non-minimal CPT-odd coupling

Inspired by the extension of the Standard Model, we investigate the effects of the space-time anisotropy caused by Lorentz symmetry violation (LSV) on a plasma column conﬁnement conﬁguration. The model of Plasma taken into account is the z-pinch model that was in the earliest efforts in fusion power research. This model comprises particles in a nonequilibrium cylindrical distribution, which remains stationary in the absence of collisions. We propose a disturbance in the distribution by a Lorentz violation environment. As proposed by Carroll, Field, and Jackiw, in a scenario of (LSV), a background ﬁeld vector could couple with the electromagnetic ﬁeld, modifying the classical behavior of the electromagnetic ﬁeld. As reported here, considering the presence of a background ﬁeld vector, the intensities of the ﬁelds and particle densities would be disturbed by the inﬂuence of the LSV. For different values of the background ﬁeld vector coupling constant, the contribution of the background vector ﬁeld could modify the intensity of the electromagnetic ﬁelds, and concentrate even more the electrons densities in the edge of the plasma column, evidencing a behavior similar to a skin effect in this plasma column.


Introduction
Thermonuclear fusion research started in the 1950s after the development of the hydrogen bomb.The development of the peaceful research counterpart appeared a few years later, with the first pinch experiments [1].The first two attempts were made using the simple schemes of θ -pinch and z-pinch.To clarify, θ and z refer to the direction of the plasma current in a e-mail: diego.n.soares@edu.ufes.br(corresponding author) b e-mail: belichjr@gmail.comc e-mail: wesley.spalenza@ifes.edu.brd e-mail: filipe.braga@ifes.edu.brterms of a cylindrical coordinate system r , θ , z.In addition to historical factors, from the point of view of technological development, both systems present simplicity in experimental implementation and theoretical treatment when compared to more complex, but more effective, forms of confinement, such as stellerators and tokamaks [2,3].Both models rest on cylindrical approximations, however, in the Z-pinch case the coiling of the confinement coils occurs in the z direction of the cylinder, and in the θ -pinch case it happens along the θ direction of the cylinder, that is, at various points of the cylindrical column.The plasma confinement configurations in an infinite Z-pinch theoretical model resemble the approximations used in large aspect ratio tokamak reactors.Among several approaches to plasma modeling in z-pinch, Weibel's work from 1959 [4] brings one of these simple and practical approaches, which was recently revised [5].
It is simple to produce plasma by ionizing hydrogen gas in a tube, and a conductive fluid is obtained by a strong current induced by discharging a capacitor bank over an external coil surrounding the gas tube.From these studies, problems related to the equilibrium of plasma columns confining using complex configurations of magnetostatic fields reported by [6,7], bicuspid arrangements investigated by [8], or self-organizing structuring as observed by [9] have been extensively studied in the literature [4,[10][11][12][13][14]].An investigation point of these investigations was related to the magnetic confinement of the fusion plasma [15].On the other hand, an investigation point that we intend to open up with this work is to use plasma to obtain background fields that can provide us with physics different from the Standard Model of Particle Physics.
The Weinberg-Salam-Glashow Standard Model (SM) results from an intense search to understand the origin and behavior of elementary particles.The large hadron collider (LHC) has made the discovery missing for the acceptance of the model, the experimental verification of the Higgs boson.
Also, physics out of the Standard Model, such as the pentaquark [16], has started to emerge at LHC.Despite its success, the Standard Model ignores gravitational interaction and the problems related to it.So the SM has nothing to clarify about the dark Matter [17][18][19] and Energy problems.Thus, it is interesting to investigate possible extensions of the SM, and this line of research is known as Physics Beyond the Standard Model.With this goal, Kostelecký and Samuel [20] have shown that the Lorentz symmetry (LS) can be broken spontaneously.The development of this proposal by the works [21,22] has established what we today know as the Standard Model Extended (SME).
In the context of electrodynamics, Carroll et al. [23][24][25] have investigated the presence of a background vector field, ξ ν , coupled to a Chern-Simons-like lagrangean in (3 + 1)dimensions, μναβ ξ ν F αβ , the dual electromagnetic tensor, which violates Lorentz symmetry [26][27][28][29][30][31].This term violates CPT-symmetry and presents interesting modified electromagnetism, such that vacuum birefringence.Thus, we propose a non-minimal coupling [24], where μνκλ is a Levi-Civita pseudo-tensor, and the 4-vector ξ ν is a fixed vector and acts on a vector field that violates the Lorentz symmetry, in the Lagrangian, such that f represents the fermions in the theory, the current matter term, where the new covariant derivative is D μ , with g as a Lorentz breaking parameter, and the current J μ carries the presence of matter.
It is studied in various branches of physics, such that in magnetic moment generation [24], in Rashba spin-orbit interaction, in Maxwell-Chern-Simons vortices [32], on vortexlike configurations [33], in Casimir effect [34,35], in cosmological constraints in electrodynamic [36,37] and on an analogy of the quantum hall conductivity [38].In the nonrelativistic limit, there are studies in an Aharonov-Bohm-Casher system, in quantum holonomies, in a Dirac neutral particle inside a two-dimensional quantum ring, on a spinorbit coupling for a neutral particle.In the relativistic case, there are studies on Einstein-Podolsky-Rosen correlations, geometric quantum phases, Landau-He-McKellar-Wilkens quantization, and bound state solutions for a Coulomb-like potential, quantum scattering [39] and in a scalar field [40][41][42][43][44][45][46].
In this article, we investigate the effective electromagnetism that a plasma undergoes in the presence of a quadrivector that violates the Lorentz symmetry coming from the non-minimum coupling.The structure of this paper is as follows: in Sect. 2 we study the influence of space-time anisotropies under a plasma environment by our non-minimal coupling, in the Sect. 3 we study the modification in electromagnetism, caused by the Lorentz breaking sector in the effects of the z-pinch model, in the Sect. 4 we present our conclusions.

Plasma with Lorentz symmetry violation
In this section, we investigated the influence of the spacetime anisotropies under a Plasma environment by our nonminimal coupling.For this, we take the Lagrangian and analyze the term where we have LSV, and then we have the modified mechanical momenta with (gξ 0 , g ξ) being a perturbation quadrivector parameter ξ μ .Under that consideration, the main Hamiltonian associated with low-energy plasma particle motions can be depicted as Using the connection between the Electromagnetic fields and the potentials we observe that E = −∇φ − ∂ A ∂t , B = ∇ × A, and given that we are interested in stationary configurations ∂ A ∂t = 0. Again, one may see that E( r ) = −∇φ( r ), B( r ) = ∇ × A( r ).Then, the Hamiltonian becomes Given the axial symmetries imposed by the Weibel model, we can observe that A = A z êz , and in a similar way we can put φ = φ(r ), then we encounter the following equations Looking more explicitly for the components obtained inside the momentum term at the Hamiltonian, we observe that ξ = ξ r êr + ξ φ êφ + ξ z êz , and considering that we are able to attribute any value to the proposed perturbation quadrivector (ξ 0 , ξ 1 , ξ 2 , ξ 3 ) only using the condition that the spatial part of it maintains its norm equal to the unity.Considering that the E.M. equations associated with the potentials for the new proposition can be written as (restrict-ing our analysis for stationary cases) 3 Non-neutral Plasma pinch effect with CPT-odd coupling Now, we will apply the considerations of the modification in electromagnetism, caused by the Lorentz breaking sector in the effects of the z-pinch model, initially described by Weibel in his paper [4], to recover this theory and verify the new effects of a given parameter g.
In other to evaluate some conclusions about J , we take the constant temperature parameter β = (k B T ) −1 , in canonical Ensemble to describes Plasma fluids F ± = a ± ex p (−βH ± ) δ( p z ), such that, where a ± are integration constants.
Introducing the space-like condition from the Lorentz vector ξ μ = (0, ξ ), where ξ = êφ , and then we choose (ξ 0 , ξ 1 , ξ 2 , ξ 3 ) = (0, 0, 1, 0), which produce and where we use J = J (r ) êz , and ρ = ρ(r ) as proposed by [4], such that ρ = en 0 (n + − n − ) and J z = J + + J − .The potential equations can be written as for the electric sector, and for the magnetic sector.Now, using the following parameterization based on Weibel's paper parameterization [4] and considering that terms with e 2 , g 2 and e g we can neglect, due to their orders of magnitude, knowing that they usually appear multiplying the potential function, generating terms of negligible values.
x, and the fields e A(r ) Being a ± represents the normalization constant of each distribution F ± .Here, y and u represent, respectively, the dimensionless electric potential and the magnetic vector potential.x is a dimensionless radial coordinate.n 0 is a typical plasma particle density in the center of the cylinder, taking into account ion and electrons, M = m + m − m + +m − is a reduced mass for the species analyzed, and β = is the dimensionless difference of mass between the ions and electrons [4,5].There is a need to redefine the derivatives depending on the new parameters dr = and Then, replacing these expressions of densities in the differential equations that we have before, we achieve this new equation system that allows us to find the fields equations and To find v(x) and u(x), we use algorithms for numeric methods to solve this ordinary differential equations system, using the Runge-Kutta (RK) method, in the Fortran-90 language, exploring the distinct values of the parameter g.In order to run the RK code, we use a computer with 2 processors Intel Quad Core Xeon E5540, 2.53 GHz, 8 M L3 cache, 5.86GT/s, and 24GB of memory, DDR3 ECC SDRAM, 1066 MHz.Note that if we take the parameter g = 0 we recover the expressions of the classic model explored in the papers [4,5].Thus, using g = 10 −15 and the same initial conditions in an algorithm, we achieve similar results for the fields and densities in Fig. 1.Considering the dimensionless terms associated with the arguments of exponential functions, it can be noted that the parameter g can be associated with the ratio between the magnetic interaction (magnetic pressure) and the electrostatic interaction (electrostatic potential energy) merged.See the demand for dimensionlessness of the product Using the input parameters from Weibel's [4] original article, g would be on the order of 10 −10 .
By changing the values of the coupling constant g, varying by one order of magnitude at a time, between the values of g = 10 −15 and g = 10 −7 , the profiles are of the same form.Although, when we achieve the value of g = 10 −6 we observe an alteration in the form of the profiles, as we see in the Fig. 2.
For the values of "g" between g = 10 −6 and g = 6×10 −6 , it was possible to see different alterations in the curves that present the field intensities and particle densities, as we see in the Fig. 3.
Analyzing the densities graphs is possible to see that in the same way that the column radius grows up, the densities of particles go to zero.This means that, in that parameter configuration, the plasma column is confined.Furthermore, it is possible to note that the parameter g can alter the profile of particle densities.As far as the coupling constant value increase, the densities of the electrons became more concentrated, assuming a profile similar to a Dirac's delta distribution function, evidencing even more the profile of skin effect [47].The phenomenological effect of electron condensation at the plasma edge, directly points to a reinforcement of the conductive properties of non-neutral plasma without collisions.These characteristics may be associated with critical conditions.For the values of "g" above that we used in the algorithm, larger than g = 6 × 10 −6 , the algorithm was unable to complete the resolution.
Exploring more possibilities related to the preferential direction of the background field vector, we investigate the perturbation contribution if the background vector had the preferential direction in (0, 0, 0, ξ 3 = 1) and (0, ξ 2 = 1, 0, 0).However, in these scenarios, the contribution of the perturbation factor is negligible, returning to the results presented in [4] and [5].
To improve our analysis, we calculate the ratio between electromagnetic energy and thermal energy, following the method used in [5].The electromagnetic energy is given by Using thermal energy as the integration of an ideal gas approximated energy.
Using the same change of variables proposed by [4], we achieve the following expressions.
As a result, of analysing the confinement configurations, the ratio between the electromagnetic energy and the thermal energy indicates a value bigger than unit E E M /E T > 1, meaning the electromagnetic field had energy enough to contain the expansion of the plasma column.Calculating this energy ratio for the different confinement configurations that we analyze and for the different values of the "g" parameter that we analyze the plasma column, the ratio E E M /E T > 1 is satisfied.It means that the behavior of thermal expansion emulating an ideal gas is suppressed by the electrostatic and magnetostatic field pressure.plasma electromagnetic confinement scenario known as zpinch.Initially, we reproduced the results of the classical model, developed by Weibel [4] and [5], using a Fortran-90 code, to solve the differential equation system from the development of the new model, considering that the values of the coupling constant "g" were null or very small (Fig. 1).
Increasing values to the coupling constant "g", we introduce the perturbation effect from the term of Carrol, Field, and Jackiw that appears in the graphs of field intensities and electronic densities.For values of the order g = 10 −6 , the fields and electronic densities begin to suffer the influence of the anisotropy caused by the LSV, modifying the aspect of the fields and electronic densities curves (Fig. 2).In the range of g = 1.10 −6 and g = 6.10 −6 , was possible to detect more modifications in the intensities of the fields e electronic densities (Fig. 3).
For values in the order of 6×10 −6 , we note an instability region, that collapses the intensity of the field, electronic densities, and the confinement configuration of the plasma column.
According to these results, in a scenario of LSV, a background field vector would couple with the electromagnetic field, disturbing the classical electromagnetic field and consequently, disturbing particle densities as well.Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made.The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material.If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder.To view a copy of this licence, visit http://creativecomm ons.org/licenses/by/4.0/.Funded by SCOAP 3 .