Static and dynamic properties of decane/water microemulsions stabilized by cetylpyridinium chloride cationic surfactant and octanol cosurfactant

Molecular dynamics simulation (MD) is used to study the static and dynamic properties of positively charged decane/water microemulsions, for various volume fractions Φ (2.8%, 6.98%, 14%, and 26.5%). An effective hybrid potential combining three potentials, namely the hard-sphere repulsive potential, the van der Waals attractive potential, and the Yukawa repulsive potential, is used to describe the microemulsion interactions. The microemulsion shape is determined using the renormalized spectra in Porod representation. The appropriate potential parameters are tested using the Ornstein–Zernike integral equation approach with the Hypernetted Chain (HNC) closure relation by a comparison between the structure factor calculated from HNC and that obtained from Small Angle Neutron Scattering experiments (SANS). Thus, the micro arrangements of microemulsions have been analyzed using the pair correlation function g(r) and the structure factor S(q) obtained from HNC, SANS, and MD simulation using these parameters. The microemulsion dynamic properties were discussed using the mean-square displacement (MSD) and the diffusion coefficient Dc calculated from MD simulations.


Introduction
Microemulsions possess interesting physicochemical properties, thermodynamic stability, low interfacial tension, low viscosity, transparency, and high solubilization power. 1-3 From an industrial viewpoint, microemulsions have a simple and lowcost preparation process. The properties of microemulsions make them useful for transdermal drug delivery; 4,5 they can provide rapid drug release due to the large interfacial area and this improves the efficacy of the drug. Also, they can solubilize both hydrophilic and lipophilic drugs. Usually, microemulsions comprise two immiscible liquids and a mixture of surfactant and cosurfactant. High amounts of surfactant and cosurfactant provide an enormous reduction in interfacial tension, 6 enabling a large interface and supporting the spontaneous formation of microemulsions. The low interfacial tension ensures excellent surface contact between the skin and the microemulsion and facilitates further skin penetration. Surfactants decrease the interfacial tension, thus facilitating dispersion during the preparation of microemulsions. They form a exible lm which provides the correct curvature for the interfacial region. Hydrophilic surfactants were found to increase percutaneous absorption whereas less hydrophilic surfactants increased skin drug deposition. 7 The cosurfactant provides the exibility required by the interfacial lm to form microemulsions with different curvatures. Short to medium alcohol chains are used as cosurfactants to decrease the interfacial tension and increase uidity. [8][9][10] Thus, microemulsions may optimize skin drug targeting and reduce further systemic absorption.
From a mesoscopic point of view, much effort has been dedicated to the development and parameterization of effective interaction potentials to study the structure and dynamics of colloidal systems. [11][12][13] This is supposed to save extensive details and yield valuable information. Thus, the physical properties, in terms of the micro arrangements and collective dynamics of colloidal systems, are determined by considering the colloids as dispersed particles interacting via an effective interaction. In this context, a series of reports have studied the theoretical modeling and numerical simulation of microemulsions, aiming, in particular, to determine their structure and dynamic properties using MD simulations, 12,13 using an effective interaction of the Sogami-Ise type. This work could make a substantial contribution to determining the dynamic laws of Pickering emulsions. However, it has not yet been conrmed by enough experimental results. In the present study, we have modeled the interactions between positively charged microemulsions using hard sphere, van der Waals, and electrostatic repulsion potentials. Fig. 1 shows a schematic representation of the main research steps. Accordingly, the proposed interaction potential model is tested by a numerical solution of the Ornstein-Zernike (OZ) integral equation with the Hypernetted Chain (HNC) approximation. 14-18 Then, we verify the validity of the potential parameters with molecular dynamics simulation and compare the results of MD, HNC approximation, and SANS experiments for the structure factor S(q). Furthermore, MD is used to determine the thermodynamical properties of the system, a microemulsion of decane in salt-water stabilized by a monolayer of cetylpyridinium chloride ionic surfactant and octanol cosurfactant.

Materials
Cetylpyridinium chloride [H 3 C-(CH 2 ) 15 ]-C 5 H 5 N + Cl À (CpCl) was purchased from Fluca and was puried by successive recrystallization in water and in acetone. Octanol [H 3 C-(CH 2 ) 7 ]-OH and decane [H 3 C-(CH 2 ) 8 CH 3 ] were obtained from Fluca and were used as received. All samples were prepared by weight in brine or deuterated brine. Brine (0.2 M NaCl) was prepared with triply distilled water or deuterated water according to the Material Safety Data Sheet (MSDS).

