0 1 2 / c 0.8 0.85 0.9 0.95 1 1.05 T/T c =7, NVT−MC =9, NVT−MC =10, NVT−MC =12, NVT−MC =15, NVT−MC =8, GC−TMMC =9, GC−TMMC =10, GC−TMMC FIG. 4: Reduced vapor-liquid coexistence curves of SR-HCAY fluid. NVT-MC data from this work and Ref. [33]. GC-TMMC data from Ref. [36]. short range HCAY fluid form a single master curve within a very good degree of accuracy, as it was reported by one of us in a previous work for longer range [35]. As the range is decreased, a slight trend towards a widening of the coexistence curve can be detected, but this effect remains small for the ranges considered here. Also, the coexistence curves reported by Singh [36] using GC-TMMC are on the master curve. The same behavior has been found for the Mie(n,m) potential [58–60]. These results confirm that, even for the short range interactions considered here, with inverse-range parameter as large as κ = 15, the vapor-liquid phase diagrams of HCAY fluid obey the law of corresponding states. This fact could be used to test theoretical approaches. V. CONCLUSIONS We have carried out a systematic study of the liquid-vapor phase diagrams for hard- core attractive Yukava potential with κ = 9 , 10 , 12 and 15, using NVT-MC simulation. Our coexistence curves for κ = 9 and 10 were compared to those reported by Singh using GC-TMMC. An excellent agreement was found between them. Besides, we have shown that reduced coexistence curves of HCAY fluid from κ = 7 to κ = 15 follow the law of corresponding states, as it was shown in a previous work [35] for longer interaction ranges. Once again, we have confirmed the linear relationship between the critical density and inverse 15