“…In the present work, we numerically solve the Heisenberg's equation of motion derived within the pseudospin formalism 29,30 in order to describe the superconducting state, and we support our interpretations and results with the derivation of the nonlin-ear susceptibility by means of a diagrammatic approach 31 . Moreover, in order to treat the pseudogap phase, characterized by incoherent pairs, we extend the pseudospin model artificially adding a phase to the Cooper pairs 1 , and solving the corresponding equations to obtain the time-evolution of the order parameter and the generated nonlinear current (see Fig. 2(b),(c)).…”