Showing 1–50 of 188 results for author: Podgornik, R

Exact moments for a run and tumble particle in a harmonic trap with a finite tumble time

Authors: Aoran Sun, Fangfu Ye, Rudolf Podgornik

Abstract: We study the problem of a run and tumble particle in a harmonic trap, with a finite run and tumble time, by a direct integration of the equation of motion. An exact 1D steady state distribution, diagram laws and a programmable Volterra difference equation are derived to calculate any order of moments in any other dimension, both for steady state as well as the Laplace transform in time for the int… ▽ More We study the problem of a run and tumble particle in a harmonic trap, with a finite run and tumble time, by a direct integration of the equation of motion. An exact 1D steady state distribution, diagram laws and a programmable Volterra difference equation are derived to calculate any order of moments in any other dimension, both for steady state as well as the Laplace transform in time for the intermediate states. We also use the moments to infer the distribution by considering a Gaussian quadrature for the corresponding measure, and from the scaling law of high order moments. △ Less

Submitted 31 August, 2024; originally announced September 2024.

Comments: 12 pages 5 figures

Sequence disorder-induced first order phase transition in confined polyelectrolytes

Authors: V. Stepanyan, A. Badasyan, V. Morozov, Y. Mamasakhlisov, R. Podgornik

Abstract: We consider a statistical mechanical model of a generic flexible polyelectrolyte, comprised of identically charged monomers with long range electrostatic interactions, and short-range interactions quantified by a disorder field along the polymer contour sequence, which is randomly quenched. The free energy and the monomer density profile of the system for no electrolyte screening are calculated in… ▽ More We consider a statistical mechanical model of a generic flexible polyelectrolyte, comprised of identically charged monomers with long range electrostatic interactions, and short-range interactions quantified by a disorder field along the polymer contour sequence, which is randomly quenched. The free energy and the monomer density profile of the system for no electrolyte screening are calculated in the case of a system composed of two infinite plates with an intervening oppositely charged polyelectrolyte chain. We show that the effect of the contour sequence disorder mediated by short-range interactions leads to an enhanced localization of the polyelectrolyte chain and a first order phase transition at a critical value of the inter-surface spacing. This phase transition results in an abrupt change of the pressure from negative to positive values, effectively eliminating polyelectrolyte mediated bridging attraction. △ Less

Submitted 9 July, 2024; v1 submitted 6 July, 2024; originally announced July 2024.

Correlation decoupling of Casimir interaction in an electrolyte driven by external electric fields

Authors: Guangle Du, David S. Dean, Bing Miao, Rudolf Podgornik

Abstract: It has been established for a long time that the long range van der Waals or thermal Casimir interaction between two semi-infinite dielectrics separated by a distance $H$ is screened by an intervening electrolyte. Here we show how this interaction is modified when an electric field of strength $E$ is applied parallel to the dielectric boundaries, leading to a non-equilibrium steady state with a cu… ▽ More It has been established for a long time that the long range van der Waals or thermal Casimir interaction between two semi-infinite dielectrics separated by a distance $H$ is screened by an intervening electrolyte. Here we show how this interaction is modified when an electric field of strength $E$ is applied parallel to the dielectric boundaries, leading to a non-equilibrium steady state with a current. The presence of the field induces a long range thermal repulsive interaction, scaling just like the thermal Casimir interaction between dielectrics without the intervening electrolyte, {\em i.e.} as $1/H^3$. At small $E$ the effect is of order $E^2$ while at large fields it saturates to an $E$ independent value. We explain the results in terms of a decoupling mechanism between the charge density fluctuations of cations and anions at large applied fields. △ Less

Submitted 9 April, 2024; originally announced April 2024.

Comments: 6 pages, 2 figures

Fabry-Pérot nanocavities controlled by Casimir forces in electrolyte solutions

Authors: Lixin Ge, Kaipeng Liu, Ke Gong, Rudolf Podgornik

Abstract: We propose a design for tuning the resonant spectra of Fabry-Pérot nanocavities mediated by the Casimir force. The system involves a suspended gold nanoplate approaching to a dielectric-coated gold substrate in a univalent electrolyte solution. The gold nanoplate can be stably suspended due to the delicate balance between repulsive and attractive components of the Casimir forces. In an electrolyte… ▽ More We propose a design for tuning the resonant spectra of Fabry-Pérot nanocavities mediated by the Casimir force. The system involves a suspended gold nanoplate approaching to a dielectric-coated gold substrate in a univalent electrolyte solution. The gold nanoplate can be stably suspended due to the delicate balance between repulsive and attractive components of the Casimir forces. In an electrolyte solution, the presence of ionic-charge fluctuations can partially or totally screen the thermal $n$=0 Matsubara term, resulting in strongly modified interactions. As a result, the separation between the gold nanoplate and the substrate experiences a significant modulation in response to variations in salt concentration. Under proper conditions, we find that the modulation of the Casimir force would strongly shift the resonances of Fabry-Pérot nanocavities at the optical frequencies, when the Debye length of the electrolyte decreases from 1000 nm to 10 nm. Finally, the temperature dependence of the thermal Casimir force would provide an additional modulation of Fabry-Pérot nanocavity resonances for their eventual fine tuning. These results open up a promising venue for general tuning of the optical resonances with potential applications in re-configurable microfluidic nanophotonics. △ Less

Submitted 3 March, 2024; originally announced March 2024.

Comments: 10 pages, 5 figures

Spool-nematic ordering of dsDNA and dsRNA under confinement

Authors: James Farrell, Jure Dobnikar, Rudolf Podgornik, Tine Curk

Abstract: The ability of double-stranded DNA or RNA to locally melt and form kinks leads to strong non-linear elasticity effects that qualitatively affect their packing in confined spaces. Using analytical theory and numerical simulation we show that kink formation entails a mixed spool-nematic ordering of double-stranded DNA or RNA in spherical capsids, consisting of an outer spool domain and an inner, twi… ▽ More The ability of double-stranded DNA or RNA to locally melt and form kinks leads to strong non-linear elasticity effects that qualitatively affect their packing in confined spaces. Using analytical theory and numerical simulation we show that kink formation entails a mixed spool-nematic ordering of double-stranded DNA or RNA in spherical capsids, consisting of an outer spool domain and an inner, twisted nematic domain. These findings explain the experimentally observed nematic domains in viral capsids and imply that non-linear elasticity must be considered to predict the configurations and dynamics of double-stranded genomes in viruses, bacterial nucleoids or gene-delivery vehicles. The non-linear elastic theory suggests that spool-nematic ordering is a general feature of strongly confined kinkable polymers. △ Less

Submitted 10 October, 2024; v1 submitted 16 December, 2023; originally announced December 2023.

Comments: 6 pages

Journal ref: Phys. Rev. Lett. 133, 148101 (2024)

On the nature of screening in charge regulated macroion solutions

Authors: Sunita Kumari, Rudolf Podgornik

Abstract: We present a derivation of the screening length for a solution containing a charge regulated macroion, e.g. protein, with its counterions. We show that it can be obtained directly from the second derivatives of the total free energy by taking recourse to the ``uncertainty relation" of the Legendre transform which connects the Hessians or the local curvatures of the free energy as a function of den… ▽ More We present a derivation of the screening length for a solution containing a charge regulated macroion, e.g. protein, with its counterions. We show that it can be obtained directly from the second derivatives of the total free energy by taking recourse to the ``uncertainty relation" of the Legendre transform which connects the Hessians or the local curvatures of the free energy as a function of density and its Legendre transform, i.e., osmotic pressure, as a function of chemical potentials. Based on the Fowler-Guggenheim-Frumkin model of charge regulation we then analyze the ``screening resonance" and the ``overscreening" of the screening properties of the charge-regulated macroion solution. △ Less

Submitted 17 November, 2023; originally announced November 2023.

Comments: 11 Pages, 5 Figures

Stability of Taylor-Couette Flow with Odd Viscosity

Authors: Guangle Du, Rudolf Podgornik

