Returning to the details in performing the simulations, let us first emphasize that the radial distribution function g(r) was calculated for both the equilibrium and nonequilibrium systems using the program RDF28 for each atom. This method characterizes the hydration behavior of each species and identifies the total number of H-bonds and the hydration number of any component . In the radial distribution profiles, the characteristic peaks corresponding to the HB distance signifies whether the hydrogen bond forms between molecules of interest, namely dodecanoic acids, water and poly(oxyethylene) ether. The coordination numbers (CN) were estimated for each species by integrating the radial density function, up to a value that corresponds to its first minimum.