BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Marco Baiesi

Marco Baiesi appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
cond-mat.stat-mechcond-mat.softphysics.geo-phBiomoleculesBiological Physicsphysics.comp-phphysics.data-anq-bio

Catalog footprint

What is connected

30works
8topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

30 published item(s)

preprint2022arXiv

Dynamic and Facilitated Binding of Topoisomerase Accelerates Topological Relaxation

How type 2 Topoisomerase (TopoII) proteins relax and simplify the topology of DNA molecules is one of the most intriguing open questions in biophysics. Most of the existing models neglect the dynamics of TopoII which is characteristics for proteins searching their targets via facilitated diffusion. Here, we show that dynamic binding of TopoII speeds up the topological relaxation of knotted substrates by enhancing the search of the knotted arc. Intriguingly, this in turn implies that the timescale of topological relaxation is virtually independent of the substrate length. We then discover that considering binding biases due to facilitated diffusion on looped substrates steers the sampling of the topological space closer to the boundaries between different topoisomers yielding an optimally fast topological relaxation. We discuss our findings in the context of topological simplification in vitro and in vivo.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Generalized equipartition for nonlinear multiplicative Langevin dynamics: application to laser-cooled atoms

The virial theorem, and the equipartition theorem in the case of quadratic degrees of freedom, are handy constraints on the statistics of equilibrium systems. Their violation is instrumental in determining how far from equilibrium a driven system might be. We extend the virial theorem to nonequilibrium conditions for Langevin dynamics with nonlinear friction and multiplicative noise. In particular, we generalize the equipartition theorem for confined laser-cooled atoms in the semi-classical regime. The resulting relation between the lowest moments of the atom position and velocity allows to measure in experiments how dissipative the cooling mechanism is. Moreover, its violation can reveal the departure from a strictly harmonic confinement or from the semi-classical regime.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Linear response of hydrodynamically-coupled particles under a nonequilibrium reservoir

A recent experiment driving colloids electromagnetically, by Bérut et al. [2014 Europhys. Lett. 107, 60004], is an ideal paradigm for illustrating a linear response theory for nonequilibrium overdamped systems including hydrodynamic interactions and, unusually, a reservoir itself out of equilibrium. Indeed, in this setup one finds a nonequilibrium environment in which the mobility and diffusivity of free particles are not simply proportional to each other. We derive both the response to a mechanical forcing and to temperature variations in terms of correlations between an observable and a path-weight action. The time-antisymmetric component of the latter turns out not to be simply proportional to the heat flowing into the environment. These results are visualized with simulations resembling conditions and protocols easily realizable in the experiment, thereby tracing a path for experimental verifications of the theory.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Nonequilibrium temperature response for stochastic overdamped systems

The thermal response of nonequilibrium systems requires the knowledge of concepts that go beyond entropy production. This is showed for systems obeying overdamped Langevin dynamics, either in steady states or going through a relaxation process. Namely, we derive the linear response to perturbations of the noise intensity, mapping it onto the quadratic response to a constant small force. The latter, displaying divergent terms, is explicitly regularized with a novel path-integral method. The nonequilibrium equivalents of heat capacity and thermal expansion coefficient are two applications of this approach, as we show with numerical examples.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Self-similar aftershock rates

In many important systems exhibiting crackling noise --- intermittent avalanche-like relaxation response with power-law and, thus, self-similar distributed event sizes --- the "laws" for the rate of activity after large events are not consistent with the overall self-similar behavior expected on theoretical grounds. This is in particular true for the case of seismicity and a satisfying solution to this paradox has remained outstanding. Here, we propose a generalized description of the aftershock rates which is both self-similar and consistent with all other known self-similar features. Comparing our theoretical predictions with high resolution earthquake data from Southern California we find excellent agreement, providing in particular clear evidence for a unified description of aftershocks and foreshocks. This may offer an improved way of time-dependent seismic hazard assessment and earthquake forecasting.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Thermal response of nonequilibrium RC-circuits

We analyze experimental data obtained from an electrical circuit having components at different temperatures, showing how to predict its response to temperature variations. This illustrates in detail how to utilize a recent linear response theory for nonequilibrium overdamped stochastic systems. To validate these results, we introduce a reweighting procedure that mimics the actual realization of the perturbation and allows extracting the susceptibility of the system from steady state data. This procedure is closely related to other fluctuation-response relations based on the knowledge of the steady state probability distribution. As an example, we show that the nonequilibrium heat capacity in general does not correspond to the correlation between the energy of the system and the heat flowing into it. Rather, also non-dissipative aspects are relevant in the nonequilbrium fluctuation response relations.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Torque-Induced Rotational Dynamics in Polymers: Torsional Blobs and Thinning

