Abstract: We investigate the coupling factor phi( mu) that quantifies the magnetic flux phi per magnetic moment mu of a point-like magnetic dipole that couples to a superconducting quantum interference device (SQUID). Representing the dipole by a tiny current-carrying (Amperian) loop, the reciprocity of mutual inductances of SQUID and Amperian loop provides an elegant way of calculating phi(mu)(r,e(mu)) vs. position r and orientation e(mu) of the dipole anywhere in space from the magnetic field B-J(r) produced by a supercurrent circulating in the SQUID loop. We use numerical simulations based on London and Ginzburg-Landau theory to calculate phi (mu) from the supercurrent density distributions in various superconducting loop geometries. We treat the far-field regime ( r greater than or similar to a= inner size of the SQUID loop) with the dipole placed on (oriented along) the symmetry axis of circular or square shaped loops. We compare expressions for phi (mu) from simple filamentary loop models with simulation results for loops with finite width w (outer size A > alpha), thickness d and London penetration depth lambda(L )and show that for thin ( d << alpha ) and narrow (w < alpha) loops the introduction of an effective loop size a(eff) in the filamentary loop-model expressions results in good agreement with simulations. For a dipole placed right in the center of the loop, simulations provide an expression phi(mu)(a,A,d,lambda(L)) that covers a wide parameter range. In the near-field regime (dipole centered at small distance z above one SQUID arm) only coupling to a single strip representing the SQUID arm has to be considered. For this case, we compare simulations with an analytical expression derived for a homogeneous current density distribution, which yields excellent agreement for lambda(L)>w,d . Moreover, we analyze the improvement of phi(mu) provided by the introduction of a narrow constriction in the SQUID arm below the magnetic dipole.

Abstract: We fabricated an asymmetric nanoscale SQUID consisting of one nanobridge weak link and one Dayem bridge weak link. The current phase relation of these particular weak links is characterized by multivaluedness and linearity. While the latter is responsible for a particular magnetic field dependence of the critical current (so-called vorticity diamonds), the former enables the possibility of different vorticity states (phase winding numbers) existing at one magnetic field value. In experiments the observed critical current value is stochastic in nature, does not necessarily coincide with the current associated with the lowest energy state and critically depends on the measurement conditions. In this paper, we unravel the origin of the observed metastability as a result of the phase dynamics happening during the freezing process and while sweeping the current. Moreover, we employ special measurement protocols to prepare the desired vorticity state and identify the (hidden) phase slip dynamics ruling the detected state of these nanodevices. In order to gain insights into the dynamics of the condensate and, more specifically the hidden phase slips, we performed time-dependent Ginzburg-Landau simulations.