Internal force distributions in 't Hooft-Polyakov monopole and Julia-Zee dyon

The energy-momentum tensor of the 't Hooft-Polyakov monopole and the Julia-Zee dyon are studied. This tensor contains important information about the pressure and the shear force distributions which define the mechanical properties of systems. Obtaining the violation of the local stability criterion for the magnetic monopole and dyon we decompose the EMTs into long- and short-range parts. This decomposition depends on the abelian field strength tensor which can not be uniquely defined. We suggest to use the modified 't Hooft definition for the tensor. Finally, the long- and short-range parts of the EMTs are computed and new equilibrium equations are obtained. Numerical values for masses, $D$-terms and various mean square radii for the monopole and the dyon are also computed.

The most interesting part of the EMT for our purpose is the ij-components, which define the stress tensor. The T ij components can be associated, according to Ref. [5], with the distribution of the shear forces s(r) and the elastic pressure p(r) inside the system. For a spherically symmetric system in static approximation in three dimensions the stress tensor is decomposed in s(r) and p(r) as follows T ij ( r) = r i r j r 2 − 1 3 δ ij s(r) + δ ij p(r).
The pressure and shear force distributions can be connected with the gravitational D(t)-form factor, which is involved in the parametrisation of the matrix element of the EMT operator, see e.g. Ref. [5]. As it follows from Ref. [5], the D-term (or the Druck-term) D ≡ D(0) can be obtained as where M is the mass of the system. The value of the D-term for a particle can be measured experimentally, like its mass, however, the value is difficult to extract from experimental data. For example, the experimental value of the D-term for the proton is D = −1.47 ± 0.06 ± 0.14, where the first error is the statistical uncertainty, and the second error is due to the systematic uncertainties, see Ref. [6]. The value of the D-term for the nucleon was first obtained theoretically in the Skyrme model in Ref. [7], the authors got the value D = −3.6 and in the bag model the value of the D-term for the nucleon is D = −1.1 [8], other values of the D-term together with the corresponding references can be found in the table 2 of Ref. [9]. The values of the D-term for the monopole and the dyon are given in Secs. III and IV, respectively.

B. Static EMT
The EMT is the conserved Noether current associated with spacetime translations, i.e. it satisfies the following condition For the static EMT the Eq. (5) turns to ∂ i T iν = 0. For the parametrisation of Eq. (3) it implies that the pressure and shear force distributions satisfy the following equation which is also called equilibrium equation ∂p(r) ∂r + 2 3 ∂s(r) ∂r + 2 r s(r) = 0. (6) Assuming that the pressure and shear force densities decay at large distances faster than ∼ 1 r 3 , multiplying on r 3 and integrating over r one obtains the von Laue stability condition given in Ref. [10] d 3 rp(r) = 0, which is necessary for stability, but not sufficient, since it is also satisfied for unstable systems, see discussion in Ref. [4]. This condition is satisfied for any system, whose EMT is conserved and implies that the pressure distribution has at least one mode. Integrating the differential equation (6), from positive value r to infinity and requiring that s(r) and p(r) vanish at the infinity, one obtains This equation describes the equilibrium of the internal forces inside a system. The combination 2 3 s(r) + p(r) describes the normal component of the total force exhibited by the system on an infinitesimal piece of area dS i , which is denoted as p r (r): F i (r) = T ij (r)dS j = 2 3 s(r) + p(r) dS i = p r (r)dS i (9) where dS i = dS r i /r. From the equilibrium equation (8) follows that for the positive shear force distribution s(r) the normal force is always positive, for the negative s(r) the normal force is always negative. In Ref. [4] it was argued that the normal force has to be directed outwards otherwise the system would collapse, this means that it has to be non-negative, i.e. 2 3 s(r) + p(r) ≥ 0. (10) Thus, according to the Eq. (8) for negative shear force distribution the system is definitely unstable. The local stability condition is also necessary, but not sufficient for stability analogously to the von Laue condition. However, due to its local character, it is stronger than the von Laue condition.

