The Magnetic Field in the Stability of the Stars

The theoretical frame developed by A. S. Eddington is used in this paper in order to study the problem of internal structure and stars stability. A new hypothesis is developed that consists of assume that at early stage of their evolution, all gaseous stars generate an intense magnetic field whose mission is to contribute to their stability. A simple methodology is shown to calculate the magnitude of the self-generated magnetic field at any point of the star, and diagrams are built for two specific cases, where the way in which this field varies from the center to the surface can be seen. Finally, the problem of the Cepheid type variables is studied obtaining an expression for the oscillation period carried out from the differential equation of the simple harmonic movement in terms of the appropriate parameters of the star.


Introduction
The theory developed by Eddington considers the case of a star composed by a fluid which from the thermodynamic viewpoint behaves as an ideal gas.To determine the state of the material, two ordinary differential equations are proposed in order to mathematically express the following conditions [1].
The mechanical equilibrium of the star.In order to fulfil this condition, it is necessary that at any internal region of the star, the pressure has the right value in order to support the weight of the material above it; and The existence of thermal equilibrium.This condition requires that temperature distribution in the star is capable of maintaining itself automatically not with standing the continual transfer of heat from one part of the star to another.
The proposed equations are integrated and conditions of the material at any point are determined.This way the distribution of pressure, mass density and temperature are obtained.As work hypothesis, it is considered that at first approximation, the scheme distribution is homologous from star to star.This means that all gaseous stars copy the same stellar model within its appropriate scale of mass, length, temperature, etc. [1].In order to, once and for all, simplify the task, a general solution to the problem is formulated and then the question is reduced to adapt it to the scale of the particular star being studied.It is important to notice that Eddington abandons J. Homer Lane's hypothesis [1] [2] when he does not consider that the thermal equilibrium in a gaseous star is due to the existence of internal convective streams.The hypothesis of convective equilibrium is replaced by the radiative equilibrium and this last basic principle is applied to the internal conditions of the sun and the stars [1].
In this work, a magnetic field as the fundamental element for the stars equilibrium will be introduced.It is proposed as a basic hypothesis, that all gaseous stars generate an intense magnetic field at an early stage of their evolution.Thus, a large mass of compressible gaseous fluid, viscous and conductor isolated in space and at very high temperature and pressure conditions that remain together by its own gravitational attraction and at dynamic equilibrium with the force produced by the sum of the pressures of radiation and the hot gases are considered.The source of power generation lies on the central part of this huge mass of gas and a little closer to the surface the so called convective zone [3] [4].The convective streams generated there are responsible of removing the heat produced and the ashes of the thermo nuclear process, feeding at the same time the power source with new nuclear fuel [5].It is assumed that these convective streams are made up by neutral atoms and a lot of electrically charged particles both positive and negative.It is also assumed that these convective streams contain a steady-state current distribution localized in some region inside the convective zone, produced by a process of maximum ionization, which by hypothesis is the generator of the magnetic field that in a first approximation, may be considered similar to the bipolar magnetic field produced by a magnetized bar.Irregularities observed in the magnetic field should be attributed to the fact that the structure of the localized steady-state current distribution is not always the same; this is because the condition of the convective streams is not always the same either, and can be dependent on the magnitude of the thermonuclear explosions produced in the oven.The solution of the problem of the magnitude and topology of the self-generated magnetic field by gaseous stars can be seen in Appendix II.It is significant to emphasize the importance that the self-generated magnetic field has for the equilibrium and the stability of this object.Additionally, this enormous concentration of matter is found distributed in a configuration that has spherical symmetry.A heavenly object with the above mentioned characteristics is a star.
It can be assured that at any point of the stellar fluid in equilibrium, there is a hydrostatic pressure which is the same in every direction.If a closed surface with in this configuration is drawn, the reaction of the external fluid to that surface over the internal fluid consists of a force per unit area exerted on the surface and along the inward normal.So that for the dynamic equilibrium existence, it is necessary for these surface forces to be balanced by body forces such as the gravitational which acts towards the inside of the star.It is also necessary to consider the contribution of the magnetic field.Its contribution to the state of equilibrium must be included in the theory in the way of a force per unit area known as the magnetic hydrostatic pressure [5]- [7].As it is well known [8], a magnetic field has the property of generating in a conducting fluid, a magnetic viscosity that confers it to certain rigidity.That rigidity in the conducting stellar fluid may be interpreted as if the star were supported by a superstructure made up by the magnetic lines of force.It will be assumed in this paper that this magnetic superstructure has the mission of keeping the shape of the star even when this object should be animated by a rotation movement, since the magnetic lines of force are frozen in the stellar fluid and move along with it [5] [9].

