Sep 21, 2012 - 1, increasing rs (increasing the plasma density) from zero-density ... Figure 3(b) within a fixed relativistic density heavier atoms crystallize into.
Phase Transition in Relativistically Degenerate Quantum Plasmas M. Akbari-Moghanjoughi
arXiv:1209.4742v1 [physics.plasm-ph] 21 Sep 2012
Azarbajan University of Shahid-Madani, Faculty of Sciences, Department of Physics, 51745-406 Tabriz, Iran In this paper we present the generalized phase separation diagram for quantum plasmas revealing the fundamental mechanism called the spontaneous core collapse (SCC) for a white dwarf star coredensity based on the extended Shukla-Eliasson force for relativistically degenerate matter taking into account the most important electron interaction features such as Fermi-Dirac statistical pressure, Coulomb, exchange-correlation and quantum tunneling effects. A critical core mass-density value is reported for white dwarf stars beyond which the SCC mechanism will be inevitable. It is revealed that the SCC critical core-density value is slightly affected by the ion charge state, Z (the core composition). Such fundamental effect which is mainly caused by the electron quantum diffraction effect can progressively lead the white dwarf star with a mass beyond the Chandrasekhar limit through SCC, core fusion and the ultimate supernova stages. PACS numbers: 52.30.-q,71.10.Ca, 05.30.-d
2 INTRODUCTION
During the past decade, a new dimension of research has initiated due to Bohmian collective interpretation of quantum systems [1–3]. Foundations of quantum hydrodynamic (QHD) formulation, earlier suggested by Madelung [4] termed as hydrodynamic formalism for quantum mechanics, was extended by Takabayasi [5–9] and others. It has been shown that in quantum system of Fermionic particles a new force due to the interaction of particle wavefunction may result in important collective plasma dynamic features. Such ne collective feature is called the Bohm quantum force [10, 11] and is the main ingredients of the quantum hydrodynamics (QHD) formulations. In dense plasmas where the inter-electron distances are comparable to or lower than the de Broglie thermal-wavelength λD = h/(2πme kB T )1/2 , quantum effects such as the Pauli exclusion [12] and the electron quantum Bohm force will be inevitable. Quantum hydrodynamic theory has been successfully employed to interprets the collective dynamics of electron-hole plasmas in a semiconductor [13–15]. Variety of applications may be foreseen for quantum plasmas other than semiconductors [16]. For instance, the physical characterization of nanofabricated devices fundamentally rely on the quantum theory of hydrodynamics. A large number of recent investigations based on QHD model confirm signicant role of quantum diffraction effect on linear and nonlinear wave dynamics of degenerate plasmas [17–22]. It has been shown that the inclusion of Bohm force can lead to many critical features [23, 24]. More recently Shukla and Eliasson by investigating the screening effects in a quantum plasma have discovered a new attractive force between ions [25] which is due to interplay between electron degeneracy, exchange-correlation and quantum diffraction effects. The existence of such attractive force among ions, if established experimentally, can mark a distinguished milestone in the history of plasma physics research with the Bohmian collective interpretation of the quantum particles. Moreover, the Wigner-Poisson formulation of hydrodynamics has been recently extended to include the electron spin magnetization effects [26, 27] or the relativistic electron momentum effects [28]. Using the quantum magnetohydrodynamic (QMHD) theory, it has been confirmed that the electron spin can play signicant role in linear and nonlinear wave dynamics of magnetized quantum plasmas [29–34]. One of the important applications of QHD theory regards the astrophysical applications of degenerate plasmas which are ubiquitous in the nature [35–40]. Chandrasekhar [41, 42] has shown that the compact stars called white dwarfs can be effectively modeled as the ideal zero-temperature quantum plasmas. It has also been shown that the large external gravitational force on stars can lead to a phenomenon known as the relativistic degeneracy [43] leading to distinct features in equation of state (EoS) of degenerate matter [44] and nonlinear wave dynamics [45] in different density regimes. However, the electron diffraction mechanism is supposed to be affected by the relativistic plasma effects, since, the relativistic effects are directly coupled to the plasma number density in the relativistically degenerate matter. Mettes and Sorg [46] and Carlos and Mahecha [47–50] have extended the fluid dynamics to include the generalized relativistic Bohm force by solving the time dependent Dirac or Klein-Gordon equations. More recently, a relativistic fluid model based on Wigner function and Dirac equation has been proposed [51] which includes the relativistic corrections to quantum diffraction effect. In the following, we consider a zero-temperature relativistically degenerate collisionless electron/ion fluid with nondegenerate ions of charge state, Z. The relativistic quantum hydrodynamic fluid model incorporating the relativistic electron motion, quantum statistical pressure, interaction effects such as Coulomb, exchange-correlation, etc. and the relativistic quantum Bohm force [51], may be written as
∂γn = 0, ∂t + ∇ · (γnu) ∂u γ ∂t + (u · ∇)u √ ∆√ n h ¯2 , = mee ∇φ − m1e n ∇PG + 2γm 2∇ n e
(1)
∆φ(r) = 4πe(γn − n0 ),
where, PG is the generalized ppressure due to relativistic electron degeneracy in addition to other electron interaction pressures [52]. Also, γ = 1/ 1 − β 2 is the relativistic factor with β = u/c. It should be noted that the additional correction to the relativistic Bohm force given in Ref. [51] is ignored for being relatively small proportional to β/c. Using the concept of generalized effective potential, same as in Ref. [53], one may obtain the plasma dielectric response to harmonic perturbations via the Fourier analysis of Eqs. (1) as follows
3
−ωn1 + k · u1 = 0, h ¯2 2 2 −ωu1 = − 4m 2 k kn1 − Lc kn1 , e h R2 a 0 R0 L = √ 0 2 − β3 R0 + 2π 3 − 3
1+R0
R0 1+R20
+
R20 sinh−1 R0 (1+R20 )3/2
i
(2) ,
It is evident that the relativistic parameter, γ, does not enter the linear response model. In the linearized model, Eq. (2), the parameter R0 = (n0 /nc )1/3 with nc = 5.9 × 1029 cm−3 is the√Chandrasekhar’s relativity parameter, a0 = e2 /¯hc ≃ 1/137 is the fine-structure constant and β = (2a0 /5)(2Z/3 π)2/3 . Therefore the static dielectric constant is then given as
2 D(0, k) = 1 + ωpe
¯ 2 k4 h + Lc2 4m2e
−1
,
(3)
p where ωpe = 4πn0 e2 /me is the electron plasmon frequency. The dielectric around a point-like stationary charge, Q = Ze, located at the origin is given by
φ(r) =
Q 2π 2
Z
exp(ik · r) 3 d k, k 2 D(0, k)
(4)
Decomposition of the integral with known solution given in [25] reads as
φ(r) =
Q (1 + b) e−k− r + (1 − b) e−k+ r . 2r
(5)
with the potential parameters defined as below 2 k∓ = k02 1∓
α=
√
1−4α , 2α
2 h ¯ 2 ωpe 4m2e c4 L2 ,
b=
k0 =
√ 1 − 4α,
ωpe √ . c L
(6)
One of the fundamental characteristics of the Shukla-Eliasson potential [25] given by Eq. (5) is that for α > 1/4 the potential is attractive and repulsive otherwise. The generalized degenerate matter phase separation then is obtained by evaluation of the parameter, α. In Fig. 1 we have shown the variation of this parameter for wide density parameter, rs = rB /r0 range, where, rB = ¯h2 /me e2 ≃ 5.3 × 10−8 cm is the Bohr-radius and r0 = (3n0 /4π)1/3 is the Wigner-Seitz cell radius defined as a sphere containing only a single electron. The relativity parameter is given as R0 ≃ 0.014rs . As it is observed form Fig. 1, increasing rs (increasing the plasma density) from zero-density first leads the pure Coulomb screening (α = +∞) to turn over to the Lennard-Jonse-like Shukla-Eliasson attractive potential of form φ(r) = (Ze/r) exp[−ℜ(k− )r]{cos[ℑ(k− )r] + b sin[ℑ(k− )r]}, where ℜ and ℑ refer to the real and imaginary parts of the parameter k− . For two different values of rs1 (Z) and rs2 (Z) we have α = 1/4. For rs1 (Z) < rs < rs2 (Z) the plasma undergoes the ordinary Debye-H¨ uckel screening and for rs > rs2 (Z) the Shukla-Eliasson attractive force reappears. It is remarked that, the value of rs2 (Z) unlike that of rs1 (Z) is not strongly sensitive to the ion charge state, Z.This value coincides exactly with the core-density of a typical white dwarf star (108 < ρc (gr/cm3 ) < 1013 ). It is remarked that rs1 = 6.8326 value for iron plasma composition extends to density of ρ ≃ 821.945gr/cm3 which is well beyond the zero-pressure metallic iron. This is the indication that the Shukla-Eliasson attractive potential leading to the plasma crystallization can be present at small-planet core density ranges. On the other hand, the rs2 ≃ 2354.8 value for the same plasma composition corresponds to the mass density of ρ ≃ 3.5 × 1010 gr/cm3 which resides in the range of a typical white dwarf core. The reappearance of Shukla-Eliasson attractive force in white-dwarf core density has some fundamental consequences for stellar evolution, as will be discussed next. Figure 2 reveals some fundamental differences of Shukla-Eliassom attraction between nonrelativistic and relativistic ¯ in degeneracy cases. Figures 2(a) and 2(b) compare the potential profiles (normalized binding energy, rB φ(r)φ/Ze) nonrelativisic and relativistic degeneracy cases against the variation in the ion number density, Z, where r¯ = r/rB . It is revealed that, the increase in the number-density leads to increase in the potential depth in both degeneracy regimes while this variation is more pronounced in the nonrelativistic degeneracy case. However, as it is remarked comparing
4 Z = 0, 2, 6, 12, 26 100
10
Α
1
0.1
0.01
0.1
1
10
100
1000
104
105
rs FIG. 1. The variation of α against the plasma density parameter, rs = aB /r0 . The critical value α = 1/4 is indicated with the horizontal line. The thickness of curves appropriately reflect the variation in the ion charge value, Z
Figs. 2(c) and 2(d) the increase in the plasma density causes reverse effects in the two density regimes. In other words, the attractive potential strength (binding energy) increases/decreases for the case of relativistic/nonrelativistic degeneracy regimes due to increase in the plasma density. The variation in the potential-minimum location also follows quite dissimilar rout in the two degeneracy regimes. Figure 3 closely inspects the variations of Shukla-Eliasson attractive potential parameters, namely, the binding energy and binding length for the relativistic degeneracy regime relevant to the white dwarf core density for three different core compositions of Helium, Carbon and Magnesium. It is observed that the unlike for the binding energy, ¯ rmin ), the binding length, r¯min , is unaffected by the change in plasma atomic-number, Z. It is further remarked φ(¯ that in the relativistic densities the binding length decreases suddenly with increase of the plasma density and reaches an asymptotic value of about 10−14 cm, which is which is two orders of magnitude smaller that the nucleon diameter which may seem controversial. However, One should note that the hydrodynamic theory used in this model may not even be valid beyond the quantum screening length. On the other hand, other processes like nuclear fusion takes place before such short ionic distances are met. Figure 3(b) within a fixed relativistic density heavier atoms crystallize into stronger lattices. Crystallization in white dwarf cores are sometimes wrongly attributed to the presence of strong gravitational inward pressure at the center of the star. However, such confinement of ions in the core may only lead to quasi-lattice formation with viscous ion fluid due to lowering of ion coupling parameter value in fusion-medium temperatures. Slattery et al. [54] have shown that for ion coupling parameter values of Γ < 178 the plasma will be in liquid state. Therefore, Shukla-Eliasson attractive force can attribute appropriately to the perfect crystallization in a wide density range relevant to both cores of small planets as well as that of a typical white dwarf star. Figure 4 reveals yet another fundamental feature of Shukla-Eliasson screening for relativistic density regimes. It depicts variation of the normalized binding strength versus the normalized binding length for different plasma compositions. It is observed that with decrease of the binding length the binding energy increases for all plasma composition with more pronounced increase for higher plasma number density. Such variations reveal that the plasma crystals with relativistic density are unstable under external compression. Any force exerted from outside such as the gravitational inward pull for the case of a compact star would decrease the binding length and lead to collapse of the configuration. Such progressive increase in the plasma binding energy due to the gravitational pressure would also result in the energy release in the core leading the nuclear fusion ignition. As the atoms fuse they will form heavier elements which are lead to form even more stronger bounds and the release of more energy. Such loop-hole may proceed infinitely until the ultimate collapse of the white dwarf core releasing enough energy to lead to a cataclysmic supernova to take place. It should be noted that such mechanism called the spontaneous core collapse is only possible in the relativistic degeneracy regime, rs > rs2 (Z) and is not possible for rs < rs1 (Z).
5 n0 = 1.3 ´ 1025 cm-3 , Z = 1, 2, 3, 4
n0 = 2 ´ 1035 cm-3 , Z = 1, 2, 3, 4
0.5
20
0.0
0
-0.5 -20
Φ HrL ´103
-1.0
Φ HrL ´103
-1.5
-40
-60 -2.0 -80
-2.5 -3.0
1.0
1.5
2.0
2.5
3.0
3.5
-100 0.004
4.0
0.006
0.008
0.010
0.012
r
r
HaL
HbL
n0 = H1.6, 5.4, 13, 25L ´ 1024 cm-3 , Z = 2
n0 = H2, 3.5, 5.5, 8.3L ´ 1035 cm-3 , Z = 2
1
0.014
0 0 -1
Φ HrL ´103
-100
-2
Φ HrL ´103 -200
-3 -4
-300 -5 -6
1
2
3
4
5
-400 0.002
0.004
0.006
0.008
r
r
HcL
Hd L
0.010
0.012
FIG. 2. (Color online) Variations of normalized electric potential φ¯ = rB φ/Ze for α > 1/4 with changes in plasma density and the ion charge-state for nonrelativistic (Figs. a and c) and relativistic (Figs. b and d) degeneracy density regimes. The distance r¯ = r/rB is measured in the Bohr-radius unit. The thickness of curves appropriately reflect the variation in the varied plasma parameter, in each plot.
We conclude that Shukla-Eliasson quantum screening potential which goes beyond the ordinary Debye or the Thomas-Fermi models can successfully explain the structure formation in solid-density matter, cores of small planets and compact stellar objects. Such theory which relies on the complete quantum interpretation of collective fermion behavior can even lead to deep understanding of the stellar formation and structure. It can also result in breakthroughs in pure as well rapidly growing nanotechnology sciences.
[1] [2] [3] [4] [5]
D. Bohm and D. Pines, Phys. Rev. 92 609(1953). D. Pines, Phys. Rev. 92 609(1953). P. Levine and O. V. Roos, Phys. Rev, 125 207(1962). E. Madelung, Z. Phys., 40 322(1926). T. Takabayasi, Prog. Theor. Phys., 8 143(1952).
6
Normalized lattice spacing, dΦHrLdr=0, for Z = 2, 6, 12
Normalized binding energy, ΦHr min L, for Z = 2, 6, 12
14
0
12 -20
10
r min ´10
8
Φ Hr minL
6
´10-3
6 4
-40
-60
2 -80
0 0
2
4
6
rs ´ 10
8
0
10
2
4
6
8
10
rs ´ 10-7
-7
HbL
HaL
FIG. 3. (Color online) Variation in the binding energy and length in the relativistic degeneracy regime for three different plasma compositions, namely, Helium, Carbon and Magnesium. The thickness of curves being a measure in the value of the varied parameter in each plot. Spontaneous Core Collapse HSCCL, Z = 2, 6, 12 0
-20
Φ Hr minL ´10-3
-40
-60
-80 0
2
4
6
8
10
12
r min´106 FIG. 4. (Color online) Variation of the normalized plasma binding energy in terms of the binding length for various plasma composition, namely, Helium, Carbon and Magnesium. The plot shows a feature specific to relativistic degeneracy regime called the spontaneous core collapse relevant to white dwarf stars. The thickness of curves being a measure of the value of the varied parameter.
[6] [7] [8] [9] [10] [11] [12] [13]
T. Takabayasi, Prog. Theor. Phys., 14 283(1955). T. Takabayasi and J. P. Vigier, Prog. Theor. Phys., 18 573(1957). T. Takabayasi, Prog. Theor. Phys., 9 187(1953). T. Takabayasi, Nuovo Cim, 7 118(1958). G. Manfredi, Fields Inst. Commun. 46 263(2005). P. K. Shukla, B. Eliasson, Nonlinear aspects of quantum plasma physics Phys. Usp. 51 53(2010). L. D. Landau and E. M. Lifshitz, Statistical Physics, Part I, Pergamon, Oxford, (1978). C. Gardner, SIAM, J. Appl. Math. 54 409(1994).
7 [14] P. A. Markowich, C. A. Ringhofer, and C. Schmeiser, Semiconductor Equations (Springer-Verlag, New York, 1990). [15] M. Bonitz, D. Semkat, A. Filinov, V. Golubnychyi, D. Kremp, D.O. Gericke, M.S. Murillo, V. Filinov, V. Fortov, W. Hoyer, J. Phys. A, 36 5921(2003). [16] H. Haug and S. W. Koch, Quantum theory of the optical and electronic properties of semi-conductors, World Scientic, (London)2004, [17] G. Manfredi and F. Haas, Phys. Rev. B 64, 075316 (2001). [18] F. Haas, L. G. Garcia, J. Goedert, and G. Manfredi, Phys. Plasmas 10, 3858(2003). [19] A. P. Misra and P. K. Shukla, Phys. Plasmas, 15 052105(2010). [20] A. P. Misra, G. Brodin, M. Marklund, and P. K. Shukla, Phys. Plasmas, 17 122306(2010). [21] Martin Stefan, Gert Brodin and Mattias Marklund, New J. Phys., 12 013006(2010). [22] A. Mushtaq, and S. V. Vladimirov, Phys. Plasmas, 17 102310(2010). [23] M. Akbari-Moghanjoughi, Phys. Plasmas, 17 052303(2010). [24] M. Akbari-Moghanjoughi, Phys. Plasmas, 17 082317(2010). [25] P. K. Shukla and B. Eliasson, Phys. Rev. Lett. 108, 165007 (2012); ibid. 108, 219902 (E) (2012); ibid. 109, 019901 (E) (2012). [26] J. Zamanian, G. Brodin and M. Marklund, New J. Phys., 11 072017(2009). [27] G. Brodin, M. Marklund, and G. Manfredi, Phys. Rev. Lett., 100 175001 (2008). [28] F. Asenjo, V. Munoz, J.A. Valdivia, S. H. Hahajan, Phys. Plas., 18 012107(2011). [29] M. Marklund and P. K. Shukla, Rev. Mod. Phys. 78, 591 (2006). [30] M. Marklund, G. Brodin, L. Stenflo, and C. S. Liu, Europhys. Lett., 84 17006 (2008). [31] G. Brodin and M. Marklund, New J. Phys., 9 277(2007). [32] G. Brodin and M. Marklund, Phys. Plasmas, 14 112107(2007). [33] Vitaly Bychkov, Mikhail Modestov, and Mattias Marklund, Phys. Plasmas, 17 112107(2010). [34] M. Akbari-Moghanjoughi, IEEE Trans. Plasma Sci., 40 1330(2012), Doi: 10.1109/TPS.2012.2189137 [35] H. R. Miller, P. J. Witta, Active Galetic Nuclei, Springer-Verlag, Berlin, (1987) pp. 202. [36] P. Goldreich, W. H. Julian, Astrophys. J., 157 869(1969). [37] F. C. Michel, Rev. Mod. Phys., 54 1(1982). [38] E. Tandberg-Hansen, A.G. Emshie, The Physics of Solar Flares, Cambridge Univ. Press, Cambridge,(1988) pp.124. [39] M. J. Rees, in: G. B. Gibbons, S. W. Hawking, S. Siklas (Eds.), The Very Early Universe, Cambridge Univ. Press, Cambridge, (1983). [40] W. Misner, K.S. Throne, J.A. Wheeler, Gravitation, Freeman, San Francisco, (1973) pp.763. [41] S. Chandrasekhar, An Introduction to the Study of Stellar Structure, Chicago, Ill. , (The University of Chicago press), (1939) pp.392. [42] S. Chandrasekhar, Mon. Not. R. Astron. Soc., 113 667(1953). [43] S. Chandrasekhar, Science, 226 4674(1984). [44] D. S. Kothari and B. N. Singh, Proceedings of the Royal Society of London. Series A, Math-ematical and Physical Sciences, Vol. 180, No. 983 (Jul. 3, 1942), pp. 414-423 [45] M. Akbari-Moghanjoughi, Phys. Plasmas, 17 092304(2010). [46] M. Mattes and M. Sorg, Nuvo Cim., 110 1323(1995). [47] C. Castro and J. Mahencha, Prog. Phys., 1 38(2006). [48] C. Castro, J. Math. Phys., 31 2633(1990). [49] C. Castro, Found. Phys. Lett., 4 81(1991). [50] C. Castro, Found. Phys., 4 569(1992). [51] J. Zhu and P. Ji, Phys. Rev. E, 81 036406(2010). [52] E. E. Salpeter, Astrophys. J. 134, 669 (1961). [53] M. Akbari-Mojhanjoughi, Shukla-Eliasson Attractive Force: Revisited, J. Plasmas Phys. 78, in press, (2012); DOI:http://dx.doi.org/S00223778120000839; E-print:arXiv:1207.5154v2 [astro-ph.EP]. [54] W. L. Slattery, G. D. Doolen, and H. E. DeWitt, Phys. Rev. A. 26, 2255 (1982).