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

Henri Orland

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

ResearcherUnclaimed source record
cond-mat.softcond-mat.stat-mechphysics.chem-phcond-mat.mtrl-sciBiomoleculesBiological Physicscond-mat.mes-hallmath-phmath.MPcond-mat.dis-nnmath.PRphysics.comp-phQuantitative Methods

Catalog footprint

What is connected

30works
13topics
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)

preprint2023arXiv

Conductivity of Concentrated Electrolytes

The conductivity of ionic solutions is arguably their most important trait, being widely used in electrochemical, biochemical, and environmental applications. The Debye-Hückel-Onsager theory successfully predicts the conductivity at very low ionic concentrations of up to a few millimolars, but there is no well-established theory applicable at higher concentrations. We study the conductivity of ionic solutions using a stochastic density functional theory, paired with a modified Coulomb interaction that accounts for the hard-core repulsion between the ions. The modified potential suppresses unphysical, short-range electrostatic interactions, which are present in the Debye-Hückel-Onsager theory. Our results for the conductivity show very good agreement with experimental data up to 3 molars, without any fit parameters. We provide a compact expression for the conductivity, accompanied by a simple analytical approximation.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Conductance of concentrated electrolytes: multivalency and the Wien effect

The electric conductivity of ionic solutions is well understood at low ionic concentrations of up to a few millimolar but becomes difficult to unravel at higher concentrations that are still common in nature and technological applications. A model for the conductivity at high concentrations was recently put forth for monovalent electrolytes at low electric fields. The model relies on applying a stochastic density-functional theory and using a modified electrostatic pair-potential that suppresses unphysical, short-range electrostatic interactions. Here, we extend the theory to multivalent ions as well as to high electric fields where a deviation from Ohm's law known as the Wien effect occurs. Our results are in good agreement with experiments and recent simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Dielectric Constant of Ionic Solutions: Combined Effects of Correlations and Excluded Volume

The dielectric constant of ionic solutions is known to reduce with increasing ionic concentrations. However, the origin of this effect has not been thoroughly explored. In this paper we study two such possible sources: long-range Coulombic correlations and solvent excluded volume. Correlations originate from fluctuations of the electrostatic potential beyond the mean-field Poisson-Boltzmann theory, evaluated by employing a field-theoretical loop expansion of the free energy. The solvent excluded-volume, on the other hand, stems from the finite ion size, accounted for via a lattice-gas model. We show that both correlations and excluded volume are required in order to capture the important features of the dielectric behavior. For highly polar solvents, such as water, the dielectric constant is given by the product of the solvent volume fraction and a concentration-dependent susceptibility per volume fraction. The available solvent volume decreases as function of ionic strength due the increasing volume fraction of ions. A similar decrease occurs for the susceptibility due to correlations between the ions and solvent, reducing the dielectric response even further. Our predictions for the dielectric constant fit well with experiments for a wide range of concentrations for different salts in different temperatures, using a single fit parameter related to the ion size.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Sampling constrained stochastic trajectories using Brownian bridges

We present a new method to sample conditioned trajectories of a system evolving under Langevin dynamics, based on Brownian bridges. The trajectories are conditioned to end at a certain point (or in a certain region) in space. The bridge equations can be recast exactly in the form of a non linear stochastic integro-differential equation. This equation can be very well approximated when the trajectories are closely bundled together in space, i.e. at low temperature, or for transition paths. The approximate equation can be solved iteratively, using a fixed point method. We discuss how to choose the initial trajectories and show some examples of the performance of this method on some simple problems. The method allows to generate conditioned trajectories with a high accuracy.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Generating stochastic trajectories with global dynamical constraints

We propose a method to exactly generate Brownian paths $x_c(t)$ that are constrained to return to the origin at some future time $t_f$, with a given fixed area $A_f = \int_0^{t_f}dt\, x_c(t)$ under their trajectory. We derive an exact effective Langevin equation with an effective force that accounts for the constraint. In addition, we develop the corresponding approach for discrete-time random walks, with arbitrary jump distributions including Lévy flights, for which we obtain an effective jump distribution that encodes the constraint. Finally, we generalise our method to other types of dynamical constraints such as a fixed occupation time on the positive axis $T_f=\int_0^{t_f}dt\, Θ\left[x_c(t)\right]$ or a fixed generalised quadratic area $\mathcal{A}_f=\int_0^{t_f}dt \,x_c^2(t)$.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A polymer model for the quantitative reconstruction of 3d chromosome architecture from Hi-C and GAM data

