Source author record

Sahin Buyukdagli

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

34works
4topics
4close collaborators

Actions

Connect this record

Log in to claim

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

34 published item(s)

preprint2022arXiv

Explicit Solvent Effects on Macromolecular Interactions From a Solvent-Augmented Contact Value Theorem

The Derjaguin-Landau-Verywey-Overbeek (DLVO) theory has been a remarkably accurate framework for the characterization of macromolecular stability in water solvent. In view of its solvent-implicit nature neglecting the electrostatics of water molecules with non-negligible charge structure and concentration, the precision of the DLVO formalism is somewhat puzzling. In order to shed light on this issue, we derive from our earlier explicit solvent formalism [S. Buyukdagli el al., Phys. Rev. E, 2013, 87, 063201] a solvent-augmented contact value theorem and assess the contribution of solvent molecules to the interaction of charged membranes. We find that in the case of hydrophobic membranes with fixed charges embedded in the membrane surface, the nearly exact cancellation of various explicit solvent effects of substantially large magnitude but opposite sign keeps the intermembrane pressure significantly close to the double layer force of the DLVO theory. Then, in the case of hydrophilic surface charge groups within the aqueous region, due to the spatial separation of the membrane substrate from the location of the fixed charges where the non-local dielectric response of the structured solvent is sharply localized, the interfacial field energy and the contact charge densities remain unaffected by the explicit solvent. As a result, the hydration of the lipid head groups suppresses the signature of the solvent molecules from the membrane interaction force.

preprint2022arXiv

Explicit Solvent Theory of Salt-Induced Dielectric Decrement

We introduce a field-theoretic electrolyte model composed of structured solvent molecules and salt ions coupled by electrostatic and hard-core interactions. Within this explicit solvent theory, we characterize the salt-driven dielectric decrement beyond weak-coupling (WC) electrostatics. The WC approximation of prior formalisms is relaxed by treating the salt charges via a virial expansion. This virial approach enables the explicit inclusion of the many-body salt-solvent interactions, and directly leads to the experimentally observed linear decay of the electrolyte permittivity with added dilute salt. The permittivity formula emerging from our approach indicates that the reduction of the solvent permittivity is induced by the salt screening of the polarization charges suppressing the dielectric response of the solvent. By comparison with experiments, we also show that the salt-dressed permittivity formula can equally reproduce the attenuation of the electrolyte permittivity with rising temperature, the thermal decay of the dielectric decrement, and its intensification with the salt valency. The consistent qualitative agreement of our theory with this wide range of experimental trends points out the electrostatic ion-solvent correlations as the primary mechanism behind the salt-induced dielectric decrement.

preprint2021arXiv

Dielectric manipulation of polymer translocation dynamics in engineered membrane nanopores

The alteration of the dielectric membrane properties by membrane engineering techniques such as carbon nanotube (CNT) coating opens the way to novel molecular transport strategies for biosensing purposes. In this article, we predict a macromolecular transport mechanism enabling the dielectric manipulation of the polymer translocation dynamics in dielectric membrane pores confining mixed electrolytes. In the giant permittivity regime of these engineered membranes governed by attractive polarization forces, multivalent ions adsorbed by the membrane nanopore trigger a monovalent ion separation and set an electroosmotic counterion flow. The drag force exerted by this flow is sufficiently strong to suppress and invert the electrophoretic velocity of anionic polymers, and also to generate the mobility of neutral polymers whose speed and direction can be solely adjusted by the charge and concentration of the added multivalent ions. These features identify the dielectrically generated transport mechanism as an efficient mean to drive overall neutral or weakly charged analytes that cannot be controlled by an external voltage. We also reveal that in anionic polymer translocation, multivalent cation addition into the monovalent salt solution amplifies the electric current signal by several factors. The signal amplification is caused by the electrostatic many-body interactions replacing the monovalent polymer counterions by the multivalent cations of higher electric mobility. The strength of this elektrokinetic charge discrimination points out the potential of multivalent ions as current amplifiers capable of providing boosted resolution in nanopore-based biosensing techniques.

preprint2020arXiv

Interactions between Zwitterionic Membranes in Complex Electrolytes

