Harmonic maps with potential from complete manifolds are considered. This is a new kind of maps more general than the usual harmonic maps relating to many interesting problems such as equilibrium system of ferromagnetic spin chain and Neumann motion. Aiming at the general properties, the author derives basic gradient estimates and then Liouville type results for these maps, which are interesting in constrast to those of the u...