Abstract: Odd viscosity can emerge in 3D hydrodynamics when the time reversal symmetry is broken and anisotropy is introduced. Its ramifications on the stability of the prototypical Taylor-Couette flow in curved geometries have remained unexplored. Here, we investigate the effects of odd viscosity on the stability of Taylor-Couette flow under axisymmetric perturbations both analytically and numerically, der… ▽ More Odd viscosity can emerge in 3D hydrodynamics when the time reversal symmetry is broken and anisotropy is introduced. Its ramifications on the stability of the prototypical Taylor-Couette flow in curved geometries have remained unexplored. Here, we investigate the effects of odd viscosity on the stability of Taylor-Couette flow under axisymmetric perturbations both analytically and numerically, deriving analytically the critical Taylor number for different odd viscosities in the narrow gap case as well as fully numerically implementing the wide gap case. We find that the odd viscosity modifies the vortex pattern by creating secondary vortices, and exerts an intriguing ``lever'' effect in the stability diagram, completely suppressing the instability of Taylor-Couette flow under axisymmetric perturbations when the odd viscosity is large, irrespective of its sign. Our findings highlight the role of odd viscosity for the rich flow patterns of the Taylor-Couette geometry and provide guidance for viscometer experiments when odd viscosity is present. △ Less

Submitted 21 September, 2023; v1 submitted 18 September, 2023; originally announced September 2023.

Electrostatic interactions between charge regulated spherical macroions

Authors: Hu Ruixuan, Arghya Majee, Jure Dobnikar, Rudolf Podgornik

Abstract: We study the interaction between two charge regulating spherical macroions with dielectric interior and dissociable surface groups immersed in a monovalent electrolyte solution. The charge dissociation is modelled via the Frumkin-Fowler-Guggenheim isotherm, which allows for multiple adsorption equilibrium states. The interactions are derived from the solutions of the mean-field Poisson-Boltzmann t… ▽ More We study the interaction between two charge regulating spherical macroions with dielectric interior and dissociable surface groups immersed in a monovalent electrolyte solution. The charge dissociation is modelled via the Frumkin-Fowler-Guggenheim isotherm, which allows for multiple adsorption equilibrium states. The interactions are derived from the solutions of the mean-field Poisson-Boltzmann type theory with charge regulation boundary conditions. For a range of conditions we find symmetry breaking transitions from symmetric to asymmetric charge distribution exhibiting annealed charge patchiness, which results in like-charge attraction even in a univalent electrolyte -- thus fundamentally modifying the nature of electrostatic interactions in charge-stabilized colloidal suspensions. △ Less

Submitted 7 December, 2023; v1 submitted 6 August, 2023; originally announced August 2023.

Comments: 18 pages, 8 figures

Journal ref: Eur. Phys. J. E (2023)46:115

Interplay between finite thickness and chirality effects on the Casimir-Lifshitz torque with nematic cholesteric liquid crystals

Authors: Wijnand Broer, Rudolf Podgornik

Abstract: We theoretically investigate the combined effects of the chirality and the finite total thickness of nematic cholesteric liquid crystals on the Casimir-Lifshitz torque. We find that, the larger the thickness, the more sinusoidal the angular dependence of the torque becomes. We use a Fourier decomposition to quantify this result. The general direction of the torque depends on whether the configurat… ▽ More We theoretically investigate the combined effects of the chirality and the finite total thickness of nematic cholesteric liquid crystals on the Casimir-Lifshitz torque. We find that, the larger the thickness, the more sinusoidal the angular dependence of the torque becomes. We use a Fourier decomposition to quantify this result. The general direction of the torque depends on whether the configuration of two cholesterics is heterochiral or homochiral. △ Less

Submitted 22 May, 2023; originally announced May 2023.

Comments: 13 pages, 5 figures

Curvature effects in interfacial acidity of amphiphilic vesicles

Authors: Petch Khunpetch, Arghya Majee, Hu Ruixuan, Rudolf Podgornik

Abstract: We analyze the changes in the vicinal acidity (pH) at a spherical amphiphilic membrane. The membrane is assumed to contain solvent accessible, embedded, dissociable, charge regulated moieties. Basing our approach on the linear Debye-Hückel as well as the non-linear Poisson-Boltzmann theory, together with the general Frumkin-Fowler-Guggenheim adsorption isotherm model of the charge regulation proce… ▽ More We analyze the changes in the vicinal acidity (pH) at a spherical amphiphilic membrane. The membrane is assumed to contain solvent accessible, embedded, dissociable, charge regulated moieties. Basing our approach on the linear Debye-Hückel as well as the non-linear Poisson-Boltzmann theory, together with the general Frumkin-Fowler-Guggenheim adsorption isotherm model of the charge regulation process, we analyse and review the dependence of the local pH on the position, as well as bulk electrolyte concentration, bulk pH and curvature of the amphiphilic single membrane vesicle. With appropriately chosen adsorption parameters of the charge regulation model, we find a good agreement with available experimental data. △ Less

Submitted 18 May, 2023; originally announced May 2023.

Comments: 17 pages, 11 figures

Continuum theory of electrostatic-elastic coupling interactions in colloidal crystals

Authors: Hao Wu, Rudolf Podgornik

Abstract: Electrostatic-elastic coupling in colloidal crystals, composed of a mobile Coulomb gas permeating a fixed background crystalline lattice of charged colloids, is studied on the continuum level in order to analyze the lattice-mediated interactions between mobile charges. The linearized, Debye-H\" uckel-like mean-field equations allow for an analytic solution incorporating the minimal coupling betwee… ▽ More Electrostatic-elastic coupling in colloidal crystals, composed of a mobile Coulomb gas permeating a fixed background crystalline lattice of charged colloids, is studied on the continuum level in order to analyze the lattice-mediated interactions between mobile charges. The linearized, Debye-H\" uckel-like mean-field equations allow for an analytic solution incorporating the minimal coupling between electrostatic and displacement fields leading to an additional effective attractive interaction between mobile charges. For small screening lengths, the interactions between mobile like charges exhibit a coupling-strength-dependent behavior, with lattice mediated effective interaction ranging from weak to strong attraction, while for large screening lengths the lattice mediated interaction remains purely attractive. By analyzing Gaussian fluctuations around the mean-field solution we identify its region of stability in terms of the electrostatic-elastic screening length. A detailed continuum theory incorporating the standard lattice elasticity and electrostatics of mobile charges augmented by the electrostatic-elastic coupling terms can serve as a baseline for more detailed microscopic models. △ Less

Submitted 16 April, 2023; originally announced April 2023.

Comments: 11 pages, 4 figures

Charge Regulation of Polyelectrolyte Gels: Swelling Transition

Authors: Bin Zheng, Yael Avni, David Andelman, Rudolf Podgornik

Abstract: We study the effects of charge-regulated acid/base equilibrium on the swelling of polyelectrolyte gels, by considering a combination of the Poisson-Boltzmann theory and a two-site charge-regulation model based on the Langmuir adsorption isotherm. By exploring the volume change as a function of salt concentration for both nano-gels and micro-gels, we identify conditions where the gel volume exhibit… ▽ More We study the effects of charge-regulated acid/base equilibrium on the swelling of polyelectrolyte gels, by considering a combination of the Poisson-Boltzmann theory and a two-site charge-regulation model based on the Langmuir adsorption isotherm. By exploring the volume change as a function of salt concentration for both nano-gels and micro-gels, we identify conditions where the gel volume exhibits a discontinuous swelling transition. This transition is driven exclusively by the charge-regulation mechanism and is characterized by a closed-loop phase diagram. Our predictions can be tested experimentally for polypeptide gels. △ Less

Submitted 4 April, 2023; originally announced April 2023.

Journal ref: Macromolecules 2023, 56, 5217

Viral RNA as a branched polymer

Authors: Domen Vaupotič, Angelo Rosa, Rudolf Podgornik, Luca Tubiana, Anže Božič

