We replace the old PACIAE 2.1 and PACIAE 2.2 program series with new ones, respectively. In the new program series, the method of creation of the partonic initial state is changed, an option for calling PYEVNT or PYEVNW in PYTHIA is added, and the few bugs in the calculation for the reduction of strange quark suppression are corrected. New version program summaryProgram Title: PACIAE version 2.1 & 2.2Program Files doi:http://dx.doi.org/10.17632/w3g68dj4d9.1Licensing provisions: CC By 4.0Programming language: FORTRAN 77Journal reference of previous version: Comput. Phys. Comm. 184 (2013) 1476, 193 (2015) 89.Does the new version supersede the previous version?: YesNature of problem: PYEVNT is the subroutine in PYTHIA to administer the generation of a high-pT event via calls a number of subroutines. PYEVNW is the new subroutine in PYTHIA to administer the generation of an event for new multiple interactions scenario. In previous PACIAE versions only PYEVNT is called, while both PYEVNT and PYEVNW can be called in this new version.Solution method: There has been added a parameter of mstp81 in PACIAE programs and input files of usux.dat (for 21a or 22a) and usu.dat (for 21b, 21c ,22b and 22c). The PACIAE will call PYEVNT when mstp81=1( default) and call PYEVNW when mstp81=21.Reasons for the new version: The method of creation of the partonic initial state is changed, an option for calling PYEVNT or PYEVNW in PYTHIA is added, and the few bugs in the calculation for the reduction of strange quark suppression are corrected in the new version.Summary of revisions: The PACIAE model [1] is a parton and hadron cascade model based on PYTHIA [2]. PACIAE model consists of four stages of the parton initiation, the parton evolution (rescattering), the hadronization, and the hadron evolution (rescattering). The stage of parton initiation is just by means of PYTHIA model and the rest are modeled by us [1]. PYTHIA model is for the high energy elementary particle collisions such as hadron–hadron and lepton–hadron collisions, but PACIAE model is also for high energy lepton–nucleus, hadron–nucleus and nucleus–nucleus collisions.So far, the PACIAE model has three series of PACIAE 2.0 (catalogue identifier: AEKI_v1_0) [1], PACIAE 2.1 (AEKI_v2_0) [3], and PACIAE 2.2 (AEKI_v2_2) [4]. The discrepancy between PACIAE 2.0 and PACIAE 2.1 is just the sampling of momentum x and y components at fixed pT. In the PACIAE 2.0 it is sampled on the circle with radius of pT, but on the circumference of ellipse with half major and minor axes of pT(1+δp) and pT(1−δp), respectively in PACIAE 2.1. Including the reaction of lepton–nucleon and lepton–nucleus is the distinction of PACIAE 2.2 series from the PACIAE 2.1. Each series of PACIAE 2.0, PACIAE 2.1 and PACIAE 2.2 consists of a, b, and c packages. The package of PACIAE 2.1a, for instance, is for the high energy elementary collisions, while the PACIAE 2.1b as well as PACIAE 2.1c for the high energy hadron–nucleus and nucleus–nucleus collisions. We refer to [1] for the difference between PACIAE 2.0b and PACIAE 2.0c (PACIAE 2.1b and PACIAE 2.1c as well as PACIAE 2.2b and PACIAE 2.2c).In this replacement the modifications are as follows: 1.The method of creation of a partonic initial state is changed from setting mstj(1)=0 (or mstp(111)=0) before “CALL PYINIT” to adding a “RETURN” statement before the statement of “CALL PYEXEC” in subroutine PYEVNT (PYEVNW) in p21b.f.2.Correcting the mistakes in the calculation for the mechanism of reduction of strange quark suppression.3.Calling PYEVNT or PYEVNW is controlled by mstp(81) (=1 or 21) instead of calling PYEVNT only originally.