Is the Spin of the Black Hole in GX 339-4 Negative?

We have studied the accreting black hole binary GX 339-4 using two highly accurate broad-band X-ray data sets in very soft spectral states from simultaneous NICER and NuSTAR observations. Simultaneous fitting of both data sets with relativistic models of the disk, its Comptonization and reflection allows us to relatively accurately determine the black-hole mass and spin, and the distance and inclination. However, we find the measured values strongly depend on the used disk model. With the widely used thin-disk Kerr models kerrbb and kerrbb2 (which employ color corrections), we find relatively low masses and strongly negative spins. Then, the models utilizing detailed disk atmospheric spectra, bhspec and slimbh, predict moderately positive spins and high masses. When adding a warm corona above the disk (as proposed before for both AGNs and accreting binaries), we find the spin is weakly constrained, but consistent with zero. In all cases, the fitted inclination is low, $\approx$30-$34^\circ$. For the spin aligned with the binary orbit, the mass function for this binary implies large values of the mass, consistent only with those obtained with either slimbh or warm corona. We also test the disk models for an assumed set of mass, distance and inclination. We find that, e.g., kerrbb yields values of the spin parameter lower than bhspec or slimbh by $\sim$0.2-0.3. Our results confirm previously found strong disk-model dependencies of the measured black-hole spin, now for a low-mass X-ray binary.