This study explores the interaction between particles in microplastic semi-coke water slurry at the molecular level using molecular simulation methods, specifically DFT calculations and MD simulations. In addition, the experiment of slurry preparation was carried out to study the viscosity and stability of the slurry. The electrostatic potential analysis shows that the interaction between microplastics and dispersant molecules occurs on atoms with large electronegativity or oxygen-containing functional groups, and the energy gap of frontier molecular orbitals indicated that PVC interacts most easily with the dispersant (0.39eV), followed by PS (1.08eV) and PET (3.65eV). In addition, it is also noted that due to the steric hindrance effect, the adsorption energy was opposite to the DFT calculation results: PET was - 213.338kcal/mol (NNO) which was highest, followed by PS (- 107.603kcal/mol, NNO), and PVC (NNO) was lowest which was - 94.808kcal/mol. And RDF shows similar results, which the probability of water molecules in the PET system was the highest, followed by PS, and finally, PVC. The MD results are consistent with the viscosity and stability characterization results of the slurry which PET has the lowest viscosity of 87.3mPa·s. Finally, this study provides new ideas for the treatment of microplastics and the improvement of the performance of semi-coke water slurry and reveals the interaction mechanism between microplastics and semi-coke water slurry. All calculations were performed using Materials Studio (MS) version 2020 software, BIOVIA Corporation. The DFT calculation was carried out through the DMol3 module. The DFT calculations include electron density, electrostatics, orbitals, and population analysis. In DMol3 module, the GGA-PBE function was selected to consider gradient changes in density in the simulated calculation. The DFT-D correction was selected, and all electrons were calculated by DNP for accurate core potentials and the DNP file was 4.4. MD simulation was performed through the Forcite module. MD simulation mainly focuses on relative concentration distribution analysis, radial distribution function, and adsorption energy calculation. All molecular geometry optimizations are performed in the Forcite module. In the molecular dynamic part, all simulations used PCFF forcefield. The NVT ensemble was adopted and using the Nosé thermostat.