C. EMT and external forces
The EMT conservation taking the form ∂ i T ij (r) = f j (r) can be interpreted according to Ref. [11] as the equilibrium equation for internal stress and external force f j (per unit of the volume). According to the parametrisation in Eq. (3) the equilibrium equation gets the following form d dr 2 3 s(r) + p(r) + 2 s(r) r = f (r), (11) where f (r) is the normal component of the external force per unit of the volume. This equation describes the balance between internal forces pushing out from center and external force pulling inwards to the center. When the corresponding forces are equal, the system is at equilibrium. The von Laue stability condition for such equilibrium equation gets the following form One can also rewrite the equilibrium equation (11) as with σ(r) = ∞ r dxf (x). The left-hand side describes the normal component of a total force of the system acting on the infinitesimal unit area dS i as it follows from below Then for systems affected by an external force, i.e. for systems described by the equilibrium equation (11), the local stability criterium of Eq. (10) can be modified as follows For a time-dependent system the conservation of EMT gives the following equation ∂ i T ij = −∂ 0 T 0j , where the left-hand side again describes internal force of the system and the right-hand side can be interpreted as an external force.

D. Mean square energy radius and mechanical radius
For a positive energy density T 00 (r) the mean square radius of the energy density can be introduced as which characterises the size of the system in which the energy is distributed. According to the local stability condition, for a stable system the radial force must be positive, so the mechanical mean square radius where the normal force is distributed can be defined as in Ref. [12], for a system with the equilibrium equation (6) it takes the following form and for the system with the equilibrium equation (11) it takes has the form: where p r (r) and σ(r) are defined above. The local stability condition has not been mathematically proven and thereby is still questioned, see e.g. criticism in Ref. [13]. Moreover, as it was mentioned in recent studies of Refs. [14][15][16], the local stability condition is not satisfied in the presence of the long-range forces. Additional to that, as pointed out in Refs. [14,15,17,18], there is another problem for systems with the long-range contribution, namely, the divergence of essential quantities, that describe mechanical properties of a system, such as the D-term and the mean square radii corresponding to the EMT densities.
At this point, it is important to mention, that the 't Hooft-Polyakov monopole is accepted to be a stable system, see arguments of Refs. [19,20]. As we will see later, the local stability condition is violated for the monopole and the dyon. On the first glance, thereby one could think that the criterium is not correct. However, the monopole and the dyon involve electromagnetic interaction, which supports the idea that the local stability condition does not apply correctly in the presence of long-range forces.

A. Equations of motion
The Grand Unified Theories combine the electromagnetic, weak, and strong forces into a single force. One of the first such theories was suggested by Howard Georgi and Sheldon Glashow in 1974 in Ref. [21]. Later Alexander Polyakov and Gerard 't Hooft have independently found that magnetic monopoles automatically appear in all Grand Unified Theories [1,2]. The most simple model, where magnetic monopole exists is the gauge SU(2) Georgi-Glashow model with Higgs triplet field ϕ a , a = 1, 2, 3, which belongs to adjoint representation. The corresponding gauge invariant action of the model is where λ is a dimensionless coupling constant, v 2 is the squared vacuum expectation value of the Higgs field; µ and ν are Lorenz indices. The covariant derivative and the non-abelian field strength tensor are defined as follows where g is the gauge coupling constant and A a µ is the gauge vector field. Choosing the vacuum expectation value as ϕ a 0 = (0, 0, v) and considering small fluctuations of the ground state in unitary gauge, one can show that the Georgi-Glashow model has one massive scalar field η 3 (x) with mass m H = √ 2λv, one massless vector field A 3 µ corresponding to the U (1) subgroup of SU (2) gauge group and two massive vector fields A 1 µ and A 2 µ both with masses m V = gv. In this model we are interested in a soliton solution, in other words, in a static solution of classical field equations with finite energy. As we want to study static soliton configuration, we consider the fields A a i ( x) and ϕ a ( x) independent of time. We also fix zero component of a vector field through the gauge A a 0 = 0 to have zero electric field. Requiring the energy to be finite the following configuration of fields can be found ϕ a = n a vh(r), with the unit vector n a = r a r , the unknown profile functions h(r) and F (r) and the following boundary conditions F (r) = 0 at r → ∞, F (r) = 1 at r → 0, h(r) = 1 at r → ∞, h(r) = 0 at r → 0.
The explicit form of the profile functions can be obtained from equations of motion which can be computed by varying the action with respect to scalar and vector fields, which for the static case reduce to the following form In terms of the profile functions the equations transform to Rescaling the argument r with dimensionless argument ρ as ρ = r R0 , where R 0 = 1 m V = 1 gv is a typical size of the solution, and introducing the new parameter we obtain the following system of equations of motion with the same boundary conditions as before. This system can be solved analytically only for the limit β = 0, otherwise one solves it numerically. However, one can find approximate solutions at small and large distances: where C F , C h , a and b are free constants. Note, that for β = 2 the constant C F = 0. The coefficients c n and d n can be expressed in terms of a and b as follows In Ref. [22] another system of differential equations is obtained and analytically solved for complex non-abelian monopole and dyon fields.