Microemulsion preparation
The surfactant lm was composed of a CpCl cationic surfactant and an octanol cosurfactant. The weight ratio of octanol to CpCl was 0.25 and the weight ratio of decane to the surfactant lm was 0.62. The samples were shaken and kept in a water bath for ten days at a constant temperature. Aer this time, phase separation was clearly observed. Then, the samples were homogenized, le to rest, and observed for four days to ensure their thermodynamic stability.
The weight ratio of octanol to CpCl is chosen to be just below the emulsication failure limit, which refers to the limit above which the microemulsion droplets are saturated with oil and coexist with excess oil. Under such conditions, the microemulsion droplets are spheres with a well-dened radius. [19][20][21] The prepared microemulsions can be diluted over a large concentration range (1-20 wt%). In order to properly study the structural and dynamic properties at different microemulsion concentrations, we performed the study on more than ten samples (1.4%, 2.8%, ., 26%). Aer analysis of these samples, we noticed that the results for the volume fractions 2.8%, 6.98%, 14%, and 25.5% describe the behavior of our system well.

Small angle neutron scattering (SANS) measurements
The SANS experiments were performed at LLB-Saclay on the PACE spectrometer. The congurations used were 1.5 m at 6 A and 4.68 m at 13 A. The range of scattering vectors covered was 0.004 A À1 < q < 0.16 A À1 and the temperature was T ¼ 293 K. The scattering data are treated according to standard procedures. They are put on an absolute scale by using water as a standard. And the obtained intensities in absolute units (cm À1 ) are found with accuracy better than 10%. All spectra models are convoluted by the instrument response function, taking into account the uncertainty of the neutron wavelength and the angular denition. 22 In order to separate the form factor P(q) from the structure factor S(q), it is necessary to measure the scattering intensity I(q). For a nite microemulsion concentration, the scattering due to one particle P(q) is isolated. Thus, the form factor, P(q), is found. Then the scattering intensity is measured in order to nd the structure factor, S(q). The size of the microemulsion is determined using SANS experiments with sphere contrast. For all studied samples, the scattering is essentially due to the microemulsions. Thus, the scattered intensity I(q) is modeled using the following expression: 23 IðqÞ ¼ APðqÞSðqÞ where S(q) is the structure factor, q ( A À1 ) is the scattering vector, F is the volume fraction of the aggregates, v (cm 3 ) is the dry volume of the aggregates, Dr ¼ 6.83 Â 10 10 cm À2 is the contrast, i.e., the difference between the aggregate and the solvent scattering length density, and P(q) is the form factor of the colloidal aggregates, such that P(q / 0) ¼ 1.

Structure and thermodynamics from the integral equation method
The physical properties of microemulsions are essentially determined by the interactions between particles which are in most cases modeled by a pair potential. However, other contributions should be taken into account. The central problem in liquid theory is the determination of the radial distribution function g(r) of the system particles. The radial distribution function is important for several reasons; its Fourier transform is the structure factor S(q), which can be measured experimentally by neutron scattering, thus allowing a direct comparison between theory and experiment. The g(r) is related to the mean interaction potential by: The average effective potential J(r) includes not only the effect of direct interactions U(r) between a particle located at the origin and other surrounding particles but also the average effect of the forces of the neighboring particles which surround each particle u(r). The structure factor S(q) is the Fourier transform of the radial distribution function g(r), which describes the three-dimensional distribution arrangement of microemulsions, where The isothermal compressibility c T is related to the structure factor S(q) at q / 0: The thermodynamic quantities of the microemulsion depend on g(r) and U(r). We write the rst relationship giving the internal energy as a summation of the kinetic and potential energies: Also, the system pressure can be expressed as a function of g(r), by connecting the 2nd virial coefficient to g(r):

