Stellar Escape from Globular Clusters. II. Clusters May Eat Their Own Tails

We apply for the first time the Monte Carlo star cluster modeling method to study tidal tail and stellar stream formation from globular clusters, assuming a circular orbit in a smooth Galactic potential. Approximating energetically unbound bodies (potential escapers; PEs) as collisionless enables this fast but spherically symmetric method to capture asymmetric tidal phenomena with unprecedented detail. Beyond reproducing known stream features, including epicyclic overdensities, we show how 'returning tidal tails' may form after the stream fully circumnavigates the Galaxy back to the cluster, enhancing the stream's velocity dispersion. While a realistically clumpy, time-dependent Galactic potential may disrupt such tails, they warrant scrutiny as potentially excellent constraints on the Galactic potential's history and substructure. Re-examining the escape timescale $\Delta t$ of PEs, we find new behavior related to chaotic scattering in the three-body problem; the $\Delta t$ distribution features sharp plateaus corresponding to distinct locally smooth patches of the chaotic saddle separating the phase space basins of escape. We study for the first time $\Delta t$ in an evolving cluster, finding that $\Delta t\sim(E_{\rm J}^{-0.1},E_{\rm J}^{-0.4})$ for PEs with (low, high) Jacobi energy $E_{\rm J}$, flatter than for a static cluster ($E_{\rm J}^{-2}$). Accounting for cluster mass loss and internal evolution -- and (roughly) for ongoing relaxation among PEs -- lowers the median $\Delta t$ from ${\sim}10\,$Gyr to ${\lesssim}100\,$Myr. We finally outline future improvements to escape physics in the Monte Carlo method intended to enable both the first large-parameter-space studies of tidal tail/stellar stream formation from full globular cluster simulations and detailed comparisons to stream observations.