Showing 1–50 of 140 results for author: Dean, D S

Singular viscoelastic perturbation to soft lubrication

Authors: Bharti Bharti, Quentin Ferreira, Aditya Jha, Andreas Carlson, David S. Dean, Yacine Amarouchene, Tak Shing Chan, Thomas Salez

Abstract: Soft lubrication has been shown to drastically affect the mobility of an object immersed in a viscous fluid in the vicinity of a purely elastic wall. In this theoretical study, we develop a minimal model incorporating viscoelasticity, carrying out a perturbation analysis in both the elastic deformation of the wall and its viscous damping. Our approach reveals the singular-perturbation nature of… ▽ More Soft lubrication has been shown to drastically affect the mobility of an object immersed in a viscous fluid in the vicinity of a purely elastic wall. In this theoretical study, we develop a minimal model incorporating viscoelasticity, carrying out a perturbation analysis in both the elastic deformation of the wall and its viscous damping. Our approach reveals the singular-perturbation nature of viscoelasticity to soft lubrication. Numerical resolution of the resulting non-linear, singular and coupled equations of motion reveals peculiar effects of viscoelasticity on confined colloidal mobility, opening the way towards the description of complex migration scenarios near realistic polymeric substrates and biological membranes. △ Less

Submitted 5 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

Uphill drift in the absence of current in single-file diffusion

Authors: Benjamin Sorkin, David S. Dean

Abstract: Single-file diffusion is a paradigmatic model for the transport of Brownian colloidal particles in narrow one-dimensional channels, such as those found in certain porous media, where the particles cannot cross each other. We consider a system where a different external uniform potential is present to the right and left of an origin. For example, this is the case when two channels meeting at the or… ▽ More Single-file diffusion is a paradigmatic model for the transport of Brownian colloidal particles in narrow one-dimensional channels, such as those found in certain porous media, where the particles cannot cross each other. We consider a system where a different external uniform potential is present to the right and left of an origin. For example, this is the case when two channels meeting at the origin have different radii. In equilibrium, the chemical potential of the particles are equal, the density is thus lower in the region with the higher potential, and by definition there is no net current in the system. Remarkably, a single-file tracer particle initially located at the origin, with position denoted by $Y(t)$, exhibits an average {\em up-hill} drift toward the region of {\em highest} potential. This drift has the late time behavior $\langle Y(t)\rangle= C t^{1/4}$, where the prefactor $C$ depends on the initial particle arrangement. This surprising result is shown analytically by computing the first two moments of $Y(t)$ through a simple and physically-illuminating method, and also via extensive numerical simulations. △ Less

Submitted 27 March, 2024; originally announced March 2024.

Comments: 17 pages, 10 figures

Self-phoretic oscillatory motion in a harmonic trap

Authors: A. Alexandre, L. Anderson, T. Collin-Dufresne, T. Guérin, D. S. Dean

Abstract: We consider the motion of a harmonically trapped overdamped particle, which is submitted to a self-phoretic force, that is proportional to the gradient of a diffusive field for which the particle itself is the source. In agreement with existing results for free particles or particles in a bounded domain, we find that the system exhibits a transition between an immobile phase, where the particle st… ▽ More We consider the motion of a harmonically trapped overdamped particle, which is submitted to a self-phoretic force, that is proportional to the gradient of a diffusive field for which the particle itself is the source. In agreement with existing results for free particles or particles in a bounded domain, we find that the system exhibits a transition between an immobile phase, where the particle stays at the center of the trap, and an oscillatory state. We perform an exact analysis giving access to the bifurcation threshold, as well as the frequency of oscillations and their amplitude near the threshold. Our analysis also characterizes the shape of two-dimensional oscillations, that take place along a circle or a straight line. Our results are confirmed by numerical simulations. △ Less

Submitted 12 March, 2024; originally announced March 2024.

Comments: 10 pages 8 figures

Single-file diffusion in spatially inhomogeneous systems

Authors: Benjamin Sorkin, David S. Dean

Abstract: We study the effect of spatially-varying potential and diffusivity on the dispersion of a tracer particle in single-file diffusion. Non-interacting particles in such a system exhibit normal diffusion at late times, which is characterised by an effective diffusion constant $D_\mathrm{eff}$. Here we demonstrate the physically appealing result that the dispersion of single-file tracers in this system… ▽ More We study the effect of spatially-varying potential and diffusivity on the dispersion of a tracer particle in single-file diffusion. Non-interacting particles in such a system exhibit normal diffusion at late times, which is characterised by an effective diffusion constant $D_\mathrm{eff}$. Here we demonstrate the physically appealing result that the dispersion of single-file tracers in this system has the same long-time behavior as that for Brownian particles in a spatially-homogeneous system with constant diffusivity $D_\mathrm{eff}$. Our results are based on a late-time analysis of the Fokker-Planck equation, motivated by the mathematical theory of homogenization. The findings are confirmed by numerical simulations for both annealed and quenched initial conditions. △ Less

Submitted 18 June, 2023; originally announced June 2023.

Comments: 14 pages, 5 figures

Interactions enhance dispersion in fluctuating channels via emergent flows

Authors: Yating Wang, David S. Dean, Sophie Marbach, Ruben Zakine

Abstract: Understanding particle motion in narrow channels is essential to guide progress in numerous applications, from filtration to vascular transport. Thermal or active fluctuations of channel walls for fluid-filled channels can slow down or increase the dispersion of tracer particles. Entropic trapping in the wall bulges slows dispersion, and hydrodynamic flows induced by wall fluctuations enhance disp… ▽ More Understanding particle motion in narrow channels is essential to guide progress in numerous applications, from filtration to vascular transport. Thermal or active fluctuations of channel walls for fluid-filled channels can slow down or increase the dispersion of tracer particles. Entropic trapping in the wall bulges slows dispersion, and hydrodynamic flows induced by wall fluctuations enhance dispersion. Previous studies primarily concentrated on the case of a single Brownian tracer either embedded in an incompressible fluid or in the ideal case where the presence of fluid is ignored. Here we address the question of what happens when there is a large ensemble of interacting Brownian tracers -- a common situation in applications. Introducing repulsive interactions between the tracer particles, while ignoring the presence of a background fluid, leads to an effective flow field. This flow field enhances tracer dispersion, a phenomenon strongly reminiscent of that seen in incompressible background fluid. We characterise the dispersion by the long-time diffusion coefficient of a single tracer, numerically and analytically with a mean-field density functional analysis. We find a surprising effect where an increased particle density enhances the diffusion coefficient, challenging the notion that crowding effects tend to reduce diffusion. Here, inter-particle interactions push particles closer to the fluctuating channel walls. Then, interactions between the fluctuating wall and the now-nearby particles drive particle mixing. Our mechanism is sufficiently general that we expect it to apply to various systems. In addition, the perturbation theory we derive quantifies dispersion in generic advection-diffusion systems with arbitrary spatiotemporal drift. △ Less