B. EMT densities
The EMT for the 't Hooft-Polyakov monopole can be obtained variating generally covariant form of Georgi-Glashow action (19) with respect to the metric For the static case one gets Using the 't Hooft-Polyakov ansatz given in Eq. (21) and the decomposition in Eq. (3) the EMT densities can be expressed in terms of the profile functions F (ρ) and h(ρ): The T 00 (ρ) component provides the information about spatial distribution of the monopole mass, p(ρ) and s(ρ) describe the pressure and shear force distributions inside the monopole. The first two terms in these densities are originating from the term ∼ F a µσ F aσ ν in action, the next two terms are coming from ∼ D µ ϕ a D ν ϕ a , and the term proportional to β 2 comes from the Higgs potential.
Before we discuss the numerical results for the EMT density distributions and their properties, we compute their asymptotic behaviour at large and small distances using the behaviour of the profile functions from Eq. (26): and where C h , C F and a, b are free constants from the asymptotic behaviour of profile functions in Eq. (26). From the large distances behaviour it is clear that the power-law decay (∼ ρ −4 ) is dominating, this behaviour corresponds to the contribution of the electromagnetic interaction. 1 This is not surprising if one remembers that the mass spectrum of the Lagrangian (19) has one massless vector particle corresponding to the U (1) group. The asymptotic behaviour of the energy density is definitely positive at the infinity as well as at the origin. This hints that we are dealing with an usual system. The outer region of the pressure distribution has definitely positive sign and that of the shear force distribution -negative one. Such behaviour leads to violation of the local stability criterion (10). However, this behaviour is related only to the electromagnetic contribution. Note, that because of this long-range contribution both mean square radii in Eqs. (16) and (17) as well as the Druck-term in Eq. (4) diverge.
In Fig. 1 the energy density for various values of β is shown. Although the energy is growing with β, according to Ref. [23] it does not diverge for β → ∞. The distributions of the shear force, pressure and normal force are shown in Fig. 2. Although the pressure is mostly negative, it changes the sign which allows the von Laue condition (7) to be satisfied, see e.g. Fig. 6. It is also interesting that the maximum of pressure is three times larger than the maximum of shear forces like in many other systems studied before, see e.g Refs. [4,5,7,12,14]. From the figure one also sees that the stability condition (10) is violated everywhere and for every choice of β, although it is proved in Refs. [19,20] that the monopole is stable. Since the local stability condition is not mathematically proved, we cannot conclude from its violation that the monopole is unstable. This violation can be also explained by the fact that the stability criterion is inapplicable for systems with the long-range force that is present in the monopole. The authors of Refs. [14,15] also obtained the violation of the stability criterion for stable systems like proton in the presence of long-range forces. All these aspects motivate us to think that the stability condition in Eq. (10) can be applied only for the systems where only the short-range interactions are present.

