Skip to main content

ORIGINAL RESEARCH article

Front. Phys., 25 March 2021
Sec. Interdisciplinary Physics
This article is part of the Research Topic 50 years of Statistical Physics in Mexico: Development, State of the Art and Perspectives View all 17 articles

Entropic Effects of Interacting Particles Diffusing on Spherical Surfaces

  • 1Departamento de Matemáticas, Universidad Autónoma Metropolitana-Iztapalapa, Avenida San Rafael Atlixco No.186 Colonia Vicentina, Delegación Iztapalapa, México City, México
  • 2Departamento de Física, Universidad Autónoma Metropolitana-Iztapalapa, Avenida San Rafael Atlixco No.186 Colonia Vicentina, Delegación Iztapalapa, México City, México
  • 3Unidad Multidiciplinaria de Docencia e Investigación-Juriquilla, Facultad de Ciencias, Universidad Nacional Autónoma de México (UNAM), Boulevard Juriquilla, Querétaro, México

We present a molecular dynamics and theoretical study on the diffusion of interacting particles embedded on the surface of a sphere. By proposing five different interaction potentials among particles, we perform molecular dynamics simulations and calculate the mean square displacement (MSD) of tracer particles under a crowded regime of high surface density. Results for all the potentials show four different behaviors passing from ballistic and transitory at very short times, to sub-diffusive and saturation behaviors at intermediary and long times. Making use of irreversible thermodynamics theory, we also model the last two stages showing that the crowding induces a sub-diffusion process similar to that caused by particles trapped in cages, and that the saturation of the MSD is due to the existence of an entropic potential that limits the number of accessible states to the particles. By discussing the convenience of projecting the motions of the particles over a plane of observation, consistent with experimental capabilities, we compare the predictions of our theoretical model with the simulations showing that these stages are remarkably well described in qualitative and quantitative terms.

1 Introduction

In several physical and biological systems the mass transport phenomena is carried out in surfaces with non vanishing curvatures. The diffusion of bio-molecules and other particles on the surface of liposomes, drops and other curved entities is of high relevance due to its potential applications in biomedicine and technology [14]. Characteristic examples of this phenomena can be found in the diffusion of proteins on curved membranes or in the surface diffusion of molecules and chemical solvents in catalytic surfaces [58].

In these curved domains, many geometric aspects relating the motion of the particles on surfaces have been studied in recent years [914]. These studies are relevant in the understanding of applied problems such as nucleation, spinodal decomposition, adsorption, and phase transitions, since the dynamics of the particles is modified due to the curvature of the surface on which the molecules are embedded [1518].

In particular, the most relevant quantity to measure in these particle motions is the mean square displacement (MSD). However, most of the previously refereed studies provide geometric information of the MSD and the effective diffusion coefficient only for the case where the interaction of the particles occurs with the surface and not among particles themselves. In this simplification, many expressions can be found in the literature for approximating the short and long time behavior of the MSD, see Ref. [9] and references therein. However, for the case of interacting particles, there are only few numerical works studying the influence of interaction and confinement in the MSD [1921].

In this work we will consider how the surface diffusion displacements are influenced by the interaction of the particles and the projection of their displacement in a plane of observation. This is because for most experimental cases of interest, the geodesic trajectories of the particles cannot be followed at all times, and only the projection in a plane of observation can be measured [11, 22, 23]. This fact emphasizes that for diffusion over closed surfaces, the observed motility depends on the curvature of the surface and on the plane of observation [2426].

We present the results of the MSD for interacting particles in a crowded spherical surface by using Molecular Dynamics (MD) simulations and implementing five different pairwise-interaction potentials between the particles in a relatively high surface density medium (ρ=0.92). The particles are free to move on the entire spherical surface of radius R. To keep the particles on the surface, we use the algorithm described by Juffer et al. [27]. The continuous interaction potentials used are different cases of the generalized form of the Mie potentials [28]: LJTS, WCA(12.6), WCA(50.49). Besides, some simple discontinuous potentials like the shoulder-square (SS) and the well-square (SW) were adapted to their continuous-mathematical form described by Munguía-Valadez et al.1

With the aim to provide a physical interpretation of the results of these numerical calculations, we formulate an analytic model based on the generalized Smoluchowski equation with time-dependent coefficients [2934]. Using this tool, we can identify how the change of perspective, confinement and crowding are incorporated in an effective diffusion coefficient that takes into account the entropic confinement and the anomalous diffusion. The model used here was derived in the context of mesoscopic non-equilibrium thermodynamics (MNET) and has been widely used in order to obtain kinetic equations for transport phenomena, like diffusion-adsorption processes, anomalous diffusion, activated processes, diffusion in pores, and diffusion in the presence of entropic barriers [3539]. In particular for our purposes, MNET has been successful in describing diffusion on other confined systems [38, 4042].

The comparison between our numerical simulations and thermodynamic-based model allows us to propose a new interpretation of the observed dynamics of the MSD for diffusion of interacting particles over a sphere, and to compare the different subdiffusive regimes associated with the hardness or softness of the particle collisions. We have chosen a spherical geometry in order to rule out the effects of different local curvatures in the surface. This will also allow us to recover previous theoretical results deduced for the sphere in absence of interaction.

Regarding the organization of work, in Section 2 we present our Molecular Dynamics simulations on the diffusion of finite sized-particles over an spherical surface for different interaction potentials among particles. Then, we study the dynamics of the diffusion and the behavior of the MSD using the results provided by a Smoluchowski description which is presented in Sections 3.1 and 3.2, for free and interacting particles, respectively. Comparison between simulations and model is presented in Section 4. Finally, discussion and concluding remarks are provided in Section 5.

2 Molecular Dynamics Simulations

We present the details of the simulations as well as a brief description of the interaction potentials used in this work.

2.1 Interaction Potentials

The continuous interaction potentials are defined through:

ϕn,m(rij)=(nnm)(nm)m/(nm)ε[(σrij)n(σrij)m],(1)

where ϕn,m(rij) is the most general form of the Mie potential [28]. For example, we can write the Truncated and Shifted Lennard-Jones (LJTS) interaction by choosing n,m=12,6, and then

