1. Introduction

JHEPGC

Journal of High Energy Physics, Gravitation and Cosmology

2380-4327

Scientific Research Publishing

10.4236/jhepgc.2022.81002

JHEPGC-113318

Articles

Physics&Mathematics

Plasmas Created in the Interaction of Antiprotons with Atomic and Ionized Hydrogen Isotopes. Suggested Fuels for Space Engines

Mohamed

Assaad Abdel-Raouf

¹^*Abdelfattah

T. Elgendy

¹Amr

Abd Al-Rahman Youssef

Physics Department, Faculty of Science, Ain Shams University, Cairo, Egypt

Mathematics Department, Faculty of Science, Al-Azhar University, Cairo, Egypt

17112021

080114247, August 202120, November 2021 23, November 2021

2014

This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/

The main objective of the present work is to investigate the properties of plasmas created by injecting a thermalized beam of antiprotons in two types of media. The first is hydrogen, deuterium, or tritium atoms localized in palladium crystals. The second medium is composed of protons, deuterons, or tritons localized in a magnetic cavity. Particularly, it is demonstrated that huge amounts of energy are released in both cases which could be used as fuels for space shuttle engines. A novel mathematical scheme is employed to calculate the energy yields in real space at different incident energies of the antiprotons.

Antiprotons Antiprotonic Hydrogens Antiprotonic Plasma Fuel for Space Engines Plasmas in Molecular Crystals Palladium as Host for Neutral Plasma Antiprotons-Ionized Hydrogen Isotopes Plasma

1. Introduction

With the discovery of antiprotons (p) by Chamberlain et al. [1] at Fermilab in the fifties of the preceding century, the ultimate proof of possible formation of antiparticles in laboratory was rigorously confirmed. Consequently, Dirac’s theory of holes [2] was acknowledged [3] as one of the fundamental theories of particle physics. The interest in the formation of cold antiprotons was demonstrated through the development of the “Low Energy Antiproton Ring (LEAR)” at CERN [4], in the early eighties of the preceding century which led in 1995 to the formation of the first antiatom in laboratory, namely the ANTIHYDROGEN by the ATHEN experiment [5] [6] (For more study about Antihydrogen see [7] [8] [9] ). Other interesting experiments were performed after the development of the ATRAP machine at CERN ( [10] [11] [12] ) most of them seeking on the one hand, the production of large number of antihydrogen atoms, protoniums and protonium ions [13]. On the other hand, efforts were made to test the behavior of antimatter from the gravity point of view [14]. Recently, several applications of matter-antimatter interactions at low energy were explored ( [15] [16] [17] ).

Experimental investigations of plasmas involving antiprotons were carried out at the AEgIS, (for a review see [18] [19] ). In this case nonneutral plasmas were studied in detail [20]. Formation of Matter-antimatter plasma, particularly high energetic electron-positron plasms, was also proposed as plausible explanation for the mysterious radiation emitted from pulsers. The motive of the present work is the experimentally confirmed fact that matter-antimatter annihilation occurring by magnetically confined plasma releases much higher energy per unit mass comparative to any other source of propulsion [20] (For a comprehensive account on plasma theory and applications see [21] - [28] ).

Our present work is concerned with two problems. The first is the formation of plasma through the interaction of antiprotons with highly populated hydrogen, deuterium or tritium atoms localized in a material host (e.g., palladium crystal). The second problem is to study plasmas formed by injecting a beam of antiprotons in a magnetic trap filled with protons, deuterons, or tritons.

The main objective of our investigations is to show that huge number of energies could be released from the plasmas created in both cases. Particularly, the employment of this plasma as fuel for engines of space shuttles is strongly emphasized.

The present paper falls in other three sections followed by a complete list of references mentioned in the text. The next section deals with the mathematical formalism of our problems. In Section 3 the results of our calculations are displayed and discussed. In Section 4 the main conclusions drawn from our investigations are presented. The paper concludes with the list of references mentioned in the text.

2. Mathematical Formalism

The present section is split into two parts, the first is devoted to the treatment of antiproton-atom plasmas and the second deals with the antiproton-nucleon plasmas.

2.1. Antiproton-Atom Plasmas

