Skip to content

Understanding combustion of H\(_2\)/O\(_2\) gases inside nanobubbles generated by water electrolysis using reactive molecular dynamics simulations

Summary

ReaxFF RMD in LAMMPS (USER-REAXC) studies spontaneous H\(_2\)/O\(_2\) chemistry inside nm cubic cells (30–120 nm) at ~300 K initial conditions and 2–80 atm equivalent pressures, with surface-assisted H/O radical scenarios. H\(_2\)O\(_2\) dominates stable oxygenated products vs H\(_2\)O under these low-T, high-P nanoconfinement conditions; net temperature rise remains small (~68 K) because heat loss to cold walls dominates. The physical picture is electrolytic nanobubble chemistry: H\(_2\) and O\(_2\) can coexist metastably in water at ambient T, but radical-seeded pathways unlock low-temperature oxidation chemistry more typical of high-T flames—here examined in closed nanoscale reactors rather than macroscopic burners.

Methods

  • Engine / code: LAMMPS with the USER-REAXC reactive molecular dynamics package implements the ReaxFF integration.
  • Potential: ReaxFF parametrization by Agrawalla for H/O chemistry emphasizing HO\(_2\)/H\(_2\)O\(_2\) pathways at high P, low T; bond-order cutoffs define species (Table 1 in paper).
  • Systems: Stoichiometric H\(_2\):O\(_2\) 2:1 mixtures (124200/62100 molecules baseline count—see article for scaling with box); cubic boxes set by target P.
  • Initial equilibration: NVT at 300 K for 0.5 ns; timestep 0.1 fs; Nosé–Hoover (100 fs relaxation).
  • Radicals: H/O added at <3% box edge thickness—either initial burst (1–6% of moles) or continuous each 0.01 ns; reference cases in Table 2.
  • Combustion integration: Hybrid NVE interior + NVT 300 K shell (<3% thickness) as isothermal boundary; dynamic atom groups; 0.1 fs timestep; thermostat 100 fs on boundary atoms.

Pressure / barostat / bias / sampling. Initial gas pressures (2–80 atm equivalents) set the cubic box side lengths; N/A — NPT barostat — the reactive stage uses fixed-volume NVE/NVT hybrid control rather than fluctuating-box NPT. N/A — applied electric field — not part of the reported protocol. N/A — umbrella sampling, metadynamics, or replica exchange — not used.

Findings

Outcomes and mechanisms. Under ~300 K initial conditions and multi-atm nanoconfinement, H\(_2\)O\(_2\) can dominate stable oxygenated products over H\(_2\)O, contrasting macroscopic flame oxidation pathways. H\(_2\) and O\(_2\) consumption accelerates with radical seeding and with higher initial pressure.

Comparisons. The authors frame the results relative to classical high-temperature combustion literature, emphasizing how low-temperature, high-pressure confinement shifts product speciation.

Sensitivity. Radical delivery (initial burst vs continuous edge injection per Table 2) changes reaction progress vs time; net temperature rise stays modest (~68 K in the summarized case) because heat loss to the cold boundary shell dominates—nanoscale “combustion” is loss-limited despite oxidation chemistry.

Limitations and PDF grounding. ReaxFF scope, idealized cubic cells, and phenomenological radical injection limit direct mapping to electrolytic nanobubbles; electrochemical polarization and explicit electrolyte ions are outside the model. Quantitative statements should track the J. Phys. Chem. A PDF (pdf_path).

Limitations

ReaxFF scope; idealized cubic nanobubble; radical injection is phenomenological vs electrolysis reality. Electrochemical overpotential, double-layer fields, and bubble nucleation kinetics at electrodes are not represented explicitly, so quantitative mapping to sonochemical or electrolytic experiments requires multiphysics extensions beyond the gas-phase ReaxFF cell studied here. Water dissociation and ion product speciation may differ if explicit electrolyte ions are included.

Relevance to group

Low-temperature high-pressure reactive MD benchmark for H/O chemistry distinct from shock or flame MD in other corpus entries.

Citations and evidence anchors

  • DOI: 10.1021/acs.jpca.8b01798.

Reader notes (navigation)

Treat this paper as a low-temperature, high-pressure complement to shock and flame ReaxFF studies elsewhere in the corpus.