Since the establishment of stochastic differential equations (SDEs) driven by Brownian motions, many scholars have contributed to study properties of SDEs, for example, [1, 3, 16, 19] and references therein. We can observe that the stochastic systems are widely applied in many fields such as biology, chemistry, finance and economy. When studying the realistic models, it is found that the real state is not only related to the present state, but also related to the past state. Stochastic differential delay equations (SDDEs) are used to describe such systems [4, 5, 20] . Moreover, if an emergency occurs, its impact on the systems must be taken into account. For example, the sudden outbreak of the new coronavirus has a huge impact on the global economy, leading to a shock in the stock market. Hence, SDEs with jumps which take both the continuous and discontinuous random effects into consideration are studied to analyze these situations [11, 21, 28] . In this paper, we will take the delay and jumps into the cosideration, i.e. we shall study SDDEs with jumps [15, 29, 30] .
Generally speaking, the true solutions of many equations can not be calculated, so it is meaningful to investigate the numerical solutions. For instance, the explicit Euler-Maruyama (EM) schemes are very popular to approximate the true solutions [16] . However, when the coefficients grow super-linearly, Hutzenthaler et al. in [13] proved that the pth moments of the EM approximations diverge to infinity for all p ∈ [1, ∞). Thus, many implicit methods have been put forward to approximate the solutions of the equations with nonlinear growing coefficients [2, 10, 24, 27] . In addition, since the explicit schemes require less computation, some modified EM methods have also been established for nonlinear stochastic equations [14, 18, 25, 26] . In particularly, the truncated EM method was originally proposed by Mao in [22] with drift and diffusion coefficients growing super-linearly. The rate of convergence of the truncated EM method was obtained in [23] . Afterward, there are many papers to study the truncated EM method for stochastic equations whose coefficients grow super-linearly, and we refer to [7, 8, 9, 12, 17] and references therein. Additionally, there are many results on the numerical solutions for SDE with jumps and SDDEs with jumps. For example, the convergence in probability of the EM method for SDDEs with jumps was discussed in [15] . The strong convergence of the EM method for SDDEs with jumps as well as the modified split-step backward Euler method approximation to the true solution was presented in [30] . The semi-implicit Euler method for SDDEs with jumps is convergent with strong order 1/2 [29] . However, there are few papers studying the numerical solutions of the super-linear SDDEs with Poisson jumps (SDDEwPJs) whose all three coefficients might grow superlinearly. Therefore, in this paper, we will investigate the strong convergence rate of the truncated EM method for super-linear SDDEwPJs in L p (p > 0) sense. This paper is organized as follows. We will introduce some necessary notation in Section 2. The rate of convergence in L p for p ≥ 2 will be discussed in Section 3. In Section 4, the rate of convergence in L p for 0 < p < 2 will be presented. Section 5 contains an example to illustrate that our main result could cover a large class of super-linear SDDEwPJs.
Throughout this paper, unless otherwise specified, we use the following notation. If A is a vector or matrix, its transpose is denoted by A T . For x ∈ R n , |x| denotes its Euclidean norm. If A is a matrix, we let |A| = trace(A T A) be its trace norm. By A ≤ 0 and A < 0, we mean A is non-positive and negative definite, respectively. If both a, b are real numbers, then a ∧ b = min{a, b} and a ∨ b = max{a, b}. Let ⌊a⌋ denote the largest integer which does not exceed a. Let R + = [0, +∞) and τ > 0. Denote by C ([−τ, 0]; R n ) the family of continuous functions ϕ from [−τ, 0] to R n with the norm ϕ = sup −τ ≤θ≤0 |ϕ(θ)|. If H is a set, denote by 1 H its indicator function; that is, 1 H (x) = 1 if x ∈ H and 1 H (x) = 0 if x / ∈ H. Let C stand for a generic positive real constant whose value may change in different appearances.
Let Ω, F , {F t } t≥0 , P be a complete probability space with a filtration {F t } t≥0 satisfying the usual conditions (i.e. it is increasing and right continuous while F 0 contains all P-null sets). Let E denote the probability expectation with respect to P. For p > 0, L p = L p (Ω, F , {F t } t≥0 , P denotes the space of random variables X with a norm |X| p :=
Brownian motion defined on the probability space. Let N(t) be a scalar Poisson process with the compensated Poisson process N(t) = N(t) − λt, where the parameter λ > 0 is the jump intensity. Furthermore, we assume that B(t) and N(t) are independent.
In this paper, we study the truncated EM method for super-linear SDDEwPJs of the form
In order to estimate the convergence rate of the truncated EM method, we assume that the initial value ξ is γ-Hölder continuous, which is a standard constraint [9, 30] .
Assumption 2.1. There exist constantsK > 0 and γ ∈ (0, 1] such that
Now, in order to obtain the rate of convergence for the truncated EM method for (2.1) in L p (p ≥ 2) sense, we need to impose the following assumptions on coefficients.
Assumption 3.1. There exist constants K 1 > 0 and β ≥ 0 such that
for any x, y,x,ȳ ∈ R n .
By Assumption 3.1, we obtain
and
for any x, y ∈ R n . Before stating the next assumption, we need more notation. Let U denote the family of continuous functions U : R n × R n → R + such that for any b > 0, there exists a constant
for any x,x ∈ R n with |x| ∨ |x| ≤ b.
Assumption 3.2. There exist constants K 2 > 0,η > 2 and U ∈ U such that
for any x, y,x,ȳ ∈ R n . Remark 1. We use an example to illustrate the necessity of setting U(·, ·). Let f (x, y) = −5x 3 + 1 8 |y| 5 4 + 2x, g(x, y) = 1 2 |x| 3 2 + y for x, y ∈ R. We could observe that there is no K 2 satisfying
but Assumption 3.2 is satisfied. The detailed proof will be provided in Section 5.
Assumption 3.3. There exist constants K 3 > 0 andp >η > 2 such that
for any x, y ∈ R n .
By using the standard method, we could derive that the moment of the true solution is bounded, i.e. E|x(t)| q < ∞, ∀T > 0. (3.8) To our best knowledge, there are few results about the strong convergence of the superlinear SDDEwPJs. On the other hand, the truncated EM method developed in [22] is a kind of useful tool to deal with the super-linear terms. Therefore, in this paper, we will show that the truncated EM solutions will converge to the true solution in L p (p > 0) sense.
To define the truncated EM scheme, we first choose a strictly increasing continuous function ϕ(r) : R + → R + such that ϕ(r) → ∞ as r → ∞ and sup |x|∨|y|≤r (|f (x, y)| ∨ |g(x, y)|) ≤ ϕ(r), ∀r ≥ 1. For example, we could choose
for some ε ∈ (0, 1/4]. For a given step size ∆ ∈ (0, 1], define the truncated mapping π ∆ : R n → R n by π ∆ (x) = |x| ∧ ϕ −1 (α(∆)) x |x| , (3.12) where we let x |x| = 0 when x = 0. It is easy to see that the truncated mapping π ∆ maps x to itself when |x| ≤ ϕ −1 (α(∆)) and to ϕ −1 (α(∆)) x |x| when |x| ≥ ϕ −1 (α(∆)). We now define the truncated functions f ∆ (x, y) = f (π ∆ (x), π ∆ (y)), g ∆ (x, y) = g(π ∆ (x), π ∆ (y)), (3.13) for any x, y ∈ R n . By definition, we could easily find that
(3.14)
Moreover, we could obtain for any x, y ∈ R n that
If Assumption 3.3 hold. Then, for any ∆ ∈ (0, 1], x, y ∈ R n , one has that
Let us now introduce the discrete-time truncated EM numerical scheme to approximate the true solution of (2.1). For some positive integer M, we take step size ∆ = τ /M. Obviously, when we choose M sufficiently large, ∆ will become sufficiently small. Define t k = k∆ for k = −M, −M + 1, −M + 2, · · · , −1, 0, 1, 2, · · · . Set X ∆ (t k ) = ξ(t k ) for k = −M, −M + 1, −M + 2, · · · , −1, 0 and then form
. As usual, there are two kinds of the continuous-time truncated EM solutions. The first one is that:
where every sample path is not continuous. The second one is defined as follows:
which the sample path is continuous. It is easy to see
is an Itô process on t ≥ 0 with its Itô differential:
We now prepare some useful lemmas. Before stating the next lemma, we define
where ⌊·⌋ is the integer part of a number.
where cp is a positive constant which is independent of ∆.
Proof. Fix anyp ≥ 2. Then by the Hölder inequality and the Burkholder-Davis-Gundy inequality, we derive from (3.14) that
By the characteristic function's argument [6] , for ∆ ∈ (0, 1], we could get
where c 0 is a positive constant independent of ∆. Then by (3.4) we have
Thus we derive that
When 0 <p < 2, an application of Jensen's inequality yields that
We complete the proof.
Lemma 3.6. Let Assumption 3.1 and 3.3 hold, then, for any 2 < q <p, we have
Proof. For any ∆ ∈ (0, 1] and t > 0, by Itô's formula and (3.17) we derive that
We now handle A 1 , A 2 , A 3 respectively. Firstly, we could see that
and
Then we get
By Lemma 3.5, (3.10), (3.14) and Young's inequality, we have
(3.31) By (3.14), we obtain that
This, together with (3.10) and Lemma 3.5, implies
Combing (3.31) and (3.33) together, we derive that 
An application of Gronwall's inequality yields that
where C is independent of ∆. Since this inequality holds for any ∆ ∈ (0, 1], the desired result follows. The proof is therefore complete.
Lemma 3.7. Let Assumption 3.1 and 3.3 hold. Then for any ∆ ∈ (0, 1] and t ∈ [0, T ], we have
and
Proof. By Lemma 3.6 and (3.22), we get (3.38) when 2 ≤ q <p. Then for any 0 < q < 2, Jensen's inequality gives that
We complete the proof.
By using the Chebyshev inequality, Lemma 3.4 and Lemma 3.6, we can immediately have the following lemma. Let us now discuss the rate of convergence in L p (p ≥ 2) sense for the truncated EM solutions to the true solution. In order to avoid the details becoming more complicated, we only calculate the rate of convergence when p = 2.
Theorem 3.9. Let Assumption 2.1, 3.1-3.3 hold. Suppose that there exists a real number q ∈ (2,p) such that q > (1 + β)η. Then, for any ∆ ∈ (0, 1], we have
and
In particular, let
with c * = 4K 1 + |f (0, 0)| + |g(0, 0)|. Then it holds for any ∆ ∈ (0, 1] that
and
We write ρ ∆,L = ρ for simplicity. Noting that forq ∈ (2,η), we have q > (1 + β)q. By Itô's formula, for any 0 ≤ t ≤ T and ∆ ∈ (0, 1], we have
(3.49) Firstly, we estimate I 1 . Noting that
Then we have
(3.50)
By Assumptions 2.1, 3.2 and (3.5), we derive that
where c 1 = 2τ K 2K 2 + 2τ κ bK 2 , κ(s) is defined as before. As for I 12 , we have 
Then by Chebyshev's inequality, Hölder's inequality and Lemma 3.6, we get
E|e ∆ (s ∧ ρ)| 2 ds + C(ϕ −1 (α(∆))) 2β+2−q .
(3.53)
We can use the same technique to handle I 122 . By Young's inequality, Hölder's inequality, Lemma 3.6, Lemma 3.7, (3.16) and the inequality 2q/(q − 2β) ≥ 2, we obtain
(3.54)
Let us now estimate I 2 . By Assumption 3.1, we obtain
Thus, using Assumption 2.1 and Lemma 3.7, we obtain
(3.55)
Combining (3.49) -(3.55) together, one can see that
An application of Gronwall's inequality yields that
The desired results (3.44) and (3.45) follow by letting L → ∞ and using Lemma 3.8 and Lemma 3.7. In particular, by the definition of varphi, we can derive (3.47) and (3.48 ). The proof is therefore complete.
4 Convergence in L p for 0 < p < 2
In this section, we will discuss the convergence and the rate of the convergence of the truncated EM method for (2.1) in L p for 0 < p < 2. To achieve this goal, we need to impose the following assumptions on coefficients.
Assumption 4.1. There exists a positive constant K R such that
for any x, y,x,ȳ ∈ R n with |x| ∨ |y| ∨ |x| ∨ |ȳ| ≤ R.
Assumption 4.2. There exist constants K 5 > 0, K 6 ≥ 0 and σ > 2 such that
for any x, y ∈ R n .
We could get the following lemma in the similar way as Lemma 3.4 was proved. In the previous section, the jump coefficient h is linear growth, but in the assumptions 4.1 and 4.2 h is allowed to grow super-linearly. Thus, we need to truncate all the three coefficients. In the same way in Section 3, we first choose a strictly increasing continuous function ϕ(r) : R + → R + such that ϕ(r) → ∞ as r → ∞ and sup |x|∨|y|≤r (|f (x, y)| ∨ |g(x, y)| ∨ |h(x, y)|) ≤ ϕ(r), ∀r ≥ 1.
(4.4)
Choose K 0 and α : (0, 1] → (0, ∞) as in (3.10) . For a given step size ∆ ∈ (0, 1], the truncated mapping π ∆ is defined as (3.12) , and the truncated functions are define as follows:
for any x, y ∈ R n . It is easy to see that
Moreover, if Assumption 4.2 holds, then it holds for any ∆ ∈ (0, 1], x, y ∈ R n that
Since |h ∆ (x, y)| ≤ α(∆), similar to the Lemma 3.5, we have the following lemma. 
and
where cp is a positive constant which is independent of ∆. As a result,
The following lemma states that the numerical solution is bounded in mean square. Proof. Since the proof is similar to that of Lemma 3.6, we only highlight how to deal with the jump term. By Itô's formula and (4.6), we derive that, for any ∆ ∈ (0, 1] and t ∈ [0, T ],
(4.11)
Moreover, by (3.10) and Lemma 4.4, we have
(4.12) Therefore, we have
We could observe that the right-hand-side term is nondecreasing in t, hence
An application of Gronwall's inequality yields that
where C is independent of ∆. We complete the proof.
Since the boundedness of the numerical solution, the estimates of stopping times in Lemma 3.8 still hold. Now, we are going to state the convergence of the truncated EM method for SDDEwPJs in L p for 0 < p < 2. Proof. Let e ∆ (t) = x(t) − x ∆ (t) for t ≥ 0 and ∆ ∈ (0, 1]. Define ρ ∆,L = τ L ∧ τ ∆,L . We write ρ ∆,L = ρ for simplicity. Obviously,
Let δ > 0 be arbitrary. By Young's inequality, we have
Hence, 
Inserting (4.18) and (4.19) into (4.17) yields that
Let ε be arbitrary. We choose δ sufficiently small such that
and choose L sufficiently large such that
Thus,
Moreover, we could use the similar technique in the proof of Theorem 3.5 in [9] to prove that
Combing (4.15), (4.21) and (4.22) together, we have
Hence, we get the desired result (4.13). Then combing (4.13) and Lemma 4.4 yields (4.14). We complete the proof.
Next, in order to estimate the rate of the convergence at time T , we have to impose an extra condition.
Assumption 4.7. There exists a positive constant K 7 such that
for any x, y,x,ȳ ∈ R n . Here, U(·, ·) is defined as before. 25) where ρ ∆,L := τ L ∧ τ ∆,L and τ L , τ ∆,L is defined in Lemma 3.8.
Proof. Let e ∆ (t) = x(t) − x ∆ (t) for t ≥ 0 and ∆ ∈ (0, 1]. We write ρ ∆,L = ρ for simplicity. For 0 ≤ s ≤ t ∧ ρ, we observe that
Recalling the definition of f ∆ , g ∆ and h ∆ , we obtain for 0 ≤ s ≤ t ∧ ρ that
By Itô's formula and Assumption 4.7, for any 0 ≤ t ≤ T we have
(4.28)
By Assumption 2.1 and Lemma 4.4, similar to the proof of Theorem 3.9, we derive that
(4.29)
These imply that
The required assertion follows by the Gronwall inequality. Then, for any T > 0, we have that
and
Proof. We use the notation of e ∆ (t) and ρ ∆,L as before. We write ρ ∆,L = ρ for simplicity. By (4.15) and (4.20) , one can see that
for any ∆ ∈ (0, 1), L > ξ and δ > 0. 
Combing Lemma 4.4 and (4.31) together, we can derive (4.32). We complete the proof.
Remark 2. If we impose an additional condition: Assume that there exist constants
for any x, y,x,ȳ ∈ R n , we could obtain better convergence rate. But our main result can cover more equations without this condition (4.37).
In this section, we give an example to illustrate our theories. Consider the super-linear scalar SDDEwPJs dx(t) = (−5x 3 (t) + 1 8 |x(t − τ )| Therefore, Assumption 3.2 is satisfied with U(x,x) = 1 4 |x −x| 2 (|x| 2 + |x| 2 ). Noting that
x T f (x, y) +p − 1 2 |g(x, y)| 2 = x(−5x 3 + 1 8 |y| 5 4 + 2x) +p − 1 2 | 1 2 |x| 3 2 + y| 2 ≤ C(1 + |x| 2 + |y| 2 ).
(5.7)
Thus, Assumption 3.3 is satisfied as well.
Additionally, it is easy to see that sup |x|∨|y|≤r (|f (x, y)| ∨ |g(x, y)| ∨ |h(x, y)|) ≤ 5r 3 , ∀r ≥ 1. (5.8) hence, we can choose ϕ(r) = 5r 3 . This means ϕ −1 (r) = ( r 5 ) 1 3 . In order for q ≥ 1+β ε to hold, we setp = 26 such that Assumption 3.3 be satisfied. Then q ∈ ((1 + β)η,p) = q ∈ (9, 26). We choose q = 25. Moreover, let ε = 1 8 , then q ≥ 1+β ε is satisfied. So α(∆) = K 0 ∆ − 1 8 . By Theorem 3.9, we have |x −x| 2 (|x| 2 + |x| 2 ) + 11 4 |y −ȳ| 2 (|y| 2 + |ȳ| 2 ) + 5λ|x −x| 2 + 3λ|y −ȳ| 2 ≤ 12(|x −x| 2 + |y −ȳ| 2 ) − 11 4 |x −x| 2 (|x| 2 + |x| 2 ) + 11 4 |y −ȳ| 2 (|y| 2 + |ȳ| 2 ). 