ϕLJTS(rij)={ϕ12,6(rij)ϕ12,6(rc),rijrc,0,rij>rc,(2)

where σ(=2a) is the particle diameter, rij=|rirj| the distance between the centers of mass of the i-th and the j-th particles, ϵ the potential well depth, and rc(=2.5σ) the cutoff radius. On the other hand, a generalization of the Weeks-Chandler-Anderson potential (WCA) [43, 44] can be written as

ϕWCA(n,m)(rij)={ϕn,m(rij)+ϵ,rij(nm)1/nmσ,0,rij>(nm)1/nmσ,(3)

where rij=(n/m)1/nmσ is the value that corresponds to the minimum potential. In our case, we use the parameters n,m=12,6, [WCA(12,6)] and n,m=50,49, [WCA(50,49)] to incorporate pseudo hard-sphere interaction models into the study, see Figure 1B [44].

FIGURE 1
www.frontiersin.org

FIGURE 1. Potentials of interacting particles used for the simulations. (A) Potential of Truncated and Shifted Lennard-Jones. (B)WCA(12,6) potential and the hard sphere pseudo-potential WCA(50,49). (C) and (D) Square-Shoulder potential continuous and continuous Square-Well respectively, given by Eq. 4.

The Square-Well (SW) and Square-Shoulder (SS) potentials can be approximated with the Generalized Continuous Multiple Step Potential (GCMS), which incorporates a contribution of excluded volume followed by multiple steps that model repulsive barriers or attractive wells as the case may be. Its simplified form can be written as:

ϕGCMS(σ,ω,q,a0,rij)=ϵ[(ωrijσ+ω)q+a01+exp[q(rijσω)/ω]],(4)

where ω defines the spatial extent of the particle core in units of σ, q is related to the hardness of the potential, and a0 is a factor that defines an attractive step (a0<0) or a repulsive one (a0>0). The parameters used to reproduce the SS and SW are: ω=0.5, q=500, a0=1 and a0=1, respectively, see Figure 1C y Figure 1D.

2.2 Simulation Method

In this work all quantities are assumed to be expressed in conventional reduced units, with m, σ and ϵ as the units of mass (set equal to 1), distance, and energy, respectively. According to this convention, the temperature (T) is in units of kBT/ϵ where kB is Boltzmanns constant, the density (ρ) in units of ρσ2, the time (δt) in units of δt(ϵ/mσ2)1/2, and the energy (ϕ) in units of ϕ/ϵ. A monodisperse system of spherical particles was studied with the restriction of moving on the surface of a sphere of radius R=9.3σ. To study the dynamic behavior of each system, Molecular Dynamics simulations were performed in the canonical assembly (NVT) for a total of N=1000 particles embedded in the spherical surface, with a surface density given by ρ(=N/A)=0.92. The R and ρ values were taken for compatibility with previous studies by Vest et al. [4547]. The system was placed in a thermal bath at a temperature T=1.0 set constant by using the isokinetic method; the integration of the equations of motion was carried out using the velocity-Verlet algorithm [48], with the restriction fi(ri)=|ri|2R2=0 (see Ref. [27]) and a time-step δt=103 for LJTS, WCA(50,49) and WCA(12,6) systems, and with δt=104 for SS and SW systems.

The simulations were run by 106 time-steps for LJTS, WCA(50,49) and WCA(12,6), and by 107 time-steps for SS and SW to reach thermodynamic equilibrium, then for 3×107 and 4×108 time-steps, respectively, for the calculation of the mean square displacement. Some representative trajectories of a tracer particle for each potential are plotted in Figure 2.

FIGURE 2
www.frontiersin.org

FIGURE 2. (A) Representative instantaneous configuration of the system with N = 1,000 particles on the spherical surface of density ρ(=N/A)=0.92 and with a radius of R=9.3σ. (B–F) Path of a tracer particle over the spherical surface for the different interaction potentials.

The results of our simulations for the MSD are plotted in Figure 3. At the right side we plot the geodesic MSD over the sphere with the saturation valor of (π24)R2/2. At the left, we plot the MSD measured when the displacements of the particles are projected the plane of observation XY, giving the saturation value of 2R2/3. These saturation values for long times only depend on the geometric configuration and therefore are the same for free and interacting diffusing particles [10].

FIGURE 3
www.frontiersin.org

FIGURE 3. MSD of interacting particles over the sphere with the potentials given in Figure 1. The geodesic arc, s2, and the projected displacement in the XY plane, r2, are plotted at left and right, respectively.

The MSD plotted in Figure 3 shows four different stages: 1) the ballistic regime where the MSD increases as t2, 2) a transition between ballistic and diffusive regime where the slope of the MSD decreases and changes of concavity, 3) the subdiffusive regime where, as we will calculate below, the MSD increases as tα with α1 and, finally, for long times, 4) the saturation regime where the MSD no longer increases since the domain of the sphere is finite. These four regimes are seen in the geodesic and projected version of the MSD plotted at right and left of Figure 3, respectively. Some or all of these four stages have been found both in numerical simulations [9, 4952] and experiments, mainly for protein and lipid diffusion on cell membranes [5356].

The form of the plots for the MSD in Figure 3 depends on the combined effect of the curvature of the surface, the interaction of the particles and, in the case of the projected displacement, the projection on the observation plane. In the next section we will deal with the connection between these different effects and their effect in the measured effective diffusion coefficient.

3 The MSD of Free and Interacting Particles

Commonly, experimental setups for measuring the diffusion of particles on curved surfaces use confocal microscopy techniques [33, 57, 58]. Due to the focal length inherent to these techniques, the measure of the motions of the particles under study is preferentially performed in terms of their projection to the focal plane [22, 23]. In view of this, in the present section we study the equivalence between the description of the free diffusion over the surface (the geodesic displacement) and compare with the evolution of its projection in the planar disk constrained by the sphere, see Figure 4. First, we will deduce an expression for the MSD of free particles moving at a flat disk by considering how the entropic restriction of the movement (imposed by the fact that the particles cannot travel beyond the radius of the disk) is reflected in the transport properties that are used to describe the motion. Then, we will modify the MNET description to include the interaction and crowding of particles.