By using the blob theory and computer simulations, we investigate the properties of a linear polymer performing a stationary rotational motion around a long impenetrable rod. In particular, in the simulations the rotation is induced by a torque applied to the end of the polymer that is tethered to the rod. Three different regimes are found, in close analogy with the case of polymers pulled by a constant force at one end. For low torques the polymer rotates maintaining its equilibrium conformation. At intermediate torques the polymer assumes a trumpet shape, being composed by blobs of increasing size. At even larger torques the polymer is partially wrapped around the rod. We derive several scaling relations between various quantities as angular velocity, elongation and torque. The analytical predictions match the simulation data well. Interestingly, we find a "thinning" regime where the torque has a very weak (logarithmic) dependence on the angular velocity. We discuss the origin of this behavior, which has no counterpart in polymers pulled by an applied force.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Inflow rate, a time-symmetric observable obeying fluctuation relations

While entropy changes are the usual subject of fluctuation theorems, we seek fluctuation relations involving time-symmetric quantities, namely observables that do not change sign if the trajectories are observed backward in time. We find detailed and integral fluctuation relations for the (time integrated) difference between "entrance rate" and escape rate in mesoscopic jump systems. Such "inflow rate", which is even under time reversal, represents the discrete-state equivalent of the phase space contraction rate. Indeed, it becomes minus the divergence of forces in the continuum limit to overdamped diffusion. This establishes a formal connection between reversible deterministic systems and irreversible stochastic ones, confirming that fluctuation theorems are largely independent of the details of the underling dynamics.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Role of trapping and crowding as sources of negative differential mobility

Increasing the crowding in an environment does not necessarily trigger negative differential mobility of strongly pushed particles. Moreover, the choice of the model, in particular the kind of microscopic jump rates, may be very relevant in determining the mobility. We support these points via simple examples and we therefore address recent claims saying that crowding in an environment is likely to promote negative differential mobility. Trapping of tagged particles enhanced by increasing the force remains the mechanism determining a drift velocity not monotonous in the driving force.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Temperature response in nonequilibrium stochastic systems

The linear response to temperature changes is derived for systems with overdamped stochastic dynamics. Holding both in transient and steady state conditions, the results allow to compute nonequilibrium thermal susceptibilities from unperturbed correlation functions. These correlations contain a novel form of entropy flow due to temperature unbalances, next to the standard entropy flow of stochastic energetics and to complementary time-symmetric dynamical aspects. Our derivation hinges on a time rescaling, which is a key procedure for comparing apparently incommensurable path weights. An interesting notion of thermal time emerges from this approach.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Knotted globular ring polymers: how topology affects statistics and thermodynamics

The statistical mechanics of a long knotted collapsed polymer is determined by a free-energy with a knot-dependent subleading term, which is linked to the length of the shortest polymer that can hold such knot. The only other parameter depending on the knot kind is an amplitude such that relative probabilities of knots do not vary with the temperature $T$, in the limit of long chains. We arrive at this conclusion by simulating interacting self-avoiding rings at low $T$ on the cubic lattice, both with unrestricted topology and with setups where the globule is divided by a slip link in two loops (preserving their topology) which compete for the chain length, either in contact or separated by a wall as for translocation through a membrane pore. These findings suggest that in macromolecular environments there may be entropic forces with a purely topological origin, whence portions of polymers holding complex knots should tend to expand at the expense of significantly shrinking other topologically simpler portions.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Models of DNA denaturation dynamics: universal properties

We briefly review some of the models used to describe DNA denaturation dynamics, focusing on the value of the dynamical exponent $z$, which governs the scaling of the characteristic time $τ\sim L^z$ as a function of the sequence length $L$. The models contain different degrees of simplifications, in particular sometimes they do not include a description for helical entanglement: we discuss how this aspect influences the value of $z$, which ranges from $z=0$ to $z \approx 3.3$. Connections with experiments are also mentioned.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Rings in Random Environments: Sensing Disorder Through Topology

In this paper we study the role of topology in DNA gel electrophoresis experiments via molecular dynamics simulations. The gel is modelled as a 3D array of obstacles from which half edges are removed at random with probability p, thereby generating a disordered environment. Changes in the microscopic structure of the gel are captured by measuring the electrophoretic mobility of ring polymers moving through the medium, while their linear counterparts provide a control system as we show they are insensitive to these changes. We show that ring polymers provide a novel non-invasive way of exploiting topology to sense microscopic disorder. Finally, we compare the results from the simulations with an analytical model for the non-equilibrium differential mobility, and find a striking agreement between simulation and theory

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Rotational dynamics of entangled polymers