Abstract: Myriad viruses use positive-strand RNA molecules as their genomes. Far from being only a repository of genetic material, viral RNA performs numerous other functions mediated by its physical structure and chemical properties. In this chapter, we focus on its structure and discuss how long RNA molecules can be treated as branched polymers through planar graphs. We describe the major results that can… ▽ More Myriad viruses use positive-strand RNA molecules as their genomes. Far from being only a repository of genetic material, viral RNA performs numerous other functions mediated by its physical structure and chemical properties. In this chapter, we focus on its structure and discuss how long RNA molecules can be treated as branched polymers through planar graphs. We describe the major results that can be obtained by this approach, in particular the observation that viral RNA genomes have a characteristic compactness that sets them aside from similar random RNAs. We also discuss how different parameters used in the current RNA folding software influence the resulting structures and how they can be related to experimentally observable quantities. Finally, we show how the connection to branched polymers can be extended to take advantage of known results from polymer physics and can be further moulded to include additional interactions, such as excluded volume or electrostatics. △ Less

Submitted 27 March, 2023; v1 submitted 1 December, 2022; originally announced December 2022.

Comments: 24 pages, 9 figures

Journal ref: Chapter 1 in: Comas-Garcia, M., Rosales-Mendoza, S. (eds) Physical Virology. Springer Series in Biophysics, vol 24. Springer, Cham (2023)

Trapping Instability of an Active Particle in Steering Potential Fields

Authors: Guangle Du, Fangfu Ye, Rudolf Podgornik

Abstract: A particle driven by active self-propulsion can be subject to inhomogeneous potential fields, steering its orientation and leading to confinement and eventual trapping. Analytical treatment of capture and/or release dynamics for general steering potentials presents a challenge due to its coupling between external potential fields and intrinsic active noise. By using the projection operator method… ▽ More A particle driven by active self-propulsion can be subject to inhomogeneous potential fields, steering its orientation and leading to confinement and eventual trapping. Analytical treatment of capture and/or release dynamics for general steering potentials presents a challenge due to its coupling between external potential fields and intrinsic active noise. By using the projection operator method we obtain the coarse-grained Dynkin equations with orientation integrated out in the large fluctuations limit, and derive explicit analytical solutions for the mean first passage time in radially symmetric point source trapping potentials. We analyze the ensuing trapping instabilities related to a critical value of the steering potential strength below which the particle either cannot be lured into the trap, or above which it is unable to leave the trap after being lured into it. △ Less

Submitted 21 June, 2022; originally announced June 2022.

A comprehensive continuum theory of structured liquids

Authors: Ralf Blossey, Rudolf Podgornik

Abstract: We develop a comprehensive continuum model capable of treating both electrostatic and structural interactions in liquid dielectrics. Starting from a two-order parameter description in terms of charge density and polarization, we derive a field-theoretic model generalizing previous theories. Our theory explicitly includes electrostatic and structural interactions in the bulk of the liquid and allow… ▽ More We develop a comprehensive continuum model capable of treating both electrostatic and structural interactions in liquid dielectrics. Starting from a two-order parameter description in terms of charge density and polarization, we derive a field-theoretic model generalizing previous theories. Our theory explicitly includes electrostatic and structural interactions in the bulk of the liquid and allows for polarization charges within a Drude model. In particular, we develop a detailed description of the boundary conditions which include the charge regulation mechanism and surface polarization. The general features for solving the saddle-point equations of our model are elucidated and future applications to predict and validate experimental results are outlined. △ Less

Submitted 22 January, 2023; v1 submitted 7 June, 2022; originally announced June 2022.

Comments: 19 pages; accepted for publication in J. Phys. A

Continuum theories of structured dielectrics

Authors: Ralf Blossey, Rudolf Podgornik

Abstract: 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 approac… ▽ More 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. △ Less

Submitted 3 June, 2022; originally announced June 2022.

Comments: 7 pages, 2 figures

Field theory of structured liquid dielectrics

Authors: Ralf Blossey, Rudolf Podgornik

Abstract: 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 treat… ▽ More 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. △ Less

Submitted 31 January, 2022; originally announced January 2022.

Comments: 13 pages, 7 figures

Curvature effects in charge-regulated lipid bilayers

Authors: Petch Khunpetch, Arghya Majee, Rudolf Podgornik

Abstract: 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 cha… ▽ More 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. △ Less

Submitted 15 February, 2022; v1 submitted 13 January, 2022; originally announced January 2022.

Comments: 25 pages, 6 figures

Journal ref: Soft Matter 18, 2597 (2022)

Qualitative chirality effects on the Casimir-Lifshitz torque with liquid crystals

Authors: Wijnand Broer, Bing-Sui Lu, Rudolf Podgornik

Abstract: We model a cholesteric liquid crystal as a planar uniaxial multilayer system, where the orientation of each layer differs slightly from that of the adjacent one. This allows us to analytically simplify the otherwise acutely complicated calculation of the Casimir-Lifshitz torque. Numerical results differ appreciably from the case of nematic liquid crystals, which can be treated like bloc birefringe… ▽ More We model a cholesteric liquid crystal as a planar uniaxial multilayer system, where the orientation of each layer differs slightly from that of the adjacent one. This allows us to analytically simplify the otherwise acutely complicated calculation of the Casimir-Lifshitz torque. Numerical results differ appreciably from the case of nematic liquid crystals, which can be treated like bloc birefringent media. In particular, we find that the torque deviates considerably from its usual sinusoidal behavior as a function of the misalignment angle. In the case of a birefringent crystal faced with a cholesteric liquid one, the Casimir-Lifshitz torque decreases more slowly as a function of distance than in the nematic case. In the case of two cholesteric liquid crystals, either in the homochiral or in the heterochiral configuration, the angular dependence changes qualitatively as a function of distance. In all considered cases, finite pitch length effects are most pronounced at distances of about 10 nm. △ Less

Submitted 30 June, 2021; originally announced June 2021.

Journal ref: Phys. Rev. Research 3, 033238 (2021)

Model of metameric locomotion in active directional filaments

Authors: Guangle Du, Sunita Kumari, Fangfu Ye, Rudolf Podgornik

Abstract: Locomotion in segmented animals, such as annelids and myriapods (centipedes and millipedes), is generated by a coordinated movement known as metameric locomotion, which can be also implemented in robots designed to perform specific tasks. We introduce a theoretical model, based on an active directional motion of the head segment and a passive trailing of the rest of the body segments, in order to… ▽ More Locomotion in segmented animals, such as annelids and myriapods (centipedes and millipedes), is generated by a coordinated movement known as metameric locomotion, which can be also implemented in robots designed to perform specific tasks. We introduce a theoretical model, based on an active directional motion of the head segment and a passive trailing of the rest of the body segments, in order to formalize and study the metameric locomotion. The model is specifically formulated as a steered Ornstein-Uhlenbeck curvature process, preserving the continuity of the curvature along the whole body filament, and thus supersedes the simple active Brownian model, which would be inapplicable in this case. We obtain the probability density by analytically solving the Fokker-Planck equation pertinent to the model. We also calculate explicitly the correlators, such as the mean-square orientational fluctuations, the orientational correlation function and the mean-square separation between the head and tail segments, both analytically either via the Fokker-Planck equation or directly by either solving analytically or implementing it numerically from the Langevin equations. The analytical and numerical results coincide. Our theoretical model can help understand the locomotion of metameric animals and instruct the design of metameric robots. △ Less

Submitted 18 June, 2021; originally announced June 2021.

Comments: 11 pages, 4 figures

Contribution of dipolar bridging to phospholipid membrane interactions: a mean-field analysis

Authors: Sahin Buyukdagli, Rudolf Podgornik

