Stopping Power Modulation by Pump Waves of Charged Particles Moving above Two-Dimensional Electron Gases

Yahong Yang; Ya Zhang; Lin Yi; Wei Jiang

doi:10.1155/2021/5580444

Stopping Power Modulation by Pump Waves of Charged Particles Moving above Two-Dimensional Electron Gases

Published online by Cambridge University Press: 01 January 2024

Lin Yi and

Yahong Yang: Affiliation:
School of Physics, Huazhong University of Science and Technology, Wuhan 430074, China
Ya Zhang*: Affiliation:
Department of Physics, Wuhan University of Technology, Wuhan 430070, China
Lin Yi: Affiliation:
School of Physics, Huazhong University of Science and Technology, Wuhan 430074, China
Wei Jiang*: Affiliation:
School of Physics, Huazhong University of Science and Technology, Wuhan 430074, China
*: Correspondence should be addressed to Ya Zhang; [email protected]
Wei Jiang; [email protected]

Article contents

Abstract
Introduction
Quantum Hydrodynamic Model
Result and Discussion
Summary
Data Availability
Conflicts of Interest
References

Rights & Permissions

Abstract

The perturbation electron density and stopping power caused by the movement of charged particles above two-dimensional quantum electron gases (2DQEG) have been studied in numerous works using the quantum hydrodynamic (QHD) theory. In this paper, the QHD is modified by introducing the two-dimensional electron exchange-correlation potential at high density Vxc2DH and the pump wave modulations. Based on the modified QHD, the perturbation electron density and stopping power are calculated for pump waves with various parameters. The results show that the stopping power values are more accurate after considering Vxc2DH. Under the modulation of pump waves with the wavelength from 0.1nm to 0.1cm, the perturbation electron density of 2DQEG and the stopping power of charged particles show periodic changes. Under the modulation of pump waves with λ = 1.76 × 10−4 cm and Φ0=2×1010e/λf, the average stopping power with respect to the time phase θ becomes negative, which means that the charged particles will gain energy and can be accelerated. This is a new phenomenon in the fields of 2DQEG and of great significance in surface physics and surface modification in nanoelectronic devices with beam matter interactions.

Type: Research Article
Information: Laser and Particle Beams , Volume 2021 , 2021 , e14

DOI: https://doi.org/10.1155/2021/5580444 [Opens in a new window]
Creative Commons: This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
Copyright: Copyright © 2021 Yahong Yang et al.

1. Introduction

Accompanied by the two-dimensional quantum electron gases (2DQEG) is confirmed to exist on the surface of the metals and semiconductor heterostructures. The interaction between charged particles and plasma targets has attracted a lot of attention in surface physics [Reference Abrahams, Kravchenko and Sarachik1–Reference Wang and Ma6]. The incident particles can also be used as probes to detect the properties of 2DQEG, or as a powerful tool of surface character modification.

In the study of the interaction between charged particles and plasma targets, people are particularly interested in the analytical calculation of energy loss [Reference Mišković, Segui, Gervasoni and Arista7] and the density distribution of the plasma targets [Reference Mišković, Liu, Goodman and Wang8]. Stopping power is used to describe the energy loss per unit length [Reference Wang and Ma6, Reference Maynard and Deutsch9–Reference Custidiano and Jakas15]. As an important physical quantity to study the interaction between charged particles and the plasma targets, a lot of experiments and theoretical studies have been done in the past decade [Reference Zaremba, Nagy and Echenique16–Reference Frank, Blažević and Grande28]. The stopping power of charged particles can be modulated by the incident velocity of charged particles [Reference Li, Song and Wang5], the amount of charge of charged particles, the density of plasma targets [Reference Li, Song and Wang5], varying altitude parameter between charged particles and 2DQEG [Reference Li, Song and Wang4, Reference Li, Song and Wang5, Reference Raphael and Gennes29], the external magnetic field [Reference Lafleur and Baalrud30, Reference Lafleur and Baalrud31], the spin effects [Reference Zhang, Jiang and Yi32], and so on.

In order to understand the interaction between charged particles and plasma targets in the presence of a laser field, a lot of theoretical studies have been carried out [Reference Arista, Galvão and Miranda33–Reference Wu, Luan and Wang39]. Within the dielectric formalism, Arist et al. [Reference Arista, Galvão and Miranda33] studied the influence of a strong laser field on the stopping power for charged particles in low-density nondegenerate plasmas, and the expression of stopping power is given when a laser field exists. The results show that the laser field reduces the stopping power of plasma and makes it anisotropic through the polarization of the electric field wave. Within the random-phase approximation theory, Nersisyan [Reference Nersisyan and Akopyan34] also investigated the stopping and acceleration effect of protons in a plasma in the presence of an intense radiation field and calculated the stopping power of charged particles in the high-frequency limit. With two-dimensional electrostatic particle-in-cell (2DPIC) simulation, Hu [Reference Hu, Song, Mišković and Wang35] studied energy loss of ion beam when plasma irradiated by a strong laser pulse and the dynamic polarization of plasma. Pertaining results make it clear that stopping power gets strongly modulated by the laser field.

However, these research methods are no longer applicable when studying the stopping power modulation by pump waves of charged particles moving above the 2DQEG. In the study of Nersisyan et al., when considering the interaction between laser and three-dimensional plasma, the expression of the vector potential A of electromagnetic waves was used [Reference Nersisyan and Akopyan34]. However, the vector potential expression of the electromagnetic wave A cannot be written for the interaction of 2DQEG and laser. Another feasible method to implement electromagnetic waves is to ignore the effect of the magnetic field and write down the expression for the electric field term of electromagnetic waves [Reference Mir, Shahid, Jamil, Rasheed and Shahbaz40]. The 2DPIC method cannot be used either to study the interaction between charged particles and 2DQEG. Because the charged particles are above the 2DQEG, a three-dimensional electrostatic particle-in-cell (3DPIC) method is required. However, the 3DPIC method is currently difficult to implement due to a large amount of calculation [Reference Hu, Song, Mišković and Wang35].

In this paper, a modified linear QHD is used to study the modulation effect of pump waves on the perturbation electron density and stopping power caused by the movement of charged particles above the 2DQEG. The outline of this article is as follows. In Section 2, the modified QHD is introduced to include the laser pump wave modulation and the two-dimensional electron exchange-correlation potential of high-density V _xc2DH effects. The general expressions of perturbation electron density of 2DQEG and stopping power of incident charged particles are derived by the modified QHD theory. In Section 3, the numerical results of perturbation electron density of 2DQEG and stopping power of incident charged particles modulated by pump waves with various physical parameters (wavelength λ, amplitude Φ₀, and time phase θ) are given. Section 4 gives a summary. Gauss system of units is adopted in this paper.

2. Quantum Hydrodynamic Model

2.1. Physic Model

The schematic diagram is given in Figure 1. 2DQEG is an idealized model, which corresponds to the electronic components on the surface or interface of metal or semiconductor in a very thin layer (only a few layers of atoms’ thickness). For the convenience of studying, it is generally considered that the thickness of such two-dimensional electron gas is taken negligible [Reference Stern and Howard41]. Consider an infinite 2DQEG composed of free electrons and motionless ions located on the plane z = 0, and use a Cartesian coordinate system to represent its position in the region z > 0 and z < 0. n ₀ is the equilibrium density of electrons and ions. At equilibrium, 2DQEG density satisfies the relation n_i = n_e = n₀ [Reference Zhang, Jiang and Yi32]. Considering that the quantum effect is more significant in the high-density electron gas, we take the equilibrium density of electrons and ions as n ₀ = 5.69 × 10¹⁵cm⁻², which is the density of aluminum surface [Reference Li, Song and Wang5].

Figure 1: The schematic diagram of the interaction between pump waves, charged particles, and 2DQEG.

Charged particles move in the x direction parallel to the 2DQEG plane, with a velocity of $v = 2 v_{B}$ and a height of z ₀ = 9λ_f away from the 2DQEG. Here, v_B = e ²/ℏ is Bohr velocity, $λ_{f} = 1 / \sqrt{2 π n_{o}} = 5.29 \times 10^{- 9} cm$ is Fermi wavelength, and e is the elementary charge. The density of charged particles with charge Z ₁e can be expressed as $ρ_{e x t} = Z_{1} e δ (r - v t) δ (z - z_{0})$ , where r = r(x, y).