It is widely believed that the folding of the chromosome in the nucleus has a major effect on genetic expression. For example co-regulated genes in several species have been shown to colocalize in space despite being far away on the DNA sequence. In this manuscript, we present a new method to model the three-dimensional structure of the chromosome in live cells, based on DNA-DNA interactions measured in high-throughput chromosome conformation capture experiments (Hi-C) and genome architecture mapping experiments (GAM). Our approach incorporates a polymer model, and directly uses the contact probabilities measured in Hi-C and GAM experiments rather than estimates of average distances between genomic loci. Specifically, we model the chromosome as a Gaussian polymer with harmonic interactions and extract the coupling coefficients best reproducing the experimental contact probabilities. In contrast to existing methods, we give an exact expression of the contact probabilities at thermodynamic equilibrium. The Gaussian effective model (GEM) reconstructed with our method reproduces experimental contacts with high accuracy. We also show how Brownian Dynamics simulations of our reconstructed GEM can be used to study chromatin organization, and possibly give some clue about its dynamics.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Multicanonical Monte Carlo Ensemble Growth Algorithm

We present a novel Ensemble Monte Carlo Growth method to sample the equilibrium thermodynamic properties of random chains. The method is based on the multicanonical technique of computing the density of states in the energy space. Such a quantity is temperature independent, and therefore microcanonical and canonical thermodynamic quantities, including the free energy, entropy, and thermal averages, can be obtained by re-weighting with a Boltzmann factor. The algorithm we present combines two approaches: the first is the Monte Carlo ensemble growth method, where a "population" of samples in the state space is considered, as opposed to traditional sampling by long random walks, or iterative single-chain growth. The second is the flat-histogram Monte Carlo, similar to the popular Wang-Landau sampling, or to multicanonical chain-growth sampling. We discuss the performance and relative simplicity of the proposed algorithm, and we apply it to known test cases.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Beyond Poisson-Boltzmann: Numerical sampling of charge density fluctuations

We present a method aimed at sampling charge density fluctuations in Coulomb systems. The derivation follows from a functional integral representation of the partition function in terms of charge density fluctuations. Starting from the mean-field solution given by the Poisson-Boltzmann equation, an original approach is proposed to numerically sample fluctuations around it, through the propagation of a Langevin like stochastic partial differential equation (SPDE). The diffusion tensor of the SPDE can be chosen so as to avoid the numerical complexity linked to long-range Coulomb interactions, effectively rendering the theory completely local. A finite-volume implementation of the SPDE is described, and the approach is illustrated with preliminary results on the study of a system made of two like-charge ions immersed in a bath of counter-ions.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Conditioned Langevin Dynamics enables efficient sampling of transition paths

We propose a novel stochastic method to generate Brownian paths conditioned to start at an initial point and end at a given final point during a fixed time $t_{f}$ under a given potential $U(x)$. These paths are sampled with a probability given by the overdamped Langevin dynamics. We show that these paths can be exactly generated by a local Stochastic Partial Differential Equation (SPDE). This equation cannot be solved in general. We present several approximations that are valid either in the low temperature regime or in the presence of barrier crossing. We show that this method warrants the generation of statistically independent transition paths. It is computationally very efficient. We illustrate the method on the two dimensional Mueller potential as well as on the Mexican hat potential.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Improved RNA pseudoknots prediction and classification using a new topological invariant

We propose a new topological characterization of RNA secondary structures with pseudoknots based on two topological invariants. Starting from the classic arc-representation of RNA secondary structures, we consider a model that couples both I) the topological genus of the graph and II) the number of crossing arcs of the corresponding primitive graph. We add a term proportional to these topological invariants to the standard free energy of the RNA molecule, thus obtaining a novel free energy parametrization which takes into account the abundance of topologies of RNA pseudoknots observed in RNA databases.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Ionic profiles close to dielectric discontinuities: Specific ion-surface interactions