Abstract: We develop a model of interacting zwitterionic membranes with rotating surface dipoles immersed in a monovalent salt, and implement it in a field theoretic formalism. In the mean-field regime of monovalent salt, the electrostatic forces between the membranes are characterized by a non-uniform trend: at large membrane separations, the interfacial dipoles on the opposing sides behave as like-charge… ▽ More We develop a model of interacting zwitterionic membranes with rotating surface dipoles immersed in a monovalent salt, and implement it in a field theoretic formalism. In the mean-field regime of monovalent salt, the electrostatic forces between the membranes are characterized by a non-uniform trend: at large membrane separations, the interfacial dipoles on the opposing sides behave as like-charge cations and give rise to repulsive membrane interactions; at short membrane separations, the anionic field induced by the dipolar phosphate groups sets the behavior in the intermembrane region. The attraction of the cationic nitrogens in the dipolar lipid headgroups leads to the adhesion of the membrane surfaces via dipolar bridging. The underlying competition between the opposing field components of the individual dipolar charges leads to the non-uniform salt ion affinity of the zwitterionic membrane with respect to the separation distance; large inter-membrane separations imply anionic excess while small, nanometer size separations, favor cationic excess. This complex ionic selectivity of zwitterionic membranes may have relevant repercussions on nanofiltration and nanofluidic transport techniques. △ Less

Submitted 18 May, 2021; v1 submitted 10 April, 2021; originally announced April 2021.

Journal ref: J. Chem Phys. 154, 224902 (2021)

Phase Separation of Polyelectrolytes: The Effect of Charge Regulation

Authors: Bin Zheng, Yael Avni, David Andelman, Rudolf Podgornik

Abstract: Complex coacervation, known as the liquid-liquid phase separation of solutions with oppositely charged polyelectrolytes, has attracted substantial interest in recent years. We study the effect of the charge regulation (CR) mechanism on the complex coacervation by including short-range interactions between the charged sites on the polymer chains as well as an association-dissociation energy paramet… ▽ More Complex coacervation, known as the liquid-liquid phase separation of solutions with oppositely charged polyelectrolytes, has attracted substantial interest in recent years. We study the effect of the charge regulation (CR) mechanism on the complex coacervation by including short-range interactions between the charged sites on the polymer chains as well as an association-dissociation energy parameter in the CR mechanism. We investigate the phase diagrams of two CR models: (i) the hopping CR model (HCR) and (ii) the asymmetric CR model (ACR). It is shown that during the phase separation that the polymers in the condensed phase are more charged than those in the dilute phase, in accordance with Le Chatelier's principle. In addition, secondary CR effects also influence the change in the volume fraction of the two phases. The latter can cause the charge difference between the two phases to change nonmonotonically as a function of the CR parameters. △ Less

Submitted 7 July, 2021; v1 submitted 5 March, 2021; originally announced March 2021.

Journal ref: J. Phys. Chem. B 125, 7863 (2021)

Embarras de richesses in non-DLVO colloid interactions

Authors: Rudolf Podgornik, David Andelman

Abstract: 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 a… ▽ More 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. △ Less

Submitted 5 January, 2021; v1 submitted 1 January, 2021; originally announced January 2021.

Comments: 4 pages

Journal ref: J. Club Cond. Matter September 02 (2020)

On the nature of screening in Voorn-Overbeek type theories

Authors: Sunita Kumari, Shikha Dwivedi, Rudolf Podgornik

Abstract: By using the recently formulated Legendre transform approach to the thermodynamics of charged systems, we explore the general form of the screening length in the Voorn-Overbeek type theories, that remains valid also in the cases where the entropy is not given by the ideal gas form as in the Debye-Huckel theory. The screening length consistent with the non-electrostatic terms in the free energy Ans… ▽ More By using the recently formulated Legendre transform approach to the thermodynamics of charged systems, we explore the general form of the screening length in the Voorn-Overbeek type theories, that remains valid also in the cases where the entropy is not given by the ideal gas form as in the Debye-Huckel theory. The screening length consistent with the non-electrostatic terms in the free energy Ansatz for the Flory-Huggins and Voorn-Overbeek type theories, derived from the local curvature properties of the Legendre transform, has distinctly different behavior than the often invoked standard Debye screening length, though it reduces to it in some limiting cases. △ Less

Submitted 8 March, 2022; v1 submitted 24 December, 2020; originally announced December 2020.

Comments: 8 pages, 4 figures

Electrostatic interaction between SARS-CoV-2 virus and charged electret fibre

Authors: Leili Javidpour, Anze Bozic, Ali Naji, Rudolf Podgornik

Abstract: While almost any kind of face mask offers some protection against particles and pathogens of different sizes, the most efficient ones make use of a layered structure where one or more layers are electrically charged. This electret layer is essential to efficient filtration of difficult-to-capture small particles, yet the exact nature of electrostatic capture with respect to both the charge on the… ▽ More While almost any kind of face mask offers some protection against particles and pathogens of different sizes, the most efficient ones make use of a layered structure where one or more layers are electrically charged. This electret layer is essential to efficient filtration of difficult-to-capture small particles, yet the exact nature of electrostatic capture with respect to both the charge on the particles and the electret fibres as well as the effect of immediate environment remains unclear. Here, we explore in detail the electrostatic interaction between the surface of a single charged electret fibre and a model of SARS-CoV-2 virus. Using Poisson-Boltzmann electrostatics coupled to a detailed spike protein charge regulation model, we show how pH and salt concentration drastically change both the scale and the sign of the interaction. Furthermore, the configuration of the few spike proteins closest to the electret fibre turns out to be as important for the strength of the interaction as their total number on the virus, a direct consequence of spike protein charge regulation. The results of our work elucidate the details of virus electrostatics and contribute to the general understanding of efficient virus filtration mechanisms. △ Less

Submitted 7 April, 2021; v1 submitted 13 December, 2020; originally announced December 2020.

Comments: 8 pages, 5 figures; 3 pages and 2 figures of SI

Journal ref: Soft Matter 17, 4296-4303 (2021)

Theory of inhomogeneous rod-like Coulomb fluids

Authors: Rudolf Podgornik

Abstract: 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 solu… ▽ More 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. △ Less

Submitted 27 January, 2021; v1 submitted 11 December, 2020; originally announced December 2020.

Comments: 17 pages, 4 figures

Charge symmetry broken complex coacervation

Authors: Arghya Majee, Markus Bier, Ralf Blossey, Rudolf Podgornik

Abstract: Liquid-liquid phase separation has emerged as one of the important paradigms in the chemical physics as well as biophysics of charged macromolecular systems. We elucidate an equilibrium phase separation mechanism based on charge regulation, i.e., protonation-deprotonation equilibria controlled by pH, in an idealized macroion system which can serve as a proxy for simple coacervation. First, a low-d… ▽ More Liquid-liquid phase separation has emerged as one of the important paradigms in the chemical physics as well as biophysics of charged macromolecular systems. We elucidate an equilibrium phase separation mechanism based on charge regulation, i.e., protonation-deprotonation equilibria controlled by pH, in an idealized macroion system which can serve as a proxy for simple coacervation. First, a low-density density-functional calculation reveals the dominance of two-particle configurations coupled by ion adsorption on neighboring macroions. Then a binary cell model, solved on the Debye-Hückel as well as the full nonlinear Poisson-Boltzmann level, unveils the charge-symmetry breaking as inducing the phase separation between low- and high-density phases as a function of pH. These results can be identified as a charge symmetry broken complex coacervation between chemically identical macroions. △ Less

Submitted 24 December, 2020; v1 submitted 1 October, 2020; originally announced October 2020.

Comments: 11 pages, 7 figures

Journal ref: Phys. Rev. Research 2, 043417 (2020)

Ordering of adsorbed rigid rods mediated by the Boussinesq interaction on a soft substrate

Authors: Sunita Kumari, Fangfu Ye, Rudolf Podgornik

Abstract: Orientational ordering driven by mechanical distortion of soft substrates plays a major role in material transformation processes such as elastocapillarity and surface anchoring. We present a theoretical model of the orientational response of anisotropic rods deposited onto a surface of a soft, elastic substrate of finite thickness. We show that anisotropic rods exhibit a continuous isotropic-nema… ▽ More Orientational ordering driven by mechanical distortion of soft substrates plays a major role in material transformation processes such as elastocapillarity and surface anchoring. We present a theoretical model of the orientational response of anisotropic rods deposited onto a surface of a soft, elastic substrate of finite thickness. We show that anisotropic rods exhibit a continuous isotropic-nematic phase transition, driven by orientational interactions between surface deposited rods. This interaction is mediated by the deformation of the underlying elastic substrate, and is quantified by the Boussinesq solution adapted to the case of slender, surface deposited rods. From the microscopic rod-rod interactions we derive the appropriate Maier-Saupe mean-field description, which includes the Boussinesq elastic free-energy contribution due to the substrate elasticity, derive the conditions for the existence of a continuous orientational ordering transition and discuss the implication of results in the soft (bio) systems context. △ Less