The Internal Structure and the Stability of a Gaseous Star
In order to make an appropriate theoretical treatment of the problem of the stability of a star and of its internal structure, it is necessary to include the self-generated magnetic field in the fundamental equation that governs the state of equilibrium [1] which in this case is magneto mechanical [5], i.e.
where p (x, t) is the whole pressure, ρ(x, t) the mass density, g the constant value of gravity acceleration and H 2 /8π the hydrostatic magnetic pressure, with ( ) being the magnetic field self-generated by the star [5].
The previous functions depend on time t and of an x distance measured from any internal region of the star to the center of it [1].The other relationship that must be considered is the equation for the radiative equilibrium [1] [10] In the previous formula p r (x, t) is the pressure of radiation, c the velocity of light in the empty space, ε the radiation energy by cm 2 and by second and k the coefficient which determines radiation absorption by stellar fluid [1].From ( 1) and ( 2) it is easy to see that Let us consider a fluid mass distribution within the star contained in a radius sphere equal to x.In a stationary state, the amount of energy released per second within the sphere is equal to L(x) and it is such that ( ) where 2 4π x ε is the whole of the radiation that flows per second through the surface of the sphere and it is clearly equal to the amount of energy released by the central thermonuclear oven [1].On the other hand, the gravitational force at the position x is only due to the mass M(x) which is found in the inside region to x [11], so that ( ) With G is the universal gravitational constant.From ( 4) and ( 5), the following result is obtained: ( ) ( ) where L(x)/M(x) is the average rate of energy release per gram from the inner region to x [1].It is assumed that this energy release is greater in the dense center of the star than in its external parts, in such way that this reason decreases as x is increased when succeeding layers of colder material are added to the average [1].Be M the total mass of the star and L the total emission of energy per second from its surface in such way that L/M is the boundary value that reaches the reason L(x)/M(x).For all of the above mentioned it is possible to write the following relation [1] ( ) ( ) There are reasons to state that η is a magnitude which is increased from 1 on the star surface up to an unknown value although not very big in the center of it [1].Its analytical form depends on the law of energy release by the thermonuclear oven [1].Nevertheless, it is possible to assume that this law is approximately of the same kind for all the stars since presumably the nuclear energy releasing mechanism is basically the same for all of them [1].If (7) is substituted in (3) the following result is obtained This last relation is an exact equation which allows establishing an upper limit for the value of the opacity for any star for which magnitudes L and M are known as a result of observation, regardless of whether the star is or not constituted by a fluid obeying the thermal equation of perfect gas [1].
In any case it is expected that the values of temperature and mass density are increased towards the inner part of the star [1].Consider that the whole pressure p is such that [1] g r where p g (x, t) is the pressure of hot gases.As on the other hand it is clear that in (8) the following fundamental relation is obtained Since the tendency of temperature and mass density is that of increasing its values towards the inner part of the star, it can be stated that the pressure of the material p g and the magnetic hydrostatic pressure H 2 /8π should increase their intensities so that the stability of the star is kept.Under these conditions the following inequality is expected to be fulfilled According to (1) Since η is always positive [1] we have that η > 1, in such way that from the previous relation the following is that with the constants numerical values is transformed into The previous numerical factor is only half of what is reported in the specialized literature [1].For the most brilliant component of Capella's binary system we have the following basic data: In that case For the Sun we have that k < 6564 cm 2 •gr −1 .The physical reason for the existence of that upper limit known by the name of Eddington limit [1] [10] is the following.So that the radiation is observed has to be emitted reaching its way through the star, and if there were too much obstruction it would blow up the star [1].However, the existence of the self-generated magnetic field is a determining factor for the star to keep its stability.In fact, the superstructure built by the magnetic lines of force is a dynamic obstacle which prevents the star from collapsing because of gravitational compression or exploding because of the added effect of the radiation pressures and of the hot gases.Ultimately, the role played by the magnetic hydrostatic pressure in this part of the theory is that of diminishing to a half the limit value of k.
This fact is important since it indicates that in general gaseous stars must be more luminous than it is believed and therefore that their age can be less than that which has been determined by Eddington's theory [12].In specialized literature [12] it is stated that the age of stars is inversely proportional to luminosity [12].Since luminosity is inversely proportional to opacity [1] as it will be seen further on, if the value of the Eddington limit of the opacity is in fact smaller than the one calculated so far, luminosity will be greater and therefore the age of the stars will be smaller.
Incidentally, an increase in luminosity is indicative of the fact that the star will live less since it is burning its fuel at a greater rate.With this result the apparent paradox of the stars which are older than the Universe itself, can be solved.In conclusion and within the present theoretic scheme, it can be stated that the age of gaseous stars is ruled by the modified Eddington limit of the coefficient of opacity.