FIGURE 4
www.frontiersin.org

FIGURE 4. Projection on an equatorial plane of the Brownian motion which takes place over the surface of a sphere of radius R. The polar radial coordinate of the particle r is measured on the plane.

3.1 The MSD of Free Punctual Particles

From a thermodynamic point of view, the problem can be tackled by estimating the equivalent force keeping the system confined to diffuse in the planar disk delimited by the surface of the sphere of radius R, see Figure 4. This equivalent force can be written in terms of a Taylor expansion. To first order, this correspond to an harmonic-like potential of the form

U(r)12κ0r2.(5)

The potential here is of an entropic origin, since it is due to the fact that the particles cannot leave the sphere projection or, equivalently, the planar disk. The maximum value of this potential, associated with the equivalent force, can be estimated by using the equipartition theorem [59]. Comparing the mean kinetic energy of the particles, K=(3/2)kT (which corresponds to one half for each degree of freedom in the three dimensional description), with the maximum potential energy available for the trapped particles, Umax=(1/2)κ0R2, we can estimate the value of the effective restorative coefficient κ0 as

κ0(R)=3kBTR2.(6)

Hence, the projected motion of the particles constrained to diffuse in a disk under the influence of the equivalent harmonic force can be described by means of the Smoluchowsky equation

Pt=D02PD0kBT(PF),(7)

where kB and T are Boltzmann constant and temperature, respectively. Since the equivalent force has radial symmetry, assuming an initial configuration with the same symmetry allow us to reduce the calculation of the MSD to the contribution in the radial coordinate r. In this case, we have

pt=D0rr[rpr+κ0r2kBTp](8)

where the force that constraint the particles motions is F=U, with U given by Eq. 5, and the gradient and Laplacian operators are represented in polar coordinates.

The mean square displacement in the constrained disk system is therefore given by

r2=0R02πr2prdr.(9)

The time evolution of the MSD is obtained by taking the derivative

dr2dt=r2pdtdA,(10)

And substituting Eq. 8. Assuming zero flux boundary condition: p/r(R)=0 in the border of the disk and the conservation normalization equation, 0Rprdr=1, we can prove, using successive integration by parts, that

dr2dt=4D02D0kBTκ0r2.(11)

The comparison of the previous result with the result of diffusion in an infinite planar surface, r2=4D0t, shows that the effect of the confinement is proportional to the ratio of the potential and kinetic energies:

dr2dt=4D0(132UK),(12)

where the average here is over the ensemble of non-interacting particles. Eq. 11 is readily solved using that r20 when t0 as

r2=2kBTκ0[1e2κ0D0t/kBT].(13)

Substituting the estimated value of κ0 given in Eq. 6, valid when particle interactions are negligible, we obtain

r2=2R23[1e6D0t/R2].(14)

This result has been obtained by using different methods in Ref. [10] and, as expected, only depends upon geometrical considerations resulting from averaging methodologies applied directly to the Laplace-Beltrami diffusion equation.

This result describing the projection of the surface to the disk can be contrasted to the result of particles really moving on a planar disk of radius R, where K=kBT (since there are only two freedom degrees), κ0=2kBT/R2, and therefore, from Eq. 11, we have

r22D=R2[1e4D0t/R2].(15)

Comparison between Eqs 14, 15 shows that the displacement in a sphere respect to a planar projection has an average reduction by a factor of 2/3 for long times. This reduction factor in the displacements can be deduced also by evaluating the ratio between the effective diffusion coefficient measured on a planar region, D0, with respect to that measured in the observation plane, D. The result is

DD0=12(1+nz2),(16)

Where nz2 is the geometric average of the quadratic vertical component of the normal vector [10, 59]. For the sphere we have nz=cosθ, from where it is straightforward deduced that nz2=1/3 giving the expected result D/D0=2/3.

3.2 The MSD of Interacting Non-punctual Particles

In the previous subsection we have shown that the entropic restrictions, present by the fact that tracer particles move over a surface with a finite number of accessible states, are responsible for the observed saturation value reached by the MSD at long times. However, as we have seen in Figure 3, the simulations show that the time behavior of the MSD at intermediate times scales with a time dependence tα with α1. Therefore, for certain interaction potentials, anomalous diffusion is observed. This fact requires that the model developed in Section 3.1 should be generalized to cope with particle interaction.

This generalization goes along the same lines of Santamaría-Holek et al. [32], where a Smoluchowski description with time-dependent coefficients was used to demonstrate that the exponent characterizing the sub-diffusion is controlled by the nature of the local cages and the free space at disposal for their motion. In our case, since the single particle MSD reported here accounts for the motion of independent particles in a crowded media, the description of the diffusion can be given also in terms of the generalized Smoluchowski equation for the single-particle distribution function [61]. This equation for the probability is the counterpart of the respective generalized Langevin equation with memory effects, commonly used to describe microrheological experiments and has the peculiarity of being equivalent in the long-time regime to the Smoluchowski equation with time dependent coefficients [29, 30]. This fact was generalized in Refs. [31, 33, 62] for the case of power-law dependent memory kernels which are equivalent to time dependent diffusion coefficients of the form

Ds(t)D0(tτ)α1,(17)

where τ is a characteristic time of the anomalous diffusion process. In this equation, the exponent α(0,1] characterizes how strong sub-diffusion is and, as it was shown in Ref. [32], has the general form

α=(1B1aξ)3kBTκ0R2.(18)

Besides the entropic effect introduced through the ratio kBT/κ0R2, the Eq. 18, through the factor 1B1(a/ξ), contains the two main ingredients mentioned previously, namely, the free-area at disposal to perform diffusion a/ξ, and the nature of the local cages through the coefficient B1. The characteristic length ξ is related with the free-area in which the finite particles can move in a surface of total area A: ξ=Afree1/2R(1Ma2/R2)1/2 where M is the number of particles forming the cage. In this case, we have

α=(1B1aRR2Ma2),(19)

