The breakup of a spherical droplet in a decaying homogeneous isotropic turbulence is studied by solving the Cahn-Hilliard-Navier-Stokes equations, using the discrete unified gas kinetic scheme combined with the free-energy-based phase-field model. We focus on the combined effects of turbulence and surface tension on the breakup process by assuming that the two fluid phases have the same density and same viscosity. The key physical parameters of the system include the volume fraction (ϕ = 6.54%), the initial Weber number (We = 21.7), and the initial Taylor microscale Reynolds number (Re λ = 58). Due to the turbulence decay, the Weber number decreases monotonically in time to a value of less than 0.01, providing a great opportunity to study the competing effects of turbulent kinetic energy and interfacial free energy on the dynamics of the two-phase system. Three distinct stages of droplet evolution are identified, namely, the deformation stage when the initially spherical droplet evolves into an irregular geometric shape with complex structures, the breakup stage when many daughter droplets are formed, and the restoration stage when the droplets relax towards spherical shape. These three stages are analyzed systematically from several perspectives: (1) a geometric perspective concerning the maximum equivalent diameter, the total number of droplets, total interface area, and probability distribution of droplet diameters, (2) a dynamic perspective concerning the evolution of local velocity and vorticity at the fluid-fluid interface, (3) a global perspective concerning the evolution of average kinetic energy / dissipation rate and their Fourier spectra, (4) spherical harmonics based energetics concerning simultaneous transfer of kinetic energy across different length scales and different radii relative to initial droplet center, and (5) the time evolution of global kinetic energy and free energy of the system. It is found that the ending time of the breakup stage can be estimated by the Hinze criterion. The kinetic energy of the twophase flow during the breakup stage is found to have a power-law decay with an exponent −1.76, compared to the exponent (−1.65) for the single-phase flow during the same time period, mainly due to the enhanced viscous dissipation generated by the daughter droplets.Energy spectra of the two-phase flow show power-law decay, with a slope between −4 and −3, at high wave numbers, both in the usual Fourier spectral space and in the spherical harmonics space.