The Magnetic Field on the Surface of a Star
On the surface of any gaseous star the magnitude of the self-generated magnetic field follows a very simple law as we shall see next.Its absolute value can be determined from the condition of magneto mechanical equilibrium (1) and considering that ρ = ρ(R, t), with R the star's radius.When that relation is integrated it is easy to obtain the following result [5] 2 constant 8π On the other hand, it is known that the pressure of hot gases ( ) satisfies the hydrostatic equation [5] [13], i.e.
( ) In that case and with the help of the relation (9) it can be proved that The preceding formula must be considered as a relation of equivalence more than an equality.From it and with the help of the expression for the pressure of the radiation, the absolute value of the magnetic field can be calculated.Since for a perfect gas it is fulfilled that where a = 7.64 × 10 −15 is Stefan's constant [1], of (18) it can be proved that with a universal constant that has the following numerical value 7 2

10 gauss K m
The relation (20) allows calculating the absolute value of the average magnetic field on the surface of a any gaseous star if its effective temperature T e is known.The effective temperature can be calculated from the radiation theory of the black body [1] [14].If E is the energy density radiated by the black body, it can be shown that The former result is known as the name of Stefan's law [1].Since for a perfect gas E = 3p r [1], the relation ( 19) is satisfied.On the other hand, the radiation emitted by a radius sphere r every second is [1] ( ) Therefore, the effective temperature of a star is defined by the following relation [1] [10] since it gives the black body temperature which produces the same amount of radiation emitted by the star [1].In (24) and (25) L is the luminosity.It is important to clarify that the effective temperature is a conventional measurement that specifies the reason of radiant flux heat per unit area.It must not be considered as the temperature at any particularly significant level in the star [1] [10] [12].Finally, it is known that the temperature of the photosphere of the Sun is equal 5.741 × 10 3 K.Consequently H  = 8.4 gauss on its surface.Another example is that of the most brilliant component of Capella's binary system for which it is known that its effective temperature is equal to 5.2 × 10 3 K [1].In that case H c = 6.8 gauss.

The Mass-Luminosity Relation and the Opacity Coefficient
Let us suppose that in the Equation ( 11) the product ηk is a constant through the star.In order to support that hypothesis it is required that the absorption coefficient k be practically constant decreasing a little towards the center to counterbalance the increase of η, in that way assuring that the product of these two quantities remains constant [1] [10] [12].It is said that there are very good reasons to believe that in general k behaves in this manner [1], wherefore it is possible to state that the constancy of this product is a good approximation [1].Be therefore with k o a constant that somehow represents the boundary value of k [1].It is important to emphasize that the value of k in the stellar photosphere could be very different from the one that it would have there k o [1].Using the approximation (26) in (11) and integrating, the following result is obtained that with the help of the relation ( 18) is transformed into what follows Clearly, it fulfills the fact that for a particular star the ratio p r /p g always maintains the same relation.Consequently At this moment it is important to introduce a constant β defined as follows [1] ( ) where p is whole pressure.β represents the ratio between hot gas pressure and whole pressure, whilst (1 − β) is the ratio between radiation pressure and whole pressure.Therefore, the ratio (1 − β)/β can be considered as a measurement of the degree of stability of a star.Substituting (30) in (29) we have that For gaseous stars the value of β is determined from the following quadric equation [1] ( ) where the symbol  indicates the mass of the Sun and µ the average molecular weight [1] [10].The relation (32) can be solved for various values of mass M in terms of the sun's mass, e.g.M = 1/4, 1/2, 1,... times ; and also for some average value of µ.It is important to point out that β only depends on the mass and average molecular weight of the material forming the star and it is independent of its radius and its opacity [1].Once obtained the value of β by solving the Equation (32) it is possible to determine the opacity k o from the relation (31).In specialized literature [1] there are tables which can be consulted to use in numerical calculations which contain values from (1 − β) for different masses and molecular weights.

Luminosity and Opacity
If the possible small changes of µ with the temperature and density [1] are neglected, it is possible to see from (32) that (1 − β) is a function only of the mass of the star.In that case and according to (31) it can be stated that for gaseous stars of the same mass, the luminosity is inversely proportional to the opacity k o [1].However the approximate constancy of k o from star to star should be distinguished from the approximate constancy of the product ηk within a single star [1].According to theoretical analysis carried out by other researchers [1] [10] [12], it can be stated that for a particular star the following relation is satisfied where k c refers to density and temperature conditions in the center of the star.On the other hand, the law of absorption of the radiation commonly accepted [1] [10] [12] has the following analytical form with k 1 is a constant.In that case it is customary to write that [1] whereα is a constant that according to (33) is equal to The constancy of ηT −1/2 depends on the ratio of energy release; release which is likewise related to the distribution of energy sources in the inner part of the star.What it normally does is to consider various degrees of concentration of those sources and to examine how the ratio of energy release per gram, ε, for different temperature powers varies [1] [10].However, even when the constancy of that factor is better for ε∼T, it is customary to explore other possibilities [1].For cases in which ηT -1/2 is not absolutely constant, the best approximation to the ratio k o /k c may be obtained from the average.From specialized literature [1] we have that 2 4 1.32 1.74 2.12 2.75 respectively.It is said that when the sources are very concentrated in the center of the star brightness diminishes.
In the case in which numerical calculations are required, the value α = 2.5 is adopted, which is indicative of an intense and uniform concentration of sources of energy release in the center of the star [1].Consequently in (31) the modified mass-luminosity relation is obtained When numerical calculations are made it is customary to assume that k c = k [1].