We investigate the electrostatic interactions of zwitterionic membranes immersed in mixed electrolytes composed of mono- and multivalent ions. We show that the presence of monovalent salt is a necessary condition for the existence of a finite electrostatic force on the membrane. As a result, the mean-field membrane pressure originating from the surface dipoles exhibits a non-uniform salt dependence, characterized by an enhancement for dilute salt conditions and a decrease at intermediate salt concentrations. Upon addition of multivalent cations to the submolar salt solution, the separate interactions of these cations with the opposite charges of the surface dipoles makes the intermembrane pressure more repulsive at low membrane separation distances and strongly attractive at intermediate distances, resulting in a discontinuous like-charge binding transition followed by the membrane binding transition. By extending our formalism to account for correlation corrections associated with large salt concentrations, we show that membranes of high surface dipole density immersed in molar salt solutions may undergo a membrane binding transition even without the multivalent cations. Hence, the tuning of the surface polarization forces by membrane engineering can be an efficient way to adjust the equilibrium configuration of dipolar membranes in concentrated salt solutions.

preprint2020arXiv

Pulling a DNA molecule through a nanopore embedded in an anionic membrane: tension propagation coupled to electrostatics

We consider the influence of electrostatic forces on driven translocation dynamics of a flexible polyelectrolyte being pulled through a nanopore by an external force on the head monomer. To this end, we augment the iso-flux tension propagation (IFTP) theory with electrostatics for a negatively charged biopolymer pulled through a nanopore embedded in a similarly charged anionic membrane. We show that for the realistic case such as a single-stranded DNA, the translocation dynamics at low salt where screening is weak and at finite negative membrane charge is unexpectedly accelerated despite the large repulsive electrostatic interactions between the polymer coil on the {\it cis} side and the charged membrane. This is due to the rapid release of the electrostatic potential energy of the coil during translocation, leading to an effectively attractive force that assists end-driven translocation. The speedup results in non-monotonic polymer length and membrane charge dependence of the exponent $α$ characterizing the translocation time $τ\propto N_0^α$ of the polymer with length $N_0$. In the regime of long polymers $N_0\gtrsim500$, the translocation exponent exceeds its upper limit $α=2$ previously observed for the same system without electrostatic interactions.

preprint2019arXiv

Schwinger-Dyson equations for composite electrolytes governed by mixed electrostatic couplings strengths

The electrostatic Schwinger-Dyson equations are derived and solved for an electrolyte mixture composed of mono- and multivalent ions confined to a negatively charged nanoslit. The closure of these equations is based on an asymmetric treatment of the ionic species with respect to their electrostatic coupling strength; the weakly coupled monovalent ions are treated within a gaussian approximation while the multivalent counterions of high coupling strength are incorporated with a strong-coupling approach. The resulting self-consistent formalism includes explicitly the interactions of the multivalent counterions with the monovalent salt. In highly charged membranes characterized by a pronounced multivalent counterion adsorption, these interactions take over the salt-membrane charge coupling. As a result, the increment of the negative membrane charge brings further salt anions into the pore and excludes salt cations from the pore into the reservoir. The corresponding like-charge attraction and opposite-charge repulsion effect is amplified by the pore confinement but suppressed by salt addition into the reservoir. The effect is particularly pronounced in high dielectric membranes where the attractive polarization forces lead to a dense multivalent cation layer at the membrane walls. These cation layers act as an effective positive surface charge, resulting in a total monovalent cation exclusion and a strong anion excess even in the case of neutral membrane walls.

preprint2016arXiv

Beyond Poisson-Boltzmann: fluctuations and fluid structure in a self-consistent theory

Poisson-Boltzmann (PB) theory is the classic approach to soft matter electrostatics which has been applied to numerous problems of physical chemistry and biophysics. Its essential limitations are the neglect of correlation effects and of fluid structure. Recently, several theoretical insights have allowed the formulation of approaches that go beyond PB theory in a systematic way. In this topical review we provide an update on the developments achieved in self-consistent formulations of correlation-corrected Poisson-Boltzmann theory. We introduce the corresponding system of coupled nonlinear equations for both continuum electrostatics with a uniform dielectric constant and a structured solvent, a dipolar Coulomb fluid, including nonlocal effects. While the approach is only approximate and also limited to corrections in the so-called weak fluctuation regime, it allows to include physically relevant effects, as we show for a range of applications of these equations.

preprint2016arXiv

Correlation-induced DNA adsorption on like-charged membranes