Submitted 6 May, 2023; originally announced May 2023.

Comments: 18 pages + Appendices + Supplemental Material

Effusion of stochastic processes on a line

Authors: David S. Dean, Satya N. Majumdar, Gregory Schehr

Abstract: We consider the problem of leakage or effusion of an ensemble of independent stochastic processes from a region where they are initially randomly distributed. The case of Brownian motion, initially confined to the left half line with uniform density and leaking into the positive half line is an example which has been extensively studied in the literature. Here we derive new results for the average… ▽ More We consider the problem of leakage or effusion of an ensemble of independent stochastic processes from a region where they are initially randomly distributed. The case of Brownian motion, initially confined to the left half line with uniform density and leaking into the positive half line is an example which has been extensively studied in the literature. Here we derive new results for the average number and variance of the number of leaked particles for arbitrary Gaussian processes initially confined to the negative half line and also derive its joint two-time probability distribution, both for the annealed and the quenched initial conditions. For the annealed case, we show that the two-time joint distribution is a bivariate Poisson distribution. We also discuss the role of correlations in the initial particle positions on the statistics of the number of particles on the positive half line. We show that the strong memory effects in the variance of the particle number on the positive real axis for Brownian particles, seen in recent studies, persist for arbitrary Gaussian processes and also at the level of two-time correlation functions. △ Less

Submitted 15 March, 2023; originally announced March 2023.

Comments: 19 pages, 2 figues

Journal ref: J. Stat. Mech. (2023) 063208

Non-Gaussian diffusion near surfaces

Authors: Arthur Alexandre, Maxime Lavaud, Nicolas Fares, Elodie Millan, Yann Louyer, Thomas Salez, Yacine Amarouchene, Thomas Guérin, David S. Dean

Abstract: We study the diffusion of particles confined close to a single wall and in double-wall planar channel geometries where the local diffusivities depend on the distance to the boundaries. Displacement parallel to the walls is Brownian as characterized by its variance, but it is non-Gaussian having a non-zero fourth cumulant. Establishing a link with Taylor dispersion, we calculate the fourth cumul… ▽ More We study the diffusion of particles confined close to a single wall and in double-wall planar channel geometries where the local diffusivities depend on the distance to the boundaries. Displacement parallel to the walls is Brownian as characterized by its variance, but it is non-Gaussian having a non-zero fourth cumulant. Establishing a link with Taylor dispersion, we calculate the fourth cumulant and the tails of the displacement distribution for general diffusivity tensors along with potentials generated by either the walls or externally, for instance gravity. Experimental and numerical studies of the motion of a colloid in the direction parallel to the wall give measured fourth cumulants which are correctly predicted by our theory. Interestingly, contrary to models of Brownian-yet-non-Gaussian diffusion, the tails of the displacement distribution are shown to be Gaussian rather than exponential. All together, our results provide additional tests and constraints for the inference of force maps and local transport properties near surfaces. △ Less

Submitted 19 December, 2022; v1 submitted 5 July, 2022; originally announced July 2022.

Nanoparticle Taylor dispersion near charged surfaces with an open boundary

Authors: Alexandre Vilquin, Vincent Bertin, Elie Raphaël, David S Dean, Thomas Salez, Joshua D. McGraw

Abstract: The dispersive spreading of microscopic particles in shear flows is influenced both by advection and thermal motion. At the nanoscale, interactions between such particles and their confining boundaries become unavoidable. We address the roles of electrostatic repulsion and absorption on the spatial distribution and dispersion of charged nanoparticles in near-surface shear flows, observed under eva… ▽ More The dispersive spreading of microscopic particles in shear flows is influenced both by advection and thermal motion. At the nanoscale, interactions between such particles and their confining boundaries become unavoidable. We address the roles of electrostatic repulsion and absorption on the spatial distribution and dispersion of charged nanoparticles in near-surface shear flows, observed under evanescent illumination. The electrostatic repulsion between particles and the lower charged surface is tuned by varying electrolyte concentrations. Particles leaving the field of vision can be neglected from further analysis, such that the experimental ensemble is equivalent to that of Taylor dispersion with absorption. These two ingredients modify the particle distribution, deviating strongly from the Gibbs-Boltzmann one at the nanoscale studied here. The overall effect is to restrain the accessible space available to particles, leading to a striking, ten-fold reduction in the spreading dynamics as compared to the non-interacting case. △ Less

Submitted 6 January, 2023; v1 submitted 15 June, 2022; originally announced June 2022.

How stickiness can speed up diffusion in confined systems

Authors: Arthur Alexandre, Matthieu Mangeat, Thomas Guérin, David S. Dean

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

Submitted 2 June, 2022; v1 submitted 10 December, 2021; originally announced December 2021.

Effective diffusivity of Brownian particles in a two dimensional square lattice of hard disks

Authors: M. Mangeat, T. Guérin, D. S. Dean

Abstract: We revisit the classic problem of the effective diffusion constant of a Brownian particle in a square lattice of reflecting impenetrable hard disks. This diffusion constant is also related to the effective conductivity of non-conducting and infinitely conductive disks in the same geometry. We show how a recently derived Green's function for the periodic lattice can be exploited to derive a series… ▽ More We revisit the classic problem of the effective diffusion constant of a Brownian particle in a square lattice of reflecting impenetrable hard disks. This diffusion constant is also related to the effective conductivity of non-conducting and infinitely conductive disks in the same geometry. We show how a recently derived Green's function for the periodic lattice can be exploited to derive a series expansion of the diffusion constant in terms of the disk's volume fraction $\varphi$. Secondly we propose a variant of the Fick-Jacobs approximation to study the large volume fraction limit. This combination of analytical results is shown to describe the behavior of the diffusion constant for all volume fractions. △ Less

Submitted 8 November, 2021; originally announced November 2021.

Journal ref: J. Chem. Phys. 152, 234109 (2020)

Steady state of overdamped particles in the non-conservative force field of a simple non-linear model of optical trap

Authors: M. Mangeat, T. Guérin, D. S. Dean

Abstract: Optically trapped particles are often subject to a non-conservative scattering force arising from radiation pressure. In this paper we present an exact solution for the steady state statistics of an overdamped Brownian particle subjected to a commonly used force field model for an optical trap. The model is the simplest of its kind that takes into account non-conservative forces. In particular, we… ▽ More Optically trapped particles are often subject to a non-conservative scattering force arising from radiation pressure. In this paper we present an exact solution for the steady state statistics of an overdamped Brownian particle subjected to a commonly used force field model for an optical trap. The model is the simplest of its kind that takes into account non-conservative forces. In particular, we present exact results for certain marginals of the full three dimensional steady state probability distribution as well as results for the toroidal probability currents which are present in the steady state, as well as for the circulation of theses currents. Our analytical results are confirmed by numerical solution of the steady state Fokker-Planck equation. △ Less

