An ''equation of motion'' for matrix elements of an arbitrary Hermitian operator with respect to nuclear coordinates is derived. In the diabatic basis, this equation expresses the smoothness of the corresponding molecular property. Its solution, which determines an adiabatic-to-diabatic transformation, is considered in the two-and three-state approximations. The relation between a smoothness of molecular property and a configurational uniformity introduced by Atchity and Ruedenberg is discussed.