This expression for the sub-diffusion exponent shows that it only depends in the radii of the particles and sphere, and on the parameter B1. This parameter indirectly depends in the nature of the interaction potentials used. In the context of hydrodynamics [63], it was related with the correction introduced by hydrodynamic interactions over the motion of a particle when this motion takes places near to a solid wall. Since these interactions are related to the potential interactions among particles, one may conjecture that the parameter B1 is a measure of the effect that different interaction potentials have on the cages’ dynamic structure and, therefore, of the different values that the sub-diffusion exponent may take.

Even in this approximation, it results difficult to estimate theoretically the parameters in Eq. 18. However, it is clear that the exponent of the subdiffusive process is linked to the shape of the interacting potential, the curvature of the surface and the projection of the forces to the observation plane. In order to find its value for the different potentials in Figure 1, we will adjust the data of the numerical simulations, using the Smoluchowski equation, but now in terms of this time dependent diffusion coefficient as we detail below.

Therefore, we will consider that the single particle Smoluchowski equation describing the dynamics of a tracer particle in an effective medium of N1 interacting non-punctual particles takes the form

pt=Ds(t)rr[rpr+1kBTκ0r2p],(20)

Where now the time dependent coefficient is given through Eqs. 17, 19. From the last equation, repeating the same procedure of the last section, the temporal evolution of the MSD is obtained to be

r2=2kBTκ0{1exp[2κ0kBT(0tDs(t)dt.)]}.(21)

Using Eq. 6 we finally have

r2=23R2{1exp[6R2(0tDs(t)dt)]},(22)

Which is the general expression for the MSD of the tracer particles valid for the whole time interval. From this equation, it is possible to deduce the temporal behavior of the diffusion coefficient for all the stages of the process in terms of the numerical data provided by the simulations as

Ds(t)=R26(dr2/dt)(t)((2/3)R2r2(t)).(23)

Notice that this equation reproduces the expected behavior for short and long times. For short times r20 and dr2/dt4D0 and therefore DsD0. For long times, it is expected that the MSD saturates and, therefore, dr2/dt0 with Ds0.

At the other hand, when (Eq. 22) is used in combination with the diffusion coefficient in Eq. 17, it accounts only for the process once the sub-diffusion stage has started. At intermediary times, it accounts for the sub-diffusion associated with particle interactions and cage diffusion effects and, for long times, it account for the fact that the particles have a restricted number of accessible states in a finite domain. In the next section we will use the numerical simulation in order to find the general behavior of the temporal diffusion coefficient and the exponents of sub-diffusion for the potentials we have considered in Section 2.

4 Comparison Between Model and Simulations

We have shown that the projection on the equatorial plane of the geodesic MSD, directly measured in the numerical simulations, is an equivalent form to characterize the diffusive dynamics of the particles. We have explained how the saturation value of the MSD rescales due to the restricted number of accessible states, giving rise to entropic constrictions on the dynamics of the particles.

For short and intermediate times, our methodology allows us to consider how different aspects modify the behavior of the MSD of interacting particles forced to diffuse in the surface of the sphere. All these aspects are clearly represented in Figure 5, where we plot the time-dependent diffusion coefficient, Ds, given in Eq. 23 for the potentials chosen, and illustrated in Figure 1. For constructing this plot, the data from the functions r2(t) and dr2/dt obtained from numerical simulations are directly substituted in (Eq. 17) for Ds(t).

FIGURE 5
www.frontiersin.org

FIGURE 5. Projected surface diffusion coefficient Ds(t) as a function of the time for the interacting particles with potentials given in Figure 1. The dashed lines represent the reference value D0 in a planar surface.

The first aspect to notice is the reference diffusion coefficient, D0. This is marked as a dashed line in Figure 5 and represents the surface diffusion coefficient in an infinite planar surface, and negligible interaction potential among particles. In Figure 5 it is seen that, for the used potentials, the value of the diffusion coefficients augments according to the order: SW, LJTS, SS, and WCA(12,6). The color key is the same as in Figure 1.

The second aspect to notice in the temporal dependence of the diffusion coefficient Ds(t) is the emergence of sub-diffusion due to the crowding effects. This results in different slopes of the curve Ds(t) at intermediary times and, therefore, corresponds to different values of the sub-diffusion exponent α. In order to compare quantitatively the sub-diffusion for the different potentials, in Figure 6, we adjust the data from numerical simulations using (Eq. 22) with the temporal dependence of the diffusion coefficient in (Eq. 17) using the least squares method. The value of the local diffusion coefficient D0 and the sub-diffusion exponent α together with the fits of the data are given in Figure 6.

FIGURE 6
www.frontiersin.org

FIGURE 6. Comparison between the projected MSD obtained from simulations (colored solid lines) and theory (black lines) according to Eq. 22. The values of the parameters used are included in each plot. The scaling factor in Eq. 17 is τ=1 for all cases.

Figure 5, show that the value of the diffusion constant D0 is smaller for hardest potentials (WCA(12,6) and SS) than for softer potentials (WCA(50,49), LJTS and SW). The potentials having wells also show this tendency since the diffusion constant for the LJTS is larger that of the SW potential. The attractive part of the potentials decreases, in turn, the kinetic energy of the particles surrounding the tracer particle and, therefore, the energy availability to perform position fluctuations. This interpretation is clear after a comparison between the potentials WCA(50,49) and WCA(12,6), which have a similar functional form but different degree of penetrability. Another aspect to be noticed is the steepest drop of the diffusion coefficients (lower α) for the potentials with a tilted repulsive part after r=1, that is, SS and WCA(12,60). A better understanding of the effect of the shape and hardness of the interaction potentials on the diffusion constants and the sub-diffusion process will require its own work.

Finally, the third aspect to notice is that both the ballistic and long-time behaviors of the MSD appear in all cases since not depend on the interaction potential. Because of this, they do not provide information of the fluid medium where the diffusion occurs. As it is seen in Figure 5, the final drop of the diffusion coefficient to zero occurs earlier as the diffusion coefficient is lower, as expected.

Our work shows that the deviations of the MSD relative to the planar behavior represent the effect that the particle interaction and crowding effects have on the dynamics. This effect is coupled with the presence of an entropic force restricting the number of accessible states. We prove that these two effects are well captured by the projected diffusion coefficient Ds, given in Figure 5 and which has been reported in previous numerical simulations and experiments [53, 64].

