Bound states in the continuum generated by supersymmetric quantum mechanics and phase rigidity of the corresponding wavefunctions

Supersymmetric quantum mechanics (SUSYQM) is a method that can be used for generating complex potentials with entirely real spectrum with bound states in the continuum (BIC). These complex potentials are isospectral with the initial one, but SUSYQM method adds discrete BIC's at selected energies. Corresponding wavefunctions created by SUSYQM are biorthogonal and complex, hence we can discuss their phase rigidity and illustrate the application of SUSYQM on the examples of three specific potential profiles (free electron, negative Dirac potential and quantum well with infinite walls).