  A key feature of integrable systems is that they can be solved to obtain exact analytical solutions. We show how new models can be constructed through generalisations of some well known nonlinear partial differential equations with PT-symmetries whilst preserving integrability. Subsequently, we develop new methods from well-known ones to obtain exact analytical soliton solutions for these new systems. The first PT-symmetric generalization we take are extensions to the complex and multicomplex fields. In agreement with the reality property present in PT-symmetric non-Hermitian quantum systems, we find PT-symmetries also play a key role in the reality of conserved charges here. We then extend our investigations to explore degenerate multi-soliton solutions for the sine-Gordon and Hirota equations. In particular, we find the usual time-delays from degenerate soliton solution scattering are time-dependent, unlike the non-degenerate multi-soliton solutions, and provide a universal formula to compute the exact time-delay values for the scattering of N-soliton solutions. Other PT-symmetric extensions of integrable systems we take are of nonlocal nature, with nonlocalities in space and/or in time, of time crystal type. Whilst developing new methods for the construction of soliton solutions for these systems, we find new types of solutions with different parameter dependence and qualitative behaviour even in the one-soliton solution cases. We exploit gauge equivalence between the Hirota system with continuous Heisenberg and Landau-Lifschitz systems to see how nonlocality is inherited from one system to another and vice versa. Extending investigations to the quantum regime, we generalize the scheme of Darboux transformations for fully time-dependent non-Hermitian quantum systems, which allows us to create an infinite tower of solvable models. 
Nonlinear Classical and Quantum Integrable Systems with PT-symmetries