C. BPS limit
The equations of motion can be solved analytically only for the limit β = 0. The solution of these equations was first obtained in Ref. [24]. In Ref. [25] Bogomolny solved these equations by reducing the system of the second order equations to the first order by considering the energy functional for the Georgi-Glashow model. By integrating the energy density over the volume one gets the energy functional Rewriting it in terms of the chromomagnetic field H a i = − 1 2 ijk F a jk , Bogomolny introduced the following inequality for the monopole energy, which is also called BPS bound where the Bianchi identity (D i H i ) a = 0 and the Gauss's theorem were used. The BPS limit gives the possible minimum of the static energy of the monopole and the equality holds for the Bogomolny equation From the definition of chromomagnetic field the following expression can be obtained Using this expression and Bogomolny equation (35) it can be shown that the spatial components of the EMT disappear: independently of the choice of the ansatz in Eq. (21). The vanishing of the pressure and shear force distributions at the BPS limit allows to assume that the main contribution for the stress tensor T ij is given by the Higgs potential. As it can be seen from Fig. 1 the largest contribution to the pressure distribution is indeed due to the term originating from the Higgs potential in Eq. (19). Although the shear force distribution does not depend explicitly on β, it is related to the pressure distribution due to the differential equation (6). We think that in all models where it is possible to construct the BPS limit, the corresponding ij components of EMT would vanish. We have also obtained such vanishing, for example, for baby Skyrme model. The vanishing of the pressure and shear force distributions for BSP limit indicates that the 't Hooft-Polyakov monopole in this limit has the isotropic matter distribution [12].

D. Abelian field strength tensor and electromagnetic EMT
In the previous subsection we obtained the EMT in non-abelian SU (2) gauge theory, however, the theory has one massless vector field corresponding to the abelian U (1) subgroup. This massless vector field is responsible for the electromagnetic contribution to the EMT and for the divergence of such mechanical properties of the monopole as the D-term, mean square energy and mechanical radii. We will subtract this long-range contribution from the EMT and study the remaining short-range structure. For this we have to determine the U (1) abelian field strength tensor. Let us denote the abelian field strength tensor as F µν . Since in the unitary gauge the massless vector field of the theory is the third component of the vector field A µ 3 (see Sec.III A), the expression for the abelian field strength tensor must be SU (2) gauge invariant and coincide with . We introduce a general definition of the abelian field strength tensor without any requirements on φ as The constant c 1 is not fixed by the above requirements for the abelian field strength tensor F µν , however, it is needed to define the magnetic charge density in such a way that it coincides with the topological charge density. G. 't Hooft suggested in Ref. [2] the tensor with the unit vector φ a =φ a = ϕ a |ϕ| and c 1 = 1, Faddeev offered in Ref. [26] the tensor with φ a = ϕ a v and c 1 = 0, where ϕ a is the ansatz from Eq. (21) and Boulware in Ref. [27] suggested the tensor with φ a =φ a = ϕ a |ϕ| and c 1 = 0. Since the under-lying theory combines the long-and short-range forces it is not possible uniquely defined the U (1) field strength tensor to separate the short-range interaction from the long one.
For the general choice of the abelian field strength tensor the magnetic charge density gets the following form One can show that for the 't Hooft's definition of the abelian field strength tensor the magnetic and topological charge densities coincide and describe single point-like particle with the magnetic charge 4π g at the origin, see Ref. [2]. Such monopole is called Dirac's monopole. In Ref. [28] Dirac derived the quantisation of the magnetic charge as In contrast to the 't Hooft's definition of the abelian field strength tensor the Faddeev's and Boulware's definitions describe the magnetic charge density smoothly without singularities at the origin. The topological and magnetic charge densities are not equal in these cases, but the quantisation of the magnetic charge is satisfied. It is also interesting to notice that the magnetic charge density in Eq. (39) for the Faddeev's and and Boulware's definitions of the abelian strength tensor coincides with the energy density in BSP limit in Eq. (34) up to some dimensional normalisation factor. For the positive magnetic charge distribution the mean square radius can be determined We define the electromagnetic EMT analogously to the EMT in electrodynamic as where F µν is the U (1) field strength tensor. As we have discussed F µν is not uniquely defined. Since we are interesting in spatial structure of the monopole it makes no sense to define T C µν according to the 't Hooft's definition of the F µν because of singularities at the origin. Faddeev's and Boulware's smooth definitions of the abelian strength tensor F µν could be a better candidate to define the T C µν , however, we suggest another definition of the U (1) field strength tensor which is also smooth as the Faddeev's and Boulware's definitions and it coincides with the 't Hooft's definition at long distances. For this definition of the abelian strength tensor the electromagnetic EMT distributions have the following form where The function Q(r) is directly related to the magnetic charge density through The electromagnetic EMT densities are shown in Fig. 3. The corresponding asymptotic behaviour can be found in App. A3. Note, that since the full spatial part of the EMT in BSP limit where β = 0 vanishes, see Sec. III C, and the electromagnetic EMT does not depend on β (and thereby does not vanish in this limit) the short-range part of the EMT in Eq. (45) has to be equal to the long-range part. It means that the scalar interaction in BSP limit becomes long-range interacting.

