This paper analyzes the semiparametric estimation of multivariate long-range dependent processes. The class of spectral densities considered is motivated by and includes those of multivariate fractionally integrated processes. The paper establishes the consistency of the multivariate Gaussian semiparametric estimator (GSE), which has not been shown in other work, and the asymptotic normality of the GSE estimator. The proposed GSE estimator is shown to have a smaller limiting variance than the two-step GSE estimator studied by Lobato (1999). Gaussianity is not assumed in the asymptotic theory. Some simulations confirm the relevance of the asymptotic results in samples of the size used in practical work.