Path integral Monte Carlo simulation of global and local superfluidity in liquid $^{4}$He reservoirs separated by nanoscale apertures

We present a path integral Monte Carlo study of the global superfluid fraction and local superfluid density in cylindrically-symmetric reservoirs of liquid $^{4}$He separated by nanoaperture arrays. The superfluid response to both translations along the axis of symmetry (longitudinal response) and rotations about the cylinder axis (transverse response) are computed, together with radial and axial density distributions that reveal the microscopic inhomogeneity arising from the combined effects of the confining external potential and the $^4$He-$^4$He interatomic potentials. We make a microscopic determination of the length-scale of decay of superfluidity at the radial boundaries of the system by analyzing the local superfluid density distribution to extract a displacement length that quantifies the superfluid mass displacement away from the boundary. We find that the longitudinal superfluid response is reduced in reservoirs separated by a septum containing sufficiently small apertures compared to a cylinder with no intervening aperture array, for all temperatures below $T_{\lambda}$. For a single aperture in the septum, a significant drop in the longitudinal superfluid response is seen when the aperture diameter is made smaller than twice the empirical temperature-dependent $^4$He healing length, consistent with the formation of a weak link between the reservoirs. Increasing the diameter of a single aperture or the number of apertures in the array results in an increase of the superfluid density toward the expected bulk value.