Submitted 29 September, 2020; v1 submitted 19 July, 2020; originally announced July 2020.

Comments: 22 pages, 9 figures

Journal ref: J. Chem. Phys. (2020)

Interactions between Zwitterionic Membranes in Complex Electrolytes

Authors: Sahin Buyukdagli, Rudolf Podgornik

Abstract: 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 dependenc… ▽ More 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. △ Less

Submitted 26 May, 2020; originally announced May 2020.

Journal ref: Phys. Rev. E 102, 012806 (2020)

Critical behavior of charge-regulated macro-ions

Authors: Yael Avni, Rudolf Podgornik, David Andelman

Abstract: Based on a collective description of electrolytes composed of charge-regulated macro-ions and simple salt ions, we analyze their equilibrium charge state in the bulk and their behavior in the vicinity of an external electrified surface. The mean-field formulation of mobile macro-ions in an electrolyte bathing solution is extended to include interactions between association/dissociation sites. We d… ▽ More Based on a collective description of electrolytes composed of charge-regulated macro-ions and simple salt ions, we analyze their equilibrium charge state in the bulk and their behavior in the vicinity of an external electrified surface. The mean-field formulation of mobile macro-ions in an electrolyte bathing solution is extended to include interactions between association/dissociation sites. We demonstrate that above a critical concentration of salt, and similar to the critical micelle concentration, a non-trivial distribution of charge states sets in. Such a charge state can eventually lead to a liquid-liquid phase separation based on charge regulation. △ Less

Submitted 9 July, 2020; v1 submitted 23 April, 2020; originally announced April 2020.

Comments: 12 pages, 7 figures

Journal ref: J. Chem. Phys. 153, 024901 (2020)

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

Authors: David S. Dean, Bing Miao, Rudolf Podgornik

Abstract: 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 condit… ▽ More 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. △ Less

Submitted 22 February, 2020; originally announced February 2020.

Comments: 18 pages, 6 figures

Surface alignment disorder and pseudo-Casimir forces in smectic-A liquid crystalline films

Authors: Fahimeh Karimi Pour Haddadan, Ali Naji, Rudolf Podgornik

Abstract: Random (disordered) components in the surface anchoring of the smectic-A liquid crystalline film in general modify the thermal pseudo-Casimir interaction. Anchoring disorder of the quenched type is in general decoupled from the thermal pseudo-Casimir force and gives rise to an additional disorder-generated interaction, in distinction to the annealed disorder, whose effect on the pseudo-Casimir for… ▽ More Random (disordered) components in the surface anchoring of the smectic-A liquid crystalline film in general modify the thermal pseudo-Casimir interaction. Anchoring disorder of the quenched type is in general decoupled from the thermal pseudo-Casimir force and gives rise to an additional disorder-generated interaction, in distinction to the annealed disorder, whose effect on the pseudo-Casimir force is non-additive. We consider the effects of the surface anchoring disorder by assuming that one of the substrates of the film is contaminated by a disorder source, resulting in a Gaussian-weighted distribution of the preferred molecular anchoring orientation (easy axes) on that substrate, having a finite mean and variance or, more generally, a homogeneous in-plane, two-point correlation function. We show that the presence of disorder, either of the quenched or annealed type, leads to a significant reduction in the magnitude of the net thermal fluctuation force between the confining substrates of the film. In the quenched case this is a direct consequence of an additive free energy dependent on the variance of the disorder, while in the annealed case, the suppression of the interaction force can be understood based on a disorder-renormalized, effective anchoring strength. △ Less

Submitted 17 September, 2019; originally announced September 2019.

Comments: 8 pages, 4 figures

Journal ref: J. Phys.: Condens. Matter 32, 325103 (2020)

Charge regulation radically modifies electrostatics in membrane stacks

Authors: Arghya Majee, Markus Bier, Ralf Blossey, Rudolf Podgornik

Abstract: Motivated by biological membrane-containing organelles in plants and photosynthetic bacteria, we study charge regulation in a model membrane stack. Considering (de)protonation as the simplest mechanism of charge equilibration between the membranes and with the bathing environment, we uncover a symmetry-broken charge state in the stack with a quasiperiodic effective charge sequence. In the case of… ▽ More Motivated by biological membrane-containing organelles in plants and photosynthetic bacteria, we study charge regulation in a model membrane stack. Considering (de)protonation as the simplest mechanism of charge equilibration between the membranes and with the bathing environment, we uncover a symmetry-broken charge state in the stack with a quasiperiodic effective charge sequence. In the case of a monovalent bathing salt solution our model predicts complex, inhomogeneous charge equilibria depending on the strength of the (de)protonation reaction, salt concentration, intermembrane separation, and their number in the stack. Our results shed light on the basic reorganization mechanism of thylakoid membrane stacks. △ Less

Submitted 6 November, 2019; v1 submitted 31 August, 2019; originally announced September 2019.

Comments: 16 pages, 7 figures

Journal ref: Phys. Rev. E 100, 050601 (2019)

Spontaneous Domain Formation in Spherically-Confined Elastic Filaments

Authors: Tine Curk, James Daniel Farrell, Jure Dobnikar, Rudolf Podgornik

Abstract: Although the free energy of a genome packing into a virus is dominated by DNA-DNA interactions, ordering of the DNA inside the capsid is elasticity-driven, suggesting general solutions with DNA organized into spool-like domains. Using analytical calculations and computer simulations of a long elastic filament confined to a spherical container, we show that the ground state is not a single spool as… ▽ More Although the free energy of a genome packing into a virus is dominated by DNA-DNA interactions, ordering of the DNA inside the capsid is elasticity-driven, suggesting general solutions with DNA organized into spool-like domains. Using analytical calculations and computer simulations of a long elastic filament confined to a spherical container, we show that the ground state is not a single spool as assumed hitherto, but an ordering mosaic of multiple homogeneously-ordered domains. At low densities, we observe concentric spools, while at higher densities, other morphologies emerge, which resemble topological links. We discuss our results in the context of metallic wires, viral DNA, and flexible polymers. △ Less

Submitted 19 March, 2024; v1 submitted 30 June, 2019; originally announced July 2019.

Comments: 5 pages

Journal ref: PRL 2019

Path integrals for higher derivative actions

Authors: David S. Dean, Bing Miao, Rudi Podgornik

Abstract: We consider Euclidean path integrals with higher derivative actions, including those that depend quadratically on acceleration, velocity and position. Such path integrals arise naturally in the study of stiff polymers, membranes with bending rigidity as well as a number of models for electrolytes. The approach used is based on the relation between quadratic path integrals and Gaussian fields and w… ▽ More We consider Euclidean path integrals with higher derivative actions, including those that depend quadratically on acceleration, velocity and position. Such path integrals arise naturally in the study of stiff polymers, membranes with bending rigidity as well as a number of models for electrolytes. The approach used is based on the relation between quadratic path integrals and Gaussian fields and we also show how it can be extended to the evaluation of even higher order path integrals. △ Less

Submitted 1 November, 2019; v1 submitted 20 June, 2019; originally announced June 2019.

Comments: Revtex 12 pages

Like-charge polymer-membrane complexation mediated by multivalent cations: one-loop-dressed strong coupling theory

Authors: Sahin Buyukdagli, Rudolf Podgornik

