Molecular　dynamics　simulation　is　made　to　study　configuration　and　total　energy　of　0.5mol%,　2.1mol%,　3.8mol%　VB　transition　metal　ions　(V5+,　Nb5+,　Ta5+)/TiO2　in　rutile　at　300K　and　101325　Pa.　As　2.1mol%　Ti4+　is　substituted　by　dopants,　configuration　remains　well　with　less　mean-square　displacements　(MSDs),　distinct　planes　of　atoms　and　higher　stability.　Nb5+　incorporates　into　TiO6　octahedra　well　and　shows　a　large　solubility　in　TiO2　compared　with　V5+　or　Ta5+　due　to　comparable　valence　and　ionic　radius　between　Nb5+　and　Ti4+.