The problem of thermal ignition in a homogeneous gas is revisited from a molecular dynamics perspective. A two-dimensional model is adopted, which assumes reactive disks of type A and B in a fixed domain that react to form type C products if an activation threshold for impact is surpassed. Such a reaction liberates kinetic energy to the product particles, representative of the heat release. The results for the ignition delay are compared with those obtained from the continuum description with the reaction rate evaluated from kinetic theory assuming local thermodynamic equilibrium and Maxwell-Boltzmann statistics, in order to assess the role played by molecular fluctuations. Ignition times obtained using molecular dynamics are ensemble averaged over 100 simulations to address the statistics of the ignition event. Results show two regimes of non-equilibrium ignition whereby ignition occurs at different times as compared to that for homogeneous ignition assuming local equilibrium. The first regime is at low activation energies, where the ignition time is found to be higher than that expected from theory for all values of heat release. The lower reaction rate is shown to occur due to a departure from local equilibrium for the different species, in agreement with predictions from Prigogine and Xhrouet. In this low activation energy regime, the ignition times from molecular dynamics are also found to be independent of domain size and there is little variance between different realizations under similar conditions, which suggests that the ignition is spa-Email addresses: nsirmas@uottawa.ca (N. Sirmas), matei@uottawa.ca (M. I. Radulescu)