Abstract. In this paper, a mass (or normalization) and magnetization conservative and energydiminishing numerical method is presented for computing the ground state of spin-1 (or F = 1 spinor) Bose-Einstein condensates (BECs). We begin with the coupled Gross-Pitaevskii equations, and the ground state is defined as the minimizer of the energy functional under two constraints on the mass and magnetization. By constructing a continuous normalized gradient flow (CNGF) which is mass and magnetization conservative and energy-diminishing, the ground state can be computed as the steady state solution of the CNGF. The CNGF is then discretized by the Crank-Nicolson finite difference method with a proper way to deal with the nonlinear terms, and we prove that the discretization is mass and magnetization conservative and energy-diminishing in the discretized level. 1. Introduction. Since its realization in dilute bosonic atomic gases [2, 13, 9], the atomic Bose-Einstein condensate (BEC) has been produced and studied extensively in the laboratory [28,29,16] and has provided a successful testing ground of theoretical studies of quantum many-body systems [28,29]. In earlier BEC experiments, atoms were spatially confined with magnetic traps, which essentially freeze the atomic internal degrees of freedom [2,13,9]. Most studies were thus focused on scalar models, i.e., single-component quantum degenerate gases [12]. One of the most important recent developments in BECs was the study of spin-1 condensates (of atoms with hyperfine quantum number F = 1) [17,27,34,10,31], and they were realized in experiments recently using both 23 Na and 87 Rb [24,35]. In fact, the emergence of the spin-1 BEC [19,20,24] has created opportunities for understanding degenerate gases with internal degrees of freedom [21,22,17,18,14,25,26,32,37].At temperature T much smaller than the critical condensate temperature T c [23], a spin-1 BEC is well described by the three-component wave function Ψ = (ψ 1 (x, t), ψ 0 (x, t), ψ −1 (x, t)) T whose evolution is governed by the coupled Gross-Pitaevskii equations (GPEs) [23,17,18,38,36]: