Antihydrogen, an atom composed of an antiproton orbited by a positron, was created at CERN in 1995 . Fundamental development of the physics of particle-antiparticle interaction at low energy was achieved by the ATHENA GROUP at CERN in 2006 . It was proved for the first time that the collision of a proton with an antihydrogen leads to the formation of a bound state composed of a proton and an antiproton, (identified by the Authors as antiprotonic hydrogen and known also as protonium), in which both particles are rotating around each other. The binding energy of the antiprotonic hydrogen was calculated as Es = −0.75 keV and the lifetime is τ ~ 1.0 × 10−6 s. Formation of protonium ions was theoretically confirmed by Abdel-Raouf . Cold fusion based on matter-antimatter interaction was proposed and studied by Abdel-Raouf   . On the other hand, the last few years have shown a great deal of interest devoted mainly to the dynamics and applications of plasmas .
An interesting application will be announced in this article in which ion acoustic solitons (IAS) are formed in a compound plasma system of hydrogen and antihydrogen. The study of analog of IAS, using a simplified hydrodynamic model was previously achieved . Since that, there has been a renewed interest in plasmas, covering different plasma modes, instabilities, and other nonlinear effects   . The most important quantities in our study are phase velocity, amplitude, and width, and the existence conditions of IAS. Considerable attention will be paid to the possible role; these quantities may play in the energy carried by the IAS. Let us recall that the interacting solitons may split, fuse, and exchange energies during their collisions. Many investigations have been made of the soliton energy for different plasma situations .
The large-scale oscillations are termed “solitons” . IAS are nonlinear and coherent plasma structures that qualitatively share many properties with the low-frequency ﬂuctuations we have measured in the system of high energetic plasma of hydrogen-antihydrogen collision . Solitons propagate in plasma near the ion acoustic speed, they are coherent and can be periodically excited. In a departure with previous linear perturbation analyses of electron-positron collision, we propose that the generation of IAS may be the source of coherent low frequency oscillations. Thus, IAS are coherent plasma oscillations born out of a balance between the nonlinearity of convection and wave dispersion. These waves have been studied extensively both experimentally and theoretically since the 1960s in double plasma devices   . To interpret these trends, we note that previous studies have shown that similar types of coherent structures propagating concurrently with ion acoustic solitons (IAS) energy. In the present work, we first review the relevant theory of IAS. We then interpret our proposed experimental measurements of hydrogen-antihydrogen collision through within the lines of this analytical theory to obtain the average value of the energy. On the other hand, in particle physics, annihilation is the process that occurs when a particle collides with its respective antiparticle to produce electromagnetic radiation, such as an electron colliding with a positron to produce two photons. High energetic electron-positron collision (more than 209 GeV) may lead to the production of particle-antiparticle pairs of larger masses, which is the basic concept of the LEP collider   at CERN composed of a circular collider the circumference of which is 27 km. The last decade has shown the establishment of the large hadron collider (LHC) based on high energetic proton-proton collisions occurring in a tunnel with 50-kilometer circumference  .
In our study, we will focus on the formation of electron-positron and proton-antiproton plasma prior the annihilation process when two thermalized beams of hydrogen and antihydrogen are injected in a palladium crystal, (see Appendix A). The released energy will be calculated and added to the fusion energy of these two fundamental atoms. Apart from this introduction, the content of this paper is distributed in other three sections. The next section deals with the mathematical model employed for describing the plasma state. Section 3 is devoted to results and discussion of our calculation in addition to the anticipated hydrogen-antihydrogen fusion energy. In Section 4 we present the main conclusions drawn from our investigations. The paper ends with the complete list of references mentioned in the text as well as two Appendices A and B, the first of which deals with the main features of palladium crystals.
2. Mathematical Model
Solitons are known to exist in a dynamic equilibrium with ion-acoustic turbulence, which has been measured in the cathode plume . The waveform for solitons can be derived from the σ fluid equations for ions and electrons and Poisson equations. In the ion frame of reference, the ion continuity for an α-particle is given by:
The equation of motion for α-particle is given by:
where, . For electron, p for positron, pr+ for proton, and pr− for antiproton.
Poisson equation is expressed as follow:
where, , such that and , and , . Introducing the following dimensionless variables , , , and the plasma frequency .
2.1. Linear Analysis
The density , velocity and potential of our novel system are introduced as:
where, and , to obtain the density and velocity of α-particle which described in terms of the electrostatic potential as:
And the neutral condition for the initial densities which is taken in the form:
And the dispersion relation between ω and k which is given by:
where, is the relative mass between a proton mass and electron mass, , , , and such as , and .
2.2. Non-Linear Analysis
In this section, stretching coordinates is introduced as: , . The perturbations of density, velocity for article and electric potential are given by , and . The propose is obtained kdv equation. The continuity and momentum Poisson Equations (1)-(3) are treated by perturbed quantities for , and in terms of and . The first perturbed density and velocity calculations can be expressed as follow:
where, . such that and that Thus, the signs of , and are always similar to . The neurtal condition of four densities and are given by:
To obtain the kdv equation, the coefficient of must be vanished which appears in Poisson equation. Thus, the coefficient of gives the equation where, , to find the phase velocity we solve the expressed equation of
which is transformed into the algebraic equation of degree six where , , , and . Then, the solution of Poisson’s equation can be expressed as:
The solution of kdv is given by for , and are expressed as follow:
where, , and . Thus. The energy of the plasma is given by:
where, , and . The total plasma energy of zero and first order perturbations calculated by
3. Renormalization of Plasma Energy Calculation
Applying to our example injected proton, antiproton, electron, and positron inside the cavity of palladium crystal. The ion Bohm acoustic velocity
, the density of one mole of electron, positron, proton and anti-proton are , The temperature ratio of electron, positron and proton respectively and , in case of free space, the plasma frequency of antiproton = 4.04 × 106 rad/s and in case of palladium crystal where, , . will be expressed in renormalized coordinate as
where, , .
4. Results and Discussion
Our objectives in this section are twofold: first to discuss the solution of the dispersion relation (8) deduced in Section 2, and second to discuss the solution of the kdv Equation (11) derived in the same section. The models of protons-antiprotons-electrons-positrons are studied for unit charge e = 1 and relative mass . The unperturbed densities are set to be , whilest the unperturbed velocities are chosen to be , and , , and .
4.1. Results of Linear Analysis
With respect to the analysis of the dispersion relation (8) for electrons, protons, positrons and anti-protons plasma, the wave number k is rewritten in terms of ω to obtain the algebraic equation of an eight degree in k. In our solution, we select the positive root that indicates the existence of ion-acoustic waves in the system and the possibility of existing soliton waves. The wave number is plotted against frequency to describe the behavior of the dispersion relation. In Figure 1, the dispersion relation has only one positive root at ω = 43.8k, which is independent of the parameters , and Λ. In this regard, the linear dispersion relation is plotted linearly by k against ω. Over ω domain from zero to 30,000 as shown in Figure 1. Hence, the reciprocal of the slope is the phase velocity and is equal to 43.8485.
4.2. Results of Nonlinear Analyses
In the present section we display the calculations of the potential and kinetic energies for nonlinear analyses. In Figure 2, the potential (14) are plotted against the stretched coordinate Y for certain values of the plasma parameters and for the following choices of the unperturbed velocity of antiproton,
Figure 1. Relation between the wave number k for , , , , , and for k against ω.
respectively , 0.009, and 0.0009. The figures show that the potential is attractive and symmetric about Y = 0, i.e., Support the possible existence different bound states between the interacting particles. Also, its range is contracted as the unperturbed velocity of antiproton decelerates at low relative velocities.
In Figure 3, the kinetic energy of the plasma corresponding to the zero and first order perturbations are sketched versus the stretched coordinate Y at the
Figure 2. Schematic diagrams for the potential energy calculated using the following parameters , , , , , at the three different values of .
Figure 3. (a)-(c) schematic diagrams for the kinetic energy of our novel system for , , , , , at different values of the antiproton velocities, .
same parameters used for calculating the potential energy presented in Figure 2. Obviously, the kinetic energy of the plasma increases when the particles are close to each other.
5. Total Energy Released from Plasma (By Renormalized Our System)
The ion Debye radius and time of plasma period in case of free space and in palladium crystal respectively , , , and . KE (J/mole) are calculate at , , and at and . The energy calculates at free space and in palladium crystal where . The energy calculations at our macroscopic study are similar.
In our novel system we could evaluate the plasma energy released of electron, positron, proton and antiproton annihilation in case of free space and in the cavity of palladium crystal. The results show that the time of existing plasma before annihilation around 13.7 time the free space. The amount of energy in case of macroscopic study of the plasma energy is the same in case of free space and palladium crystal. By, changing the velocity of antiproton up to the ion Bohm acoustic wave we could control the amount of energy released from the system. At high speed, the energy of the plasma is very high and within time scale 0.11 s the amount of energy decay very fast. At very low velocity the energy very low but the system is more stable.
6. Total Energy Released from the Proposed Annihilates Mechanism
As we know the electron-positron pair annihilates into two gamma photons the total energy of which is 1.022 MeV. Also, the protonium annihilation yields energy of the order of 2 × 0.938, i.e. 1.876 GeV. Thus, it is anticipated that the yield of the hydrogen-antihydrogen system is 1886.22 MeV (=3.02205748 × 10−13 Joule). Which means that we have 3.02205748 × 10−13 × 6.02214 × 1023 = 18.1992532386 × 1010 Joule/mole = 1.882 × 1016 Joule/kg. Comparing this result with yields of chemical reaction, nuclear fission and nuclear fusion which are of the order of 1 × 107, 8 × 1012, 1 × 1014 Joule/kg, respectively, we not that our suggested experiment is superior.
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 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 ×1020 J/mole. On the other hand, if the velocity of the anti-proton 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.
In this Appendix we present a brief account on the properties of palladium crystals. The palladium could absorb more than 900 times its own volume of hydrogen and antihydrogen. It expands and becomes harder, stronger, and less ductile in the process. Palladium is a metal in group VIII of the periodic table with atomic number 46, an atomic weight of 106.4, and a density of 12 Mg/m3. Its melting point is 1552 C, and it boils at 3980 C. The electronic configuration of is (Kr) (4d10) (5s0). Its atomic radius is 0.137 nm. At room temperature, Palladium has a face-centered cubic crystal structure with a = 0.388 nm. The linear coefficient of thermal expansion is 12.29 m/m/K and the thermal conductivity is 0.718 W/cm/K.
In the Appendix, we present closed forms of the final analyses of the terms A₁, B₁ and C₁ appeared in Section 2.
 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.
 Abdel-Raouf, M.A. (1989) Possibility of Producing an Intensive Beam of Antihydrogens and Its Consequences. 12th International Conference on Few Body Problems in Physics, Vancouver, 2-8 July 1989, A17.
 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.
 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.
 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.
 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.
 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.
 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.
 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.
 Myers, S. and Picasso, E. (1990) The Design, Construction and Commissioning of the CERN Large Electron-Positron Collider. Contemporary Physics, 31, 387-403.
 Georgin, M.P., Jorns, B.A. and Gallimore, A.D. (2018) Experimental Evidence for Ion Acoustic Solitons in the Plume of a Hollow Cathode. 2nd Space Propulsion Conference, Sevilla, 14-18 May 2018, 12 p.