Bimodal distribution function of a three-dimensional wormlike chain with a fixed orientation of one end

Grafted Semiflexible Nunchucks with a Magnetic Bead Attached to the Free End

Semiflexible nunchucks are block copolymers, which consist of two long blocks of high bending stiffness jointed together by a short block of low bending stiffness. Semiflexible nunchucks that consist of two DNA nanorods jointed by a short segment of double-stranded (ds) DNA and confined in two dimensions have been used in recent experiments by Fygenson and coworkers as a tool to magnify the bending fluctuations of the linking dsDNA, which in turn are used to deduce the persistence length of dsDNA. In a recent theoretical analysis, we showed that in a semiflexible nunchuck with one end grafted, the fluctuations of the position of the free end that is transverse to the grafting direction exhibit a pronounced bimodality, provided that the bending stiffness of the hinge is not very large. In this article, we theoretically analyse a grafted semiflexible nunchuck with a magnetic bead attached to its free end. We show that a transverse magnetic field induces an asymmetry in the bimodal distribution of the transverse fluctuations of the free end. This asymmetry is very sensitive to interactions with a magnetic field and, in principle, could be used in magnetometry (the measurement of a magnetic field or the magnetic moment of the bead). We also investigate how the response of the bimodal distribution of the transverse fluctuations of the free end to a magnetic field depends on the bending stiffness of the nunchuck hinge. In addition, we analyse the closely related systems of a single filament and two filaments jointed at a kink point with one end grafted and the other end attached to a magnetic bead.

Orientational Fluctuations and Bimodality in Semiflexible Nunchucks

Semiflexible nunchucks are block copolymers consisting of two long blocks with high bending rigidity jointed by a short block of lower bending stiffness. Recently, the DNA nanotube nunchuck was introduced as a simple nanoinstrument that mechanically magnifies the bending angle of short double-stranded (ds) DNA and allows its measurement in a straightforward way [Fygenson et al., Nano Lett. 2020, 20, 2, 1388-1395]. It comprises two long DNA nanotubes linked by a dsDNA segment, which acts as a hinge. The semiflexible nunchuck geometry also appears in dsDNA with a hinge defect (e.g., a quenched denaturation bubble or a nick), and in end-linked stiff filaments. In this article, we theoretically investigate various aspects of the conformations and the tensile elasticity of semiflexible nunchucks. We analytically calculate the distribution of bending fluctuations of a wormlike chain (WLC) consisting of three blocks with different bending stiffness. For a system of two weakly bending WLCs end-jointed by a rigid kink, with one end grafted, we calculate the distribution of positional fluctuations of the free end. For a system of two weakly bending WLCs end-jointed by a hinge modeled as harmonic bending spring, with one end grafted, we calculate the positional fluctuations of the free end. We show that, under certain conditions, there is a pronounced bimodality in the transverse fluctuations of the free end. For a semiflexible nunchuck under tension, under certain conditions, there is bimodality in the extension as a function of the hinge position. We also show how steric repulsion affects the bending fluctuations of a rigid-rod nunchuck.

Emergence of stationary multimodality under two-timescaled dichotomic noise

We study a linear Langevin dynamics driven by an additive non-Markovian symmetrical dichotomic noise. It is shown that when the statistics of the time intervals between noise transitions is characterized by two well differentiated timescales, the stationary distribution may develop multimodality (bi- and trimodality). The underlying effects that lead to a probability concentration in different points include intermittence and also a dynamical locking of realizations. Our results are supported by numerical simulations as well as by an exact treatment obtained from a Markovian embedding of the full dynamics, which leads to a third-order differential equation for the stationary distribution.

Bimodal probability density characterizes the elastic behavior of a semiflexible polymer in 2D under compression

We explore the elastic behavior of a wormlike chain under compression in terms of exact solutions for the associated probability densities. Strikingly, the probability density for the end-to-end distance projected along the applied force exhibits a bimodal shape in the vicinity of the critical Euler buckling force of an elastic rod, reminiscent of the smeared discontinuous phase transition of a finite system. These two modes reflect the almost stretched and the S-shaped configuration of a clamped polymer induced by the compression. Moreover, we find a bimodal shape of the probability density for the transverse fluctuations of the free end of a cantilevered polymer as fingerprint of its semiflexibility. In contrast to clamped polymers, free polymers display a circularly symmetric probability density and their distributions are identical for compression and stretching forces.

End-to-end distribution for a wormlike chain in arbitrary dimensions

We construct an efficient methodology for calculating wormlike chain statistics in arbitrary D dimensions over all chain rigidities, from fully rigid to completely flexible. The structure of our exact analytical solution for the end-to-end distribution function for a wormlike chain in arbitrary D dimensions in Fourier-Laplace space (i.e., Fourier-transformed end position and Laplace-transformed chain length) adopts the form of an infinite continued fraction, which is advantageous for its compact structure and stability for numerical implementation. We then proceed to present a step-by-step methodology for performing the Fourier-Laplace inversion in order to make full use of our results in general applications. Asymptotic methods for evaluating the Laplace inversion (power-law expansion and Rayleigh-Schrödinger perturbation theory) are employed in order to improve the accuracy of the numerical inversions of the end-to-end distribution function in real space. We adapt our results to the evaluation of the single-chain structure factor, rendering simple, closed-form expressions that facilitate comparison with scattering experiments. Using our techniques, the accuracy of the end-to-end distribution function is enhanced up to the limit of the machine precision. We demonstrate the utility of our methodology with realizations of the chain statistics, giving a general methodology that can be applied to a wide range of biophysical problems.