The adsorption of DNA or other polyelectrolyte molecules on charged membranes is a recurrent motif in soft matter and bionanotechnological systems. Two typical situations encountered are the deposition of single DNA chains onto substrates for further analysis, e.g., by force microscopy, or the pulling of polyelectrolytes into membrane nanopores, as in sequencing applications. In this paper, we present a theoretical analysis of such scenarios based on the self-consistent field theory approach, which allows us to address the important effect of charge correlations. We calculate the grand potential of a stiff polyelectrolyte immersed in an electrolyte in contact with a negatively charged dielectric membrane. For the sake of conciseness, we neglect conformational polymer fluctuations and model the molecule as a rigid charged line. At strongly charged membranes, the adsorbed counterions enhance the screening ability of the interfacial region. In the presence of highly charged polymers such as double-stranded DNA molecules close to the membrane, this enhanced interfacial screening dominates the mean-field level DNA-membrane repulsion and results in the adsorption of the DNA molecule to the surface. This picture provides a simple explanation for the recently observed DNA binding onto similarly charged substrates [G. L.-Caballero et al., Soft Matter 10, 2805 (2014)] and points out charge correlations as a non-negligible ingredient of polymer-surface interactions.

preprint2016arXiv

Electrostatics of polymer translocation events in electrolyte solutions

We develop an analytical theory that accounts for the image and surface charge interactions between a charged dielectric membrane and a DNA molecule translocating through the membrane. Translocation events through neutral carbon-based membranes are driven by a competition between the repulsive DNA-image-charge interactions and the attractive coupling between the DNA segments on the trans and the cis sides of the membrane. The latter effect is induced by the reduction of the coupling by the dielectric membrane. In strong salt solutions where the repulsive image-charge effects dominate the attractive trans-cis coupling, the DNA molecule encounters a translocation barrier of ~10 kBT. In dilute electrolytes, the trans-cis coupling takes over image-charge forces and the membrane becomes a metastable attraction point that can trap translocating polymers over long time intervals. This mechanism can be used in translocation experiments in order to control the DNA motion by tuning the salt concentration of the solution.

preprint2015arXiv

Dielectric anisotropy in polar solvents under external fields

We investigate dielectric saturation and increment in polar liquids under external fields. We couple a previously introduced dipolar solvent model to a uniform electric field and derive the electrostatic kernel of interacting dipoles. This procedure allows an unambiguous definition of the liquid dielectric permittivity embodying non-linear dielectric response and correlation effects.We find that the presence of the external field results in a dielectric anisotropy characterized by a two-component dielectric permittivity tensor. The increase of the electric field amplifies the permittivity component parallel to the field direction, i.e. dielectric increment is observed along the field. However, the perpendicular component is lowered below the physiological permittivity, indicating dielectric saturation perpendicular to the field. By comparison with Molecular Dynamics simulations from the literature, we show that the mean-field level dielectric response theory underestimates dielectric saturation. The inclusion of dipolar correlations at the weak-coupling level intensify the mean-field level dielectric saturation and improves the agreement with simulation data at weak electric fields. The correlation-corrected theory predicts as well the presence of a metastable configuration corresponding to the antiparallel alignment of dipoles with the field. This prediction can be verified by solvent-explicit simulations where solvent molecules are expected to be trapped transiently in this metastable state.

preprint2015arXiv

Electrostatic energy barriers from dielectric membranes upon approach of translocating DNA molecules

We probe the electrostatic cost associated with the approach phase of DNA translocation events. Within an analytical theory at the Debye-Huckel level, we calculate the electrostatic free energy of a rigid DNA molecule interacting with a dielectric membrane. For carbon or silicon based low permittivity neutral membranes, the DNA molecule experiences a repulsive energy barrier between 10 kBT and 100 kBT. In the case of engineered membranes with high dielectric permittivity, the membrane surface attracts the DNA with an energy of the same magnitude. Both the repulsive and attractive interactions result from image-charge effects and their magnitude survive even for the thinnest graphene-based membranes of size d~6 A. For weakly charged membranes, the electrostatic free energy is always attractive at large separation distances but switches to repulsive close to the membrane surface. We also characterise the polymer length dependence of the interaction energy. For specific values of the membrane charge density, low permittivity membranes repel short polymers but attract long polymers. Our results can be used to control the strong electrostatic free energy of DNA-membrane interactions prior to translocation events by chemical engineering of the relevant system parameters.

preprint2015arXiv

Electrostatic interactions in charged nanoslits within an explicit solvent theory

