The work attempted to combine two forms of diffusion starting with phenome- nological approach motivated by Fick’s laws and their mathematical solutions and the “random walk” of diffusing particles suggested by Robert Brown. This random walk of microscopic particles in suspension in a fluid has since been adopted as the “Brownian motion” in honour of Brown. The mathematical form of the Brownian motion was derived by Einstein  . The Fick’s laws follow from Fourier’s law of thermal conduction     .
The research encouraging modelling of random fluctuations found in PK/PD (Pharmacokinetic and Pharmacodynamic) relations has been suggested. It has been noted that there is an increasing need to extend the deterministic models which are currently favoured to models including a stochastic component in modelling pharmacological processes  .
The transport diffusion is the concentration gradient dependent driven movement and describes the change in concentration of a solution particle. The transport diffusion can further be studied from Fick’s laws by use of partial differential equations, however reseachers have used different variable spaces to study this flow    . Furthermore, in this work we tracked the transport diffusion by making use of the one to one function relationship between secondary saturation and concentration  . It is also noted that concentration is not a movement parameter as compared to secondary saturation. The self diffu- sion is a concentration gradient independent driven movement and describes the non-interacting movement of solution particle in the volume space. The observation of self diffusion component has been attributed to Brown  .
The work has allowed for a proposition on the possibility to study the two movements in the 24 h dosing period that aid the process of diffusion by considering stochastic differential equations   . It describes the movement of an interacting particle. Inferences were proposed on the possible relationship between the two movements (interacting and non-interacting). Numerical solutions on diffusional constituent movement characterisation were derived for a patient on efavirenz. This work presents a stochastic coupled diffusion model. This enables us to study the relationship between transport and self forms of diffusion. Stochastic coupled models have been applied on combined PK/PD processes  . In addition, in this case the coupling is in modelling the relation- ship of the two separate diffusional movements. A measure is developed to estimate the extent of gradient driven diffusion with the aid of a coupled stochastic differential equation. The Verhulst model is proposed to estimate the logistic relation obtained from modeling the self and transport forms of diffusion  .
Simulated projected data on secondary saturation movement, time and concentration was taken from pharmacokinetic projections made on Patient P on 600 mg dose considered in Nemaura   . Stochastic differential equations, Ordinary differential equations and linear and non-linear regression models are used to model diffusion in plasma. The following softwares, MATLAB and R were used for development of models. The solution of the Stochastic Differential Equations are approximated with the Euler-Maruyama scheme.
2.1. Deterministic Models for Concentration and Gradient Driven Movements
Initially there was modelling of the concentration profile of a solution particle of patient P in time. Letting model concentration of solution particle at time t and was given by,
for some . The process was sufficiently modelled by the following equation,
We modelled the gradient driven movement by initially considering, the derived secondary saturation movement , for some patient P, in Nemaura (2014) from relations found which were informed from Nemaura (2015). The diffusant was the drug efavirenz in blood plasma. We gave the estimation of the secondary saturation movement in patient P with respect to time and concentration,
The following was noted,
2.2. Stochastic Model
The secondary saturation movement followed a one to one relation with concen- tration and was thus used as compared to other forms of secondary movement that is convection and advection  . The process modelled both the interaction and non-interacting diffusion movement of the solution particle in the volume space. The process modelled these two diffusion forms of the solution particle and was thus given by,
where the process modelled the independent standard Brownian motion (Weiner process), was the relative acceleration drift of a solution particle that allowed for the exchange of concentration material to its neighbourhood, -auxilliary concentration gradient dependant driven movement, -auxilliary concentration gradient independant driven movement and -as the random term in . Additionally,
where at for . The parameter, m is the least upper bound of drift for possible positive gradient driven process to occur. Thus the following:
Furthermore, for and and and some the following conditions are satistified,
Conditions 13 and 14 allowed for the existence and uniqueness of solutions for Equation (11)  .
3.1. Modelling Concentration and the Concentration Gradient Driven Movement in the Absence of Self Diffusion
We estimated the parameters in Equation (1) for the concentration-time curve (Table 1).
Table 1. Parameter estimates in modelling deterministic solution particle concentration.
Table 2. Parameter estimates in modelling saturation movement rate with respect to t (Model (3(i))).
Table 3. Parameter estimates in modelling saturation movement rate with respect to x (Model (3(ii))).
There was consideration of how the concentration and movement processes informed potential relationships. The relationship of the concentration gradient driven movement and concentration rate of change was established by observing the relationship between equation 1 and 7. The following was initially considered where and with and There were two types of movement considered that were postulated to contribute to diffusion. The part modelled the auxilliary transport constituent of diffusion and measures the auxilliary self aspect of diffusion. A high correlation value was obtained of (estimated the extent of transport diffusion) in the 24 h dose interval, where were solutions of Equations (1) and (7) (Figure 1).
Informed by results herewith there was an adoption that transport diffusion was driven by concentration. Thus concentration could be used as a substitute parameter for transport diffusion.
3.2. Stochastic Model: Modelling the Solution Particle Movement with Variable Relative Acceleration (Drift)
The following case was considered where and However, it was projected that at the solution particle movement had zero correlation with the concentration process . In the interval where , a solution particle movement was projected to have no transport diffusion. A solution particle with increasing relative acceleration drift to its neighbour would have less magnitude of extent with respect to transport diffusion. However, a solution particle with low relative drift to its neighbour had a high potential of transport diffusion. That is transport diffusion was more pronounced in neighbouring particles with relatively low acceleration drift. A relationship of the correlation between and , and the drift in self diffusion was proposed (Figure 2).
An equation for the logistic decay relationship  was proposed,
where k and l (decay) are constants (Table 4). We obtained an analogous parameter to the Malthusian parameter and had a value of
where m generally defines the `carrying capacity’ however in this case it was the value of supremum of acceleration drift that enables transfer of concentration material  .
Figure 1. (a) Projected concentration and solution particle movement(transport diffusion with no self diffusion); (b) projected Estimated solution particle movement from the fitted equation and RMS value of 1.9554.
Table 4. Parameter estimates in modelling solution particle concentration and solution particle movement correlation relationship in relation to relative acceleration drift.
Figure 2. Projected solution particle concentration and solution particle movement correlation in relation to relative acceleration drift with and 1.
It is important to note that the relationship between concentration change and solution particle movement was inferred to be predominantly affected by relative acceleration drift.
Other researchers have been able to correlate self and transport forms of diffusion from experimental observations  . Transport diffusion has been shown to thrive in neighbouring solution particles that have zero relative drift. The transport diffusion in a neighbourhood of a solution particle was projected to be affected by drift. Considering an increased level of self diffusion of a solution particle relative to its active neighbourhood, showed that it will give out less in terms of concentration material to be exchanged through gradient. A logistic relation was projected between this correlation and drift. An elaborate study of the dynamics of the logistic decay curve proposed was done in Gonze, and considering the “Malthusian parameter value” resulted in an unstable steady state and unstable cycles for the system derived for Patient P since it had a value greater than 1  . Relatively slow movement of neighbouring solution particles was projected to enable optimum drug transfer through transport diffusion in the plasma.
Throughout this work, there was an inherent assumption that the particles were already in orientation that enabled transport diffusion. This was because for transport diffusion to occur the particles should be in such an orientation which enables progression of that process  . The extent of the transport diffusion of a solution particle was projected to decrease with increasing relative drift to an active neighbourhood.
The author would like to thank the following; C. Nhachi, C. Masimirembwa, and G. Kadzirange, AIBST and The College of Health Sciences, University of Zimbabwe.