Branch structure of IV-characteristics in the capacitively coupled Josephson junctions model with the diffusion current

We have solved numerically a system of dynamical equations for the gauge-invariant phase differences between superconducting layers for a stack of N intrinsic junctions and obtained a total branch structure in the current-voltage characteristics (IVC) of the stack. The coupling dependence of the branch's slopes is investigated and demonstrated that the equidistance of the branch structure in capacitively coupled Josephson junctions (CCJJ) model is broken at small values of coupling parameter. Changes in the parameters of the boundary conditions and the use of periodic boundary conditions do not affect this result. In the framework of the CCJJ model with the diffusion current we simulate an experiment and obtain the IV-characteristic with equidistant branch structure at different values of model parameters