Some recent results on the rotational dynamics of polymers are reviewed and extended. We focus here on the relaxation of a polymer, either flexible or semiflexible, initially wrapped around a rigid rod. We also study the steady polymer rotation generated by a constant torque on the rod. The interplay of frictional and entropic forces leads to a complex dynamical behavior characterized by non-trivial universal exponents. The results are based on extensive simulations of polymers undergoing Rouse dynamics and on an analytical approach using force balance and scaling arguments. The analytical results are in general in good agreement with the simulations, showing how a simplified approach can correctly capture the complex dynamical behavior of rotating polymers.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Thermal response in driven diffusive systems

Evaluating the linear response of a driven system to a change in environment temperature(s) is essential for understanding thermal properties of nonequilibrium systems. The system is kept in weak contact with possibly different fast relaxing mechanical, chemical or thermal equilibrium reservoirs. Modifying one of the temperatures creates both entropy fluxes and changes in dynamical activity. That is not unlike mechanical response of nonequilibrium systems but the extra difficulty for perturbation theory via path-integration is that for a Langevin dynamics temperature also affects the noise amplitude and not only the drift part. Using a discrete-time mesh adapted to the numerical integration one avoids that ultraviolet problem and we arrive at a fluctuation expression for its thermal susceptibility. The algorithm appears stable under taking even finer resolution.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Unwinding Dynamics of a Helically Wrapped Polymer

We study the rotational dynamics of a flexible polymer initially wrapped around a rigid rod and unwinding from it. This dynamics is of interest in several problems in biology and constitutes a fundamental instance of polymer relaxation from a state of minimal entropy. We investigate the dynamics of several quantities such as the total and local winding angles and metric quantities. The results of simulations performed in two and three dimensions, with and without self-avoidance, are explained by a theory based on scaling arguments and on a balance between frictional and entropic forces. The early stage of the dynamics is particularly rich, being characterized by three coexisting phases.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Unwinding relaxation dynamics of polymers

The relaxation dynamics of a polymer wound around a fixed obstacle constitutes a fundamental instance of polymer with twist and torque and it is of relevance also for DNA denaturation dynamics. We investigate it by simulations and Langevin equation analysis. The latter predicts a relaxation time scaling as a power of the polymer length times a logarithmic correction related to the equilibrium fluctuations of the winding angle. The numerical data support this result and show that at short times the winding angle decreases as a power-law. This is also in agreement with the Langevin equation provided a winding-dependent friction is used, suggesting that such reduced description of the system captures the basic features of the problem.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

An update on nonequilibrium linear response

The unique fluctuation-dissipation theorem for equilibrium stands in contrast with the wide variety of nonequilibrium linear response formulae. Their most traditional approach is "analytic", which, in the absence of detailed balance, introduces the logarithm of the stationary probability density as observable. The theory of dynamical systems offers an alternative with a formula that continues to work when the stationary distribution is not smooth. We show that this method works equally well for stochastic dynamics, and we illustrate it with a numerical example for the perturbation of circadian cycles. A second "probabilistic" approach starts from dynamical ensembles and expands the probability weights on path space. This line suggests new physical questions, as we meet the frenetic contribution to linear response, and the relevance of the change in dynamical activity in the relaxation to a (new) nonequilibrium condition.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Fibril elongation mechanisms of HET-s prion-forming domain: Topological evidence for growth polarity

The prion-forming C-terminal domain of the fungal prion HET-s forms infectious amyloid fibrils at physiological pH. The conformational switch from the non-prion soluble form to the prion fibrillar form is believed to have a functional role, since HET-s in its prion form participates in a recognition process of different fungal strains. Based on the knowledge of the high-resolution structure of HET-s(218-289) (the prion forming-domain) in its fibrillar form, we here present a numerical simulation of the fibril growth process which emphasizes the role of the topological properties of the fibrillar structure. An accurate thermodynamic analysis of the way an intervening HET-s chain is recruited to the tip of the growing fibril suggests that elongation proceeds through a dock and lock mechanism. First, the chain docks onto the fibril by forming the longest $β$-strands. Then, the re-arrangement in the fibrillar form of all the rest of molecule takes place. Interestingly, we predict also that one side of the HET-s fibril is more suitable for substaining its growth with respect to the other. The resulting strong polarity of fibril growth is a consequence of the complex topology of HET-s fibrillar structure, since the central loop of the intervening chain plays a crucially different role in favouring or not the attachment of the C-terminus tail to the fibril, depending on the growth side.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Universal properties of knotted polymer rings

