Involution factorizations of Ewens random permutations

An involution is a bijection that is its own inverse. Given a permutation $σ$ of $[n],$ let $\mathsf{invol}(σ)$ denote the number of ways $σ$ can be expressed as a composition of two involutions of $[n].$ We prove that the statistic $\mathsf{invol}$ is asymptotically lognormal when the symmetric groups $\mathfrak{S}_n$ are each equipped with Ewens Sampling Formula probability measures of some fixed positive parameter $θ.$ This paper strengthens and generalizes previously determined results about the limiting distribution of $\log(\mathsf{invol})$ for uniform random permutations, i.e. the specific case of $θ= 1$. We also investigate the first two moments of $\mathsf{invol}$ itself, detailing the phase transition in asymptotic behavior at $θ= 1,$ and provide a functional refinement and a convergence rate for the Gaussian limit law which is demonstrably optimal when $θ= 1.$