The excellent agreement in Figure 6 among simulation results and the theory proposed indicates that particle interactions control the cage effect and the sub-diffusion regime. The delay induced by the crowded scenario appears since the transient stage, and its influence is captured by the memory effects that the effective diffusion coefficient Ds(t) incorporates. We have found that the value of the parameters used for fitting the MSD curves depend upon the used potentials and further numerical and theoretical studies are necessary in order to show the relation between the form of the potential interaction and the surface diffusion coefficient. Similar commentaries can be done concerning the specific value for the radii of the particles and the surface density. These values allowed us to capture the effect of the surface density and the form of the potential on the sub-diffusion observed. Nonetheless, a more detailed study on the dependence of the diffusion properties on this factor will require a future work.

5 Conclusion

We have studied the statistical and dynamical properties of a system of interacting particles constituting a fluid embedded on a sphere’s surface. To understand the link between the microscopic properties of the fluid and the observed MSD, we have proposed a diffusion model incorporating memory and entropic effects that allowed us to derive a simple and powerful expression for the projected time-dependent diffusion coefficient. This expression allows us to characterize the dependence of the diffusive dynamics on the particle interactions and surface curvature. This approach takes into account the change of perspective inherent to the experimental measure of diffusion of particles on closed surfaces. The formalism proposed in this work reproduces also some limiting cases previously studied in the literature and exemplifies, with the use of ideas emerging from irreversible thermodynamics, how we can provide a more flexible description that helps to interpret complex phenomena occurring at the surface.

In this work, we have shown in Figure 3 that the dynamics of the MSD reflects four stages. Two, at short times t101s, are associated with the ballistic regime and how it ceases; and other two (t101s) associated with the diffusion regime in presence of cages at the intermediate and long time regimes, corresponding to anomalous diffusion and saturation, respectively. Using the ideas of the MNET we use the Smoluchowski equation for describing the process once the diffusion regime starts and this is shown Figure 6 for the intermediary and long time behaviors of the MSD. In contrast, modeling the regimes at short times will require another kind of kinetic mechanism. However, it is worth mentioning that, under the perspective of biological applications of surface diffusion in crowded membrane environments we discuss below, only the subdiffusive regime is accessible to measurements with the actual experimental time-resolution capabilities and therefore, this does not constitute a limitation to our conclusions.

In the realm of applications, this study allows us to establish a connection with models and experiments studying the diffusion of several tracers in the cytoplasm of living cells which often exhibit heterogeneous distribution of macromolecular crowders [6567]. This crowding is extremely important in the surface membrane and might affect the properties of anomalous diffusion. It is known, for example, that going from a less to a more crowded region will slow down the dynamics, both in terms of the exponent and diffusivity [6870]. This results very important since the slow transport in the cell membrane of lipids and proteins is linked with protein cluster formation, phase segregation, lipid droplet formation, signal propagation and other crucial functions occurring on the cell surface which can be enhanced by the presence of sub-diffusion [71, 72].

Data Availability Statement

The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author.

Author Contributions

IS-H, AL-D and JM-R designed the research. AL-D and JM-V performed the calculations and methology. IS-H, JM-R, AL-D and JM-V performed analysis. IS-H, JM-R, AL-D and SH wrote the original draft. All authors contributed to the analysis and discussion of the data and the writing of the paper.

Funding

SIH is grateful to projects UNAMDGAPA-PAPIIT IN114721 and LANCAD-UNAM-DGTIC-276. ISH acknowledges UNAM-DGAPA for financial support under grant number IN117419.

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

Authors appreciate the technical support of Carlos Sair Flores Bautista, Alejandro de León Cuevas, Jair Santiago García Sotelo and Luis Alberto Aguilar Bautista from Laboratorio Nacional de Visualización Científica Avanzada (LAVIS-UNAM)

Footnotes

1Munguía-Valadez J, Chávez-Rojo MA, Moreno-Razo JA, Sambriski EJ (under review). The generalized continuous multiple step potential: model systems and related properties. J Chem Phys.

References

1. Moreno-Razo JA, Sambriski EJ, Abbott NL, Hernández-Ortiz JP, de Pablo JJ. Liquid-crystal-mediated self-assembly at nanodroplet interfaces. Nature (2012) 485(7396):86. doi:10.1038/nature11084

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Holly C. Gaede and Klaus Gawrisch. Lateral diffusion rates of lipid, water, and a hydrophobic drug in a multilamellar liposome. Biophys J (2003) 85(3):1734–40. doi:10.1016/S0006-3495(03)74603-7

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Nuytten N, Hakimhashemi M, Ysenbaert T, Defour L, Trekker J, Soenen SJ, et al. Pegylated lipids impede the lateral diffusion of adsorbed proteins at the surface of (magneto)liposomes. Colloids Surf B Biointerfaces (2010) 80(2):227–31. doi:10.1016/j.colsurfb.2010.06.009

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Danelian E, Karlén A, Karlsson R, Winiwarter S, Hansson A, Löfâs S, et al. Spr biosensor studies of the direct interaction between 27 drugs and a liposome surface: correlation with fraction absorbed in humans. J Med Chem (2000) 43(11):2083–6. doi:10.1021/jm991156g

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Ledesma-Durán A, Hernández SI, Santamaría-Holek I. Effect of surface diffusion on adsorption-desorption and catalytic kinetics in irregular pores. I. Local kinetics. J Phys Chem C (2017) 121(27):14544. doi:10.1021/acs.jpcc.7b03652

CrossRef Full Text | Google Scholar

6. Ledesma-Durán A, Hernández SI, Santamaría-Holek I. Effect of surface diffusion on adsorption-desorption and catalytic kinetics in irregular pores. II. Macro-kinetics. J Phys Chem C (2017) 121(27):14557. doi:10.1021/acs.jpcc.7b03653

CrossRef Full Text | Google Scholar

7. Faraudo J. Diffusion equation on curved surfaces. I. Theory and application to biological membranes. J Chem Phys (2002) 116(13):5831. doi:10.1063/1.1456024

CrossRef Full Text | Google Scholar

