In this study, we investigate the effect of structural disorder on structural, elastic, and magnetic properties of Fe \(_{100-x}\) Ga \(_x\) ( \(x = 3.7, 5.5, 9.3\) at.%) alloys in the A2 phase by means of the density functional theory and Monte Carlo simulations. For each composition of Ga, ten configurations with a random arrangement of non-magnetic atoms in the A2 structure were generated. Calculations performed using the density functional theory have shown that, depending on the localization of atoms in the calculation cell, the calculated equilibrium volume of the unit cell per atom can either increase or decrease, which affects the magnitude of the magnetic moment per Fe atom. The elastic moduli, obtained in this work, has the same trend as observed in experiment. The magnetic exchange interaction between nearest neighbor Fe pair is ferromagnetic and about 22–25 meV. By using the long-range exchange interactions between magnetic atoms in the classical Heisenberg Hamiltonian, we have simulated the temperature dependencies of magnetization for Fe \(_{100-x}\) Ga \(_x\) alloys at zero magnetic field and estimated Curie temperatures. By using the long-range exchange interactions between magnetic atoms in the classical Heisenberg Hamiltonian, we have simulated the temperature dependencies of magnetization for Fe \(_{100-x}\) Ga \(_x\) alloys at zero magnetic field. It was shown that the Curie temperature is directly related to the values of the magnetic exchange parameter and demonstrates an increasing dependence. The following results were obtained for the studied compounds: 1138 K for Fe \(_{96.3}\) Ga \(_{3.7}\) , 1185 K for Fe \(_{94.5}\) Ga \(_{5.5}\) , and 1267 K for Fe \(_{90.7}\) Ga \(_{9.3}\) .