Consider a plasma which is composed of antiprotons (P^-), electrons (e⁻) and hydrogen ions (H⁺), (deuteron ions or triton ions, i.e., D⁺ or T⁺, respectively). For simplicity, it is assumed that the fluid of particles is forming an ideal gas.

Hence, the pressure of each stream is defined by p α = 3 2 k B T α n , where p α ; T α ;

n α are, respectively, pressure, the temperature, and density of the particle gas, whilst k_B is Boltzmann constant. The model is explained by the fluid equations.

The continuity equation for particles is given by.

∂ n α ∂ t + ∂ ∂ x ( n α u α ) = 0. (1)

The equation of motion for particle has the form.

∂ u α ∂ t + u α ∂ u α ∂ x + 1 m α n α ∂ p α ∂ x + q α m α ∂ ϕ ∂ x = 0 , (2)

where α stands for e⁻; P^-; H⁺; D⁺; or T⁺.

Poisson’s equation is expressed by

∂ 2 ϕ ∂ x 2 = e ϵ ( n P − + n e − n i ) (3)

where m_α; u_α; and q_α denote, respectively, the mass, velocity and charge of the particle α in the stream. ϕ is the electrostatic potential and the index i refers to H⁺; D⁺; or T⁺. The fluid system of Equations (1)-(3) are normalized by the dimensionless variables x → ω i x / C s ; t → ω i x / C s ; N = n / n 0 ; U α = u α / C s and ϕ = e ϕ / 2 k B T e ; where ω i is the plasma frequency of particle and n₀ is the

equilibrium density, and C s = 2 k B T e m i is the acoustic wave speed and T_e is the

electron temperature respectively. Thus, the normalized system is described by the equations.

∂ N α ∂ t + ∂ ∂ x ( N α U α ) = 0 (4)

∂ U α ∂ t + U α ∂ U α ∂ x + ν α 0 N α 1 ∂ N α ∂ x + ν α 1 ∂ Φ ∂ x = 0 ∂ 2 Φ ∂ x 2 = N P − + N e − N i (5)

where ν α 0 = 3 4 M i α τ e α , ν α 1 = Q α M i 0 ; M i α = m i / m α , τ e α = T α / T e and

Q α = q α / e with Q e , p − = − 1 ; and Q i = 1 . Stretching coordinates are defined by ξ = ε 1 / 2 ( x − v t ) ; and τ = ε 3 / 2 v t ; and

N α = N α 0 + ε N α 1 + ε 2 N α 2 + ⋯ , U α = U α 0 + ε U α 1 + ε 2 U α 2 + ⋯ and Φ = ε Φ 1 + ε 2 Φ 2 + ⋯ (6)

N i 0 = N P − 0 + N e 0 (7)

∂ Φ 1 ∂ τ + 2 A 1 Φ 1 ∂ Φ 1 ∂ ξ + B 1 ∂ 3 Φ 1 ∂ ξ 3 = 0 (8)

and

h p − 0 + h e 0 − h i 0 = 0. (9)

Hence, the solution of kdv Equation (8) and the first order perturbation of velocity and density of particle are given by

Φ 1 = 3 v 2 A 1 sech 2 ( v / 4 B 1 Y ) U α 1 = 3 ν α 1 v 2 A 1 ( v − U α 0 ) η α 0 sech 2 ( v / 4 B 1 Y ) N α 1 = 3 ν α 1 v 2 A 1 N α 0 η α 0 sech 2 ( v / 4 B 1 Y ) (10)

where, A 1 = h i 2 − h p − 2 − h e 2 h i 1 − h p − 1 − h e 1 , B 1 = 1 h i 1 − h p r − 1 − h e 1 , h α 1 = 2 v η a 0 ( v − U α 0 ) h α 0 , h α 2 = ν α 1 2 η α 0 2 h α 0 ( 3 ( v − U α 0 ) 2 − ν α 0 ) , h α 0 = ν α 1 N α 0 ( v − U α 0 ) 2 − ν α 0 N α 0 2 , η α 0 = 1 ( v − U 00 ) 2 − ν 00 and Y = ξ − v τ . The normalized plasma kinetic energy of zero and first perturbations for velocities and densities is.

2.2. Antiproton-Nucleon Plasma