Within a dipolar Poisson-Boltzmann theory including electrostatic correlations, we consider the effect of explicit solvent structure on solvent and ion partition confined to charged nanopores. We develop a relaxation scheme for the solution of this highly non-linear integro-differential equation for the electrostatic potential. The scheme is an extension of the approach previously introduced for simple planes (S. Buyukdagli and Ralf Blossey, J. Chem. Phys. 140, 234903 (2014)) to nanoslit geometry. We show that the reduced dielectric response of solvent molecules at the membrane walls gives rise to an electric field significantly stronger than the field of the classical Poisson-Boltzmann equation. This peculiarity associated with non-local electrostatic interactions results in turn in an interfacial counterion adsorption layer absent in continuum theories. The observation of this enhanced counterion affinity in the very close vicinity of the interface may have important impacts on nanofludic transport through charged nanopores. Our results indicate the quantitative inaccuracy of solvent implicit nanofiltration theories in predicting the ionic selectivity of membrane nanopores.

preprint2015arXiv

Ionic current inversion in pressure-driven polymer translocation through nanopores

We predict streaming current inversion with multivalent counterions in hydrodynamically driven polymer translocation events from a correlation-corrected charge transport theory including charge fluctuations around mean-field electrostatics. In the presence of multivalent counterions, electrostatic many-body effects result in the reversal of the DNA charge. The attraction of anions to the charge-inverted DNA molecule reverses the sign of the ionic current through the pore. Our theory allows for a comprehensive understanding of the complex features of the resulting streaming currents. The underlying mechanism is an efficient way to detect DNA charge reversal in pressure-driven translocation experiments with multivalent cations.

preprint2014arXiv

Controlling polymer translocation and ion transport via charge correlations

We develop a correlation-corrected transport theory in order to predict ionic and polymer transport properties of membrane nanopores in physical conditions where mean-field electrostatics breaks down. The experimentally observed low KCl conductivity of open alpha-Hemolysin pores is quantitatively explained by the presence of surface polarization effects. Upon the penetration of a DNA molecule into the pore, these polarization forces combined with the electroneutrality of DNA sets a lower boundary for the ionic current, explaining the weak salt dependence of blocked pore conductivities at dilute ion concentrations. The addition of multivalent counterions into the solution results in the reversal of the polymer charge and the direction of the electroosmotic flow. With trivalent spermidine or quadrivalent spermine molecules, the charge inversion is strong enough to stop the translocation of the polymer and to reverse its motion. This mechanism can be used efficiently in translocation experiments in order to improve the accuracy of DNA sequencing by minimizing the translocation velocity of the polymer.

preprint2014arXiv

Dipolar correlations in structured solvents under nanoconfinement

We study electrostatic correlations in structured solvents confined to nanoscale systems. We derive variational equations of Netz-Orland type for a model liquid composed of finite size dipoles. These equations are solved for both dilute solvents and solvents at physiological concentrations in a slit nanopore geometry. Correlation effects are of major importance for the dielectric reduction and anisotropy of the solvent resulting from dipole image interactions and also lead to a reduction of van der Waals attractions between low dielectric bodies. Finally, by comparison with other recently developed self-consistent theories and experiments, we scrutinize the effect of solvent-membrane interactions on the differential capacitance of the charged liquid in contact with low dielectric substrates. The interfacial solvent depletion driven by solvent-image interactions plays the major role in the observed low values of the experimental capacitance data, while non-locality associated with the extended charge structure of solvent molecules only brings a minor contribution.

preprint2014arXiv

Electrostatic correlations on the ionic selectivity of cylindrical membrane nanopores

We characterize the role of electrostatic fluctuations on the charge selectivity of cylindrical nanopores confining electrolyte mixtures. To this end, we develop an extended one-loop theory that can account for correlation effects induced by the surface charge, nanoconfinement of the electrolyte, and interfacial polarization charges associated with the low permittivity membrane. We validate the quantitative accuracy of the theory by comparisons with previously obtained Monte-Carlo simulation data from the literature, and scrutinize in detail the underlying forces driving the ionic selectivity of the nanopore. In the biologically relevant case of electrolytes with divalent cations such as CaCl2 in negatively charged nanopores, electrostatic correlations associated with the dense counterion layer in the channel result in an increase of the pore coion density with the surface charge. This peculiarity analogous to the charge inversion phenomenon remains intact for dielectrically inhomogeneous pores, which indicates that the effect should be observable in nanofiltration membranes or DNA-blocked nanopores characterized by a low membrane permittivity. Our results show that a quantitatively accurate consideration of correlation effects is necessary to determine the ionic selectivity of nanopores in the presence of electrolytes with multivalent counterions.

preprint2014arXiv

Nonlocal and nonlinear electrostatics of a dipolar Coulomb fluid