Submitted 8 October, 2021; originally announced October 2021.

Comments: 14 pages, 3 figures

Journal ref: Journal of Statistical Mechanics: Theory and Experiment (2021) 113205

Generalized Taylor dispersion for translationally invariant microfluidic systems

Authors: Arthur Alexandre, Thomas Guérin, David S. Dean

Abstract: We consider Taylor dispersion for tracer particles in micro-fluidic planar channels with strong confinement. In this context, the channel walls modify the local diffusivity tensor and also interactions between the tracer particles and the walls become important. We provide a simple and general formula for the effective diffusion constant along the channel as well as the first non-trivial finite ti… ▽ More We consider Taylor dispersion for tracer particles in micro-fluidic planar channels with strong confinement. In this context, the channel walls modify the local diffusivity tensor and also interactions between the tracer particles and the walls become important. We provide a simple and general formula for the effective diffusion constant along the channel as well as the first non-trivial finite time correction for arbitrary flows along the channel, arbitrary interaction potentials with the walls and arbitrary expressions for the diffusion tensor. The formula are in particular amenable to a straightforward numerical implementation, rendering them extremely useful for comparison with experiments. We present a number of applications, notably for systems which have parabolically varying diffusivity profiles, to systems with attractive interactions with the walls as well as electroosmotic flows between plates with differing surface charges within the Debye-Hückel approximation. △ Less

Submitted 11 August, 2021; v1 submitted 13 May, 2021; originally announced May 2021.

Comments: 33 pages, 7 figures

Journal ref: Physics of Fluids 33, 082004 (2021)

Wigner function for noninteracting fermions in hard wall potentials

Authors: Benjamin De Bruyne, David S. Dean, Pierre Le Doussal, Satya N. Majumdar, Gregory Schehr