Consider a system of fluid plasma consisting of two kinds of oppositely charged particles, namely negatively charged antiprotons (P^-) and positively charged Protons (P), Deuterons (D) or Tritons (T). Furthermore, assume that all plasmas are ideal gases. Consequently, the pressure of each plasma is calculated by p α = 3 2 k B T α n , where p α , T α , and n α stand, respectively, for the pressure, temperature, and density of particle stream, whilst k_B is Boltzmann constant. The model is explained by the following system of equations.

The continuity equation

∂ n α ∂ t + ∂ ∂ x ( n α u α ) = 0 , (13)

the equation of motion

∂ u α ∂ t + u α ∂ u α ∂ x + 1 m α n α ∂ p α ∂ x + q α m α ∂ ϕ ∂ x = 0 , (14)

with α = P^-, P, D, T) and

Poisson’s equation

∂ 2 ϕ ∂ x 2 = e ϵ ( n P − − n i ) (15)

where m_α, u_α, and q_α are the mass, velocity, and charge of the particle α. ϕ is the electrostatic potential and the index i refers to P, D, or T.

The Plasma system of Equations (1)-(3) are normalized by the dimensionless variables x → ω i x / C s , t → ω i x / C s , N = n / n 0 , U α = u α / C s and

ϕ = e ϕ / 2 k B T e , where ω i is the plasma frequency of the particlei and n₀ is the equilibrium density, and C s = 2 k B T e m i is the acoustic wave speed and T_e is the electron temperature respectively. Thus, the normalized system is described by.

∂ N α ∂ t + ∂ ∂ x ( N α U α ) = 0 ∂ U α ∂ t + U α ∂ U α ∂ x + ν α 0 N α − 1 ∂ N α ∂ x + ν α 1 ∂ Φ ∂ x = 0 ∂ 2 Φ ∂ x 2 = N P − − N i (16)

where ν α 0 = 3 4 M i α τ e α , ν α 1 = Q α M i 0 ; M i α = m i / m α , τ e α = T α / T e and Q α = q α / e with Q e , p − = − 1 , and Q i = 1 . The stretching coordinates are defined by.

ξ = ε 1 / 2 ( x − v t ) ,

and

τ = ε 3 / 2 v t ,

and

N α = N α 0 + ε N α 1 + ε 2 N α 2 + ⋯ , U α = U α 0 + ε U α 1 + ε 2 U α 2 + ⋯ and Φ = ε Φ 1 + ε 2 Φ 2 + ⋯

N i 0 = N P − 0

∂ Φ 1 ∂ τ + 2 A 1 Φ 1 ∂ Φ 1 ∂ ξ + B 1 ∂ 3 Φ 1 ∂ ξ 3 = 0 (20)

and

h P − 0 − h i 0 = 0 . (21)

Hence, the solution of kdv Equation (8) and the first order perturbation of the velocity and density of the particle are given by

Φ 1 = 3 v 2 A 1 sech 2 ( v / 4 B 1 Y ) U α 1 = 3 ν α 1 v 2 A 1 ( v − U a 0 ) η α 0 sech 2 ( v / 4 B 1 Y ) N α 1 = 3 ν α 1 v 2 A 1 N α 0 η α 0 sech 2 ( v / 4 B 1 Y ) (22)

A 1 = h i 2 − h P − 2 h i 1 − h P − 1 , B 1 = 1 h i 1 − h P − 1 , h α 1 = 2 v η α 0 ( v − U a 0 ) h α 0 .

h α 2 = ν α 1 2 η a 0 2 h a 0 ( 3 ( v − U a 0 ) 2 − ν o 0 ) , h α 0 = ν α 1 N a 0 ( v − U a 0 ) 2 − ν α 0 N α 0 2 ,

η α 0 = 1 ( v − U α 0 ) 2 − ν α 0 and Y = ξ − v τ .

The normalized kinetic energies of the plasma corresponding to the zero and first order perturbations of the velocities and densities are determined by

K E = 1 2 ∑ α ( N α 0 + N α 1 ) ( U α 0 + U α 1 ) 2 K E = 1 2 ∑ m α ( n α 0 + n α 1 ) ( u o 0 + u α 1 ) 2 .

3. Results and Discussion

As shown in tables. The following values of particle masses have been employed in the present work.

