Full capacitance-matrix effects in driven Josephson-junction arrays

We study the dynamic response to external currents of periodic arrays of Josephson junctions, in a resistively capacitively shunted junction (RCSJ) model, including full capacitance-matrix effects}. We define and study three different models of the capacitance matrix $C_{vec{r},vec{r}}$: Model A includes only mutual capacitances; Model B includes mutual and self capacitances, leading to exponential screening of the electrostatic fields; Model C includes a dense matrix $C_{vec{r},vec{r}}$ that is constructed approximately from superposition of an exact analytic solution for the capacitance between two disks of finite radius and thickness. In the latter case the electrostatic fields decay algebraically. For comparison, we have also evaluated the full capacitance matrix using the MIT fastcap algorithm, good for small lattices, as well as a corresponding continuum effective-medium analytic evaluation of a finite voltage disk inside a zero-potential plane. In all cases the effective $C_{vec{r},vec{r}}$ decays algebraically with distance, with different powers. We have then calculated current voltage characteristics for DC+AC currents for all models. We find that there are novel giant capacitive fractional steps in the I-Vs for Models B and C, strongly dependent on the amount of screening involved. We find that these fractional steps are quantized in units inversely proportional to the lattice sizes and depend on the properties of $C_{vec{r},vec{r}}$. We also show that the capacitive steps are not related to vortex oscillations but to localized screened phase-locking of a few rows in the lattice. The possible experimental relevance of these results is also discussed.