8. Reister-Gottfried E, Leitenberger S, Seifert U. Diffusing proteins on a fluctuating membrane: analytical theory and simulations. Phys Rev E - Stat Nonlinear Soft Matter Phys (2010) 81(3):031903. doi:10.1103/PhysRevE.81.031903

CrossRef Full Text | Google Scholar

9. Castro-Villarreal P, Villada-Balbuena A, Méndez-Alcaraz JM, Castañeda-Priego R, Estrada-Jiménez S. A brownian dynamics algorithm for colloids in curved manifolds. J Chem Phys (2014) 140(21):214115. doi:10.1063/1.4881060

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Castro-Villarreal P. Intrinsic and extrinsic measurement for brownian motion. J Stat Mech (2014) 2014(5):P05017. doi:10.1088/1742-5468/2014/05/p05017

CrossRef Full Text | Google Scholar

11. Reister-Gottfried E, Leitenberger SM, Seifert U. Hybrid simulations of lateral diffusion in fluctuating membranes. Phys Rev E - Stat Nonlinear Soft Matter Phys (2007) 75(1):011908. doi:10.1103/PhysRevE.75.011908

CrossRef Full Text | Google Scholar

12. Castro-Villarreal P. Brownian motion meets riemann curvature. J Stat Mech (2010) 2010(8):P08006. doi:10.1088/1742-5468/2010/08/p08006

CrossRef Full Text | Google Scholar

13. Apaza L, Sandoval M. Brownian self-driven particles on the surface of a sphere. Phys Rev E (2017) 96(2):022606. doi:10.1103/PhysRevE.96.022606

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Petersen CF, Schrack L, Franosch T. Static properties of quasi-confined hard-sphere fluids. J Stat Mech (2019) 2019(8):083216. doi:10.1088/1742-5468/ab3342

CrossRef Full Text | Google Scholar

15. Jinnai H, Koga T, Nishikawa Y, Hashimoto T, Hyde ST. Curvature determination of spinodal interface in a condensed matter system. Phys Rev Lett (1997) 78(11):2248. doi:10.1103/physrevlett.78.2248

CrossRef Full Text | Google Scholar

16. Saito A, Foley HC. Curvature and parametric sensitivity in models for adsorption in micropores. AIChE J (1991) 37(3):429. doi:10.1002/aic.690370312

CrossRef Full Text | Google Scholar

17. Tvergaard V. Effect of yield surface curvature and void nucleation on plastic flow localization. J Mech Phys Solid (1987) 35(1):43. doi:10.1016/0022-5096(87)90027-5

CrossRef Full Text | Google Scholar

18. Roux A, Cuvelier D, Nassoy P, Prost J, Bassereau P, Goud B. Role of curvature and phase transition in lipid sorting and fission of membrane tubules. EMBO J (2005) 24(8):1537. doi:10.1038/sj.emboj.7600631

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Resende FJ, Costa BV. Molecular-dynamics study of the diffusion coefficient on a crystal surface. Phys Rev B (2000) 61(19):12697. doi:10.1103/physrevb.61.12697

CrossRef Full Text | Google Scholar

20. Wang D, He C, Stoykovich MP, Schwartz DK. Nanoscale topography influences polymer surface diffusion. ACS Nano (2015) 9(2):1656. doi:10.1021/nn506376n

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Bulnes FM, Pereyra VD, Riccardo JL. Collective surface diffusion: n-fold way kinetic Monte Carlo simulation. Phys Rev E (1998) 58(1):86. doi:10.1103/physreve.58.86

CrossRef Full Text | Google Scholar

22. Benavides-Parra JC, Jacinto-Méndez D, Brotons G, Carbajal-Tinoco MD. Brownian motion near a liquid-gas interface. J Chem Phys (2016) 145(11):114902. doi:10.1063/1.4962746

CrossRef Full Text | Google Scholar

23. Sandoval-Jiménez IM, Jacinto-Méndez D, Toscano-Flores LG, Carbajal-Tinoco MD. Brownian-particle motion used to characterize mechanical properties of lipid vesicles. J Chem Phys (2020) 152(1):014901. doi:10.1063/1.5133092

PubMed Abstract | CrossRef Full Text | Google Scholar

24. King MR. Apparent 2-d diffusivity in a ruffled cell membrane. J Theor Biol (2004) 227(3):323–6. doi:10.1016/j.jtbi.2003.11.010

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Naji A, Brown FL. Diffusion on ruffled membrane surfaces. J Chem Phys (2007) 126(23):235103. doi:10.1063/1.2739526

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Gov NS. Diffusion in curved fluid membranes. Phys Rev E-Stat Nonlinear Soft Matter Phys (2006) 73(4):041918. doi:10.1103/PhysRevE.73.041918

CrossRef Full Text | Google Scholar

27. Juffer AH, Berendsen HJC. Dynamic surface boundary conditions. Mol Phys (1993) 79(3):623. doi:10.1080/00268979300101501

CrossRef Full Text | Google Scholar

28. Mie G. Zur kinetischen Theorie der einatomigen Körper. Ann Phys (1903) 316(8):657. doi:10.1002/andp.19033160802

CrossRef Full Text | Google Scholar

29. Adelman SA. Fokker-Planck equations for simple non‐Markovian systems. J Chem Phys (1976) 64(1):124–30. doi:10.1063/1.431961

CrossRef Full Text | Google Scholar

30. Tokuyama M, Mori H. Statistical-mechanical theory of random frequency modulations and generalized brownian motions. Prog Theor Phys (1976) 55(2):411–29. doi:10.1143/ptp.55.411

CrossRef Full Text | Google Scholar

31. Metzler R, Barkai E, Klafter J. Deriving fractional fokker-planck equations from a generalised master equation. Europhys Lett (1999) 46(4):431. doi:10.1209/epl/i1999-00279-7

CrossRef Full Text | Google Scholar

32. Santamaría-Holek I, Rubi JM. Finite-size effects in microrheology. J Chem Phys (2006) 125(6):64907. doi:10.1063/1.2241190

PubMed Abstract | CrossRef Full Text | Google Scholar