Proton mass m_P = 1.6726231 × 10²⁷ kg, Deuteron mass m_D = 3.3476 × 10⁻²⁷ kg; and triton mass m_T = 5.0225 × 10⁻²⁷ kg, M p ¯ = M_P = 1.0073 × 10⁻⁶ g, M_D = 2.0160 × 10⁻⁶ g, M_T = 3.0246 × 10⁻⁶ g, M e − = 5.4858 × 10⁻² g, M P P e − = 3.02245 × 10⁻⁶ g; M D P e − = 5.03985 × 10⁻⁶ g; M T P e − = 7.05705 × 10⁻⁶ g.

Normalized kinetic energy (KE) against Y for electrons, antiprotons and hydrogen, deuterium, or tritium are demonstrated in Figures 1-3 and for antiprotons and protons, deuterons, and tritons are displayed in Figures 4-6.

Data Availability Statements

The data of the calculated energy that support the findings of this study are derivative in at https.//doi.org/10.1016/j.physleta.2014.10.006, reference number [22] [23].

Table 1 Stretching velocity for protons, deuterons, or tritons with electrons and antiprotons for different velocities of antiprotons, ue = Cs, ui = 9 × 107Cs; Ne = 1, NP = 1, and Ni = 2

u_P₁	0.9C_s	0.09C_s	0.009C_s
v_H	1.2243118163856073	1.181085869585997	0.2761483671130386
v_D	1.2243118376157038	1.1834150339194047	0.33189855601735063
v_T	1.2243118459019127	1.1841406167744444	0.35936943838330293

Table 2 Frame stretching velocity for protons, deuterons, or tritons with electrons and for different velocities of antiprotons ui = 9 × 107Cs; NP = 1; and Ni = 1

u_P₁	0.9C_s	0.09C_s	0.009C_s
v_H	1.2243118163856073	1.181085869585997	0.2761483671130386
v_D	1.2243118376157038	1.1834150339194047	0.33189855601735063
v_T	1.2243118459019127	1.1841406167744444	0.35936943838330293

Table 3 The non-normalized kinetic energy KE (J) of Plasma in Palladium for n0 = 6.0221367 × 11023 mol−1; ε = 13.7; ui = 9 × 107Cs; ue = Cs, (H; D; or T, with P; e) Y = 0 at Te (K) = 273

u_P	KE_H	KE_D	KE_T
0.9C_s	3.043700	3.042580	3.042210
0.09C_s	13.571500	13.92480	14.04160
0.009C_s	0.000672722	0.00146563	0.00220382

Table 4 For (P; D; or T, with P) Y = 0

u_P	KE_P	KE_D	KE_T
0.9C_s	1.02086 × 10⁶	1.49712 × 10⁶	1.75179 × 10⁶
0.09C_s	11.4938	1.6.6024	1.9.3171
0.009C_s	0.000249961	0.000333362	0.000376228

Table 5 For Te (K) = 300, 320, 340, 360

T_e (K) = 300	(H; D; or T, with P; e) Y = 0			(P; D; or T, with P) Y = 0
u_P	KE_H	KE_D	KE_T	KE_P	KE_D	KE_T
0.9C_s	3.344720	3.343500	3.343090	1.1218 × 10⁶	1.64518 × 10⁶	1.92504 × 10⁶
0.09C_s	1.491380	1.530200	1.543030	12.6305	18.2444	21.2276
0.009C_s	0.000739254	0.00161059	0.00242178	0.000274682	0.000366332	0.000413437
T_e (K) = 320
0.9C_s	3.567700	3.566400	3.565970	1.1967 × 10⁶	1.7549 × 10⁶	2.05338 × 10⁶
0.09C_s	1590800	1632210	1645900	11.9661	19.461	2.2.6423
0.009C_s	788.538	1717.96	2583.23	0.000292995	0.000390754	0.000441
T_e (K) = 320
0.9C_s	3790690	3789300	3788840	1.2714 × 10⁶	1.86454 × 10⁶	2.18171 × 10⁶
0.09C_s	1690230	1734230	1748770	14.315	20.1	24.05
0.009C_s	837.822	1825.33	2744.68	0.0003113	0.0004157	0.0004686
T_e (K) = 360
0.9C_s	4.013670	4.012200	4.011710	1.3462 × 10⁶	1.9742 × 10⁶	2.31 × 10⁶
0.09C_s	1.789650	1.836240	1.851640	1.5156	2.18	2.547
0.009C_s	887.105	1932.7	2906.14	0.000329	0.000439	0.0004961