The Central Temperature
In order to calculate the central temperature in any gaseous star, relations ( 16) and ( 17) are used as well as the fact that the relation of equivalence (18) is always fulfilled as much on the surface as in any inner region of the star.Since for gaseous stars the thermal equation of ideal gas is fulfilled where and µ the molecular weight which is generally taken as a constant numerically equal to 2.2 [1] [10] [12], of (9) and from relations (30) the following result is obtained The previous relation is equivalent to the equation where results (38) and (40) were used.
In general, ratio T 3 /ρ which is a constant through any star is also the same for all the stars that have the same mass; whenever it is possible small differences in the average molecular weight µ that may exist among them, should be ignored [1] [10] [12].Consequently in stars of the same mass, the temperature at homologous points in the interior varies as the cube root of mean density measured at these points does [1] [10].As it is easy to see from paragraph 2, effective temperature is subject to a different law.

The Problem of Variable Stars of the Cepheid Type
The theory commonly accepted about the variable stars of the Cepheid type [1] [10] attributes the variation in their brightness to a regular pulsation they experiment.
In order to make an adequate theoretical analysis of the problem, it will be assumed that they are gaseous objects where the self-generated magnetic field has lost much of its original intensity but, that in the position of maximum compression keeps enough of it, so as to reduce the oscillation preventing the final collapse.Next, and with the help of radiation pressure and that of hot gases, that diminished magnetic field bounces, starting the subsequent expansion.If it is admitted that the Cepheids are in a stage of their evolution such that the fluid that forms them is no longer totally supported by the superstructure formed by the magnetic lines of force, it can be stated that the important dynamic agents acting on them are the huge gravitational force as well as the combined pressures of radiation and hot gases.Subject to these dynamic conditions, the Cepheids variables oscillate around some equilibrium position losing and recovering brightness alternately, as they expand and contract themselves with a noticeable regularity and with a perfectly determined period.
In accordance to the laws of gases, when the star collapses because of the effect of its huge weight and reaching even an extreme position or minimum size, the gaseous fluid heats up and the brightness of the Cepheid is increased.Right after that the hot gases, the radiation pressure and the residual magnetic field which behaves as a spring, act against gravity distending the star as far as another extreme position of maximum amplitude, causing the stellar fluid to cool and the Cepheid to lose brightness.
The cycle is repeated once and once again with regularity very much alike that observed in the movement of a simple harmonic oscillator [15].To obtain the differential equation governing the phenomenon and also an expression for the period of harmonic oscillation, Equations ( 11), ( 18) and (26) must be considered in order to obtain the following relation: Integrating the above equation the following is obtained ( ) where o r p and o g p are integration constants with unities of force per unit area.Besides p g satisfied the hy- drostatic Equation (17).Let us now consider a radius sphere r concentric with the star and with its surface almost coinciding with that of the star in such a way that the amount of mass M(r) contained in it is practically equal to the total mass M of the star.Under these conditions, the gravity acceleration on the surface of the star is It is proposed that as much on the surface of the sphere as on that of the star, o g p and the constant in (17) be equal to zero.However, there is 0 o r p ≠ because the pressure of radiation does not disappear on the surface [1].According to the first of the definitions (30), in (43) the following result is obtained ( ) where p is whole pressure such that 2 4π with f the magnitude of total force, 4πR 2 the area of the surface of the star and R its radius.Be in such way that in (44) the following is obtained .
In r = R we have that ( ) According to Newton's second law, of (47) the following result is obtained [15] . The above relation is the differential equation that governs the harmonic oscillations observed on the Cepheid variable stars.In it, ( ) is a constant which depends on some basic parameters of the particular star being studied.The oscillations period is given by the relation τ = 2π(M/K) 1/2 , namely ( ) where the relation (35) was used.Since M and L are total mass and luminosity respectively, in (50) average density ρ m must be used.From the mass-luminosity relation (38) and given that 4πG = 8.3 × 10 −7 ergs•cm•gr −2 , it is easy to see that ( ) ( ) As it was to be expected [3] [16], the product of the period in days and the square root of the mean density is equal to a constant that can be calculated from the theory.Next the periods of three known Cepheids are calculated and each result is compared to the one measured by direct observation.

δ Cepheid
For this intrinsic variable there are the following data (See Appendix) The period directly measured is of 5.366 days.As it is easy to see, theoretical calculation and direct measuring are practically equal.

Polaris
For this case we have that [1]  The period measured for this variable is equal to 3.968 days.

