Beyond the mean-field theory, a new model of the Gross-Pitaevskii equation (GPE) that describes the dynamics of Bose-Einstein condensates (BECs) is derived using an appropriate phase-imprint on the old wavefunction. This modified version of the GPE in addition to the two-body interactions term, also takes into account effects of the threebody interactions. The three-body interactions consist of a quintic term and the delayed nonlinear response of the condensate system term. Then, the modulational instability (MI) of the new GPE confined in an attractive harmonic potential is investigated. The analytical study shows that the three-body interactions destabilize more the condensate system while the external potential alleviates the instability. Numerical results confirm the theoretical predictions. Further numerical investigations of the behavior of solitons reveal that the three-body interactions enhance the appearance of solitons, increase the number of solitons generated and deeply change the lifetime of solitons. Moreover, the external potential delays the appearance of solitons. Besides, a new initial condition is introduced which enables to increase the number of solitons created and deeply affects the trail of chains of solitons generated. Moreover, the MI of a condensate without the external potential, and in a repulsive potential is also investigated.