The solubility data of epsilon-2,4,6,8,10,12-hexanitro-2,4,6,8,10,12-hexaazaiso-wurtzitane (?-HNIW) in five pure solvents (ethyl acetate, N,N-dimethylformamide, acetonitrile, cyclohexanone, N-methylpyrrolidone) at temperatures ranging from 283.15 to 333.15 K at atmospheric pressure were measured by the isothermal method; meanwhile, the solubility behavior of epsilon-HNIW in four binary solvents (dichloromethane, chloroform, tetrachloromethane, and 1,2-dichloroethane mixed with various mole fractions of ethyl acetate, respectively) was experimentally studied by the equilibrium method at temperatures ranging from 283.15 to 323.15 K. The experimental data in pure solvents were correlated by the vant Hoff equation and modified Apelblat equation. It was found that both of them could obtain good correlation results according to relative average deviations, root-mean-square deviations, and correlation coefficients (R-2). The solubilities of epsilon-HNIW in binary solvents were correlated with the modified Apelblat equation, the combined nearly ideal binary solvent/Redlich-Kister (CNIBS/R-K) equation, and the Jouyban-Acree model. For the binary solvents studied, all of the models provide a satisfactory correlation.