We study a model Coulomb fluid consisting of dipolar solvent molecules of finite extent generalizing the point-like Dipolar Poisson-Boltzmann model (DPB) previously introduced by Coalson and Duncan (J. Phys. Chem 100, 2612 (1996)) and Abrashkin et al. (Phys. Rev. Lett. 99, 077801 (2007)). We formulate a nonlocal Poisson-Boltzmann equation (NLPB) and study both linear and nonlinear dielectric response in this model for the case of a single plane geometry. Our results shed light on the relevance of nonlocal vs nonlinear effects in continuum models of material electrostatics.

preprint2013arXiv

Alteration of gas phase ion polarizabilities upon hydration in high dielectric liquids

We investigate the modification of gas phase ion polarizabilities upon solvation in polar solvents and ionic liquids. To this aim, we develop a classical electrostatic theory of charged liquids composed of solvent molecules modeled as finite size dipoles, and embedding polarizable ions that consist of Drude oscillators. In qualitative agreement with ab-initio calculations of polar solvents and ionic liquids, the hydration energy of a polarizable ion in both type of dielectric liquid is shown to favor the expansion of its electronic cloud. Namely, the ion carrying no dipole moment in the gas phase acquires a dipole moment in the liquid environment, but its electron cloud also reaches an enhanced rigidity. We find that the overall effect is an increase of the gas phase polarizability upon hydration. In the specific case of ionic liquids, it is shown that this hydration process is driven by a collective solvation mechanism where the dipole moment of a polarizable ion induced by its interaction with surrounding ions self-consistently adds to the polarization of the liquid, thereby amplifying the dielectric permittivity of the medium in a substantial way. We propose this self-consistent hydration as the underlying mechanism behind the high dielectric permittivities of ionic liquids composed of small charges with negligible gas phase dipole moment. Hydration being a correlation effect, the emerging picture indicates that electrostatic correlations cannot be neglected in polarizable liquids.

preprint2013arXiv

Microscopic formulation of non-local electrostatics in polar liquids embedding polarizable ions

Non-local electrostatic interactions associated with the finite solvent size and ion polarizability are investigated within the mean-field linear response theory. To this end, we introduce a field theoretic model of a polar liquid composed of linear multipole solvent molecules and embedding polarizable ions modeled as Drude oscillators. Unlike previous dipolar Poisson-Boltzmann formulations treating the solvent molecules as point dipoles, our model is able to qualitatively reproduce the non-local dielectric response behavior of polar liquids observed in Molecular Dynamics simulations and Atomic Force Microscope experiments for water solvent at charged interfaces. The present theory explains the formation of the associated interfacial hydration layers in terms of a cooperative dipolar response mechanism driven by the reaction of the solvent molecules to their own polarization field. We also incorporate into the theory the relative multipole/dipole moments of water molecules obtained from quantum mechanical calculations, and show that the multipolar contributions to the dielectric permittivity are largely dominated by the dipolar one. We find that this stems from the mutual cancellation of the first two interfacial hydration layers of opposite net charge for multipolar liquids. Within the same non-local dielectric response theory, we show that the induced ion polarizability reverses the interfacial ion density trends predicted by the Poisson-Boltzmann theory, resulting in a surface affinity of coions and exclusion of counterions. The results indicate that the consideration of the discrete charge composition of solvent molecules and ions is the key step towards a microscopic understanding of non-local electrostatic effects in polar solvents.

preprint2012arXiv

Dipolar depletion effect on the differential capacitance of carbon based materials

The remarkably low experimental values of the capacitance data of carbon based materials in contact with water solvent needs to be explained from a microscopic theory in order to optimize the efficiency of these materials. We show that this experimental result can be explained by the dielectric screening deficiency of the electrostatic potential, which in turn results from the interfacial solvent depletion effect driven by image dipole interactions. We show this by deriving from the microscopic system Hamiltonian a non-mean-field dipolar Poisson-Boltzmann equation. This can account for the interaction of solvent molecules with their electrostatic image resulting from the dielectric discontinuity between the solvent medium and the substrate. The predictions of the extended dipolar Poisson-Boltzmann equation for the differential capacitance are compared with experimental data and good agreement is found without any fitting parameters.

preprint2012arXiv

Electrostatic correlations in inhomogeneous charged fluids beyond loop expansion