β Cepheid
The data we have are as follows [1]  in the two last mentioned cases are due to the fact that there are not more accurate data for these stars.However, the magnitude order is adequate.The long-period variables have, with some exceptions, periods ranging from 100 to 500 days, with a strong preference for periods near 300 days.The periodicity in not perfect and the stars may be some weeks behind or ahead of the predicted phase, the amplitude of the variation is not always the same in successive periods.Nevertheless, there is enough evidence that long-period variation and Cepheid variation are essentially the same phenomenon.The very low mass density and temperatures of the long-period variables exaggerates and renders more erratic the effects of the same kinds of pulsation as in the Cepheids [1].In general, the longer the period of a red variable, the less regular is its pulsation; as in the cases of Betelgeuse and Mira.So, for that class of variables the period-mass density relation is ( ) This is so, because if 1 sec = 0.1157 × 10 −4 days, and the period of the red variable star which will be use as the unit measure is equal to 3 × 10 2 days, it has that Hence from (50), 6.87 × 10 3 × 0.347 × 10 −2 = 24.Then, and given that the long-period variation and Cepheid variation are the same phenomenon, it is fulfilled that in such o way that 0.24 = 0.08 × 3. It is easy to see that is enough to multiply the numerical factor of the relation (51) by 3, in orden to obtain (52).Next, the results for the red variable Betelgeuse are given: Omicron Ceti, also called Mira, is a red variable star which has a period equal to 3.3 × 10 2 days, that is to say, its period is less that the one of Betelgeuse and hence, its average mass density must be greater than the Betelgeuse's average mass density.In consequence and given that τ o = 0.9 τ B , it follows that ρ mo = 1.11ρ mB ; and then, for the case of Omicron Ceti the following results are obtained.