We study, by incorporating short-range ion-surface interactions, ionic profiles of electrolyte solutions close to a non-charged interface between two dielectric media. In order to account for important correlation effects close to the interface, the ionic profiles are calculated beyond mean-field theory, using the loop expansion of the free energy. We show how it is possible to overcome the well-known deficiency of the regular loop expansion close to the dielectric jump, and treat the non-linear boundary conditions within the framework of field theory. The ionic profiles are obtained analytically to one-loop order in the free energy, and their dependence on different ion-surface interactions is investigated. The Gibbs adsorption isotherm, as well as the ionic profiles are used to calculate the surface tension, in agreement with the reverse Hofmeister series. Consequently, from the experimentally-measured surface tension, one can extract a single adhesivity parameter, which can be used within our model to quantitatively predict hard to measure ionic profiles.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Effective Langevin equations for constrained stochastic processes

We propose a novel stochastic method to exactly generate Brownian paths conditioned to start at an initial point and end at a given final point during a fixed time $t_{f}$. These paths are weighted with a probability given by the overdamped Langevin dynamics. We show how these paths can be exactly generated by a local stochastic differential equation. The method is illustrated on the generation of Brownian bridges, Brownian meanders, Brownian excursions and constrained Ornstein-Uehlenbeck processes. In addition, we show how to solve this equation in the case of a general force acting on the particle. As an example, we show how to generate constrained path joining the two minima of a double-well. Our method allows to generate statistically independent paths, and is computationally very efficient.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Phase behaviour of DNA in presence of DNA-binding proteins

To characterize the thermodynamical equilibrium of DNA chains interacting with a solution of non-specific binding proteins, a Flory-Huggins free energy model was implemented. We explored the dependence on DNA and protein concentrations of the DNA collapse. For physiologically relevant values of the DNA-protein affinity, this collapse gives rise to a biphasic regime with a dense and a dilute phase; the corresponding phase diagram was computed. Using an approach based on Hamiltonian paths, we show that the dense phase has either a molten globule or a crystalline structure, depending on the DNA bending rigidity, which is influenced by the ionic strength. These results are valid at the thermodynamical equilibrium and should therefore be consistent with many biological processes, whose characteristic timescales range typically from 1 ms to 10 s. Our model may thus be applied to biological phenomena that involve DNA-binding proteins, such as DNA condensation with crystalline order, which occurs in some bacteria to protect their chromosome from detrimental factors; or transcription initiation, which occurs in clusters called transcription factories that are reminiscent of the dense phase characterized in this study.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The unknotted strands of life: knots are very rare in RNA structures

The ongoing effort to detect and characterize physical entanglement in biopolymers has so far established that knots are present in many globular proteins and also abound in viral DNA packaged inside bacteriophages. RNA molecules, on the other hand, have not yet been systematically screened for the occurrence of physical knots. We have accordingly undertaken the systematic profiling of the ~6,000 RNA structures present in the protein data bank. The search identified no more than three deeply-knotted RNA molecules. These are ribosomal RNAs solved by cryo-em and consist of about 3,000 nucleotides. Compared to the case of proteins and viral DNA, the observed incidence of RNA knots is therefore practically negligible. This suggests that either evolutionary selection, or thermodynamic and kinetic folding mechanisms act towards minimizing the entanglement of RNA to an extent that is unparalleled by other types of biomolecules. The properties of the three observed RNA knotting patterns provide valuable clues for designing RNA sequences capable of self-tying in a twist-knot fold.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Coherent States Formulation of Polymer Field Theory

We introduce a stable and efficient complex Langevin (CL) scheme to enable the first numerical simulations of the coherent-states (CS) formulation of polymer field theory. In contrast with Edwards' well known auxiliary-field (AF) framework, the CS formulation does not contain an embedded non-linear, non-local functional of the auxiliary fields, and the action of the field theory has a fully explicit, finite-order and semi-local polynomial character. In the context of a polymer solution model, we demonstrate that the new CS-CL dynamical scheme for sampling fluctuations in the space of coherent states yields results in good agreement with now-standard AF simulations. The formalism is potentially applicable to a broad range of polymer architectures and may facilitate systematic generation of trial actions for use in coarse-graining and numerical renormalization-group studies.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Dipolar Poisson-Boltzmann Approach to Ionic Solutions: A Mean Field and Loop Expansion Analysis

