Abstract:Two knots on a string can either be separated or intertwined, and may even pass through each other. At the microscopic scale, such transitions may occur spontaneously, driven by thermal fluctuations, and can be associated with a topological free energy barrier. In this manuscript, we study the respective location of a trefoil (3 1 ) and a figure-eight (4 1 ) knot on a semiflexible polymer, which is parameterized to model dsDNA in physiological conditions. Two cases are considered: first, end monomers are grafted to two confining walls of varying distance. Free energy profiles and transition barriers are then compared to a subset of free chains, which contain exactly one 3 1 and one 4 1 knot. For the latter, we observe a small preference to form an intertwined state, which can be associated with an effective entropic attraction. However, the respective free energy barrier is so small that we expect transition events to occur spontaneously and frequently in polymers and DNA, which are highly knotted for sufficient strain lengths.

Heterogeneous multiscale methods (HMM) combine molecular accuracy of particle-based simulations with the computational efficiency of continuum descriptions to model flow in soft matter liquids. In these schemes, molecular simulations typically pose a computational bottleneck, which we investigate in detail in this study. We find that it is preferable to simulate many small systems as opposed to a few large systems, and that a choice of a simple isokinetic thermostat is typically sufficient while thermostats such as Lowe-Andersen allow for simulations at elevated viscosity. We discuss suitable choices for time steps and finite-size effects which arise in the limit of very small simulation boxes. We also argue that if colloidal systems are considered as opposed to atomistic systems, the gap between microscopic and macroscopic simulations regarding time and length scales is significantly smaller. We propose a novel reducedorder technique for the coupling to the macroscopic solver, which allows us to approximate a non-linear stress-strain relation efficiently and thus further reduce computational effort of microscopic simulations. best attributes of both parts: the molecular accuracy with the computational efficiency of continuum models.