Solution of the Ornstein-Zernike integral equation
The second step consists of determining the RDF by solving the well-known Ornstein-Zernike (OZ) integral equation. The latter is solved by the total-correlation-function, h(r) ¼ g(r) À 1, 24 Here c(r) represents the direct-correlation-function. Since the considered system is homogeneous and isotropic, the associated functions g(r), h(r) and c(r) depend only on the particle distance, r, where one particle is taken as the origin. Then, r is the position-vector of the second particle. Fourier transforming the above integral equation and taking into account relation (3) yields the following general formal expression for S(q).
wherec(q) is the Fourier transform of the direct-correlationfunction, c(r). The OZ equation (eqn (7)) contains two unknown functions which are c(r) and h(r). Therefore, to get h(r) (and thus g(r)), we need an additional equation (closure relation) connecting h(r) and c(r) to the pair potential, U(r). Within the framework of HMSA, the adopted closure relation is as follows 24 with the associated bridge-function Generally, such a function is dened as an innite series of correlation diagrams. We note that HMSA interpolates between "So Core Mean Spherical" and HNC approximations. In eqn (9) and (10), the function G(r) is the difference between the total and direct-correlation-functions, i.e., G(r) ¼ h(r) À g(r), and the interaction potential, U(r), is divided into a short-range part, U 0 (r), and a long-range attractive tail, U 1 (r), as prescribed by Weeks et al. 25 We have with and Here r m denotes the position of the minimum of the overall potential, U(r). The above decomposition was proposed by Weeks, Chandler, and Anderson. 25 In eqn (10), the quantity f(r) represents the mixing-function. 26 A new form of this function was suggested by Bretonnet and Jakse. 27 Its advantage is that it ensures thermodynamic consistency in the calculation of the isothermal compressibility by two different ways. These authors have proposed the following form for the mixing-function where a À1 is some characteristic length and f(0) is the interpolation constant. The latter is an adjustable parameter such that 0 # f 0 # 1. This constant that serves to eliminate thermodynamic inconsistency can be xed by equating the compressibility deduced from the virial pressure, using eqn (6), to that calculated from the zero-scattering-angle limit of the structure factor using eqn (4).

Molecular dynamics simulation (MD)
2.6.1 Interaction potential. Microemulsions are approximated as non-overlapping hyperspherical uniformly charged particles, interacting by an effective potential based on the DLVO model, 28,29 which is the sum of a hard sphere repulsive potential, a van der Waals attractive term 28 and a Yukawa type screened coulombic repulsive potential: 29 with: This journal is © The Royal Society of Chemistry 2020 RSC Adv., 2020, 10, 36155-36163 | 36157 Paper RSC Advances where r is the center-to-center distance and R is the microemulsion radius. A H is the Hamaker constant, with A H ¼ 1.1k B T; this value is appropriate for decane microemulsions interacting through water. 30 R c is the cut off distance for the interaction potential, R c ¼ 12R.
where l B is the Bjerrum length, T is the temperature, and k À1 represents the Debye length, k À1 ¼ 6.7 A, which corresponds to the concentration of small ions added to the solution. Z eff is the number of effective charges per microemulsion, Z eff ¼ 130 for each microemulsion. Fig. 2 presents the interaction potential plotted using these parameters. 2.6.2 Simulation details. Molecular dynamic simulations are performed in the NVT statistical ensemble, where N ¼ 10 6 is the number of simulated microemulsions, V is the volume of the simulation box and T is the temperature. Periodic boundary conditions in three dimensions are applied to remove the surface effect. 31 Dimensionless units are used to optimize the calculation time. 31,32 Accordingly, the distance is reduced by the microemulsion radius R, the energies are reduced by k B T, the time is reduced by R ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi ðm=k B TÞ p ; the pressure P is reduced by k B T/R 3 , the isothermal compressibility c T is reduced by R 3 /k B T, the diffusion coefficient D c is reduced by R ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi ðk B T=mÞ p ; the viscosity h is reduced by ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi ffi ðððk B TÞ=mÞ=R 2 Þ p and the dimensionless temperature is set to T* ¼ 1. The equations of motion are integrated using the Verlet algorithm 33 for one million steps with a time step dt ¼ 0.001. All molecular dynamics simulations were performed using the LAMMPS simulation package (Fig. 3). 34 3 Results and discussions

