Knotting probabilities after a local strand passage in unknotted self-avoiding polygons

Abstract

We investigate, both theoretically and numerically, the knotting probabilities after a local strand passage is performed in an unknotted self-avoiding polygon (SAP) on the simple cubic lattice. In the polygons studied, it is assumed that two polygon segments have already been brought close together for the purpose of performing a strand passage. This restricts the polygons considered to those that contain a specific pattern called Θ at a fixed location; an unknotted polygon containing Θ is called a Θ-SAP. It is proved that the number of n-edge Θ-SAPs grows exponentially (with n) at the same rate as the total number of n-edge unknotted SAPs (those with no prespecified strand passage structure). Furthermore, it is proved that the same holds for subsets of n-edge Θ-SAPs that yield a specific after-strand-passage knot-type. Thus, the probability of a given after-strand-passage knot-type does not grow (or decay) exponentially with n. Instead, it is conjectured that these after-strand-passage knot probabilities approach, as n goes to infinity, knot-type dependent amplitude ratios lying strictly between 0 and 1. This conjecture is supported by numerical evidence from Monte Carlo data generated using a composite (aka multiple) Markov chain Monte Carlo BFACF algorithm developed to study Θ-SAPs. A new maximum likelihood method is used to estimate the critical exponents relevant to this conjecture. We also obtain strong numerical evidence that the after-strand-passage knotting probability depends on the local structure around the strand-passage site. If the local structure and the crossing sign at the strand-passage site are considered, then we observe that the more ‘compact’ the local structure, the less likely the after-strand-passage polygon is to be knotted. This trend for compactness versus knotting probability is consistent with results obtained for other strand-passage models; however, we are the first to note the influence of the crossing-sign information. We use two measures of ‘compactness’: one involves the size of a smallest polygon that contains the structure and the other is in terms of an ‘opening’ angle. The opening angle definition is consistent with one that is measurable from single molecule DNA experiments. The theoretical and numerical approaches presented here are more broadly applicable to other SAP models.