We study the variation of the dielectric response of ionic aqueous solutions as function of their ionic strength. The effect of salt on the dielectric constant appears through the coupling between ions and dipolar water molecules. On a mean-field level, we account for any internal charge distribution of particles. The dipolar degrees of freedom are added to the ionic ones and result in a generalization of the Poisson-Boltzmann (PB) equation called the Dipolar PB (DPB). By looking at the DPB equation around a fixed point-like ion, a closed-form formula for the dielectric constant is obtained. We express the dielectric constant using the "hydration length" that characterizes the hydration shell of dipoles around ions, and thus the strength of the dielectric decrement. The DPB equation is then examined for three additional cases: mixture of solvents, polarizable medium and ions of finite size. Employing field-theoretical methods we expand the Gibbs free-energy to first order in a loop expansion and calculate self-consistently the dielectric constant. For pure water, the dipolar fluctuations represent an important correction to the mean-field value and good agreement with the water dielectric constant is obtained. For ionic solutions we predict analytically the dielectric decrement that depends on the ionic strength in a non-linear way. Our prediction fits rather well a large range of concentrations for different salts using only one fit parameter related to the size of ions and dipoles. A linear dependence of the dielectric constant on the salt concentration is observed at low salinity, and a noticeable deviation from linearity can be seen for ionic strength above 1\,M, in agreement with experiments.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Dynamics of Polymers: a Mean-Field Theory

We derive a general mean-field theory of inhomogeneous polymer dynamics; a theory whose form has been speculated and widely applied, but not heretofore derived. Our approach involves a functional integral representation of a Martin-Siggia-Rose type description of the exact many-chain dynamics. A saddle point approximation to the generating functional, involving conditions where the MSR action is stationary with respect to a collective density field $ρ$ and a conjugate MSR response field $ϕ$, produces the desired dynamical mean-field theory. Besides clarifying the proper structure of mean-field theory out of equilibrium, our results have implications for numerical studies of polymer dynamics involving hybrid particle-field simulation techniques such as the single-chain in mean-field method (SCMF).

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Block Copolymer at Nano-Patterned Surfaces

We present numerical calculations of lamellar phases of block copolymers at patterned surfaces. We model symmetric di-block copolymer films forming lamellar phases and the effect of geometrical and chemical surface patterning on the alignment and orientation of lamellar phases. The calculations are done within self-consistent field theory (SCFT), where the semi-implicit relaxation scheme is used to solve the diffusion equation. Two specific set-ups, motivated by recent experiments, are investigated. In the first, the film is placed on top of a surface imprinted with long chemical stripes. The stripes interact more favorably with one of the two blocks and induce a perpendicular orientation in a large range of system parameters. However, the system is found to be sensitive to its initial conditions, and sometimes gets trapped into a metastable mixed state composed of domains in parallel and perpendicular orientations. In a second set-up, we study the film structure and orientation when it is pressed against a hard grooved mold. The mold surface prefers one of the two components and this set-up is found to be superior for inducing a perfect perpendicular lamellar orientation for a wide range of system parameters.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Block Copolymer Films with Free Interfaces: Ordering by Nano-Patterned Substrates

We study block copolymers (BCP) on patterned substrates, where the top polymer film surface is not constrained but is a free interface that adapts its shape self-consistently. In particular, we investigate the combined effect of the free interface undulations with the wetting of the BCP film as induced by nano-patterned substrates. For a finite volume of BCP material, we find equilibrium droplets composed of coexisting perpendicular and parallel lamellar domains. The self-assembly of BCP on topographic patterned substrates was also investigated and it was found that the free interface induces mixed morphologies of parallel and perpendicular domains coupled with a non-flat free interface. In both cases, the free interface relaxes the strong constraints that would otherwise be imposed by a fixed top boundary (which is commonly used in simulations), and affects strongly the BCP ordering. Our study has some interesting consequences for experimental setups of graphoepitaxy and nanoimprint lithography.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Tailoring Nanostructures Using Copolymer Nanoimprint Lithography

