Monte Carlo simulation of 3-star lattice polymers pulled from an adsorbing surface

We study uniform 3-star polymers with one branch tethered to an attractive surface and another branch pulled by a force away from the surface. Each branch of the 3-star lattice is modelled as a self-avoiding walk on the simple cubic lattice with one endpoint of each branch joined at a common node. Recent theoretical work (Ref.~1) found four phases for this system: free, fully adsorbed, ballistic and mixed. The mixed phase occurs between the ballistic and fully adsorbed phase. We investigate this system by using the flatPERM Monte Carlo algorithm with special restrictions on the endpoint moves to simulate 3-stars up to branch length 128. We provide numerical evidence of the four phases and in particular that the ballistic-mixed and adsorbed-mixed phase boundaries are first-order transitions. The position of the ballistic-mixed and adsorbed-mixed boundaries are found at the expected location in the asymptotic regime of large force and large surface-monomer interaction energy. These results indicate that the flatPERM algorithm is suitable for simulating star lattice polymers and opens up new avenues for numerical study of non-linear lattice polymers.