The Magnetic Field in the Inner Part of a Gaseous Star
Just as temperature follows different laws on the surface and in the inner part of a star, something similar happens to the self-generated magnetic field.On the surface of the star its behavior is ruled by the relation (20) whereas in the inner part of it, it follows another different law as it will next be seen.In order to calculate the magnitude of the magnetic field at the center of the star as well as at any other inner point of it, it is necessary to use the polytropic gas sphere theory [1] [10].In terms of the gravitational potential φ(x), the acceleration of gravity is by definition [1]  Now we have the relations (54) and (55) in order to be able to determine the three following unknown functions of x: P, ρ y φ.Then, a third equation is also required in order to take into account the thermodynamic state of the star.In general and regardless of whether the stellar gas is perfect or not, it is always possible to make any value of P correspond to a given mass density if temperature is adequately fixed.What is usual [1] is to use as a third relation the following whereκ and γ are disposable constants.Thus and for different values of κ and γ it is possible to investigate a va-riety of temperature distributions.It is said that distribution is polytropic [1] if it obeys an equation such as (56).The problem is reduced to redoing the analysis carried out by Eddington [1], task that will not be repeated in this paper.However, and following that methodology, we have that the whole pressure p in terms of gravitational potential has the following ( ) where n is a positive integer number.Given that the thermal equation of the state of the ideal gas (39) is fulfilled, from (57) we have that ( ) As it is easy to see, T = constant φ.Be φ o the gravitational potential in the center of the sphere and T o the central temperature so that T o = constantφ o ; where the constant is the same as before.It can be seen that according to (58) and with the results of specialized literature [1], it is easy to see that for a particular star ( ) where M' y R' are parameters calculated from the polytropic gas sphere theory [1] and The relation ( 59) is used to calculate the central temperature of any gaseous star for which its radius and mass are known [1].Apart from this and according to the first of the relations (30), from (57) the following result is obtained For gaseous stars the relation (56) takes the following form [1] 4 3 in such a way that n = 3 [1].With this result and with the relation of equivalence (18), it is easy to see from (58) that for a particular star where Consequently, for stars of the same mass as well as for homologous points in the inner part of a given star, the self-generated magnetic field varies like the square root of the product of the mass density and the gravitational potential; both calculated at those points.
For any inner region of the star, the magnitude of the self-generated magnetic field can be calculated from (63), namely It is customary for numerical calculations to use the values given in Table 1 [1].
Usually the problem is reduced to finding the internal distribution for the density as well as for the pressure, in a star for which its mass and its radius or its mean density are known.The expressions used are the following [1].
The condition u = 0 indicates the boundary of the star, whereas the numerical values of the parameters R' and M' may be consulted in the last line of the table.Another important relation is the one for the ratio between the mean and central densities Its numerical value can be found at the bottom of the sixth column.On the other hand and since φ/φ o = u, it is where T o is the central temperature which can be calculated from the relation (58) for a constant value of the average molecular weight µ equal to 2.2 [1].The numerical value of u can be looked up on the second column of the table for different points within the star.For density we have the following relation [1] With ρ is the density at some inner region of the star.It is calculated from ρ o and using the numerical values of u n for different points within the star which are recorded on the third column.In order to show how the formulas above and the table should be used, consider as the first example that of Capella's brightest component.
Amongst the many data known of her we have that M = 8.3 × 10 33 gr and R = 9.55 × 10 11 cm.Its mean density is ρ m = 0.00227 and in column sixth it can be seen that the ratio ρ o /ρ m is equal to 54.36; so that in the center of the According to (60) the gravitational potential in the center is [ On the other hand it is known that [1] 1 0.283 0.717; in such way that the magnitude of the self-generated magnetic field in the center of the star is, according to (65) ( ) Consider another point inside the star, like for example the line z = 3.5 of the table.From (66) it is easy to see that 3.5 0.507 6.901 The point considered is found situated a little further from half the center of the star.From column sixth of the table we have that In other words, with that choice 90%of the star mass is being taken into account.The gravitational potential at that point is  The magnetic field at that point has the following absolute value ( ) which means that the magnitude of the magnetic field at the middle part of the star has been reduced to 11% of the value that it has in the center.Since due to boundary condition φ = 0 is always taken on the surface, in that region the magnetic field must be calculated with the Formula (20).That way it is obtained that on the 6.8gaussIn that position the gravitational potential has the following value  As it was expected, the magnitude of solar magnetic field diminishes in the same proportion as Capella's magnetic field.Indeed, its value is again equal to the 11% of the one it has in the center.Consequently the corresponding graphic is equal to that of the previous example, as it is easy to see in Figure 2.
The behavior of the magnetic field is in general very similar to the behavior observed in the density of the stellar mass in A.S. Eddington's theory [1] [4].
From the theoretical scheme developed in the above paragraph, it is possible to derive an expression to calculate the average value of the residual magnetic field which acts in the Cepheid type variable stars, in the maximum gravitational compression as well as in that of maximum expansion.If the absolute value of total force is considered, of the equations ( 45) and (47) as well as the mass-luminosity relation (38), it is possible to demonstrate that in r = R ( ) ( ) ( ) where the bar indicates the average values of those quantities; H is the residual magnetic field absolute value, ρ the mass density and R the star radius.Besides, the first of the definitions (30) was used.On the other hand, from the result (51) it is easy to see that ( ) ( ) in such way that in (70) the following is obtained where ( ) is a constant which depends on the basic parameters of the particular star being studied.Therefore, it is possible to state that for each particular Cepheid, the residual magnetic field is inversely proportional to the product of the period and the square root of the radius; so that for the situation of maximum compression we have that whereas for the maximum expansion The maximum and minimum values of the radius reached by the star in its extreme volumes can be calculated from the relation (25) and from the observational data obtained for the luminosity and the effective temperature at the position of maximum and minimum amplitude of the harmonic oscillations carried out by the particular Cepheid being observed.

Conclusions
In the present paper, a fundamental hypothesis is made which consists of assuming that all gaseous stars at some early stage of their evolution, and in a natural way, generate an intense magnetic field whose magnitude diminishes from the center of the star to the surface of it in a very peculiar manner.This field generates in the conducting stellar fluid, a magnetic viscosity which gives it a great rigidity that can be interpreted as a rigid superstructure formed by the magnetic lines of force, whose mission is to keep the shape of the star.
That internal magnetic superstructure prevents the star from collapsing gravitationally or exploding because of the added effect of radiation pressures and those of hot gases.In some way, the magnetic lines of force of the self-generated magnetic field behave in the inner part of the star as a powerful shock absorber acting in both ways: against gravitational compression as well as against hot gases and radiation expansion.
When the magnetic field is introduced in the theory in the form of a magnetic hydrostatic pressure [5]- [7], the results previously obtained by Eddington are slightly but substantially modified [1], particularly, it is reduced to half the value of Eddington limit for k [1].
With the help of this result, it is easy to point out the important role the modified limit has in the problem of the stars age and at the same time to indicate a possible solution to the paradox consisting of the existence of stars that are apparently older than the Universe itself.The present theoretical frame derives a simple formula to calculate the magnitude of the average magnetic field on the surface of any gaseous star for which its effective temperature is known.
Since the behavior of the self-generated magnetic field follows a different law in the inner part of the star, a formula is deduced which shows that this field varies like the square root of the product of mass density and the gravitational potential, in both, the center of the star as well as in inner homologous points.Numerical calculations are carried out and the results are presented and then drawn for two specific cases: that of Capella's brightest component and that of the sun.
As it is easy to see, the graphs are exactly alike; which shows that the magnetic field self-generated by the gaseous stars, has the same general behavior for all of them.
Another interesting result refers to the solution of the problem of the Cepheid type variable stars.For those cases, it is proposed that the magnetic field is much weakened due to the fact that those stars are in their last evolutionary stage.However, it is not so much so as not to be able to stop the gravitational collapse acting as a spring which is compressed when the star reaches a certain minimal size.From that minimal volume, the residual magnetic field and the sum of radiation pressures and those of hot gases, inflate the star up to a certain maximum volume from which the compression-expansion cycle is again initiated.
Given the regularity of the oscillations, it is assumed that the process is simple harmonic with a perfectly determined period.Integrating the fundamental result of the theory, which is the relation (11), an ordinary differential equation of simple harmonic oscillator and the expression for the oscillations period in terms of the basic parameters of the star are obtained.In perfect agreement with the observations carried out by other researchers [3] [16] and from the mass-luminosity relation (38), it is shown that the product of the period and the square root of the mean density is equal to a constant.
The numerical value of the constant depends on the star being studied, and it is easily obtained from the theory.It is important to make clear that the above relation is empirically proposed and recorded in the specialized literature as a fact perfectly well established in the different models that on the stars internal structure have been proposed [1] [3] [16].
Suppose that the gravitational collapse suffered by a star that is in the last stage of its evolution is so fast and vigorous that the magnetic field noticeably weakened does not have any more an intensity powerful enough so as to stop it initiating the rebound.Under these conditions, the star collapses, and due to the combined effect of material and radiation pressures as well as the sudden hydrogen combustion still present in large amounts in the regions of the star atmosphere, a huge explosion is produced which can send into the space up to 9/10 parts of the star total mass.In that case, it is said that a nova or a supernova has risen depending on the magnitude of the sidereal catastrophe.Whether one thing or another should occur possibly depends on how massive the star is.
The final outcome of the catastrophe will also depend on that last condition.It is possible to have a recurrent nova, a white dwarf, a pulsar, or even a black hole.Next, to calculate its mass, the mass-luminosity relation (38) will be used as follows where the relations (34), ( 40) and (59) were used.From the equation [1]

R
a constant, mass M can be eliminated.In that case, it is easy to see that ( )

R
, in such way that ( ) With constant values we have that ( ) For k 1 the value found for the brightest component of Capella's binary system will be used [1] and it will be taken for the average molecular weight µ the usual numerical value 2.11 [1].Since β is the only unknown quantity in the above relation; the equation obtained finally is ( )

=× ⋅
To conclude, it is easy to prove that using the relation of equivalence (18) and the formulas (30) and (26), the following equation is obtained ( ) by means of which it is also possible to calculate the magnitude of the self-generated magnetic field, in both, the center as well as the inner homologous points of any gaseous star.
From the results of paragraph 7 and with the help of numerical data from the table, from relations (59) and (67) it is obtained for the Capella case that In such away that in the center of the star, the magnitude of the self-generated magnetic field is again From the second column of the table we have that for z = 3.5, u = 0.27629 so that at this point  It is easy to see that for the calculated values of H and H c with this last method, it is also fulfilled that in both examples the ratio H/H c always keeps the same proportion.In other words, it can be proved that ( ) The numerical values of u n+1 can be looked up on the fourth column of the table for different points within the star.
From the Equations ( 20), ( 22), (25) and from the mass-luminosity relation (38), the following expression to calculate the mass for any gaseous star, using some of their basic parameters is obtained ( ) With the numerical values of the constants, the next result is obtained 1.12 10 gr cm gauss 4 Consequently, for a particular star is fulfilled that ( ) where 6 constant 1.12 10 In other words, the mass of a gaseous star is proportional to the square of the product of its radio and the magnitude that the self-generated magnetic field has on its surface.From the data we have from the Sun it is easy to verify by direct calculations that M  = 2.0 × 10 33 gr; whereas for Capella we have that M = 8.3 × 10 33 gr.
where p(x, t) is the total pressure, ρ(x, t) the mass density and η, ζ the coefficients of viscosity [6] [9].The relation (B-2) is the momentum balance equation of magnetohidro dynamics (MHD).For any gaseous star, the mass density is in general, a function of time t and the radius R * , such that, ρ * = ρ(R * , t).
It can be directly shown that (B-2) takes the following equivalent form ( ) ( ) where ( ) is the viscosity stress tensor [6] [9], the relation of equivalence (18) and the following well known formula in vector analysis H rot H were used; after integration by parts and taking into account the condition (B-1).Suppose that the star revolves in such a way that the flux of the stellar fluid is steady, and then 0 t * ∂ ∂ = v . In this case, it is fulfilled that ( ) ( ) Using this last result in (B-3) and with the help of Equations ( 9), ( 19) and (39), it can be directly obtained that . This last result is valid for any gaseous star with radius R * , that revolves with a rotational velocity ( ) . On the other hand is easy to see that in the inner regions of a given star, the first term is greatest than the second one, and therefore that term can be neglected.Consequently and for any inner region of the star, the magnitude of the magnetic field can be calculated from the following equation At the surface of the star, the mass density is practically zero, and then, the magnitude of the magnetic field can be obtained from the next relation  20) of the text.Let's, consider that H * is the magnitude of the self-generated magnetic field by gaseous stars.Then, it is easy to see that in the inner regions of gaseous stars, the self-generated magnetic field varies like the square root of the product of the mass density and the absolute temperature; both calculated at those regions.On the other hand, at the surface, it depends on its effective temperature; according to different laws perfectly established.
In the theoretical frame of MHD and because the magnetic permeability of the media differs only slightly from unity, and the difference is unimportant [6] [9], it is consider that, * = H B and then, magnetic field fulfill the basic laws of magneto statics; that in their differential form are the condition (B-1) and the following relation where j is the steady-state current distribution localized in some region of the convective zone and c is the velocity of light in empty space.According to (B-1), ( ) * B x must be the curl of some vector field ( ) * A x , called the vector potential [7]; that is, For a steady-state current distribution localized in a relative small region of space, the vector potential is given by the following expression [7] ( ) ( ) where ′ x is a distance measured relative to a suitable origin in the localized current distribution [7].Starting with (B-8) it expands the denominator in powers of ′ x , until the lowest order of approximation.Then for a given component of the vector potential will have the next expansion [7] ( ) ( ) ( ) This can be established by an integration by parts of the second term, followed by expansion of ( ) With f = 1, and i g x′ = , the last expression establishes that ( ) The first term in (B-9) corresponding to the monopole term in the electrostatic expansion is therefore absent.The integral in the second term of (B-9) can therefore be written in the following form ( ) It is customary to define the magnetic moment density or magnetization as [7] ( ) ( ) and its integral as the magnetic moment m; that is Then, the vector potential from the second term in (B-9) is the magnetic dipole vector potential ( ) This is the lowest nonvanishing term in the expansion of * A for a localized steady-state current distribution.The corresponding magnetic induction * B can be calculated directly by evaluating the curl of the last equation [7]; that is ( ) where n is a unit vector in the direction x .The magnetic induction * B has exactly the form of the field of adipole.Faraway from any localized steady-state current distribution, the magnetic induction is that of a magnetic dipole of dipole momentum given by (B-11).Additionally, from Equation (B-2) it can be shown that where T is the total torque and only the magnetic part of the force was used and taking into account the Equation (B-6).From this last result and writing out the triple vector product, it can be directly obtained that ( ) ( ) The second integral vanishes for a localized steady-state current distribution, as can be seen from equation before (B-9) with f = g = x' [12].The first integral is therefore [12] This is the familiar expression for the torque on a dipole; it's one of the ways of defining the magnitude and direction of the magnetic induction * B [7].Consequently, the self-generated magnetic field of gaseous stars, is produced by some special kind of mechanism.In fact, according to density and temperature conditions, it can exist some region into the convective zone that has a maximum of ionization.
The electrically charged particles are moved by the convective streams across that region, making its contribution to the steady-state current distribution localized and goon; and they are continuously replaced by other particles.
That is an effect likely in form to the one that produce the rainbow.The drops of water moved across the region where the rainbow is produced, making its contribution to the phenomenon and continuing with its fall; and they are continuously replaced by other drops.
It is for that effect, that the rainbow is a steady-state and localized optical phenomenon; like is too, the proposed steady-state current distribution localized.
Since this current distribution is produced by the high ionization of the region and the process of ionization depends on density and temperature conditions of the region, the self-generated magnetic field of gaseous stars is a function of these variables, as is easy to see form Equation (B-5).
To conclude, it is important to mention what follows: concerning the elaboration of an alternative theoretical scheme on the origin and structure of the magnetic field self-generated by gaseous stars, many researches have engaged themselves, throughout 40 or more years, to the self-exited dynamo models [19].
Unfortunately, the theoretical results obtained by then are far from satisfactory even now [20].The model was initially proposed in 1919 by J. Larmor [21] with the purpose of giving an explanation to the phenomenon of Sunspots.That suggestion was quickly rejected for being inadequate and inconsistent to the astronomical observations about the phenomenon [22].
Also, by means of this model, the idea is to explain the origin and structure of the magnetic field self-generated by the Earth.Neither in this case it has been possible to give any explanation which is satisfactory concerning the basic characteristic of the above mentioned magnetic field [23].
τ (measured) = 0.190 days.It is possible that discrepancies observed P = p − H 2 /8π, in such a way that in (1) the following result is obtained ∇ = grad the laplacian operator[9], we have that for spherical symmetry the above relation takes the following form, from the third column of the table that With these three points it is possible to build Figure1showing the general behaviour of the field.Consider the case of the Sun, for which we have the following of the Sun the gravitational potential has the following value
H to H o we have that for the case of the Sun 0.11 o H H = .

3 0
Since it is also known that ρ m =0.00227, we have that As second example consider the case of the Sun for which we have the following data hand, from the absolute bolometric magnitude determined for the intrinsic variable δ Cepheid [ Using those values in the quadric Equation (32) it is easily obtained that of the intrinsic variable δ Cepheid has the following value It can be proven that in the point z = 3.5, the coordinate of a point at a great distance of the localized current distribution.The fact that j is a localized divergenceless current distribution, allows simplification and transformation of expansion (B-9).Let ( )f ′ x and ( ) g ′x be well-behaved functions of ′x .Then, if ( ) ′ j x is localized and has zero divergence, it has[7]