Source author record

David S. Dean

David S. Dean 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

43works
14topics
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

43 published item(s)

preprint2022arXiv

How stickiness can speed up diffusion in confined systems

The paradigmatic model for heterogeneous media used in diffusion studies is built from reflecting obstacles and surfaces. It is well known that the crowding effect produced by these reflecting surfaces slows the dispersion of Brownian tracers. Here, using a general adsorption desorption model with surface diffusion, we show analytically that making surfaces or obstacles attractive can accelerate dispersion. In particular, we show that this enhancement of diffusion can exist even when the surface diffusion constant is smaller than that in the bulk. Even more remarkably, this enhancement effect occurs when the effective diffusion constant, when restricted to surfaces only, is lower than the effective diffusivity with purely reflecting boundaries. We give analytical formulas for this intriguing effect in periodic arrays of spheres as well as undulating micro-channels. Our results are confirmed by numerical calculations and Monte Carlo simulations.

preprint2020arXiv

Kernels for noninteracting fermions via a Green's function approach with applications to step potentials

The quantum correlations of $N$ noninteracting spinless fermions in their ground state can be expressed in terms of a two-point function called the kernel. Here we develop a general and compact method for computing the kernel in a general trapping potential in terms of the Green's function for the corresponding single particle Schrödinger equation. For smooth potentials the method allows a simple alternative derivation of the local density approximation for the density and of the sine kernel in the bulk part of the trap in the large $N$ limit. It also recovers the density and the kernel of the so-called {\em Airy gas} at the edge. This method allows to analyse the quantum correlations in the ground state when the potential has a singular part with a fast variation in space. For the square step barrier of height $V_0$, we derive explicit expressions for the density and for the kernel. For large Fermi energy $μ>V_0$ it describes the interpolation between two regions of different densities in a Fermi gas, each described by a different sine kernel. Of particular interest is the {\em critical point} of the square well potential when $μ=V_0$. In this critical case, while there is a macroscopic number of fermions in the lower part of the step potential, there is only a finite $O(1)$ number of fermions on the shoulder, and moreover this number is independent of $μ$. In particular, the density exhibits an algebraic decay $\sim 1/x^2$, where $x$ is the distance from the jump. Furthermore, we show that the critical behaviour around $μ= V_0$ exhibits universality with respect with the shape of the barrier. This is established (i) by an exact solution for a smooth barrier (the Woods-Saxon potential) and (ii) by establishing a general relation between the large distance behavior of the kernel and the scattering amplitudes of the single-particle wave-function.

preprint2020arXiv

Noninteracting trapped Fermions in double-well potentials: inverted parabola kernel

We study a system of $N$ noninteracting spinless fermions in a confining, double-well potential in one dimension. When the Fermi energy is close to the value of the potential at its local maximum we show that physical properties, such as the average density and the fermion position correlation functions, display a universal behavior that depends only on the local properties of the potential near its maximum. This behavior describes the merging of two Fermi gases, which are disjoint at sufficiently low Fermi energies. We describe this behavior in terms of a new correlation kernel that we compute analytically and we call it the inverted parabola kernel". As an application, we calculate the mean and variance of the number of particles in an interval of size $2L$ centered around the position of the local maximum, for sufficiently small $L$. Finally, we discuss the possibility of observing our results in experiments, as well as the extensions to nonzero temperature and to higher space dimensions.

preprint2020arXiv

Position distribution in a generalised run and tumble process