Finding affordable ways of generating high-density ordered nanostructures that can be transferred to a substrate is a major challenge for industrial applications like memories or optical devices with high resolution features. In this work, we report on a novel technique to direct self-assembled structures of block copolymers by NanoImprint Lithography. Surface energy of a reusable mold and nanorheology are used to organize the copolymers in defect-free structures over tens of micrometers in size. Versatile and controlled in-plane orientations of about 25 nm half-period lamellar nanostructures are achieved and, in particular, include applications to circular tracks of magnetic reading heads.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The Dielectric Constant of Ionic Solutions: A Field-Theory Approach

We study the variation of the dielectric response of a dielectric liquid (e.g. water) when a salt is added to the solution. Employing field-theoretical methods we expand the Gibbs free-energy to first order in a loop expansion and calculate self-consistently the dielectric constant. We predict analytically the dielectric decrement which depends on the ionic strength in a complex way. Furthermore, a qualitative description of the hydration shell is found and is characterized by a single length scale. Our prediction fits rather well a large range of concentrations for different salts using only one fit parameter related to the size of ions and dipoles.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

A field-theoretic approach to nonequilibrium work identities

We study nonequilibrium work relations for a space-dependent field with stochastic dynamics (Model A). Jarzynski's equality is obtained through symmetries of the dynamical action in the path integral representation. We derive a set of exact identities that generalize the fluctuation-dissipation relations to non-stationary and far-from-equilibrium situations. These identities are prone to experimental verification. Furthermore, we show that a well-studied invariance of the Langevin equation under supersymmetry, which is known to be broken when the external potential is time-dependent, can be partially restored by adding to the action a term which is precisely Jarzynski's work. The work identities can then be retrieved as consequences of the associated Ward-Takahashi identities.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Generating Transition Paths by Langevin Bridges

We propose a novel stochastic method to generate paths conditioned to start in an initial state and end in a given final state during a certain time $t_{f}$. These paths are weighted with a probability given by the overdamped Langevin dynamics. We show that these paths can be exactly generated by a non-local stochastic differential equation. In the limit of short times, we show that this complicated non-solvable equation can be simplified into an approximate stochastic differential equation. For longer times, the paths generated by this approximate equation can be reweighted to generate the correct statistics. In all cases, the paths generated by this equation are statistically independent and provide a representative sample of transition paths. In case the reaction takes place in a solvent (e.g. protein folding in water), the explicit solvent can be treated. The method is illustrated on the one-dimensional quartic oscillator.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Organization of Block Copolymers using NanoImprint Lithography: Comparison of Theory and Experiments

We present NanoImprint lithography experiments and modeling of thin films of block copolymers (BCP). The NanoImprint lithography is used to align perpendicularly lamellar phases, over distances much larger than the natural lamellar periodicity. The modeling relies on self-consistent field calculations done in two- and three-dimensions. We get a good agreement with the NanoImprint lithography setups. We find that, at thermodynamical equilibrium, the ordered BCP lamellae are much better aligned than when the films are deposited on uniform planar surfaces.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Secondary structure formation of homopolymeric single-stranded nucleic acids including force and loop entropy: implications for DNA hybridization

Loops are essential secondary structure elements in folded DNA and RNA molecules and proliferate close to the melting transition. Using a theory for nucleic acid secondary structures that accounts for the logarithmic entropy c ln m for a loop of length m, we study homopolymeric single-stranded nucleic acid chains under external force and varying temperature. In the thermodynamic limit of a long strand, the chain displays a phase transition between a low temperature / low force compact (folded) structure and a high temperature / high force molten (unfolded) structure. The influence of c on phase diagrams, critical exponents, melting, and force extension curves is derived analytically. For vanishing pulling force, only for the limited range of loop exponents 2 < c < 2.479 a melting transition is possible; for c <= 2 the chain is always in the folded phase and for 2.479 < c always in the unfolded phase. A force induced melting transition with singular behavior is possible for all loop exponents c < 2.479 and can be observed experimentally by single molecule force spectroscopy. These findings have implications for the hybridization or denaturation of double stranded nucleic acids. The Poland-Scheraga model for nucleic acid duplex melting does not allow base pairing between nucleotides on the same strand in denatured regions of the double strand. If the sequence allows these intra-strand base pairs, we show that for a realistic loop exponent c ~ 2.1 pronounced secondary structures appear inside the single strands. This leads to a lower melting temperature of the duplex than predicted by the Poland-Scheraga model. Further, these secondary structures renormalize the effective loop exponent c^, which characterizes the weight of a denatured region of the double strand, and thus affect universal aspects of the duplex melting transition.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Accelerated Stochastic Sampling of Discrete Statistical Systems