33. Santamaría-Holek I, Rubí JM, Gadomski A. Thermokinetic approach of single particles and clusters involving anomalous diffusion under viscoelastic response. J Phys Chem B (2007) 111(9):2293. doi:10.1021/jp0675375

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Sevilla FJ, Rodríguez RF, Gómez-Solano JR. Generalized ornstein-uhlenbeck model for active motion. Phys Rev E (2019) 100:032123. doi:10.1103/PhysRevE.100.032123

PubMed Abstract | CrossRef Full Text | Google Scholar

35. Reguera D, Rubí JM, Vilar JM. The mesoscopic dynamics of thermodynamic systems. J Phys Chem B (2005) 109(46):21502. doi:10.1021/jp052904i

PubMed Abstract | CrossRef Full Text | Google Scholar

36. Pan W, Vekilov PG, Lubchenko V. Origin of anomalous mesoscopic phases in protein solutions. J Phys Chem B (2010) 114(22):7620. doi:10.1021/jp100617w

PubMed Abstract | CrossRef Full Text | Google Scholar

37. Dekker H. Nonisothermal activation: nonequilibrium thermodynamics of metastable mesoscopic systems. Phys Rev, A (1991) 43(8):4224. doi:10.1103/physreva.43.4224

PubMed Abstract | CrossRef Full Text | Google Scholar

38. Santamaría-Holek I, Hernández S, García-Alcántara C, Ledesma-Durán A. Review on the macro-transport processes theory for irregular pores able to perform catalytic reactions. Catalysts (2019) 9(3):281. doi:10.3390/catal9030281

CrossRef Full Text | Google Scholar

39. Santamaría-Holek I, Grzywna ZJ, Rubi JM. Entropic effects in diffusion-adsorption processes in micropores. Eur Phys J Spec Top (2013) 222(1):129. doi:10.1140/epjst/e2013-01831-2

CrossRef Full Text | Google Scholar

40. Reguera D, Rubí JM. Kinetic equations for diffusion in the presence of entropic barriers. Phys Rev E-Stat Nonlinear Soft Matter Phys (2001) 64(6):061106. doi:10.1103/PhysRevE.64.061106

CrossRef Full Text | Google Scholar

41. Chacón-Acosta G, Pineda I, Dagdug L. Diffusion in narrow channels on curved manifolds. J Chem Phys (2013) 139(21):214115. doi:10.1063/1.4836617

PubMed Abstract | CrossRef Full Text | Google Scholar

42. Santamaría-Holek I, Ledesma-Durán A, Hernández SI, García-Alcántara C, Andrio A, Compañ V. Entropic restrictions control the electric conductance of superprotonic ionic solids. Phys Chem Chem Phys (2020) 22(2):437. doi:10.1039/c9cp05486c

PubMed Abstract | CrossRef Full Text | Google Scholar

43. Andersen HC, Weeks JD, Chandler D. Relationship between the hard-sphere fluid and fluids with realistic repulsive forces. Phys Rev A (1971) 4(4):1597. doi:10.1103/physreva.4.1597

CrossRef Full Text | Google Scholar

44. Jover J, Haslam AJ, Galindo A, Jackson G, Müller EA. Pseudo hard-sphere potential for use in continuous molecular-dynamics simulation of spherical and chain molecules. J Chem Phys (2012) 137(14):144505. doi:10.1063/1.4754275

PubMed Abstract | CrossRef Full Text | Google Scholar

45. Julien-Piera V, Gilles T, Pascal V. Dynamics of a monodisperse Lennard-Jones system on a sphere. Mol Phys (2014) 112(9–10):1330. doi:10.1080/00268976.2014.901568

CrossRef Full Text | Google Scholar

46. Julien-Piera V, Gilles T, Pascal V. Mode-coupling approach for the slow dynamics of a liquid on a spherical substrate. J Chem Phys (2015) 143(8):084505. doi:10.1063/1.4928513

PubMed Abstract | CrossRef Full Text | Google Scholar

47. Julien-Piera V, Gilles T, Pascal V. Glassy dynamics of dense particle assemblies on a spherical substrate. J Chem Phys (2018) 148(16):164501. doi:10.1063/1.5027389

PubMed Abstract | CrossRef Full Text | Google Scholar

48. Swope WC, Andersen HC, Berens PH, Wilson KR. A computer simulation method for the calculation of equilibrium constants for the formation of physical clusters of molecules: application to small water clusters. J Chem Phys (1982) 76(1):637. doi:10.1063/1.442716

CrossRef Full Text | Google Scholar

49. Yamamoto E, Akimoto T, Yasui M, Yasuoka K. Origin of subdiffusion of water molecules on cell membrane surfaces. Sci Rep (2014) 4:4720. doi:10.1038/srep04720

PubMed Abstract | CrossRef Full Text | Google Scholar

50. Chipot C, Comer J. Subdiffusion in membrane permeation of small molecules. Sci Rep (2016) 6:35913. doi:10.1038/srep35913

PubMed Abstract | CrossRef Full Text | Google Scholar

51. Tan P, Liang Y, Xu Q, Mamontov E, Li J, Xing X, et al. Gradual crossover from subdiffusion to normal diffusion: a many-body effect in protein surface water. Phys Rev Lett (2018) 120(24):248101. doi:10.1103/PhysRevLett.120.248101

PubMed Abstract | CrossRef Full Text | Google Scholar

52. Ghosh SK, Cherstvy AG, Grebenkov DS, Metzler R. Anomalous, non-Gaussian tracer diffusion in crowded two-dimensional environments. New J Phys (2016) 18(1):013027. doi:10.1088/1367-2630/18/1/013027

CrossRef Full Text | Google Scholar

53. Schwille P, Korlach J, Webb WW. Fluorescence correlation spectroscopy with single-molecule sensitivity on cell and model membranes. Cytometry (1999) 36(3):176–82. doi:10.1002/(sici)1097-0320(19990701)36:3<176::aid-cyto5>3.0.co;2-f

PubMed Abstract | CrossRef Full Text | Google Scholar

54. Golan Y, Sherman E. Resolving mixed mechanisms of protein subdiffusion at the t cell plasma membrane. Nat Commun (2017) 8(1):15851–15. doi:10.1038/ncomms15851

