We employ a functional integral technique to calculate the macroscopic quantum tunneling rate at zero temperature of a current-biased Josephson junction weakly coupled to a resonator. We allow for the effects of environmental dissipation on both the junction and the resonator via baths of harmonic oscillators, and we consider both cases of weak and strong junction damping. The resonator can be, in principle, any quantum harmonic oscillator that couples to the junction bilinearly in either the coordinate or the velocity. The Josephson phase difference $\ensuremath{\phi}$ is the tunneling variable in this system, and the low-lying bound states in the junction's potential energy well are weakly metastable due to tunneling through a finite-width energy barrier. There has been some interest in using such biased junctions as qubits with the resonator providing multiqubit coupling. The main result of this work is that coupling to the resonator has a suppressive effect on the junction's tunneling: the stronger the coupling strength between the junction and resonator, the greater the reduction of the tunneling rate. Including damping to the junction also suppresses tunneling, a well-known result, but damping of the resonator actually reduces the magnitude of suppression; i.e., damping the resonator partially counteracts the suppression provided directly by the junction-resonator coupling. Details of the junction-resonator coupling yield interesting variations on this theme that may be useful for qubit design. For example, the coupling
Uint between a current-biased junction and an AlN dilatational resonator should depend on the angular frequency of the resonator ${\ensuremath{\omega}}_{R}$ in a power-law fashion [A. N. Cleland and M. R. Geller, Phys. Rev. Lett. 93, 070501 (2004)], specifically ${U}_{\mathit{int}}\ensuremath{\propto}{({\ensuremath{\omega}}_{R})}^{n}$, where
n=1∕2. We find that for any value of
n in the range
0<n<1 the junction's tunneling rate exhibits a nonmonotonic dependence on ${\ensuremath{\omega}}_{R}$, and for
n=1∕2, the tunneling rate is maximally suppressed for ${\ensuremath{\omega}}_{R}∕{\ensuremath{\omega}}_{J}\ensuremath{\approx}1$, where ${\ensuremath{\omega}}_{J}$ is the bias current-dependent plasma frequency of the junction.
The WKB approach to quantum tunneling from a metastable well is generalized to include a weak oscillatory force. Both one-dimensional decay and multidimensional decay, appropriate for a system coupled to a dissipative environment, are considered. In both cases, general expressions are derived for the enhancement of the decay rate linear in the applied force. The enhancement is expressed exclusively in terms of the instanton trajectories of the unperturbed system. For a cubic potential with no damping, the decay rate exhibits sharp resonances when the applied frequency coincides with the energy spacing between the ground and excited states in the metastable well. In the strong-damping limit, this resonance structure is completely suppressed.