By performing Monte Carlo sampling of $N$-steps self-avoiding polygons embedded on different Bravais lattices we explore the robustness of universality in the entropic, metric and geometrical properties of knotted polymer rings. In particular, by simulating polygons with $N$ up to $10^5$ we furnish a sharp estimate of the asymptotic values of the knot probability ratios and show their independence on the lattice type. This universal feature was previously suggested although with different estimates of the asymptotic values. In addition we show that the scaling behavior of the mean squared radius of gyration of polygons depends on their knot type only through its correction to scaling. Finally, as a measure of the geometrical self-entanglement of the SAPs we consider the standard deviation of the writhe distribution and estimate its power-law behavior in the large $N$ limit. The estimates of the power exponent do depend neither on the lattice nor on the knot type, strongly supporting an extension of the universality property to some features of the geometrical entanglement.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The modified Sutherland--Einstein relation for diffusive nonequilibria

There remains a useful relation between diffusion and mobility for a Langevin particle in a periodic medium subject to nonconservative forces. The usual fluctuation-dissipation relation easily gets modified and the mobility matrix is no longer proportional to the diffusion matrix, with a correction term depending explicitly on the (nonequilibrium) forces. We discuss this correction by considering various simple examples and we visualize the various dependencies on the applied forcing and on the time by means of simulations. For example, in all cases the diffusion depends on the external forcing more strongly than does the mobility. We also give an explicit decomposition of the symmetrized mobility matrix as the difference between two positive matrices, one involving the diffusion matrix, the other force--force correlations.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Elastic Lattice Polymers

We study a model of "elastic" lattice polymer in which a fixed number of monomers $m$ is hosted by a self-avoiding walk with fluctuating length $l$. We show that the stored length density $ρ_m = 1 - <l>/m$ scales asymptotically for large $m$ as $ρ_m=ρ_\infty(1-θ/m + ...)$, where $θ$ is the polymer entropic exponent, so that $θ$ can be determined from the analysis of $ρ_m$. We perform simulations for elastic lattice polymer loops with various sizes and knots, in which we measure $ρ_m$. The resulting estimates support the hypothesis that the exponent $θ$ is determined only by the number of prime knots and not by their type. However, if knots are present, we observe strong corrections to scaling, which help to understand how an entropic competition between knots is affected by the finite length of the chain.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The entropic cost to tie a knot

We estimate by Monte Carlo simulations the configurational entropy of $N$-steps polygons in the cubic lattice with fixed knot type. By collecting a rich statistics of configurations with very large values of $N$ we are able to analyse the asymptotic behaviour of the partition function of the problem for different knot types. Our results confirm that, in the large $N$ limit, each prime knot is localized in a small region of the polygon, regardless of the possible presence of other knots. Each prime knot component may slide along the unknotted region contributing to the overall configurational entropy with a term proportional to $\ln N$. Furthermore, we discover that the mere existence of a knot requires a well defined entropic cost that scales exponentially with its minimal length. In the case of polygons with composite knots it turns out that the partition function can be simply factorized in terms that depend only on prime components with an additional combinatorial factor that takes into account the statistical property that by interchanging two identical prime knot components in the polygon the corresponding set of overall configuration remains unaltered. Finally, the above results allow to conjecture a sequence of inequalities for the connective constants of polygons whose topology varies within a given family of composite knot types.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Correlated earthquakes in a self-organized model

Motivated by the fact that empirical time series of earthquakes exhibit long-range correlations in space and time and the Gutenberg-Richter distribution of magnitudes, we propose a simple fault model that can account for these types of scale-invariance. It is an avalanching process that displays power-laws in the event sizes, in the epicenter distances as well as in the waiting-time distributions, and also aftershock rates obeying a generalized Omori law. We thus confirm that there is a relation between temporal and spatial clustering of the activity in this kind of models. The fluctuating boundaries of possible slipping areas show that the size of the largest possible earthquake is not always maximal, and the average correlation length is a fraction of the system size. This suggests that there is a concrete alternative to the extreme interpretation of self-organized criticality as a process in which every small event can cascade to an arbitrary large one: the new picture includes fluctuating domains of coherent stress field as part of the global self-organization. Moreover, this picture can be more easily compared with other scenarios discussing fluctuating correlations lengths in seismicity.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Nonequilibrium linear response for Markov dynamics, I: jump processes and overdamped diffusions

