Eccentricity dynamics of wide binaries -- II. The effect of stellar encounters and constraints on formation channels

GAIA wide stellar binaries (separations $\sim 10^3-10^{4.5}$ AU) are observed to have a superthermal eccentricity distribution function (DF), well-fit by $P(e) \propto e^\alpha$ with $\alpha \sim 1.2$. In Modak \& Hamilton (2023), we proved that this DF cannot have been produced by Galactic tidal torques starting from any realistic DF that was not already superthermal. Here, we consider the other major dynamical effect on wide binaries: encounters with passing stars. We derive and solve the Fokker-Planck equation governing the evolution of binaries in semimajor axis and eccentricity under many weak, impulsive, penetrative stellar encounters. We show analytically that these encounters drive the eccentricity DF towards thermal on the same timescale as they drive the semimajor axes $a$ towards disruption, $t_\mathrm{ion} \sim 4\,\mathrm{Gyr}\,(a/10^4\,\mathrm{AU})^{-1}$. We conclude that the observed superthermal DF must derive from an even more superthermal (i.e. higher $\alpha$) birth distribution. This requirement places strong constraints on the dominant binary formation channels. A testable prediction of our theory is that $\alpha$ should be a monotonically decreasing function of binary age.