We propose a method to reduce the relaxation time towards equilibrium in stochastic sampling of complex energy landscapes in statistical systems with discrete degrees of freedom by generalizing the platform previously developed for continuous systems. The method starts from a master equation, in contrast to the Fokker-Planck equation for the continuous case. The master equation is transformed into an imaginary-time Schrödinger equation. The Hamiltonian of the Schrödinger equation is modified by adding a projector to its known ground state. We show how this transformation decreases the relaxation time and propose a way to use it to accelerate simulated annealing for optimization problems. We implement our method in a simplified kinetic Monte Carlo scheme and show an acceleration by an order of magnitude in simulated annealing of the symmetric traveling salesman problem. Comparisons of simulated annealing are made with the exchange Monte Carlo algorithm for the three-dimensional Ising spin glass. Our implementation can be seen as a step toward accelerating the stochastic sampling of generic systems with complex landscapes and long equilibration times.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

TT2NE: A novel algorithm to predict RNA secondary structures with pseudoknots

We present TT2NE, a new algorithm to predict RNA secondary structures with pseudoknots. The method is based on a classification of RNA structures according to their topological genus. TT2NE guarantees to find the minimum free energy structure irrespectively of pseudoknot topology. This unique proficiency is obtained at the expense of the maximum length of sequence that can be treated but comparison with state-of-the-art algorithms shows that TT2NE is a very powerful tool within its limits. Analysis of TT2NE's wrong predictions sheds light on the need to study how sterical constraints limit the range of pseudoknotted structures that can be formed from a given sequence. An implementation of TT2NE on a public server can be found at http://ipht.cea.fr/rna/tt2ne.php.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Accelerated Sampling of Boltzmann distributions

The sampling of Boltzmann distributions by stochastic Markov processes, can be strongly limited by the crossing time of high (free) energy barriers. As a result, the system may stay trapped in metastable states, and the relaxation time to the equilibrium Boltzmann distribution may be very large compared to the available computational time. In this paper, we show how, by a simple modification of the Hamiltonian, one can dramatically decrease the relaxation time of the system, while retaining the same equilibrium distribution. The method is illustrated on the case of the one-dimensional double-well potential.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

Dipolar Poisson-Boltzmann Equation: Ions and Dipoles Close to Charged Surfaces

We present an extension to the Poisson-Boltzmann model where the dipolar features of solvent molecules are taken explicitly into account. The formulation is derived at mean-field level and can be extended to any order in a systematic expansion. It is applied to a two-plate system with oppositely charged surfaces. The ion distribution and profiles in the dipolar order parameter are calculated and can result in a large correction to the inter-plate pressure.

0 citations0 reviewsNo code linkedOpen access
preprint1999arXiv

Beyond Poisson-Boltzmann: Fluctuations and Correlations

We formulate the non-linear field theory for a fluctuating counter-ion distribution in the presence of a fixed, arbitrary charge distribution. The Poisson-Boltzmann equation is obtained as the saddle-point, and the effects of fluctuations and correlations are included by a loop-wise expansion around this saddle point. We show that the Poisson equation is obeyed at each order in the loop expansion and explicitly give the expansion of the Gibbs potential up to two loops. We then apply our formalism to the case of an impenetrable, charged wall, and obtain the fluctuation corrections to the electrostatic potential and counter-ion density to one-loop order without further approximations. The relative importance of fluctuation corrections is controlled by a single parameter, which is proportional to the cube of the counter-ion valency and to the surface charge density. We also calculate effective interactions between charged particles, which reflect counter-ion correlation effects.

0 citations0 reviewsNo code linkedOpen access