4. Conclusion

The one-dimensional hydrodynamic model of linear and nonlinear analyses in the term of electron-positron and proton-antiproton are studied. The linear analyses show a positive root that indicates ion-acoustic waves in the system and the possibility of existing soliton waves. Also, we obtain the calculations of the nonlinear analysis of the deformed kdv equation and the energy released of the formed plasma. The results, see Tables 1-5 and Figures 1-6, show that the lifetime of formed plasma increases inside the palladium crystal 13.7 times concerning the free space. We have the same energy of formed plasma in the case of free space and palladium crystal. The renormalized example study assumes that when the velocity of the antiproton is around the ion-acoustic wave the released energy of the plasma could be reached to 1.15723 × 10⁶ J/mole. On the other hand, if the velocity of the antiproton decreases, the total energy of the plasma system decreases. This result supports very much the conclusion that we suggested model for obtaining huge amount of energy by constructing a plasma state between thermalized hydrogen and antihydrogen in molecular crystals with approximate velocity less than the velocity of ion acoustic wave of the system. This energy could find interesting application in cold fusion and building up engines for space shuttles.

Conflicts of Interest

The authors declare no conflicts of interest regarding the publication of this paper.

Cite this paper

Abdel-Raouf, M.A., Elgendy, A.T. and Youssef, A.A.A.-R. (2022) Plasmas Created in the Interaction of Antiprotons with Atomic and Ionized Hydrogen Isotopes. Suggested Fuels for Space Engines. Journal of High Energy Physics, Gravitation and Cosmology, 8, 14-24. https://doi.org/10.4236/jhepgc.2022.81002

References1

Chamberlain, O., Segrè, E., Wiegand, C. and Ypsilantis, T. (1955) Observation of Antiprotons. Physical Review Journals Archive, 100, 947-950. https://doi.org/10.1103/PhysRev.100.947

Dirac, P.A.M. (1930) A Theory of Electrons and Protons. Proceedings of the Royal Society A, 126, 360-365. https://doi.org/10.1098/rspa.1930.0013

Dirac, P.A.M. (1931) Quantized Singularities in the Electromagnetic Field. Proceedings of the Royal Society A, 133, 60-72. https://doi.org/10.1098/rspa.1931.0130

Hardt, W., Lefevre, P., Mohl, D. and Plass, G. (1980) The CERN Low Energy Antiproton Ring (LEAR) Project. 11th International Conference on High-Energy Accelerators, Geneva, 7-11 July 1980, 819-824. https://doi.org/10.1007/978-3-0348-5540-2_119

Baur, G., et al. (1996) Production of Antihydrogen. Physics Letters B, 368, 251-258. https://doi.org/10.1016/0370-2693(96)00005-6

Blanford, G., et al. (1998) Observation of Atomic Antihydrogen. Physical Review Letters, 80, 3037-3040. https://doi.org/10.1103/PhysRevLett.80.3037

Amoretti, M., et al. (2002) Production and Detection of Cold Antihydrogen Atoms. Nature, 419, 456-459. https://doi.org/10.1038/nature01096

Gabrielse, G., et al. (2002) Background-Free Observation of Cold Antihydrogen with Field-Ionization Analysis of Its States. Physical Review Letters, 89, Article ID: 213401. https://doi.org/10.1103/PhysRevLett.89.213401

Amoretti, M., et al. (2002) Positron Plasma Diagnostics and Temperature Control for Antihydrogen Production. Physical Review Letters, 91, Article ID: 055001. https://doi.org/10.1103/PhysRevLett.91.055001

Gabrielse, G., et al. (2002) Stacking of Cold Antiprotons. Physics Letters B, 548(3-4), 140-145. https://doi.org/10.1016/S0370-2693(02)02850-2

Bonomi, G., et al. (2005) Antihydrogen Production Mechanisms in ATHENA. Nuclear Physics A, 752, 97-100.