We assume that the pump waves are collimated, monochromatic, continuous electromagnetic waves with wavelength λ propagating in the +x direction. Along the +x direction, the pump wave with the wavelength from 0.1nm (γ rays) to 0.1cm (microwave) irradiates the 2DQEG. The electric field component of the pump wave acts on the 2DQEG as an external high-frequency electric field in the plane of x0y, and the reaction of the 2DQEG on the pump wave is without the scope of our consideration [Reference Hu, Song, Mišković and Wang35, Reference Zhang and Jiang42]. The electrons in the 2DQEG are affected by the pump wave electric field, which changes the spatial distribution of electrons in the 2DQEG. Then, the exciting electric field by electrons in 2DQEG acts on the incident particle, which affects the velocity of the incident particle, thus affecting the stopping power of incident particles. In terms of electric field, the pump waves can be expressed as $E = - i k_{0} Φ_{0} \exp (- i (ω_{0} t - k_{0} x))$ , where k₀ is the wave vector of pump waves [Reference Mir, Shahid, Jamil, Rasheed and Shahbaz40]. The time variable of the pump wave is t. The wavelength of the pump wave is λ = 1/k ₀. ω ₀ is the frequency of the pump wave, where $ω_{0} = k_{0} c$ , in which c = 2.998 × 10¹⁰ is the speed of light. The intensity of the pump wave can be calculated by $W_{p} = c {(k_{0} Φ_{0})}^{2} / 8 π$ [Reference Nersisyan and Akopyan34], and the maximum intensity of the pump wave is 10²³W cm⁻² in this paper. When k ₀ is constant, the strength of the pump wave is varying with the pump amplitude Φ₀.

2.2. Mathematical Model

As a powerful research method, QHD theory is widely used in plasma research [Reference Zhang, Song, Zhao and Wang43–Reference Kumar and Ahmad45]. Li et al. [Reference Li, Song and Wang4, Reference Li, Song and Wang5, Reference Zhang, Zhao, Qi, Duan, Xiao and Yang46] conducted a series of studies on the stopping power of charged particles above the 2DQEG using linearized QHD theory, which is consistent with dielectric response theory and experimental results. In our research, the equilibrium electron density of 2DQEG will be disturbed by pump waves and injecting charged particles, generating charged fluid velocity field ue (r, t) and electron gas density n_e (r, t). According to the linearized QHD theory, the density and velocity of electrons on the 2DQEG surface can be described by the continuity equation

(1)

\begin{matrix} \frac{\partial n_{e}}{\partial t} + \nabla_{∥} \cdot (n_{e} u_{e}) = 0, \end{matrix}

the momentum-balance equation

(2)

\begin{matrix} m_{e} (\frac{\partial u_{e}}{\partial t} + u_{e} \cdot \nabla_{∥} u_{e}) = e \nabla_{∥} {Φ |}_{z = 0} - \nabla_{∥} w_{e} + \frac{ℏ^{2}}{2 m_{e}} \nabla_{∥} (\frac{1}{\sqrt{n_{e}}} \nabla_{∥}^{2} \sqrt{n_{e}}) - γ m_{e} u_{e} - e E - V_{x c 2 DH} \frac{\nabla n_{e 1}}{n_{e}}, \end{matrix}

and Poisson’s equation

(3)

\begin{matrix} \nabla^{2} Φ = 4 π e (n_{e} δ (z) - n_{0} δ (z) - Z_{1} δ (r - v t) δ (z - z_{0})) . \end{matrix}

Here, $\nabla_{∥} = \partial / \partial x e_{x} + \partial / \partial y e_{y}$ , whereas in equation (3), the differentiation $\nabla = \partial / \partial x e_{x} + \partial / \partial y e_{y} + \partial / \partial z e_{z}$ is unrestricted. Φ denotes total electrostatic potential built on wake potential Φ_ind caused by the collective excitation of electron gas and the charged particle’s external potential. On the right-hand side of equation (2), the first term is the force of electrons on the 2DQEG surface due to the tangential electric field, and the second term is the force caused by the quantum statistical effect. Here, $w_{e} = π ℏ^{2} n_{e} / m_{e}$ refers to the quantum statistical pressure term acting on electrons in 2DQEG, the third term is the force due to the quantum diffraction effects, the fourth term is the frictional force caused by the interaction between electrons and positive charge background where γ is the friction coefficient, the fifth term is the force on electrons owing to the electric field of pump wave, and the last term is the force on electrons due to electron exchange-correlation potential V _xc2DH. The validity of QHD equations (1) and (2) implies that the off-diagonal components of the pressure tensor and viscosity remain negligible.

2.3. Exchange-Correlation Potential

Exchange-correlation potential is an important quantum effect term. Mir et al. considered the exchange-correlation potential when studying degenerate plasmas [Reference Mir, Shahid, Jamil, Rasheed and Shahbaz40]. Haas derived a quantum fluid equation containing exchange correlation, which is described by an effective potential [Reference Haas47]. Pollack nd Perdew have made an in-depth study on the exchange-correlation potential of two-dimensional electron gas [Reference Pollack and Perdew48]. However, the exchange-correlation potential is not considered in the study of the interaction between incident particles and 2DQEG using the quantum hydrodynamic method so far. The three-dimensional electron exchange-correlation potential is expressed as follows:

(4)

