Order parameter quantum fluctuations in a two-dimensional system of mesoscopic Josephson junctions

The boson lattice Hubbard model is used to study the role of quantum fluctuations of the phase and local density of the superfluid component in establishing a global superconducting state for a system of mesoscopic Josephson junctions or grains. The quantum Monte Carlo method is used to calculate the density of the superfluid component and fluctuations in the number of particles at sites of the two-dimensional lattice for various average site occupation numbers n0 (i.e., number of Cooper pairs per grain). For a system of strongly interacting bosons, the phase boundary of the ordered superconducting state lies above the corresponding boundary for its quasiclassical limit—the quantum XY-model—and approaches the latter as n0 increases. When the boson interaction is weak in the boson Hubbard model (i.e., the quantum fluctuations of the phase are small), the relative fluctuations of the order parameter modulus are significant when n0<10, while quantum fluctuations in the phase are significant when n0<8; this determines the region of mesoscopic behavior of the system. Comparison of the results of numerical modeling with theoretical calculations show that mean-field theory yields a qualitatively correct estimate of the difference between the phase diagrams of the quantum XY-model and the Hubbard model. For a quantitative estimate of this difference the free energy and thermodynamic averages of the Hubbard model are expanded in powers of 1/n0 using the method of functional integration.