Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userStatistical Mechanics of Semiflexible Chains: A Meanfield Variational Approach
65
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
We describe a simple meanfield variational approach to study a number of properties of intrinsically stiff chains which are appropriate models for a large class of biopolymers. We present the calculation of the distribution of end-to-end distance and the elastic response of stiff chains under tension using this approach. In the former example we find that the simple expression almost quantitatively fits the results of computer simulation. For the case of the stiff chain under tension we recover analytically all the known limits. We obtain quantitative agreement with recent experiments on the stretching of DNA. The limitations of our approach are also discussed.
rate research
Read More
We propose a lattice model for RNA based on a self-interacting two-tolerant trail. Self-avoidance and elements of tertiary structure are taken into account. We investigate a simple version of the model in which the native state of RNA consists of just one hairpin. Using exact arguments and Monte Carlo simulations we determine the phase diagram for this case. We show that the denaturation transition is first order and can either occur directly or through an intermediate molten phase.
The basic notions of statistical mechanics (microstates, multiplicities) are quite simple, but understanding how the second law arises from these ideas requires working with cumbersomely large numbers. To avoid getting bogged down in mathematics, one can compute multiplicities numerically for a simple model system such as an Einstein solid -- a collection of identical quantum harmonic oscillators. A computer spreadsheet program or comparable software can compute the required combinatoric functions for systems containing a few hundred oscillators and units of energy. When two such systems can exchange energy, one immediately sees that some configurations are overwhelmingly more probable than others. Graphs of entropy vs. energy for the two systems can be used to motivate the theoretical definition of temperature, $T= (partial S/partial U)^{-1}$, thus bridging the gap between the classical and statistical approaches to entropy. Further spreadsheet exercises can be used to compute the heat capacity of an Einstein solid, study the Boltzmann distribution, and explore the properties of a two-state paramagnetic system.
Thermally fluctuating sheets and ribbons provide an intriguing forum in which to investigate strong violations of Hookes Law: large distance elastic parameters are in fact not constant, but instead depend on the macroscopic dimensions. Inspired by recent experiments on free-standing graphene cantilevers, we combine the statistical mechanics of thin elastic plates and large-scale numerical simulations to investigate the thermal renormalization of the bending rigidity of graphene ribbons clamped at one end. For ribbons of dimensions $Wtimes L$ (with $Lgeq W$), the macroscopic bending rigidity $kappa_R$ determined from cantilever deformations is independent of the width when $W<ell_textrm{th}$, where $ell_textrm{th}$ is a thermal length scale, as expected. When $W>ell_textrm{th}$, however, this thermally renormalized bending rigidity begins to systematically increase, in agreement with the scaling theory, although in our simulations we were not quite able to reach the system sizes necessary to determine the fully developed power law dependence on $W$. When the ribbon length $L > ell_p$, where $ell_p$ is the $W$-dependent thermally renormalized ribbon persistence length, we observe a scaling collapse and the beginnings of large scale random walk behavior.
The equations of hydrodynamics including mass, linear momentum, angular momentum, and energy are derived by coarse-graining the microscopic equations of motion for systems consisting of rotary dumbbells driven by internal torques.
Using a recently developed bead-spring model for semiflexible polymers that takes into account their natural extensibility, we report an efficient algorithm to simulate the dynamics for polymers like double-stranded DNA (dsDNA) in the absence of hydrodynamic interactions. The dsDNA is modelled with one bead-spring element per basepair, and the polymer dynamics is described by the Langevin equation. The key to efficiency is that we describe the equations of motion for the polymer in terms of the amplitudes of the polymers fluctuation modes, as opposed to the use of the physical positions of the beads. We show that, within an accuracy tolerance level of $5%$ of several key observables, the model allows for single Langevin time steps of $approx1.6$, 8, 16 and 16 ps for a dsDNA model-chain consisting of 64, 128, 256 and 512 basepairs (i.e., chains of 0.55, 1.11, 2.24 and 4.48 persistence lengths) respectively. Correspondingly, in one hour, a standard desktop computer can simulate 0.23, 0.56, 0.56 and 0.26 ms of these dsDNA chains respectively. We compare our results to those obtained from other methods, in particular, the (inextensible discretised) WLC model. Importantly, we demonstrate that at the same level of discretisation, i.e., when each discretisation element is one basepair long, our algorithm gains about 5-6 orders of magnitude in the size of time steps over the inextensible WLC model. Further, we show that our model can be mapped one-on-one to a discretised version of the extensible WLC model; implying that the speed-up we achieve in our model must hold equally well for the latter. We also demonstrate the use of the method by simulating efficiently the tumbling behaviour of a dsDNA segment in a shear flow.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