Abstract: The Wigner function $W_N({\bf x}, {\bf p})$ is a useful quantity to characterize the quantum fluctuations of an $N$-body system in its phase space. Here we study $W_N({\bf x}, {\bf p})$ for $N$ noninteracting spinless fermions in a $d$-dimensional spherical hard box of radius $R$ at temperature $T=0$. In the large $N$ limit, the local density approximation (LDA) predicts that… ▽ More The Wigner function $W_N({\bf x}, {\bf p})$ is a useful quantity to characterize the quantum fluctuations of an $N$-body system in its phase space. Here we study $W_N({\bf x}, {\bf p})$ for $N$ noninteracting spinless fermions in a $d$-dimensional spherical hard box of radius $R$ at temperature $T=0$. In the large $N$ limit, the local density approximation (LDA) predicts that $W_N({\bf x}, {\bf p}) \approx 1/(2 π\hbar)^d$ inside a finite region of the $({\bf x}, {\bf p})$ plane, namely for $|{\bf x}| < R$ and $|{\bf p}| < k_F$ where $k_F$ is the Fermi momentum, while $W_N({\bf x}, {\bf p})$ vanishes outside this region, or "droplet", on a scale determined by quantum fluctuations. In this paper we investigate systematically, in this quantum region, the structure of the Wigner function along the edge of this droplet, called the Fermi surf. In one dimension, we find that there are three distinct edge regions along the Fermi surf and we compute exactly the associated nontrivial scaling functions in each regime. We also study the momentum distribution $\hat ρ_N(p)$ and find a striking algebraic tail for very large momenta $\hat ρ_N(p) \propto 1/p^4$, well beyond $k_F$, reminiscent of a similar tail found in interacting quantum systems (discussed in the context of Tan's relation). We then generalize these results to higher $d$ and find, remarkably, that the scaling function close to the edge of the box is universal, i.e., independent of the dimension~$d$. △ Less

Submitted 30 April, 2021; v1 submitted 11 April, 2021; originally announced April 2021.

Comments: 31 pages, 14 figures

Journal ref: Phys. Rev. A 104, 013314 (2021)

Impurities in systems of noninteracting trapped fermions

Authors: David S. Dean, Pierre Le Doussal, Satya N. Majumdar, Gregory Schehr

Abstract: We study the properties of spin-less non-interacting fermions trapped in a confining potential in one dimension but in the presence of one or more impurities which are modelled by delta function potentials. We use a method based on the single particle Green's function. For a single impurity placed in the bulk, we compute the density of the Fermi gas near the impurity. Our results, in addition to r… ▽ More We study the properties of spin-less non-interacting fermions trapped in a confining potential in one dimension but in the presence of one or more impurities which are modelled by delta function potentials. We use a method based on the single particle Green's function. For a single impurity placed in the bulk, we compute the density of the Fermi gas near the impurity. Our results, in addition to recovering the Friedel oscillations at large distance from the impurity, allow the exact computation of the density at short distances. We also show how the density of the Fermi gas is modified when the impurity is placed near the edge of the trap in the region where the unperturbed system is described by the Airy gas. Our method also allows us to compute the effective potential felt by the impurity both in the bulk and at the edge. In the bulk this effective potential is shown to be a universal function only of the local Fermi wave vector, or equivalently of the local fermion density. When the impurity is placed near the edge of the Fermi gas, the effective potential can be expressed in terms of Airy functions. For an attractive impurity placed far outside the support of the fermion density, we show that an interesting transition occurs where a single fermion is pulled out of the Fermi sea and forms a bound state with the impurity. This is a quantum analogue of the well-known Baik-Ben Arous-Péché (BBP) transition, known in the theory of spiked random matrices. The density at the location of the impurity plays the role of an order parameter. We also consider the case of two impurities in the bulk and compute exactly the effective force between them mediated by the background Fermi gas. △ Less

Submitted 24 February, 2021; v1 submitted 27 December, 2020; originally announced December 2020.

Comments: 35 pages, 8 figures. Revised version: references added and typos corrected

Journal ref: SciPost Phys. 10, 082 (2021)

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

Authors: David S. Dean, Pierre Le Doussal, Satya N. Majumdar, Gregory Schehr, Naftali R. Smith

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

Submitted 27 September, 2020; originally announced September 2020.

Comments: 35 pages, 9 figures

Journal ref: J. Phys. A: Math. Theor. 54 084001 (2021)

Position distribution in a generalised run and tumble process

Authors: David S. Dean, Satya N. Majumdar, Hendrik Schawe

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

Submitted 3 September, 2020; originally announced September 2020.

Comments: 24 pages, 4 Figures

Journal ref: Phys. Rev. E 103, 012130 (2021)

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

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

Authors: Naftali R. Smith, David S. Dean, Pierre Le Doussal, Satya N. Majumdar, Gregory Schehr

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

Submitted 17 May, 2020; v1 submitted 31 January, 2020; originally announced January 2020.

Comments: 16 pages, 8 figures

Journal ref: Phys. Rev. A 101, 053602 (2020)

The effect of driving on model C interfaces

Authors: David S. Dean, Paul Gersberg, Peter C. W. Holdsworth

Abstract: We consider the effect of uniform driving on the interface between two phases which are described by model C dynamics. The non-driven system has a classical Gaussian interface described by capillary wave theory. The model under driving retains Gaussian statistics but the interface statistics are modified by driving, notably the height fluctuations are suppressed and the correlation length of the f… ▽ More We consider the effect of uniform driving on the interface between two phases which are described by model C dynamics. The non-driven system has a classical Gaussian interface described by capillary wave theory. The model under driving retains Gaussian statistics but the interface statistics are modified by driving, notably the height fluctuations are suppressed and the correlation length of the fluctuations is increased. The model we introduce can also be used as a model for the effect of activity on interface dynamics. △ Less

Submitted 13 January, 2020; originally announced January 2020.

Comments: 10 pages, Revtex

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

Nonequilibrium dynamics of noninteracting fermions in a trap

Authors: David S. Dean, Pierre Le Doussal, Satya N. Majumdar, Gregory Schehr

Abstract: We consider the real time dynamics of $N$ noninteracting fermions in $d=1$. They evolve in a trapping potential $V(x)$, starting from the equilibrium state in a potential $V_0(x)$. We study the time evolution of the Wigner function $W(x,p,t)$ in the phase space $(x,p)$, and the associated kernel which encodes all correlation functions. At $t=0$ the Wigner function for large $N$ is uniform in phase… ▽ More We consider the real time dynamics of $N$ noninteracting fermions in $d=1$. They evolve in a trapping potential $V(x)$, starting from the equilibrium state in a potential $V_0(x)$. We study the time evolution of the Wigner function $W(x,p,t)$ in the phase space $(x,p)$, and the associated kernel which encodes all correlation functions. At $t=0$ the Wigner function for large $N$ is uniform in phase space inside the Fermi volume, and vanishes at the Fermi surf over a scale $e_N$ being described by a universal scaling function related to the Airy function. We obtain exact solutions for the Wigner function, the density, and the correlations in the case of harmonic and inverse square potentials, for several $V_0(x)$. In the large $N$ limit, near the edges where the density vanishes, we obtain limiting kernels (of the Airy or Bessel types) that retain the form found in equilibrium, up to a time dependent rescaling. For non-harmonic traps the evolution of the Fermi volume is more complex. Nevertheless we show that, for intermediate times, the Fermi surf is still described by the same equilibrium scaling function, with a non-trivial time and space dependent width which we compute analytically. We discuss the multi-time correlations and obtain their explicit scaling forms valid near the edge for the harmonic oscillator. Finally, we address the large time limit where relaxation to the Generalized Gibbs Ensemble (GGE) was found to occur in the "classical" regime $\hbar \sim 1/N$. Using the diagonal ensemble we compute the Wigner function in the quantum case (large $N$, fixed $\hbar$) and show that it agrees with the GGE. We also obtain the higher order (non-local) correlations in the diagonal ensemble. △ Less

Submitted 18 March, 2019; v1 submitted 7 February, 2019; originally announced February 2019.

Comments: Main text: 8 pages, 1 figure + Supplementary Material: 28 pages, 5 figures

Journal ref: EPL, 126 (2019) 20006

The role of non-conservative scattering forces and damping on Brownian particles in optical traps

Authors: Matthieu Mangeat, Yacine Amarouchene, Yann Louyer, Thomas Guérin, David S. Dean

Abstract: We consider a model of a particle trapped in a harmonic optical trap but with the addition of a non-conservative radiation induced force. This model is known to correctly describe experimentally observed trapped particle statistics for a wide range of physical parameters such as temperature and pressure. We theoretically analyse the effect of non-conservative force on the underlying steady state d… ▽ More We consider a model of a particle trapped in a harmonic optical trap but with the addition of a non-conservative radiation induced force. This model is known to correctly describe experimentally observed trapped particle statistics for a wide range of physical parameters such as temperature and pressure. We theoretically analyse the effect of non-conservative force on the underlying steady state distribution as well as the power spectrum for the particle position. We compute perturbatively the probability distribution of the resulting non-equilibrium steady states for all dynamical regimes, underdamped through to overdamped and give expressions for the associated currents in phase space (position and velocity). We also give the spectral density of the trapped particle's position in all dynamical regimes and for any value of the non-conservative force. Signatures of the presence non-conservative forces are shown to be particularly strong for in the underdamped regime at low frequencies. △ Less

Submitted 21 December, 2018; originally announced December 2018.

Comments: 25 pages, 5 figures

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

Nonequilibrium dynamics induced by scattering forces for optically trapped nanoparticles in strongly inertial regimes

Authors: Yacine Amarouchene, Matthieu Mangeat, Benjamin Vidal Montes, Lukas Ondic, Thomas Guérin, David. S. Dean, Yann Louyer

Abstract: As a first approximation, the forces acting on optically trapped particles are commonly assumed to be conservative. The influence of the nonconservative force has been shown to be negligible in overdamped liquid environments. However, its impact in the underdamped regime remains unexplored. Here, we experimentally study the combined effects of gradient and scattering forces on the nonlinear dynami… ▽ More As a first approximation, the forces acting on optically trapped particles are commonly assumed to be conservative. The influence of the nonconservative force has been shown to be negligible in overdamped liquid environments. However, its impact in the underdamped regime remains unexplored. Here, we experimentally study the combined effects of gradient and scattering forces on the nonlinear dynamics of trapped particles at various air pressures, leading to significant low-frequency noise excess along the optical axis and the emergence of toroidal Brownian vortices. Our results supported, by numerical simulations and a theoretical model, provide fundamental insights into the functioning of optical tweezers and a means for investigating nonequilibrium steady states in the presence of nonconservative forces. △ Less

Submitted 9 April, 2019; v1 submitted 14 December, 2018; originally announced December 2018.

Journal ref: Phys. Rev. Lett. 122, 183901 (2019)

Noninteracting fermions in a trap and random matrix theory

Authors: David S. Dean, Pierre Le Doussal, Satya N. Majumdar, Gregory Schehr

Abstract: We review recent advances in the theory of trapped fermions using techniques borrowed from random matrix theory (RMT) and, more generally, from the theory of determinantal point processes. In the presence of a trap, and in the limit of a large number of fermions $N \gg 1$, the spatial density exhibits an edge, beyond which it vanishes. While the spatial correlations far from the edge, i.e. close t… ▽ More We review recent advances in the theory of trapped fermions using techniques borrowed from random matrix theory (RMT) and, more generally, from the theory of determinantal point processes. In the presence of a trap, and in the limit of a large number of fermions $N \gg 1$, the spatial density exhibits an edge, beyond which it vanishes. While the spatial correlations far from the edge, i.e. close to the center of the trap, are well described by standard many-body techniques, such as the local density approximation (LDA), these methods fail to describe the fluctuations close to the edge of the Fermi gas, where the density is very small and the fluctuations are thus enhanced. It turns out that RMT and determinantal point processes offer a powerful toolbox to study these edge properties in great detail. Here we discuss the principal edge universality classes, that have been recently identified using these modern tools. In dimension $d=1$ and at zero temperature $T=0$, these universality classes are in one-to-one correspondence with the standard universality classes found in the classical unitary random matrix ensembles: soft edge (described by the "Airy kernel") and hard edge (described by the "Bessel kernel") universality classes. We further discuss extensions of these results to higher dimensions $d\geq 2$ and to finite temperature. Finally, we discuss correlations in the phase space, i.e., in the space of positions and momenta, characterized by the so called Wigner function. △ Less

Submitted 30 October, 2018; originally announced October 2018.

Comments: 32 pages, 5 figures, contribution to a special issue of Journal of Physics A on Random Matrices: the first 90 years

Journal ref: J. Phys. A: Math. Theor. 52 144006 (2019)

Dispersion in two-dimensional periodic channels with discontinuous profiles

Authors: M. Mangeat, T. Guérin, D. S. Dean

Abstract: The effective diffusivity of Brownian tracer particles confined in periodic micro-channels is smaller than the microscopic diffusivity due to entropic trapping. Here, we study diffusion in two-dimensional periodic channels whose cross-section presents singular points, such as abrupt changes of radius or the presence of thin walls, with openings, delimiting periodic compartments composing the chann… ▽ More The effective diffusivity of Brownian tracer particles confined in periodic micro-channels is smaller than the microscopic diffusivity due to entropic trapping. Here, we study diffusion in two-dimensional periodic channels whose cross-section presents singular points, such as abrupt changes of radius or the presence of thin walls, with openings, delimiting periodic compartments composing the channel. Dispersion in such systems is analyzed using the Fick-Jacobs' approximation. This approximation assumes a much faster equilibration in the lateral than in the axial direction, along which the dispersion is measured. If the characteristic width $a$ of the channel is much smaller than the period $L$ of the channel, i.e. $\varepsilon = a/L$ is small, this assumption is clearly valid for Brownian particles. For discontinuous channels, the Fick-Jacobs' approximation is only valid at the lowest order in $\varepsilon$ and provides a rough, though on occasions rather accurate, estimate of the effective diffusivity. Here we provide formulas for the effective diffusivity in discontinuous channels that are asymptotically exact at the next-to-leading order in $\varepsilon$. Each discontinuity leads to a reduction of the effective diffusivity. We show that our theory is consistent with the picture of effective {\em trapping rates} associated with each discontinuity, for which our theory provides explicit and asymptotically exact formulas. Our analytical predictions are confirmed by numerical analysis. Our results provide a precise quantification of the kinetic entropic barriers associated with profile singularities. △ Less

Submitted 14 July, 2018; originally announced July 2018.

Comments: 13 pages, 7 figures

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

Wigner function of noninteracting trapped fermions

Authors: David S. Dean, P. Le Doussal, Satya N. Majumdar, Gregory Schehr

Abstract: We study analytically the Wigner function $W_N({\bf x},{\bf p})$ of $N$ noninteracting fermions trapped in a smooth confining potential $V({\bf x})$ in $d$ dimensions. At zero temperature, $W_N({\bf x},{\bf p})$ is constant over a finite support in the phase space $({\bf x},{\bf p})$ and vanishes outside. Near the edge of this support, we find a universal scaling behavior of… ▽ More We study analytically the Wigner function $W_N({\bf x},{\bf p})$ of $N$ noninteracting fermions trapped in a smooth confining potential $V({\bf x})$ in $d$ dimensions. At zero temperature, $W_N({\bf x},{\bf p})$ is constant over a finite support in the phase space $({\bf x},{\bf p})$ and vanishes outside. Near the edge of this support, we find a universal scaling behavior of $W_N({\bf x},{\bf p})$ for large $N$. The associated scaling function is independent of the precise shape of the potential as well as the spatial dimension $d$. We further generalize our results to finite temperature $T>0$. We show that there exists a low temperature regime $T \sim e_N/b$ where $e_N$ is an energy scale that depends on $N$ and the confining potential $V({\bf x})$, where the Wigner function at the edge again takes a universal scaling form with a $b$-dependent scaling function. This temperature dependent scaling function is also independent of the potential as well as the dimension $d$. Our results generalize to any $d\geq 1$ and $T \geq 0$ the $d=1$ and $T=0$ results obtained by Bettelheim and Wiegman [Phys. Rev. B ${\bf 84}$, 085102 (2011)]. △ Less

Submitted 8 January, 2018; originally announced January 2018.

Comments: 16 pages, 4 figures

Journal ref: Phys. Rev. A 97, 063614 (2018)

Stresses in non-equilibrium fluids: Exact formulation and coarse grained theory

Authors: M. Krüger, A. Solon, V. Démery, C. M. Rohwer, D. S. Dean

Abstract: Starting from the stochastic equation for the density operator, we formulate the exact (instantaneous) stress tensor for interacting Brownian particles, whose average value agrees with expressions derived previously. We analyze the relation between the stress tensor and forces on external potentials, and observe that, out of equilibrium, particle currents give rise to extra forces. Next, we derive… ▽ More Starting from the stochastic equation for the density operator, we formulate the exact (instantaneous) stress tensor for interacting Brownian particles, whose average value agrees with expressions derived previously. We analyze the relation between the stress tensor and forces on external potentials, and observe that, out of equilibrium, particle currents give rise to extra forces. Next, we derive the stress tensor for a Landau-Ginzburg theory in non-equilibrium situations, finding an expression analogous to that of the exact microscopic stress tensor, and discuss the computation of out-of-equilibrium (classical) Casimir forces. We use these relations to study the spatio-temporal correlations of the stress tensor in a Brownian fluid, which we derive exactly to leading order in the interaction potential strength. We observe that, after integration over time, the spatial correlations generally decay as power laws in space. These are expected to be of importance for driven confined systems. We also show that divergence-free parts of the stress tensor do not contribute in the Green-Kubo relation for the viscosity. △ Less

Submitted 14 December, 2017; originally announced December 2017.

Comments: 12 pages

Dispersion in two dimensional channels - the Fick-Jacobs approximation revisited

Authors: M. Mangeat, T. Guérin, D. S. Dean

Abstract: We examine the dispersion of Brownian particles in a symmetric two dimensional channel, this classical problem has been widely studied in the literature using the so called Fick-Jacobs' approximation and its various improvements. Most studies rely on the reduction to an effective one dimensional diffusion equation, here we drive an explicit formula for the diffusion constant which avoids this redu… ▽ More We examine the dispersion of Brownian particles in a symmetric two dimensional channel, this classical problem has been widely studied in the literature using the so called Fick-Jacobs' approximation and its various improvements. Most studies rely on the reduction to an effective one dimensional diffusion equation, here we drive an explicit formula for the diffusion constant which avoids this reduction. Using this formula the effective diffusion constant can be evaluated numerically without resorting to Brownian simulations. In addition a perturbation theory can be developed in $\varepsilon = h_0/L$ where $h_0$ is the characteristic channel height and $L$ the period. This perturbation theory confirms the results of Kalinay and Percus (Phys. Rev. E 74, 041203 (2006)), based on the reduction, to one dimensional diffusion are exact at least to ${\cal O}(\varepsilon^6)$. Furthermore, we show how the Kalinay and Percus pseudo-linear approximation can be straightforwardly recovered. The approach proposed here can also be exploited to yield exact results an appropriate limit $\varepsilon \to \infty$, we show that here the diffusion constant remains finite and show how the result can be obtained with a simple physical argument. Moreover we show that the correction to the effective diffusion constant is of order $1/\varepsilon$ and remarkably has a some universal characteristics. Numerically we compare the analytical results obtained with exact numerical calculations for a number of interesting channel geometries. △ Less

Submitted 8 November, 2017; v1 submitted 7 October, 2017; originally announced October 2017.

Comments: 21 pages, 5 figures

Journal ref: J. Stat. Mech. (2017) 123205

Geometry controlled dispersion in periodic corrugated channels

Authors: M. Mangeat, T. Guérin, D. S. Dean

Abstract: The effective diffusivity $D_e$ of tracer particles diffusing in periodically corrugated axisymmetric two and three dimensional channels is studied. The majority of previous studies of this class of problems are based on perturbative analyses about narrow channels, where the problem can be reduced to an effectively one dimensional one. Here we show how to analyze this class of problems using a muc… ▽ More The effective diffusivity $D_e$ of tracer particles diffusing in periodically corrugated axisymmetric two and three dimensional channels is studied. The majority of previous studies of this class of problems are based on perturbative analyses about narrow channels, where the problem can be reduced to an effectively one dimensional one. Here we show how to analyze this class of problems using a much more general approach which even includes the limit of infinitely wide channels. Using the narrow and wide channel asymptotics, we provide a Pad{é} approximant scheme that is able to describe the dispersion properties of a wide class of channels. Furthermore, we systematically identify all the exact asymptotic scaling regimes of $D_e$ and the accompanying physical mechanisms that control dispersion, clarifying the distinction between smooth channels and compartmentalized ones, and identifying the regimes in which $D_e$ can be linked to first passage problems. △ Less

Submitted 12 September, 2017; originally announced September 2017.

Comments: Main text (6 pages) + Supplementary Information (5 pages)

Journal ref: EPL, 118, 40004 (2017)

A Gaussian theory for fluctuations in simple liquids

Authors: Matthias Krüger, David S. Dean

Abstract: Assuming an effective quadratic Hamiltonian, we derive an approximate, linear stochastic equation of motion for the density-fluctuations in liquids, composed of overdamped Brownian particles. From this approach, time dependent two point correlation functions (such as the intermediate scattering function) are derived. We show that this correlation function is exact at short times, for any interacti… ▽ More Assuming an effective quadratic Hamiltonian, we derive an approximate, linear stochastic equation of motion for the density-fluctuations in liquids, composed of overdamped Brownian particles. From this approach, time dependent two point correlation functions (such as the intermediate scattering function) are derived. We show that this correlation function is exact at short times, for any interaction and, in particular, for arbitrary external potentials so that it applies to confined systems. Furthermore, we discuss the relation of this approach to previous ones, such as dynamical density functional theory as well as the formally exact treatment. This approach, inspired by the well known Landau-Ginzburg Hamiltonians, and the corresponding "Model B" equation of motion, may be seen as its microscopic version, containing information about the details on the particle level. △ Less

Submitted 26 January, 2017; originally announced January 2017.

Comments: 8 pages

Universal time dependent dispersion properties for diffusion in a one-dimensional critically tilted potential

Authors: T. Guérin, D. S. Dean

Abstract: We consider the time dependent dispersion properties of overdamped tracer particles diffusing in a one dimensional periodic potential under the influence of an additional constant tilting force $F$. The system is studied in the region where the force is close to the critical value $F_c$ at which the barriers separating neighboring potential wells disappear. We show that, when $F$ crosses the criti… ▽ More We consider the time dependent dispersion properties of overdamped tracer particles diffusing in a one dimensional periodic potential under the influence of an additional constant tilting force $F$. The system is studied in the region where the force is close to the critical value $F_c$ at which the barriers separating neighboring potential wells disappear. We show that, when $F$ crosses the critical value, the shape of the Mean-Square Displacement (MSD) curves is strongly modified. We identify a diffusive regime at intermediate time scales, with an effective diffusion coefficient which is much larger than the late time diffusion coefficient for $F>F_c$, whereas for $F<F_c$ the late time and intermediate time diffusive regimes are indistinguishable. Explicit asymptotic regimes for the MSD curves are identified at all time scales. △ Less

Submitted 9 January, 2017; originally announced January 2017.

Journal ref: Phys. Rev. E, 95, 012109 (2017)

Statistics of the maximal distance and momentum in a trapped Fermi gas at low temperature

Authors: David S. Dean, Pierre Le Doussal, Satya N. Majumdar, Gregory Schehr

Abstract: We consider $N$ non-interacting fermions in an isotropic $d$-dimensional harmonic trap. We compute analytically the cumulative distribution of the maximal radial distance of the fermions from the trap center at zero temperature. While in $d=1$ the limiting distribution (in the large $N$ limit), properly centered and scaled, converges to the squared Tracy-Widom distribution of the Gaussian Unitary… ▽ More We consider $N$ non-interacting fermions in an isotropic $d$-dimensional harmonic trap. We compute analytically the cumulative distribution of the maximal radial distance of the fermions from the trap center at zero temperature. While in $d=1$ the limiting distribution (in the large $N$ limit), properly centered and scaled, converges to the squared Tracy-Widom distribution of the Gaussian Unitary Ensemble in Random Matrix Theory, we show that for all $d>1$, the limiting distribution converges to the Gumbel law. These limiting forms turn out to be universal, i.e., independent of the details of the trapping potential for a large class of isotropic trapping potentials. We also study the position of the right-most fermion in a given direction in $d$ dimensions and, in the case of a harmonic trap, the maximum momentum, and show that they obey similar Gumbel statistics. Finally, we generalize these results to low but finite temperature. △ Less

Submitted 12 December, 2016; originally announced December 2016.

Comments: 28 pages, 4 figures

Journal ref: J. Stat. Mech. (2017) 063301

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

Authors: David S. Dean, Pierre Le Doussal, Satya N. Majumdar, Gregory Schehr

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

Submitted 14 September, 2016; originally announced September 2016.

Comments: 55 pages, 8 figures

Journal ref: Phys. Rev. A 94, 063622 (2016)

Sample-to-sample fluctuations of power spectrum of a random motion in a periodic Sinai model

Authors: David S. Dean, Antonio Iorio, Enzo Marinari, Gleb Oshanin

Abstract: The Sinai model of a tracer diffusing in a quenched Brownian potential is a much studied problem exhibiting a logarithmically slow anomalous diffusion due to the growth of energy barriers with the system size. However, if the potential is random but periodic, the regime of anomalous diffusion crosses over to one of normal diffusion once a tracer has diffused over a few periods of the system. Here… ▽ More The Sinai model of a tracer diffusing in a quenched Brownian potential is a much studied problem exhibiting a logarithmically slow anomalous diffusion due to the growth of energy barriers with the system size. However, if the potential is random but periodic, the regime of anomalous diffusion crosses over to one of normal diffusion once a tracer has diffused over a few periods of the system. Here we consider a system in which the potential is given by a Brownian Bridge on a finite interval $(0,L)$ and then periodically repeated over the whole real line, and study the power spectrum $S(f)$ of the diffusive process $x(t)$ in such a potential. We show that for most of realizations of $x(t)$ in a given realization of the potential, the low-frequency behavior is $S(f) \sim {\cal A}/f^2$, i.e., the same as for standard Brownian motion, and the amplitude ${\cal A}$ is a disorder-dependent random variable with a finite support. Focusing on the statistical properties of this random variable, we determine the moments of ${\cal A}$ of arbitrary, negative or positive order $k$, and demonstrate that they exhibit a multi-fractal dependence on $k$, and a rather unusual dependence on the temperature and on the periodicity $L$, which are supported by atypical realizations of the periodic disorder. We finally show that the distribution of ${\cal A}$ has a log-normal left tail, and exhibits an essential singularity close to the right edge of the support, which is related to the Lifshitz singularity. Our findings are based both on analytic results and on extensive numerical simulations of the process $x(t)$. △ Less

Submitted 30 June, 2016; originally announced July 2016.

Comments: 8 pages, 5 figures

Journal ref: Phys. Rev. E 94, 032131 (2016)

Nonequilibrium Tuning of the Thermal Casimir Effect

Authors: David S. Dean, Bing-Sui Lu, A. C. Maggs, Rudolf Podgornik

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

Submitted 15 June, 2016; originally announced June 2016.

Comments: 7 pages (including Supplementary Material), 2 figures

Journal ref: Phys. Rev. Lett. 116, 240602 (2016)

Phase diagram of a bulk 1d lattice Coulomb gas

Authors: V. Démery, R. Monsarrat, D. S. Dean, R. Podgornik

Abstract: The exact solution, via transfer matrix, of the simple one dimensional lattice Coulomb gas (1d LCG) model can reproduce peculiar features of ionic liquid capacitors, such as overscreening, layering, and camel- and bell-shaped capacitance curves. Using the same transfer matrix method, we now compute the bulk properties of the 1d LCG in the constant voltage ensemble. We unveil a phase diagram with r… ▽ More The exact solution, via transfer matrix, of the simple one dimensional lattice Coulomb gas (1d LCG) model can reproduce peculiar features of ionic liquid capacitors, such as overscreening, layering, and camel- and bell-shaped capacitance curves. Using the same transfer matrix method, we now compute the bulk properties of the 1d LCG in the constant voltage ensemble. We unveil a phase diagram with rich structure exhibiting a low density disordered and high density ordered phases, separated by a first order phase transition at low temperature; the solid state at full packing can be ordered or not, depending on the temperature. This phase diagram, which is strikingly similar to its three dimensional counterpart, also sheds light on the behaviour of the confined system. △ Less

Submitted 23 November, 2015; originally announced November 2015.

Comments: 6 pages, 6 figures

Journal ref: EPL 113 18008 (2016)

The conductivity of strong electrolytes from stochastic density functional theory

Authors: Vincent Démery, David S Dean

Abstract: Stochastic density functional theory is applied to analyze the conductivity of strong two species electrolytes at arbitrary field strengths. The corresponding stochastic equations for the density of the electrolyte species are solved by linearizing them about the mean density of ionic species, yielding an effective Gaussian theory. The non-equilibrium density-density correlation functions are comp… ▽ More Stochastic density functional theory is applied to analyze the conductivity of strong two species electrolytes at arbitrary field strengths. The corresponding stochastic equations for the density of the electrolyte species are solved by linearizing them about the mean density of ionic species, yielding an effective Gaussian theory. The non-equilibrium density-density correlation functions are computed and the conductivity of the electrolyte is deduced. In the bulk, our results give a simple derivation of the results of Onsager and coworkers who used very different methods. The method developed here can also be used to study electrolytes confined in one and two dimensions and interacting via either the three dimensional Coulomb interaction or the Coulomb interaction corresponding to that dimension of space. △ Less

Submitted 20 October, 2015; originally announced October 2015.

Comments: 16 pages, 3 figures

Journal ref: J. Stat. Mech. (2016) 023106

Kubo formulas for dispersion in heterogeneous periodic non-equilibrium systems

Authors: T. Guérin, D. S. Dean

Abstract: We consider the dispersion properties of tracer particles moving in non-equilibrium heterogeneous periodic media. The tracer motion is described by a Fokker-Planck equation with arbitrary spatially periodic (but constant in time) local diffusion tensors and drift, eventually with the presence of obstacles. We derive a Kubo-like formula for the time dependent effective diffusion tensor valid in any… ▽ More We consider the dispersion properties of tracer particles moving in non-equilibrium heterogeneous periodic media. The tracer motion is described by a Fokker-Planck equation with arbitrary spatially periodic (but constant in time) local diffusion tensors and drift, eventually with the presence of obstacles. We derive a Kubo-like formula for the time dependent effective diffusion tensor valid in any dimension. From this general formula, we derive expressions for the late time effective diffusion tensor and drift in these systems. In addition, we find an explicit formula for the late finite time corrections to these transport coefficients. In one dimension, we give closed analytical formula for the transport coefficients. The formulas derived here are very general and provide a straightforward method to compute the dispersion properties in arbitrary non-equilibrium periodic advection-diffusion systems. △ Less

Submitted 9 September, 2015; originally announced September 2015.

Journal ref: Phys. Rev. E 92, 062103 (2015)

Out of equilibrium thermal Casimir effect between Brownian conducting plates

Authors: B. -S. Lu, D. S. Dean, R. Podgornik

Abstract: We study the thermal fluctuation induced interactions between two surfaces containing Brownian charges which are held at different temperatures. Using a dynamical form of Debye-Hückel theory implemented within the stochastic equation for the density of mobile Brownian charges, we derive expressions for the average force between the two surfaces as well as its variance. The latter is found both for… ▽ More We study the thermal fluctuation induced interactions between two surfaces containing Brownian charges which are held at different temperatures. Using a dynamical form of Debye-Hückel theory implemented within the stochastic equation for the density of mobile Brownian charges, we derive expressions for the average force between the two surfaces as well as its variance. The latter is found both for the normal, of finite mean, as well as the lateral force, of zero mean, between the surfaces. △ Less

Submitted 10 December, 2015; v1 submitted 27 August, 2015; originally announced August 2015.

Comments: 6 pages, 3 figures, EPL style

Journal ref: EPL 112 (2015) 20001

Force-induced dispersion in heterogeneous media

Authors: Thomas Guérin, David S Dean

Abstract: The effect of a constant applied external force, induced for instance by an electric or gravitational field, on the dispersion of Brownian particles in periodic media with spatially varying diffusivity, and thus mobility, is studied. We show that external forces can greatly enhance dispersion in the direction of the applied force and also modify, to a lesser extent and in some cases non-monotonica… ▽ More The effect of a constant applied external force, induced for instance by an electric or gravitational field, on the dispersion of Brownian particles in periodic media with spatially varying diffusivity, and thus mobility, is studied. We show that external forces can greatly enhance dispersion in the direction of the applied force and also modify, to a lesser extent and in some cases non-monotonically, dispersion perpendicular to the applied force. Our results thus open up the intriguing possibility of modulating the dispersive properties of heterogeneous media by using externally applied force fields. These results are obtained via a Kubo formula which can be applied to any periodic advection diffusion system in any spatial dimension. △ Less

Submitted 16 July, 2015; originally announced July 2015.

Comments: Main text (4 pages, 2 figures) + Supplemental Material (4 pages)

Journal ref: Physical Review Letters, 115, 020601 (2015)

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

Authors: David S. Dean, Pierre Le Doussal, Satya N. Majumdar, Gregory Schehr

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

Submitted 7 January, 2016; v1 submitted 6 May, 2015; originally announced May 2015.

Comments: 5 pages + 8 pages of supplementary material, 2 figures, published version

Journal ref: Europhys. Lett. 112, 60001 (2015)

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

Authors: David S. Dean, Pierre Le Doussal, Satya N. Majumdar, Gregory Schehr

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

Submitted 9 December, 2015; v1 submitted 4 December, 2014; originally announced December 2014.

Comments: 5 pages + 6 pages of supplementary material, 3 figures, published version

Journal ref: Phys. Rev. Lett. 114, 110402 (2015)

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

Authors: David S. Dean, Thomas Guérin

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

Submitted 3 December, 2014; originally announced December 2014.

Comments: 8 pages, 2 figures, revtex

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

Authors: David S. Dean, V. Adrian Parsegian, R. Podgornik

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

Submitted 5 August, 2014; originally announced August 2014.

Comments: 25 pages Latex IOP Style, 2 figures

Diffusion in periodic, correlated random forcing landscapes

Authors: David S. Dean, Shamik Gupta, Gleb Oshanin, Alberto Rosso, Gregory Schehr

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

Submitted 29 August, 2014; v1 submitted 10 June, 2014; originally announced June 2014.

Comments: 5 pages (main text) + 2 pages (supplemental material); v2: 9 pages, 3 figures, published version

Journal ref: J. Phys. A: Math. Theor. 47, 372001 (2014)

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

Authors: David S. Dean, Gleb Oshanin

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

Submitted 15 May, 2014; originally announced May 2014.

Comments: 11 pages RevTex, 3 figure

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

Authors: David S. Dean, Rudolf Podgornik

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

Submitted 27 March, 2014; v1 submitted 17 December, 2013; originally announced December 2013.

Comments: 10 pages, RevTex, I figure. Final Published version, title change with respect to initial preprint, figure added, protocol for model slightly modified

Journal ref: Phys. Rev. E 89, 032117 (2014)

On the fluctuation of thermal van der Waals forces due to dipole fluctuations

Authors: D. S. Dean, V. A. Parsegian, R. Podgornik

Abstract: Fluctuations of the thermal or classical component of the van der Waals force between two dielectric slabs, modelled as an ensemble of polarizable dipoles which interact via the usual electrostatic dipole-dipole interaction, are evaluated. In the model the instantaneous force is a deterministic function of the dipole configurations in the slabs and its fluctuations are purely due to dipole fluctua… ▽ More Fluctuations of the thermal or classical component of the van der Waals force between two dielectric slabs, modelled as an ensemble of polarizable dipoles which interact via the usual electrostatic dipole-dipole interaction, are evaluated. In the model the instantaneous force is a deterministic function of the dipole configurations in the slabs and its fluctuations are purely due to dipole fluctuations (no background thermal fluctuations of the electromagnetic field are considered). The average of the force and its variance are computed. The fluctuations of the force exhibit normal thermodynamic scaling in that they are proportional to the area of the two plates, and even more importantly, do not depend on any microscopic cut-off in the theory. The average and the variance of the thermal van der Waals forces give a unique fingerprint of these fluctuation interactions. △ Less

Submitted 11 April, 2013; originally announced April 2013.

Comments: 5 pages Revtex

Journal ref: Phys. Rev. A 87, 032111 (2013)

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

Authors: Denis Boyer, David S. Dean, Carlos Mejía-Monasterio, Gleb Oshanin

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

Submitted 18 March, 2013; originally announced March 2013.

Comments: Mini-review paper presented at the 4-th Conference on Statistical Physics: Modern Trends and Applications, July 3-6, 2012 Lviv, Ukraine. It contains new results on the optimal L-S estimator of the diffusion constant for diffusion in the presence of a random potential

Journal ref: Eur. Phys. J. Special Topics 216, 57-71 (2013)