The supersymmetry (SUSY), applied to treat fermions and bosons in a unified way in elementary particle physics since the concept first arose in 1971 by Ramond, Golfand and Likhtman, has been researched extensively during past four decades  - . The starting point of SUSY is the supersymmetric versions of well known KdV equation first by Kupershmidt in 1984 (a simple fermionic but not supersymmetric extension) and later found independently of the work of Manin-Radual on super KP hierarchy    . It was pointed out afterwards that the latter is indeed a truly N = 1 (N refers to the number of supersymmetries, for N = 1 standard, for N > 1 extended) sKdV equation which is invariant under supersymmetric transformation . Since then, various properties have been established for its supersymmetric versions, such as Lax representation, bi-Hamiltonian structures, Backlund transformation (BT), Painleve analysis, N-soliton solutions, (non) local conserved quantities, etc.  - .
The integrability of sKdV equation can be established in the way to supersymmetrize the unique Gardner transformation. For the well known KdV equation
we extend the classical spacetime to a super-spacetime , where is a Grassmannian odd number . Now we write the N = 1 supersymmetric KdV (sKdV) equation accompanied with a fermionic super-variable under the compact form
Here the covariant super-derivative D is defined by . Mathieu found that a unique extension
of establishing the integrability maps a solution of the superfield Gardner equation
into a solution of the sKdV equation   . The Gardner equation is also called the extended KdV equation with the variable-sign cubic non-linear term or the combined KdV and mKdV (KdV-mKdV) equation. It is widely used in various branches of physics, such as plasma physics, fluid physics, nonlinear phenomena and quantum field theory, etc., and it also describes a variety of wave phenomena in plasma and solid state     .
This map was also used to recover an infinite number of conservation laws for the sKdV equation, and construct interesting BT . It is easy to show that such super equation is invariant under the supersymmetry transformation: , ( is an anticommuting parameter.). The component form of the above equation with the superfield can be rewritten as
Note that and u are new setted fermionic and bosonic functions, the usual Gardner equation is recovered by setting the fermionic variable to be absent and the sKdV equation is the limiting case where .
Nonlinear partial differential equations play an important role in nonlinear physics, even nonlinear science. Various effective methods have been proposed to derive explicit or formal solutions. Recently, a simple but powerful bosonization approach, which main idea is to consider fields in a Grassmannian algebra and rewrite a system in a basis of this algebra to arrive at a system of ordinary (commutative) evolution equations, can effectively simplify such systems containing anti-communicating fermionic fields   . In , B. Ren et al. used this approach in the N = 1 supersymmetric Burgers (SB) system, and the exact solutions of the usual pure bosonic systems are obtained with the mapping and deformation method and Lie point symmetries theory. In , the Lie point symmetries of the supersymmetric KdV-a system are considered and similarity reductions of it are conducted. Several types of similarity reduction solutions of the coupled bosonic equations are also simply obtained. The motivation and purpose of this paper is to show this procedure and outcome of the method by taking the superfield Gardner equation and to acquire novel traveling wave solutions of this equation.
The rest of this paper is organized as follows. Section 2 and 3 are brief reminders of fairly basic illustrations of the bosonization approach of the superfield Gardner equation with two and three fermionic parameters. In Section 4 we present the N fermionic parameters bosonization case. In Section 5 we give the N = 1 supersymmetric KdV and mKdV equations using parallel procedure in two particular cases, and we will also give a short summary.
2. Two-Fermionic-Parameter Bosonization
To get direct comprehension and fixed notations of the bosonizaion approach for superfield Gardner equation with multi-fermionic parameters, we first concentrate on a linear space . Mathematically, such uncomplicated method used for vanishing fermionic fields is based on direct sum of superspace , where and represents subspace containing even and odd elements, respectively . Here and in the following we omit the exterior algebra sign and denote it briefly by ordinary multiplication. Moreover, relevant solutions may involve rich (super) symmetries physically. An exact example is the cases of two and three fermionic parameters bosonization, thereby we can directly derive the usual Gardner equation and coupled equations appear below. For the case of two fermionic parameters and with , let the two component fields u and be expanded as
, , (6)
here , are all usual bosonic functions with respect to spacetime variable x and t, thus we get nonlinear Partial Differential Equations (PDEs) in the component form by using (5)
where and .
Next we introduce the traveling wave variable along with the constants of wavenumber k, angular frequency and phase , therefore, above equations would be changed to a system consisting of ordinary differential equations(ODEs):
Note that we denote in here.
The traveling waves we discuss are only in the usual spacetime but not in the super-spacetime , for example, with Grassmannian constant is different from those in the usual spacetime. In addition to a directly integrable ODE in , the solution of the residual system which are related to third-order linear (non)homogeneous ODEs in , and , can be obtained through the variable transformation from ordinary coordinates space to phase space on the base of periodic wave solutions of usual Gardner equation. We first solve out from Equation (8a), and the result reads
where, and are two arbitrary integral constants, here and below the new variable z stands for function.
To get the mapping relations of, and, we introduce the variable transformations as follows, and . Applying the transformation via Equation (8a), the linear ODEs (8b)-(8c) are reduced to mapping and deformation relations between the traveling wave solutions of the classical Gardner equation and its supersymmetric equation by exploiting the known solutions of classical Gardner equation
where linear operators read
and the nonhomogeneous term is
while is just an arbitrary integral constant. On this basis, the mapping and deformation relations are obtained as
where and are arbitrary constants, and with auxiliary function . Coefficients for used are defined as , and.
Thus, we have constructed the general two-fermionic-parameter traveling wave solutions of the supersymmetric version of Gardner system
with the known solution or Equation (8a) of the usual Gardner equation which have been extensively studied in many literatures  . For a special case which takes to eliminate nonhomogeneous terms, the above traveling wave solution becomes
It is interesting to see that the expression (9) is a trivial type of the symmetries or conservation quantity of standard Gardner equation ;,. In fact, for any given being a solution of the usual Gardner equation, a certain type of solutions of the bosonic-looking equation can be constructed as follows:, , , where represents any symmetry of the usual Gardner equation, and we have much freedom to choose such that it can construct solutions without restricting to the traveling wave type solutions and infinitely many symmetries of the superfield Gardner equations.
It is clear that the solution (9) can be expressed by the form of the Jacobian elliptic sine functions, i.e.,
where, modulus m and constants are related to other known constants through following relations:
Therefore, we derive special type solutions of the superfield Gardner equation:
and are all arbitrary integral constants, the Jacobian elliptic functions are denoted as and, the incomplete elliptic integrals are denoted as and. We also note that for the critical modulus m = 1, above solutions read:
3. Three-Fermionic-Parameter Bosonization
For the case of three fermionic parameters with, let the two component fields u and be expressed as
, ( 22)
here and are all usual bosonic functions with respect to spacetime variable x and t. The symbol is the third-order Levi-Civita tensor. Then from (5), we get PDEs in component form as
where the somewhat complex nonhomogeneous terms read
Introducing the traveling wave variable with constants k, and and variable transformations and, similar to the previous case, the above bosonization system becomes following ODEs
where nonhomogeneous terms are
with arbitrary integral constants. Mapping and deformation relations with respect to via traveling wave variable are constructed as
where and are some arbitrary integral
constants, and functioin with integral constant.
Therefore, we have obtained the three-fermionic-parameter traveling wave solutions of the superfield Gardner system. While one of the Grassmann numbers tends to zero, the solution turns back to above section. Similar to the two-fermionic-parameter case, we write a special type solution
where are some constants, are arbitrary symmetries of the usual Gardner equation, and is an arbitrary solution of the usual Gardner equation.
4. N-Fermionic-Parameter Bosonization
Motivated by above sections, we repeat same procedure to get traveling wave solutions of the superfield Gardner equation via bosonizaion approach with N fermionic parameters here. The component fields u and can be expressed as
Here and below we denote by the set of multi-indices which satisfies, and the upper bound of
summations. The elements and are bosonic smooth real or complex valued functions defined in the commutative algebra. Thus the super Gardner model (5) is transformed to a new pure bosonic-looking system with coupled nonlinear PDEs:
Operators related to Gâteaux derivative of Equation (4) or Equations (5) with an operator decomposition read
and nonhomogeneous terms are
We have also used the shorthand notations
and for even case; for odd case. The parity function is defined in . The symbol is the generalized Levi-Civita anti-symmetric tensor in m dimensions. Equation (30a) apparently suggests that it is the standard Gardner equation. Each solution in the case of fermionic parameters is a portion of solutions of the system consisting of n homologous Grassmannian variables. Introducing the traveling wave variable and variable transformations, mapping transformations are arranged as follows:
1) Equation of: (33a)
2) Equations of: (33b)
3) Equations of: (33c)
Here the nonhonogenerous terms are, each usual bosonic functions and are represented by and, respectively.
Finally, the mapping and deformation relations of and with N fermionic parameters can be formally rewritten as
in which auxiliary functions are
, and are some arbitrary integral constants.
5. Discussion and Summary
With relationships at hand, in case of preserving the square and cubic non-linear terms in Equations (5) or taking limitation, one can conveniently acquire the corresponding result of the integrable N = 1 sKdV regardless of some integral constants. Similarly to Section 4, the general traveling wave solutions of the N = 1 sKdV Equation (2) with N fermionic parameters are
Due to the case of taking limitation and the absence of non-linear quadratic term in Equation (4), traveling wave solutions of the usual integrable N = 1 smKdV Equation  
with N-fermionic-parametric Bosonization procedure can be derived as
In summary, the bosonization approach with multi-fermionic parameters to deal with supersymmetric systems is developed in the super Gardner equation with the role of traveling wave solution. The procedure and technique are also available for N = 1 sKdV and smKdV equations derived from two particular cases. We expect this procedure exhibited in our paper could be successfully applied or formulated in the N = 1 supersymmetric sine-Gordon equation, especially in the N = 2 version of KdV (SKdVa) equations     . For
example, in the case of two fermionic parameters, letting, , in bosonic field, the
simple traveling wave solutions of the N = 1 supersymmetric Sine-Gordon equation are
where and are the usual super derivatives, S represents and, E, , , and are integral constants, is a solution of the standard sine-Gordon equation.
We thank the Editor and the referee for their comments. This research is funded by the National Natural Science Foundation of China (Grant nos. 31570364) and the National Key Research and Development Program of China (Grant No. 2018YFE0103700). This support is greatly appreciated.
 Tian, K. (2010) Symmetry Approach to N=1 Supersymmetric Systems. China-Japan Joint Workshop on Integrable Systems on the Occasion of Prof. Ryogo Hirota’s 77th Birthday, Shaoxing, 7-10 January 2010.
 Liu, Q.P. and Hu, X.-B. (2005) Bilinearization of N=1 Supersymmetric Korteweg-de Vries Equation Revisited. Journal of Physics A: Mathematical and General, 38, 6371-6378.
 Carstea, A.S., Ramani, A. and Grammaticos, B. (2001) Constructing the Soliton Solutions for the N=1 Supersymmetric KdV Hierarchy. Nonlinearity, 14, 1419-1423.
 Nakoulima, O., Zahibo, N., Pelinovsky, E., Talipova, T., Slunyaev, A. and Kurkin, A. (2004) Analytical and Numerical Studies of the Variable-Coefficient Gardner Equation. Applied Mathematics and Computation, 152, 449-471.
 Zhang, L.-H., Dong, L.-H. and Yan, L.-M. (2008) Construction of Non-Travelling Wave Solutions for the Generalized Variable-Coefficient Gardner Equation. Applied Mathematics and Computation, 203, 784-791.
 Lu, D.C. and Liu, C.L. (2010) A Sub-ODE Method for Generalized Gardner and BBM Equation with Nonlinear Terms of Any Order. Applied Mathematics and Computation, 217, 1404-1407.
 Andrea, S., Restuccia, A. and Sotomayor, A. (2001) An Operator Valued Extension of the Super Korteweg Cde Vries Equations. Journal of Mathematical Physics, 42, 2625-2634.
 Ren, B., Gao, X.-N., Yu, J. and Liu, P. (2015) Exact Solutions of Supersymmetric Burgers Equation with Bosonization Procedure. Open Mathematics, 13, 502-509.
 Gao, X.N., Yang, X.-D. and Lou, S.-Y. (2012) Exact Solutions of Supersymmetric KdV-a System via Bosonization Approach. Communications in Theoretical Physics. Communications in Theoretical Physics, 58, 617-622.
 Fan, E.G. and Hon, Y.C. (2010) Quasiperiodic Wave Solutions of N=2 Supersymmetric KdV Equation in Superspace. Studies in Applied Mathematics, 125, 343-371.
 Bai, C.-L., Bai, C.-J. and Zhao, H. (2005) A New Generalized Algebraic Method and Its Application in Nonlinear Evolution Equations with Variable Coefficients. Zeitschrift für Naturforschung, 60, 211-220.
 Grammaticos, B., Ramani, A. and Carstea, A.S. (2001) Bilinearization and Soliton Solutions of the N=1 Supersymmetric Sine-Gordon Equation. Journal of Physics A: Mathematical and General, 34, 4881-4886.
 Popowicz, Z. (1998) Integrable Extensions of N=2 Supersymmetric KdV Hierarchy Associated with the Nonuniqueness of the Roots of the Lax Operator. Physics Letters A, 249, 204-208.