“…Expressing the components d 12 and d 22 of the matrix D from the conditions (2.6) where we assume that k 12 = 0 and then changing coordinates as x 1 = z 1 + iz 2 , x 2 =x 1 , x 3 =ẋ 1 , x 4 =ẋ 2 , where i = √ −1 and the overbar denotes complex conjugation, we transform equation (2.1) to iẋ = Hx, where For example, when k 11 = k 22 , we have PH T = HP, where P = diag(1, −1, −1, 1) [26], that is, in this particular case, H is a PT -symmetric non-Hermitian Hamiltonian [8,9] that contains, for example, the model considered by Schindler et al [6]. Is H in (2.11) PT -symmetric also when k 11 = k 22 ?…”