E. Short-range part of EMT and its consequences
We denote the short-range part of the EMT as T SR µν and define it as follows where the full EMT T µν is obtained in Eq. (30) and electromagnetic EMT -in Eq. (43). After some simple algebraic calculation the following expressions for the short-range part of the EMT densities can be found The corresponding asymptotic behaviour can be found in App. A4. The short-range part of the pressure and shear force distributions are presented in Fig. 4. The short-range part of the pressure distribution is negative and the shortrange part of the shear force distribution is positive for every choice of β in considered range of r. The total normal force distribution is positive for any value of β, which satisfies the stability criterium (15). It is also interesting to notice that the pressure distribution of the monopole does not change drastically after the subtraction of the long-range part, however, the shear force distribution does. As we have already discussed in the first section, the static EMT according to the Noether theorem must be conserved. The conservation of EMT couples pressure and shear force dinsities by the differential equation (6). However, after the decomposition of the EMT into the short-and long-range parts, the separate EMTs are not conserved and the corresponding short-and long-range parts of the pressure and of the shear force densities couple due to the new equilibrium equations: d dr where Q M (r) = | x|<r d 3 xρ M ( x) = 4π g Q(r) is magnetic charge contained in a sphere of radius r. So the right-hand side describes "Coulomb force" of the magnetically charged sphere acting on the magnetic charge density. 2 Thereby the first equation is the equation of magnetostatic equilibrium between the "Coulomb stress" pushing the monopole outwards and the magnetic "Coulomb force" pulling the monopole inward to the center. In contrast, the second equation describes the balance between the "short-range stress" pulling the monopole inward to the center and the repulsive magnetic "Coulomb force" pushing the monopole outward. We notice, that since the right-hand sides of the equilibrium equations are associated with the magnetic charge density, it is not possible to define the long-and short-range parts of the EMT using the ambiguity of the F µν in such a way, that the decomposed EMTs would be conserved separately, i.e. the right-hand sides would vanish, except the case with zero magnetic charge density. As we discussed, the D-term of the monopole diverges because of the long-range contribution. After subtracting the long-range part from the EMT, the corresponding D-term still diverges for small values of β. We illustrate this divergence in the following. The D-term is defined through the shear force distribution as in Eq. (4). Since for small values of β 1 the main contribution to the short-range part of the shear force density is provided by the asymptotic behaviour at large-distances, see App. A4, the D-term diverges for the large enough R as The mean square magnetic charge radius in Eq. (40), the mean square energy radius of the short-range part in Eq. (16) and the mean square mechanical radius in Eq.  (4), the mean square radius of the magnetic charge density is computed due to the definition in Eq. (40), the mean square radius of the short-range part of the energy distribution is computed according to the definition in Eq. (16) and the short-range part of the mechanical square radius -according to Eq. (18). It is remarkable that for the small values of β the main contribution to the monopole mass gives the short range part and for the large values of β 5 -the long range part. It would be interesting to find the analytic form of the pressure and shear force densities depending on the parameter β like it was done in Refs. [23,29,30] for the static energy.