Microemulsion shape
In all studied cases, the polydispersity of the microemulsion size is described by a Gaussian distribution with mean radius R, and standard deviation DR, assuming that microemulsions are spherical particles. Then, the average form factor of the polydisperse size of the studied microemulsions is: with the form factor Here, R is the microemulsion radius. The Porod representation (q 4 $I(q)) as a function of q is used in order to amplify the P(q) oscillations, in the limit of large q. 35,37 According to eqn (20) and Fig. 4, for the rst and second maxima we nd qR ¼ 2.73 and qR ¼ 6.12, and for the rst and second minima, we nd qR ¼ 0 and qR ¼ 4.99, respectively. Thus, for weakly polydisperse microemulsions, R has been estimated using eqn (19). The re-normalized spectra are examined in the Porod representation (q 4 $I(q)/F) for different volume fractions ranging from F ¼ 2.8% to F ¼ 26.5%. All spectra display well-dened oscillations at intermediate and large q for all examined cases. These spectra are superimposed in the large q range, which indicates that the size and shape of the microemulsion remain constant for various F. This behavior is observed in  many experimental reports investigating the microemulsion structure using the Porod representation. [35][36][37] Assuming that decane/water microemulsions are spherical shaped, their radius is deduced from the maximum and minimum positions of the Porod representation. To investigate the size polydispersity, we assume a Gaussian distribution of the microemulsion radius. The spectra calculated using eqn (19) and (20) give a reasonable t with mean radius R ¼ 62 A and polydispersity p ¼ DR/ R ¼ 10%. For simplicity, the microemulsions are considered as monodisperse spheres of mean radius R. 38 The g(r) is calculated at various volume fractions (F ¼ 2.8%, 6.98%, 14%, 26.5%). For the dilute systems (F ¼ 2.8%, 6.98%), the depicted g(r) curves obtained from HNC and MD are in good agreement ( Fig. 5A and B). For high concentrations ( Fig. 5C and D), the HNC approximation slightly underestimates the principal peak height, g(r m ), at the position r m , and the discrepancy between HNC and MD is noticeable for F ¼ 26.5%. The microemulsion correlation increases on increasing the volume fraction F from 2.8% to 26.5%, and the height of the rst g(r) peak increases from 1.07 to 2.36. Furthermore, the position of the principal peak is shied to smaller values. This indicates that the interactions between microemulsions become more repulsive and the organization of the suspension structure is more ordered. Thus, the separation distance between microemulsions is shied to smaller values, as can be observed in Fig. 5. Fig. 6 presents a comparison of the structure factors S SANS (q), S HNC (q), and S MD (q) obtained from the SANS experiment, HNC approximation, and molecular dynamics simulation, respectively. An excellent agreement is found for the studied spectra, proving that the proposed potential model is satisfactory for all studied volume fractions. The slight increase in the rst peak height for F ¼ 14% and F ¼ 26.5% does not affect the t quality as long as molecular dynamics simulation reproduces the peaks with the same q max . We notice that the principal peak becomes narrower and more pronounced with increasing F. Indeed, going from 2.8% to 26.5%, S max varies from 1.012 to 1.833. On the other hand, the position of the main peak is shied to higher values when F increases. This fact is not surprising because q is inversely proportional to the mean distance between microemulsions d (see Table 1). The present ndings in terms of g(r) and S(q) are in good agreement with the results of other experimental and numerical approaches for similar colloidal systems. [11][12][13][39][40][41] The results presented in Table 2 indicate that the pressure and the potential energy increase monotonically with increasing F, and the osmotic compressibility decreases with increasing F. Fig. 4 Renormalized spectra in Porod representation (q 4 $I(q)) at different volume fractions (F ¼ 2.8%, F ¼ 6.98%, F ¼ 14%, and F ¼ 26.5%). The solid line is the computed spectrum for R¼ 62 A and DR ¼ 6.2 A.  where m stands for the mass of the microemulsion, v(t) stands for its velocity, z stands for the friction coefficient, and F s (t) represents a random force. The friction coefficient z is related to the microemulsion viscosity, h, and the microemulsion radius, R, by the classical Stokes relation z ¼ 6phR. The viscosity of the solution depends on the number density of the dispersed

