Spectral variability in NGC 1042 ULX1

We report X-ray spectral variability in an ultraluminous X-ray source NGC 1042 ULX1, using archival XMM-Newton and recent NuSTAR observations. In long-term evolution, the source has shown a trend of variation in spectral hardness. The variability in different XMM-Newton observations is prominent above ~1 keV. Cool thermal disc component with a characteristic temperature of ~0.2 keV manifests that the spectral state of NGC 1042 ULX1 in all epochs is similar to that of the ultraluminous state sources. An apparent anti-correlation between luminosity and power-law index demonstrates that the source becomes spectrally harder when it is in a brighter state. That is conceivably related to variation in accretion rate, strength of Comptonization, wind/outflow in the system or a manifestation of varying disc occultation. Typical hard ultraluminous-type spectra indicate that NGC 1042 ULX1 is a low inclination system in general. Spectral properties suggest that, like many other ULXs which show spectral curvature around ~6-10 keV, NGC 1042 ULX1 could be another stellar-mass super-Eddington accretor.