We study a class of stochastic processes of the type $\frac{d^n x}{dt^n}= v_0\, σ(t)$ where $n>0$ is a positive integer and $σ(t)=\pm 1$ represents an `active' telegraphic noise that flips from one state to the other with a constant rate $γ$. For $n=1$, it reduces to the standard run and tumble process for active particles in one dimension. This process can be analytically continued to any $n>0$ including non-integer values. We compute exactly the mean squared displacement at time $t$ for all $n>0$ and show that at late times while it grows as $\sim t^{2n-1}$ for $n>1/2$, it approaches a constant for $n<1/2$. In the marginal case $n=1/2$, it grows very slowly with time as $\sim \ln t$. Thus the process undergoes a {\em localisation} transition at $n=1/2$. We also show that the position distribution $p_n(x,t)$ remains time-dependent even at late times for $n\ge 1/2$, but approaches a stationary time-independent form for $n<1/2$. The tails of the position distribution at late times exhibit a large deviation form, $p_n(x,t)\sim \exp\left[-γ\, t\, Φ_n\left(\frac{x}{x^*(t)}\right)\right]$, where $x^*(t)= v_0\, t^n/Γ(n+1)$. We compute the rate function $Φ_n(z)$ analytically for all $n>0$ and also numerically using importance sampling methods, finding excellent agreement between them. For three special values $n=1$, $n=2$ and $n=1/2$ we compute the exact cumulant generating function of the position distribution at all times $t$.

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

Non-interacting fermions at finite temperature in a $d$-dimensional trap: universal correlations

We study a system of $N$ non-interacting spin-less fermions trapped in a confining potential, in arbitrary dimensions $d$ and arbitrary temperature $T$. The presence of the trap introduces an edge where the average density of fermions vanishes. Far from the edge, near the center of the trap (the so called "bulk regime"), physical properties of the fermions have traditionally been understood using the Local Density Approximation. However, this approximation drastically fails near the edge where the density vanishes. In this paper we show that, even near the edge, novel universal properties emerge, independently of the details of the confining potential. We show that for large $N$, these fermions in a confining trap, in arbitrary dimensions and at finite temperature, form a determinantal point process. As a result, any $n$-point correlation function can be expressed as an $n \times n$ determinant whose entry is called the kernel. Near the edge, we derive the large $N$ scaling form of the kernels. In $d=1$ and $T=0$, this reduces to the so called Airy kernel, that appears in the Gaussian Unitary Ensemble (GUE) of random matrix theory. In $d=1$ and $T>0$ we show a remarkable connection between our kernel and the one appearing in the $1+1$-dimensional Kardar-Parisi-Zhang equation at finite time. Consequently our result provides a finite $T$ generalization of the Tracy-Widom distribution, that describes the fluctuations of the rightmost fermion at $T=0$. In $d>1$ and $T \geq 0$, while the connection to GUE no longer holds, the process is still determinantal whose analysis provides a new class of kernels, generalizing the $1d$ Airy kernel at $T=0$ obtained in random matrix theory. Some of our finite temperature results should be testable in present-day cold atom experiments, most notably our detailed predictions for the temperature dependence of the fluctuations near the edge.

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

Universal ground state properties of free fermions in a $d$-dimensional trap

The ground state properties of $N$ spinless free fermions in a $d$-dimensional confining potential are studied. We find that any $n$-point correlation function has a simple determinantal structure that allows us to compute several properties exactly for large $N$. We show that the average density has a finite support with an edge, and near this edge the density exhibits a universal (valid for a wide class of potentials) scaling behavior for large $N$. The associated edge scaling function is computed exactly and generalizes to any $d$ the edge electron gas result of Kohn and Mattsson in $d=3$ [Phys. Rev. Lett. 81, 3487 (1998)]. In addition, we calculate the kernel (that characterizes any $n$-point correlation function) for large $N$ and show that, when appropriately scaled, it depends only on dimension $d$, but has otherwise universal scaling forms, at the edges. The edge kernel, for higher $d$, generalizes the Airy kernel in one dimension, well known from random matrix theory.

preprint2015arXiv

Finite temperature free fermions and the Kardar-Parisi-Zhang equation at finite time

We consider the system of $N$ one-dimensional free fermions confined by a harmonic well $V(x) = mω^2 {x^2}/{2}$ at finite inverse temperature $β= 1/T$. The average density of fermions $ρ_N(x,T)$ at position $x$ is derived. For $N \gg 1$ and $β\sim {\cal O}(1/N)$, $ρ_N(x,T)$ is given by a scaling function interpolating between a Gaussian at high temperature, for $β\ll 1/N$, and the Wigner semi-circle law at low temperature, for $β\gg N^{-1}$. In the latter regime, we unveil a scaling limit, for $β{\hbar ω}= b N^{-1/3}$, where the fluctuations close to the edge of the support, at $x \sim \pm \sqrt{2\hbar N/(mω)}$, are described by a limiting kernel $K^{\rm ff}_b(s,s')$ that depends continuously on $b$ and is a generalization of the Airy kernel, found in the Gaussian Unitary Ensemble of random matrices. Remarkably, exactly the same kernel $K^{\rm ff}_b(s,s')$ arises in the exact solution of the Kardar-Parisi-Zhang (KPZ) equation in 1+1 dimensions at finite time $t$, with the correspondence $t= b^3$.

preprint2014arXiv

Approach to asymptotically diffusive behavior for Brownian particles in media with periodic diffusivities

We analyze the mean squared displacement of a Brownian particle in a medium with a spatially varying local diffusivity which is assumed to be periodic. When the system is asymptotically diffusive the mean squared displacement, characterizing the dispersion in the system, is, at late times, a linear function of time. A Kubo type formula is given for the mean squared displacement which allows the recovery of some known results for the effective diffusion constant $D_e$ in a direct way, but also allows an understanding of the asymptotic approach to the diffusive limit. In particular, as well as as computing the slope of a linear fit to the late time mean squared displacement, we find a formula for the constant where the fit intersects the y axis.

preprint2014arXiv

Approach to asymptotically diffusive behavior for Brownian particles in periodic potentials : extracting information from transients

A Langevin process diffusing in a periodic potential landscape has a time dependent diffusion constant which means that its average mean squared displacement (MSD) only becomes linear at late times. The long time, or effective diffusion constant, can be estimated from the slope of a linear fit of the MSD at late times. Due to the cross over between a short time microscopic diffusion constant, which is independent of the potential, to the effective late time diffusion constant, a linear fit of the MSD will not in general pass through the origin and will have a non-zero constant term. Here we address how to compute the constant term and provide explicit results for Brownian particles in one dimension in periodic potentials. We show that the constant is always positive and that at low temperatures it depends on the curvature of the minimum of the potential. For comparison we also consider the same question for the simpler problem of a symmetric continuous time random walk in discrete space. Here the constant can be positive or negative and can be used to determine the variance of the hopping time distribution.

preprint2014arXiv

Diffusion in periodic, correlated random forcing landscapes

We study the dynamics of a Brownian particle in a strongly correlated quenched random potential defined as a periodically-extended (with period $L$) finite trajectory of a fractional Brownian motion with arbitrary Hurst exponent $H \in (0,1)$. While the periodicity ensures that the ultimate long-time behavior is diffusive, the generalised Sinai potential considered here leads to a strong logarithmic confinement of particle trajectories at intermediate times. These two competing trends lead to dynamical frustration and result in a rich statistical behavior of the diffusion coefficient $D_L$: Although one has the typical value $D^{\rm typ}_L \sim \exp(-βL^H)$, we show via an exact analytical approach that the positive moments ($k>0$) scale like $\langle D^k_L \rangle \sim \exp{[-c' (k βL^{H})^{1/(1+H)}]}$, and the negative ones as $\langle D^{-k}_L \rangle \sim \exp(a' (k βL^{H})^2)$, $c'$ and $a'$ being numerical constants and $β$ the inverse temperature. These results demonstrate that $D_L$ is strongly non-self-averaging. We further show that the probability distribution of $D_L$ has a log-normal left tail and a highly singular, one-sided log-stable right tail reminiscent of a Lifshitz singularity.

preprint2014arXiv

Fluctuation mediated interactions due to rigidity mismatch and their effect on miscibility of lipid mixtures in multicomponent membranes

We consider how membrane fluctuations can modify the miscibility of lipid mixtures, that is to say how the phase diagram of a boundary-constrained membrane is modified when the membrane is allowed to fluctuate freely in the case of zero surface tension. In order for fluctuations to have an effect, the different lipid types must have differing Gaussian rigidities. We show, somewhat paradoxically, that fluctuation-induced interactions can be treated approximately in a mean-field type theory. Our calculations predict that, depending on the difference in bending and Gaussian rigidity of the lipids, membrane fluctuations can either favor or disfavor mixing.

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

Distribution of the least-squares estimators of a single Brownian trajectory diffusion coefficient

In this paper we study the distribution function $P(u_α)$ of the estimators $u_α \sim T^{-1} \int^T_0 \, ω(t) \, {\bf B}^2_{t} \, dt$, which optimise the least-squares fitting of the diffusion coefficient $D_f$ of a single $d$-dimensional Brownian trajectory ${\bf B}_{t}$. We pursue here the optimisation further by considering a family of weight functions of the form $ω(t) = (t_0 + t)^{-α}$, where $t_0$ is a time lag and $α$ is an arbitrary real number, and seeking such values of $α$ for which the estimators most efficiently filter out the fluctuations. We calculate $P(u_α)$ exactly for arbitrary $α$ and arbitrary spatial dimension $d$, and show that only for $α= 2$ the distribution $P(u_α)$ converges, as $ε= t_0/T \to 0$, to the Dirac delta-function centered at the ensemble average value of the estimator. This allows us to conclude that only the estimators with $α= 2$ possess an ergodic property, so that the ensemble averaged diffusion coefficient can be obtained with any necessary precision from a single trajectory data, but at the expense of a progressively higher experimental resolution. For any $α\neq 2$ the distribution attains, as $ε\to 0$, a certain limiting form with a finite variance, which signifies that such estimators are not ergodic.

preprint2013arXiv

On ergodic least-squares estimators of the generalized diffusion coefficient for fractional Brownian motion

We analyse a class of estimators of the generalized diffusion coefficient for fractional Brownian motion $B_t$ of known Hurst index $H$, based on weighted functionals of the single time square displacement. We show that for a certain choice of the weight function these functionals possess an ergodic property and thus provide the true, ensemble-averaged, generalized diffusion coefficient to any necessary precision from a single trajectory data, but at expense of a progressively higher experimental resolution. Convergence is fastest around $H\simeq0.30$, a value in the subdiffusive regime.

preprint2013arXiv

Optimal least-squares estimators of the diffusion constant from a single Brownian trajectory

Modern developments in microscopy and image processing are revolutionising areas of physics, chemistry, and biology as nanoscale objects can be tracked with unprecedented accuracy. However, the price paid for having a direct visualisation of a single particle trajectory with high temporal and spatial resolution is a consequent lack of statistics. This naturally calls for reliable analytical tools which will allow one to extract the properties specific to a statistical ensemble from just a single trajectory. In this article we briefly survey different analytical methods currently used to determine the ensemble average diffusion coefficient from single particle data and then focus specifically on weighted least-squares estimators, seeking the weight functions for which such estimators are ergodic. Finally, we address the question of the effects of disorder on such estimators.

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

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

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

Optimal estimates of the diffusion coefficient of a single Brownian trajectory

Modern developments in microscopy and image processing are revolutionizing areas of physics, chemistry and biology as nanoscale objects can be tracked with unprecedented accuracy. The goal of single particle tracking is to determine the interaction between the particle and its environment. The price paid for having a direct visualization of a single particle is a consequent lack of statistics. Here we address the optimal way of extracting diffusion constants from single trajectories for pure Brownian motion. It is shown that the maximum likelihood estimator is much more efficient than the commonly used least squares estimate. Furthermore we investigate the effect of disorder on the distribution of estimated diffusion constants and show that it increases the probability of observing estimates much smaller than the true (average) value.

preprint2012arXiv

Optimal fits of diffusion constants from single time data points of Brownian trajectories

Experimental methods based on single particle tracking (SPT) are being increasingly employed in the physical and biological sciences, where nanoscale objects are visualized with high temporal and spatial resolution. SPT can probe interactions between a particle and its environment but the price to be paid is the absence of ensemble averaging and a consequent lack of statistics. Here we address the benchmark question of how to accurately extract the diffusion constant of one single Brownian trajectory. We analyze a class of estimators based on weighted functionals of the square displacement. For a certain choice of the weight function these functionals provide the true ensemble averaged diffusion coefficient, with a precision that increases with the trajectory resolution.

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

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

Electromediated formation of DNA complexes with cell membranes and its consequences for gene delivery

Electroporation is a physical method to induce the uptake of therapeutic drugs and DNA, by eukaryotic cells and tissues. The phenomena behind electro-mediated membrane permeabilization to plasmid DNA have been shown to be significantly more complex than those for small molecules. Small molecules cross the permeabilized membrane by diffusion whereas plasmid DNA first interacts with the electropermeabilized part of the cell surface, forming localized aggregates. The dynamics of this process is still poorly understood because direct observations have been limited to scales of the order of seconds. Here, cells are electropermeabilized in the presence of plasmid DNA and monitored with a temporal resolution of 2 ms. This allows us to show that during the first pulse application, plasmid complexes, or aggregates, start to form at distinct sites on the cell membrane. FRAP measurements show that the positions of these sites are remarkably immobile during the application of further pluses. A theoretical model is proposed to explain the appearance of distinct interaction sites, the quantitative increase in DNA and also their immobility leading to a tentative explanation for the success of electro-mediated gene delivery.

preprint2011arXiv

Insights into the mechanisms of electromediated gene delivery and application to the loading of giant vesicles with negatively charged macromolecules

We present experimental results regarding the electrotransfer of plasmid DNA into phosphatidylcholine giant unilamellar vesicles (GUVs). Our observations indicate that a direct entry is the predominant mechanism of electrotransfer. A quantitative analysis of the DNA concentration increments inside the GUVs is also performed, and we find that our experimental data are very well described by a simple theoretical model in which DNA entry is mostly driven by electrophoresis. Our theoretical framework allows for the prediction of the amount of transfered DNA as a function of the electric field parameters, and thus paves the way towards a novel method for encapsulating with high efficiency not only DNA, but any negatively charged macromolecule into GUVs.

preprint2011arXiv

On the distribution of estimators of diffusion constants for Brownian motion

We discuss the distribution of various estimators for extracting the diffusion constant of single Brownian trajectories obtained by fitting the squared displacement of the trajectory. The analysis of the problem can be framed in terms of quadratic functionals of Brownian motion that correspond to the Euclidean path integral for simple Harmonic oscillators with time dependent frequencies. Explicit analytical results are given for the distribution of the diffusion constant estimator in a number of cases and our results are confirmed by numerical simulations.

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

Diffusion of active tracers in fluctuating fields

The problem of a particle diffusion in a fluctuating scalar field is studied. In contrast to most studies of advection diffusion in random fields we analyze the case where the particle position is also coupled to the dynamics of the field. Physical realizations of this problem are numerous and range from the diffusion of proteins in fluctuating membranes and the diffusion of localized magnetic fields in spin systems. We present exact results for the diffusion constant of particles diffusing in dynamical Gaussian fields in the adiabatic limit where the field evolution is much faster than the particle diffusion. In addition we compute the diffusion constant perturbatively, in the weak coupling limit where the interaction of the particle with the field is small, using a Kubo-type relation. Finally we construct a simple toy model which can be solved exactly.

preprint2010arXiv

Drag forces in classical fields

Inclusions, or defects, moving at constant velocity through free classical fields are shown to be subject to a drag force which depends on the field dynamics and the coupling of the inclusion to the field. The results are used to predict the drag exerted on inclusions, such as proteins, in lipid membranes due to their interaction with height and composition fluctuations. The force, measured in Monte Carlo simulations, on a point like magnetic field moving through an Ising ferromagnet is also well explained by these results.

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.

preprint2010arXiv

The out of equilibrium behavior of Casimir type fluctuation induced forces for free classical fields

We present a general method to study the non-equilibrium behavior of Casimir type fluctuation induced forces for classical free scalar field theories. In particular we analyze the temporal evolution of the force towards its equilibrium value when the field dynamics is given by a general class of over damped stochastic dynamics (including the model A and model B class). The steady state force is also analyzed for systems which have non-equilibrium steady states, for instance where they are driven by colored noise. The key to the method is that out of equilibrium force is computed by specifying an energy of interaction between the field and the surfaces in the problem. In general we find that there is a mapping of the dynamical problem onto a corresponding static one, and in the case where the latter can be solved the full dynamical behavior of the force can be extracted. The method is used how to compute the non-equilibrium Casimir force induced between two parallel plates by a fluctuating field, in the cases of Dirichlet, Neumann and mixed boundary conditions. Various other examples, such as the fluctuation induced force between inclusions in fluctuating media are discussed.

preprint2009arXiv

Dipole diffusion in a random electrical potential

We study the Langevin dynamics of a dipole diffusing in a random electrical field E derived from a quenched Gaussian potential. We show that in a suitable adiabatic limit (where the dynamics of the dipole moment is much faster than the dynamics of its position), one can reduce the coupled stochastic equations to an effective Langevin equation for a particle diffusing in an effective potential with a spatially varying and anisotropic local diffusivity $κ_{ij}$. Analytic results, close to the adiabatic limit, for the diffusion constant $κ_e$ are found in one dimension and a finite temperature dynamical transition is found. The system is also studied numerically. In particular, we study the anomalous diffusion exponent in the low temperature regime. Our findings strongly support the conclusion that the location of the dynamical transition and the anomalous diffusion exponents are determined by purely static considerations, i.e. they are independent of the relative values of the diffusion constants of the particle position and its dipole moment.

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.

preprint2004arXiv

Equilibrium statistics of a slave estimator in Langevin processes

We analyze the statistics of an estimator, denoted by xi_t and referred to as the slave, for the equilibrium susceptibility of a one dimensional Langevin process x_t in a potential phi(x). The susceptibility can be measured by evolving the slave equation in conjunction with the original Langevin process. This procedure yields a direct estimate of the susceptibility and avoids the need, when performing numerical simulations, to include applied external fields explicitly. The success of the method however depends on the statistical properties of the slave estimator. The joint probability density function for x_t and xi_t is analyzed. In the case where the potential of the system has a concave component the probability density function of the slave acquires a power law tail characterized by a temperature dependent exponent. Thus we show that while the average value of the slave, in the equilibrium state, is always finite and given by the fluctuation dissipation relation, higher moments and indeed the variance may show divergences. The behavior of the power law exponent is analyzed in a general context and it is calculated explicitly in some specific examples. Our results are confirmed by numerical simulations and we discuss possible measurement discrepancies in the fluctuation dissipation relation which could arise due to this behavior.

preprint2000arXiv

Slow Logarithmic Decay of Magnetization in the Zero Temperature Dynamics of an Ising Spin Chain: Analogy to Granular Compaction

We study the zero temperature coarsening dynamics in an Ising chain in presence of a dynamically induced field that favors locally the `-' phase compared to the `+' phase. At late times, while the `+' domains still coarsen as $t^{1/2}$, the `-' domains coarsen slightly faster as $t^{1/2}\log (t)$. As a result, at late times, the magnetization decays slowly as, $m(t)=-1 +{\rm const.}/{\log (t)}$. We establish this behavior both analytically within an independent interval approximation (IIA) and numerically. In the zero volume fraction limit of the `+' phase, we argue that the IIA becomes asymptotically exact. Our model can be alternately viewed as a simple Ising model for granular compaction. At late times in our model, the system decays into a fully compact state (where all spins are `-') in a slow logarithmic manner $\sim 1/{\log (t)}$, a fact that has been observed in recent experiments on granular systems.

preprint1995arXiv

Full dynamical solution for a spherical spin-glass model

We present a detailed analysis for the Langevin dynamics of a spherical spin-glass model (the spherical Sherrington-Kirkpatrick model). All the spins in the system are coupled by pairs via a random interaction matrix taken from the Gaussian ensemble. One finds that for a general initial configuration the system never reaches an equilibrium state and the theorems associated to `equilibrium dynamics' are violated. Only very particular initial conditions drive the system to equilibrium. The weak ergodicity breaking scenario is demonstrated for general `non-equilibrium' initial conditions. Two-time quantities such as the auto-correlation function explicitly depend on both times. When the time difference is short compared to the smaller time one finds stationary dynamics with time translation invariance and the fluctuation-dissipation theorem satisfied. Instead, when the time difference is of the order of the smaller time one finds non-stationary dynamics with aging phenomena, the system {\em remembers} the time spent after the initial time (the quench below the critical temperature). Interestingly enough the short time-difference dynamics (FDT regime) for non-equilibrium initial conditions is identical to the relaxation within the equilibrium states obtained for the particular `equilibrium' initial conditions. This points to a self-similaririty of the energy landascape. In addition we analyse the effect of temperature variations on the behaviour of the correlation function and energy density. In this way we are able to make some comparisons between experimentally observed results and exact calculations for the model. Somewhat surprisingly, this simple model captures some of the effects seen in laboratory spin-glasses.