PubMed Abstract | CrossRef Full Text | Google Scholar

55. Ehrig J, Petrov EP, Schwille P. Near-critical fluctuations and cytoskeleton-assisted phase separation lead to subdiffusion in cell membranes. Biophys J (2011) 100(1):80–9. doi:10.1016/j.bpj.2010.11.002

PubMed Abstract | CrossRef Full Text | Google Scholar

56. Weiss M, Elsner M, Kartberg F, Nilsson T. Anomalous subdiffusion is a measure for cytoplasmic crowding in living cells. Biophys J (2004) 87(5):3518–24. doi:10.1529/biophysj.104.044263

PubMed Abstract | CrossRef Full Text | Google Scholar

57. Zhong Y, Zhao L, Tyrlik PM, Wang G. Investigating diffusing on highly curved water-oil interface using three-dimensional single particle tracking. J Phys Chem C (2017) 121(14):8023–32. doi:10.1021/acs.jpcc.7b01721

CrossRef Full Text | Google Scholar

58. Sarmiento-Gomez E, Santamaría-Holek I, Castillo R. Mean-square displacement of particles in slightly interconnected polymer networks. J Phys Chem B (2014) 118(4):1146. doi:10.1021/jp4105344

PubMed Abstract | CrossRef Full Text | Google Scholar

59. Mendoza CI, Santamaría-Holek I, Pérez-Madrid A. Effective temperatures and the breakdown of the Stokes-einstein relation for particle suspensions. J Chem Phys (2015) 143(10):104506. doi:10.1063/1.4930550

PubMed Abstract | CrossRef Full Text | Google Scholar

60. Gustafsson S, Halle B. Diffusion on a flexible surface. J Chem Phys (1997) 106(5):1880–7. doi:10.1063/1.473326

CrossRef Full Text | Google Scholar

61. Mason TG, Ganesan K, van Zanten JH, Wirtz D, Kuo SC. Particle tracking microrheology of complex fluids. Phys Rev Lett (1997) 79(17):3282. doi:10.1103/physrevlett.79.3282

CrossRef Full Text | Google Scholar

62. Rodríguez RF, Salinas-Rodríguez E, W Dufty J. Fokker-Planck and Langevin descriptions of fluctuations in uniform shear flow. J Stat Phys (1983) 32(2):279–98. doi:10.1007/BF01012711

CrossRef Full Text | Google Scholar

63. Happel J, Brenner H. Low Reynolds number hydrodynamics: with special applications to particulate media, Vol. 1. Berlin, Germany: Springer Science & Business Media (2012).

64. Saxton MJ. A biological interpretation of transient anomalous subdiffusion. i. qualitative model. Biophys J (2007) 92(4):1178–91. doi:10.1529/biophysj.106.092619

PubMed Abstract | CrossRef Full Text | Google Scholar

65. Szymanski J, Weiss M. Elucidating the origin of anomalous diffusion in crowded fluids. Phys Rev Lett (2009) 103(3):038102. doi:10.1103/PhysRevLett.103.038102

PubMed Abstract | CrossRef Full Text | Google Scholar

66. Sharonov A, Bandichhor R, Burgess K, Petrescu AD, Schroeder F, Kier AB, et al. Lipid diffusion from single molecules of a labeled protein undergoing dynamic association with giant unilamellar vesicles and supported bilayers. Langmuir (2008) 24(3):844–50. doi:10.1021/la702600w

PubMed Abstract | CrossRef Full Text | Google Scholar

67. Javanainen M, Hammaren H, Monticelli L, Jeon JH, Miettinen MS, Martinez-Seara H, et al. Anomalous and normal diffusion of proteins and lipids in crowded lipid membranes. Faraday Discuss (2013) 161:397–59. doi:10.1039/c2fd20085f

PubMed Abstract | CrossRef Full Text | Google Scholar

68. Guigas G, Weiss M. Sampling the cell with anomalous diffusion - the discovery of slowness. Biophys J (2008) 94(1):90–4. doi:10.1529/biophysj.107.117044

PubMed Abstract | CrossRef Full Text | Google Scholar

69. Jeon JH, Monne HM, Javanainen M, Metzler R. Anomalous diffusion of phospholipids and cholesterols in a lipid bilayer and its origins. Phys Rev Lett (2012) 109(18):188103. doi:10.1103/PhysRevLett.109.188103

PubMed Abstract | CrossRef Full Text | Google Scholar

70. Goose JE, Sansom MS. Reduced lateral mobility of lipids and proteins in crowded membranes. PLoS Comput Biol (2013) 9(4):e1003033. doi:10.1371/journal.pcbi.1003033

PubMed Abstract | CrossRef Full Text | Google Scholar

71. Metzler R, Jeon JH, Cherstvy AG. Non-brownian diffusion in lipid membranes: experiments and simulations. Biochim Biophys Acta (2016) 1858(10):2451–67. doi:10.1016/j.bbamem.2016.01.022

PubMed Abstract | CrossRef Full Text | Google Scholar

72. Guigas G, Weiss M. Effects of protein crowding on membrane systems. Biochim Biophys Acta (2016) 1858(10):2441–50. doi:10.1016/j.bbamem.2015.12.021

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: molecular simulations, mean square displacement, surface sphere, diffusion coefficient, curved surfaces

Citation: Ledesma-Durán A, Munguía-Valadez J, Moreno-Razo JA, Hernández SI and Santamaría-Holek I (2021) Entropic Effects of Interacting Particles Diffusing on Spherical Surfaces. Front. Phys. 9:634792. doi: 10.3389/fphy.2021.634792

Received: 28 November 2020; Accepted: 08 January 2021;
Published: 25 March 2021.

Edited by:

Ramon Castañeda-Priego, University of Guanajuato, Mexico

Reviewed by:

Leonardo Dagdug, Metropolitan Autonomous University, Mexico
Pavel Castro-Villarreal, Autonomous University of Chiapas, Mexico
Andrey Cherstvy, University of Potsdam, Germany

Copyright © 2021 Ledesma-Durán, Munguía-Valadez, Moreno-Razo, Hernández and Santamaría-Holek. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: I. Santamaria Holek, isholek.fc@gmail.com

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.