Abstract: We probe the electrostatic mechanism driving adsorption of polyelectrolytes onto like-charged membranes upon the addition of tri- and tetravalent counterions to a bathing monovalent salt solution. We develop a one-loop-dressed strong coupling theory that treats the monovalent salt at the electrostatic one-loop level and the multivalent counterions within a strong-coupling approach. It is shown tha… ▽ More We probe the electrostatic mechanism driving adsorption of polyelectrolytes onto like-charged membranes upon the addition of tri- and tetravalent counterions to a bathing monovalent salt solution. We develop a one-loop-dressed strong coupling theory that treats the monovalent salt at the electrostatic one-loop level and the multivalent counterions within a strong-coupling approach. It is shown that the adhesive force of the multivalent counterions mediating the like-charge adsorption arises from their strong condensation at the charged membrane. The resulting interfacial counterion excess locally maximizes the screening ability of the electrolyte and minimizes the electrostatic polymer grand potential. This translates into an attractive force that pulls the polymer to the similarly charged membrane. We show that the high counterion valency enables this adsorption transition even at weakly charged membranes. Additionally, strongly charged membranes give rise to salt-induced correlations and intensify the interfacial multivalent counterion condensation, strenghtening the complexation of the polymer with the like-charged membrane, as well as triggering the orientational transition of the molecule prior to its adsorption. Finally, our theory provides two additional key features as evidenced by previous adsorption experiments: first, the critical counterion concentration for polymer adsorption decreases with the rise of the counterion valency, and second, the addition of monovalent salt enhances the screening of the membrane charges and suppresses salt correlations. This weakens the interfacial multivalent counterion condensation and results in the desorption of the polymer from the substrate. △ Less

Submitted 10 May, 2019; originally announced May 2019.

Journal ref: J. Chem Phys. 151, 094902 (2019)

Orientational transition and complexation of DNA with anionic membranes: weak and intermediate electrostatic coupling

Authors: Sahin Buyukdagli, Rudolf Podgornik

Abstract: We characterize the role of charge correlations in the adsorption of a short, rod-like anionic polyelectrolyte onto a similarly charged membrane. Our theory reveals two different mechanisms driving the like-charge polyelectrolyte-membrane complexation: in weakly charged membranes, repulsive polyelectrolyte-membrane interactions lead to the interfacial depletion and a parallel orientation of the po… ▽ More We characterize the role of charge correlations in the adsorption of a short, rod-like anionic polyelectrolyte onto a similarly charged membrane. Our theory reveals two different mechanisms driving the like-charge polyelectrolyte-membrane complexation: in weakly charged membranes, repulsive polyelectrolyte-membrane interactions lead to the interfacial depletion and a parallel orientation of the polyelectrolyte with respect to the membrane; while in the intermediate membrane charge regime, the interfacial counterion excess gives rise to an attractive "salt-induced" image force. This furthermore results in an orientational transition from a parallel to a perpendicular configuration and a subsequent short-ranged like-charge adsorption of the polyelectrolyte to the substrate. A further increase of the membrane charge engenders a charge inversion, originating from surface-induced ionic correlations, that act as a separate mechanism capable of triggering the like-charge polyelectrolyte-membrane complexation over an extended distance interval from the membrane surface. The emerging picture of this complexation phenomenon identifies the interfacial "salt-induced" image forces as a powerful control mechanism in polyelectrolyte-membrane complexation. △ Less

Submitted 3 June, 2019; v1 submitted 14 March, 2019; originally announced March 2019.

Journal ref: Phys. Rev. E 99, 062501 (2019)

Pseudo-Casimir interactions and surface anchoring duality in bookshelf geometry of smectic-A liquid crystals

Authors: Fahimeh Karimi Pour Haddadan, Ali Naji, Rudolf Podgornik

Abstract: We analyze the transverse intersubstrate pseudo-Casimir force, arising as a result of thermal fluctuations of the liquid crystalline layers of a smectic-A film confined between two planar substrates in a bookshelf geometry, in which the equidistant smectic layers are placed perpendicular to the bounding surfaces. We discuss the variation of the interaction force as a function of the intersubstrate… ▽ More We analyze the transverse intersubstrate pseudo-Casimir force, arising as a result of thermal fluctuations of the liquid crystalline layers of a smectic-A film confined between two planar substrates in a bookshelf geometry, in which the equidistant smectic layers are placed perpendicular to the bounding surfaces. We discuss the variation of the interaction force as a function of the intersubstrate separation in the presence of surface anchoring to the substrates, showing that the force induced by confined fluctuations is attractive and depends on the penetration length as well as the layer spacing. The strongest effect occurs for tightly confined fluctuations, in which the surface anchoring energy is set to infinity, where the force per area scales linearly with the thermal energy and inversely with the fourth power of the intersubstrate separation. By reducing the strength of the surface anchoring energy, the force first becomes weaker in magnitude but then increases in magnitude as the surface anchoring strength is further reduced down to zero, in which case the force tends to that obtained in the limit of strong anchoring. △ Less

Submitted 27 November, 2018; originally announced November 2018.

Journal ref: Soft Matter 15, 2216 (2019)

Density--nematic coupling in isotropic linear polymers

Authors: Aleksandar Popadić, Daniel Svenšek, Rudolf Podgornik, Matej Praprotnik

Abstract: Linear polymers and other connected "line liquids" exhibit a coupling between density and equilibrium nematic order on the macroscopic level that gives rise to a Meyer-de Gennes vectorial conservation law. Nevertheless, isotropic linear polymer melts/solutions exhibit fluctuations of the density and of the nematic order that are not coupled by this vectorial constraint, just like for isotropic liq… ▽ More Linear polymers and other connected "line liquids" exhibit a coupling between density and equilibrium nematic order on the macroscopic level that gives rise to a Meyer-de Gennes vectorial conservation law. Nevertheless, isotropic linear polymer melts/solutions exhibit fluctuations of the density and of the nematic order that are not coupled by this vectorial constraint, just like for isotropic liquids composed of disconnected non-spherical particles. It takes the proper tensorial description of the nematic order in linear polymer liquids, leading to a tensorial conservation law connecting density and orientational order, that finally implicates coupled density and nematic order fluctuations, already in the isotropic system and not subject to the existence of an orientational phase transition. This coupling implies that a spatial variation of density or a local concentration gradient will induce nematic order and thereby an acoustic or osmotic optical birefringence even in an otherwise isotropic polymer melt/solution. We validate the theoretical conceptions by performing detailed Monte Carlo simulations of isotropic melts of "soft" worm-like chains with variable length and flexibility, and comparing the numerically determined orientation correlation functions with predictions of the macroscopic theory. The methodology drawn sets forth a means of determining the macroscopic parameters by microscopic simulations to yield realistic continuum models of specific polymeric materials. △ Less

Submitted 4 December, 2018; v1 submitted 13 November, 2018; originally announced November 2018.

Comments: 15 pages, 6 figures

Journal ref: Adv. Theory Simul. 2 (2019)

On the voltage-controlled assembly of NP arrays at electrochemical solid/liquid interfaces

Authors: Cristian Zagar, Ryan-Rys Griffith, Rudolf Podgornik, Alexei A. Kornyshev

Abstract: Research in the field of nanoplasmonic metamaterials is moving towards more and more interesting and, potentially useful, applications. The present work tackles the problem of nanoparticle self-assembly at an electrochemical solid-liquid interface from a purely theoretical perspective. We perform a simplified, comprehensive analysis of the stability of a nanoparticle arrays under different conditi… ▽ More Research in the field of nanoplasmonic metamaterials is moving towards more and more interesting and, potentially useful, applications. The present work tackles the problem of nanoparticle self-assembly at an electrochemical solid-liquid interface from a purely theoretical perspective. We perform a simplified, comprehensive analysis of the stability of a nanoparticle arrays under different conditions and assembly. From the Poisson-Bjerrum model of electrostatic interactions between a metallic nanoparticle and the electrode and between the nanoparticles at the electrode, as well the Hamaker-Lifshitz model of the corresponding van der Waals interactions, we reach some conclusions regarding the possibility to build arrays of charged nanoparticles on electrodes and disassemble them, subject to variation of applied voltage. Since system of this type have been shown, recently, to provide nontrivial electrotuneable optical response, such analysis is crucial for answering the question whether such scenarios of electrochemical plasmonics are feasible. △ Less

Submitted 11 October, 2018; originally announced October 2018.