A. Equation of motion
In the previous section we considered monopole solution carrying magnetic charge. In this part we will consider a dyon, that is a hypothetical particle carrying both electric and magnetic charges, such particle was first suggested by Julian Schwinger in 1969 in Ref. [31]. The dyon solution in Georgi-Glashow model was first obtained in 1975 by Anthony Zee and Bernard Julia in Ref. [3]. We will use the same Georgi-Glashow model as we did in the first part with the same action (19). However, we will not fix the zero component of the vector field as we did it for the monopole case, A a 0 ( x) = 0. We are again interested in the static soliton solution with finite energy. Thereby, we choose the same spherically-symmetric ansatz as we chose for the 't Hooft-Polyakov monopole in Eq. (21) with the same boundary condition as in Eq. (22). In contrast to monopole case, for the dyon the A a 0 is non-zero and one chooses the ansatz which follows from requiring the energy to be finite and definition of the abelian field strength tensor as in Eq. (42). The boundary conditions for the function J(r) are where Q D is an electric charge of a dyon, Q M = 4π g is a monopole charge and m is a free constant, that has dimension of mass. From the variation of the action follow the equations of motion The equations for the profile functions in dimensionless variable ρ = gvr are Note, that the boundary condition for J(ρ) also changes: The approximate solution for small and large distances can be found: The a, b, c are free constants, all other constants for the small r behaviour can be expressed in terms of these constants, for example, Since we searched the asymptotic behaviour at large distances for real and falling functions we define a new function J(ρ) =J(ρ)ρ and find that the constant C has restricted region: 0 ≤ C < 1. From numerical analysis of the solutions we obtain that the charge ratio is restricted in the following range 0 ≤ Q D Q M ≤ 1. More detailed analysis can be found in Ref. [32], where also the dependence of the parameter C on the charge ratio Q D Q M is found.

B. EMT densities
The EMT for the Julia-Zee dyon can be obtained analogically to the 't Hooft-Polyakov monopole by variating generally covariant form of the Georgi-Glashow action (19) with respect to the metric and it can be decomposed into magnetic and electric parts: where the magnetic part equals to the EMT of monopole which we have already computed in the previous section, see Eqs. (29), (30) and (31). The electric part of EMT is Note, that same as for the monopole case the spatial components of the EMT for the dyon vanish in BSP limit, where β = 0. For electric part of the energy density, pressure and shear force distributions of the dyon in dimensionless variable ρ the following expressions are obtained With the help of asymptotic behaviour of profile functions in Eq. (53) the asymptotic behaviour of the EMT for the dyon can be obtained From this behaviour one sees that the long-range contribution presented in the dyon is even stronger than in the monopole case, see Eq. (31). In App. A 2 one can find the asymptotic behaviour of the EMT near the origin as well as behaviour of the electric part of the EMT only. Again already from the asymptotic behaviour it is clear that the stability condition of Eq. (10) is violated. In Fig. 5 the full pressure and shear force distributions of the dyon are presented. The shear force and pressure distributions have very similar form as it was in the monopole case, see Fig. 2. These distributions are mostly negative for every choice of parameters. However, pressure distribution changes its sing, what allows the von Laue condition (7) to be satisfied, see e.g. Fig. 6. It is also interesting to notice that the main contribution to the dyon energy comes from the monopole part as it can be seen in Fig. 7. Moreover, for the growing β this contribution is increasing. The pressure distribution has one mode for any choice of the parameters.