Electrostatic correlation effects in inhomogeneous symmetric electrolytes are investigated within a previously developed electrostatic self-consistent (SC) theory (R.R. Netz and H. Orland, Eur. Phys.J. E 11, 301 (2003)). To this aim, we introduce two computational approaches that allow to solve the SC equations beyond the loop expansion. Both approaches can handle the case of dielectrically discontinuous boundaries where the one-loop theory is known to fail. By comparing the theoretical results obtained from these schemes with the results of the MC simulations that we ran for ions at neutral single dielectric interfaces as well as with previous MC data for charged interfaces, we first show that the weak coupling (WC) Debye-Huckel (DH) theory remains quantitatively accurate up to the bulk ion density rhob=0.01 M, whereas the SC theory exhibits a good quantitative accuracy up to rhob=0.2 M. Then, we derive from the perturbative SC scheme the one-loop theory of asymmetrically partitioned salt systems around a dielectrically homogeneous charged surface. It is shown that correlation effects originate in these systems from a competition between the salt screening loss at the interface driving the ions to the bulk region, and the interfacial counterion screening excess attracting them towards the surface. In the case of weak surface charges, the interfacial salt screening loss is the dominant effect. As a result, correlations decrease the MF density of both coions and counterions. With increasing surface charge, the surface-attractive counterion screening excess starts to dominate, and correlation effects amplify in this regime the MF density of both type of ions. We also show that at a characteristic value of the electrostatic coupling parameter, electrostatic correlations result in a charge inversion effect.

preprint2012arXiv

Excluded volume effects in macromolecular forces and ion-interface interactions

A charged Yukawa liquid confined in a slit nanopore is studied in order to understand excluded volume effects in the interaction force between the pore walls. A previously developed self-consistent scheme (S. Buyukdagli el al., J. Stat. Mech. P05033 (2011)) and a new simpler variational procedure that self-consistently couple image forces, surface charge induced electric field and pore modified core interactions are used to this aim. For neutral pores, it is shown that with increasing pore size, the theory predicts a transition of the interplate pressure from an attractive to a strongly repulsive regime associated with an ionic packing state, an effect observed in previous Monte Carlo simulations for hard core charges. The role of the range of core interactions in the ionic rejection and interplate pressure is thoroughly analyzed. We show that the physics of the system can be split into two screening regimes. The ionic packing effect takes place in the regime of moderately screened core interactions characterized with the bare screening parameter of the Yukawa potential b<3/lB, where lB is the Bjerrum length. In the second regime of strongly screened core interactions b>3/lB, solvation forces associated with these interactions positively contribute to the ionic rejection driven by electrostatic forces and enhance the magnitude of the attractive pressure. The pronounced dependence of the interplate pressure and ionic partition coefficients on the magnitude and the range of core interactions indicates excluded volume effects as an important ion specificity and a non-negligible ingredient for the stability of macromolecules in electrolyte solutions.

preprint2011arXiv

Ion size effects at ionic exclusion from dielectric interfaces and slit nanopores

A previously developed field-theoretic model [R.D. Coalson et al., J. Chem. Phys. 102, 4584 (1995)] that treats core collisions and Coulomb interactions on the same footing is investigated in order to understand ion size effects on the partition of neutral and charged particles at planar interfaces and the ionic selectivity of slit nanopores. We introduce a variational scheme that can go beyond the mean-field (MF) regime and couple in a consistent way pore modified core interactions, steric effects, electrostatic solvation and image-charge forces, and surface charge induced electrostatic potential. We show that in the dilute limit, the MF and the variational theories agree well with MC simulation results, in contrast to a recent RPA method. The partition of charged Yukawa particles at a neutral dielectric interface (e.g air-water or protein-water interface) is investigated. It is shown that as a result of the competition between core collisions that push the ions towards the surface, and repulsive solvation and image forces that exclude them from the interface, a concentration peak of finite size ions sets in close to the dielectric interface. We also characterize the role played by the ion size on the ionic selectivity of neutral slit nanopores. We show that the complex interplay between electrostatic forces, excluded volume effects induced by core collisions and steric effects leads to an unexpected reversal in the ionic selectivity of the pore with varying pore size: while large pores exhibits a higher conductivity for large ions, narrow pores exclude large ions more efficiently than small ones.

preprint2011arXiv

Ionic exclusion phase transition in neutral and weakly charged cylindrical nanopores

