Molecular dynamics study on surface adsorption of hexyl trimethylammonium bromide, octyl trimethylammonium bromide,...

Molecular dynamics simulations were performed on four surfactants (hexyl trimethylammonium bromide, octyl trimethylammonium bromide, sodium hexyl sulfate, and sodium heptyl sulfate) to investigate their surface adsorption behavior in aqueous solution. In this work of simulations, surfactant molecules were initially distributed in both the surface and bulk regions of the solution, enabling them to form micelles spontaneously within the aqueous solution. The results show that the solute surface concentration and bulk concentration for both cationic surfactants (hexyl trimethylammonium bromide and octyl trimethylammonium bromide) and anionic surfactants (sodium hexyl sulfate and sodium heptyl sulfate) can be accurately described by a Padé approximant, which is consistent with the results of our liquid surface adsorption theory. The type of hydrophilic head groups and the length of hydrophobic tail groups affect the speed of surface adsorption. Molecular dynamics simulation demonstrated that the thickness of the surface layer in surfactant solution increases as the temperature and concentration increase. This is due to the weak intermolecular interaction between solvent molecules.

Molecular dynamics simulations were carried out in the NVT ensemble using the GROMACS software package. The interaction potential was chosen to be the Lennard-Jones 12-6 model, with long-range electrostatic interactions treated using the Particle Mesh Ewald (PME) method. A cutoff distance of 1.0 nm was taken for both non-bonded and electrostatic interactions. The V-rescale thermostat was used to realize a thermal equilibrium state at 298.15 K. The time steps were all taken as 2 fs. A typical total simulation time was 35 ns, of which the initial period was used to bring the system to equilibrium, and the remaining steps were used to dump the trajectories for subsequent analysis via VMD.