A quark coalescence model, based on semi-relativistic molecular dynamics with color interactions among quarks, is presented and applied to pp collisions. A phenomenological potential with two tunable parameters is introduced to describe the color interactions between quarks and antiquarks. The interactions drive the process of hadronization that finally results in different color neutral clusters, which can be identified as hadrons based on some criteria. A Monte Carlo generator PYTHIA is used to generate quarks in the initial state of hadronization, and different values of tunable parameters ...