Abstract
The deoxyribonucleic acid (DNA) dynamical equation, which emerges from the oscillator chain known as the Peyrard–Bishop (PB) model for abundant optical soliton solutions, is presented, along with a novel fractional derivative operator. The Kudryashov expansion method and the extended hyperbolic function (HF) method are used to construct novel abundant exact soliton solutions, including light, dark, and other special solutions that can be directly evaluated. These newly formed soliton solutions acquired here lead one to ask whether the analytical approach could be extended to deal with other nonlinear evolution equations with fractional space–time derivatives arising in engineering physics and nonlinear sciences. It is noted that the newly proposed methods’ performance is most reliable and efficient, and they will be used to construct new generalized expressions of exact closed-form solutions for any other NPDEs of fractional order.