Effects of the two-color laser on relativistic muon motion

Effects of the two-color laser on relativistic muon motion 1 corresp first-author Corresponding author, shi-chunhua@163.com Corresponding author, shi-chunhua@163.com shi-chunhua@163.com 2 1 Department of Biomedical Engineering, Changzhi Medical College, Changzhi046000, China Department of Biomedical Engineering, Changzhi Medical College, Changzhi046000, China Department of Physics, Changzhi Medical College, Changzhi046000, China Department of Physics, Changzhi Medical College, Changzhi046000, China Considering the mixture after muon-catalyzed fusion ( μ CF) reaction as overdense plasma, we analyze muon motion in the plasma induced by a linearly polarized two-colour laser, particularly, the effect of laser parameters on the muon momentum and trajectory. The results show that muon drift along the propagation of laser and oscillation perpendicular to the propagation remain after the end of the laser pulse. Under appropriate parameters, muon can go from the skin layer into field-free matter in a time period of much less than the pulse duration. The electric-field strength ratio or frequency ratio of the fundamental to the harmonic has more influence on muon oscillation. The laser affects little on other particles in the plasma. Hence, in theory, this work can avoid muon sticking to α effectively and reduce muon-loss probability in μ CF. Muon Two-color laser Overdense plasma Muon sticking to alpha I. INTRODUCTION 1 The muon catalyzed fusion ( μ CF) of hydrogen isotopes, especially D-T fusion, has been studied as a realizable candidate of an energy source as thermal energies [1, 2, 3, 4, 5, 6, 7]. In the liquid mixture of D 2 and T 2 , the muon plays the role of a catalyzer through the formation of a muonic molecule. After the fusion process, the muon is released normally and again it is utilized for a subsequent fusion until the muon decays. However, in fact, after the μ CF, some muon may be strapped by α particles and form into ( αμ ) + . This has become a major limitation on the number of fusion cycles catalyzed by each muon [ 8 8 , 9 9 ]. Because of muon-loss, μ CF is far from the aim to attain energy plus. In the past decade, researchers found sticking phenomena and new insights in μ CF [ 10 10 , 11 11 , 12 12 , 13 13 , 14 14 ], and in experiments μ CF numbers per muon extend to 150 during the lifetime of ∼ 2.2×10 -6 s [ 15 15 ]. Althoght this catalyzed efficiency is not enough for atomic energy application, muon regeneration from the α -sticking has been increasingly studied. In this work, we were trying to avoid muon- α sticking by inducting an external force. Here, the laser is a good choice. It is well-known that ultra-short, ultra-intense laser pulses are potentially useful for a variety of applications. Recently, ultra-short lasers have been applied extensively to studies of electron acceleration schemes due to their potential of offering compact and low-cost setups through accelerator gradients several orders of magnitude higher than those of conventional RF linear accelerators [ 16 16 , 17 17 , 18 18 ]. Unlike one-color laser, a two-color laser beam is the sum of two copropagating pulsed radially polarized laser beams, with central angular frequencies ω and nω , of unequal pulse duration and peak power. Also, the relative phase of the two-color laser pulses adds one more degree of freedom [ 19 19 ]. Previous simulations have shown that a one-color pulsed radially polarized laser beam can accelerate an initially stationary electron only up to 40% of the theoretical energy gain limit [ 20 20 , 21 21 ], while a two- color pulsed beam can accelerate an electron by over 90% of the one-color beam’s theoretical gain limit, for a given total energy and pulse duration [ 22 22 ]. So, in order to make the muon well out of the mixture, the choice of the two-color laser becomes especially important [ 23 23 , 24 24 ]. In this paper, the two-color laser is induced in the D-D fusion, and muon can be moved at a certain track so that μ is not captured by α . To study the muon drift of D-D fusion in the two-color laser, it is necessary to choose some appropriate laser parameters. Here, we explore the possibility of the muon drift and trajectory, particularly; study influences of various laser-pulse parameters (i.e., frequency, intensity, and relative phase) on the muon drift and investigate coherent effects on other particles in a two-level field. II. COMPUTATIONAL MODEL AND METHOD 1 When a cycle of μ CF(d-d): μ +d+d → dd μ → α +n+ μ is just finished, the α and μ particles are existing independently, i.e. the μ ions are not captured by the α particles. Then, at the moment, there is the mixture which consists of α , neutron (n), electron (e), d, and μ in the fusion reaction target. It is known that the use of a plasma medium is an attractive way of achieving laser-driven electron acceleration [ 17 17 ]. So, we consider the above mixture as the plasma with a frequency ω p , ( <math id="M1"><mrow><msubsup><mi>ω</mi><mtext>p</mtext><mn>2</mn></msubsup><mo>=</mo><msubsup><mi>ω</mi><mrow><mtext>p</mtext><mi>α</mi></mrow><mn>2</mn></msubsup><mo>+</mo><msubsup><mi>ω</mi><mrow><mtext>pe</mtext></mrow><mn>2</mn></msubsup><mo>+</mo><msubsup><mi>ω</mi><mrow><mtext>p</mtext><mi>μ</mi></mrow><mn>2</mn></msubsup></mrow></math> , ω p α , ω pe and ω p μ are frequency of α , e and μ , respectively). Also, as ω pe > > ω p μ > > ω p α , we judge <math id="M2"><mrow><msubsup><mi>ω</mi><mtext>p</mtext><mn>2</mn></msubsup><mo>=</mo><msubsup><mi>ω</mi><mrow><mtext>pe</mtext></mrow><mn>2</mn></msubsup></mrow></math> . At the same time of the above mixture formation, the two-color laser pulse is turned on and off adiabatically, i.e. <math display="block" id="M3"><mrow><mi>E</mi><mo stretchy="false">(</mo><mi>t</mi><mo stretchy="false">)</mo><mo>=</mo><mi>f</mi><mo stretchy="false">(</mo><mi>t</mi><mo stretchy="false">)</mo><mo>⋅</mo><mo>{</mo><msub><mi>E</mi><mtext>f</mtext></msub><mi>exp</mi><mo stretchy="false">(</mo><mo>−</mo><mi>i</mi><msub><mi>ω</mi><mtext>f</mtext></msub><mi>t</mi><mo stretchy="false">)</mo><mo>+</mo><msub><mi>E</mi><mtext>h</mtext></msub><mi>exp</mi><mo stretchy="false">[</mo><mo>−</mo><mi>i</mi><mo stretchy="false">(</mo><msub><mi>ω</mi><mtext>h</mtext></msub><mi>t</mi><mo>+</mo><mi>φ</mi><mo stretchy="false">)</mo><mo stretchy="false">]</mo><mo>}</mo><mn>,</mn></mrow></math> where, subscript f denotes the fundamental part in the two-color laser pulse and subscript h denotes the harmonic part; E f and E h are amplitude of the electric-field strength of the fundamental laser field and the harmonic field, respectively; for the frequency ω f and ω h , the n th harmonic field means ω h = nω f ; φ is the relative phase between the fundamental field and the harmonic; and f ( t ) is the shape of the laser pulse, i.e. <math display="block" id="M4"><mrow><mi>f</mi><mo stretchy="false">(</mo><mi>t</mi><mo stretchy="false">)</mo><mo>=</mo><mi>sec</mi><msup><mi>h</mi><mn>2</mn></msup><mo stretchy="false">(</mo><mi>t</mi><mo>/</mo><mi>τ</mi><mo stretchy="false">)</mo><mn>,</mn></mrow></math> where, τ is the pulse duration. When the plasma frequency ω p meets the relation of ω p > > ω h > ω ( ω =( ω f + ω h )/2 is the equivalent frequency of a two-color laser), the laser beam can be considered propagating in overdense plasma [ 25 25 ]. For the sake of simplicity, we assume that the linearly polarized laser beam propagates perpendicular to the plane surface of the target. Let us direct the axis X along the propagation of the laser pulse, the axis Y along the electric-field strength, and the axis Z along the magnetic-field strength. Throughout the paper, the system of units m e = e = c = 1 is used. A laser beam propagates in overdense plasma with a plasma frequency ω p greater than the laser frequency ω in the skin layer having a depth δ = 1/ ω p . The boundary conditions for the electric and magnetic fields on the surface vacuum-overdense plasma are of the well-known form: <math display="block" id="M5"><mrow><mtable columnalign="left"><mtr columnalign="left"><mtd columnalign="left"><mrow><msup><mi>F</mi><mrow><mtext>in</mtext></mrow></msup></mrow></mtd><mtd columnalign="left"><mrow><mo>=</mo><mfrac><mn>2</mn><mrow><mn>1</mn><mo>+</mo><msqrt><mi>ε</mi></msqrt></mrow></mfrac><mi>F</mi><mo stretchy="false">(</mo><mtext>t</mtext><mo stretchy="false">)</mo><mo>≈</mo><mfrac><mn>2</mn><mrow><msqrt><mi>ε</mi></msqrt></mrow></mfrac><mi>F</mi><mo stretchy="false">(</mo><mtext>t</mtext><mo stretchy="false">)</mo><mn>,</mn></mrow></mtd></mtr><mtr columnalign="left"><mtd columnalign="left"><mrow><msup><mi>B</mi><mrow><mtext>in</mtext></mrow></msup></mrow></mtd><mtd columnalign="left"><mrow><mo>=</mo><mfrac><mrow><mn>2</mn><msqrt><mi>ε</mi></msqrt></mrow><mrow><mn>1</mn><mo>+</mo><msqrt><mi>ε</mi></msqrt></mrow></mfrac><mi>F</mi><mo stretchy="false">(</mo><mtext>t</mtext><mo stretchy="false">)</mo><mo>≈</mo><mn>2</mn><mi>F</mi><mo stretchy="false">(</mo><mtext>t</mtext><mo stretchy="false">)</mo><mn>,</mn></mrow></mtd></mtr></mtable></mrow></math> where, F in is the electric-field amplitude inside the plasma, F ( t ) meets the relation of F ( t )=| E ( t )| and varepsilon is the dielectric constant of plasma, i.e. <math id="M6"><mrow><mi>ε</mi><mo>=</mo><mn>1</mn><mo>−</mo><msubsup><mi>ω</mi><mtext>p</mtext><mn>2</mn></msubsup><mo>/</mo><msup><mi>ω</mi><mn>2</mn></msup><mo>≈</mo><mo>−</mo><msubsup><mi>ω</mi><mtext>p</mtext><mn>2</mn></msubsup><mo>/</mo><msup><mi>ω</mi><mn>2</mn></msup></mrow></math> . Then, the Newton classical equations for muon momentum P x and P y inside the skin layer ( x>0 ) are of the form: <math display="block" id="M7"><mrow><mrow><mo>{</mo><mrow><mtable columnalign="left"><mtr columnalign="left"><mtd columnalign="left"><mrow><mi>k</mi><mfrac><mrow><mtext>d</mtext><msub><mi>P</mi><mi>x</mi></msub></mrow><mrow><mtext>d</mtext><mi>t</mi></mrow></mfrac><mo>=</mo></mrow></mtd><mtd columnalign="left"><mrow><mn>2</mn><msub><mi>ν</mi><mi>y</mi></msub><mi>f</mi><mo stretchy="false">(</mo><mi>t</mi><mo stretchy="false">)</mo><mtext>exp</mtext><mo stretchy="false">(</mo><mo>−</mo><mi>x</mi><mo>/</mo><mi>δ</mi><mo stretchy="false">)</mo><mo>{</mo><msub><mi>E</mi><mi>f</mi></msub><mi>sin</mi><mo stretchy="false">(</mo><msub><mi>ω</mi><mi>f</mi></msub><mi>t</mi><mo stretchy="false">)</mo></mrow></mtd></mtr><mtr columnalign="left"><mtd columnalign="left"><mrow/></mtd><mtd columnalign="left"><mrow><mo>+</mo><msub><mi>E</mi><mi>h</mi></msub><mi>sin</mi><mo stretchy="false">(</mo><msub><mi>ω</mi><mi>h</mi></msub><mi>t</mi><mo>+</mo><mi>φ</mi><mo stretchy="false">)</mo><mo>}</mo></mrow></mtd></mtr><mtr columnalign="left"><mtd columnalign="left"><mrow><mi>k</mi><mfrac><mrow><mtext>d</mtext><msub><mi>P</mi><mi>y</mi></msub></mrow><mrow><mtext>d</mtext><mi>t</mi></mrow></mfrac><mo>=</mo></mrow></mtd><mtd columnalign="left"><mrow><mn>2</mn><mfrac><mi>ω</mi><mrow><msub><mi>ω</mi><mi>p</mi></msub></mrow></mfrac><mi>f</mi><mo stretchy="false">(</mo><mi>t</mi><mo stretchy="false">)</mo><mtext>exp</mtext><mo stretchy="false">(</mo><mo>−</mo><mi>x</mi><mo>/</mo><mi>δ</mi><mo stretchy="false">)</mo><mo>{</mo><msub><mi>E</mi><mi>f</mi></msub><mi>cos</mi><mo stretchy="false">(</mo><msub><mi>ω</mi><mi>f</mi></msub><mi>t</mi><mo stretchy="false">)</mo></mrow></mtd></mtr><mtr columnalign="left"><mtd columnalign="left"><mrow/></mtd><mtd columnalign="left"><mrow><mo>+</mo><msub><mi>E</mi><mi>h</mi></msub><mi>cos</mi><mo stretchy="false">(</mo><msub><mi>ω</mi><mi>h</mi></msub><mi>t</mi><mo>+</mo><mi>φ</mi><mo stretchy="false">)</mo><mo>}</mo></mrow></mtd></mtr><mtr columnalign="left"><mtd columnalign="left"><mrow/></mtd><mtd columnalign="left"><mrow><mo>−</mo><mn>2</mn><msub><mi>ν</mi><mi>x</mi></msub><mi>f</mi><mo stretchy="false">(</mo><mi>t</mi><mo stretchy="false">)</mo><mtext>exp</mtext><mo stretchy="false">(</mo><mo>−</mo><mi>x</mi><mo>/</mo><mi>δ</mi><mo stretchy="false">)</mo><mo>{</mo><msub><mi>E</mi><mi>f</mi></msub><mi>sin</mi><mo stretchy="false">(</mo><msub><mi>ω</mi><mi>f</mi></msub><mi>t</mi><mo stretchy="false">)</mo></mrow></mtd></mtr><mtr columnalign="left"><mtd columnalign="left"><mrow/></mtd><mtd columnalign="left"><mrow><mo>+</mo><msub><mi>E</mi><mi>h</mi></msub><mi>sin</mi><mo stretchy="false">(</mo><msub><mi>ω</mi><mi>h</mi></msub><mi>t</mi><mo>+</mo><mi>φ</mi><mo stretchy="false">)</mo><mo>}</mo></mrow></mtd></mtr></mtable></mrow></mrow></mrow></math> where k = m μ / m e ≈ 207, δ = 1/ ω p is the depth of the skin layer, and x =0 is the surface of the target. The velocity ν x and ν y are given by <math display="block" id="M8"><mrow><mtable columnalign="left"><mtr columnalign="left"><mtd columnalign="left"><mrow><msub><mi>ν</mi><mi>x</mi></msub></mrow></mtd><mtd columnalign="left"><mrow><mo>=</mo><mfrac><mrow><mtext>d</mtext><mi>x</mi></mrow><mrow><mtext>d</mtext><mi>t</mi></mrow></mfrac><mo>=</mo><mfrac><mrow><msub><mi>P</mi><mi>x</mi></msub></mrow><mrow><msqrt><mrow><msup><mn>1</mn><mn>2</mn></msup><mo>+</mo><msubsup><mi>P</mi><mi>x</mi><mn>2</mn></msubsup><mo>+</mo><msubsup><mi>P</mi><mi>y</mi><mn>2</mn></msubsup></mrow></msqrt></mrow></mfrac><mo>;</mo></mrow></mtd></mtr><mtr columnalign="left"><mtd columnalign="left"><mrow><msub><mi>ν</mi><mi>y</mi></msub></mrow></mtd><mtd columnalign="left"><mrow><mo>=</mo><mfrac><mrow><mtext>d</mtext><mi>y</mi></mrow><mrow><mtext>d</mtext><mi>t</mi></mrow></mfrac><mo>=</mo><mfrac><mrow><msub><mi>P</mi><mi>y</mi></msub></mrow><mrow><msqrt><mrow><msup><mn>1</mn><mn>2</mn></msup><mo>+</mo><msubsup><mi>P</mi><mi>x</mi><mn>2</mn></msubsup><mo>+</mo><msubsup><mi>P</mi><mi>y</mi><mn>2</mn></msubsup></mrow></msqrt></mrow></mfrac><mn>.</mn></mrow></mtd></mtr></mtable></mrow></math> III. NUMERICAL RESULTS AND ANALYSIS 1 We assume that there is no muon motion along the magnetic-field strength, i.e. the axis Z . Besides, collisions of muons are neglected in our consideration. Eqs. (4) and (5) can be solved with the initial conditions that at the time instance t = -∞ a muon is at rest when x = y = 0. Of course, the average muon velocity along the Y axis is zero. After the end of a laser pulse, a muon has zero velocity in the Y direction and nonzero velocity in the X direction. A. The effect of electric-field strength 2 The relativistic plasma frequency is ω p = (4 πn e / γ ) 1/2 , γ =[1+( F in / ω ) 2 ] 1/2 is the relativistic factor. From Eq. (1) to Eq. (5), we fixed the value ω p / ω f = 21/2 and ω f τ / π =40 in our numerical computation which correspond to the pulse duration τ =83 fs and the electron number density n e =8×10 22 cm -3 . For the two-color laser, we take the 2 th harmonic field ω h =2 ω f and the relative phase φ =0. Based on the above values, by solving the Eqs. (4) and (5), the muon motion is determined only by one dimensionless parameter g = 4 E f /(3 ω f ) and the amplitude ratio of the electric-field strength E f / E h . The value g =1 corresponds to the peak laser intensity of 5×10 17 W/cm 2 (at the laser photon energy ħ ω = 1.5 eV). Figure 1 Figure 1 shows the dependence of muon momentum P x (in units of m e c , in this paper, the units of the momentum are all m e c .) on the value of g and the different ratio of E f to E h , where we take E h = E f and 3 E h = E f , respectively. The value g from 0 to 100 corresponds to the peak laser intensity of 0∼5×10 21 W/cm 2 . The results show that the two curves change in a similar way basically. The momentum increases with the laser intensity, but it is a nonrelativistic quantity ( P x < < m e c ) even at the relativistic peak laser intensity 5×10 21 W/cm 2 . Fig. 1 Fig. 1 also shows the different ratio of the fundamental electric-field strength to the harmonic has effects on the muon momentum P x , but the effect is obvious at g =0∼70 and is not obvious at g > 70. Figure 2 Figure 2 demonstrates the muon momentum P x and P y at the leading edge of the laser pulse for the typical value of g =80 (the peak laser intensity of the fundamental wave is 3.2×10 21 W/cm 2 according to Eqs. (4) and (5), with E h = E f in Fig. 2(a1–a3) and 2Eh=Ef Fig. 2(a 1 –a 3 ) and 2 E h = E f in Fig. 2(b1–b3) Fig. 2(b 1 –b 3 ) . As can be seen in the figures, a muon goes from the skin layer inside field-free matter before the laser pulse (∝ sec h 2 ( t / τ )) reaches its maximum because of the small depth of the skin layer. i.e. a muon stops acceleration on the rear border of the skin layer. It is also obvious that P y → 0 when a muon leaves the skin layer, because of adiabaticity of the laser pulse. In addition, the different E f / E h ratio affects little on the muon momentum in the X axis, but obviously in the Y axis ( Figs. 2(a3) and 2(b3) Figs. 2(a 3 ) and 2(b 3 ) ). The muon trajectory in the plane for g =80 is shown in Fig. 3 Fig. 3 . The dimensionless muon coordinates x and y are given in units of the skin depth δ . The amplitude of muon oscillations in the transverse direction Y is much less than its drift motion in the direction of laser pulse propagation (the X axis). This indicates that a muon penetrates many skin depths into field-free matter even before the laser pulse reaches its maximum. The amplitude and shape of muon oscillations in the transverse direction Y in Figs. 3(a) and 3(b) Figs. 3(a) and 3(b) differ from each other, due to different harmonic wave electric-field strength. At a greater m = E f / E h , the amplitude of muon oscillations in the transverse direction Y reduces and the shape is more complicated. To explain that a muon goes from the skin layer inside field-free matter or penetrates many skin depths into field-free matter before the laser pulse reaches its maximum, we make an estimate qualitatively. We consider that the muon drift velocity in the direction of laser pulse propagation has the estimate ν x = N ( E f / ω p ) 2 from Eq. (4). Then, the time t p for penetration of the skin layer having the depth δ is approximately, i.e. <math id="M9"><mrow><msub><mi>t</mi><mtext>p</mtext></msub><mo>=</mo><mi>δ</mi><mo>/</mo><msub><mi>ν</mi><mi>x</mi></msub><mo>=</mo><msub><mi>ω</mi><mtext>p</mtext></msub><mo>/</mo><mo stretchy="false">(</mo><mi>N</mi><msubsup><mi>E</mi><mtext>f</mtext><mn>2</mn></msubsup><mo stretchy="false">)</mo><mo>≈</mo><mn>48</mn><mtext> fs</mtext></mrow></math> . The result shows t p < τ , i.e. the time t p at which a muon goes into field-free matter from the skin layer is much less than the laser duration τ . Thus, a muon goes quickly from the skin layer into field-free matter. Besides, in order to better illustrate effects of the two-color laser intensity on muon motion, we calculate the muon momentum and trajectory in one-color laser ( Fig. 4 Fig. 4 , in a peak laser intensity of 6.4×10 21 W/cm 2 and other parameters being the same value as the above. In Fig. 4(c) Fig. 4(c) the maximum of muon momentum P x in the one-color laser field is less than in Fig. 2(a1) Fig. 2(a 1 ) . It shows that the influence of the two-color laser intensity on the momentum is not a simple addition of two one-color laser beams, and a two-color pulsed beam can accelerate muon to the same maximum with lower intensity. In Figs. 4(a) and 4(b) Figs. 4(a) and 4(b) , the oscillations of muon momentum P y in the one-color laser field are different from Fig. 2(a2) and 2(a3) Fig. 2(a 2 ) and 2(a 3 ) . The muon oscillations have just one peak in every cycle in one-color laser field ( Fig. 4(b) Fig. 4(b) ) but two peaks in two-color laser field ( Fig. 2(a3) Fig. 2(a 3 ) ). The muon trajectory in the one-color laser field is shown in Fig. 4(d) Fig. 4(d) , where the amplitude of muon oscillations in the transverse direction Y is obviously less than that in the two-color laser field ( Fig. 3(a) Fig. 3(a) ). In Fig. 3(a) Fig. 3(a) , the muon oscillation in each cycle has two peaks (red and black arrows) in the two-color laser field, but in Fig. 4(d) Fig. 4(d) it has a single peak (red arrow). Certainly, in the one-color laser field, a muon goes from the skin layer inside field-free matter before the laser pulse reaches its maximum, and P y → 0 when a muon leaves the skin layer. However, a muon takes a longer time to go into field-free matter from the skin layer in the one-color field than in the two-color field. B. The effect of frequency 2 To evaluate the frequency effect, we still fix the pulse duration τ = 83 fs and the electron density n e =8×10 22 cm -3 . For the two-color laser, we take 2 E h = E f , E f =8.0×10 20 W/cm 2 and φ =0. Therefore, the muon motion is determined only by the dimensionless parameter ω p / ω ( ω =( ω f + ω h )/2) and the frequency ratio of the fundamental laser field ω f to the harmonic field ω h . Figure 5 Figure 5 shows the dependence of muon momentum P x on ω p / ω , and on ω f / ω h , at ω h = ω f , ω h = 2 ω f and ω h = 4 ω f , respectively. The ω p / ω changes from 5 to 20 so that the plasma frequency ω p is far larger than the laser frequency ω . In Fig. 5 Fig. 5 , the three curves change in a similar way basically, but different frequency ratio affects the muon momentum P x , especially at ω p / ω = 6 14. Besides, the three curves have different peak positions. Figure 6 Figure 6 shows the muon momentum P x and P y at the leading edge of the laser pulse for typical value of the peak laser intensity of the fundamental wave being 8.0×10 20 W/cm 2 and ω p / ω =8 according to Fig. 5 Fig. 5 . We take ω h =2 ω f in Figs. 6(a1) and 6(a2) Figs. 6(a 1 ) and 6(a 2 ) and ω h = 4 ω f in Figs. 6(b1) and 6(b2) Figs. 6(b 1 ) and 6(b 2 ) . From the figures, a muon goes from the skin layer inside field-free matter before the laser pulse reaches its maximum, and P y → 0 when a muon leaves the skin layer. In addition, the frequency ratio affects the muon momentum obviously in the X axis and the Y axis. At a higher frequency ratio, i.e. a greater n = ω h / ω f , the maximum of the muon momentum P x increase, and so does the amplitude of the muon oscillations of P y in the Y axis, with a more complicated shape. Figure 7 Figure 7 shows the muon trajectory in the plane for the peak laser intensity of fundamental wave 8.0×10 20 W/cm 2 and ω p / ω =8. The dimensionless muon coordinates x and y are given in units of the skin depth δ . The amplitude of muon oscillations in the transverse direction Y is much less than its drift motion in the direction of laser pulse propagation (the X axis). This indicates that a muon penetrates the skin depths into field-free matter even before the laser pulse reaches its maximum. At a greater n = ω h / ω f , with an increased muon momentum in the X and the Y axises, the amplitude of muon oscillations in the transverse direction Y increases while the frequency of muon oscillations reduces. Thus, the muon oscillations in the Y axis reduce and at the same time the muon leaves the skin layer by a faster speed in the X axis. In the same way, based on above parameters, the time t p at which a muon goes into field-free matter from the skin layer is about 44 fs. A smaller-than- τ t p means that a muon penetrates the skin depth into field-free matter before the laser pulse reaches its maximum. So, the muon leaves the skin layer faster to enter the free field, at ω p / ω =8 and peak laser intensity of E f =8.0×10 20 W/cm 2 , than the results in Sec. III. A at much larger parameters of g and E f . C. The effect of relative phase 2 For the two-color laser with the parameters of ω p / ω =8, E f = 8.0×10 20 W/cm 2 , electron number density n e =8.0×10 22 cm -3 , 2 E h = E f and ω h =2 ω f , the muon motion is determined only by the pulse duration τ and the relative phase φ . The relative phase in two-color lasers means the propagation distance in the plasma media. It can be converted into the media length L = φλ h /(2 π ). Figure 8 Figure 8 shows the muon momentum P x at the relative phase of φ = π /6 ( L = λ h /12) and the pulse duration τ =10 100 fs. At τ < 50 fs, the muon momenturm P x is less than those in Sec. III. A and Sec. III. B, hence a longer t p , the time for penetrating the skin layer of the depth δ . Then, probably, once the laser pulse duration is less than 50 fs, the muon cannot penetrate the skin layer at the end of a pulse. Thus, to study the effect of the relative phase on the muon trajectory, we take τ =50 fs. Figure 9 Figure 9 shows the muon momentum P y at the leading edge of the laser pulse, with peak laser intensity of the fundamental wave being 8.0×10 20 W/cm 2 at ω p / ω =8, τ =50 fs, and the relative phase φ =0 ( L =0), π /4 ( L = λ h /8), π /2 ( L = λ h /4) and 3 π /4 ( L =3 λ h /8). We also find that the relative phase φ has obvious effects on the muon momentum in just the Y axis. The muon oscillations shape and the peak change with the relative phase. The red and the black arrows indicate peaks that change constantly with the relative phase. Because the relative phase has less effect on the muon momentum P x , the muon trajectories are similar to the results above. The muon speed in the X axis is not affected and the muon still can penetrate skin depths into field-free matter. However, as shown in Fig. 10 Fig. 10 , in the transverse direction Y , the muon oscillations are in different amplitudes and shapes at different relative phases. In Figs. 9 Figs. 9 and 10 10 , at greater the relative phases, the muon oscillation peaks marked with red arrows increase while the muon oscillation peaks marked with black arrows reduce consistently. D. The effects of muon initial velocity 2 For all above analysis, initial velocity of muon is neglected in our consideration, i.e. Eqs. (4) and (5) can be solved with the initial conditions that ν ( t = -∞) = 0. In fact, after every μ CF reaction, velocity of muon is not zero but its value is not unique. Most muons are free, and their velocity depends on the energy release of nuclear fusion [ 13 13 ]. Some muons stick to alpha particles and form muonic helium ions, and the initial sticking probability is 0.912 [ 26 26 ]. Where muonic helium ions are formed with energy of 3:47 MeV, they are slowed down toward thermal energy by collision with the surrounding D2 and DT molecules [ 27 27 ]. So, muonic helium ions can be stripped as a result of collisions and some muons regain free, probability of which is 25% [ 28 28 ]. Their velocities are much less. Other muonic helium ions may be stripped by the external force [ 29 29 ] and velocity of free muon are related to the force. Based on the energy release of nuclear fusion, we take two value of initial velocity and solve the Eqs. (4) and (5). Figure 11 Figure 11 shows the muon momentum P x and P y at the leading edge of the laser pulse and where values of relevant parameters are the same with Fig. 2(a1-a3) Fig. 2(a 1- a 3 ) . The only different is the value of muon initial velocity. According to Eqs. (4) and (5), we take the initial velocity ν x = ν y = 0.001 in Fig. 11(a1–a3) Fig. 11(a 1 –a 3 ) , and ν x = ν y = 0.0001 in Fig. 11(b1–b3) Fig. 11(b 1 –b 3 ) , instead of ν x = ν y = 0 in Fig. 2(a1-a3) Fig. 2(a 1- a 3 ) , the unit of velocity is c =3×10 8 m/s. The results shows that a muon goes from the skin layer inside field-free matter before the laser pulse reaches its maximum, and P y →0 when a muon leaves the skin layer. Comparing Figs. 11(a3) and 11(b3) Figs. 11(a 3 ) and 11(b 3 ) , the detail changes of P y , the muon initial velocity has obvious effect on the muon momentum in the Y axis, while it has a little effect on the muon momentum in the X axis. In the case of superintense laser pulses, the initial velocity in the X axis plays a little role in muon drift, for muon has relatively high drift velocities induced by the laser in the X axis. IV. DISCUSSION AND CONCLUSION 1 In addition, with the induction of the two-color laser, we have to consider the influence of the laser on other particles in the overdense plasma. According to similar methods, we study the other particles trajectory in different muzzle velocities. Because the muzzle velocity of other particles is small, the laser has little influence on them, i.e. the motion of those particles produced by the superintense laser pulse is far less than the motion of muons. Also, the duration of laser pulse is so short that the work does not introduce the disadvantage into μ CF reaction in the superintense laser, but make the muon penetrate skin depths into field-free matter to avoid the α -sticking. In this paper, we consider the mixture after μ CF (d-d) reaction as overdense plasma with the frequency ω p , and the two-color lasers with the frequency ω h and ω f are induced. We study the effect of two-color laser parameters on the muon momentum and trajectory, including the peak intensity, ratio of the fundamental electric-field strength to the harmonic, frequency ratio, relative phase and pulse duration. We find that the time at which a muon penetrates into field-free matter from the skin layer is much less than the pulse duration τ under the right conditions. This penetration occurs at the leading edge of the laser pulse. The muon drift in overdense plasma along the propagation of laser radiation produced by a magnetic part of a laser field remains after the end of the laser pulse. Hence, after μ CF (d-d) reaction is finished and before the muon is captured by α , inducting a suitable superintense laser pulse can avoid the α -sticking efficiently. References Gershtein S S , Petrov I V , Ponomarev L I , et al . Zh Eksp Teor Fiz+ , 1980 , 78 : 2099 - 2111 . Zh Eksp Teor Fiz+ Alexander S A , Froelich P , Monkhorst H J . Phys Rev A , 1990 , 41 : 2854 - 2857 . Phys Rev A Ponomarev L I . Contemp Phys , 1990 , 31 : 219 - 245 . Contemp Phys Petitjean C , Atchison F , Heidenreich G , et al . Fusion Technol , 1994 , 25 : 437 - 449 . Fusion Technol Nagamine K , Matsuzaki T , Ishida K , et al . Hyperfine Interact , 1999 , 119 : 273 - 280 . Hyperfine Interact Nakamura S N , Nagamine K , Matsuzaki T , et al . Phys Lett B , 2000 , 473 : 226 - 232 . Phys Lett B Nagamine K . Nucl Energy , 2005 , 3B : 555 - 602 . Nucl Energy Ceperley D and Alder B J . Phys Rev A , 1985 , 31 : 1999 - 2004 . Phys Rev A Davies J D , England J B A , Pyle G J , et al . J Phys G Nucl Partic , 1990 , 10 : 1529 - 1537 . J Phys G Nucl Partic Ishida K , Okamoto K , Kawasaki Y , et al . Phys Rev Lett , 2002 , 89 : 107202 . Phys Rev Lett Nagamine K and Ponomarev L . Nucl Phys A , 2003 , 721 : c863 - c866 . Nucl Phys A Ishida K , Nagamine K , Matsuzaki T . Nucl Phys B , 2005 , 149 : 348 - 350 . Nucl Phys B Kimura S and Bonasera A . Radiat Eff Defect S , 2008 , 163 : 287 - 292 . Radiat Eff Defect S Balin D V , Ganzha V A , Kozlov S M , et al . Phys Part Nuclei , 2011 , 42 : 185 - 214 . Phys Part Nuclei He J T . Physics , 2003 , 18 : 461 - 466 . (in Chinese) Physics Malka V , Faure J , Gauduel Y A , et al . Nat Phys , 2008 , 4 : 447 - 453 . Nat Phys Esarey E , Schroeder C B , Leemans W P . Rev Mod Phys , 2009 , 81 : 1229 - 1285 . Rev Mod Phys Salamin Y I , Hu S X , Hatsagortsyan K Z , et al . Phys Rep , 2006 , 427 : 41 - 155 . Phys Rep Shi C H , Qiu X J , Li R X . J Mod Optic , 2007 , 54 : 1435 - 1444 . J Mod Optic Wong L J and Kärtner F X , Opt Express , 2010 , 18 : 25035 . Opt Express Fortin P L , Piché M , Varin C . J Phys B , 2010 , 43 : 025401 . J Phys B Liang J W and Franz X K . Opt Lett , 2011 , 36 : 957 - 959 . Opt Lett Zhao K and Chu T S . Chem Phys Lett , 2011 , 511 : 166 - 171 . Chem Phys Lett Johnson L A , Palastro J P , Antonsen T M , et al . Phys Rev A , 2013 , 88 : 063804 . Phys Rev A Rastunkov V S and Krainov V P . Phys Rev E , 2004 , 69 : 037402 . Phys Rev E Hu C Y , Hale G M , Cohen J S . Phys Rev A , 1994 , 49 : 4481 - 4488 . Phys Rev A Jones S E . Nature , 1986 , 321 : 127 - 133 . Nature Zeldovich Y B and Gershtein S S . Sov Phys Uspekhi , 1961 , 3 : 593 - 623 . Sov Phys Uspekhi Shi C H , Qiu X J , An W K , et al . Acta Phys Sin , 2005 , 54 : 4087 - 4091 . (in Chinese) Acta Phys Sin 10.13538/j.1001-8042/nst.25.020201 Supported by Natural Science Foundation of Scientific and Technological Innovation Programs of Higher Education Institutions in Shanxi (No.20091025) 2013-10-24 2014-12-17 2014-04-20 2014-04-20 2021-04-11 2014 Nuclear Science and Techniques 2 25