A field theoretic variational approach is introduced to study ion penetration into water-filled cylindrical nanopores in equilibrium with a bulk reservoir. It is shown that (i) an ion located in a neutral pore undergoes two opposing mechanisms: the deformation of its surrounding ionic cloud of opposite charge, with respect to the reservoir, which increases the surface tension and tends to exclude ions form the pore, and (ii) the attractive contribution to the ion self-energy of the repulsive image forces associated with the dielectric jump between the solvent and the pore wall, which becomes more and more screened when ions enter the pore. For pore radii around 1 nm and bulk concentrations lower than 0.2 mol/L, this mechanism leads to a first-order phase transition, similar to capillary "evaporation", from an ionic-penetration state to an ionic-exclusion state. The discontinuous phase transition survives within the biological concentration range (~ 0.15 mol/L) for small enough membrane dielectric permittivities (< 5). In the case of a weakly charged pore, counterion penetration exhibits a non-monotonic behaviour and is characterized by two regimes: at low reservoir concentration or small pore radii, coions are excluded and the energy barrier for counterions is reduced but image forces remain strong enough, so that the counterion partition coefficient in the pore decreases with increasing reservoir concentration up to a characteristic value; for larger reservoir concentrations, image forces are screened and the partition coefficient of counterions increases with the reservoir electrolyte concentration.

preprint2010arXiv

Ionic Capillary Evaporation in Weakly Charged Nanopores

Using a variational field theory, we show that an electrolyte confined to a neutral cylindrical nanopore traversing a low dielectric membrane exhibits a first-order ionic liquid-vapor pseudo-phase-transition from an ionic-penetration "liquid" phase to an ionic-exclusion "vapor" phase, controlled by nanopore-modified ionic correlations and dielectric repulsion. For weakly charged nanopores, this pseudotransition survives and may shed light on the mechanism behind the rapid switching of nanopore conductivity observed in experiments.

preprint2010arXiv

Mapping between the order of thermal denaturation and the shape of the critical line of mechanical unzipping in 1-dimensional DNA models

In this Letter, we investigate the link between thermal denaturation and mechanical unzipping for two models of DNA, namely the Dauxois-Peyrard-Bishop model and a variant thereof we proposed recently. We show that the critical line that separates zipped from unzipped DNA sequences in mechanical unzipping experiments is a power-law in the temperature-force plane. We also prove that for the investigated models the corresponding critical exponent is proportional to the critical exponent alpha, which characterizes the behaviour of the specific heat in the neighbourhood of the critical temperature for thermal denaturation.

preprint2010arXiv

On the application of Canonical Perturbation Theory up to the dissociation threshold

We investigate a model system consisting of a Morse oscillator strongly coupled to a doubly-degenerate bending degree of freedom and show that Canonical Perturbation Theory is able to provide a fairly precise, though not exact, approximation of the Hamiltonian up to the dissociation threshold. Quantum mechanical results and classical ones are discussed in this Letter.

preprint2010arXiv

Origin of negative differential thermal resistance in a chain of two weakly coupled nonlinear lattices

Negative differential resistance in electronic conduction has been extensively studied, but it is not the case for its thermal counterpart, namely, negative differential thermal resistance (NDTR). We present a classical Landauer formula in which the nonlinearity is incorporated by the self-consistent phonon theory in order to study the heat flux across a chain consisting of two weakly coupled lattices. Two typical nonlinear models of hard and soft on-site potentials are discussed, respectively. It is shown that the nonlinearity has strong impacts on the occurring of NDTR. As a result, a transition from the absence to the presence of NDTR is observed. The origin of NDTR consists in the competition between the temperature difference, which acts as an external field, and the temperature-dependent thermal boundary conductance. Finally, the onset of the transition is clearly illustrated for this model. Our analytical calculation agrees reasonably well with numerical simulations.

preprint2010arXiv

Thermal conductivity of anharmonic lattices: Effective phonons and quantum corrections

We compare two effective phonon theories, which have both been applied recently to study heat conduction in anharmonic lattices. In particular, we study the temperature dependence of the thermal conductivity of the Fermi-Pasta-Ulam model via the Debye formula, showing the equivalence of both approaches. The temperature for the minimum of the thermal conductivity and the corresponding scaling behavior are analytically calculated, which agree well with the result obtained from nonequilibrium simulations. We also give quantum corrections for the thermal conductivity from quantum self-consistent phonon theory. The vanishing behavior at the low temperature regime and the existence of an umklapp peak are qualitatively consistent with experimental studies.

preprint2010arXiv

Variational approach for electrolyte solutions: from dielectric interfaces to charged nanopores

