QGSjet II and EPOS hadronic interaction models: comparison with the Yakutsk EAS array data

Various hadronic interaction models were used in extensive air showers simulations. This resulted in ambiguous estimation of primary energy, cosmic ray flux intensity, mass composition, etc. Several revisions of models have been made recently; for example, third major version of QGSjet II (QGSjet II-03) model was released, new models based on actual accelerator data appeared (EPOS). Employment of newer models always leads to new comprehension of experimental results. Nevertheless, in this case there still is some ambiguity. It is a matter of how correct does the model extrapolate characteristics of primary particle interaction with nuclei of the air from high energies to ultra-high.