Molecular dynamics (MD) simulation has been applied to study of pull-out of Halloysite nanotubes (HNTs) from a polyurethane (PU) matrix. First, the Machine learning (ML) particle swarm optimization (PSO) method was used to obtain force field parameters for MD from data of density functional theory (DFT) calculations. The current study shows the possibility of using a PSO technique to modify the force field with DFT data with less than 5 kcal/mol discrepancy. Second, we considered the influence of atomic interface on pulling out of HNT from PU. Energy variation has been proposed as the cohesion strength between matrix and nanoparticle. In addition, the best Lennard Jones parameters in the MD simulation make good agreement with an experimental sample stress-strain response.