\begin{matrix} V_{x c 3 D^{'}} = \frac{0.985}{3} \frac{e^{2}}{ε} {(n_{0})}^{1 / 3} + \frac{0.985 \times 0.034}{3 a_{B}} \frac{e^{2}}{ε} \frac{18.37 a_{B} {(n_{0})}^{1 / 3}}{1 + 18.37 a_{B} {(n_{0})}^{1 / 3}}, \end{matrix}

where $a_{B} = ε ℏ^{2} / m_{e} e^{2}$ is effective Bohr atomic radius, in which ε = 1 is the relative dielectric constant and m_e is the electron mass [Reference Mir, Shahid, Jamil, Rasheed and Shahbaz40]. We adapt the dimension of the three-dimensional electron exchange-correlation potential to the dimension of the two-dimensional electron exchange-correlation potential based on [Reference Mir, Shahid, Jamil, Rasheed and Shahbaz40].

(5)

\begin{matrix} V_{x c 3 D} = \frac{0.985}{3} \frac{1}{r_{s} {(2 π)}^{1 / 2}} + \frac{0.985 \times 0.034}{3} \frac{18.37 / r_{s} {(2 π)}^{1 / 2}}{1 + 18.37 / r_{s} {(2 π)}^{1 / 2}} . \end{matrix}

The two-dimensional electron exchange-correlation potential of low density is then expressed as follows [Reference Mir, Shahid, Jamil, Rasheed and Shahbaz40]:

(6)

\begin{matrix} V_{x c 2 DL} = 0.5058 (\frac{1.3311}{r_{s}^{2}} (\sqrt{1 + 1.5026 r_{s}} - 1) - \frac{1}{r_{s}}), \end{matrix}

and the two-dimensional electron exchange-correlation potential of high density is expressed as follows [Reference Mir, Shahid, Jamil, Rasheed and Shahbaz40]:

(7)

\begin{matrix} V_{x c 2 DH} = \frac{- 0.6002}{r_{s}}, \end{matrix}

where $r_{s} = {(2 π n_{0} a_{B}^{2})}^{- 1 / 2}$ is the Wigner–Seitz radius associated with the 2DQEG density [Reference Pollack and Perdew48].

The V _xc3D (solid line), the V _xc2DL (dashed line), and the V _xc2DH (dotted line) are shown in Figure 2. Here, the value of V _{xc3 D} is positive, and the values of V _xc2DL and V _xc2DH are negative. It is reasonable to assume that two-dimensional electron exchange-correlation potential has completely different properties from its three-dimensional electron homologue. The Wigner–Seitz radius of 2DQEG that we are considering is r_s = 1, the corresponding density is n ₀ = 5.69 × 10¹⁵cm⁻², and the corresponding exchange-correlation potentials are V _{xc3 D} = 0.1362, V _xc2DL = −0.114, and V _xc2DH = −0.6002 (these are dimensionless values). In previous studies of the interaction between charged particles and 2DQEG, the exchange-correlation potential was ignored due to its small value. Here, the absolute value of V _xc2DH is the largest compared with V _xc3 D and V _xc2DL and the absolute value of V _xc2DH is almost six times as much as the absolute value of V _xc3D and V _xc2DL. It is necessary to consider V _xc2DH in the research of interaction between incident charged particles and 2DQEG by using the QHD theory due to the high density in the quantum electron gas. This indicates that increasing density could enhance the exchange-correlation effect, as a function of the density. The exchange-correlation effect has been considered in quantum wells, GaAs/GaAlAs heterostructures, and silicon inversion layers, with the local density approximation method [Reference Stern and Das Sarma49, Reference Brey, Dempsey, Johnson and Halperin50].

Figure 2: The change of electron exchange-correlation potential with Wigner–Seitz radius of 2DQEG.

2.4. Linearized QHD Model

Let $n_{e} (r, t) = n_{0} + n_{e 1} (r, t), u_{e} (r, t) = u_{e 1} (r, t), and Φ (R, t) = Φ_{1} (R, t)$ , where $n_{e 1} ≪ n_{0}$ . $u_{e 1} (r, t)$ , $n_{e 1} (r, t)$ , and $Φ_{1} (R, t)$ are the first-order perturbed quantities of velocity, density, and potential, respectively. We linearize the above equation to obtain the perturbation electron density n _e1 and the induced potential Φ_ind [Reference Li, Song and Wang5]. Here are the linearized continuity equation

(8)

\begin{matrix} \frac{\partial n_{e 1}}{\partial t} + n_{0} \nabla_{∥} \cdot u_{e 1} = 0, \end{matrix}

the momentum-balance equation

(9)

\begin{matrix} m_{e} \frac{\partial u_{e 1}}{\partial t} = e \nabla_{∥} {Φ_{1} |}_{z = 0} \frac{π ℏ^{2}}{m_{e}} \nabla_{∥} n_{e 1} + \frac{ℏ^{2}}{4 m_{e} n_{0}} \nabla_{∥} (\nabla_{∥}^{2} n_{e 1}) - γ m_{e} u_{e 1} - e E - V_{x c 2 DH} \frac{\nabla n_{e 1}}{n_{0}}, \end{matrix}

and Poisson’s equation

(10)

\begin{matrix} \nabla^{2} Φ_{1} = 4 π e (n_{e 1} δ (z) - z_{1} δ (r - v t) δ (z - z_{0})), \end{matrix}

where Φ₁ = Φ_ext + Φ_ind. Φ_ext is the external potential of the incident particle. Φ_ind is the induced potential by the perturbation electron density. Take the Fourier transform of the above linearized equations in space-time, and obtain the perturbation electron density n _e1 and the induced potential Φ_ind:

(11)

\begin{matrix} n_{e 1} (r, t) = n_{0} \frac{z_{1} e^{2}}{2 π m_{e}} \int_{- \infty}^{+ \infty} - \frac{k}{D (k, ω)} e^{- k z_{0}} e^{i k \cdot (r - v t)} d k \\ - \frac{e k_{0}^{2} Φ_{0} n_{0}}{m_{e} {(D_{0} (k_{0}, ω_{0}))}^{2}} (Re (D_{0} (k_{0}, ω_{0})) \cos (k_{0} x - ω_{0} t) + Im (D_{0} (k_{0}, ω_{0})) \sin (k_{0} x - ω_{0} t)), \end{matrix}

(12)

\begin{matrix} Φ_{ind} (R, t) = n_{0} \frac{z_{1} e^{3}}{m_{e}} \int_{- \infty}^{+ \infty} \frac{1}{D (k, ω)} e^{- k z_{0} - k z} e^{i k \cdot (r - v t)} d k \\ + \frac{2 π e^{2} n_{0} k_{0} Φ_{0}}{m_{e} {(D_{0} (k_{0}, ω_{0}))}^{2}} (Re (D_{0} (k_{0}, ω_{0})) \cos (k_{0} x - ω_{0} t) + Im (D_{0} (k_{0}, ω_{0})) \sin (k_{0} x - ω_{0} t)), \end{matrix}

where

(13)

\begin{matrix} D (k, ω) = ω \tilde{ω} - ω_{p}^{2} a_{B} k - \frac{V_{x c 2 DH}}{m_{e}} k^{2} - {(k v_{F})}^{2} \frac{(1 + k^{2} / (2 k_{F}^{2}))}{2}, \\ D_{0} (k_{0}, ω_{0}) = ω_{0} {\tilde{ω}}_{0} - ω_{p}^{2} a_{B} k_{0} - \frac{V_{x c 2 DH}}{m_{e}} k_{0}^{2} - {(k_{0} v_{F})}^{2} \frac{(1 + k_{0}^{2} / (2 k_{F}^{2}))}{2}, \end{matrix}

with $\tilde{ω} = ω + i γ$ . Bohr radius $a_{B} = ℏ^{2} / m_{e} e^{2}$ , Fermi velocity $v_{F} = ℏ k_{F} / m_{e}$ , Fermi wave-number $k_{F} = \sqrt{2 π n_{0}}$ , and electron plasma frequency $ω_{p} = {(2 π n_{0} e^{2} / m_{e} a_{B})}^{1 / 2}$ . k = {k_x, k_y} is a two-dimensional wave vector. Consider that $ω = k \cdot v$ and the projectile velocity v and the x axis are going in the same direction, hence ω = k_xv.

For convenience, we introduce the dimensionless variables: $ω / ω_{p} ⟶ ω$ , $ω_{0} / ω_{p} ⟶ ω_{0}$ , $k_{y} / k_{F} ⟶ k_{y}$ , $k / k_{F} ⟶ k$ , $k_{0} / k_{F} ⟶ k_{0}$ , $v / v_{B} ⟶ v$ , $c / v_{B} ⟶ c$ , $(x - v t) / λ_{f} ⟶ x$ , $y / λ_{f} ⟶ y$ , $z / λ_{f} ⟶ z$ , $z_{0} / λ_{f} ⟶ z_{0}$ , $γ / ω_{p} ⟶ γ$ , and $Φ_{0} / ({e / λ}_{f}) ⟶ Φ_{0}$ . Thus, equation (11) can be reduced to

(14)

\begin{matrix} \frac{n_{e 1} (r, t)}{n_{0}} = A_{n} \int_{- \infty}^{+ \infty} d k_{y} \int_{- \infty}^{+ \infty} d ω B_{n} \cos (y k_{y}) \times (Re (\tilde{D}) \cos (x ω / v) + Im (\tilde{D}) \sin (x ω / v)) \\ + C_{n} (Re (\tilde{D_{0}} (k_{0}, ω_{0})) \cos (x k_{0} + θ) + Im (\tilde{D_{0}} (k_{0}, ω_{0})) \sin (x k_{0} + θ)), \end{matrix}

where $A_{n} = (- Z_{1}) / (2 π r_{s} v)$ , $B_{n} = (k e^{- k z_{0}}) / ({(\tilde{D})}^{2})$ , $C_{n} = (Φ_{0} k_{0}^{2}) / ({(\tilde{D_{0}})}^{2} r_{s})$ ,

(15)

\begin{matrix} \tilde{D} = ω^{2} - \frac{k}{r_{s}} - (\frac{k^{2}}{(2 r_{s}^{2})}) (\frac{1 + k^{2}}{2}) - V_{x c 2 DH} k^{2}, \\ \tilde{D_{0}} = ω_{0}^{2} - \frac{k_{0}}{r_{s}} - (k_{0}^{2} / 24 r_{s}^{2}) (\frac{1 + k_{0}^{2}}{2}) - V_{x c 2 DH} k_{0}^{2}, \end{matrix}

Here, t denotes time variable, ω ₀ is constant when the wavelength λ is fixed, and θ = −ω ₀t is used to modulate the phase, and it highlights the sine function factor oscillation of pump waves. Note that the introduced V _xc2DH affects the electron density distribution by modifying the value of $\tilde{D}$ , determined through the Wigner–Seitz radius r_s thus by the equilibrium density n ₀. The second term of the reduced perturbation electron density is induced by the electric field of the pump wave in equation (14). Because the pump wave satisfies the relation k ₀ = 1/λ, the second term of the reduced perturbation electron density is modulated by the wavelength λ, the amplitude Φ₀, and the time phase θ. Furthermore, the stopping power S(v) determined by wake force can be given by the value of the gradient of the induced potential Φ_ind [Reference Zhang, Jiang and Yi32]:

(16)

\begin{matrix} S (v) = e Z_{1} {\frac{\partial Φ_{ind}}{\partial x} |}_{z = z_{0}}^{r = v t} = n_{0} \frac{z_{1}^{2} e^{4}}{m_{e}} \int_{- \infty}^{+ \infty} \frac{i k_{x}}{D (k, ω)} e^{- k z_{0} - k z} e^{i k \cdot (r - v t)} d k + \frac{2 π z_{1} e^{3} n_{0} k_{0}^{2} Φ_{0} e^{- k_{0} z}}{m_{e} {(D_{0} (k_{0}, ω_{0}))}^{2}} (Im (D_{0} (k_{0}, ω_{0})) \cos (k_{0} x - ω_{0} t) - Re (D_{0} (k_{0}, ω_{0})) \sin (k_{0} x - ω_{0} t)) . \end{matrix}

In the calculation above, the residue theorem is used. Using dimensionless variables, equation (16) can be reduced to

(17)

\begin{matrix} \frac{S (v)}{S_{0}} = A_{s} \int_{- \infty}^{+ \infty} d k_{y} \int_{- \infty}^{+ \infty} d ω B_{s} Im ({\tilde{D}}^{- 1}) + C_{s} (Im (\tilde{D_{0}} (k_{0}, ω_{0})) \cos (k_{0} (v - c) + θ) - Re (\tilde{D_{0}} (k_{0}, ω_{0})) \sin (k_{0} (v - c) + θ)), \end{matrix}

where

(18)

\begin{matrix} S_{0} = \frac{{(Z_{1} e / a_{B})}^{2}}{2 π}, \\ A_{s} = \frac{(v_{B}^{2})}{(2 π r_{s}^{3} v^{2})}, \\ B_{s} = - ω e^{- 2 k z_{0}}, \\ C_{s} = \frac{Φ_{0} k_{0}^{2} e^{- k_{0} z_{0}}}{(r_{s}^{3} {(\tilde{D_{0}})}^{2})} . \end{matrix}

The stopping power is related to the velocity v of the incident particle. When the velocity variable v is extracted from the sine and cosine functions, the time variable t is also exposed. Similarly, to study the influence of the time variable t on the sine and cosine functions, we introduce the time phase θ, which means that the time variable t is changing the sine and cosine functions with $- γ_{0} t$ in equation (16) while the time phase θ is changing the sine and cosine functions through θ in equation (17). Note that the second term of the stopping power is induced by the excited electric field of the pump wave in equation (16). The pump wave satisfies the relation ω₀ = k ₀c and 3k ₀ = 1/λ, and hence, the second term of the stopping power is modulated by the wavelength λ, the amplitude Φ₀, and the time phase θ.

In the next section, we will give the result of the effect of V _xc2DH on the stopping power of the incident particle. To investigate the effect of pump wave on the perturbation electron density of 2DQEG and the stopping power of the incident particle, we calculate the perturbation electron density and the stopping power under different wavelength λ, amplitude Φ₀, and time phase θ. In what follows, we take r_s = 1, the collision frequency γ = 0.02ω_P, incident particle as a proton Z₁ = 1, and the distance between the incident particle and 2DQEG z ₀ = 9λ_f, while the wavelength λ, the amplitude Φ₀, and the time phase θ of the pump wave are varied.

3. Result and Discussion

3.1. Two-Dimensional Electron Exchange-Correlation Potential of High Density

In the past, Li et al. and Zhang et al. have done a lot of theoretical research studies on the stopping power of incident particles in the study of the charged particles moving above two-dimensional electron gases [Reference Li, Song and Wang4, Reference Li, Song and Wang5, Reference Zhang, Jiang and Yi32, Reference Zhang, Song and Wang38, Reference Zhang, Song, Zhao and Wang43]. On the basis of the work of Li et al. and Zhang et al., we investigate the effect of exchange-correlation potential on the stopping power of incident particles. In this paper, the stopping power of incident particles calculated without considering the exchange-correlation potential is consistent with the research results of Li et al. and Zhang et al.

As it can be assumed, V _xc2DH modified the stopping power of the incident particle because V _xc2DH has been considered in the momentum-balance equation, thereby changing the value of D and D ₀ in equation (16) and modifying the stopping power of the incident particle. In order to verify the necessity of introducing V _xc2DH to modifying QHD, Figure 3 presents the stopping power of the incident particle versus particle velocity under different conditions.

Figure 3: The stopping power of incident particles, respectively, with V _xc3D = 0.1362, V _xc2DH = −0.6002, and V _xc2DL = −0.114. Here, Z ₁ = 1, r_s = 1, γ = 0.02, and Z ₀ = 9.

Comparing the curve of stopping power without exchange-correlation potential (solid line), the curve of stopping power with V _xc3D = 0.1362 (dashed line), the curve of stopping power with V _xc2DL = −0.114 (dotted line), and the curve of stopping power with V _xc2DH = −0.6002 (dash-dotted line) are shown in Figure 3. Considering the peak of the stopping power of the incident particle near V = 3.5V_B in Figure 3, the value of stopping power without exchange-correlation potential exhibits an obvious difference with respect to stopping power featured by V _{xc3 D}, V _xc2DL, and V _xc2DH. The value of stopping power without exchange-correlation potential is 0.01517S ₀. The values of stopping power with V _xc3D, V _xc2DL, and V _xc2DH, respectively, are 0.01506S ₀, 0.01525S ₀, and 0.01559S ₀, which, respectively, are changed by 0.725%, 0.527%, and 2.77% relative to stopping power without V_xc (these are dimensionless values). The results show that the two-dimensional exchange-correlation potential, especially V _xc2DH, influences the stopping power of incident particles. V _xc2DH is adapted throughout the paper within subsequent results.

3.2. Amplitude of Pump Wave

In order to study the modulation of pump wave amplitude Φ₀ on the perturbation electron density and the stopping power, Figures 4 and 5 show the spatial distribution of the perturbation electron density as a function of x and the stopping power of the incident particle versus particle velocity under three different intensity conditions of the pump wave, where the wavelength is λ = 1.67 × 10⁻⁸ cm. The comparison of the perturbation electron density without the pump wave and with the pump wave for $Φ_{0} = 1 \times 10^{2} e / λ_{f}$ , $Φ_{0} = 2 \times 10^{2} e / λ_{f}$ , and $Φ_{0} = 3 \times 10^{2} e / λ_{f}$ , respectively, shows that the oscillation amplitude of the perturbation electron density increases when considering the pump wave in most cases and the oscillation amplitude of the perturbation electron density becomes larger as the amplitude of the pump wave increases in Figure 4. This is because the value of the second term in equation (14) is proportional to the factor of Φ₀. But the oscillation amplitude of the perturbation electron density displays specific features near x = −29λ_f. This is due to the value of the incident particle external potential and induced potential by the perturbation electron density. The comparison of the stopping power without the pump wave and with the pump wave for $Φ_{0} = 1 \times 10^{2} e / λ_{f}$ , $Φ_{0} = 2 \times 10^{2} e / λ_{f}$ , and $Φ_{0} = 3 \times 10^{2} e / λ_{f}$ , respectively, shows that the pump wave yields obvious modulation effects on the stopping power, and the intensity of this modulation is proportional to the amplitude of the pump wave in Figure 5. This is because the value of the second term in equation (17) is proportional to Φ₀. In the low-velocity region (v < 8.3v_B), the stopping power that peaks at v = 3.5v_B with the pump wave for $Φ_{0} = 1 \times 10^{2} e / λ_{f}$ , $Φ_{0} = 2 \times 10^{2} e / λ_{f}$ , and $Φ_{0} = 3 \times 10^{2} e / λ_{f}$ , respectively, increases by 12.4%, 24.7%, and 37.1% relative to the case without the pump wave. Moreover, in the high-velocity region (v > 8.3v_B), it is easy to see that the pump wave dominates the stopping power as the particle velocity increases compared to the stopping power quickly vanishing without the pump wave. The value of stopping power at v = 13.7v_B with pump wave for $Φ_{0} = 3 \times 10^{2} e / λ_{f}$ decreases to −0.003772S ₀. It is surprising that a negative stopping power [Reference Nersisyan and Akopyan34] is found due to the existence of the pump wave, different from our observation without the pump wave. The appearance of negative stopping power means that the particle drains energy from the 2DQEG.

Figure 4: 2DQEG perturbation electron density. Here, Z ₁ = 1, r_s = 1, v = 2, γ = 0.02, Z ₀ = 9, λ = 1.67 × 10⁻⁸ cm, and θ = 0.

Figure 5: Stopping power of the incident charged particle. Here, Z ₁ = 1, r_s = 1, γ = 0.02, Z ₀ = 9, λ = 1.67 × 10⁻⁸ cm, and θ = 0.

Generally, when the pump wave does not exist, due to the disturbance of the electron density in 2DQEG caused by charged particles, negative charges will accumulate behind the charged particles in 2DQEG to excite the wake field, thereby reducing the speed of the charged particles. When the pump wave exists, in addition to the density disturbance of 2DQEG caused by charged particles, we also need to consider the density disturbance caused by the electric field component of the pump wave. Since the electric field of the pump wave changes periodically with time, the density disturbance caused by the electric field component of the pump wave will also change periodically. The deceleration and acceleration effects of the electric field excited by the perturbation density caused by the pump wave on the charged particles will change periodically with time. When the electric field excited by the perturbation density caused by the pump wave in 2DQEG has a deceleration effect on the charged particles, the stopping power of the charged particles will increase. When the electric field excited by the perturbation density caused by the pump wave accelerates the charged particles, the stopping power of the charged particles will decrease. When the acceleration of the electric field excited by the perturbation density caused by the pump wave has a stronger effect on the charged particles than the deceleration of wake field excited by the charged particles, the velocity of the charged particles will increase and negative stopping power appears. In other words, the polarity of wake field can be tuned in the presence of the pump wave, yielding the negative stopping power.

In the case of λ = 1.56 × 10⁻⁷ cm, the comparison of the perturbation electron density and stopping power without the pump wave and with the pump wave shows that the modulation effect on the oscillation amplitude of the perturbation electron density and stopping power is proportional to the amplitude of the pump wave in Figures 6 and 7. As the value of the second term in equations (14) and (17) increases with increase in Φ₀, it is easy to see that the oscillation amplitude of the perturbation electron density at x = −5λ_f with the pump wave for $Φ_{0} = 1 \times 10^{2} e / λ_{f}$ , $Φ_{0} = 2 \times 10^{2} e / λ_{f}$ , and $Φ_{0} = 3 \times 10^{2} e / λ_{f}$ , respectively, increases by 22.2%, 44.5%, and 66.8% in Figure 6. At x > 10λ_f or x < −46λ_f, the pump wave with wavelength λ = 1.56 × 10⁻⁷ cm dominates the oscillation amplitude of the perturbation electron density compared to the perturbation electron density quickly vanishing without the pump wave. The stopping power at v = 3.5v_B with the pump wave for Φ₀ = 1 × 10²e/λ_f, $Φ_{0} = 2 \times 10^{2} e / λ_{f}$ , and $Φ_{0} = 3 \times 10^{2} e / λ_{f}$ , respectively, increases to 0.03991S ₀, 0.06422S ₀, and 0.08854S ₀ in Figure 7.

Figure 6: 2DQEG perturbation electron density. Here, Z ₁ = 1, r_s = 1, v = 2, γ = 0.02, Z ₀ = 9, λ = 1.56 × 10⁻⁷ cm, and θ = 0.

Figure 7: Stopping power of the incident charged particle. Here, Z ₁ = 1, r_s = 1, γ = 0.02, Z ₀ = 9, λ = 1.56 × 10⁻⁷ cm, and θ = 0.

In the case of λ = 1.76 × 10⁻⁶ cm, the comparison of the perturbation electron density and stopping power without the pump wave and with the pump wave shows that the pump wave has an obvious modulation effect on the oscillation amplitude of the perturbation electron density and stopping power in Figures 8 and 9, and the modulation effect of pump wave strengthens as the amplitude of pump wave increases, which again arises from equations (14) and (17) second term proportional to Φ₀. The modulation effect of the pump wave only changes the amplitude, not the shape of the perturbation electron density and stopping power. It is easy to see that the oscillation amplitude of the perturbation electron density at x = −5λ_f with the pump wave for $Φ_{0} = 1 \times 10^{2} e / λ_{f}$ , $Φ_{0} = 2 \times 10^{2} e / λ_{f}$ , and $Φ_{0} = 3 \times 10^{2} e / λ_{f}$ , respectively, increases by 22.9%, 45.8%, and 68.7% in Figure 8. The oscillation amplitude of the perturbation electron density increases with the pump wave amplitude. The amplitude of the stopping power at v = 3.5v_B with pump wave for $Φ_{0} = 1 \times 10^{2} e / λ_{f}$ , $Φ_{0} = 2 \times 10^{2} e / λ_{f}$ , and $Φ_{0} = 3 \times 10^{2} e / λ_{f}$ , respectively, decreases to 0.004182S ₀, −0.0007226S ₀, and −0.01862S ₀ in Figure 9. The oscillation amplitude of the stopping power rises as the amplitude of the pump wave increases. Interestingly, the negative stopping power is also found in Figure 9.

Figure 8: 2DQEG perturbation electron density. Here, $Z_{1} = 1, r_{s} = 1, v = 2, γ = 0.02, Z_{0} = 9, λ = 1.76 \times 10^{- 6}$ cm, and θ = 0.

Figure 9: Stopping power of the incident charged particle. Here, $Z_{1} = 1, r_{s} = 1, γ = 0.02, Z_{0} = 9, λ = 1.76 \times 10^{- 6}$ cm, and θ = 0.

Figures 10–15 show the perturbation electron density and stopping power with the pump wave for λ = 1.76 × 10⁻⁴ cm, λ = 1.76 × 10⁻² cm, and λ = 1.76 × 10⁻¹ cm. The results of Figures 12 and 14 are similar to that of Figure 10. And the results of Figures 13 and 15 are similar to that of Figure 11 because the variations of the second term in equations (14) and (17) remain small in the domain of $- 120 λ_{f} < x > 80 λ_{f}$ and $0 < v > 20 v_{B}$ . The oscillation amplitude of the perturbation electron density increases as the amplitude of pump wave rises in Figures 10, 12, and 14. The effect of the pump wave on the stopping power strengthens as the amplitude of pump wave increases in Figures 11, 13, and 15.

Figure 10: 2DQEG perturbation electron density. Here, $Z_{1} = 1, r_{s} = 1, v = 2, γ = 0.02, Z_{0} = 9, λ = 1.76 \times 10^{- 4}$ cm, and θ = 0.

Figure 11: Stopping power of the incident charged particle. Here, $Z_{1} = 1, r_{s} = 1, γ = 0.02, Z_{0} = 9, λ = 1.76 \times 10^{- 4}$ cm, and θ = 0.

Figure 12: 2DQEG perturbation electron density. Here, $Z_{1} = 1, r_{s} = 1, v = 2, γ = 0.02, Z_{0} = 9, λ = 1.76 \times 10^{- 2}$ cm, and θ = 0.

Figure 13: Stopping power of the incident charged particle. Here, $Z_{1} = 1, r_{s} = 1, γ = 0.02, Z_{0} = 9, λ = 1.76 \times 10^{- 2}$ cm, and θ = 0.

Figure 14: 2DQEG perturbation electron density. Here, $Z_{1} = 1, r_{s} = 1, v = 2, γ = 0.02, Z_{0} = 9, λ = 1.76 \times 10^{- 1}$ cm, and θ = 0.

Figure 15: Stopping power of the incident charged particle. Here, $Z_{1} = 1, r_{s} = 1, γ = 0.02, Z_{0} = 9, λ = 1.76 \times 10^{- 1}$ cm, and θ = 0.

Therefore, the amplitude of the pump wave has a significant effect on the oscillation amplitude of the perturbation electron density and stopping power, while the oscillation amplitude of the perturbation electron density and stopping power rises with the pump wave amplitude. The amplitudes $Φ_{0} \leq 5 \times 10^{5} e / λ_{f}$ are realistic values in current experimental technologies and consistent with those used in [Reference Zhang, Zhai and Yi51].

3.3. Wavelength of the Pump Wave

The wavelength of the pump wave is an essential parameter for the modulation of the pump wave on the perturbation electron density and the stopping power. Comparing Figures 14 and 15, it is found that the influence of wavelength of the pump wave on the perturbation electron density and the stopping power is mainly reflected in the modulation of the shape of the perturbation electron density curve and the stopping power curve.

Let us consider the perturbation electron density curve and the stopping power curve for a short wavelength of the pump wave of λ = 1.67 × 10⁻⁸ cm. In the space range from −120λ _f to 80λ _f, the perturbation electron density shows a violent periodic oscillation, contrasted to the no pump wave case, as shown in Figure 4. The stopping power also has been modulated by the pump wave in the velocity range from 0v_B to 20v_B in Figure 5. The reason is that the second term contains the factor k ₀ in equations (14) and (17), which is inversely proportional to the wavelength λ. So the perturbation electron density and the stopping power are strongly modulated by the pump wave.

Observe the perturbation electron density curve and the stopping power curve with a wavelength of λ = 1.56 × 10⁻⁷ cm. The number of oscillating periods of the perturbation electron density decreases over a range of −120λ _f to 80λ _f in Figure 6. The oscillation period of the stopping power is smaller than one period in the velocity range from 0v_B to 20v_B in Figure 7. This is because the k ₀ value decreases in equations (14) and (17) as the wavelength λ is increasing. Therefore, the effect of the pump wave on the perturbation electron density and stopping power modulation is reduced compared to the wavelength λ = 1.67 × 10⁻⁸ cm.

With the wavelength of the pump wave λ = 1.76 × 10⁻⁶ cm, the perturbation electron density oscillates less than one period in the space range of −120λ _f to 80λ _f in Figure 8. The stopping power oscillation period is also shorter than one period in the velocity range from 0v_B to 20v_B in Figure 9. As the wavelength λ becomes longer, the value of k ₀ becomes smaller. The pump wave modulation to perturbation electron density and stopping power is further reduced in the region of $- 120 λ_{f} < x < 80 λ_{f}$ and $0 < v < 20 v_{B}$ .

When the amplitude of the pump wave is $Φ_{0} = 3 \times 10^{2} e / λ_{f}$ , the maximum value of perturbation electron density curves with wavelength λ = 1.67 × 10⁻⁸ cm, λ = 1.56 × 10⁻⁷ cm, and λ = 1.76 × 10⁻⁶ cm, respectively, is 0.2082n ₀, 0.2469n ₀, and 0.2497n ₀ in Figures 4, 6, and 8. Those results show that the maximum value of the perturbation electron density increases in the region of $- 120 λ_{f} < x < 80 λ_{f}$ with the wavelength which means that the wavelength of the pump wave not only modifies the shape of the curve but also the amplitude of the perturbation electron density. Comparing the stopping power curves when the amplitude of pump wave is $Φ_{0} = 3 \times 10^{2} e / λ_{f}$ in Figures 5, 7, and 9, the maximum value of stopping power with wavelength λ = 1.67 × 10⁻⁸ cm, λ = 1.56 × 10⁻⁷ cm, and λ = 1.76 × 10⁻⁶ cm, respectively, is 0.02137S ₀, 0.08854S ₀, and −0.01863S ₀. The value of stopping power dramatically changes because k ₀ varied when the wavelength is increasing.

As the wavelength of the pump wave increases, the pump wave modulates the shape of the perturbation electron density curve more and more weakly in the space range of −120λ_f to 80λ_f in Figures 10, 12, and 14. Similarly, as the wavelength of the pump wave increases, the shape of the stopping power curve is modulated more weakly by the pump waves in the 0v_B to 20v_B speed range in Figures 11, 13, and 15. The modulation of the wavelength of the pump wave to the shape of the perturbation electron density curve and the stopping power curve decreases as the wavelength increases, arising from the tiny second term change in equations (14) and (17).

3.4. Phase of Pump Wave

Finally, in order to study the effect of time phase θ on the perturbation electron density and stopping power, Figures 16–19 exhibit the perturbation electron density and stopping power with wavelength λ = 1.67 × 10⁻⁸ cm and λ = 1.76 × 10⁻⁴ cm for θ = 0 (solid line), θ = π/2 (dashed line), θ = π (dotted line), and $θ = 3 π / 2$ (dash-dotted line). Here, $Φ_{0} = 1 \times 10^{2} e / λ_{f}$ . It is easy to see that the amplitude of perturbation electron density and stopping power oscillates periodically with different time phase θ in Figures 16–19. This is because the second term includes the parameter time phase θ in equations (14) and (17). Therefore, the time phase θ obviously impacts the perturbation electron density and stopping power.

Figure 16: 2DQEG perturbation electron density. Here, Z ₁ = 1, r_s = 1, v = 2, γ = 0.02, Z ₀ = 9, λ = 1.67 × 10⁻⁸ cm, and $Φ_{0} = 1 \times 10^{2} e / λ_{f}$ .

Figure 17: Stopping power of the incident charged particle. Here, Z ₁ = 1, r_s = 1, γ = 0.02, Z ₀ = 9, λ = 1.67 × 10⁻⁸ cm, and $Φ_{0} = 1 \times 10^{2} e / λ_{f}$ .

Figure 18: 2DQEG perturbation electron density. Here, Z ₁ = 1, r_s = 1, v = 2, γ = 0.02, Z ₀ = 9, λ = 1.76 × 10⁻⁴ cm, and $Φ_{0} = 1 \times 10^{2} e / λ_{f}$ .

Figure 19: Stopping power of the incident charged particle. Here, Z ₁ = 1, r_s = 1, γ = 0.02, Z ₀ = 9, λ = 1.76 × 10⁻⁴ cm, and $Φ_{0} = 1 \times 10^{2} e / λ_{f}$ .

It is worth pointing out that the negative stopping power in Figures 5, 9, and 19 we found is instantaneous, which is difficult to verify experimentally, and it is also difficult to apply. So we consider averaging the stopping power with respect to the time phase θ. After averaging the stopping power with respect to the time phase θ from 0 to 2π, we found that the pump wave with $Φ_{0} = 1 \times 10^{2} e / λ_{f}$ decreases the value of stopping power in Figure 20. This is because the value of averaged second term is negative in equation (17) with respect to the time phase θ. But amplitude reduction shows up as a relatively small one in connection to the too low pump wave intensity (Figure 20). This is because the intensity of pump wave is too low. Linearized QHD is valid only for low intensity of pump wave. The comparison of the stopping power without the pump wave and the average stopping power with pump wave for $Φ_{0} = 2 \times 10^{10} e / λ_{f}$ and λ = 1.76 × 10⁻⁴ cm demonstrates that the average stopping power decreases in the region of $0 < v < 20 V_{B}$ in Figure 20.

Figure 20: The average value of stopping power for time phase. Here, Z ₁ = 1, r_s = 1, γ = 0.02, Z ₀ = 9, and λ = 1.76 × 10⁻⁴ cm.

It is an exciting result that the average stopping power of the incident particles is negative in the region $v < 2.3 v_{B}$ and $v > 5.8 v_{B}$ when λ = 1.76 × 10⁻⁴ cm and $Φ_{0} = 2 \times 10^{10} e / λ_{f}$ . The incident particles will gain energy and can be accelerated within the velocity range of $v < 2.3 v_{B}$ and $v > 5.8 v_{B}$ . Although this result is not strictly obtained by using the linearized QHD due to the amplitude of pump wave which is too large, this could predict qualitatively that the time averaged stopping power has indeed been affected by the pump waves.

In conclusion, the time phase of pump wave θ can adjust the value of perturbation electron density effectively and stopping power, and the averaging of the stopping power with respect to the time phase θ also has been impacted. Therefore, one can see that the perturbation electron density and stopping power depend on λ, Φ₀, and θ, indicating that the pump wave using different parameters can modulate the perturbation electron density and stopping power of 2DQEG.

4. Summary

In this paper, considering the electron exchange-correlation potential V _xc2DH in QHD, we have calculated the perturbation electron density and stopping power with the modulation of pump wave for a charged particle with a constant distance and speed moving above 2DQEG. In our calculation, we compare the stopping power with and without V _xc2DH and show that it is necessary to consider V _xc2DH in 2DQEG. Then, we compare the perturbation electron density and stopping power with different wavelength λ, amplitude Φ₀, and time phase θ of the pump wave and show that it is meaningful to investigate the modulation effect of the pump wave on the perturbation electron density and the stopping power. The wavelength λ, the amplitude Φ₀, and time phase θ of the pump wave can effectively modulate the value of the perturbation electron density and stopping power. The temporal averaged stopping power with respect to the time phase θ also has been influenced by pump waves.

The possible experimental equipment to verify the prediction includes three components: pump wave, charged particle beam, and 2DQEG. The pump wave can be provided by laser or microwave device, the plasma wake field [Reference Yao, Li and Cao52] can produce the charged particle beam, and 2DQEG has already been observed in experiments [Reference Qi and Zhang2, Reference Di Pietro, Ortolani and Limaj3]. All of them can be realized with present technology. At present, the intensity of laser systems provided by the laboratory has arrived 10²⁴ ∼ 10²⁶Wcm⁻² [Reference Kemp, Fiuza and Debayle53]. Thus, the parameter of the pump wave considered here can be easily satisfied. We can verify the negative stopping power by measuring the velocity of the incident particle. For instance, charged particles move in the x direction parallel to the 2DQEG plane, with a velocity of v and a height of z ₀ away from the 2DQEG, and the pump wave irradiates the 2DQEG along the +x direction. The wavelength of pump wave is λ = 1.76 × 10⁻⁴ cm and $Φ_{0} = 2 \times 10^{10} e / λ_{f}$ . The density of two-dimensional electron gases is n ₀ = 5.69 × 10¹⁵cm⁻². The velocity of incident particles is v = 2v_B (v_B = 2.18 × 10⁸ cm/s), and the distance between incident particles and two-dimensional electron gases is z ₀ = 9λ_f = 4.76 × 10⁻⁸cm. The measured particles flying out the right side of the 2DQEG will gain energy from the 2DQEG and get accelerated.

It is appropriate to investigate the stopping power modulation by low-intensity pump waves of the charged particle moving above two-dimensional electron gases using the linearized QHD, but linearized QHD is not strictly applicable to high-intensity pump waves. In future work, we attempt to further study the effect of high-intensity pump waves on the stopping power with nonlinear QHD.

Data Availability

The data that support the findings of this study are available from the corresponding author upon reasonable request.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This work was supported by the National Natural Science Foundation of China (11975174 and 11775090), the National Magnetic Confinement Fusion Energy Research Project (2017YFE0301805 and 2017YFE0301803), and the Fundamental Research Funds for the Central Universities (WUT: 2020IB023 and 2018IB011).

References

Abrahams, E., Kravchenko, S. V., and Sarachik, M. P., “Metallic behavior and related phenomena in two dimensions,” Reviews of Modern Physics, vol. 73, no. 2, p. 251, 2001.10.1103/RevModPhys.73.251CrossRef Google Scholar

Qi, X.-L. and Zhang, S.-C., “Topological insulators and superconductors,” Reviews of Modern Physics, vol. 83, no. 4, p. 1057, 2011.10.1103/RevModPhys.83.1057CrossRef Google Scholar

Di Pietro, P., Ortolani, M., Limaj, O. et al., “Observation of Dirac plasmons in a topological insulator,” Nature Nanotechnology, vol. 8, no. 8, p. 556, 2013.10.1038/nnano.2013.134CrossRef Google Scholar

Li, C.-Z., Song, Y.-H., and Wang, Y.-N., “Electrical properties of epitaxial SrTiO3 tunnel barriers on (001) Pt/SrTiO3 substrates,” Physical Review A, vol. 79, Article ID 062903, 2009.Google Scholar

Li, C.-Z., Song, Y.-H., and Wang, Y.-N., “Wake effects and stopping power for a charged particle moving above two-dimensional quantum electron gases,” Physics Letters A, vol. 372, no. 24, p. 4500, 2008.10.1016/j.physleta.2008.04.034CrossRef Google Scholar

Wang, Y.-N. and Ma, T.-C., “Consistent calculation of the stopping power for slow ions in two-dimensional electron gases,” Physical Review A, vol. 55, no. 3, p. 2087, 1997.10.1103/PhysRevA.55.2087CrossRef Google Scholar

Mišković, Z. L., Segui, S., Gervasoni, J. L., and Arista, N. R., “Plasmons in graphene traversed by relativistic electronic beam with arbitrary trajectories,” Physics Letters B, vol. 94, Article ID 125414, 2016.Google Scholar

Mišković, Z. L., Liu, W.-K., Goodman, F., and Wang, Y.-N., “Spatial distribution of ion charges in fast, partially stripped clusters traversing solid targets,” Nuclear Instruments and Methods in Physics Research Section B: Beam Interactions with Materials and Atoms, vol. 193, p. 221, 2002.10.1016/S0168-583X(02)00754-1CrossRef Google Scholar

Maynard, G. and Deutsch, C., “Energy loss and straggling of ions with any velocity in dense plasmas at any temperature,” Physical Review A, vol. 26, no. 1, p. 665, 1982.10.1103/PhysRevA.26.665CrossRef Google Scholar

Horing, N. J. M., Tso, H. C., and Gumbs, G., “Fast-particle energy loss in the vicinity of a two-dimensional plasma,” Physical Review B, vol. 36, no. 3, p. 1588, 1987.10.1103/PhysRevB.36.1588CrossRef Google Scholar PubMed

Bret, A. and Deutsch, C., “Dielectric response function and stopping power of a two-dimensional electron gas,” Physical Review E, vol. 48, no. 4, p. 2994, 1993.10.1103/PhysRevE.48.2994CrossRef Google Scholar PubMed

Bret, A. and Deutsch, C., “Ion stopping in two-dimensional electron layers,” Europhysics Letters (EPL), vol. 25, no. 4, p. 291, 1994.10.1209/0295-5075/25/4/009CrossRef Google Scholar

Lindhard, J. and Rensen, A. H., “Relativistic theory of stopping for heavy ions,” Physical Review A, vol. 53, no. 4, p. 2443, 1996.10.1103/PhysRevA.53.2443CrossRef Google Scholar PubMed

Zwicknagel, G., Toepffer, C., and Reinhard, P.-G., “Stopping of heavy ions in plasmas at strong coupling,” Physics Reports, vol. 309, no. 3, p. 117, 1999.10.1016/S0370-1573(98)00056-8CrossRef Google Scholar

Custidiano, E. R. and Jakas, M. M., “Classical-trajectory Monte Carlo calculations of the electronic stopping cross section for keV protons and antiprotons impinging on hydrogen atoms,” Physical Review A, vol. 72, Article ID 022708, 2005.10.1103/PhysRevA.72.022708CrossRef Google Scholar

Zaremba, E., Nagy, I., and Echenique, P. M., “Nonlinear screening and stopping power in two-dimensional electron gases,” Physical Review B, vol. 71, Article ID 125323, 2005.10.1103/PhysRevB.71.125323CrossRef Google Scholar

Patel, P. K., Mackinnon, A. J., Key, M. H. et al., “Sequence dependent DNA-mediated conduction,” Physical Review Letters, vol. 91, Article ID 125004, 2003.Google Scholar

Chen, S. N., Atzeni, S., Gauthier, M. et al., “Nuclear instruments and methods in physics research section A: accelerators, spectrometers, detectors and associated equipment,” Proceedings of the First European Advanced Accelerator Concepts Workshop, vol. 740, p. 105, 2013.Google Scholar

Zylstra, A., Frenje, J., Grabowski, P. et al., “Non-equilibrium thermodynamics of gravitational screens,” Physical Review Letters, vol. 114, Article ID 215002, 2015.Google Scholar

Frenje, J., Grabowski, P., Li, C. et al., “Negative tension branes as stable thin shell wormholes,” Physical Review Letters, vol. 115, Article ID 205001, 2015.Google Scholar

Hayes, A. C., Jungman, G., Schulz, A. E. et al., “Reaction-in-flight neutrons as a test of stopping power in degenerate plasmas,” Physics of Plasmas, vol. 22, Article ID 082703, 2015.10.1063/1.4928104CrossRef Google Scholar

Gus’kov, S. Y., Solyanikova, M. S., and Korneev, P. A., “Specifics of powerful shock initialization by energetic ion beam,” Plasma Physics and Controlled Fusion, vol. 61, Article ID 045006, 2019.10.1088/1361-6587/ab00daCrossRef Google Scholar

Jacoby, J., Hoffmann, D. H. H., Laux, W. et al., “Stopping of Heavy Ions in a Hydrogen Plasma,” Physical Review Letters, vol. 74, no. 9, p. 1550, 1995.10.1103/PhysRevLett.74.1550CrossRef Google Scholar

Cambridge, R., “LPB volume 36 issue 4 Cover and Front matter,” Laser and Particle Beams, vol. 36, pp. 98–104, 2018.Google Scholar

Jieru Ren, Z. D., “Pitfalls in identifying active catalyst species,” Nature Communications, vol. 11, p. 5157, 2020.Google Scholar

Hoffmann, D. H. H., Tahir, N. A., Udrea, S. et al., “High energy density physics with heavy ion beams and related interaction phenomena,” Contributions to Plasma Physics, vol. 50, no. 1, p. 7, 2010.10.1002/ctpp.201010004CrossRef Google Scholar

Hoffmann, D. H. H., Fortov, V. E., Kuster, M. et al., “High energy density physics generated by intense heavy ion beams,” Astrophysics and Space Science, vol. 322, no. 1–4, p. 167, 2009.10.1007/s10509-009-0001-2CrossRef Google Scholar

Frank, A., Blažević, A., Grande, P. L. et al., “Phase field modeling of the evolution of partical interface shape distribution during coarsening,” Physical Review E, vol. 81, Article ID 026401, 2010.Google Scholar

Raphael, E. and Gennes, P.-G. d., “Effective mass of a charged particle travelling above a dielectric fluid surface,” Europhysics Letters (EPL), vol. 31, no. 5-6, p. 293, 1995.10.1209/0295-5075/31/5-6/008CrossRef Google Scholar

Lafleur, T. and Baalrud, S. D., “Transverse force induced by a magnetized wake,” Plasma Physics and Controlled Fusion, vol. 61, no. 12, Article ID 125004, 2019.10.1088/1361-6587/ab45d4CrossRef Google Scholar

Lafleur, T. and Baalrud, S. D., “Friction in a strongly magnetized neutral plasma,” Plasma Physics and Controlled Fusion, vol. 62, no. 9, Article ID 095003, 2020.10.1088/1361-6587/ab9beaCrossRef Google Scholar

Zhang, Y., Jiang, W., and Yi, L., “Stopping power of two-dimensional spin quantum electron gases,” Nuclear Instruments and Methods in Physics Research Section B: Beam Interactions with Materials and Atoms, vol. 349, p. 72, 2015.10.1016/j.nimb.2015.02.016CrossRef Google Scholar

Arista, N. R., Galvão, R. M. O., and Miranda, L. C. M., “Influence of a strong laser field on the stopping power for charged test particles in nondegenerate plasmas,” Journal of the Physical Society of Japan, vol. 59, p. 544, 1990.10.1143/JPSJ.59.544CrossRef Google Scholar

Nersisyan, H. B. and Akopyan, E. A., “Stopping and acceleration effect of protons in a plasma in the presence of an intense radiation field,” Physics Letters A, vol. 6, 1999.Google Scholar

Hu, Z.-H., Song, Y.-H., Mišković, Z. L., and Wang, Y.-N., “Energy dissipation of ion beam in two-component plasma in the presence of laser irradiation,” Laser and Particle Beams, vol. 29, no. 3, p. 299, 2011.10.1017/S0263034611000334CrossRef Google Scholar

Zhang, Y., Song, Y.-H., and Wang, Y.-N., “Nonlinear wake potential and stopping power for charged particles interacting with a one-dimensional electron gas,” Physics of Plasmas, vol. 18, no. 11, Article ID 112705, 2011.Google Scholar

Hu, Z.-H., Song, Y.-H., and Wang, Y.-N., “Critical properties of the S4 model for Sierpinski carpet,” Physical Review E, vol. 85, Article ID 016402, 2012.Google Scholar

Zhang, Y., Song, Y.-H., and Wang, Y.-N., “Simulations of interactions of high-energy proton beam with high dense matter based on two-dimensional quantum hydrodynamic model,” Laser and Particle Beams, vol. 31, no. 2, p. 345, 2013.10.1017/S0263034613000086CrossRef Google Scholar

Wu, D., Luan, S. X., Wang, J. W. et al., “Plasma detachment in linear devices,” Plasma Physics and Controlled Fusion, vol. 59, Article ID 065004, 2017.Google Scholar

Mir, Z., Shahid, M., Jamil, M., Rasheed, A., and Shahbaz, A., “Nonlinear beat excitation of low frequency wave in degenerate plasmas,” Physics Plasmas, vol. 8, 2018.Google Scholar

Stern, F. and Howard, W. E., “Properties of semiconductor surface inversion layers in the electric quantum limit,” Physical Review, vol. 163, no. 3, p. 816, 1967.10.1103/PhysRev.163.816CrossRef Google Scholar

Zhang, Y. and Jiang, W., “Pseudomagnetic field modulation of stopping power for a charged particle moving above graphene,” Physics of Plasmas, vol. 25, Article ID 072107, 2018.10.1063/1.5039588CrossRef Google Scholar

Zhang, Y., Song, Y.-H., Zhao, Y.-T., and Wang, Y.-N., “Two-dimensional quantum hydrodynamic model for the heating of a solid target using a Gaussian cluster,” Laser and Particle Beams, vol. 30, no. 4, pp. 671–677, 2012.10.1017/S0263034612000559CrossRef Google Scholar

Kumar, P., Singh, S., and Ahmad, N., “Conversion efficiency of even harmonics of whistler pulse in quantum magnetoplasma,” Laser and Particle Beams, vol. 37, no. 1, pp. 5–11, 2019.10.1017/S0263034619000089CrossRef Google Scholar

Kumar, P. and Ahmad, N., “Surface plasma wave in spin-polarized semiconductor quantum plasma,” Laser and Particle Beams, vol. 38, no. 2, pp. 159–164, 2020.10.1017/S026303462000018XCrossRef Google Scholar

Zhang, L.-Y., Zhao, X.-Y., Qi, X., Duan, W.-S., Xiao, G. Q., and Yang, L., “Collective effects on the wakefield and stopping power of an ion beam pulse in plasmas,” Physics of Plasmas, vol. 22, Article ID 053109, 2015.10.1063/1.4921133CrossRef Google Scholar

Haas, F., “Kinetic theory derivation of exchange-correlation in quantum plasma hydrodynamics,” Plasma Physics and Controlled Fusion, vol. 61, Article ID 044001, 2019.10.1088/1361-6587/aaffe1CrossRef Google Scholar

Pollack, L. and Perdew, J. P., “Evaluating density functional performance for the quasi-two-dimensional electron gas,” Journal of Physics: Condensed Matter, vol. 12, no. 7, pp. 1239–1252, 2000.Google Scholar

Stern, F. and Das Sarma, S., “Electron energy levels in GaAs-Ga1−xAlxAsheterojunctions,” Physical Review B, vol. 30, no. 2, p. 840, 1984.10.1103/PhysRevB.30.840CrossRef Google Scholar

Brey, L., Dempsey, J., Johnson, N. F., and Halperin, B. I., “Infrared optical absorption in imperfect parabolic quantum wells,” Physical Review B, vol. 42, no. 2, p. 1240, 1990.10.1103/PhysRevB.42.1240CrossRef Google Scholar PubMed

Zhang, Y., Zhai, F., and Yi, L., “Study of spin-polarized plasma driven by spin force in a two-dimensional quantum electron gas,” Physics Letters A, vol. 380, no. 46, p. 3908, 2016.10.1016/j.physleta.2016.09.050CrossRef Google Scholar

Yao, W., Li, B., Cao, L. et al., “Generation of monoenergetic proton beams by a combined scheme with an overdense hydrocarbon target and an underdense plasma gas irradiated by ultra-intense laser pulse,” Laser and Particle Beams, vol. 32, no. 4, p. 583, 2014.10.1017/S0263034614000561CrossRef Google Scholar

Kemp, A. J., Fiuza, F., Debayle, A. et al., “Laser–plasma interactions for fast ignition,” Nuclear Fusion, vol. 54, Article ID 054002, 2014.10.1088/0029-5515/54/5/054002CrossRef Google Scholar

Figure 1: The schematic diagram of the interaction between pump waves, charged particles, and 2DQEG.

Figure 2: The change of electron exchange-correlation potential with Wigner–Seitz radius of 2DQEG.

Figure 3: The stopping power of incident particles, respectively, with Vxc3D = 0.1362, Vxc2DH = −0.6002, and Vxc2DL = −0.114. Here, Z1 = 1, rs = 1, γ = 0.02, and Z0 = 9.

Figure 4: 2DQEG perturbation electron density. Here, Z1 = 1, rs = 1, v = 2, γ = 0.02, Z0 = 9, λ = 1.67 × 10−8 cm, and θ = 0.

Figure 5: Stopping power of the incident charged particle. Here, Z1 = 1, rs = 1, γ = 0.02, Z0 = 9, λ = 1.67 × 10−8 cm, and θ = 0.

Figure 6: 2DQEG perturbation electron density. Here, Z1 = 1, rs = 1, v = 2, γ = 0.02, Z0 = 9, λ = 1.56 × 10−7 cm, and θ = 0.

Figure 7: Stopping power of the incident charged particle. Here, Z1 = 1, rs = 1, γ = 0.02, Z0 = 9, λ = 1.56 × 10−7 cm, and θ = 0.

Figure 8: 2DQEG perturbation electron density. Here, Z1=1,rs=1,v=2,γ=0.02,Z0=9,λ=1.76×10−6 cm, and θ = 0.

Figure 9: Stopping power of the incident charged particle. Here, Z1=1,rs=1,γ=0.02,Z0=9,λ=1.76×10−6 cm, and θ = 0.

Figure 10: 2DQEG perturbation electron density. Here, Z1=1,rs=1,v=2,γ=0.02,Z0=9,λ=1.76×10−4 cm, and θ = 0.

Figure 11: Stopping power of the incident charged particle. Here, Z1=1,rs=1,γ=0.02,Z0=9,λ=1.76×10−4 cm, and θ = 0.

Figure 12: 2DQEG perturbation electron density. Here, Z1=1,rs=1,v=2,γ=0.02,Z0=9,λ=1.76×10−2 cm, and θ = 0.

Figure 13: Stopping power of the incident charged particle. Here, Z1=1,rs=1,γ=0.02,Z0=9,λ=1.76×10−2 cm, and θ = 0.

Figure 14: 2DQEG perturbation electron density. Here, Z1=1,rs=1,v=2,γ=0.02,Z0=9,λ=1.76×10−1 cm, and θ = 0.

Figure 15: Stopping power of the incident charged particle. Here, Z1=1,rs=1,γ=0.02,Z0=9,λ=1.76×10−1 cm, and θ = 0.

Figure 16: 2DQEG perturbation electron density. Here, Z1 = 1, rs = 1, v = 2, γ = 0.02, Z0 = 9, λ = 1.67 × 10−8 cm, and Φ0=1×102e/λf.

Figure 17: Stopping power of the incident charged particle. Here, Z1 = 1, rs = 1, γ = 0.02, Z0 = 9, λ = 1.67 × 10−8 cm, and Φ0=1×102e/λf.

Figure 18: 2DQEG perturbation electron density. Here, Z1 = 1, rs = 1, v = 2, γ = 0.02, Z0 = 9, λ = 1.76 × 10−4 cm, and Φ0=1×102e/λf.

Figure 19: Stopping power of the incident charged particle. Here, Z1 = 1, rs = 1, γ = 0.02, Z0 = 9, λ = 1.76 × 10−4 cm, and Φ0=1×102e/λf.

Figure 20: The average value of stopping power for time phase. Here, Z1 = 1, rs = 1, γ = 0.02, Z0 = 9, and λ = 1.76 × 10−4 cm.

Article contents

Stopping Power Modulation by Pump Waves of Charged Particles Moving above Two-Dimensional Electron Gases

Abstract

1. Introduction

2. Quantum Hydrodynamic Model

2.1. Physic Model

2.2. Mathematical Model

2.3. Exchange-Correlation Potential

2.4. Linearized QHD Model

3. Result and Discussion

3.1. Two-Dimensional Electron Exchange-Correlation Potential of High Density

3.2. Amplitude of Pump Wave

3.3. Wavelength of the Pump Wave

3.4. Phase of Pump Wave

4. Summary

Data Availability

Conflicts of Interest

Acknowledgments

References

Save article to Kindle

Save article to Dropbox

Save article to Google Drive

Reply to: Submit a response

Your details

You have entered the maximum number of contributors

Conflicting interests