To perform a single Monte Carlo step, we first randomly select one of the considered reactions, each of them with equal probability. Next, we choose whether to carry out the reaction move in the forward or backward direction, also with equal probabilities. For a chosen direction, a trial reaction move is carried out by removing particles, by changing the particle identities, inserting new particles at random positions or deleting old ones, as prescribed by the stoichiometry of the corresponding reaction. The proposed new trial state (n) is accepted with the probability Ultimately, the input of the method is the set of equilibrium constants of all reactions, Γ. The output is the composition of the system, determined by the average concentration of the exchangeable species in the simulation box. One can perceive this method as a combination of Monte Carlo in reaction ensemble, and semi-grandcanonical ensemble, where we fix the linear combinations of chemical potentials of species because the chemical potentials of individual ionic species are ill-defined. 105