Mike John Edwards (Majid Farzin)*
Received: December 01, 2023; Published: December 14, 2023
*Corresponding author: Mike John Edwards (Majid Farzin), Independent research Physicist, Germany
DOI: 10.26717/BJSTR.2023.54.008506
The problem of polymer brush bilayer under stationary shear is studied by using the DFT, the scaling theory and MD simulations. Both theory and simulations confirm that the shear stress follows the universal power law γ1.46 for the brush bilayers with in- terpenetration and in the absence of the interpenetration, the shear stress scales linearly with the shear rate. It is also revealed that the presence of explicit solvent molecules prevents the brushes to form an interpenetration zone, therefore with explicit solvents the shear stress scales linearly with shear rate. Hence, this study strongly confirms that there is no sublinear regime in the world of polymer brush bilayer, neither by solvents nor by hydrodynamic effects. As long as there is an interpenetration zone, the superlinear regime dominates and in the absence of the interpenetration zone the linear regime dominates. Thus, polymer brushes are not a good candidate for lubrication and all works suggesting that this system is a super lubricant are completely wrong.
Abbreviations: FJC: Freely Jointed Chain Model; PET: Perturbation Expansion Theory; ST: Scaling Theory; MFT: Mean Field Theory; DFT: Density Functional Theory; CSHL: Cold Spring Harbor Laboratory; PBC: Periodic Boundary Conditions
Polymers are linear macromolecular structures that are already known as building blocks of life. They are composed of repeating units of atoms or molecules. Each repeating unit of polymers are known as monomers. The connectivity throughout monomers is established via covalent bonds in which the valence electrons of atoms or molecules are shared together. The covalent bonds are among the strongest bonds in nature. The Brownian motion of monomers which arise from collisions by interstitial water molecules, causes monomers to fluctuate. The fluctuations of monomers makes the whole chainto undergo all possible configurations in the long time measurement. The average chain extension plus other time-averaged quantities are key properties of polymers that physicists seek. The end-to-end distance Re and the radius of gyration Rg [1-3] are two candidates for average chain extension, however, Rg is more accurate.Simple calculations show that
In this system, the steric repulsion among the monomers stretch the chains drastically in the perpendicular direction. The density functional theory (DFT) reveals that the average perpendicular chain extension in perpendicular direction or brush height is
In this article, I approach the PBBs theoretically and numerically. The theoretical part is done by a combination of the DFT framework and scaling theory. The numerical part is done by MD simulations. In the coming two sections, I briefly describe the theoretical part and in one section I describe the numerical simulation part. Finally, in the last section, I present the results and I discuss the results and conclude this research.
One of the useful methods to study mechanics of a many body system is the DFT. In the DFT, the Euler-Lagrange equation of the system is calculated from a grand potential functional which is primarily based on a guess. There are also other methods such as Newton’s equation of motion and Hamilton-Jacobi method. In 1988, Hirz in his Master’s thesis applied the DFT method to polymer brush problem and obtained density profiles as well as brush height. [6] Following the same procedure that Hirz did, the PBBs grand potential functional could be given as follows,
Where n(z) and m(z) are the density profiles of the bottom and up brushes respectively, [8] h the brush height, µ the chemical potential, σ the grafting density and D the distance between flat surfaces. Three Euler-Lagrange equations could be obtained plus one constraint which all are given as follows,
The first three equations say that equilibrium density profiles and height is calculated by setting functional derivatives of grand potential with respect to those quantities to zero. The last equation says that the total number of particles are fixed. To obtain equilibrium n(z), m(z), h and µ that are primarily unknown quantities, the four equations above need to be solved simultaneously because they are coupled to each other. The solutions of the above set of coupled equations are given as follows,
where χ0 and χ1 are volumes defined as follows,
By plotting h as a function of each parameter, it turns out that h follows the following similarity relation,
To calculate the equation of state of the PBBs, all calculated quantities are put back into Equation (1), and the pressure via the Gibbs-Dohem relation Φ = −PV is calculated. Thus, the equation of state of the PBBs is calculated as follows,
where the following volumes are defined,
By plotting the pressure equation, it could be simply inferred that the pressure follows the following similarity relation,
It turns out that, the equation [17] of state scales as (N/D)4 which is very important outcome and violates all previous calculations. [18] It reveals that pressure is strongly sensitive to N and D unlike single brush. The quadratic dependence of PBBs pressure on the degree of polymerization as well as distance between surfaces originates from complex interactions between monomers. The phenomenological arguments approaches a physical system by means of simple as- sumptions and obtains the general behavior of that system. The scaling theory is based on the phenomenological arguments and is capable of solving complex physical problems in a simple way. The problem of the PBB under shear seems to be very complicated at first sight. However, by using the following argument, it gets easier to solve. Upon shearing the PBB, above a critical shear rate, the chains start stretching in the shear direction. This is similar to the equilibrium condition but with chains stretched laterally. We can use this concept to build our scaling theory. I introduce a scaling function in the form of
Again from statistical Physics, it is known that for each system
For simplicity, it is useful to work with a dimensionless quantity W ÿ=τ Rγ instead of shearrate. This dimensionless quantity is called Weissenberg number. Now, in an step-by-step fashion, I calculate the physical quantities both in linear and non-linear regimes. The chain extensions at linear regime are given as follows,
The diffusion constant of center of mass monomer in the Rouse’s dynamics with excluded volume interactions is given as
The viscosity in the linear regime is
Now, I utilize the scaling argument to calculate the physical quantities at non-linear regime i.e. when the shear forces dominate the fluctuation forces. When the shear rate exceeds the critical shear rate
In the case of stress tensor, the scaling argument implies that the shear stress at non- linear regime, scales as N4. This way, I determine the scaling function for stresses as
And finally, the kinetic friction coefficient at non-linear response regime is calculated via dividing shear stress by normal stress,
This is a very important result that we will see in the next sections that is in agreement with MD simulations results.
At high shear rate, the chains continue stretching in the shear direction up to a moment when the top and bottom brushes do not interpenetrate anymore. This is the moment when the normal extension equals half the wall distance. It means that when
At γ˙s the equation of state by which we had build our scaling theory is changed. This is completely different scenario with respect to previous one. So we need to do everything from the scratch. The equation of state in this condition is the equation of state of a single brush which is given as follows,
The equation above is calculated by using the grand potential of a single brush
The mean free time is then calculated as follows,
So we can now calculate the viscosity in this scenario as follows,
Thus the shear force in linear regime is given as follows,
Now, we can use the scaling theory to check shear stress at high shear rates in this scenario. As before, we build an scaling function
The normal stress is obtained by using the fact that it must become proportional to N3 at high shear rate. This leads us to the equation which clearly gives us α = −1/ (2ν + 1). Therefore, the normal stress in this scenario is given as follows,
Hence, the kinetic friction coefficient in this scenario becomes the following relation,
Repeating the same procedure for the chain extensions, we get the following equations
for the shear and normal chain extensions,
The molecular dynamic simulation is based upon calculation of trajectories of each particle in a discrete manner. The trajectories of each particle is calculated from the Newton’s equation of motion. In the present research, I use the Leap-Frog algorithm to discrete the equations of motion of monomers. This scheme calculates trajectories of monomers with a accuracy of the order of magnitude ∆ t4 which is very precise. It also allows for using large time step and faster calculations. In the simulations, the system is composed of two surfaces located at z = 0 and z = D. The surfaces are build by arranging atoms in a 2D simple cubic lattice. The area of the surfaces is Lx = 20σ times Ly = 20σ . The linear polymer chains of degree N are grafted to randomly chosen wall atoms. The periodic boundary conditions (PBC) is employed to mimic the bulk properties in the system. The shear motion is produced by moving the bottom wall atoms in the +y direction and the top wall atoms in the −y direction with speed v/2. The Lennard-Jones potential which make the hard-core excluded volume interactions between particles is given as follows,
And the finitely-extensible-nonlinear-elastic (FENE) potential connects the monomers in the chains is given as follows,
The Langevin thermostat is given as follows,
In the coarse-grained simulation all quantum mechanical degrees of freedom are inte- grated. It is called Kremer-Grest model. [20,21] The simulations have run for 105 time steps and to capture the most accurate results. Table 1 shows the numerical values of simulation parameters in my simulation.
The problem of the PBBs under stationary shear motion is studied by means of the DFT framework, the scaling theory and the MD simulations. In Figure 2(i), the shear force exerted on the top wall in terms of the shear rate is shown. Here, the similarity relation Equation. (9) for the shear force in non-linear regime of the low shear scenario, Equation. (14) for the shear force in non-linear regime of the high shear scenario and the MD simulations results are shown together. As it is seen, the shear force scales as γ˙1.46 in low shear rate scenario where the brushes interpenetrate and it scales as γÿ−0.46 .˙ in high shear rate scenario where the brushes do not interpenetrate. In Figure 2(i), we see a perfect agreement between the DFT and the MD simulations. Each plot shows the DFT and the MD simulations results at low and high shear rate scenarios. In Figure 2(ii), the normal force exerted on the top wall in terms of the shear rate is shown. Here, the similarity relation Equation. (9) for the normal stress of the low shear scenario, Equation (15) for the normal stress in high shear scenario and MD simulations results are shown together. It is shown that the normal stress at both low and high shear scenario scales as . Here again, we see a great agreement between DFT and MD simulations. In Figure 2(iii), the kinetic friction coefficient for both scenarios are shown. Here, the similarity relation Equation. (10), Equation. (16) and the MD results are shown. It turns out that the kinetic friction coefficient at low shear scenario scales as γÿ1.92 and at high shear scenario scales as γÿ1.46 . We find a great agreement between the DFT and the MD simulations. In Figure 2(iv), the shear viscosity in terms of the shear rate is shown. Note that the shear viscosity could be calculated by dividing the shear force in each scenario by the shear rate.
It turns out that the shear viscosity scales as γÿ0.46 at low shear scenario and it is independent of the shear rate at high shear scenario. These universal power laws are well confirmed by the MD simulations as well. In Figure 2(v), we see the shear chain extension in terms of the shear rate. Here, the MD simulations results reveal an anomalous behavior of chains in the shear direction. It turns out that by increasing the shear rate, the shear chain extension increases, then decreases and again increases. Apparently, when we increase the shear rate, there is a critical shear in which the chains have minimum extension in the shear direction. This critical shear rate could be the topic of more researches in the future. And finally Figure 2(vi) shows the normal extension of the chains in terms of the shear rate. Here, we see Equation. (8), Equation. (17) and the MD simulations results together. It turns out that the normal extension at low shear scenario scales as γÿ−0.34 while at high shear scenario, it scales as γÿ−0.18 As it is seen, we again get a great agreement between the DFT and the MD simulations results. Without any complicated calculations or extensive numerical simulations, one immediately infers that when two opposing polymer brushes intermediately compressed such that their chains interpenetrate into each other, the friction between them is larger than simple fluids. The reason is clear the interpenetrating chains produce larger friction.
However, to find out the power laws we must use theoretical tools. So, my theoretical approach predicts that when two interpenetrating brushes are sheared, the shear stress scales as γÿ1.46 , of course, as long as the brushes interpenetrate. This universal power law is very important because it correctly shows that if we use the equation of state of the polymer brush bilayer with interpenetrating chains, and apply the scaling theory to this equation of state, we get the correct power law. This is really amazing because molecular dynamic simulations confirm this power law as well. The second scenario takes place when the chains stretch drastically in the shear direction such that eventually the normal chain extension reaches D/2 and at this moment the interpenetration zone vanishes. At this moment, the equation of state is changed and we have to apply the scaling theory to the equation of state of a single brush or non-interpenetrating brushes which are equal in meaning. When we do the same procedure
on this equation of state, we surprisingly get the power law γÿ for the shear stress. This is really really amazing because both theory and simulations strongly agree on this scenario. In conclusion, the interpenetration between top and bottom brushes play the key role in the behavior of the system. If we take into account the explicit solvent molecules, it is always and in any conditions the solvent molecules concentrate in the middle of brushes and push the brush monomers out of the middle zone. It means that with explicit solvent molecules we never ever have an interpenetration zone between brushes. Note that you cannot reach any interpenetration by compressing the brushes over each other because the solvents are always there. Therefore, with explicit solvents, the shear stress always scales linearly with the shear rate. It means that there is no sublinear regime in the PBB system. Hence, at the end of this article I would say that the works already published [5,13-15] are completely wrong.
It is my honor to acknowledge the Biomedical Journal of Scientific and Technical Research (BJSTR) for supporting my research. Also, I am honored to acknowledge the Cold Spring Harbor Laboratory (CSHL) and the Regeneron Pharmaceuticals Inc. [22-26] for giving credit to my research. I also, thank Weizmann institute for the XMGrace software, the Wolfram Research for the Mathematica software and the GNU Inc. for the GFortran without which my research was not come into reality.