C. Electromagnetic and short-range parts of the EMT
According to the definition of the electromagnetic EMT in Eq. (41) the long-range part of the EMT for dyon is Using the modified 't Hooft abelian field strength tensor in Eq. (42) one gets the following expressions for the electromagnetic EMT densities of a dyon where the function Q is defined in Eq. (43) and is related to the magnetic charge density in Eq. (44) and the functioñ Q is related to the electric charge density of a dyon Comparing these expressions with the expressions for the monopole case in Eq. (43), it is clear that both have similar behaviour. Hence, it is again not possible to define the various mean square radii or D-terms, because they diverge. We will exclude the long-range contribution given in Eq. (60) from the EMT of dyon in Eq. (55) in the same way as we did it for the monopole case using Eq. (45). After simple algebraic calculations we obtain the following expression for the short-range part of the EMT of the dyon Here The equilibrium equation of the long-and short-range parts of the EMT analogously to Eq. (47) is given as d dr where Q M (r) andQ D (r) are magnetic and electric charges contained in the sphere of radius r: The first equilibrium equation describes the balance between two types of "Coulomb forces", namely magnetic and electric ones, pulling the dyon inwards to the center and the "Coulomb stress" pushing the dyon outwards. The second equation describes the balance between the two types of the repulsive "Coulomb forces" pushing the dyon outward and the "short-range stress" pulling the dyon inwards to the center. Same as for the monopole case it is not possible to define the short-range part of the EMT using the ambiguity of the F µν in such a way that the long-and short-range parts of the EMT are conserved separately unless the electric and magnetic charge densities both vanish.  In table II the full mass of dyon is denoted as M , the monopole contribution to the full dyon mass as M M , the M C and M SR are the long-and short-range contributions to the full mass, correspondingly, the D-terms are computed with the help of the short-range part of the shear force distribution in Eq. (62). The following quantities can be also found in the table: the energy mean square radius r 2 E , the magnetic and electric charge mean square radii r 2

V. SUMMARY AND CONCLUSIONS
In this work the EMTs of the 't Hooft-Polyakov monopole and the Julia-Zee dyon are studied. These EMTs contain long-range contributions which present analogies to the "Coulomb interaction". The local stability condition containing the pressure and shear force distributions is violated for both cases, the monopole and the dyon. Thereby, the applicability of the condition in the presence of long-range contributions is questioned. Moreover, such important quantities which give information about mechanical properties of a system as the D-term, energy and charge mean square radii of the monopole and the dyon can not be computed due to the presence of the long-range interaction.
To shed more light on the mechanical properties of the monopole and the dyon and on the local stability condition we exclude the long-range contribution from the EMT of the monopole and the dyon. The difficulty of such calculation is that this contribution can not be uniquely defined. We suggest the modified 't Hooft definition of the abelian field strength tensor for this purpose. Dealing with the separate long-and short-range parts of the EMTs we obtained the equilibrium equations which couple the pressure, shear force distributions and the external force acting on the system. In this line we have also modified the local mechanical stability condition for systems in the presence of external forces. According to this condition the short-range part of the 't Hooft-Polyakov monopole as well as of the Julia-Zee dyon is stable for every choice of the parameter β.
Further, in this paper numerous figures describing mechanical properties of the 't Hooft-Polyakov monopole and the Julia-Zee dyon can be found as well as tables with varied contributions to the masses, D-terms and mean square energy, magnetic and electric charge radii.
The asymptotic behaviour of the magnetic charge density is The "Coulomb EMT" behaves asymptotically as (A3) The asymptotic behaviour of the short-range part of the EMT is

For dyon part
The functionQ(ρ) has the following asymptotic behaviour The asymptotic behaviour of the electric charge density is The asymptotic behaviour of the full EMT of the dyon near the origin is The asymptotic behaviour of the electric part of the full EMT for the dyon: The asymptotic behaviour of the long-range part of the EMT is The asymptotic behaviour of the short-range part of the EMT of the dyon is as follows:

(A10)
Appendix B: Divergence of radii As it follows from the Eqs. (A2) and (A4) the main contribution to the magnetic charge density and to the shortrange part of the energy density, respectively, is provided by the asymptotic behaviour at large-distances. Then the corresponding mean square radii for the small values of β and for the large enough R diverge as According to the modified definition of the mechanical radius in Eq. (18), to the external force in equilibrium equation (47) for the short-range and to the asymptotic behaviour in Eqs. (A1) and (A4) the mean square mechanical radius of the monopole for small values of β diverges as follows Proceeding analogously to the monopole case, we obtain for the dyon the same divergences as in Eqs. (48), (B1), (B2) and (B3). Namely, the D-term diverges as in Eq. (48), the mean square electric and magnetic charge radii as in Eq. (B1), the mean square energy radius of the short-range part as in Eq. (B2) and the mean square mechanical radius as in Eq. (B3).