Charge Regulation with Fixed and Mobile Charges

Authors: Yael Avni, David Andelman, Rudolf Podgornik

Abstract: Uncompensated charges do not occur in Nature and any local charge should be a result of charge separation. Dissociable chemical groups at interfaces in contact with ions in solution, whose chemical equilibrium depends both on short-range non-electrostatic and long-range electrostatic interactions, are the physical basis of this charge separation, known as charge regulation phenomena. The charged g… ▽ More Uncompensated charges do not occur in Nature and any local charge should be a result of charge separation. Dissociable chemical groups at interfaces in contact with ions in solution, whose chemical equilibrium depends both on short-range non-electrostatic and long-range electrostatic interactions, are the physical basis of this charge separation, known as charge regulation phenomena. The charged groups can be either fixed and immobile, as in the case of solvent-exposed solid substrate and soft bounding surfaces, (e.g., molecularly smooth mica surfaces and soft phospholipid membranes), or free and mobile, as in the case of charged macro-ions, (e.g., protein or other biomolecules). Here, we review the mean-field formalism used to describe both cases, with a focus on recent advances in the modeling of mobile charge-regulated macro-ions in an ionic solution. The general form of the screening length in such a solution is derived, and is shown to combine the concept of intrinsic capacitance, introduced by Lund and Jönsson, and bulk capacitance, resulting from the mobility of small ions and macro-ions. The advantages and disadvantages of different formulations, such as the cell model vs. the collective approach, are discussed, along with several suggestions for future experiments and modeling. △ Less

Submitted 21 August, 2018; originally announced August 2018.

Journal ref: Curr. Opin. Electrochem. 13, 70 (2019)

Role of metallic core for the stability of virus-like particles in strongly coupled electrostatics

Authors: Leili Javidpour, Anže Lošdorfer Božič, Rudolf Podgornik, Ali Naji

Abstract: We investigate the osmotic (electrostatic) pressure acting on the proteinaceous shell of a generic model of virus-like particles (VLPs), comprising a charged outer shell and a metallic nanoparticle core, coated by a charged layer and bathed in an aqueous electrolyte (salt) solution. Motivated by the recent studies accentuating the role of multivalent ions for the stability of VLPs, we focus on the… ▽ More We investigate the osmotic (electrostatic) pressure acting on the proteinaceous shell of a generic model of virus-like particles (VLPs), comprising a charged outer shell and a metallic nanoparticle core, coated by a charged layer and bathed in an aqueous electrolyte (salt) solution. Motivated by the recent studies accentuating the role of multivalent ions for the stability of VLPs, we focus on the effects of multivalent cations and anions in an otherwise monovalent ionic bathing solution. We perform extensive Monte-Carlo simulations based on appropriate Coulombic interactions that consistently take into account the effects of salt screening, the dielectric polarization of the metallic core, and the strong-coupling electrostatics due to the presence of multivalent ions. We specifically study the intricate roles these factors play in the electrostatic stability of the model VLPs. It is shown that while the insertion of a metallic nanoparticle by itself can produce negative, inward-directed, pressure on the outer shell, addition of only a small amount of multivalent counterions can robustly engender negative pressures, enhancing the VLP stability across a wide range of values for the system parameters. △ Less

Submitted 16 August, 2018; originally announced August 2018.

Comments: 15 pages, 7 figures

Journal ref: Sci. Rep. 9, 3884 (2019)

Rodlike counterions at heterogeneously charged surfaces

Authors: Ali Naji, Kasra Hejazi, Elnaz Mahgerefteh, Rudolf Podgornik

Abstract: We study the spatial and orientational distribution of rodlike counterions (such as mobile nanorods) as well as the effective interaction mediated by them between two plane-parallel surfaces that carry fixed (quenched) heterogeneous charge distributions. The rodlike counterions are assumed to have an internal charge distribution, specified by a multivalent monopolar moment and a finite quadrupolar… ▽ More We study the spatial and orientational distribution of rodlike counterions (such as mobile nanorods) as well as the effective interaction mediated by them between two plane-parallel surfaces that carry fixed (quenched) heterogeneous charge distributions. The rodlike counterions are assumed to have an internal charge distribution, specified by a multivalent monopolar moment and a finite quadrupolar moment, and the quenched surface charge is assumed to be randomly distributed with equal mean and variance on the two surfaces. While equally charged surfaces are known to repel within the traditional mean-field theories, the presence of multivalent counterions has been shown to cause attractive interactions between uniformly charged surfaces due to the prevalence of strong electrostatic couplings that grow rapidly with the counterion valency. We show that the combined effects due to electrostatic correlations (caused by the coupling between the mean surface field and the multivalent, monopolar, charge valency of counterions) as well as the disorder-induced interactions (caused by the coupling between the surface disorder field and the quadrupolar moment of counterions) lead to much stronger attractive interactions between two randomly charged surfaces. The interaction profile turns out to be a nonmonotonic function of the intersurface separation, displaying an attractive minimum at relatively small separations, where the ensuing attraction can exceed the maximum strong-coupling attraction (produced by multivalent monopolar counterions between uniformly charged surfaces) by more than an order of magnitude. △ Less

Submitted 7 September, 2018; v1 submitted 12 June, 2018; originally announced June 2018.

Comments: 15 pages, 4 figures

Journal ref: J. Chem. Phys. 149, 134702 (2018)

General theory of charge regulation and surface capacitance

Authors: Rudolf Podgornik

Abstract: A generalization of the mean-field approach will be derived that will take into account the ion-ion as well as ion-surface non-electrostatic effects on an equal footing, being based on the bulk and surface equations of state in the absence of electrostatic interactions. This approach will be applied to the calculation of the differential capacitance for several models of the specific ion-surface n… ▽ More A generalization of the mean-field approach will be derived that will take into account the ion-ion as well as ion-surface non-electrostatic effects on an equal footing, being based on the bulk and surface equations of state in the absence of electrostatic interactions. This approach will be applied to the calculation of the differential capacitance for several models of the specific ion-surface non-electrostatic interactions providing a thermodynamic insight into the characteristics of the differential capacitance. The packing considerations due to surface adsorbed ions will be shown to be more relevant then the bulk packing constraints for ions vicinal to the surface, and to set in prior to the conditions where the bulk packing constraints would become relevant. △ Less

Submitted 18 June, 2018; v1 submitted 6 June, 2018; originally announced June 2018.

Mean-field theory of active electrolytes: dynamic adsorption and overscreening

Authors: Derek Frydel, Rudolf Podgornik

Abstract: We investigate active electrolytes within the mean-field level of description. The focus is on how the double-layer structure of passive, thermalized charges is affected by active dynamics of all constituting ions. One feature of active dynamics is that particles adhere to hard-surfaces, regardless of chemical properties of a surface and specifically in complete absence of any chemipsorption or ph… ▽ More We investigate active electrolytes within the mean-field level of description. The focus is on how the double-layer structure of passive, thermalized charges is affected by active dynamics of all constituting ions. One feature of active dynamics is that particles adhere to hard-surfaces, regardless of chemical properties of a surface and specifically in complete absence of any chemipsorption or physisorption. To carry out the mean-field analysis of the system that is out of equilibrium, we develop the "mean-field simulation" technique, where the simulated system consists of charged parallel sheets moving on a line and obeying active dynamics, with the interaction strength rescaled by the number of sheets. The mean-field limit becomes exact in the limit of an infinite number of movable sheets. △ Less

Submitted 19 May, 2018; v1 submitted 30 April, 2018; originally announced May 2018.

Journal ref: Phys. Rev. E 97, 052609 (2018)

Anomalous Multipole Expansion: Charge Regulation of Patchy, Inhomogeneously Charged Spherical Particles

Authors: Anže Lošdorfer Božič, Rudolf Podgornik