A variational theory is developed to study electrolyte solutions, composed of interacting point-like ions in a solvent, in the presence of dielectric discontinuities and charges at the boundaries. Three important and non-linear electrostatic effects induced by these interfaces are taken into account: surface charge induced electrostatic field, solvation energies due to the ionic cloud, and image charge repulsion. Our variational equations thus go beyond the mean-field theory. The influence of salt concentration, ion valency, dielectric jumps, and surface charge is studied in two geometries. i) A single neutral air-water interface with an asymmetric electrolyte. A charge separation and thus an electrostatic field gets established due to the different image charge repulsions for coions and counterions. Both charge distributions and surface tension are computed and compared to previous approximate calculations. For symmetric electrolyte solutions close to a charged surface, two zones are characterized. In the first one, with size proportional to the logarithm of the coupling parameter, strong image forces impose a total ion exclusion, while in the second zone the mean-field approach applies. ii) A symmetric electrolyte confined between two dielectric interfaces as a simple model of ion rejection from nanopores. The competition between image charge repulsion and attraction of counterions by the membrane charge is studied. For small surface charge, the counterion partition coefficient decreases with increasing pore size up to a critical pore size, contrary to neutral membranes. For larger pore sizes, the whole system behaves like a neutral pore. The prediction of the variational method is also compared with MC simulations and a good agreement is observed.

preprint2007arXiv

1/f fluctuations of DNA temperature at thermal denaturation

We theoretically investigated the temperature fluctuations of DNA close to denaturation and observed a strong enhancement of these fluctuations at the critical temperature. Although in a much lower frequency range, such a sharp increase was also reported in the recent experimental work of Nagapriya et al [Phys. Rev. Lett. 96, 038102 (2006)]. We showed that there is instead no enhancement of temperature fluctuations when the dissipation coefficient in Langevin equations is assumed to be larger than a few tens of ps-1, and pointed out the possible role of the solvent in real experiments. We sought for a possible correlation between the growth of large bubbles and the enhancement of temperature fluctuations but found no direct evidence thereof. Finally, we showed that neither the enhancement of fluctuations nor the 1/f dependence are observed at the scale of a single base pair, while these properties show up when summing the contributions of a large number of base pairs. We therefore conclude that both effects result from collective motions that are facilitated by the divergence of the correlation length at denaturation.

preprint2007arXiv

Theoretical investigation of finite size effects at DNA melting

We investigated how the finiteness of the length of the sequence affects the phase transition that takes place at DNA melting temperature. For this purpose, we modified the Transfer Integral method to adapt it to the calculation of both extensive (partition function, entropy, specific heat, etc) and non-extensive (order parameter and correlation length) thermodynamic quantities of finite sequences with open boundary conditions, and applied the modified procedure to two different dynamical models. We showed that rounding of the transition clearly takes place when the length of the sequence is decreased. We also performed a finite-size scaling analysis of the two models and showed that the singular part of the free energy can indeed be expressed in terms of an homogeneous function. However, both the correlation length and the average separation between paired bases diverge at the melting transition, so that it is no longer clear to which of these two quantities the length of the system should be compared. Moreover, Josephson's identity is satisfied for none of the investigated models, so that the derivation of the characteristic exponents which appear, for example, in the expression of the specific heat, requires some care.

preprint2006arXiv

On scaling laws at the phase transition of systems with divergent order parameter and/or internal length : the example of DNA denaturation

We used the Transfer-Integral method to compute, with an uncertainty smaller than 5%, the six fundamental characteristic exponents of two dynamical models for DNA thermal denaturation and investigate the validity of the scaling laws. Doubts concerning this point arise because the investigated systems (i) have a divergent internal length, (ii) are described by a divergent order parameter, (iii) are of dimension 1. We found that the assumption that the free energy can be described by a single homogeneous function is robust, despite the divergence of the order parameter, so that Rushbrooke's and Widom's identities are valid relations. Josephson's identity is instead not satisfied. This is probably due to the divergence of the internal length, which invalidates the assumption that the correlation length is solely responsible for singular contributions to thermodynamic quantities. Fisher's identity is even wronger. We showed that this is due to the d=1 dimensionality and obtained an alternative law, which is well satisfied at DNA thermal denaturation.

preprint2005arXiv

Towards more realistic dynamical models for DNA secondary structure

We propose a dynamical model for the secondary structure of DNA, which is based on the finite stacking enthalpies used in thermodynamics calculations. In this model, the two strands can separate and the bases are allowed to rotate perpendicular to the sequence axis. We show, through molecular dynamics simulations, that the model has the correct behaviour at the denaturation transition.