  The dissolution of CO$_2$ in water followed by the subsequent hydrolysis
reactions is of great importance to the global carbon cycle and carbon capture
and storage. Despite enormous previous studies, the reaction pathways and
kinetics are still not fully understood at the atomistic scale. Here, we
constructed Markov state models based on $ab\ initio$ molecular dynamics
simulations to elucidate the reaction mechanisms and kinetics of dissolved
carbon in supercritical water both in the bulk and nanoconfined states. Unlike
many previous molecular simulations using enhanced sampling methods, our method
can automatically identify complex reaction coordinates and pathways with
multiple intermediates using unsupervised machine learning techniques instead
of $a\ priori$ human speculation. In the bulk solution, CO$_2$(aq) tends to
directly react with H$_2$O or OH$^-$ to generate HCO$_3^-$ and H$_2$CO$_3$(aq),
whereas under graphene nanoconfinement the dissolution of CO$_2$(aq) involves
the pyrocarbonate (C$_2$O$_5^{2-}$(aq)) ion as an intermediate state. While it
is known that pyrocarbonate does not exist in aqueous solutions, our study
suggests that the extreme hydrophobic confinement may greatly enhance the
stability of pyrocarbonate ions in water. We found that carbonation reactions
involve collective proton transfers along transient water wires, which exhibit
concerted behavior in the bulk solution but proceed stepwise under
nanoconfinement. Our study provides valuable insights into the reaction kinetic
network of aqueous carbon, yielding significant implications for the deep
carbon cycle and the sequestration of CO$_2$. Combining quantum molecular
dynamics with Markov state models shows great potential to elucidate complex
reaction pathways and kinetics.