Abstract: Charge regulation is an important aspect of electrostatics in biological and colloidal systems, where the charges are generally not fixed, but depend on the environmental variables. Here, we analyze the charge regulation mechanism in patchy, inhomogeneously charged spherical particles, such as globular proteins, colloids, or viruses. Together with the multipole expansion of inhomogeneously charged… ▽ More Charge regulation is an important aspect of electrostatics in biological and colloidal systems, where the charges are generally not fixed, but depend on the environmental variables. Here, we analyze the charge regulation mechanism in patchy, inhomogeneously charged spherical particles, such as globular proteins, colloids, or viruses. Together with the multipole expansion of inhomogeneously charged spherical surfaces, the charge regulation mechanism on the level of linear approximation is shown to lead to a mixing between different multipole moments depending on their capacitance -- the response function of the charge distribution with respect to the electrostatic potential. This presents an additional anomalous feature of molecular electrostatics in the presence of ionic screening. We demonstrate the influence of charge regulation on several examples of inhomogeneously charged spherical particles, showing that it leads to significant changes in their multipole moments. △ Less

Submitted 31 May, 2018; v1 submitted 20 April, 2018; originally announced April 2018.

Comments: 9 pages, 4 figures

Journal ref: J. Chem. Phys. 149, 163308 (2018)

Charge Regulating Macro-ions in Salt Solutions: Screening Properties and Electrostatic Interactions

Authors: Yael Avni, David Andelman, Tomer Markovich, Rudi Podgornik

Abstract: We revisit the charge-regulation mechanism of macro-ions and apply it to mobile macro-ions in a bathing salt solution. In particular, we examine the effects of correlation between various adsorption/desorption sites and analyze the collective behavior in terms of the solution effective screening properties. We show that such a behavior can be quantified in terms of the charge {\em asymmetry} of th… ▽ More We revisit the charge-regulation mechanism of macro-ions and apply it to mobile macro-ions in a bathing salt solution. In particular, we examine the effects of correlation between various adsorption/desorption sites and analyze the collective behavior in terms of the solution effective screening properties. We show that such a behavior can be quantified in terms of the charge {\em asymmetry} of the macro-ions, defined by their preference for a non-zero effective charge, and their {\em donor/acceptor} propensity for exchanging salt ions with the bathing solution. Asymmetric macro-ions tend to increase the screening, while symmetric macro-ions can in some cases decrease it. Macro-ions that are classified as donors display a rather regular behavior, while those that behave as acceptors exhibit an anomalous non-monotonic Debye length. The screening properties, in their turn, engender important modifications to the disjoining pressure between two charged surfaces. Our findings are in particular relevant for solutions of proteins, whose exposed amino acids can undergo charge dissociation/association processes to/from the bathing solution, and can be considered as a solution of charged regulated macro-ions, as analyzed here. △ Less

Submitted 31 July, 2018; v1 submitted 9 April, 2018; originally announced April 2018.

Comments: 13 pages, 8 figures, minor revisions

Journal ref: Soft Matter 14, 6058 (2018)

Splay--density coupling in semiflexible main-chain nematic polymers with hairpins

Authors: Aleksandar Popadić, Daniel Svenšek, Rudolf Podgornik, Kostas Ch. Daoulas, Matej Praprotnik

Abstract: We establish a macroscopic description of the splay--density coupling in semiflexible main-chain nematic polymers with hairpins, using a vectorial continuity constraint for the "recovered" polar order of the chain tangents and introducing chain backfolds (hairpins) as its new type of sources besides chain ends. We treat both types of sources on a unified basis as a mixture of two ideal gases with… ▽ More We establish a macroscopic description of the splay--density coupling in semiflexible main-chain nematic polymers with hairpins, using a vectorial continuity constraint for the "recovered" polar order of the chain tangents and introducing chain backfolds (hairpins) as its new type of sources besides chain ends. We treat both types of sources on a unified basis as a mixture of two ideal gases with fixed composition. Performing detailed Monte Carlo simulations of nematic monodomain melts of "soft" worm-like chains with variable length and flexibility, we show via their structure factors that the chain backfolding weakens the splay--density coupling, and demonstrate how this weakening can be consistently quantified on the macroscopic level. We also probe and discuss the deviations from the noninteracting gas idealization of the chain ends and backfolds. △ Less

Submitted 14 November, 2018; v1 submitted 21 February, 2018; originally announced February 2018.

Comments: 6 pages, 5 figures

Journal ref: Soft Matter, 2018, 14, 5898

Hydrodynamic stress correlations in fluid films driven by stochastic surface forcing

Authors: Masoud Mohammadi-Arzanagh, Saeed Mahdisoltani, Rudolf Podgornik, Ali Naji

Abstract: We study hydrodynamic fluctuations in a compressible and viscous fluid film confined between two rigid, no-slip, parallel plates, where one of the plates is kept fixed, while the other one is driven in small-amplitude, translational, displacements around its reference position. This jiggling motion is assumed to be driven by a stochastic, external, surface forcing of zero mean and finite variance.… ▽ More We study hydrodynamic fluctuations in a compressible and viscous fluid film confined between two rigid, no-slip, parallel plates, where one of the plates is kept fixed, while the other one is driven in small-amplitude, translational, displacements around its reference position. This jiggling motion is assumed to be driven by a stochastic, external, surface forcing of zero mean and finite variance. Thus, while the transverse (shear) and longitudinal (compressional) hydrodynamic stresses produced in the film vanish on average on either of the plates, these stresses exhibit fluctuations that can be quantified through their equal-time, two-point, correlation functions. For transverse stresses, we show that the correlation functions of the stresses acting on the same plate (self-correlators) as well as the correlation function of the stresses acting on different plates (cross-correlators) exhibit universal, decaying, power-law behaviors as functions of the inter-plate separation. At small separations, the exponents are given by -1, while at large separations, the exponents are found as -2 (self-correlator on the fixed plate), -4 (excess self-correlator on the mobile plate) and -3 (cross-correlator). For longitudinal stresses, we find much weaker power-law decays in the large separation regime, with exponents -3/2 (excess self-correlator on the mobile plate) and -1 (cross-correlator). The self-correlator on the fixed plate increases and levels off upon increasing the inter-plate separation, reflecting the non-decaying nature of the longitudinal forces acting on the fixed plate. △ Less

Submitted 30 June, 2018; v1 submitted 13 December, 2017; originally announced December 2017.

Comments: 14 pages, 10 figures

Journal ref: Phys. Rev. Fluids 3, 064201 (2018)

Varieties of charge distributions in coat proteins of ssRNA+ viruses

Authors: Anze Losdorfer Bozic, Rudolf Podgornik

Abstract: A major part of the interactions involved in the assembly and stability of icosahedral, positive-sense single-stranded RNA (ssRNA+) viruses is electrostatic in nature, as can be inferred from the strong $pH$- and salt-dependence of their assembly phase diagrams. Electrostatic interactions do not act only between the capsid coat proteins (CPs), but just as often provide a significant contribution t… ▽ More A major part of the interactions involved in the assembly and stability of icosahedral, positive-sense single-stranded RNA (ssRNA+) viruses is electrostatic in nature, as can be inferred from the strong $pH$- and salt-dependence of their assembly phase diagrams. Electrostatic interactions do not act only between the capsid coat proteins (CPs), but just as often provide a significant contribution to the interactions of the CPs with the genomic RNA, mediated to a large extent by positively charged, flexible N-terminal tails of the CPs. In this work, we provide two clear and complementary definitions of an N-terminal tail of a protein, and use them to extract the tail sequences of a large number of CPs of ssRNA+ viruses. We examine the $pH$-dependent interplay of charge on both tails and CPs alike, and show that -- in contrast to the charge on the CPs -- the net positive charge on the N-tails persists even to very basic $pH$ values. In addition, we note a limit to the length of the wild-type genomes of those viruses which utilize positively charged tails, when compared to viruses without charged tails and similar capsid size. At the same time, we observe no clear connection between the charge on the N-tails and the genome lengths of the viruses included in our study. △ Less

Submitted 29 November, 2017; v1 submitted 21 September, 2017; originally announced September 2017.

Comments: 12 pages, 7 figures; 24 pages, 10 figures, and 5 tables of Supplementary Material; accepted for publication in JP:CM

Journal ref: J. Phys. Condens. Matter 30, 024001 (2018)