Source author record

Rudolf Podgornik

Rudolf Podgornik 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

72works
13topics
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

72 published item(s)

preprint2022arXiv

Continuum theories of structured dielectrics

Aequous dielectrics are ubiquitous in soft- and bio-nano matter systems. The theoretical description of such systems in terms of continuum (`macroscopic') theory remains a serious challenge. In this perspective we first review the existing continuum phenomenological approaches that have been developed in the past decades. In order to describe a path to advance continuum theory beyond these approaches we then take recourse to the Onsager-Dupuis theory of the dielectric behaviour of ice, which, for the case of a solid dielectric, exemplified important conceptual issues we deem relevant for the development of a more fundamental continuum theory of liquid dielectrics. Subsequently, we discuss our recently proposed continuum field theory of structured dielectrics, which provides a generalized approach to the dielectric behavior of such systems.

preprint2022arXiv

Curvature effects in charge-regulated lipid bilayers

We formulate a theory of electrostatic interactions in lipid bilayer membranes where both monolayer leaflets contain dissociable moieties that are subject to charge regulation. We specifically investigate the coupling between membrane curvature and charge regulation of a lipid bilayer vesicle using both the linear Debye-Hückel (DH) and the non-linear Poisson-Boltzmann (PB) theory. We find that charge regulation of an otherwise symmetric bilayer membrane can induce charge symmetry breaking, non-linear flexoelectricity and anomalous curvature dependence of free energy. The pH effects investigated go beyond the paradigm of electrostatic renormalization of the mechano-elastic properties of membranes.

preprint2022arXiv

Field theory of structured liquid dielectrics

We develop of a field-theoretic approach for the treatment of both the non-local and the non-linear response of structured liquid dielectrics. Our systems of interest are composed of dipolar solvent molecules and simple salt cations and anions. We describe them by two independent order parameters, the polarization field for the solvent and the charge density field for the ions, including and treating the non-electrostatic part of the interactions explicitly and consistently. We show how to derive functionals for the polarization and the electrostatic field of increasingly finer scales and solve the resulting mean-field saddle-point equations in the linear regime. We derive criteria for the character of their solutions that depend on the structural lengths and the polarity of the solvent. Our approach provides a systematic way to derive generalized polarization theories.

preprint2021arXiv

Embarras de richesses in non-DLVO colloid interactions

In its original formulation, the seminal Deryaguin-Landau-Verwey-Overbeek (DLVO) theory of colloidal stability seemed like a simple but realistic description of the world of colloid interactions in electrolyte solutions. It is based on a straightforward superposition of the mean-field Poisson-Boltzmann (PB) electrostatics with the electrodynamic van der Waals (vdW) interactions driven by thermal and quantum fluctuations. However, subsequent developments continued to reveal a much richer and deeper structure of fundamental interactions on the nano- and micro-scale: the granularity and structure of the solvent, charging equilibria of dissociable charge groups, inhomogeneous charge distributions, the finite size of the ions, non-mean-field electrostatics, ion-ion correlations, and more. Today, the original simplicity is gone and we are left with an embarrassingly rich variety of interactions that defy simple classification and reduction to a few fundamental mechanisms. In this mini-review, we comment on the contemporary state-of-the-art picture of colloidal interactions, in view of some recent progress in experiments.

preprint2021arXiv

Theory of inhomogeneous rod-like Coulomb fluids

A field theoretic representation of the classical partition function is derived for a system composed of a mixture of anisotropic and isotropic mobile charges that interact {\sl via} long range Coulomb and short range nematic interactions. The field theory is then solved on a saddle-point approximation level, leading to a coupled system of Poisson-Boltzmann and Maier-Saupe equations. Explicit solutions are finally obtained for a rod-like counterion-only system in proximity of a charged planar wall, generalizing the standard Gouy-Chapman results. The nematic order parameter profile, the counterion density profile and the electrostatic potential profile are interpreted within the framework of a nematic wetting layer with a (Donnan) potential difference.

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

Thermal Casimir interactions for higher derivative field Lagrangians: generalized Brazovskii models

We examine the Casimir effect for free statistical field theories which have Hamiltonians with second order derivative terms. Examples of such Hamiltonians arise from models of non-local electrostatics, membranes with non-zero bending rigidities and field theories of the Brazovskii type that arise for polymer systems. The presence of a second derivative term means that new types of boundary conditions can be imposed, leading to a richer phenomenology of interaction phenomena. In addition zero modes can be generated that are not present in standard first derivative models, and it is these zero modes which give rise to long range Casimir forces. Two physically distinct cases are considered: (i) unconfined fields, usually considered for finite size embedded inclusions in an infinite fluctuating medium, here in a two plate geometry the fluctuating field exists both inside and outside the plates, (ii) confined fields, where the field is absent outside the slab confined between the two plates. We show how these two physically distinct cases are mathematically related and discuss a wide range of commonly applied boundary conditions. We concentrate our analysis to the critical region where the underlying bulk Hamiltonian has zero modes and show that very exotic Casimir forces can arise, characterised by very long range effects and oscillatory behavior that can lead to strong metastability in the system.

preprint2016arXiv

Chirality modifies the interaction between knots

In this study we consider an idealization of a typical optical tweezers experiment involving a semiflexible double-knotted polymer, with steric hindrance and persistence length matching those of dsDNA in high salt concentration, under strong stretching. Using exhaustive Molecular Dynamics simulations we show that not only does a double-knotted dsDNA filament under tension possess a free energy minimum when the two knots are intertwined, but also that the depth of this minimum depends on the relative chirality of the two knots. We rationalize this dependence of the effective interaction on the chirality in terms of a competition between chain entropy and bending energy.

preprint2016arXiv

Effects of RNA branching on the electrostatic stabilization of viruses

Many single-stranded (ss) RNA viruses self assemble from capsid protein subunits and the nucleic acid to form an infectious virion. It is believed that the electrostatic interactions between the negatively charged RNA and the positively charged viral capsid proteins drive the encapsidation, although there is growing evidence that the sequence of the viral RNA also plays a role in packaging. In particular the sequence will determine the possible secondary structures that the ssRNA will take in solution. In this work, we use a mean field theory to investigate how the secondary structure of the RNA combined with electrostatic interactions affects the efficiency of assembly and stability of the assembled virions. We show that the secondary structure of RNA may result in negative osmotic pressures while a linear polymer causes positive osmotic pressures for the same conditions. This may suggest that the branched structure makes the RNA more effectively packaged and the virion more stable.

preprint2016arXiv

From polymers to proteins: effect of side chains and broken symmetry in the formation of secondary structures within a Wang-Landau approach

We study the equilibrium properties of a single flexible homopolymers where consecutive monomers are represented by impenetrable hard spherical beads tangential to each other, and non-consecutive monomers interact via a square-well potential. To this aim, we use both replica exchange canonical simulations and micro-canonical Wang-Landau techniques. We perform a close comparative analysis of the corresponding results, and find perfect agreement between the two methods as well as the past results. The model is then refined in two different directions. By allowing partial overlapping between consecutive beads, we break the spherical symmetry and thus provide a severe constraint on the possible conformations of the chain. This leads to a single helix or a double helix ground state, depending on the range of the interactions. Alternatively, we introduce additional spherical beads at specific positions to represent the steric hindrance of the side chains in real proteins. The ground state in this case is found to be again a globule, thus showing that contrary to some recent suggestions, side chains alone are not sufficient to insure a protein-like structure. Finally, we show that a combination of the above two effects increases the stability of the obtaining secondary structures, in agreement with previous results. The fundamental role played by the range of the square-well attraction is highlighted, and it is shown that they play a role similar to that found in simple liquids and polymers. Perspectives in terms of protein folding are finally discussed.

preprint2016arXiv

Modulation of elasticity and interactions in charged lipid multibilayers: monovalent salt solutions

We have studied the electrostatic screening effect of NaCl solutions on the interactions between anionic lipid bilayers in the fluid lamellar phase using a Poisson-Boltzmann based mean-field approach with constant charge and constant potential limiting charge regulation boundary conditions. The full DLVO potential, including the electrostatic, hydration and van der Waals interactions, was coupled to thermal bending fluctuations of the membranes via a variational Gaussian Ansatz. This allowed us to analyze the coupling between the osmotic pressure and the fluctuation amplitudes and compare them both simultaneously with the measured dependence on the bilayer separation, determined by the small-angle X-ray scattering experiments. High-structural resolution analysis of the scattering data revealed no significant changes of membrane structure as a function of salt concentration. Parsimonious description of our results is consistent with the constant charge limit of the general charge regulation phenomenology, with fully dissociated lipid charge groups, together with a four-fold reduction of the membranes' bending rigidity upon increasing NaCl concentration.

preprint2016arXiv

Nonequilibrium Tuning of the Thermal Casimir Effect

In net-neutral systems correlations between charge fluctuations generate strong attractive thermal Casimir forces and engineering these forces to optimize nanodevice performance is an important challenge. We show how the normal and lateral thermal Casimir forces between two plates containing Brownian charges can be modulated by decorrelating the system through the application of an electric field, which generates a nonequilibrium steady state with a constant current in one or both plates, reducing the ensuing fluctuation-generated normal force while at the same time generating a lateral drag force. This hypothesis is confirmed by detailed numerical simulations as well as an analytical approach based on stochastic density functional theory.

preprint2016arXiv

Pseudo-Casimir stresses and elasticity of a confined elastomer film

Investigations of the elastic behavior of bulk elastomers have traditionally proceeded on the basis of classical rubber elasticity, which regards chains as thermally fluctuating but disregards the thermal fluctuations of the cross-links. Here, we consider an incompressible and flat elastomer film of an axisymmetric shape confined between two large hard co-planar substrates, with the axis of the film perpendicular to the substrates. We address the impact that thermal fluctuations of the cross-links have on the free energy of elastic deformation of the system, subject to the requirement that the fluctuating elastomer cannot detach from the substrates. We examine the behavior of the deformation free energy for one case where a rigid pinning boundary condition is applied to a class of elastic fluctuations at the confining surfaces, and another case where the same elastic fluctuations are subjected to soft "gluing" potentials. We find that there can be significant departures (both quantitative and qualitative) from the prediction of classical rubber elasticity theory when elastic fluctuations are included. Finally, we compare the character of the attractive part of the elastic fluctuation-induced, or pseudo-Casimir, stress with the standard thermal Casimir stress in confined but non-elastomeric systems, finding the same power law decay behavior when a rigid pinning boundary condition is applied, for the case of the gluing potential, we find that the leading order correction to the attractive part of the fluctuation stress decays inversely with the inter-substrate separation.

preprint2016arXiv

RNA topology remolds electrostatic stabilization of viruses

Simple RNA viruses efficiently encapsulate their genome into a nano-sized protein shell: the capsid. Spontaneous co-assembly of the genome and the capsid proteins is driven predominantly by electrostatic interactions between the negatively charged RNA and the positively charged inner capsid wall. Using field theoretic formulation we show that the inherently branched RNA secondary structure allows viruses to maximize the amount of encapsulated genome and make assembly more efficient, allowing viral RNAs to out-compete cellular RNAs during replication in infected host cells.

preprint2016arXiv

Titratable Macroions in Multivalent Electrolyte Solutions: Strong Coupling Dressed Ion Approach

We present a theoretical description of the effect of polyvalent ions on the interaction between titratable macro-ions. The model system consists of two point-like macro-ions with dissociable sites, immersed in an asymmetric ionic mixture of monovalent and polyvalent salts. We formulate a {\em dressed ion strong coupling theory}, based on the decomposition of the asymmetric ionic mixture into a weakly electrostatically coupled monovalent salt, and into polyvalent ions that are strongly electrostatically coupled to the titratable macro-ions. The charge of the macroions is not considered as fixed, but is allowed to respond to local bathing solution parameters (electrostatic potential, $pH$ of the solution, salt concentration) through a simple {\em charge regulation} model. The approach presented, yielding an effective polyvalent-ion mediated interaction between charge-regulated macro-ions at various solution conditions, describes the strong coupling equivalent of the Kirkwood-Schumaker interaction.

preprint2016arXiv

Van der Waals torque and force between dielectrically anisotropic layered media

We analyse van der Waals interactions between a pair of dielectrically anisotropic plane-layered media interacting across a dielectrically isotropic solvent medium. We develop a general formalism based on transfer matrices to investigate the van der Waals torque and force in the limit of weak birefringence and dielectric matching between the ordinary axes of the anisotropic layers and the solvent. We apply this formalism to study the following systems: (i) a pair of single anisotropic layers, (ii) a single anisotropic layer interacting with a multilayered slab consisting of alternating anisotropic and isotropic layers, and (iii) a pair of multilayered slabs each consisting of alternating anisotropic and isotropic layers, looking at the cases where the optic axes lie parallel and/or perpendicular to the plane of the layers. For the first case, the optic axes of the oppositely facing anisotropic layers of the two interacting slabs generally possess an angular mismatch, and within each multilayered slab the optic axes may either be the same, or undergo constant angular increments across the anisotropic layers. In particular, we examine how the behaviors of the van der Waals torque and force can be "tuned" by adjusting the layer thicknesses, the relative angular increment within each slab, and the angular mismatch between the slabs.

preprint2015arXiv

Bending rigidities and interdomain forces in membranes with coexisting lipid domains

In order to precisely quantify the fundamental interactions between heterogeneous lipid membranes with coexisting liquid-ordered (Lo) and liquid-disordered (Ld) domains, we performed detailed osmotic stress SAXS experiments by exploiting the domain alignment in raft-mimicking lipid multibilayers. Performing a Monte Carlo (MC) based analysis allowed us to determine with high reliability the magnitude and functional dependence of interdomain forces concurrently with the bending elasticity moduli. In contrast to previous methodologies, this approach enabled us to consider the entropic undulation repulsions on a fundamental level, without having to take recourse to crudely justified mean-field like additivity assumptions. Our detailed Hamaker coefficient calculations indicated only small differences in the van der Waals attractions of coexisting Lo and Ld phases. In contrast, the repulsive hydration and undulation interactions differed significantly, with the latter dominating the overall repulsions in the Ld phase. Therefore, alignment of like domains in multibilayers appears to originate from both, hydration and undulation repulsions.

preprint2015arXiv

Effective interactions between fluid membranes

A self-consistent theory is proposed for the general problem of interacting undulating fluid membranes subject to the constraint that they do not interpenetrate. We implement the steric constraint via an exact functional integral representation, and through the use of a saddle-point approximation transform it into a novel effective steric potential. The steric potential is found to consist of two contributions: one generated by zero mode fluctuations of the membranes, and the other by thermal bending fluctuations. For membranes of cross-sectional area $S$, we find that the bending fluctuation part scales with the inter-membrane separation $d$ as $d^{-2}$ for $d \ll \sqrt{S}$, but crosses over to $d^{-4}$ scaling for $d \gg \sqrt{S}$, whereas the zero mode part of the steric potential always scales as $d^{-2}$. For membranes interacting exclusively via the steric potential, we obtain closed-form expressions for the effective interaction potential and for the rms undulation amplitude $σ$, which becomes small at low temperatures $T$ and/or large bending stiffnesses $κ$. Moreover, $σ$ scales as $d$ for $d \ll \sqrt{S}$, but saturates at $\sqrt{k_{\rm B} T S/κ}$ for $d \gg \sqrt{S}$. In addition, using variational Gaussian theory, we apply our self-consistent treatment to study inter-membrane interactions subject to three different types of potential: (i)~the Moreira-Netz potential for a pair of strongly charged membranes with an intervening solution of multivalent counterions, (ii)~an attractive square well, (iii)~the Morse potential, and (iv)~a combination of hydration and van der Waals interactions.

preprint2015arXiv

Ion-mediated interactions between net-neutral slabs: Weak and strong disorder effects

We investigate the effective interaction between two randomly charged but otherwise net-neutral, planar dielectric slabs immersed in an asymmetric Coulomb fluid containing a mixture of mobile monovalent and multivalent ions. The presence of charge disorder on the apposed bounding surfaces of the slabs leads to substantial qualitative changes in the way they interact, as compared with the standard picture provided by the van der Waals and image-induced, ion-depletion interactions. While, the latter predict purely attractive interactions between strictly neutral slabs, we show that the combined effects from surface charge disorder, image depletion, Debye (or salt) screening and also, in particular, their coupling with multivalent ions, give rise to a more diverse behavior for the effective interaction between net-neutral slabs. Disorder effects show large variation depending on the properly quantified strength of disorder, leading either to non-monotonic effective interaction with both repulsive and attractive branches when the surface charges are weakly disordered (small disorder variance) or to a dominating attractive interaction that is larger both in its range and magnitude than what is predicted from the van der Waals and image-induced, ion-depletion interactions, when the surfaces are strongly disordered (large disorder variance).

preprint2015arXiv

Phase diagram of the ground states of DNA condensates

Phase diagram of the ground states of DNA in a bad solvent is studied for a semi-flexible polymer model with a generalized local elastic bending potential characterized by a nonlinearity parameter $x$ and effective self-attraction promoting compaction. $x=1$ corresponds to the worm-like chain model. Surprisingly, the phase diagram as well as the transition lines between the ground states are found to be a function of $x$. The model provides a simple explanation for the results of prior experimental and computational studies and makes predictions for the specific geometries of the ground states. The results underscore the impact of the form of the microscopic bending energy at macroscopic observable scales.

preprint2015arXiv

Repulsive Casimir interaction: Boyer oscillators at nanoscale

We study the effect of temperature on the repulsive Casimir interaction between an ideally permeable and an ideally polarizable plate {\it in vacuo}. At small separations or for low temperatures the quantum fluctuations of the electromagnetic field give the main contribution to the interaction, while at large separations or for high temperatures the interaction is dominated by the classical thermal fluctuations of the field. At intermediate separations or finite temperatures both the quantum and thermal fluctuations contribute. For a system composed of one infinitely permeable plate between two ideal conductors at a finite temperature, we identify a {\it stable mechanical equilibrium} state, if the infinitely permeable plate is located in the middle of the cavity. For small displacements the restoring force of this {\it Boyer oscillator} is linear in the deviation from the equilibrium position, with a spring constant that depends inversely on the separation between the two conducting plates and linearly on temperature. Furthermore, an array of such oscillators presents an ideal Einsteinian crystal that displays a fluctuation force between its outer boundaries stemming from the displacement fluctuations of the Boyer oscillators.

preprint2015arXiv

Solvent Effects in the Helix-Coil Transition Model can explain the Unusual Biophysics of Intrinsically Disordered Proteins

We analyze a model statistical description of the polypeptide chain helix-coil transition, where we take into account the specificity of its primary sequence, as quantified by the phase space volume ratio of the number of all accessible states to the number corresponding to a helical conformation. The resulting transition phase diagram is then juxtaposed with the unusual behavior of the secondary structures in Intrinsically Disordered Proteins (IDPs) and a number of similarities are observed, even if the protein folding is a more complex transition than the helix-coil transition. In fact, the deficit in bulky and hydrophobic amino acids observed in IDPs, translated into larger values of phase space volume, allows us to locate the region in parameter space of the helix-coil transition that would correspond to the secondary structure transformations that are intrinsic to conformational transitions in IDPs and that is characterized by a modified phase diagram when compared to globular proteins. Here, we argue how the nature of this modified phase diagram, obtained from a model of the helix-coil transition in a solvent, would illuminate the turned-out response of IDPs to the changes in the environment conditions that follow straightforwardly from the re-entrant (cold denaturation) branch in their folding phase diagram.

preprint2015arXiv

Surface Tension of Electrolyte Interfaces: Ionic Specificity within a Field-Theory Approach

We study the surface tension of ionic solutions at air/water and oil/water interfaces. By using field-theoretical methods and including a finite proximal surface-region with ionic-specific interactions. The free energy is expanded to first-order in a loop expansion beyond the mean-field result. We calculate the excess surface tension and obtain analytical predictions that reunite the Onsager-Samaras pioneering result (which does not agree with experimental data), with the ionic specificity of the Hofmeister series. We derive analytically the surface-tension dependence on the ionic strength, ionic size and ion-surface interaction, and show consequently that the Onsager-Samaras result is consistent with the one-loop correction beyond the mean-field result. Our theory fits well a wide range of salt concentrations for different monovalent ions using one fit parameter per electrolyte, and reproduces the reverse Hofmeister series for anions at the air/water and oil/water interfaces.

preprint2015arXiv

Van der Waals interaction between polymers with sequence-specific polarizabilities: Stiff polymers and flexible Gaussian coils

We consider the van der Waals interaction between a pair of polymers with quenched heterogeneous sequences of local polarizabilities along their backbones, and study the effective pairwise interaction energy for both stiff polymers and flexible Gaussian coils. In particular, we focus on the cases where the pair of polarizability sequences are (i)~distinct and (ii)~identical. We find that the pairwise interaction energies of distinct and identical Gaussian coils are both isotropic and exhibit the same decay behavior for separations larger than the gyration radius, in contradistinction to the orientationally anisotropic interaction energies of distinct and identical stiff polymers. For both Gaussian coils and stiff polymers, the attractive interaction between identical polymers is enhanced if the polarizability sequence is more heterogeneous.

preprint2014arXiv

A unified description of solvent effects in the helix-coil transition

We analyze the problem of the helix-coil transition in explicit solvents analytically by using spin-based models incorporating two different mechanisms of solvent action: explicit solvent action through the formation of solvent-polymer hydrogen bonds that can compete with the intrinsic intra-polymer hydrogen bonded configurations (competing interactions) and implicit solvent action, where the solvent-polymer interactions tune biopolymer configurations by changing the activity of the solvent (non-competing interactions). The overall spin Hamiltonian is comprised of three terms: the background \emph{in vacuo} Hamiltonian of the "Generalized Model of Polypeptide Chain" type and two additive terms that account for the two above mechanisms of solvent action. We show that on this level the solvent degrees of freedom can be {\sl explicitly} and {\sl exactly} traced over, the ensuing effective partition function combining all the solvent effects in a unified framework. In this way we are able to address helix-coil transitions for polypeptides, proteins, and DNA, with different buffers and different external constraints. Our spin-based effective Hamiltonian is applicable for treatment of such diverse phenomena as cold denaturation, effects of osmotic pressure on the cold and warm denaturation, complicated temperature dependence of the hydrophobic effect as well as providing a conceptual base for understanding the behavior of Intrinsically Disordered Proteins and their analogues.

preprint2014arXiv

Asymmetric Coulomb fluids at randomly charged dielectric interfaces: Anti-fragility, overcharging and charge inversion

We study the distribution of multivalent counterions next to a dielectric slab, bearing a quenched, random distribution of charges on one of its solution interfaces, with a given mean and variance, both in the absence and in the presence of a bathing monovalent salt solution. We use the previously derived approach based on the dressed multivalent-ion theory that combines aspects of the strong and weak coupling of multivalent and monovalent ions in a single framework. The presence of quenched charge disorder on the charged surface of the dielectric slab is shown to substantially increase the density of multivalent counterions in its vicinity. In the counterion-only model (with no monovalent salt ions), the surface disorder generates an additional logarithmic attraction potential and thus an algebraically singular counterion density profile at the surface. This behavior persists also in the presence of a monovalent salt bath and results in significant violation of the contact-value theorem, reflecting the anti-fragility effects of the disorder that drive the system towards a more ordered state. In the presence of an interfacial dielectric discontinuity, depleting the counterion layer at the surface, the charge disorder still generates a much enhanced counterion density further away from the surface. Likewise, the charge inversion and/or overcharging of the surface occur more strongly and at smaller bulk concentrations of multivalent counterions when the surface carries quenched charge disorder. Overall, the presence of quenched surface charge disorder leads to sizable effects in the distribution of multivalent counterions in a wide range of realistic parameters and typically within a distance of a few nanometers from the charged surface.

preprint2014arXiv

Charge regulation in ionic solutions: thermal fluctuations and Kirkwood-Schumaker interaction

We study the behavior of two macroions with dissociable charge groups, regulated by local variables such as pH and electrostatic potential, immersed in a mono-valent salt solution, considering cases where the net charge can either change sign or remain of the same sign depending on these local parameters. The charge regulation, in both cases, is described with the proper free energy function for each of the macroions, while the coupling between the charges is evaluated on the approximate Debye-Hückel level. The charge correlation functions and the ensuing charge fluctuation forces are calculated analytically and numerically. Strong attraction between like-charged macroions is found close to the point of zero charge, specifically due to asymmetric, anticorrelated charge fluctuations of the macroion charges. The general theory is then implemented for a system of two protein-like macroions, generalizing the form and magnitude of the Kirkwood-Schumaker interaction.

preprint2014arXiv

Field theoretic description of charge regulation interaction

In order to find the exact form of the electrostatic interaction between two proteins with dissociable charge groups in aqueous solution, we have studied a model system composed of two macroscopic surfaces with charge dissociation sites immersed in a counterion-only ionic solution. Field-theoretic representation of the grand canonical partition function is derived and evaluated within the mean-field approximation, giving the Poisson-Boltzmann theory with the Ninham-Parsegian boundary condition. Gaussian fluctuations around the mean-field are then analyzed in the lowest order correction that we calculate analytically and exactly, using the path integral representation for the partition function of a harmonic oscillator with time-dependent frequency. The first order (one loop) free energy correction gives the interaction free energy that reduces to the zero-frequency van der Waals form in the appropriate limit but in general gives rise to a mono-polar fluctuation term due to charge fluctuation at the dissociation sites. Our formulation opens up the possibility to investigate the Kirkwood-Shumaker interaction in more general contexts where their original derivation fails.

preprint2014arXiv

Fluctuation-induced interactions in nematics with disordered anchoring energy

We examine fluctuation-induced (pseudo-Casimir) interactions in nematic liquid-crystalline films confined between two surfaces, where one of the surfaces imposes a strong homeotropic anchoring (ensuring a uniform mean director profile), while the other one is assumed to be a chemically disordered substrate exhibiting an annealed, random distribution of anchoring energies. We employ a saddle-point approximation to evaluate the free energy of interaction mediated between the two surfaces and investigate how the interaction force is influenced by the presence of disordered surface anchoring energy. It is shown that the disorder results in a renormalization of the effective surface anchoring parameter in a way that it leads to quantitative and qualitative changes (including a change of sign at intermediate inter-surface separations) in the pseudo-Casimir interaction force when compared with the interaction force in the absence of disorder.

preprint2014arXiv

From toroidal to rod-like condensates of semiflexible polymers

The competition between toroidal and rod-like conformations as possible ground states for DNA condensation is studied as a function of the stiffness, the length of the DNA and the form of the long-range interactions between neighboring molecules, using analytical theory supported by Monte Carlo simulations. Both conformations considered are characterized by a local nematic order with hexagonal packing symmetry of neighboring DNA molecules, but differ in global configuration of the chain and the distribution of its curvature as it wraps around to form a condensate. The long-range interactions driving the DNA condensation are assumed to be of the form pertaining to the attractive depletion potential as well as the attractive counterion induced soft potential. In the stiffness-length plane we find a transition between rod-like to toroid condensate for increasing stiffness at a fixed chain length $L$. Strikingly, the transition line is found to have a $L^{1/3}$ dependence irrespective of the details of the long-range interactions between neighboring molecules. When realistic DNA parameters are used, our description reproduces rather well some of the experimental features observed in DNA condensates.

preprint2014arXiv

Kirkwood-Shumaker interactions and general thermal fluctuation forces

We present an extension of the Kirkwood-Shumaker (KS) theory of proton-fluctuation interactions to situations where the perturbation theory, usually invoked to derive these interactions, fails. In order to do that we formulate a generalized theory of fluctuation interactions with non-linear macroion (surface) free energy term that naturally leads to the long-range KS interactions, but can also be straightforwardly generalized to situations where the perturbation expansion ceases to be valid. It also allows for other surface free energy terms to be treated within the same formalism, enabling us to derive the complete fluctuation interaction with monopolar and dipolar fluctuations included, leading to the KS theory as a limiting form.

preprint2014arXiv

Relaxation of the thermal Casimir force between net neutral plates containing Brownian charges

We investigate the dynamics of thermal Casimir interactions between plates described within a living conductor model, with embedded mobile anions and cations, whose density field obeys a stochastic partial differential equation which can be derived starting from the Langevin equations of the individual particles. This model describes the thermal Casimir interaction in the same way that the fluctuating dipole model describes van der Waals interactions. The model is analytically solved in a Debye-Hückel-like approximation. We identify several limiting dynamical regimes where the time dependence of the thermal Casimir interactions can be obtained explicitly. Most notably we find a regime with diffusive scaling, even though the charges are confined to the plates and do not diffuse into the intervening space, which makes the diffusive scaling difficult to anticipate and quite unexpected on physical grounds.

preprint2013arXiv

Electrostatic stability and encapsidation of charged nano-droplets

We investigate electrostatic stability of charged droplets, modeled as permeable, charged spheres, and their encapsidation in thin, arbitrarily charged nano-shells, immersed in a neutralizing asymmetric electrolyte background. The latter consists of a small concentration of mobile multivalent counterions in a bathing solution of monovalent (positive and negative) ions. We use extensive Monte-Carlo simulations to investigate the spatial distribution of multivalent counterions and the electrostatic component of their osmotic pressure on the bounding surface of the spherical nano-shell. The osmotic pressure can be negative (inward pressure), positive (outward pressure) or zero, depending on the system parameters such as the charge density of the droplet, the charge density of the shell, and the electrolyte screening, which thus determine the stability of the nano-container. The counter-intuitive effects of multivalent counterions comprise the increased stability of the charged droplet with larger charge density, increased stability in the case of encapsidating shell of charge density of the same sign as the charged droplet, as well as the possibility to dispense altogether with the encapsidating shell, its confining effect taken over by the multivalent counterions. These dramatic effects are in stark contrast to the conventional mean-field picture, which in particular implies that a more highly charged spherical droplet should be electrostatically less stable because of its larger (repulsive) self-energy.

preprint2013arXiv

Interaction of a point charge with the surface of a uniaxial dielectric

We analyze the force on a point charge moving at relativistic speeds parallel to the surface of a uniaxial dielectric. Two cases are examined: a lossless dielectric with no dispersion and a dielectric with a plasma type response. The treatment focuses on the peculiarities of the strength and direction of the interaction force as compared to the isotropic case. We show that a plasma type dielectric can, under specific conditions, repel the point charge.

preprint2013arXiv

Perspective: Coulomb fluids -- weak coupling, strong coupling, in between and beyond

We present a personal view on the current state of statistical mechanics of Coulomb fluids with special emphasis on the interactions between macromolecular surfaces, concentrating on the weak and the strong coupling limits. Both are introduced for a (primitive) counterion-only system in the presence of macroscopic, uniformly charged boundaries, where they can be derived systematically. Later we show how this formalism can be generalized to the cases with additional characteristic length scales that introduce new coupling parameters into the problem. These cases most notably include asymmetric ionic mixtures with mono- and multivalent ions that couple differently to charged surfaces, ions with internal charge (multipolar) structure and finite static polarizability, where weak and strong coupling limits can be constructed by analogy with the counterion-only case and lead to important new insights into their properties that can not be derived by any other means.

preprint2013arXiv

Pseudo-Casimir interactions across nematic films with disordered anchoring axis

We consider the effective pseudo-Casimir interaction forces mediated by a nematic liquid-crystalline film bounded by two planar surfaces, one of which imposes a random (disordered) distribution of the preferred anchoring axis in the so-called easy direction. We consider both the case of a quenched as well as an annealed disorder for the easy direction on the disordered surface and analyze the resultant fluctuation-induced interaction between the surfaces. In the case of quenched disorder, we show that the disorder effects appear additively in the total interaction and are dominant at intermediate inter-surface separations. Disorder effects are shown to be unimportant at both very small and very large separations. In the case of annealed disorder its effects are non-additive in the total inter-surface interaction and can be rationalized in terms of a renormalized extrapolation length.

preprint2013arXiv

Statistical analysis of sizes and shapes of virus capsids and their resulting elastic properties

From the analysis of sizes of approximately 130 small icosahedral viruses we find that there is a typical structural capsid protein, having a mean diameter of 5 nm and a mean thickness of 3 nm, with more than two thirds of the analyzed capsid proteins having thicknesses between 2 nm and 4 nm. To investigate whether, in addition to the fairly conserved geometry, capsid proteins show similarities in the way they interact with one another, we examined the shapes of the capsids in detail. We classified them numerically according to their similarity to sphere and icosahedron and an interpolating set of shapes in between, all of them obtained from the theory of elasticity of shells. In order to make a unique and straightforward connection between an idealized, numerically calculated shape of an elastic shell and a capsid, we devised a special shape fitting procedure, the outcome of which is the idealized elastic shape fitting the capsid best. Using such a procedure we performed statistical analysis of a series of virus shapes and we found similarities between the capsid elastic properties of even very different viruses. As we explain in the paper, there are both structural and functional reasons for the convergence of protein sizes and capsid elastic properties. Our work presents a specific quantitative scheme to estimate relatedness between different proteins based on the details of the (quaternary) shape they form (capsid). As such, it may provide an information complementary to the one obtained from the studies of other types of protein similarity, such as the overall composition of structural elements, topology of the folded protein backbone, and sequence similarity.

preprint2013arXiv

Symmetry effects in electrostatic interactions between two arbitrarily charged spherical shells in the Debye-Hückel approximation

Inhomogeneous charge distributions have important repercussions on electrostatic interactions in systems of charged particles but are often difficult to examine theoretically. We investigate how electrostatic interactions are influenced by patchy charge distributions exhibiting certain point group symmetries. We derive a general form of the electrostatic interaction energy of two permeable, arbitrarily charged spherical shells in the Debye-Hückel approximation and apply it to the case of particles with icosahedral, octahedral, and tetrahedral inhomogeneous charge distributions. We analyze in detail how charge distribution symmetry modifies the interaction energy and find that local charge inhomogeneities reduce the repulsion of two overall equally charged particles, while sufficient orientational variation in the charge distribution can turn the minimum interaction energy into an attraction. Additionally we show that larger patches and thus lower symmetries and wave numbers result in bigger attraction given the same variation.

preprint2013arXiv

Tensorial conservation law for nematic polymers

We derive the "conservation law" for nematic polymers in tensorial form valid for quadrupolar orientational order in contradistinction to the conservation law in the case of polar orientational order. Due to microscopic differences in the coupling between the orientational field deformations and the density variations for polar and quadrupolar order, we find that respective order parameters satisfy fundamentally distinct constraints. Being necessarily scalar in its form, the tensorial conservation law is obtained straightforwardly from the gradients of the polymer nematic tensor field and connects the spatial variation of this tensor field with density variations. We analyze the differences between the polar and the tensorial forms of the conservation law, present some explicit orientational fields that satisfy this new constraint and discuss the role of singular "hairpins", which do not affect local quadrupolar order of polymer nematics, but nevertheless influence its gradients.

preprint2012arXiv

Aspects of One-Dimensional Coulomb Gases

In this short review, we discuss recent advances in exact solutions of models based on a one- dimensional (1D) Coulomb gas by means of field-theoretic functional integral methods. The exact solutions can be used to assess the accuracy of various approximations such as the weak coupling Poisson-Boltzmann theory as well as the strong coupling theory of Coulomb gases. We consider three different 1D models: the Coulomb fluid configuration in the case of the soap film model consisting of positively and negatively charged particles between adsorbing boundaries, counterions between two charged surfaces, and an ionic liquid lattice capacitor with positively and negatively charged particles on a lattice between one positive and one negative bounding surface.

preprint2012arXiv

Attraction between Neutral Dielectrics Mediated by Multivalent Ions in an Asymmetric Ionic Fluid

We study the interaction between two neutral plane-parallel dielectric bodies in the presence of a highly asymmetric ionic fluid, containing multivalent as well as monovalent (salt) ions. Image charge interactions, due to dielectric discontinuities at the boundaries, as well as effects from ion confinement in the slit region between the surfaces are taken fully into account, leading to image-generated depletion attraction, ion correlation attraction and steric-like repulsive interactions. We investigate these effects by employing a combination of methods including explicit-ion and implicit-ion Monte-Carlo simulations, as well as an effective interaction potential analytical theory. The latter incorporates strong ion-image charge correlations, which develop in the presence of high valency ions in the mixture. We show that the implicit-ion simulations and the proposed analytical theory can describe the explicit simulation results on a qualitative level, while excellent quantitative agreement can be obtained for sufficiently large monovalent salt concentrations. The resultant attractive interaction between the neutral surfaces is shown to be significant, as compared with the usual van der Waals interactions between semi-infinite dielectrics, and can thus play a significant role at the nano scale.

preprint2012arXiv

Chirality-Dependent Properties of Carbon Nanotubes. Electronic Structure: Optical Dispersion Properties, Hamaker Coefficients and van der Waals - London dispersion interactions

Optical dispersion spectra at energies up to 30 eV play a vital role in understanding the chirality-dependent van der Waals London dispersion interactions of single wall carbon nanotubes (SWCNTs). We use one-electron theory based calculations to obtain the band structures and the frequency dependent dielectric response function from 0-30 eV for 64 SWCNTs differing in radius, electronic structure classification, and geometry. The resulting optical dispersion properties can be categorized over three distinct energy intervals (M, pi, and sigma, respectively representing 0-0.1, 0.1-5, and 5-30 eV regions) and over radii above or below the zone-folding limit of 0.7 nm. While pi peaks vary systematically with radius for a given electronic structure type, peaks are independent of tube radius above the zone folding limit and depend entirely on SWCNT geometry. Based on these calculated one-electron dielectric response functions we compute and review Van derWaals - London dispersion spectra, full spectral Hamaker coefficients, and van derWaals - London dispersion interaction energies for all calculated frequency dependent dielectric response functions. Our results are categorized using a new optical dielectric function classification scheme that groups the nanotubes according to observable trends and notable features (e.g. the metal paradox) in the 0-30 eV part of the optical dispersion spectra. While the trends in these spectra begin to break down at the zone folding diameter limit, the trends in the related van derWaals - London dispersion spectra tend to remain stable all the way down to the smallest single wall carbon nanotubes in a given class.

preprint2012arXiv

Coil-globule transition of a homopolymer chain in a square-well potential: Comparison between Monte Carlo canonical replica exchange and Wang-Landau sampling

We study the equilibrium properties of a flexible homopolymer where consecutive monomers are represented by impenetrable hard spheres that are tangent to each other, and non-consecutive monomers interact via a square-well potential. To this aim, we use both replica exchange canonical simulations and micro-canonical Wang-Landau techniques for relatively short chains, and perform a close comparative analysis of the corresponding results. These investigations are then further exploited to reproduce, at a much shorter scale and, hence, computational effort, the phase diagram previously studied with much longer chains. This opens up the possibility of improving the model and introduce specificities typical, among other examples, of protein folding.

preprint2012arXiv

Confined chiral polymer nematics: ordering and spontaneous condensation

We investigate condensation of a long confined chiral nematic polymer inside a spherical enclosure, mimicking condensation of DNA inside a viral capsid. The Landau-de Gennes nematic free energy {\sl Ansatz} appropriate for nematic polymers allows us to study the condensation process in detail with different boundary conditions at the enclosing wall that simulate repulsive and attractive polymer-surface interactions. Increasing the chirality, we observe a transformation of the toroidal condensate into a closed surface with an increasing genus, akin to the ordered domain formation observed in cryo-microscopy of bacteriophages.

preprint2012arXiv

Coulomb interactions between disordered charge distributions

We review some of the recent results in the context of the Coulomb interaction between dielectric surfaces which are randomly charged but remain net-neutral on the average. Such surfaces are found to interact in vacuum with attractive long-range normal forces due to the presence of the quenched charge disorder which can compete with the standard Casimir-Van der Waals forces. They can also interact with random lateral forces and torques, which exhibit zero mean values but finite and long-range root-mean-square values. These effects can play an important role in Casimir experiments as well as in the interaction between solid surfaces as well as biomolecular surfaces that are often covered by disordered charge distributions.

preprint2012arXiv

Electromagnetic fluctuation-induced interactions in randomly charged slabs

Randomly charged net-neutral dielectric slabs are shown to interact across a featureless dielectric continuum with long-range electrostatic forces that scale with the statistical variance of their quenched random charge distribution and inversely with the distance between their bounding surfaces. By accounting for the whole spectrum of electromagnetic field fluctuations, we show that this long-range disorder-generated interaction extends well into the retarded regime where higher-order Matsubara frequencies contribute significantly. This occurs even for highly clean samples with only a trace amount of charge disorder and shows that disorder effects can be important down to the nano scale. As a result, the previously predicted non-monotonic behavior for the total force between dissimilar slabs as a function of their separation distance is substantially modified by higher-order contributions, and in almost all cases of interest, we find that the equilibrium inter-surface separation is shifted to substantially larger values compared to predictions based solely on the zero-frequency component. This suggests that the ensuing non-monotonic interaction is more easily amenable to experimental detection. The presence of charge disorder in the intervening dielectric medium between the two slabs is shown to lead to an additional force that can be repulsive or attractive depending on the system parameters and can, for instance, wash out the non-monotonic behavior of the total force when the intervening slab contains a sufficiently large amount of disorder charges.

preprint2012arXiv

Electrostatic interactions mediated by polarizable counterions: weak and strong coupling limits

We investigate the statistical mechanics of an inhomogeneous Coulomb fluid composed of charged particles with static polarizability. We derive the weak- and the strong-coupling approximations and evaluate the partition function in a planar dielectric slab geometry with charged boundaries. We investigate the density profiles and the disjoining pressure for both approximations. Comparison to the case of non-polarizable counterions shows that polarizability brings important differences in the counterion density distribution as well as the counterion mediated electrostatic interactions between charged dielectric interfaces.

preprint2012arXiv

How simple can a model of an empty viral capsid be? Charge distributions in viral capsids

We investigate and quantify salient features of the charge distributions on viral capsids. Our analysis combines the experimentally determined capsid geometry with simple models for ionization of amino acids, thus yielding the detailed description of spatial distribution for positive and negative charge across the capsid wall. The obtained data is processed in order to extract the mean radii of distributions, surface charge densities and dipole moment densities. The results are evaluated and examined in light of previously proposed models of capsid charge distributions, which are shown to have to some extent limited value when applied to real viruses.

preprint2012arXiv

One dimensional Potts model with many-body interactions and the Generalized Model of Polypeptide Chain for the helix-coil transition

Helix-coil transition in polypeptides is an example of a spin model with a preferred spin direction, in the sense that a theoretical formulation of this problem requires to assign a preferred value of spin to the helical conformation in order to account for different symmetries of the helical {\sl vs.} the coil states. This leads to the spin Hamiltonian of the {\sl Generalized Model of Polypeptide Chain} (GMPC) variety as opposed to the Potts model variety, both with many-body interactions. We compare the explicit solution of the Potts model and the solution of the GMPC within the transfer-matrix formalism. Comparison of both secular equations reveals that the largest eigenvalue of the Potts model with $Δ$ many-body interactions is identical to the largest eigenvalue of the GMPC model with $Δ-1$ many-body interactions, indicating the equivalence of both free energies. In distinction, the second largest eigenvalues do not coincide, leading to different thermal behavior of the spatial correlation length, related to the helix-coil transition interval. Spin models with built-in spin anisotropy thus engender different physical properties in the thermodynamic limit that we explore in detail.

preprint2012arXiv

Osmotic pressure induced coupling between cooperativity and stability of a helix-coil transition

Most helix-coil transition theories can be characterized by a set of three parameters: energetic, describing the (free) energy cost of forming a helical state in one repeating unit; entropic, accounting for the decrease of entropy due to the helical state formation; and geometric, indicating how many repeating units are affected by the formation of one helical state. Depending on their effect on the helix-coil transition, solvents or co-solutes can be classified with respect to their action on these parameters. Solvent interactions that alter the entropic cost of helix formation by their osmotic action can affect both the stability (transition temperature) and the cooperativity (transition interval) of the helix-coil transition. A consistent inclusion of osmotic pressure effects in a description of helix-coil transition for poly(L-glutamic acid) in solution with polyethylene glycol can offer an explanation of the experimentally observed linear dependence of transition temperature on osmotic pressure as well as the concurrent changes in the cooperativity of the transition.

preprint2012arXiv

Out of equilibrium thermal Casimir effect in a model polarizable material

Relaxation of the thermal Casimir or van der Waals force for a model dielectric medium is investigated. We start with a model of interacting polarization fields with a dynamics that leads to a frequency dependent dielectric constant of the Debye form. In the static limit the usual zero frequency Matsubara mode component of the Casimir force is recovered. We then consider the out of equilibrium relaxation of the van der Waals force to its equilibrium value when two initially uncorrelated dielectric bodies are brought into sudden proximity. It is found that the spatial dependence of the out of equilibrium force is the same as the equilibrium one but it has a time dependent amplitude, or Hamaker coefficient, which increases in time to its equilibrium value. The final relaxation to the equilibrium value is exponential in systems with a single or finite number of polarization field relaxation times. However, in systems, such as those described by the Havriliak-Negami dielectric constant, with a broad distribution of relaxation times, we observe a much slower power law decay to the equilibrium value.

preprint2012arXiv

Polymers pushing Polymers: Polymer Mixtures in Thermodynamic Equilibrium with a Pore

We investigate polymer partitioning from polymer mixtures into nanometer size cavities by formulating an equation of state for a binary polymer mixture assuming that only one (smaller) of the two polymer components can penetrate the cavity. Deriving the partitioning equilibrium equations and solving them numerically allows us to introduce the concept of "polymers-pushing-polymers" for the action of non-penetrating polymers on the partitioning of the penetrating polymers. Polymer partitioning into a pore even within a very simple model of a binary polymer mixture is shown to depend in a complicated way on the composition of the polymer mixture and/or the pore-penetration penalty. This can lead to enhanced as well as diminished partitioning, due to two separate energy scales that we analyse in detail.

preprint2012arXiv

Sample-to-sample torque fluctuations in a system of coaxial randomly charged surfaces

Polarizable randomly charged dielectric objects have been recently shown to exhibit long-range lateral and normal interaction forces even when they are effectively net neutral. These forces stem from an interplay between the quenched statistics of random charges and the induced dielectric image charges. This type of interaction has recently been evoked to interpret measurements of Casimir forces in vacuo, where a precise analysis of such disorder-induced effects appears to be necessary. Here we consider the torque acting on a randomly charged dielectric surface (or a sphere) mounted on a central axle next to another randomly charged surface and show that although the resultant mean torque is zero, its sample-to-sample fluctuation exhibits a long-range behavior with the separation distance between the juxtaposed surfaces and that, in particular, its root-mean-square value scales with the total area of the surfaces. Therefore, the disorder-induced torque between two randomly charged surfaces is expected to be much more pronounced than the disorder-induced lateral force and may provide an effective way to determine possible disorder effects in experiments, in a manner that is independent of the usual normal force measurement.

preprint2012arXiv

The one dimensional Coulomb lattice fluid capacitor

The one dimensional Coulomb lattice fluid in a capacitor configuration is studied. The model is formally exactly soluble via a transfer operator method within a field theoretic representation of the model. The only interactions present in the model are the one dimensional Coulomb interaction between cations and anions and the steric interaction imposed by restricting the maximal occupancy at any lattice site to one particle. Despite the simplicity of the model, a wide range of intriguing physical phenomena arise, some of which are strongly reminiscent of those seen in experiments and numerical simulations of three dimensional ionic liquid based capacitors. Notably we find regimes where over-screening and density oscillations are seen near the capacitor plates. The capacitance is also shown to exhibit strong oscillations as a function of applied voltage. It is also shown that the corresponding mean field theory misses most of these effects. The analytical results are confirmed by extensive numerical simulations.

preprint2012arXiv

Wrapping Transition and Wrapping-Mediated Interactions for Discrete Binding along an Elastic Filament: An Exact Solution

The wrapping equilibria of one and two adsorbing cylinders are studied along a semi-flexible filament (polymer) due to the interplay between elastic rigidity and short-range adhesive energy between the cylinder and the filament. We show that statistical mechanics of the system can be solved exactly using a path integral formalism which gives access to the full effect of thermal fluctuations, going thus beyond the usual Gaussian approximations which take into account only the contributions from the minimal energy configuration and small fluctuations about this minimal energy solution. We obtain the free energy of the wrapping-unwrapping transition of the filament around the cylinders as well as the effective interaction between two wrapped cylinders due to thermal fluctuations of the elastic filament. A change of entropy due to wrapping of the filament around the adsorbing cylinders as they move closer together is identified as an additional source of interactions between them. Such entropic wrapping effects should be distinguished from the usual entropic configuration effects in semi-flexible polymers. Our results may be applicable to the problem of adsorption of proteins as well as synthetic nano-particles on semi-flexible polymers such as DNA.

preprint2011arXiv

Dressed Counterions: Poly- and monovalent ions at charged dielectric interfaces

We investigate the ion distribution and overcharging at charged interfaces with dielectric inhomogeneities in the presence of asymmetric electrolytes containing polyvalent and monovalent ions. We formulate an effective "dressed counterion" approach by integrating out the monovalent salt degrees of freedom and show that it agrees with results of explicit Monte-Carlo simulations. We then apply the dressed counterion approach within the framework of the strong-coupling theory, valid for polyvalent ions at low concentrations, which enables an analytical description for salt effects as well as dielectric inhomogeneities in the limit of strong Coulomb interactions on a systematic level. Limitations and applicability of this theory are examined by comparing the results with simulations.

preprint2011arXiv

Electrostatic self-energy of a partially formed spherical shell in salt solution: application to stability of tethered and fluid shells -- viruses and vesicles

We investigate the electrostatics of a partially formed, charged spherical shell in a salt solution. We solve the problem numerically at the Poisson-Boltzmann level and analytically in the Debye-Huckel regime. From the results on energetics of partially formed shells we examine the stability of tethered (crystalline) and fluid shells towards rupture. We clearly delineate different regimes of stability towards rupture, where, for fluid shells, we also include the effects of bending elasticity of the shells. Our analysis shows how charging of the shell induces its instability towards rupture but also provides insight regarding growth of charged shells.

preprint2011arXiv

Energies and pressures in viruses: contribution of nonspecific electrostatic interactions

We summarize some aspects of electrostatic interactions in the context of viruses. A simplified but, within well defined limitations, reliable approach is used to derive expressions for electrostatic energies and the corresponding osmotic pressures in single-stranded RNA viruses and double-stranded DNA bacteriophages. The two types of viruses differ crucially in the spatial distribution of their genome charge which leads to essential differences in their free energies, depending on the capsid size and total charge in a quite different fashion. Differences in the free energies are trailed by the corresponding characteristics and variations in the osmotic pressure between the inside of the virus and the external bathing solution.

preprint2011arXiv

Many-body effects in van der Waals-Casimir interaction between graphene layers

Van der Waals-Casimir dispersion interactions between two apposed graphene layers, a graphene layer and a substrate, and in a multilamellar graphene system are analyzed within the framework of the Lifshitz theory. This formulation hinges on a known form of the dielectric response function of an undoped or doped graphene sheet, assumed to be of a random phase approximation form. In the geometry of two apposed layers the separation dependence of the van der Waals-Casimir interaction for both types of graphene sheets is determined and compared with some well known limiting cases. In a multilamellar array the many-body effects are quantified and shown to increase the magnitude of the van der Waals-Casimir interactions.

preprint2011arXiv

Ordering of anisotropic polarizable polymer chains on the full many-body level

We study the effect of dielectric anisotropy of polymers on their equilibrium ordering within mean-field theory but with a formalism that takes into account the full n-body nature of van der Waals forces. Dielectric anisotropy within polymers is to be expected as the electronic properties of the polymer will typically be different along the polymer than across its cross section. It is therefore physically intuitive that larger charge fluctuations can be induced along the chain than perpendicular to it. We show that this dielectric anisotropy leads to n-body interactions which can induce an isotropic--nematic transition. The two body and three body components of the full van der Waals interaction are extracted and it is shown how the two body term behaves like the phenomenological self-aligning-pairwise nematic interaction. At the three body interaction level we see that the nematic phase that is energetically favorable is discotic, however on the full n-body interaction level we find that the normal axial nematic phase is always the stable ordered phase. The n-body nature of our approach also shows that the key parameter driving the nematic-isotropic transition is the bare persistence length of the polymer chain.

preprint2011arXiv

Overscreening in 1D lattice Coulomb gas model of ionic liquids

Overscreening in the charge distribution of ionic liquids at electrified interfaces is shown to proceed from purely electrostatic and steric interactions in an exactly soluble one dimensional lattice Coulomb gas model. Being not a mean-field effect, our results suggest that even in higher dimensional systems the overscreening could be accounted for by a more accurate treatment of the basic lattice Coulomb gas model, that goes beyond the mean field level of approximation, without any additional interactions.

preprint2010arXiv

Confined polymer nematics: order and packing in a nematic drop

We investigate the tight packing of nematic polymers inside a confining hard sphere. We model the polymer {\sl via} the continuum Frank elastic free energy augmented by a simple density dependent part as well as by taking proper care of the connectivity of the polymer chains when compared with simple nematics. The free energy {\sl ansatz} is capable of describing an orientational ordering transition within the sample between an isotropic polymer solution and a polymer nematic phase. We solve the Euler-Lagrange equations numerically with the appropriate boundary conditions for the director and density field and investigate the orientation and density profile within a sphere. Two important parameters of the solution are the exact locations of the beginning and the end of the polymer chain. Pending on their spatial distribution and the actual size of the hard sphere enclosure we can get a plethora of various configurations of the chain exhibiting different defect geometry.

preprint2010arXiv

Counterion-Mediated Weak and Strong Coupling Electrostatic Interaction between Like-Charged Cylindrical Dielectrics

We examine the effective counterion-mediated electrostatic interaction between two like-charged dielectric cylinders immersed in a continuous dielectric medium containing neutralizing mobile counterions. We focus on the effects of image charges induced as a result of the dielectric mismatch between the cylindrical cores and the surrounding dielectric medium and investigate the counterion-mediated electrostatic interaction between the cylinders in both limits of weak and strong electrostatic couplings (corresponding, e.g., to systems with monovalent and multivalent counterions, respectively). The results are compared with extensive Monte-Carlo simulations exhibiting good agreement with the limiting weak and strong coupling results in their respective regime of validity.

preprint2010arXiv

Dressed Counterions: Strong Electrostatic Coupling in the Presence of Salt

We reformulate the theory of strong electrostatic coupling in order to describe an asymmetric electrolyte solution of monovalent salt ions and polyvalent counterions using field-theoretical techniques and Monte-Carlo simulations. The theory is based on an asymmetric treatment of the different components of the electrolyte solution. The weak coupling Debye-Huckel approach is used in order to describe the monovalent salt ions while a strong coupling approach is used to tackle the polyvalent counterions. This combined weak-strong coupling approach effectively leads to dressed interactions between polyvalent counterions and thus directly affects the correlation attraction mediated by polyvalent counterions between like-charged objects. The general theory is specifically applied to a system composed of two uniformly charged plane-parallel surfaces in the presence of salt and polyvalent counterions. In the strong coupling limit for polyvalent counterions the comparison with Monte-Carlo simulations shows good agreement for large enough values of the electrostatic coupling parameter. We delineate two limiting laws that in fact encompass all the Monte-Carlo data.

preprint2010arXiv

Exotic Electrostatics: Unusual Features of Electrostatic Interactions between Macroions

We present an overview of our understanding of electrostatic interactions between charged macromolecular surfaces mediated by mobile counter- and coions. The dichotomy between the weak and the strong coupling regimes is described in detail and the way they engender repulsive and attractive interactions between nominally equally charged macroions. We also introduce the concept of dressed counterions in the case of many-component Coulomb fluids that are partially weakly and partially strongly coupled to local electrostatic fields leading to non-monotonic interactions between equally charged macroions. The effect of quenched surface charge disorder on the counterion-mediated electrostatic interactions is analyzed within the same conceptual framework and shown to lead to unexpected and extraordinary electrostatic interactions between randomly charged surfaces with equal mean surface charge densities or even between effectively neutral macroion surfaces. As a result, these recent developments challenge some cherished notions of pop culture.

preprint2010arXiv

Non-monotoic fluctuation-induced interactions between dielectric slabs carrying charge disorder

We investigate the effect of monopolar charge disorder on the classical fluctuation-induced interactions between randomly charged net-neutral dielectric slabs and discuss various generalizations of recent results (A. Naji et al., Phys. Rev. Lett. 104, 060601 (2010)) to highly inhomogeneous dielectric systems with and without statistical disorder correlations. We shall focus on the specific case of two generally dissimilar plane-parallel slabs, which interact across vacuum or an arbitrary intervening dielectric medium. Monopolar charge disorder is considered to be present on the bounding surfaces and/or in the bulk of the slabs, may be in general quenched or annealed and may possess a finite lateral correlation length reflecting possible `patchiness' of the random charge distribution. In the case of quenched disorder, the bulk disorder is shown to give rise to an additive long-range contribution to the total force, which decays as the inverse distance between the slabs and may be attractive or repulsive depending on the dielectric constants of the slabs. We show that in the case of two dissimilar slabs the net effect due to the interplay between the disorder-induced and the pure van der Waals interactions can lead to a variety of unusual non-monotonic interaction profiles between the dielectric slabs. In particular, when the intervening medium has a larger dielectric constant than the two slabs, we find that the net interaction can become repulsive and exhibit a potential barrier, while the underlying van der Waals force is attractive. On the contrary, when the intervening medium has a dielectric constant in between that of the two slabs, the net interaction can become attractive and exhibit a free energy minimum, while the pure van der Waals force is repulsive. Therefore, the charge disorder, if present, can drastically alter the effective interaction between net-neutral objects.

preprint2010arXiv

Sample-to-sample fluctuations of electrostatic forces generated by quenched charge disorder

It has been recently shown that randomly charged surfaces can exhibit long range electrostatic interactions even when they are net neutral. These forces depend on the specific realization of charge disorder and thus exhibit sample to sample fluctuations about their mean value. We analyze the fluctuations of these forces in the parallel slab configuration and also in the sphere-plane geometry via the proximity force approximation. The fluctuations of the normal forces, that have a finite mean value, are computed exactly. Surprisingly, we also show that lateral forces are present, despite the fact that they have a zero mean, and that their fluctuations have the same scaling behavior as the normal force fluctuations. The measurement of these lateral force fluctuations could help to characterize the effects of charge disorder in experimental systems, leading to estimates of their magnitudes that are complementary to those given by normal force measurements.

preprint2009arXiv

Dispersion interactions in stratified anisotropic and optically active media at all separations

We propose a method to calculate dispersion interactions in a system composed of a one dimensional layering of finite thickness anisotropic and optically active slabs. The result is expressed within the algebra of 4x4 matrices and is demonstrated to be equivalent to the known limits of isotropic, nonretarded and uniaxial dispersion interactions. The method is also capable of handling dielectric media with smoothly varying anisotropy axes.

preprint2009arXiv

Fluctuation-Induced Interaction between Randomly Charged Dielectrics

Monopolar charge disorder effects are studied in the context of fluctuation-induced interactions between neutral dielectric slabs. It is shown that quenched bulk charge disorder gives rise to an additive contribution to the net interaction force which decays as the inverse distance between the slabs and may thus completely mask the standard Casimir--van der Waals force at large separations. By contrast, annealed (bulk or surface) charge disorder leads to a net interaction force whose large-distance behavior coincides with the universal Casimir force between perfect conductors, which scales as inverse cubic distance, and the dielectric properties enter only in subleading corrections.

preprint2009arXiv

The effects of dielectric disorder on van der Waals interactions in slab geometries

We study the thermal Casimir effect between two thick slabs composed of plane-parallel layers of random dielectric materials interacting across an intervening homogeneous dielectric. It is found that the effective interaction at long distances is self averaging and its value is given by a that between non-random media with the effective dielectric tensor of the corresponding random media. The behavior at short distances becomes random (sample dependent) and is dominated by the local values of the dielectric constants proximal to each other across the homogeneous slab. These results are extended to the regime of intermediate slab separations by using perturbation theory for weak disorder and also by extensive numerical simulations for a number of systems where the dielectric function has a log-normal distribution.

preprint2009arXiv

Thermodynamics of nano-spheres encapsulated in virus capsids

We investigate the thermodynamics of complexation of functionalized charged nano-spheres with viral proteins. The physics of this problem is governed by electrostatic interaction between the proteins and the nano-sphere cores (screened by salt ions), but also by configurational degrees of freedom of the charged protein N-tails. We approach the problem by constructing an appropriate complexation free energy functional. On the basis of both numerical and analytical studies of this functional we construct the phase diagram for the assembly which contains the information on the assembled structures that appear in the thermodynamical equilibrium, depending on the size and surface charge density of the nano-sphere cores. We show that both the nano-sphere core charge as well as its radius determine the size of the capsid that forms around the core.

preprint2006arXiv

Comparison of the DFT and field approaches to Van der Waals interactions in planparallel geometry

We establish a general equivalence between van der Waals interaction energies within the formalism of the non-local van der Waals functional of the density functional theory and within the formalism of the field approach based on the secular determinants of the electromagnetic field modes. We then compare the two methods explicitly in the case of a planparallel geometry with a continuously varying dielectric response function and show that their respective numerical implementations are not equivalent. This allows us to discuss the merits of the two approaches and possible advantages of either method in a simple model calculation.