Systems out of equilibrium, in stationary as well as in nonstationary regimes, display a linear response to energy impulses simply expressed as the sum of two specific temporal correlation functions. There is a natural interpretation of these quantities. The first term corresponds to the correlation between observable and excess entropy flux yielding a relation with energy dissipation like in equilibrium. The second term comes with a new meaning: it is the correlation between the observable and the excess in dynamical activity or reactivity, playing an important role in dynamical fluctuation theory out-of-equilibrium. It appears as a generalized escape rate in the occupation statistics. The resulting response formula holds for all observables and allows direct numerical or experimental evaluation, for example in the discussion of effective temperatures, as it only involves the statistical averaging of explicit quantities, e.g. without needing an expression for the nonequilibrium distribution. The physical interpretation and the mathematical derivation are independent of many details of the dynamics, but in this first part they are restricted to Markov jump processes and overdamped diffusions.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Nonequilibrium Linear Response for Markov Dynamics, II: Inertial Dynamics

We continue our study of the linear response of a nonequilibrium system. This Part II concentrates on models of open and driven inertial dynamics but the structure and the interpretation of the result remain unchanged: the response can be expressed as a sum of two temporal correlations in the unperturbed system, one entropic, the other frenetic. The decomposition arises from the (anti)symmetry under time-reversal on the level of the nonequilibrium action. The response formula involves a statistical averaging over explicitly known observables but, in contrast with the equilibrium situation, they depend on the model dynamics in terms of an excess in dynamical activity. As an example, the Einstein relation between mobility and diffusion constant is modified by a correlation term between the position and the momentum of the particle.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Computation of current cumulants for small nonequilibrium systems

We analyze a systematic algorithm for the exact computation of the current cumulants in stochastic nonequilibrium systems, recently discussed in the framework of full counting statistics for mesoscopic systems. This method is based on identifying the current cumulants from a Rayleigh-Schrodinger perturbation expansion for the generating function. Here it is derived from a simple path-distribution identity and extended to the joint statistics of multiple currents. For a possible thermodynamical interpretation, we compare this approach to a generalized Onsager-Machlup formalism. We present calculations for a boundary driven Kawasaki dynamics on a one-dimensional chain, both for attractive and repulsive particle interactions.

0 citations0 reviewsNo code linkedOpen access
preprint2004arXiv

Complex networks of earthquakes and aftershocks

We invoke a metric to quantify the correlation between any two earthquakes. This provides a simple and straightforward alternative to using space-time windows to detect aftershock sequences and obviates the need to distinguish main shocks from aftershocks. Directed networks of earthquakes are constructed by placing a link, directed from the past to the future, between pairs of events that are strongly correlated. Each link has a weight giving the relative strength of correlation such that the sum over the incoming links to any node equals unity for aftershocks, or zero if the event had no correlated predecessors. A correlation threshold is set to drastically reduce the size of the data set without losing significant information. Events can be aftershocks of many previous events, and also generate many aftershocks. The probability distribution for the number of incoming and outgoing links are both scale free, and the networks are highly clustered. The Omori law holds for aftershock rates up to a decorrelation time that scales with the magnitude, $m$, of the initiating shock as $t_{\rm cutoff} \sim 10^{βm}$ with $β\simeq 3/4$. Another scaling law relates distances between earthquakes and their aftershocks to the magnitude of the initiating shock. Our results are inconsistent with the hypothesis of finite aftershock zones. We also find evidence that seismicity is dominantly triggered by small earthquakes. Our approach, using concepts from the modern theory of complex networks, together with a metric to estimate correlations, opens up new avenues of research, as well as new tools to understand seismicity.

0 citations0 reviewsNo code linkedOpen access
preprint2004arXiv

Scaling and precursor motifs in earthquake networks

A measure of the correlation between two earthquakes is used to link events to their aftershocks, generating a growing network structure. In this framework one can quantify whether an aftershock is close or far, from main shocks of all magnitudes. We find that simple network motifs involving links to far aftershocks appear frequently before the three biggest earthquakes of the last 16 years in Southern California. Hence, networks could be useful to detect symptoms typically preceding major events.

0 citations0 reviewsNo code linkedOpen access
preprint2002arXiv

A simple model of DNA denaturation and mutually avoiding walks statistics

Recently Garel, Monthus and Orland (Europhys. Lett. v 55, 132 (2001)) considered a model of DNA denaturation in which excluded volume effects within each strand are neglected, while mutual avoidance is included. Using an approximate scheme they found a first order denaturation. We show that a first order transition for this model follows from exact results for the statistics of two mutually avoiding random walks, whose reunion exponent is c > 2, both in two and three dimensions. Analytical estimates of c due to the interactions with other denaturated loops, as well as numerical calculations, indicate that the transition is even sharper than in models where excluded volume effects are fully incorporated. The probability distribution of distances between homologous base pairs decays as a power law at the transition.

0 citations0 reviewsNo code linkedOpen access