Power GCD Matrices for a UFD