def METHOD_NAME(self, order=1):
    """Calculate phonon state moment of specified order."""
    if self._eigenvectors is None:
        self._get_moment(order)
    else:
        self._get_projected_moment(order)