Amoretti, M., et al. (2004) Antihydrogen Production Temperature Dependence. Physics Letters B, 583(1-2), 59-67.

Zurlo, N., et al. (2002) Production of Slow Protonium in Vacuum. Hyperfine Interactions, 172, 97-105. https://doi.org/10.1007/s10751-007-9529-0

Rizzini, L. (2012) Further Evidence for Low Energy Protonium Production in Vacuum. The European Physical Journal Plus, 127, Article No. 124. https://doi.org/10.1140/epjp/i2012-12124-9

Alpha Collaboration (2020) Investigation of the Fine Structure of Antihydrogen. Nature, 578, 375-380. https://doi.org/10.1038/s41586-020-2006-5

Abdel-Raouf, M.A. (2009) Matter-Antimatter Physics at Low Energy. Journal of Physics: Conference Series, 1253, Article ID: 12006.

Abdel-Raouf, M.A. (2020) Novel Consequences of Coexistence of Matter and Antimatter in Nature. Journal of High Energy Physics, Gravitation and Cosmology, 6, 251-258. https://doi.org/10.4236/jhepgc.2020.62019

Aghion, S., et al. (2018) Compression of a Mixed Antiproton and Electron Non-Neutral Plasma to High Densities. The European Physical Journal D, 72, Article No. 76.

Higaki, H. (2005) On the Possibility of Non-Neutral Antiproton Plasmas and Antiproton-Positron Plasmas. AIP Conference Proceedings, 793, Article No. 351.https://doi.org/10.1063/1.2121997

Lapointe, M. (1991) Antiproton Powered Propulsion with Magnetically Confined Plasma Engines. Journal of Propulsion and Power, 7, 749-759.https://doi.org/10.2514/3.23388

Elgendy, A.T. (2021) A Global Model of the Collisional Plasma Boundary Sheath Using Step Model. IEEE Transactions on Plasma Science, 49, 711-717.https://doi.org/10.1109/TPS.2021.3049302

El-Tantawy, S.A., Elgendy, A.T. and Ismail, S. (2017) Cylindrical Freak Waves in a Non-Maxwellian Dusty Bulk-Sheath Plasma. An Approximate Solution for the Cylindrical Nonlinear Schrödinger Equation. Physics Letters A, 381, 3465-3471. https://doi.org/10.1016/j.physleta.2017.08.054

El-Tantawy, S.A., Ali Shan, S., Akhtar, N. and Elgendy, A.T. (2018) Impact of Electron Trapping in Degenerate Quantum Plasma on the Ion-Acoustic Breathers and Super Freak Waves. Chaos, Solitons & Fractals, 113, 356-364.https://doi.org/10.1016/j.chaos.2018.04.037

Elgendy, A.T. (2019) Plasma Boundary of Nonlinear Sheath Dynamics for Arbitrary Waveforms in Capacitive Discharge. Journal of Physics: Conference Series, 1253, Article ID: 012010. https://doi.org/10.1088/1742-6596/1253/1/012010

Kouser, S., Shah, K.H., Qureshi, M.N.S. and Shah, H.A. (2020) Nonlinear Ion-Acoustic Waves in e-p-i Plasmas with (r, q) Distributed Electrons and Positrons. AIP Advances, 10, Article ID: 55123. https://doi.org/10.1063/5.0011128

Malik, H.K., Singh, S. and Dahiya, R.P. (1994) Kadomtsev-Petviashvili Solitons in Inhomogeneous Plasmas with Finite Temperature Drifting Ions. Physics Letters A, 195, 369-372. https://doi.org/10.1016/0375-9601(94)90044-2

Singh, S. and Honzawa, T. (1993) Kadomtsev-Petviashivili Equation for an Ion-Acoustic Soliton in a Collisionless Weakly Relativistic Plasma with Finite Ion Temperature. Physics of Fluids B: Plasma Physics, 5, 2093-2097.https://doi.org/10.1063/1.860745

Petrov, A.E., Sarksyan, K.A., Skvortsova, N.N. and Kharchev, N.K. (2001) Observation of Nonlinear Coupling between Drift and Ion-Acoustic Oscillations in Low-Frequency Plasma Turbulence. Plasma Physics Reports, 27, 56-61. https://doi.org/10.1134/1.1339484