  Autocatalytic reaction system with a small number of molecules is studied
numerically by stochastic particle simulations. A novel state due to
fluctuation and discreteness in molecular numbers is found, characterized as
extinction of molecule species alternately in the autocatalytic reaction loop.
Phase transition to this state with the change of the system size and flow is
studied, while a single-molecule switch of the molecule distributions is
reported. Relevance of the results to intracellular processes are briefly
discussed.