Static properties
: For low volume fractions, the effective viscosity obeys the well-known law 42 Here n accounts for the reduced microemulsion number density, h M is the solvent viscosity, and h N is the viscosity of the microemulsions. We have h(n) > h M , since h N > h M . The stochastic force is considered to be white noise with: where d(t) is the Dirac distribution. The brackets h.i mean an average over time. The mean-square displacement (MSD) h(r À r 0 ) 2 i h hD 2 ri is expressed as: where g ¼ z/m is known as the relaxation rate. For times much longer than the inverse relaxation rate, that is for t [ g À1 h s, the MSD increases linearly with time, and it is dependant on the diffusion coefficient D c , dened as the compromise between the thermal uctuations and the dissipation such that D c ¼ k B T z ; according to: For t ( s, the Taylor expansion of e Àgt gives: e Àgt z 1 À gt þ 1 2 ðgtÞ 2 ; and then we nd: 3.3.2 Theoretical prediction of the microemulsion dynamic properties for high volume fractions. When the number density of the microemulsions reaches a critical density n c , their dynamics are more complex than those at low number density. We recall that, for low density, particles move ballistically for short times. Thus the mean-square displacement hD 2 ri $ t 2 , which is followed by a crossover to Fickian diffusion, characterized by hD 2 ri $ t for long times. But, in dense colloidal systems, a caging effect, where the particles are trapped by their neighbors, causes a subdiffusion regime intermediate between the ballistic and the diffusive regimes. This regime is characterized by hD 2 ri $ t a , with 0 # a # 1. For the time being, the theoretical prediction of the value of a in the presence of the cage effect remains poorly understood. But, in fact there are some theoretical approaches based on some complex memory functions developed to discuss the subdiffusion laws. From the literature, Zwanzig, in a series of studies, developed a theoretical approach based on a generalized Langevin equation (GLE), 42,43 and this approach can be adopted to discuss the subdiffusion observed in microemulsions which can be considered as three dimensional colloidal particles swimming in a solvent. From a mathematical viewpoint this is simply an extension of the standard Langevin equation for simple liquids, where the friction is assumed to be determined by the instantaneous velocity of particles. It is noted that the difficulty of understanding the subdiffusion phenomenon lies in mathematically representing the cage effect that, physically speaking, depends on various physical parameters, namely, the temperature, the number density and the nature of the solvent. In this context, the GLE is expressed as where v is the velocity of a moving microemulsion, g ¼ z m is the relaxation rate, where z is the friction coefficient and m is the mass of the tracer, k is the memory function that expresses the retardation due to friction, and F s (t) is a random force felt by the moving microemulsion due to its collisions with the solvent molecules. Thus, the random force is obtained from the equation For these considerations, the VACF solves the following differential equation: To solve the GLE, we adopt a memory function, recently proposed by Flenner, developed to study the dynamic aspect of lipid atoms. Flenner's approach is based on the Zwanzig-Mori projection method for modeling hD 2 ri. The starting point is to consider the equation of motion for the density autocorrelation function F s (q,t) ¼ hn(Àq,0)n(q,t)i of a selected particle at wave vector q, where v ¼ ffiffiffiffiffiffiffiffi k B T p =m; k B stands for the Boltzmann constant, T denotes the temperature and m is the microemulsion mass.
The equation of motion for the MSD can be obtained from where k 0 ðtÞ ¼ lim Mðq; tÞ: Thus, Flenner et al. propose a memory function with explicit crossover time scales which reproduce the three diffusion regimes, namely, ballistic, normal, and subdiffusive 44 where d(t) is the Dirac delta function, B is a dimensionless parameter, s 1 is the characteristic time for the crossover from subdiffusion to normal diffusion, s 2 is the onset time of the subdiffusion regime, and s 3 is the characteristic time for the crossover from the ballistic regime to the subdiffusive regime. Analysis of this memory function revealed that, for B ¼ 0, the Brownian diffusion laws are recovered and hD 2 ri behaves as The subdiffusion phenomenon is represented mathematically in this memory function by the power law. The analytical determination of the diffusion exponent a seems to be complex. On the other hand, a numerical computation, taking into account the power law, 44 revealed that the dynamics pass slowly from the ballistic regime to the subdiffusive regime with a < 1. The region of subdiffusion separating the ballistic and the normal regimes is due to a cage effect. In this region the meansquare displacement, hD 2 ri, behaves as hD 2 ri $ t a (a < 1), t > s 2 (37) 3.3.3 MD simulation results. The inuence of the volume fraction on the dynamic properties of decane/water microemulsions is investigated using NVT molecular dynamics simulations. For this study, the microemulsion radius, the effective charge of the microemulsion, and the temperature are xed to the values: R ¼ 62, Z eff ¼ 130, and T* ¼ 1. But, the volume fraction varies over a large range. Fig. 7 shows the reduced MSD, Dr CM 2 (t), as a function of reduced time, t* ¼ t=ðR ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi ðm=k B TÞ p Þ; on a log-log scale, for many values of volume fraction. For early times, we note that the curves are superimposable and the random-walker exhibits ballistic diffusion, and we have Dr CM 2 (t) f t 2 . Thus, the microemulsion diffusion is practically independent of the volume fraction. For later times, we note that the MSD depends on F. Accordingly, for low and medium volume fractions (F ¼ 2.8%, 6.98%, 14%, 26.5%), for a reduced time below the transition time (s*), the microemulsions experience ballistic diffusion. But aer this transition time, the microemulsions exhibit normal diffusion. The obtained results are in good agreement with the theoretical predictions discussed previously (eqn (26) and (27)). This can be explained by the fact that at early times, the random-walker, referring to an individual microemulsion, does not feel the presence of its close neighbors yet. This diffusion regime is well known as Fickian diffusion. 32,45 For high density microemulsions, ballistic diffusion is also observed for early times. However, for intermediate times beyond s * 2 ; the tracer, representing the considered microemulsion, is trapped in a cage formed by its close neighbors and then moves slowly. The associated MSD then behaves as Dr CM 2 (t) f t a , with subdiffusion exponent a ¼ 0.47. This result conforms to the theoretical predictions we described above (eqn Fig. 7 Log-log plot of MSD versus time step from MD simulation. This journal is © The Royal Society of Chemistry 2020 RSC Adv., 2020, 10, 36155-36163 | 36161 (37)). If the number of surrounding neighbors in the cage where the tracer is trapped is denoted by N c , the stay duration in this cage, Dt c , scales as Dt c fN c 1 a: Therefore, the stay duration in the cage becomes longer as the volume fraction is increased, because of the attraction that favors contact between the microemulsions and then blocks their movement. For high volume fractions, the available space for the microemulsions becomes smaller because of the strong attraction. This attraction then leads to slowing-down of the microemulsion scattering, which is localized for a certain time (stay duration in the cage). The latter increases with the volume fraction. Aer the subdiffusive regime, for t . s * 2 ; the long-time diffusive regime (where Dr CM 2 (t) f t a , a ¼ 1) is established. Generally, the ballistic and subdiffusive regimes do not last for a long time and a crossover to normal diffusion is always observed. 46,47 In Table 3, we give the normal diffusion coefficient, D * c ; for the studied microemulsions as a function of the volume fraction. We observe that, when F increases, D c decreases exponentially. This can be explained by the fact that as F increases, the correlations between the microemulsions become larger. Similar dynamic behavior has been recently reported in MD simulations of a large number of colloidal particles. [48][49][50][51][52] According to Fig. 8, D * c decreases following an Arrhenius-like law, where: with D * 0 ¼ 88:67 and F 0 ¼ 0.0156. This behavior was observed in a recent MD simulation of Pickering microemulsions. 13

Conclusion
In this work, we investigated the static and dynamic properties of positively charged spherical decane/water microemulsions using SANS experiments, the HNC numerical method, and NVT molecular dynamics simulations, for microemulsions with different volume fractions F. The calculated properties from the different approaches are found to be in good agreement, which conrms the appropriate choice and the correct parameterization of the interaction potential. It is noted that the static and dynamic properties of microemulsions strongly depend on F. On the one hand, the interaction between the microemulsions increases proportionally to F. On the other hand, the average distance between the microemulsions, the isothermal compressibility, and the diffusion coefficient decrease with increasing F.
It should be noted that the determination and parameterization of an effective interaction potential for studying colloidal particles is a real challenge. Although the DLVO model is complete, the parameterization of this potential model requires, at least, an experimental adjustment in terms of structural properties. Thus, the present procedure used to parameterize the interaction potential between the decane/ water microemulsions stabilized by the CpCl surfactant and the octanol cosurfactant can be carried out for different types of colloidal particles. Then, molecular dynamics simulations can be used to determine several physical properties based on a proven effective interaction.
Finally, we propose that the present study can be extended to study microemulsions with biocompatible gra polymers on their surfaces, assuming that the incorporation of polymers can increase their stability and their diffusivity and thus increase the efficiency of microemulsions as drug carriers. This hypothesis will be the subject of our next investigations.

Conflicts of interest
There are no conicts to declare.   Fig. 8 Diffusion coefficient D * c as a function of the volume fraction F.