Biohydrogen Production via Genetically Optimized Dual-Stage Fermentation with Directed CO₂ Capture and Biological Recycling: A Theoretical Framework

1. Introduction

The transition to clean and sustainable energy sources remains one of the most urgent scientific challenges of the 21st century. Hydrogen has emerged as a promising candidate due to its high energy content and zero-emission profile at the point of use. However, current industrial hydrogen production methods—most notably steam methane reforming (SMR)—are carbon-intensive, heavily reliant on fossil fuels, and contribute substantially to greenhouse gas emissions. While water electrolysis offers a cleaner alternative, it is hindered by high energy requirements and dependence on renewable electricity infrastructure, limiting its scalability—particularly in regions with limited access to renewable energy systems.

Biological hydrogen production through microbial fermentation has gained increasing attention as a low-energy, renewable pathway. Specific strains such as Clostridium butyricum and Rhodobacter sphaeroides can convert organic substrates into hydrogen gas under anaerobic conditions. Dark fermentation using C. butyricum can rapidly degrade sugars like glucose into hydrogen and volatile fatty acids, while photofermentation by R. sphaeroides enables the further conversion of these acids—particularly butyric and acetic acid—into additional hydrogen using light energy. In addition to carbon-based substrates, water molecules serve as critical electron donors during photofermentation, providing the oxygen atoms released as carbon dioxide and the protons required for sustained hydrogen generation.

Despite their promise, single-stage microbial systems are typically limited by incomplete substrate utilization, accumulation of metabolic byproducts, and suboptimal hydrogen yields. This paper proposes a dual-stage hydrogen production model that integrates genetically optimized strains of C. butyricum and R. sphaeroides in a vertically structured bioreactor. By coupling dark fermentation with downstream photofermentation, the system aims to maximize hydrogen recovery from both glucose and its metabolic derivatives. Genetic modifications are introduced to enhance hydrogenase activity, reduce competing metabolic fluxes, and improve light harvesting under low-intensity conditions.

In addition to yield optimization, the system incorporates a membrane-based gas separation unit to extract and isolate carbon dioxide produced during fermentation. Rather than storing or releasing this CO₂, the model proposes redirecting it to external biological systems—such as algae bioreactors, sugarcane plantations, or bamboo forests—for natural carbon sequestration and glucose regeneration. This combined approach aims to support the development of a decentralized, low-emission hydrogen production platform that aligns with global sustainability goals.

This study not only proposes a theoretical biohydrogen system but also introduces a time-resolved kinetic modelling framework to simulate microbial reaction rates and predict real-time behaviour. Additionally, it incorporates a carbon sequestration strategy via high-efficiency C₄ crop integration, allowing passive redirection of emitted CO₂ into a sugarcane-based sink. These elements collectively enhance the model's viability for real-world deployment.

2. Methodology and Theoretical Modelling Framework

This study employs a layered theoretical modelling approach to simulate microbial hydrogen production via a genetically optimized dual-stage fermentation system. As the model is conceptual and not experimentally validated, the methodology emphasizes stoichiometric yield estimation, metabolic pathway design, kinetic approximation, and energy budgeting under idealized constraints. The following components define the modelling framework:

2.1. Stoichiometric Modelling

Hydrogen yield predictions are derived from balanced metabolic equations for Clostridium butyricum (dark fermentation) and Rhodobacter sphaeroides (photofermentation). Mass- and mole-based hydrogen outputs are calculated using standard molar masses under complete substrate conversion assumptions.

2.2. Genetic Engineering Assumptions

Proposed genetic modifications are modelled as fully functional and chromosomally integrated. Edits such as overexpression of hydrogenases and metabolic enzymes or gene knockouts are based on documented literature. The model assumes no expression instability, mutation, or metabolic burden, representing best-case functional performance.

2.3. Bioreactor Design Assumptions

The reactor is conceptualized as a vertically integrated, gravity-assisted dual-chamber system with scaffold-supported microbial biofilms. Mass transfer of VFAs between stages is passive. Each chamber operates under fixed environmental conditions (anaerobic, temperature-stabilized), with light access provided for the photofermentative stage. Uniform biofilm structure, metabolite diffusion, and chamber isolation are assumed.

2.4. Kinetic Approximations

Microbial activity is modelled using Monod kinetics for glucose consumption and light-saturation kinetics for photofermentation. Empirical constants are adapted from published literature. These models provide time-based performance insights but are not calibrated to experimental data.

Important Note on Hydrogen Feedback Inhibition:

While this model simulates microbial kinetics under ideal extraction conditions, it does not account for enzymatic inhibition due to hydrogen accumulation in the reactor headspace. Literature reports show that hydrogenase and nitrogenase activity in C. butyricum and R. sphaeroides respectively are suppressed at H₂ partial pressures >0.6–0.8 atm. This inhibition can significantly reduce gas output in poorly ventilated systems. Future kinetic models should incorporate dynamic headspace gas modelling and feedback inhibition thresholds for accurate time-course prediction.

2.5. Energy and Cost Calculations

Energy requirements for illumination, gas extraction, and basic control systems are estimated using power consumption values from typical lab-scale equipment. Hydrogen output energy is calculated using the combustion enthalpy of H₂ (286 kJ/mol). Cost analysis is based on U.S. bulk pricing for glucose, amortized component use, and power rates, reflecting theoretical minimums.

(All pricing references are based on pre-2020 U.S. bulk rates, before tariff-induced market distortions introduced under Trump’s presidency. This ensures pricing reflects commodity baselines unaffected by policy-driven inflation.)

2.6. Gas Separation and CO₂ Routing

A vacuum-assisted membrane system is modelled with perfect hydrogen/CO₂ separation efficiency, no fouling, and lossless operation. Separated CO₂ is routed to sugarcane-based sinks via passive or root-zone injection. Biological sinks are assumed to operate at full fixation capacity with no lag or underperformance.

Author’s Note on Electrochemical Capture Alternatives:

During the modelling phase, I reviewed electrochemical CO₂ capture approaches such as the fluoflavine-based FFDS redox system [Wei et al., 2025], which offers oxygen-tolerant operation and reversible capture via pH-swing chemistry. While conceptually interesting, I found the overall energy-to-capture ratio and CO₂ managing cycle to be suboptimal for integration into this system. In particular, the reliance on high energy input without a built-in utilization path renders it inefficient when compared to biological fixation strategies.

Therefore, this model intentionally favours photosynthetically active CO₂ sinks—namely sugarcane plantations and algae bioreactors—which not only sequester CO₂ passively but regenerate glucose substrates for further biohydrogen production. Unless future work demonstrates a tightly coupled valorisation pathway for electrochemically captured CO₂, I consider this class of sorbents unsuitable for the present system.

2.7. Scope Limitations

This framework assumes ideal substrate purity, enzyme functionality, membrane integrity, and microbial behaviour. Stochastic variability, environmental drift, and hardware degradation are not included. These simplifications serve to define a theoretical upper-bound system performance.

This section consolidates the modelling strategies underlying the theoretical framework, ensuring clarity and reproducibility for future studies. While certain sections of the paper include personal reflections to contextualize design constraints, the scientific modelling remains rigorous and idealized for conceptual benchmarking.

2.8. Assumptions Tiers and Realism Factors

To clearly delineate the scope of this theoretical model, all system assumptions are categorized by domain and impact severity. While idealized conditions are necessary to establish a performance ceiling, this section outlines the tiers of assumptions applied and proposes adjustments for more realistic modelling.

Assumptions are grouped into four tiers—structural, biological, engineering, and environmental—each with distinct implications for system behaviour and model reliability.

Tier	Assumption Type	Examples	Treatment
Tier 1 – Structural	Gravity-fed VFA transfer, scaffold-stable biofilms, anaerobic integrity	Accepted for conceptual clarity	Supported by bench-scale analogues, but flow obstruction and gradient formation likely in real systems
Tier 2 – Biological	100% enzyme expression, no mutation, full substrate conversion, no biomass loss	Idealized; adjustable using Biological Realism Factor (BRF)	Literature supports 10–25% expression drop over time; real-world BRF = 0.75–0.90
Tier 3 – Engineering	Perfect membrane separation (100% H₂/CO₂ split), zero fouling, no vacuum loss	Accepted for theoretical benchmarking	Real systems show 85–95% separation efficiency with membrane aging and backpressure; fouling mitigation required
Tier 4 – Environmental	Continuous CO₂ fixation by sugarcane, zero plant stress, no uptake lag	Assumed for CO₂ neutrality model	Real-world uptake is climate- and growth-phase-dependent; sequestration efficiency highly variable

To account for biological and engineering variability, the Biological Realism Factor (BRF) is introduced as a global modifier. This factor adjusts theoretical hydrogen yields and gas separation performance downward based on expected losses in enzyme expression, cofactor depletion, metabolic burden, and system wear.

Typical BRF values based on published microbial expression data range from 0.75 to 0.90, depending on strain, expression system, and cultivation duration. All future yield, energy, and cost calculations in this paper can be scaled using this factor to reflect more realistic operating envelopes.

All core values can be scaled by the BRF for realism-adjusted estimates. Tables in Section 9 and Section 14 reflect uncorrected theoretical ceilings; adjusted versions appear alongside them for realistic comparison.

The tiered structure also supports future scenario modelling—allowing researchers to evaluate system resilience under degraded conditions, such as partial membrane failure, incomplete VFA conversion, or temporary CO₂ sink inoperability.

2.8.1. Realism Adjustment Factors and Modelling Scope

To quantitatively bridge the gap between idealized assumptions and real-world system behaviour, this model introduces a set of realism adjustment factors. Each factor targets a specific constraint domain—biological, genetic, engineering, or metabolic—allowing theoretical outputs to be scaled toward plausible operating ranges. These factors serve as modular modifiers for sensitivity analysis, scenario modelling, and future experimental calibration. This table summarizes the core factors used throughout this study.

Factor	Full Name	Purpose	Typical Range	Applied In
BRF	Biological Realism Factor	Adjusts for enzyme expression loss, mutation, and metabolic burden	0.75 – 0.90	H₂ yield, cofactor availability
GSCF	Gene Stack Correction Factor	Reduces yield due to combinatorial gene expression overload	0.6 – 0.9	Total H₂ yield per cycle
MFCF	Membrane Flux Correction Factor	Compensates for real-world limitations in gas extraction throughput	0.7 – 0.85	Effective H₂ recovery vs. inhibition
BDF	Biomass Diversion Factor	Represents substrate carbon diverted to biomass instead of gas	0 – 30%	Net CO₂ & H₂ output, system efficiency

2.9. Core Parameter Reference Table

This table consolidates key assumptions, constants, and system parameters used throughout the model for ease of cross-reference.

Parameter	Value	Units	Context/Notes
Molar mass of glucose	180.16	g/mol	Primary substrate
Molar mass of water	18.015	g/mol	Photofermentation donor
Molar mass of H₂	2.02	g/mol	For mass yield
H₂ energy content	286	kJ/mol	Based on combustion enthalpy
H₂ yield per cycle	12	mol	Theoretical max (6 from glucose, 6 from water)
CO₂ output per cycle	6	mol	2 from dark, 4 from photofermentation
Membrane separation efficiency	92	%	Real-world polymer-MOF hybrid estimate
LED light intensity	30	µmol/m²/s	Minimum needed for nitrogenase activation
Nitrogenase ATP cost	20	mol ATP/mol H₂	Average demand estimate
Vacuum pump energy	216	kJ/cycle	30 W for 2 hours
Total input energy	960	kJ/cycle	Including LED + mixing + control
BRF range	0.75–0.90	-	Adjusts for biological inefficiencies
BDF range	0–30	%	Represents substrate diverted to biomass
Areal light requirement	≥0.02	m² per L	Minimum surface exposure for photofermentation
Photofermentative H₂ yield	10	Mol	From 1 mol butyric acid + 6 mol water
Effective efficiency (realistic)	40–60	%	Total system, biologically corrected
H₂ volume per cycle	~89.7	-	At STP, for 12 mol H₂

Consistency Notes on STP Assumptions

All gas volumes reported in this paper are assumed to be measured under Standard Temperature and Pressure (STP) conditions (1 atm, 273.15 K), unless otherwise specified. This standardization ensures comparability across reported hydrogen yields, especially in volumetric metrics.

Terminology Reference Table

To ensure clarity throughout this paper, the following standardized terminology is adopted:

Term	Definition	Usage Context
Theoretical Ceiling	The maximum yield or performance achievable under idealized biological and engineering conditions. Assumes perfect enzyme expression, complete substrate conversion, and zero system losses.	Used for defining best-case outputs in stoichiometric and energy models.
Idealized Model	A conceptual model that excludes stochastic effects, degradation, or regulatory limitations. Used to benchmark the system's theoretical potential.	Applied in early-stage modelling (Section 2, Section 3, Section 4 and Section 5).
Realistic Scenario	A model adjusted using Biological Realism Factor (BRF) and Biomass Diversion Factor (BDF) to account for expression loss, cofactor limitations, and carbon retention.	Applied in Section 8.3–9 for grounded yield and efficiency estimates.
Yield by Mass	Percentage of hydrogen mass relative to input substrate mass. Contextualized against either glucose-only, glucose + water, or butyrate + water.	Standard in biochemical literature (Sections 3, 8).
Energy Yield	Total energy output from H₂ combustion (in kJ or MJ), based on molar hydrogen production.	Used for benchmarking system viability and efficiency (Section 9).
Volumetric Productivity	Hydrogen output per litre of reactor volume per day (L/L/day).	Used for scaling comparisons in fermentation systems (Section 10).
Areal Productivity	Hydrogen output per square meter of light-exposed surface area per day (L/m²/day).	Used to evaluate light-limited photofermentation stages (Section 5.2.1).

3. Biochemical Pathways and Theoretical Yield Modelling

Hydrogen production through microbial metabolism is driven by a series of enzymatically regulated biochemical pathways, with stoichiometric constraints that define the theoretical maximum yields achievable from specific substrates. This section models the metabolic pathways of two key microorganisms—Clostridium butyricum and Rhodobacter sphaeroides—within the context of a dual-stage fermentation system. Yield calculations are based on molar balances and mass conversions to assess the energy potential and gas output of the proposed design.

3.1. Dark Fermentation Pathway (Clostridium butyricum)

Dark fermentation is an anaerobic process wherein glucose is catabolized to yield hydrogen gas, carbon dioxide, and volatile fatty acids (VFAs), primarily butyric and acetic acids. The simplified stoichiometric reaction is as follows:

C₆H₁₂O₆ → C₄H₈O₂ + 2CO₂ + 2H₂

3.1.1. Theoretical Hydrogen Yield Calculation:

• Molar mass of glucose (C₆H₁₂O₆): 180.16 g/mol
• Molar mass of hydrogen (H₂): 2.02 g/mol
• 2 mol H₂ = 4.04 g hydrogen per mole glucose
• Yield by mass = $\left({\displaystyle \frac{4.04}{180.16}}\right)\times 100\approx 2.24\mathit{\%}$

• Key Enzymes:

• Glycolysis: Glucose → Pyruvate
• Pyruvate:ferredoxin oxidoreductase: Pyruvate → Acetyl-CoA + CO₂ + reduced ferredoxin
• [FeFe]-Hydrogenase: Reduced ferredoxin → H₂

3.2. Photofermentation Pathway (Rhodobacter sphaeroides)

In the second stage, R. sphaeroides photoferments VFAs under anaerobic, light-activated conditions to produce additional hydrogen. A representative reaction using butyric acid is:

C₄H₈O₂ + 6H₂O + Light → 10H₂ + 4CO₂

The representative stoichiometric reaction for photofermentation used in this model—C₄H₈O₂ + 6H₂O + Light → 10H₂ + 4CO₂—reflects the complete oxidation of butyric acid under idealized conditions. This formulation aligns with theoretical yields reported in prior photofermentation studies, where yields of up to 10 mol H₂ per mol of butyrate have been proposed based on full hydrogen atom recovery from both the substrate and water-derived protons. Notably, studies by Ghosh and Hallenbeck (2010) and Nath & Das (2004) support theoretical ceilings in the 8–10 mol H₂ range when nitrogenase is fully activated and substrate inhibition is minimal. These values serve as a theoretical upper limit and are not typically achieved in practical systems, where reported yields range from 5.6 to 7.8 mol H₂/mol butyrate due to cofactor limitations, metabolic branching, and photon attenuation.

3.2.1. Theoretical Hydrogen Yield Calculation:

• Molar mass of butyric acid (C₄H₈O₂): 88.11 g/mol
• 10 mol H₂ = 20.2 g hydrogen per mole butyric acid
• Total mass input = 88.11 + 108.09 = 196.2 g
• Yield percentage = $\left({\displaystyle \frac{20.2}{196.2}}\right)\times 100\approx 10.29\mathit{\%}$ by mass

• Key Components:

• Light-harvesting complexes (LH1, LH2)
• Nitrogenase enzyme (light activated, ATP-dependant)
• Electron donors: reduced Ferredoxin and NAD(P)H

Side Note: While the yield of 10 mol H₂ per mol butyric acid represents the full theoretical recovery of all hydrogen atoms, experimental photofermentation systems commonly achieve yields between 5.6 and 7.8 mol H₂ per mol substrate due to energy losses, biomass generation, and metabolic branch pathways [2,3]. This model assumes full oxidation as a theoretical ceiling and excludes these real-world constraints.

• Real-World Photofermentation Ceiling Adjustment

While the 10 mol H₂ per mol butyrate value used in this model aligns with theoretical maximums reported in prior photofermentation studies [Ghosh & Hallenbeck, 2010; Nath & Das, 2004], real-world yields are consistently lower due to kinetic constraints, light attenuation, cofactor scarcity, and regulatory bottlenecks. Most experimental setups report practical ranges of 5.6–7.8 mol H₂/mol butyrate under optimized but realistic conditions.

To reflect this discrepancy, a corrected yield range is proposed below. These values can be used to recalibrate energy output, system efficiency, and CO₂ generation when evaluating real-world deployment scenarios:

Metric	Ideal	Realistic (7 mol)	Realistic (5.6 mol)
H₂ per mol butyrate	10 mol	7 mol	5.6 mol
Total H₂ per cycle	12 mol	9 mol	7.6 mol
Energy Output	3432 kJ	2574 kJ	2174 kJ
Net Energy Gain	2472 kJ	1614 kJ	1214 kJ
Efficiency (EROIE)	3.58	2.68	2.26

This yield gradient reinforces the importance of integrating Biological Realism Factor (BRF) and cofactor correction layers into all performance metrics, rather than relying solely on stoichiometric ceilings. Future reactor designs should benchmark against the 7 mol case unless continuous light optimization and genetic stabilization are proven experimentally.

3.3. Overall Hydrogen Yield from One Cycle

Assuming complete substrate conversion in both stages:

From 1 mol of glucose and 6 mol of water, the system yields a total of 12 mol of hydrogen—recovering 6 mol from glucose-derived metabolites and 6 mol from water-derived protons via photofermentation. This distribution reflects full hydrogen recovery from both the substrate and external electron donor sources. This corresponds to mass-based yields of 13.45% relative to glucose alone, and 8.41% when water input is included. Presenting both molar and mass-based values provides clarity on yield interpretation across biochemical and engineering contexts.

Side note: This yield assumes full constructive collaboration across all genetic edits. Applying a GSCF of 0.75–0.9 reduces this to 9–10.8 mol under more realistic gene stack limits.

H2

sub>2 Volume Clarification:

All hydrogen volume calculations throughout this model assume Standard Temperature and Pressure (STP) conditions—specifically, 1 atm pressure and 273.15 K (0°C). Under these parameters, one mole of an ideal gas occupies 22.414 litres. Thus, the 12 mol H₂ yield per cycle corresponds to approximately 268.97 litres of gas at STP. This value scales linearly with reactor operation and can be recalculated for non-STP conditions using the ideal gas law (PV = nRT). All real-world volume estimates should incorporate system temperature and backpressure to avoid yield inflation.

H2 Source	Origin Compound	Hydrogen Atoms Contributed	Mol H2 Equivalent	Note
Glucose	Direct Substrate	12	6	Via dark fermentation (partial) and VFA-based photofermentation
Water	Electron Donor in Photofermentation	12	6	Split via light-driven nitrogenase pathway
Total		24	12	Mass balance validated

For yield definitions (mass yield, energy yield, productivity), see the table in Section 8.0.2.

3.3.1. Hydrogen Atom Origin Clarification

This breakdown aligns with mass balance and theoretical atom conservation. This ensures mass conservation is maintained and avoids misinterpretation of hydrogen recovery calculations.

Reaction Stage	Input Molecule(s)	Process	Hydrogen Atoms Contributed	Mol H2 Generated
Dark Fermentation	1 mol Glucose	Glycolysis → VFAs + H₂	4 atoms	2 mol
Photofermentation	1 mol Butyric Acid (from glucose) + 6 mol	Light-driven nitrogenase activity	20 atoms	10 mol
Total			24 atoms	12 mol

3.4. CO₂ Generation Summary

• Dark fermentation: 2 mol CO₂
• Photofermentation: 4 mol CO₂
• Total = 6 mol CO₂ per mol glucose
• Molar mass of CO₂ = 44.01 g/mol
• Total CO₂ mass = 264.06 g per cycle
• CO₂-to-glucose mass ratio = 146.56%
• Clarification on Carbon Flow and CO₂ Mass Output:

The CO₂-to-glucose mass ratio of 146.56% reflects not an increase in carbon atoms, but rather the combined effects of complete substrate oxidation and the incorporation of oxygen atoms from water during photofermentation. Under ideal conditions, all six carbon atoms from glucose are oxidized—two during dark fermentation and four during photofermentation—yielding 6 mol of CO₂.

However, when biomass diversion occurs (modelled via the Biomass Diversion Factor, BDF), a fraction of glucose-derived carbon is retained in microbial biomass instead of being oxidized. Since CO₂ production is tightly coupled to substrate oxidation—particularly in photofermentation—this diversion results in a direct reduction in CO₂ output and, by extension, a proportional drop in H₂ yield.

Therefore, this 146.56% ratio applies only under idealized, full-oxidation scenarios. In realistic systems, both CO₂ evolution and hydrogen production will decline according to the degree of biomass carbon retention.

3.5. Process Assumption

Parameter	Value	Justification
Glucose purity	100%	Ideal condition
Reactor temperature	37°C	Optimal for both organisms
Anaerobic Conditions	Strict	Necessary for both fermentation stages; oxygen presence inhibits key enzymes
Light intensity	30 µmol/m2/s	Low, continuous light (~30 µmol/m²/s) ensures constant nitrogenase activation
Compete substrate conversion	Yes	Assumes full conversion of glucose to VFAs, and VFAs to H₂, under ideal conditions
Water availability	Excess	Water is available and maintained at 37°C, optimal for both strains
Genetic enhancements	Functional	Assumes overexpression is successful and stable, with no mutation or loss
No microbial cross-contamination	True	Each microbial population functions independently without inhibitory interactions

3.5.1. Biomass Diversion Factor (BDF):

In realistic systems, 5–30% of substrate carbon is diverted to microbial biomass, reducing theoretical gas yield. This model introduces a BDF (Biomass Diversion Factor) to represent the proportion of substrate unavailable for gas production due to cell growth and maintenance. BDF-adjusted yield scenarios are presented in Section 8.2.2.

The model assumes full conversion of VFAs to hydrogen, excluding the effects of accumulation, substrate inhibition, or partial VFA routing into biomass. Future kinetic models should incorporate substrate conversion efficiency factors (e.g., 80–90%) for more realistic projections. These genetic modifications are designed with a balance between yield optimization and metabolic stability in mind. While theoretical in nature, all proposed edits are based on prior microbial engineering literature and assumed to be achievable using contemporary gene editing platforms. Further modelling of pathway fluxes and gene expression regulation would be required in future experimental stages.

• Biomass Diversion and Carbon Retention Note:

In practical systems, microbial biomass formation diverts a portion of carbon and reducing equivalents away from gas production. Each glucose molecule contains six carbon atoms, all of which are fully oxidized to CO₂ under ideal conditions. However, when cells retain carbon for growth, fewer carbon atoms are available for CO₂ formation—and by extension, for the hydrogen-producing oxidation reactions that release them. This coupling between carbon flow and hydrogen yield is critical; even small diversions (e.g., 10–20% of carbon) can reduce final hydrogen output by a comparable margin. These effects are quantified in Section 8.2.2 using the Biomass Diversion Factor (BDF).

4. Genetic Modifications and System Burden Analysis

To maximize hydrogen output, specific genetic modifications are proposed for the microbial strains Clostridium butyricum and Rhodobacter sphaeroides. These edits aim to enhance hydrogenase activity, suppress competing pathways, expand substrate flexibility, and increase cofactor and ATP availability. All modifications are assumed to be integrated chromosomally via CRISPR-Cas9 or homologous recombination techniques, with stable expression under anaerobic or light-inducible promoters. Potential metabolic burden, redox imbalances, and long-term evolutionary responses are excluded from the scope of this model.

4.1. Genetic Enhancements in Clostridium butyricum

Gene	Modification	Purpose	Rationale
hydA	Overexpression	Increase [FeFe]-hydrogenase activity	Boosts conversion of reduced ferredoxin to H₂
fruA	Insertion	Enable fructose uptake	Expands substrate range beyond glucose
ptsF	Insertion	Support fructose transport via PTS system	Enhances sugar transport efficiency
pta	Insertion	Increase acetyl-CoA → acetate flux	Raises ATP output during fermentation
ackA	Insertion	Complements pta for ATP-linked acetate production	Further improves energy yield and redox balance

4.2. Genetic Enhancements in Rhodobacter sphaeroides

Gene	Modification	Purpose	Rationale
hydG	Knockout	Prevent hydrogenase assembly suppression	Avoids inhibition of hydrogenase maturation
hupL	Knockout	Eliminate hydrogen uptake hydrogenase	Prevents H₂ consumption by native uptake systems
bchP	Overexpression	Improve bacteriochlorophyll synthesis	Enhances light harvesting for photofermentation
crtB	Overexpression	Support carotenoid biosynthesis	Stabilizes light-harvesting complexes
ackA	Insertion	Boost ATP yield from acetate metabolism	Supplies energy for nitrogenase-dependent H₂ production

4.3. Genetic Integration Challenges and Model Assumptions

• Stack Compatibility and Synergy Limitations

While each proposed genetic edit in isolation has been documented to improve hydrogen yield or metabolic efficiency, the assumption of full functional compatibility across a multi-gene stack presents a significant risk of overestimation. In real-world systems, combinatorial gene stacking beyond 3–4 edits frequently result in metabolic crosstalk, resource exhaustion, and unexpected feedback inhibition.

For facultative anaerobes like Clostridium butyricum and Rhodobacter sphaeroides, regulatory networks often suppress overexpression to maintain redox balance and cellular homeostasis. Literature suggests that above 4–5 simultaneous edits, constructive collaboration tends to plateau—and in some cases, becomes negative due to ribosomal bottlenecking, cofactor depletion, or transcript interference.

Therefore, while this model includes a 6–8 gene modification set for theoretical benchmarking, it is explicitly treated as a ceiling scenario. Actual implementation may require modular optimization or sequential engineering to preserve system stability. Future modelling efforts should incorporate a constructive collaboration correction factor or simulate dynamic trade-offs using flux balance or gene regulatory network models to determine yield inflection points under escalating metabolic burden.

This model assumes these modifications are functionally expressed and do not trigger regulatory shutdowns or metabolic collapse. These assumptions are made to isolate the theoretical yield ceiling of optimized microbial strains. Future iterations of the model will need to incorporate dynamic flux balance analysis (dFBA), metabolic burden modelling, or resource allocation frameworks to refine viability assessments for these engineered systems.

Similar enhancements in Clostridium spp. and Rhodobacter spp. have achieved 1.8x–3.5x increases in H₂ yield in batch systems under lab conditions [e.g., Hallenbeck et al., 2009; McKinlay & Harwood, 2010; Thaiwong et al., 2015]. This model builds upon those pathways with additional combinatorial edits. While constructive collaboration between these edits has not been empirically validated as a unified stack, this model assumes compatibility under optimized cellular conditions.

Based on literature, multi-gene overexpression systems in facultative anaerobes often experience 10–25% loss in target expression over extended growth due to mutation, regulatory repression, or resource depletion [ref Hallenbeck, McKinlay]. While the current model assumes full expression, a Biological Realism Factor (BRF) of 0.75–0.9 can be used to adjust theoretical rates in future kinetic simulations.

These genetic assumptions define the theoretical design envelope for engineered biohydrogen systems. While this paper explores the idealized ceiling, future work should define a performance corridor bounded by expected losses in expression fidelity, cofactor availability, and adaptive mutagenesis. Such a corridor would allow scenario-based system planning those accounts for both synthetic potential and biological resistance.

• Engineering Trade-Offs and Metabolic Cost Estimates

Genetic Edit	Intended Effect	Metabolic Cost	Stability Risk	Impact on Yield
hydA↑	↑ H₂ production	ATP, Fe-S cluster draw	Moderate	High
hydG∆	↑ Maturation	Low	Low	Medium
bchP↑	↑ Light absorption	NADPH usage	High	Medium
hupL∆	↓ H2 Loss	None	Very Low	Medium
crtB↑	↑ Light Stability	Membrane strain	Medium	Medium

4.3.1. Stack Compatibility and Synergy Limitations

While each proposed genetic edit in isolation has been documented to improve hydrogen yield or metabolic efficiency, the assumption of full functional compatibility across a multi-gene stack presents a significant risk of overestimation. In real-world systems, combinatorial gene stacking beyond 3–4 edits frequently result in metabolic crosstalk, resource exhaustion, and unexpected feedback inhibition.

For facultative anaerobes like Clostridium butyricum and Rhodobacter sphaeroides, regulatory networks often suppress overexpression to maintain redox balance and cellular homeostasis. Literature suggests that above 4–5 simultaneous edits, constructive collaboration tends to plateau—and in some cases, becomes negative due to ribosomal bottlenecking, cofactor depletion, or transcript interference.

Therefore, while this model includes a 6–8 gene modification set for theoretical benchmarking, it is explicitly treated as a ceiling scenario. Actual implementation may require modular optimization or sequential engineering to preserve system stability. Future modelling efforts should incorporate a constructive collaboration correction factor or simulate dynamic trade-offs using flux balance or gene regulatory network models to determine yield inflection points under escalating metabolic burden.

4.4. Gene Stack Resource Load and Metabolic Burden Quantification

While the previous sections described genetic modifications in Clostridium butyricum and Rhodobacter sphaeroides to enhance hydrogen production, this subsection addresses the metabolic and biochemical resource requirements imposed by these edits. Quantifying these demands is essential to assess whether the engineered strains could sustain high-yield hydrogen production without exceeding physiological tolerances or collapsing under resource depletion.

• Gene Stack Limit Curve

Genes Added	Expected Yield Boost (%)	Stability Risk
1–2	+20–40%	Low
3–4	+40–70%	Medium
5–6	+70–90%	High
7+	Marginal or Negative	Very High

• Note on Enzyme Yield Non-Linearity:

While this model assumes that overexpression of hydrogenases and associated enzymes leads to enhanced hydrogen yield, it is important to note that enzyme activity does not scale linearly with gene expression levels. In practice, increasing transcription may not proportionally increase catalytic turnover due to substrate limitation, cofactor saturation, and feedback regulation. For example, hydrogenase activity may plateau at high expression levels if reduced ferredoxin or NAD(P)H becomes limiting. Additionally, overexpression may lead to improper protein folding or inclusion body formation, particularly under metabolic stress. Therefore, the relationship between gene dosage and hydrogen yield is subject to diminishing returns beyond certain thresholds, and this model's assumptions should be interpreted as reflecting best-case catalytic availability rather than guaranteed performance scaling.

• Gene Stack Correction Factor (GSCF) Introduction

While this model assumes stable expression of a 6–8 gene enhancement set, real-world systems often exhibit diminished returns beyond 4 simultaneous edits due to transcriptional overload, ribosomal crowding, and cofactor scarcity. To bridge the gap between theoretical performance and practical expression limits, we introduce a Gene Stack Correction Factor (GSCF) as a heuristic yield modifier. This factor reflects the non-linear yield penalties associated with stacked edits and can be applied multiplicatively to the uncorrected hydrogen yield. Based on published microbial engineering studies, a GSCF of 0.6–0.9 is proposed depending on stack complexity and metabolic burden. For example, if the system theoretically produces 12 mol H₂ per cycle, a GSCF of 0.75 would reduce the adjusted yield to 9 mol—a more biologically plausible outcome under high-expression load. This correction aligns the model with empirical observations while preserving its utility as a theoretical ceiling. All future system designs should benchmark gene stack size against yield plateau thresholds to avoid metabolic collapse or regulatory interference.

4.4.1. ATP Demand from Nitrogenase Activity

The photofermentation stage, driven by R. sphaeroides, relies on nitrogenase to catalyse hydrogen evolution from water and VFAs. Nitrogenase is ATP-intensive, requiring approximately 16–24 ATP per mole of H₂ produced. This model adopts a midpoint estimate of 20 ATP/mol H₂, resulting in a total demand of 200 ATP per mole of glucose processed (given 10 mol H₂ from photofermentation).

While some of this ATP is regenerated via acetate metabolism (enhanced by ackA insertion) and cyclic photophosphorylation, sustained operation requires a robust intracellular ATP supply and efficient cofactor cycling.

4.4.1.1. ATP Demand vs. ATP Supply — Balancing Constraints

Nitrogenase-driven hydrogen evolution imposes a heavy ATP cost—estimated here at 20 ATP per mole of H₂, resulting in a demand of 200 mol ATP per cycle to produce 10 mol H₂ in the photofermentative stage. While the model assumes that ATP is regenerated via acetate metabolism (enhanced through ackA insertion) and cyclic photophosphorylation, this assumption represents a best-case energy budget that is rarely achievable without active metabolic steering.

Acetate metabolism contributes an estimated 1–2 ATP per mol acetate, and photophosphorylation offers 2–4 ATP per mol H₂ at optimal photon flux. However, real-world systems often experience suboptimal light exposure, cofactor scarcity, or carbon flux diversion—resulting in partial ATP regeneration. A conservative integration of these yields suggests a probable output of only 100–150 mol ATP per cycle, leaving a shortfall of 25–50%.

This ATP deficit would likely throttle nitrogenase activity, delaying or capping photofermentative hydrogen yield. If ATP availability drops below critical thresholds, not only does H₂ evolution slow, but redox balancing and cofactor recycling are also disrupted—introducing systemic instability and potential feedback inhibition.

A corrected model scenario assuming only 150 mol ATP available would reduce photofermentative H₂ production proportionally—from 10 mol down to ~7.5 mol, decreasing the total cycle yield from 12 mol to ~9.5 mol.

Future iterations of this model should integrate dynamic ATP flux simulations, using flux balance analysis (FBA) or constraint-based modelling to quantify enzyme activity throttling under resource-limited conditions. Until then, this section presents an upper-bound yield scenario, with clearly defined ATP bottlenecks as a realism modifier.

Clarification on Cyclic Photophosphorylation ATP Output:

The estimated ATP output from cyclic photophosphorylation (2–4 mol ATP per mol H₂) is based on indirect analogues from purple non-sulphur bacteria and phototrophic growth models under nitrogen-limited conditions. Empirical studies, such as those by Gest and Kamen (1993) and Madigan et al. (2010), suggest that cyclic electron flow through the photosynthetic apparatus in Rhodobacter sphaeroides can regenerate ATP at rates sufficient to support nitrogenase activity when paired with acetate metabolism. While exact stoichiometries vary due to differences in proton motive force efficiency and ATP synthase coupling ratios, ATP yields in this range have been inferred from nitrogen-fixation-associated photophosphorylation experiments. Therefore, the model adopts this bracketed estimate as a biologically plausible range under ideal light intensity and uninterrupted electron flow conditions. Further refinement would require kinetic simulations of electron transport rates and intracellular ATP flux.

References available on request or in supplementary materials.

A conservative ATP accounting is as follows:

Source	Estimated ATP Yield	Comments
Acetate pathway (ackA/pta)	~1–2 mol ATP per mol acetate	Limited by carbon flux and available acetate; assumes full routing from glucose
Cyclic photophosphorylation	Variable (2–4 ATP/mol H₂ at best)	Highly dependent on light intensity and electron transport efficiency
Total Expected ATP	~100–150 mol (upper-bound)	Potentially insufficient under realistic light or metabolic conditions
Required ATP	200 mol	For 10 mol H₂ at 20 ATP/mol H₂

Gap: Up to 50–100 mol ATP shortfall may occur if substrate-level phosphorylation or light-driven ATP regeneration is suboptimal.

Implication: Without precise control of metabolic routing and light efficiency, nitrogenase activity may be throttled due to intracellular ATP scarcity. This would reduce hydrogen yield and/or trigger stress responses that disrupt expression stability.

A full ATP budget under dynamic conditions should be incorporated into future models using flux balance analysis (FBA) or kinetic simulations. For now, the system’s ATP sufficiency is treated as a limiting factor, and actual hydrogen yields may fall below the theoretical ceiling unless additional ATP-regenerating strategies are employed (e.g., enhanced cyclic electron flow, acetate supplementation, or heterologous ATP synthase integration).

Realistic ATP-Limited Yield Adjustment

Based on current pathway estimates, the photofermentation stage demands ~200 mol ATP per cycle to drive nitrogenase-catalysed hydrogen evolution. However, the model's ATP supply—via acetate metabolism and cyclic photophosphorylation—realistically maxes out between 100–150 mol ATP per cycle under optimal conditions. This leaves a projected ATP deficit of 25–50%, which proportionally throttles nitrogenase activity. Assuming a linear ATP-to-H₂ coupling (20 ATP/mol H₂), this implies a revised photofermentative hydrogen yield of 7.5–9 mol H₂, depending on metabolic routing and light conditions.

Accordingly, the system’s total hydrogen yield per cycle drops from 12 mol to ~9.5–11 mol, once ATP availability is considered. Future models should integrate ATP-coupled flux penalties directly into kinetic simulations to capture this constraint more dynamically. Until then, the current model's full-yield scenario should be interpreted as a theoretical ceiling contingent on ideal energy regeneration, uninterrupted light flux, and acetate availability. Any deviation from these assumptions risks substantial yield collapse due to ATP bottlenecks.

4.4.2. Fe-S Cluster and Metal Cofactor Requirements

Both [FeFe]-hydrogenase (overexpressed via hydA↑) and nitrogenase require iron-sulphur (Fe-S) clusters, with nitrogenase additionally requiring molybdenum. The biosynthesis of these cofactors is ATP-dependent and sensitive to intracellular metal availability.

Enzyme	Fe Requirement	Key Cofactor
[Fe-Fe]-hydrogenase	2-3 Fe atoms per unit	Fe-S clusters
Nitrogenase (NifDK)	~8 Fe, 7 S, 1 Mo	Fe-S MoFe cofactor

Without sufficient Fe²⁺ and MoO₄²⁻ supplementation, expression levels will bottleneck. Real-world microbial cultures often experience yield collapse when metal availability becomes rate-limiting, particularly in continuous systems.

4.4.2.1. Cofactor Limitations as a Performance Bottleneck

The model assumes sufficient supplementation of trace metals such as Fe²⁺ and MoO₄²⁻ to support the biosynthesis of Fe-S clusters and the FeMo cofactor required for [FeFe]-hydrogenase and nitrogenase, respectively. However, this oversimplifies the bioavailability kinetics and competitive uptake dynamics that occur in microbial cultures.

In continuous bioreactors or batch systems with high expression burdens, trace metals often become rate-limiting not due to absolute absence—but due to poor solubility, ion competition, or metabolic uptake saturation. For instance, Fe²⁺ must compete with Mn²⁺, Zn²⁺, and Cu²⁺ at membrane transporters, while molybdate (MoO₄²⁻) can be competitively inhibited by sulphate or phosphate anions.

Additionally, Fe-S cluster assembly is an ATP-dependent, protein-mediated process requiring coordination with HydE/F/G maturation pathways. Any disruption—whether from cofactor scarcity, oxidative stress, or pH drift—can stall hydrogenase/nitrogenase assembly and sharply reduce yield.

This model does not yet simulate:

• Transporter affinity thresholds for metal uptake
• Chelation kinetics in complex media
• Metal precipitation risks under suboptimal pH or redox conditions

Future models should integrate trace metal mass balances, transporter saturation kinetics, and bioavailability correction factors—especially for Fe and Mo—to predict system yield more accurately under industrial or continuous-use scenarios. For now, metal availability is flagged as a potential hidden bottleneck, capable of undermining even fully expressed enzyme systems.

Both enzymes rely on complex and resource-intensive metalloprotein assembly pathways:

Enzyme	Cofactor Requirement	Biosynthetic Complexity
[FeFe]-Hydrogenase	2–3 Fe atoms per unit; Fe-S clusters	Dependent on functional HydEFG maturation proteins and iron homeostasis
Nitrogenase (NifDK)	~8 Fe, 7 S, 1 Mo per unit	Requires full Fe-S cluster biogenesis and molybdenum incorporation via NifEN pathway

Cofactor Availability and Economic Feasibility Note: Molybdenum (Mo) is an essential cofactor for nitrogenase activity, typically incorporated into the FeMo-cofactor cluster at a rate of ~1 mol Mo per mol of active nitrogenase. For the modelled system—assuming full activation of nitrogenase to produce 10 mol H₂ per cycle—this corresponds to a theoretical requirement of approximately 1 mmol Mo per mole of glucose processed. At current bulk pricing for molybdate salts (e.g., ammonium molybdate), this equates to a cost of less than $0.50 per batch (based on $40 per 100 g, molar mass ≈ 123.9 g/mol). Given its industrial availability, low cost, and widespread use in agriculture, molybdenum does not present a logistical or economic bottleneck for system deployment. However, its bioavailability must be managed carefully to avoid toxicity, and real-world systems may benefit from trace metal recycling strategies to maintain sustainable operation.

4.4.3. NAD(P)H Demand and Redox Balancing

Overexpression of bchP and crtB in R. sphaeroides boosts photosystem formation and light-harvesting efficiency but increases NADPH consumption. Biosynthesis of bacteriochlorophylls and carotenoids requires multiple NADPH molecules per pigment unit, particularly during biofilm establishment phases. This can cause early-stage redox imbalance and may necessitate transhydrogenase activity or metabolic reallocation.

Estimated NADPH diversion during initial pigment biosynthesis may approach 20–30% of the available pool, competing with core metabolic functions if unregulated.

4.4.4. Expression Load and Transcriptional Capacity

Stacking 6–8 genes under high-strength promoters increases the demand for RNA polymerase, ribosomes, and processing enzymes. Even with chromosomal integration, excessive constitutive expression can saturate transcriptional machinery, triggering unintended feedback loops or silencing effects. Promoter crowding, mRNA degradation limits, and translational burden must be considered.

A gene burden summary is provided below:

Gene	Resource Burden	ATP Demand	Cofactor Required	Risk Summary
hydA↑	High	Moderate	Fe-S	Critical but unstable under low Fe
crtB↑	Medium	Low	NADPH	Affects light-harvesting onset
bchP↑	High	High	NADPH	Major redox sink; timing-sensitive
ackA↑	Low	Net ATP gain	Acetate flux	Energy positive
hupLΔ	None	None	N/A	Eliminates H₂ reuptake
hydGΔ	Low	None	N/A	Enables hydA maturation

4.4.5. Systemic Load Implications

To support these expression demands, the model assumes:

• Supplemented growth media with Fe²⁺ (≥20 µM), trace Mo, and buffered NAD⁺/NADP⁺ regeneration capacity
• Stage-specific promoter design (e.g., light-inducible for bchP, stress-inducible for hydA) to stagger gene expression
• Potential inclusion of ATP-regenerating pathways or cyclic photophosphorylation enhancers

Failure to meet these conditions in practical systems could result in:

• Loss of expression fidelity over time
• Early plateauing of hydrogen yields due to cofactor exhaustion.
• Activation of stress response pathways or unintended mutations in high-burden operons

Author’s Note: While these burdens are quantified to the extent possible within a theoretical framework, their exact thresholds are highly context-dependent and would require dynamic flux balance analysis (dFBA) and proteomic data to model precisely. The intent here is to flag these risks for consideration in future experimental or in silico simulations.

• Personal Limitation Acknowledgement

I want to be fully transparent about this part of the model. While I have done my best to design and rationalize the proposed genetic modifications based on existing literature, I do not currently have the expertise—or the computational tools—to model the deeper layers of metabolic burden, cofactor imbalance, or regulatory drift that these changes would realistically trigger.

These things are real and complex, and I do not want to pretend otherwise. Modelling gene expression stability, metabolic flux distributions, or system-wide energy budgets would require advanced simulations like flux balance analysis (FBA), and that is well beyond my current skill set.

So, for now, I am treating these genetic edits as a best-case hypothetical scenario—not because biology is this clean, but because I wanted to explore the ceiling of what such a system might achieve if everything worked optimally.

• Future Requirement: Dynamic Flux and Cofactor Modelling

While this framework establishes a clear theoretical ceiling, it does not include dynamic modelling of intracellular resource allocation, cofactor competition, or redox balancing. Biological systems are governed by non-linear flux networks where enzyme activity, ATP/NAD(P)H levels, and trace metal availability form tightly regulated, feedback-coupled systems. The overexpression of hydrogenase, nitrogenase, and light-harvesting proteins places enormous demand on Fe-S clusters, ATP pools, and NAD(P)H regeneration, all of which interact in ways that static stoichiometric models cannot resolve.

To capture these constraints, future iterations of this model should incorporate flux balance analysis (FBA) or constraint-based modelling (CBM) using genome-scale metabolic reconstructions of C. butyricum and R. sphaeroides. Such tools would allow time-resolved prediction of cofactor depletion, ATP bottlenecks, and gene stack performance under varying substrate, light, and nutrient conditions. Until then, this paper offers an upper-bound performance envelope under the assumption of stable expression, non-limiting cofactors, and uninterrupted energy flux. These assumptions are flagged for future experimental validation or computational refinement.

While cofactor and redox dynamics are acknowledged as critical, they are beyond the current model’s static scope and are proposed for future constraint-based refinement.

5. Two-Stage Bioreactor Design and Structural Considerations

The proposed system is a two-stage, vertically integrated microbial hydrogen production reactor designed to optimize substrate conversion and gas recovery efficiency. It combines dark fermentation and photofermentation processes in a unified chassis, with integrated membrane-based gas separation positioned above the biological chambers. The system is conceptually treated as a synthetic microbial consortium, where cross-species metabolic routing and spatial compartmentalization are leveraged to maximize efficiency. The reactor model is developed at bench-scale for a single-cycle processing of 1 mol glucose (~180 g), representing a conceptual prototype volume of approximately 1–2 litres.

5.1. Structural Layout

The upper chamber houses a scaffold-supported biofilm of Clostridium butyricum, maintained under strict anaerobic and dark conditions. This chamber receives glucose-rich substrate and facilitates dark fermentation, generating hydrogen, carbon dioxide, and volatile fatty acids (VFAs)—primarily butyrate and acetate. VFAs then flow downward via gravity through a perforated partition into the lower chamber, minimizing the need for pumping. To prevent VFA stagnation and nutrient layering, the system incorporates intermittent low-power peristaltic mixing, allowing fluid recirculation without disturbing biofilm structure. This gentle agitation also redistributes thermal gradients and ensures substrate availability across chamber depth.

System Flow Schematic (Verbal Representation):

• [Glucose Input]
• ↓
• [Dark Fermentation Chamber – C. butyricum]
• ↓ VFAs via gravity flow
• [Photofermentation Chamber – R. sphaeroides]
• ↓
• [Shared Headspace with Gas Collection]
• ↓
• [Vacuum-Assisted Membrane Separator]
• ↙ ↘
• [H₂ Tank] [CO₂ Redirection Line]
• ↓
• [Sugarcane / Algae Sink]

(This simplified schematic outlines the metabolic and physical routing within the system. A visual diagram based on this structure is recommended for future iterations.)

The lower chamber contains a photofermentative biofilm of Rhodobacter sphaeroides, grown on light-permeable support structures. This zone is illuminated continuously with low-intensity light (~30 µmol/m²/s), activating nitrogenase-dependent photofermentation of VFAs into additional hydrogen and carbon dioxide. To prevent uneven light distribution across biofilm surfaces, the lower chamber incorporates multi-angle internal LED arrays and limits biofilm panel thickness to <1 mm. This configuration reduces shadow zones, maintains consistent nitrogenase activation, and mitigates attenuation effects common in optically dense biofilms. LED units are housed in sealed, transparent enclosures integrated within the chamber walls to preserve anaerobic integrity. Thermal insulation prevents local convection currents and protects microbial gradients, while transparent polycarbonate barriers allow uninterrupted light transmission. Both chambers are equipped with separate nutrient ports and pH stabilization controls.

5.1.1. Thermal Stability and Heat Load Control

Although the reactor is designed for short-cycle operation (≤72 hours), metabolic heat generation and continuous photofermentation lighting introduce non-negligible thermal stress. Both microbial systems—Clostridium butyricum and Rhodobacter sphaeroides—exhibit peak hydrogen productivity at temperatures between 28–32°C, with yield drop-offs and enzymatic instability observed above ~35°C.

This model assumes passive heat dissipation through structural conduction and ambient thermal diffusion. However, realistic implementation may require active thermal stabilization to maintain optimal internal temperatures. Potential thermal management solutions include:

• Passive Measures:

  ○

  High-conductivity aluminium or copper heat-sinks

  ○

  Surface-mounted thermal paste layers or graphene cooling films

  ○

  Natural convection via side-ventilation grills

• Active Measures:

  ○

  Thermoelectric (Peltier) cooling plates affixed to outer reactor walls.

  ○

  External water-jacket circulation systems for continuous heat removal

  ○

  Embedded thermistor array with automated fan/pump triggers

A simplified heat generation estimate based on LED power input and microbial activity suggests potential internal gains of 1.5–3.5°C per 24 hours if unmanaged. While this system remains theoretically stable over one batch cycle, future scale-ups or continuous operation will necessitate formal thermal modelling and fail-safes to prevent enzyme denaturation and biofilm instability.

5.1.2. Thermal Load Estimate (Bench-Scale)

To quantify potential internal temperature rise, a simplified thermal load estimate was performed using input power values from key reactor components and known heat capacity values of water.

Assumptions:

• LED system = 5 Watts continuous (assume 100% of power converted to heat)
• Control systems and minor mixing = 2 Watts average thermal load
• Total thermal input power ≈ 7 W
• System runtime: 24 hours = 86,400 seconds
• Reactor fluid volume = 2 L = 2000 g (assuming density ≈ water)
• Specific heat capacity of water (c) = 4.18 J/g°C
• Step 1: Total Heat Energy Generated

$$Q=P\times t=7W\times \mathrm{86,400}s=\mathrm{604,800}J$$

• Step 2: Estimate Temperature Rise

$$∆T={\displaystyle \frac{Q}{m\times c}}={\displaystyle \frac{\mathrm{604,800}J}{2000g\times 4.18J/g°C}}\approx 72.3°C$$

This represents maximum uncontrolled heating if no heat escapes. In practice, natural convection, conduction through reactor walls, and evaporative cooling remove >95% of that heat at bench scale.

Conservative estimate of retained thermal energy: 2–5%

$$∆T_{realistic}\approx 1.4-3.6℃/day$$

This supports the prior claim of ~1.5–3.5°C rise per 24 hours, aligning with published values in similar low-power biofilm systems.

5.2. Biofilm Limitations and Realism Considerations

While this model assumes biofilm layers remain structurally stable and spatially uniform across both fermentation stages, real-world systems experience significant challenges in maintaining long-term biofilm integrity. Issues such as detachment, overgrowth, nutrient gradient formation, and light attenuation (particularly in the photofermentation stage) are well-documented. Biofilm clogging may disrupt fluid and VFA flow between chambers, while irregular microbial growth could reduce gas exchange efficiency. Additionally, biofilms often require scheduled cleaning or regenerative protocols to maintain performance in continuous systems.

These limitations are not included in this theoretical model to preserve stoichiometric clarity and kinetic tractability. However, future system design must address biofilm management through structural scaffolding, flow control strategies, or cyclic biofilm replacement. This remains one of the most critical engineering challenges in translating microbial hydrogen systems from bench-scale to continuous industrial operation.

The reactor’s design supports modular biofilm panel exchange, allowing saturated or detached panels to be manually or mechanically rotated out on a regular schedule. This approach minimizes downtime and supports continuous operation, especially in pilot-scale or industrial implementations.

5.2.1. Biofilm Degradation Timing and Panel Replacement Rates:

Although this model assumes spatially uniform and stable biofilm performance throughout each fermentation cycle, practical systems must contend with gradual biofilm degradation, metabolic fatigue, and physical detachment. Biofilm productivity typically follows a bell-curve pattern over time—rising during initial colonization, peaking at full metabolic activity, and declining as nutrient gradients, waste accumulation, or quorum signalling trigger senescence or lysis. In Clostridium and Rhodobacter systems, peak productivity is often sustained for 48–72 hours before measurable drops in hydrogen output occur.

Based on literature from continuous-flow and batch bioreactors, biofilm panels may require replacement or regeneration every 2–4 operational cycles (approx. 6–10 days) to maintain consistent gas yields. Productivity loss is typically inferred from declining gas evolution rates or increasing metabolite residue buildup. In future implementations, real-time monitoring systems such as inline gas flowmeters or capacitive biofilm thickness probes could be used to trigger automated panel rotation schedules. These interventions would minimize downtime and reduce the risk of clogging or metabolic collapse in long-term setups. For this theoretical model, the system is assumed to operate within the optimal productivity window, with biofilm exchange occurring post-cycle.

5.2.2. Biofilm Surface Area Exposure and Photonic Limitation

While the reactor model emphasizes volumetric throughput (e.g., mol H₂ per litre per day), the true limiting factor in photofermentation is not volume but surface area exposed to light. Unlike dark fermentation, photofermentative hydrogen production by Rhodobacter sphaeroides depends on the availability of photons to activate nitrogenase and light-harvesting complexes.

Even under continuous illumination, light attenuation in microbial biofilms becomes a serious constraint beyond ~1 mm thickness due to absorption, scattering, and shadowing effects.

Estimated Photon Requirements:

Based on empirical data:

• Nitrogenase activation requires 30–50 µmol photons/m²/s for optimal function.
• Commercial LED panels produce ~100–200 µmol/m²/s at a proximity.
• If photofermentation is to yield 10 mol H₂ per glucose as proposed, then:

Approximate light demand = 300–500 µmol photons per second per mole of substrate

That translates to:

• Photon exposure area = 0.04–0.08 m² per cycle at minimum, assuming tight efficiency
• For your 2-litre system, this means a minimum of 0.02–0.04 m² per litre of reactor volume.
• Surface Area:Volume Ratio (SA:V) Targeting

To avoid photonic bottlenecking, the following design constraint should be integrated into reactor planning:

Target Surface Area ≥ 0.02 m² per litre of reactor volume, preferably ≥ 0.03 m² for safety margin.

For a 2-litre reactor, that’s:

∙

40–60 cm² of light-exposed biofilm area, ideally distributed evenly across the inner surface via panel or spiral scaffolds

∙

This value can be met via:

∙

Vertical light-permeable baffles

∙

Spiral microfilm scaffolds

∙

Rotating thin-panel arrays submerged in the photoreactor zone.

∙

Photon Flux and Surface Area Clarification

Photofermentation rates depend on surface-exposed biofilm area rather than total volume. The model uses a minimum light exposure requirement of 0.02–0.03 m² per litre of reactor volume to avoid nitrogenase inhibition.

Added standard metric:

Areal light exposure: 0.04 m² per 2 L reactor (meets required 0.02 m²/L threshold)

Ensure areal productivity is always reported alongside volumetric values to avoid misleading interpretations.

Mitigation of Shadow Zones

To ensure consistent light activation:

∙

Multi-angle LED arrays should be positioned to eliminate angular blind spots.

∙

Panel spacing must allow reflected light and convection to pass through.

∙

Light diffusers or frosted panels may help even out photon flux.

∙

Implication on Reported Yield Metrics

If the surface area per litre drops below threshold:

• Nitrogenase activation becomes suboptimal.
• Effective H₂ yield drops, even if VFA supply and reactor volume are sufficient.
• Reported 10 mol H₂ per mol glucose becomes non-physical due to light starvation.

To reflect this:

Yield reporting should always be normalized per m² of light-exposed surface for photofermentation stages. For example:

Metric	Value
H₂ Productivity (volumetric)	44.85 L H₂/L/day
H₂ Productivity (areal)	~1120 L H₂/m²/day (at 0.04 m² per 2 L)

• Photonic Constraint and Surface Area-Linked Yield Limitations

While volumetric hydrogen productivity (e.g., L H₂ per L reactor per day) remains a common benchmarking metric, it can be misleading in photofermentation systems, where light availability—not reactor volume—is the true limiting factor. Nitrogenase activity in Rhodobacter sphaeroides is directly dependent on continuous photon flux. In systems where surface area is insufficient, biofilm zones become light-starved, reducing effective hydrogen yield despite substrate availability.

To ensure photon flux supports the theoretical 10 mol H₂ per mole of glucose, a minimum light-exposed surface area of 0.02–0.03 m² per litre is required. If surface area falls below this threshold, nitrogenase activation becomes suboptimal and the actual H₂ yield may drop well below volumetric predictions. This model assumes light saturation is achieved across all productive surfaces, but future reporting should always include:

• Volumetric yield (L H₂/L/day)
• Areal yield (L H₂/m²/day)
• SA:V ratio (m² per L)

These constraints prevent overestimation of system output and align the theoretical ceiling with physically enforceable design parameters. A surface area shortfall should be considered a hard cap on hydrogen productivity, regardless of substrate excess or reactor scale.

Author's Note:

While volumetric normalization is traditional in fermentation literature, it can conceal photofermentation failure modes. By including a surface-area constraint, this model avoids a common pitfall in hydrogen system design—if light just “gets in there” without mapping where, how, or if it lands.

5.3. Gas Capture and Flow Control

Gas production discussed here follows biochemical conversion within the structured reactor system detailed in Section 5.1 and 5.2.

Gas produced in both chambers accumulates in a shared headspace connected to a vacuum-assisted membrane separation system. The membrane system is assumed to be positioned externally above the bioreactor chambers and isolated from direct microbial contact to avoid fouling or biofilm disruption. Gas is extracted from the headspace without disturbing pressure conditions inside the bioreactor. Hydrogen is separated and collected as product, while carbon dioxide is directed toward external fixation systems as detailed in the CO₂ redirection strategy.

5.3.1. Hydrogen Accumulation and Feedback Inhibition

While this model acknowledges that Clostridium butyricum and Rhodobacter sphaeroides exhibit feedback inhibition at H₂ partial pressures exceeding 0.6–0.8 atm, it does not currently simulate the rate of gas accumulation or the membrane system's throughput capacity. This omission creates a dangerous assumption: that gas removal is effectively instantaneous and continuous.

The kinetic mismatch between hydrogen generation and membrane diffusion is a critical design vulnerability. Based on a production rate of 0.167 mol H₂/hr and a 0.5 L gas headspace, the partial pressure threshold for inhibition is reached in under 5 minutes without active extraction. This translates to a required membrane throughput of at least 2.7 mmol/min just to maintain sub-inhibitory headspace conditions.

Yet the model does not specify:

• Membrane surface area (m²)
• Permeability rate (Nm³/m²/h/bar)
• Vacuum cycle timing or capacity.

Without these values, it is impossible to guarantee that the system avoids enzyme suppression. Real membranes exhibit finite flux governed by Fick's Law, and polymer-MOF hybrids typically achieve 0.1–0.5 Nm³/m²/h/bar at best. To match the required flow, even a bench-scale system may need 10–20 cm² of active membrane area under steady vacuum—a nontrivial engineering constraint.

Dynamic Headspace Model and Membrane Throughput Limitations

To simulate enzyme viability more accurately during gas accumulation, a time-dependent headspace pressure model should be integrated into future iterations. Under standard fermentation rates (~0.167 mol H₂/hr), hydrogen accumulation in a 0.5 L headspace reaches inhibitory partial pressures (>0.6 atm) in under five minutes. This imposes a minimum membrane throughput requirement of ~2.7 mmol/min to maintain sub-inhibitory levels. Real membrane flux rates, governed by Fick’s law, typically range between 0.1–0.5 Nm³/m²/h/bar for MOF-polymer hybrids, implying a required membrane area of 10–20 cm² under vacuum-assisted flow. The current model assumes continuous, lossless gas extraction—an unrealistic constraint at bench or pilot scale.

Therefore, we propose a Membrane Flux Correction Factor (MFCF) to account for suboptimal gas clearance. For conservative benchmarking, yields in future models should be scaled by 0.7–0.85 depending on membrane permeability, fouling, and system response time. Real-time pressure regulation and multi-stage gas extraction loops may be required to prevent inhibitory backpressure during high-output periods. Until such a kinetic gas module is developed, all uncorrected H₂ productivity estimates in this paper should be treated as theoretical ceilings contingent on perfect ventilation and membrane performance.

Furthermore, extraction delay or vacuum cycling lag could create pressure pulses, briefly pushing partial pressure above the enzymatic inhibition point, thereby throttling yield or triggering stress responses.

• Model Extension Recommendation: Headspace Accumulation Kinetics

Future iterations of this model should incorporate a dynamic gas accumulation module, simulating pressure rise curves and matching them against membrane flow capacities. Until then, current H₂ yield estimates may be significantly inflated due to unmodeled feedback inhibition risks during periods of poor headspace ventilation.

• Representative Membrane Materials and Flux Parameters for Hydrogen Separation

This table presents example membrane types, flux characteristics, and performance constraints under vacuum-assisted operation. Values are normalized to STP conditions and 1 atm pressure differential. These benchmarks inform the Membrane Flux Correction Factor (MFCF) applied to theoretical gas clearance rates.

Membrane Type	Material System	Flux (Nm³/m²/h/bar)	Selectivity (H₂/CO₂)	Notes
MOF-Polymer Composite	ZIF-8 in Pebax matrix	0.25 – 0.50	~20 – 30	High-performance lab-scale membranes; sensitive to fouling
Polyimide	Commercial polymer	0.15 – 0.30	~10 – 15	Widely available, thermally stable, moderate selectivity
Silica-Based	Microporous glass	0.10 – 0.20	~5 – 8	Chemically stable, lower selectivity, brittle
PTFE / ePTFE	Expanded Teflon	0.05 – 0.10	~2 – 5	Durable, biofouling-resistant, low separation capacity

5.3.2. Control Systems and Feedback Architecture

Real-time control of environmental parameters is essential for maintaining the functional stability of the dual-stage microbial hydrogen production system. While earlier sections describe a theoretically ideal bioreactor with constant substrate diffusion, light penetration, and anaerobic integrity, real systems require feedback loops to compensate for internal fluctuations, metabolic drift, and external stressors.

This section outlines a modular control strategy designed to monitor and correct deviations in temperature, pH, light intensity, gas accumulation, and CO₂ sequestration rates. All controls are assumed to operate on low-power microcontroller-based platforms (e.g., Arduino, Raspberry Pi) using standard biosensor interfaces.

5.3.3. Quantified Hydrogen Accumulation and Inhibition Thresholds:

The inhibitory effects of hydrogen accumulation on enzymatic activity are well-documented in microbial hydrogen systems. Hydrogenase in Clostridium butyricum and nitrogenase in Rhodobacter sphaeroides exhibit feedback suppression at elevated hydrogen partial pressures, typically around 0.6–0.8 atm. Under the current reactor model, which assumes approximately 12 mol of hydrogen produced over a 72-hour cycle (~0.167 mol/hr), this threshold is reached rapidly if no gas removal occurs.

For a 2-litre reactor with a 0.5-litre gas headspace—a standard design assumption—hydrogen accumulation reaches 0.6 atm partial pressure in less than 5 minutes of continuous production. This corresponds to the accumulation of just ~13.4 mmol of H₂, based on ideal gas calculations at STP (22.4 L/mol). Without active gas extraction or continuous vacuum cycling, the system would quickly enter a regime of enzymatic inhibition, reducing hydrogen yield and potentially inducing redox or energy stress in the microbial population.

This analysis highlights the critical importance of real-time gas removal or pressure-controlled headspace venting in batch and continuous systems. Future iterations of this model should incorporate hydrogen accumulation kinetics and headspace pressure dynamics to simulate enzyme behaviour more accurately under varying operational conditions.

Supporting Calculation: Time to Reach Inhibitory H₂ Pressure

Given:

∙

Hydrogen production rate: 12 mol / 72 hr = 0.167 mol/hr

∙

Reactor headspace volume: 0.5 L

∙

Inhibitory partial pressure threshold: 0.6 atm

∙

Molar volume at STP: 22.4 L/mol

∙

Required moles of H₂ to reach 0.6 atm in 0.5 L headspace:

$$n={\displaystyle \frac{P\times V}{RT}}\approx {\displaystyle \frac{0.6\times 0.5}{22.4}}\approx 0.0134mol$$

∙

Time to reach 0.6 atm at constant production:

$$t={\displaystyle \frac{0.0134}{0.167}}\approx 0.08hour=4.8minutes$$

Thus, if hydrogen is not extracted or diffused out of the headspace, the system will exceed the inhibitory threshold in under 5 minutes of operation.

Side Note: This model assumes headspace management sufficient to prevent hydrogen accumulation above inhibitory thresholds. Future work will address dynamic gas control strategies to maintain enzymatic activity over extended operation.

5.3.3.1. Core Monitoring Parameters

The following table summarizes the critical control parameters, their sensor types, operating ranges, and programmed responses:

Parameter	Sensor Type	Target Range	Trigger Action
pH (both chambers)	pH probe (glass electrode)	6.2–7.4	<6.0 → alkaline buffer injection; >7.8 → mild acid pulse
Temperature	Thermistor array / RTD	28–32°C	>33°C → Peltier cooling plate or fan activation
Light Intensity (photo chamber)	Photodiode array	25–40 µmol/m²/s	<20 → LED boost; >50 → dimming circuit
H₂ Pressure (headspace)	MEMS pressure sensor	<1.1 atm	>1.2 atm → vacuum pump activation
CO₂ Accumulation (in pod)	NDIR gas sensor	<900 ppm	>1000 ppm → alert or temporary gas flow halt
Biofilm Thickness (optional)	Capacitive displacement probe	<1mm	>1.2 mm → alert for manual panel rotation

These sensors ensure system stability without requiring full laboratory automation. All thresholds are programmable and can be adapted per strain tolerance or reactor scale.

5.3.3.2. Control Logic Flow and Safety Behaviour

Each parameter operates within a hierarchical control loop, where low-severity deviations trigger passive adjustments (e.g., LED dimming), while high-severity deviations initiate active fail-safes (e.g., CO₂ rerouting or system shutdown).

Example Control Chain:

1.

Normal: Temp = 30°C → stable

2.

Warning: Temp = 33°C → Fan or TEC turns on

3.

Critical: Temp = 36°C → LED cutoff + halt fermentation → log error

3.

Fail-safe behaviour includes:

∙

Auto-throttling of gas flow if CO₂ sensors detect oversaturation in the sequestration pod.

∙

Redundant pH control with dual dosing pumps to prevent sensor drift-induced pH crash.

∙

Real-time log file generation for offline diagnostics

∙

Optional software modules can include:

∙

PID (proportional–integral–derivative) control for smooth system corrections

∙

Simple logic-based control (if/then) for robustness in off-grid deployments

∙

GSM or Wi-Fi module for remote alerting and telemetry (for field setups)

5.3.3.3. Suggested Control System Schematic (Textual)

Here is a verbal schematic you can translate into a diagram later:

[Biofilm Reactor]

↓

[Sensor Hub]

↓

[Microcontroller Unit]

↓

[Actuator Array]

↓

[Threshold Logic + Timer]

↓

[Data Log]

↑

Feedback Loop

Component Descriptions

∙

Sensor Hub: pH, temperature, gas (H₂/CO₂), light, biofilm thickness

∙

Microcontroller: Control firmware (e.g., Arduino, Pi)

∙

Actuator Array: Pumps, solenoids, LEDs, cooling fans, TEGs

∙

Threshold Logic + Timer: Evaluates state changes; executes correction routines.

∙

Data Log: Time-stamped CSV or onboard flash storage

5.3.3.4. Deployment Note

While these control systems are low-cost and modular at bench scale, scaling up would require industrial-grade PLCs (programmable logic controllers), multi-channel sensor arrays, and real-time process monitoring dashboards. For field or developing-world setups, simplicity and resilience are prioritized over digital optimization.

5.4. Assumptions and Limitations

• Uniform substrate diffusion and metabolic activity are assumed across both biofilms. Nutrient gradients, quorum dynamics, and spatial heterogeneity are excluded from this model.
• No cross-contamination or metabolic inhibition is assumed between microbial zones.
• Light distribution in the photofermentation chamber is considered sufficient and non-limiting.
• Vacuum extraction system is modelled as non-intrusive and energetically lossless, though this remains a significant engineering challenge.

Author's Note:

A theoretical system cannot be called "complete" if it assumes perfect environmental stability and ignores control infrastructure. These feedback systems not only maintain biological integrity but also function as ethical enforcers for the CO₂ redirection requirement, ensuring that industrial misuse cannot bypass sustainability protocols without triggering visible alarms or system throttling.

6. Gas Separation and Post Processing

The dual-stage microbial hydrogen production process yields a gaseous mixture composed primarily of hydrogen (H₂) and carbon dioxide (CO₂), which must be separated efficiently to preserve product purity and enable sustainable downstream processing. This section outlines the theoretical gas separation mechanism employed in the model and its associated assumptions.

6.1. Membrane Separation

The system uses a vacuum-assisted membrane separation module designed to selectively extract hydrogen from the fermentation headspace. Hydrogen is extracted first due to its high diffusivity and favourable permeation properties in most polymer membranes. While molecular size plays a role, effective separation depends on the selectivity of the membrane material and the partial pressure differential across it. For this model, two classes of membranes are considered:

• Polyimide-based membranes: Common in gas separation, offering moderate H₂/CO₂ selectivity and high mechanical stability.
• Metal-Organic Frameworks (MOFs): Emerging materials with tuneable pore sizes and potentially higher selectivity, modelled here as a theoretical performance benchmark [5].

Vacuum pressure is applied across the membrane surface to enhance gas flux. The system is assumed to operate under thermally stable, anaerobic conditions to prevent oxidation and enzymatic inhibition. While modelled here as perfectly staged and non-overlapping, each membrane step in real systems results in partial gas mixing. Recirculation loops, pressure swing adsorption, or cryogenic polishing may be required to achieve the purity assumed in the current idealized model.

6.1.1. Membrane Lifespan, Degradation, and Fouling Considerations

The gas separation system in this model utilizes semi-permeable vacuum-assisted membranes for hydrogen and CO₂ extraction. While theoretical performance assumes >95% separation efficiency, practical implementations face performance drift due to material degradation, physical stress, and fouling from microbial byproducts.

Key degradation and fouling risks include:

• Protein Mist Fouling: Fine droplets from active fermentation can condense on membranes, blocking pores and reducing flow rate.
• Salt Crystallization: Trace ionic species in the broth may vaporize and deposit on membranes, especially in long-cycle or high-CO₂ environments.
• Vacuum Fatigue: Repeated vacuum cycling induces microfractures and pore deformation, particularly in polymer-based membranes.
• Acidic Gas Corrosion: CO₂-rich gas streams can lower pH and erode polymeric layers over multiple cycles if not buffered.

Stressor	Primary Effect	Degradation Onset (Est.)	Mitigation
Biofilm mist / protein fog	Pore clogging	24–48 h	Pre-filtration mesh, vertical flow path
CO₂ acidification	Polymer erosion	5+ cycles	Use of fluoropolymer or ceramic layers
Vacuum cycling	Micro-fracturing	~30–50 cycles	Lower vacuum pulsing frequency
Salts / metabolites	Crystallization/fouling	48–72 h	Gas polishing (desiccant + filter)

To simplify modelling, this system assumes no membrane replacement within a 72-hour cycle, with a conservative performance loss margin of ≤5% factored into gas recovery estimates. For continuous or extended-cycle setups, performance decay should be modelled dynamically, and membrane health sensors or differential flow detection can be integrated for predictive maintenance.

Future iterations may benefit from pre-filtration chambers, dedicated gas polishing columns, and switchable membrane banks for redundancy.

6.1.2. Membrane Throughput Constraints and Gas Flux Limitations

While membrane-based gas separation systems offer scalable selectivity and moderate energy requirements, real-world performance is ultimately constrained by membrane area, gas flux rate, and partial pressure differentials. In the current model, gas extraction is treated as continuous and instantaneous, with all hydrogen produced (12 mol per cycle) assumed to be removed in real time. However, without quantifying membrane throughput (typically expressed in Nm³/m²/h/bar), the system risks hydrogen accumulation in the headspace—potentially surpassing the inhibitory threshold for hydrogenase and nitrogenase activity (~0.6 atm).

To prevent enzymatic suppression, future designs must benchmark membrane surface area, permeability coefficients, and operational pressure differentials to ensure adequate volumetric flow rates. If extraction lags gas production, internal feedback inhibition may compromise yield regardless of stoichiometric potential.

6.1.3. Gas Purity and Post-Separation Polishing

Although the current model assumes hydrogen purity suitable for energy applications, typical membrane separations alone do not meet the ≥99.97% purity threshold required for polymer electrolyte membrane fuel cells (PEMFCs). Even with a 92% selectivity rate, trace CO₂, water vapor, and residual metabolites may remain in the product stream. Realistic systems may therefore require secondary polishing units—such as pressure swing adsorption (PSA), catalytic oxidation beds, or activated carbon filters—to achieve fuel-grade standards. These components are omitted from this model but should be considered in lifecycle design for application-specific deployment.

6.1.4. Backpressure and Vacuum-Induced Stress Considerations

Vacuum-assisted separation introduces negative pressure differentials across the membrane and reactor headspace. While the model assumes this extraction occurs without disrupting internal microbial dynamics, real-world systems may experience shear stress on biofilm surfaces, altered gas solubility, or microbially relevant pressure fluctuations. These shifts could destabilize biofilm attachment or impair metabolic output. Future iterations should consider pressure dampeners, buffer chambers, or variable-flow regulators to mitigate these risks and preserve system integrity during extended operation.

6.2. CO₂ Routing and Sink Options

For theoretical yield benchmarking, this model assumes best-case performance for membrane gas separation. However, to contextualize feasibility, a realistic performance corridor is also described below.

• Hydrogen/CO₂ separation efficiency:Modelled at 100% for theoretical ceiling, but real-world membrane recovery rates range from 85–95% depending on material, pore stability, and backpressure control [Bernardo et al., 2009].
• Membrane lifespan:Idealized as indefinite, but real systems experience biofouling, pore compaction, and selectivity loss over time. Cleaning cycles or replacements every 300–500 hours are common.
• Vacuum extraction loss:Assumed lossless, but actual systems incur 5–15% efficiency penalties due to pump inefficiency, heat loss, or pulsed extraction behaviour.
• Product purity:Fuel-grade H₂ is assumed, but in practice polishing steps like pressure swing adsorption (PSA), activated carbon filters, or moisture traps may be required.

These constraints are not simulated in this model but must be considered for realistic deployment. Future extensions may incorporate fouling coefficients, separation efficiency decay curves, and energy cost scaling with purity thresholds.

6.2.1. Yield Impact from Separation Losses

To quantify the effect of real-world membrane inefficiencies on hydrogen recovery, a comparative energy yield analysis was conducted.

Assumptions:

• Theoretical H₂ production per cycle = 12 mol
• Combustion energy per mol H₂ = 286 kJ/mol
• Ideal separation efficiency = 100%
• Realistic separation efficiency = 92% (typical for polymer/MOF hybrid membranes)
• Energy loss from H₂ capture inefficiency:

$$Captured{H}_2=12mol\times 0.92=11.04mol$$

$$RecoveredEnergy=11.04mol\times 286kJ/mol\approx 3157.4kJ$$

• Compared to the ideal yield of 3432 kJ, this is a loss of 274.6 kJ, or ~8% energy penalty purely from imperfect gas separation.
• Impact on Net Energy Gain: Original net gain = 3432 − 960 = 2472 kJ
• Adjusted net gain = 3157.4 − 960 = 2197.4 kJ
• → Net energy drop: ~11.1%

This demonstrates that even modest membrane inefficiencies can meaningfully reduce system energy balance. For scaled systems, such losses may justify staged membrane banks or polishing units to preserve fuel-grade output.

6.3. System Losses and Gas Fate

This section assumes best-case membrane function for modelling clarity. In real-world systems, gas separation faces multiple engineering constraints including:

• Partial pressure gradients that limit flux across membranes.
• Membrane aging and fouling from bioaerosols or organic residues.
• Backpressure that could affect gas accumulation rates or microbial equilibrium.
• Vacuum system inefficiencies, such as pulsed flow, heat buildup, and non-linear extraction rates.

While these effects are not modelled here, they present clear avenues for future investigation. For instance, hydrogen recovery efficiencies in real systems often fall between 85–95%, and energy requirements for vacuum pumps scale non-linearly with extraction volume and purity targets. Additionally, pressure differential stress across biofilms could lead to microbial adhesion loss or diffusion barriers, requiring pressure dampening designs in future prototypes. These challenges reinforce the need for experimental membrane characterization, CFD simulation of gas flow, and pressure-regulated biofilm integration in scaled systems. Until such data becomes available, the present model defines an upper-bound, zero-loss gas separation baseline.

To further clarify the performance gap between the idealized membrane model and real-world system behaviour, this table summarizes critical membrane and vacuum gas handling parameters, with estimated realistic performance ranges and their impact on system operation.

Parameter	Ideal Model	Realistic Range	Impact
H₂ Separation Efficiency	100%	85-95%	Lowers product yield, energy recovery
Vacuum Pump Efficiency	100%	85-90%	Increases power draw
Membrane Lifetime	Infinite	300-500 Cycles	Affects maintenance and cost
Biofouling	None	Mild-Moderate	Requires cleaning cycles
Purity (H₂%)	100%	93-98%	May need PSA or polishing

7. CO₂ Fixation and Redirection Strategy

The dual-stage microbial system yields significant volumes of carbon dioxide as a metabolic byproduct—estimated at 6 mol CO₂ per mole of glucose, or approximately 264 g CO₂ per 180.16 g glucose, resulting in a CO₂-to-glucose mass ratio of 146.56%. This value exceeds the substrate mass due to oxygen atoms derived from water during photofermentation reactions. To maintain environmental sustainability, the model emphasizes redirecting this CO₂ to biological carbon sinks.

7.1. Candidate Biological Sinks

Three primary biological systems are proposed for CO₂ assimilation:

7.1.1. Sugarcane Plantations (C₄ Crops)

• High CO₂ fixation rates
• Dual purpose: sequestration + potential glucose recovery
• Proven viability in tropical and subtropical zones

7.1.2. Bamboo Forests

• Among the fastest-growing terrestrial plants
• High year-round carbon uptake per hectare
• Suitable for industrial land reclamation or buffer zones

7.1.3. Algae Bioreactors [8]

• High surface area to volume ratio for CO₂ uptake
• Suitable for urban, modular, or space-limited deployment
• Biomass can be harvested for nutrient recovery or biofertilizer.
• Other systems like aquatic plant ponds (e.g., duckweed) were explored but deprioritized due to acidification risks and high maintenance needs for pH stabilization.

7.2. Storage and Transport Assumptions

Separated CO₂ is assumed to be stored at ambient pressure in intermediate buffer tanks before redirection. No liquefaction or compression is modelled. Transport to biological sinks is assumed to occur via insulated pipelines or short-range pressure vessels depending on site layout and distance. These assumptions are idealized and exclude energy or infrastructure costs.

7.3. Scientific Frustrations and Realism Acknowledgement

While the paper outlines multiple CO₂ fixation strategies, achieving a truly closed-loop carbon system is far from straightforward. Each sink presents trade-offs:

• Land-based sinks (sugarcane, bamboo) demand significant space, irrigation, and land-use planning.
• Algal systems require strict environmental control (light, nutrients, mixing) and capital investment.
• Aquatic sinks, though attractive, suffer from ecosystem acidification without alkaline buffering.
• Moreover, the actual CO₂ uptake efficiency of each system is highly variable and context dependent. Climate, crop health, microbial symbiosis, and system maintenance all contribute to sequestration inconsistency.

This section represents a combination of plausible strategies and educated speculation. While redirecting CO₂ biologically is possible, it is also clear that implementation would require site-specific optimization, lifecycle cost analysis, and multi-disciplinary collaboration. The paper emphasizes this limitation as one of the most significant engineering and ecological challenges in making microbial hydrogen production truly carbon neutral.

7.3.1. Temporal and Lifecycle Constraints in CO₂ Fixation

While this model proposes redirecting CO₂ emissions into high-efficiency biological sinks such as sugarcane plantations, algae bioreactors, and bamboo forests, it currently treats these systems as instantaneous and lossless carbon absorbers. In practice, CO₂ fixation is delayed, variable, and biologically constrained by multiple factors:

Factor	Effect on CO2 Fixation
Growth Phase	Young or senescent plants fix less CO₂ than those in peak growth
Climate & Light	Reduced fixation under low light, drought, or temperature extremes
Nighttime Respiration	Plants respire CO₂ during darkness, offsetting net capture
Soil and Root Limitations	Poor soil health or root stress can limit nutrient uptake, slowing CO₂ assimilation
Algal Turnover	Bioreactor systems require regular harvesting or risk CO₂ saturation and stalling

Implication: Even with optimal gas routing, real-time CO₂ absorption may lag emission, resulting in temporary atmospheric accumulation unless properly buffered.

To address this, future models should include:

• Lag curves to model sequestration delay vs. emission rate.
• CO₂ buffer tanks or regulated injection pacing
• Lifecycle-based tracking of plant biomass carbon retention vs. turnover or decay

Net-zero status should be assessed over longer periods (e.g., weekly/monthly) rather than per-cycle, acknowledging that biological sinks accumulate CO₂ gradually, not on demand.

7.3.2. Lifecycle Carbon Accounting and Sink Saturation Limitations

While this model treats CO₂ redirection as an environmentally positive strategy, a complete life-cycle carbon budget remains outside the current scope. Key variables such as plant respiration, post-harvest biomass decay, transport emissions, and land-use change are not modelled. As such, carbon neutrality is evaluated at the point of CO₂ injection into the biological sink, not the full assimilation and retention pathway. This creates a temporal and spatial mismatch between emission and sequestration, which may be acceptable in short-cycle systems but problematic at industrial scale.

Moreover, all candidate sinks—sugarcane, bamboo, and algae—have upper limits on their CO₂ fixation capacity. Algae bioreactors, for instance, can become saturated within 24–48 hours if biomass is not harvested regularly, while sugarcane plantations exhibit declining CO₂ uptake as atmospheric concentrations plateau near the plant’s physiological maximum (~500 ppm in localized zones). None of these dynamic limitations are currently simulated in this version of the model.

Future implementations should include saturation response curves, storage buffer capacities, and CO₂ pacing algorithms to prevent over-saturation, backflow, or temporal CO₂ buildup. This is particularly important in closed loop or semi-autonomous installations where manual oversight is limited. Additionally, routing infrastructure—assumed here to be lossless and energetically neutral—should be evaluated for compression costs, leak risk, and maintenance overhead in long-distance or industrial deployments.

Clarification on Temporal Lag in CO₂ Sequestration:

While this model assumes that carbon dioxide emissions from microbial metabolism are routed to external biological sinks such as sugarcane plantations or algae bioreactors, it currently treats sequestration as temporally immediate and lossless. In real ecosystems, CO₂ uptake is delayed by biological growth cycles, seasonal variation, and environmental stressors. For example, sugarcane exhibits peak fixation rates only during active vegetative growth, with significantly reduced carbon assimilation during early growth or senescence phases. Likewise, algae systems may saturate with CO₂ under suboptimal light or nutrient depletion, introducing lag between gas delivery and biomass conversion.

Therefore, while the model achieves theoretical carbon neutrality on a per-cycle basis, actual atmospheric offsetting will occur over extended timescales—days to weeks for algae systems, and weeks to months for land-based crops. This temporal mismatch may result in transient CO₂ accumulation unless buffered via intermediate gas storage or throttled injection rates. Future models should integrate sequestration lag curves and dynamic sink response profiles to represent net-zero trajectories more accurately over realistic operational timelines.

Terminology Note: Distinction Between Biomass Types

In this paper, the term biomass is used in two distinct biological contexts. To avoid confusion, the following terminology is applied throughout:

• Microbial Biomass (also referred to as cellular biomass): Refers to the accumulation of bacterial or phototrophic cell mass during fermentation. This includes protein, lipids, DNA/RNA, and structural polysaccharides formed from substrate carbon.
• Plant Biomass (also referred to as sugarcane biomass or CO₂ sink biomass): Refers to vegetative tissue formed by C₄ plants like sugarcane as part of the system’s carbon fixation loop.

When yield metrics refer to "biomass diversion" or "Biomass Diversion Factor (BDF)", they pertain to microbial biomass only—not plant tissue.

This distinction ensures clarity in carbon accounting, resource budgeting, and life-cycle assessments across both fermentation and sequestration components.

8. Theoretical Results and Yield Modelling

This section consolidates the theoretical stoichiometric outcomes of the dual-stage microbial hydrogen production system. Calculations are based on idealized assumptions of complete substrate conversion, lossless gas extraction, and unrestricted enzyme performance.

• Hydrogen quantities are provided in both moles (mol) and grams (g). The molar mass of hydrogen is taken as 2.02 g/mol.
• Carbon dioxide (CO₂) outputs are reported in grams (g), based on a molar mass of 44.01 g/mol.
• Gas volumes are calculated at Standard Temperature and Pressure (STP), defined in this model as 0°C and 1 atm, with a molar volume of 22.4 L/mol. Deviations due to temperature and pressure are not modelled in this version.
• Cost values are reported in U.S. dollars (USD) and are based on industrial utility rates and bulk glucose pricing within the United States.

All tables and quantitative outputs from Section 8 onward adopt these units unless otherwise specified. This standardization supports comparative interpretation, external benchmarking, and ease of conversion across related studies.

Metric	Definition	Example
Mass Yield (%)	Percentage of hydrogen mass relative to input substrate mass.	24.24 g H₂ from 180.16 g glucose = 13.45%
Energy Yield (kJ / MJ)	Total energy content of hydrogen produced based on 286 kJ/mol H₂.	12 mol H₂ = 3432 kJ = 3.43 MJ
Energy Efficiency (%)	Ratio of energy output to total system input. Two forms are used: NEF and EROIE.	NEF = 78.1%; EROIE = 3.58
Yield by Input Basis	Reference used for input mass: glucose-only, glucose + water, or butyrate + water.	Mass yield = 13.45% (glucose), 8.41% (glucose + water)
Volumetric Productivity	Hydrogen volume produced per litre of reactor volume per day (L/L/day).	44.85 L/L/day
Areal Productivity	Hydrogen volume produced per square meter of light-exposed surface area per day (L/m²/day).	~1120 L/m²/day
Carbon-to-CO₂ Ratio	Total CO₂ output mass relative to input glucose mass, reflecting oxidation extent.	264.06 g CO₂ from 180.16 g glucose = 146.56%
Net Energy Gain (kJ)	Energy output minus system input, accounting for ideal or real-world losses.	2472 kJ net gain at theoretical max

• Clarification on Energy Efficiency Metrics:

• Net Energy Fraction (NEF) = Output Energy / (Output + Input)
• Energy Return on Input Energy (EROIE) = Output Energy / InputThese definitions provide complementary perspectives on system performance, reflecting both absolute and relative energy returns.

8.1. Yield Summary

Starting from 1 mol glucose, the system yields 12 mol H₂ (2 from dark fermentation and 10 from photofermentation), with 6 mol water consumed. This equals 13.45% mass yield (glucose only) or 8.41% (glucose + water).

Stage	H2 Produced	Yield by Mass (Input Basis)
Dark Fermentation	2 mol	2.24% (per mol glucose)
Photofermentation	10 mol	10.29% (per mol butyrate + water)
Total (per cycle)	12 mol	13.45% (glucose only) / 8.41% (glucose + water)

Substrate Contribution to Energy Yield (Not Mass-Based)

Input Component	Input Mass (g)	Role	Energy Contribution
Glucose	180.16	Primary fermentable substrate; donor of electrons and carbon	~100% of caloric energy input
Water	108.09	Electron donor during photofermentation; contributes hydrogen atoms but no caloric energy	~0% (non-energetic, non-combustible)

Clarification on Water’s Role:

Water contributes hydrogen atoms, not caloric energy. During photofermentation, water serves as an electron donor whose hydrogen atoms are split via light-driven nitrogenase activity. However, this process does not provide chemical energy—only electrons and protons. Therefore, the system’s entire energy yield originates from glucose oxidation, even though the total hydrogen yield includes atoms from both glucose and water.

Side Note: While water provides hydrogen atoms for photofermentation via light-driven nitrogenase activity, it does not contribute any usable energy in the form of enthalpy. The entire system’s energy output derives from glucose oxidation alone. Any reference to “yield including water” reflects mass-based input comparisons, not energy sources.

While the total mass input includes both glucose and water, their contributions to energy output are not equivalent. Glucose serves as the primary energy substrate, whereas water functions as a hydrogen donor with negligible caloric value but serving as a critical electron donor in photofermentation. This table clarifies the input components' roles in the system’s energy yield.

Clarification on Yield Metrics:

While mass yield is a convenient metric for comparing biochemical conversion efficiency, hydrogen is primarily valued as an energy carrier. Therefore, energy-normalized yields (e.g., MJ per mol substrate) provide a more meaningful basis for system benchmarking and economic comparison. Mass-based values are still included to enable direct comparison with microbial and enzymatic yield literature, where this format remains common.

• Molar mass of H₂: 2.02 g/mol
• Total H₂ mass: 12 mol × 2.02 g/mol = 24.24 g
• Total input mass: 180.16 g (glucose) + 108.09 g (water) = 288.25 g

Yield Basis Clarification: Hydrogen yield can be expressed using multiple mass-based references, depending on whether input mass includes only the primary substrate or all reactants. This model provides three key perspectives for clarity:

Yield Basis	Input Mass (g)	H2 Produced (g)	Mass Yield (%)	Notes
Glucose only	180.16	24.24	13.45	Traditional metric in microbial literature
Glucose + Water	180.16 + 108.09 = 288.25	24.24	8.41	Reflects total input mass for full-cycle fermentation
Butyric Acid + Water	88.11 + 108.09 = 196.20	20.20	10.29	Applies to photofermentation stage only

Each format is useful depending on whether the analysis focuses on:

• Substrate-to-product conversion (glucose-only)
• Total resource consumption (glucose + water)
• Stage-specific yield (butyrate + water)

These yield values are based on complete conversion assumptions and are recalculated in adjusted form later using BDF and BRF values in Section 8.3.2.

8.1.1. H₂ Yield Sensitivity

To illustrate system variability under suboptimal conditions, this table presents a range of hydrogen yields and their associated energy and efficiency impacts. These values help define a realistic operating envelope for future system iterations.

Note on Efficiency Units:

All efficiency values in the following table are calculated based on a fixed input energy of 960 kJ per cycle, reflecting the theoretical system energy requirement under ideal conditions. These values represent apparent system efficiency, assuming consistent operation time and energy expenditure per batch.

In practical systems, energy input may scale dynamically with process duration, hydrogen yield, or system throttling. For this reason, a separate table modelling scaled input efficiency is provided in Section 8.3.2. This distinction ensures clarity between static efficiency benchmarking and adaptive, yield-dependent energy analysis.

H2 Yield (mol)	H2 Energy Output (kJ)	Net Energy Gain (kJ)	Apparent System Efficiency (%)	H2 Yield by Mass (%)
12	3432 (3.43 MJ)	2472	78.14	13.45
10	2860	1900	74.87	11.21
8	2288	1328	70.44	8.97
6	1716	756	64.13	6.73

• Efficiency Under Scaled Input Conditions (Projected)

This table models dynamic input scaling based on hydrogen output, assuming the system runs for shorter durations or uses less power when yield is reduced. It shows that efficiency does not drop sharply if energy input is adjusted proportionally.

H2 Yield Sensitivity	H2 Energy Output	Scaled Input Energy (kJ)	Net Energy Gain (kJ)	Scaled System Efficiency (%)
12	3432	960	2472	78.1
10	2860	800	2060	72.0
8	2288	640	1648	71.8
6	1716	480	1236	71.6

8.1.2. CO₂ Output Analysis

• CO₂ Produced: 6 mol per mol glucose.
• Molar mass of CO₂: 44.01 g/mol
• Total CO₂ mass: 6 × 44.01 = 264.06 g
• CO₂-to-glucose mass ratio: (264.06 / 180.16) × 100 ≈ 146.56%

Side Note: For earlier discussion of oxygen contribution from water molecules, see Section 3.4.

8.2. Modelling Assumptions

• Full conversion of substrate to desired products
• Ideal anaerobic conditions, optimal temperature (37°C), and sufficient light intensity
• No gas loss, contamination, or membrane inefficiency
• Biofilm performance is spatially uniform.

These theoretical outcomes represent upper-bound estimates. Real-world yields may vary due to substrate limitations, microbial kinetics, enzyme inhibition, or incomplete gas capture, all of which are excluded from this model for simplicity and clarity.

8.3. Sensitivity Considerations and Yield Variability

While this model reports a theoretical hydrogen yield of 12 mol H₂ per glucose cycle under ideal conditions, real-world systems rarely achieve perfect conversion. Key process variables—including substrate concentration, light intensity, microbial activity, and gas separation efficiency—can introduce significant variability.

For instance, a 20–30% reduction in photofermentation efficiency (e.g., due to suboptimal lighting or VFA buildup) would lower total hydrogen yield to ~9–10 mol H₂ per cycle. Similarly, incomplete CO₂ separation or membrane fouling could reduce effective gas collection rates, impacting the net energy recovery and skewing cost estimates.

In energy terms, a drop from 12 mol H₂ to 9 mol would reduce gross energy output from ~3.43 MJ to ~2.57 MJ—cutting net energy gain by nearly 25%.

These outcomes highlight the importance of system robustness, adaptive control mechanisms, and experimental optimization. While this paper establishes a theoretical upper-bound, future models should incorporate probabilistic ranges, confidence intervals, or empirical yield curves to better represent operational variability.

8.3.1. Yield Realism Gradient

To contextualize the theoretical yield of 12 mol H₂ per mole of glucose, a performance gradient is presented below. This stratification illustrates hydrogen yield performance under progressively realistic constraints, ranging from idealized models to conservative real-world outcomes.

Scenario	H2 Yield (mol)	Description
Theoretical Ceiling	12	Assumes complete glucose and VFA conversion, no biomass generation, no energy loss, and perfect enzyme expression.
Optimized System	10	Partial VFA conversion (~85–90%), minimal diversion into biomass, high metabolic efficiency.
Practical High-End	8	Includes biomass formation (10–15%), moderate kinetic losses, and less efficient light or nutrient use.
Conservative Estimate	6	Reflects experimental averages with significant VFA retention, non-optimized strains, and real-world kinetic bottlenecks.

8.3.2. Biomass Corrected Yield Estimates

Scenario	BDF (%)	BRF (%)	Effective H2 Yield (mol)	Energy Output (kJ)	Efficiency (%)	Yield by mass (%)
Theoretical Max	0	1.00	12	3432	78.14	13.45
Moderate Biomass	10	0.90	10.8	3080	74.2	12.11
Realistic Field	20	0.80	9.6	2746	70.0	10.76
Conservative	30	0.75	9.0	2574	67.3	10.08

These estimates assume proportional diversion of glucose- and VFA-derived carbon into non-gas biomass structures. This includes amino acids, nucleotides, and structural polysaccharides. Although the degree of biomass formation depends on nutrient conditions, reactor cycling, and microbial growth phase, empirical studies support a diversion range of 10–30% in batch systems [3].

Side Note: Adjusted hydrogen yields in this table are derived by estimating carbon diversion from glucose into microbial biomass. As biomass retains significant substrate carbon, this diversion directly reduces the number of carbon atoms oxidized to CO₂ during photofermentation. Since hydrogen production in R. sphaeroides is stoichiometrically linked to CO₂ evolution, yield losses are approximated proportionally based on carbon retention assumptions. Further refinement using empirical carbon flux analysis is planned in future models.

8.3.2.1. Carbon Diversion and Impact on H₂ Yield

In the dual-stage model, hydrogen yield is tightly coupled to the oxidation of glucose-derived carbon atoms. Each mole of glucose (C₆H₁₂O₆) contains six carbon atoms. Under ideal, complete oxidation:

• 2 carbon atoms are released as CO₂ during dark fermentation.
• The remaining 4 carbon atoms are released during photofermentation of butyric acid.

Because photofermentative hydrogen production by Rhodobacter sphaeroides is stoichiometrically linked to CO₂ evolution—with an approximate 2:1 molar ratio of H₂ to CO₂—any carbon retained in biomass instead of being oxidized results in a proportional loss of H₂ yield.

Hydrogen Yield Equation (BDF-Corrected):

To account for the diversion of substrate carbon into microbial biomass, the following formula calculates the total hydrogen yield per mole of glucose based on the Biomass Diversion Factor (BDF).

Let x represent the fraction of carbon atoms from glucose that are diverted to biomass. This value ranges from 0 (no diversion) to 0.5 (half of all carbon lost to biomass).

In the system:

• Dark fermentation consistently yields 2 mol of hydrogen, unaffected by biomass diversion.
• Photofermentation ideally yields 10 mol of hydrogen, assuming all 4 remaining carbon atoms from glucose (after dark fermentation) are used to form butyric acid and are fully oxidized.

However, when biomass diversion occurs, fewer carbon atoms are available for oxidation in the photofermentation stage. Since 6 carbon atoms exist per glucose molecule and 2 are always released in dark fermentation, the remaining 4 are subject to diversion.

The number of available carbon atoms for photofermentation becomes: 4 – (6 × x)

This means the actual photofermentation yield is: 10 × (available carbon ÷ 4) = 10 × (1 – 1.5x)

Total hydrogen yield, combining both stages, is then: H₂ yield = 2 + 10 × (1 – 1.5x)

Which simplifies to: H₂ yield = 12 – 15x

Final equation:

Adjusted hydrogen yield = 12 – 15 × BDF

This equation assumes dark fermentation is constant, and photofermentation yield declines linearly with increasing biomass diversion.

8.3.2.2. Biological Realism Factor (BRF) Adjustment

These values can be combined with the Biomass Diversion Factor (BDF) for layered correction in scenario modelling. While BRF reflects metabolic inefficiencies even in optimized strains, BDF represents substrate diversion into biomass structures. Together, they define a more biologically constrained performance envelope.

8.4. Energy Yield Estimate

Each mole of hydrogen releases ~286 kJ of energy upon combustion. Based on 12 mol H₂ per glucose cycle:

• Energy yield per mol glucose: 12 × 286 = 3432 kJ, or ~3.43 MJ

This figure assumes complete combustion of hydrogen and no energy loss during storage or transfer.

9. Net Energy Balance and System Efficiency

To evaluate the practical viability of the proposed dual-stage microbial hydrogen production system, it is essential to quantify not only the theoretical energy output but also the energy required to operate the system components. This section outlines a preliminary energy balance based on stoichiometric hydrogen yield and estimated energy inputs under idealized bench-scale conditions.

9.1. Hydrogen Energy Output

Under ideal conditions, the system produces 12 mol of H₂ per mole of glucose, resulting in a gross energy output of:

12

>12 mol × 286 kJ/mol = 3432 kJ (or ~3.43 MJ)

To reflect biological inefficiencies in enzyme expression and metabolic function, a Biological Realism Factor (BRF) is applied. BRF-scaled energy values are summarized in Table 8.3.2, alongside yield and efficiency adjustments. These corrected values define the system’s expected energy ceiling under non-ideal operating conditions.

9.2. Estimated Energy Inputs

Energy inputs were estimated based on typical laboratory-scale setups and published literature. All values assume one full processing cycle (~24 hours) for 1 mole of glucose.

Component	Estimated Input (kJ)	Assumptions
LED Illumination (For Photofermentation)	600 kJ	~5–10 Watts LEDs over 24 hours
Vacuum-Assisted Gas Extraction	216 kJ	~30 Watts pump for 2 hours
Stirring, Control and pH regulation	144 kJ	6 Watts average load for 6 hours
Total Input Energy	960 kJ

9.2.1. Net Energy Gain and Efficiency

Under theoretical conditions, the system produces 3432 kJ of energy from hydrogen combustion and requires 960 kJ of energy to operate per cycle.

• Net Energy Gain

$$NetGain=3432kJ\left(ouput\right)-960kJ\left(input\right)=2472kJ$$

To quantify system performance, two different efficiency metrics are introduced:

1.

>1. Energy Return on Input Energy (EROIE): This reflects the sum of energy is recovered per unit of input energy.

$$EROIE={\displaystyle \frac{OutputEnergy}{InputEnergy}}={\displaystyle \frac{3432}{960}}\approx 3.575$$

2.

>2. Net Energy Fraction (NEF): This represents the proportion of the total energy cycle that results in useful output.

$$NEF={\displaystyle \frac{OuputEnergy}{Ouput+InputEnergy}}={\displaystyle \frac{3432}{3432+960}}\approx 0.7814(or78.14\%)$$

Clarification: The term "system efficiency" is often misused in biological models. Here, both EROIE and NEF are reported to clarify different performance aspects:

∙

EROIE reflects viability from an energy investment perspective.

∙

NEF reflects how much of the total system energy remains in usable form.

9.3. Real-World Expectations

While the theoretical efficiency approaches 80%, practical efficiencies are expected to be lower. A biomass yield penalty of 5–10% is typically observed in practical systems due to resource allocation toward microbial maintenance and division. This was excluded from the current energy model to isolate system performance potential. Based on operational inefficiencies, microbial kinetics, and engineering limitations, realistic system efficiency may range between 40–60%. Factors influencing actual performance include:

• Partial pressure and selectivity losses in membrane separation
• Biofilm detachment or overgrowth
• Incomplete substrate conversion
• Heat loss and energy dissipation in control systems.
• LED degradation and photon scattering inefficiencies.
• These must be addressed in future experimental trials and lifecycle assessments.

9.3.1. Summary and Comparative Outlook

Metric	Value
H2 Output (12 mol)	3432 kJ (3.43 MJ)
Estimated Input Energy	960 kJ
Net Energy Gain	2472 kJ
Idealized System Efficiency	~78.1%
Expected Real Efficiency	40–60% (projected)

This energy model provides a foundational benchmark for the future development of a scalable, energy-positive, microbial hydrogen platform. Further refinements will require integration with kinetic modelling, pilot-scale data, and dynamic energy budgeting across varying reactor designs.

9.3.2. Clarifying Efficiency Metrics: Separation Loss vs. Total System Efficiency

The paper currently references two distinct ranges for practical energy efficiency:

• 67–72% efficiency — reflecting only gas separation penalties and membrane inefficiencies under ideal microbial performance.
• 40–60% efficiency — reflecting full-system inefficiencies including microbial losses, biomass diversion, expression instability, and light/mass transfer limitations.

Clarification: These ranges are not contradictory but represent different layers of system constraint.

To resolve potential confusion, the following delineation is made:

Efficiency Scope	Description	Estimated Range
Gas-Adjusted Efficiency	Accounts for membrane inefficiencies and gas separation penalties (e.g., 8% loss from 92% H₂ recovery)	67–72%
Biology-Adjusted Efficiency	Includes gas losses plus microbial inefficiencies such as ATP shortfalls, cofactor scarcity, biomass diversion, and yield throttling	40–60%
Theoretical Max (NEF)	Net energy fraction under idealized operation	~78.1%

The higher value (67–72%) assumes microbes operate near their theoretical ceiling, and the only losses occur in gas handling.

The lower value (40–60%) reflects real-world biological constraints such as:

• BRF-adjusted expression losses
• Biomass carbon diversion (BDF)
• ATP or cofactor bottlenecks
• Incomplete substrate conversion
• Conclusion: When quoting system efficiency, it is essential to specify which domain the losses originate from. Future empirical models should separately quantify:
• Engineering penalties (e.g., membrane loss, pump inefficiency)
• Biological penalties (e.g., yield suppression, carbon retention) to provide a more precise picture of overall energy performance.

Gas Separation Penalty

Real-world membrane systems rarely achieve perfect separation. With a hydrogen recovery efficiency of ~92% and vacuum pump efficiency of ~90%, two energy penalties must be considered:

• Hydrogen Loss (~275 kJ): Due to unrecovered hydrogen molecules during membrane passage, reducing usable gas output.
• Vacuum Operation Cost (~216 kJ): Represents the power required to run the vacuum pump, even under efficient cycling.

Together, these factors reduce the theoretical energy output of 3432 kJ per cycle by a total of ~491 kJ, resulting in an adjusted net energy gain of ~1960–2200 kJ, depending on membrane and vacuum system efficiency. This corresponds to a practical energy efficiency of 67–72%, compared to the theoretical maximum of 78.1%.

These adjustments reinforce the importance of staged membrane arrays, fouling resistance, and supplemental polishing technologies (e.g., PSA or desiccant traps) for maintaining fuel-grade output.

• Clarified Statement on Gas Separation Losses and Energy Penalty

In the original model, the energy penalty from gas separation was expressed as a drop in usable hydrogen energy—from 3432 kJ (theoretical) to ~3157 kJ (realistic), amounting to an 8% loss. This value includes both the loss of hydrogen due to imperfect membrane selectivity (typically 5–10% unrecovered gas) and the energy required to operate the vacuum extraction system (216 kJ per cycle based on 30W pump usage for 2 hours).

It is important to distinguish between these two contributors:

• Gas Loss Penalty (~274 kJ): Represents the portion of hydrogen that diffuses inefficiently or is unrecovered during membrane separation.
• Vacuum Energy Cost (216 kJ): Represents the energetic overhead of powering the extraction system.
• Clarification on Energy Penalty Attribution:

The total energy penalty from the gas separation system is composed of two distinct components that must be treated separately:

• Hydrogen Loss Penalty (~275 kJ): This represents the reduction in usable energy output due to membrane inefficiency—i.e., hydrogen gas that escapes capture and is therefore unavailable for combustion.
• Vacuum Pump Energy Input (~216 kJ): This is an operational energy cost, representing the power required to run the membrane extraction system during gas separation.

These values should not be combined into a single energy loss from the output stream, as they affect different sides of the energy balance equation:

• The hydrogen loss reduces gross energy output:3432 kJ → 3157 kJ usable H₂ energy.
• The vacuum pump cost adds to system input energy:Total input = 960 kJ (including 216 kJ from vacuum system)

Therefore, the net energy gain remains:

• Net energy = 3157 − 960 = 2197 kJ
• Real-world system efficiency ≈ 64% (2197 / 3432)

This refined breakdown avoids double-counting and accurately distinguishes between energy losses due to system inefficiency and energy costs due to system operation.

Side Note: Yield values assuming 10 mol H₂ from photofermentation reflect a best-case ATP supply scenario. Under energy-limited conditions, ATP throttling may reduce this to 7.5–9 mol, resulting in cycle yields of ~9.5–11 mol H₂ and a corresponding drop in energy output and system EROIE.

10. Production Rate and Daily Throughput

To contextualize the system's practicality and operational output, this section presents a normalized daily breakdown of production, energy output, and carbon emissions based on a complete reactor cycle duration of 72 hours (3 days).

10.1. Cycle Duration Justification

The proposed system operates through a sequential two-stage microbial process:

10.1.1. Stage 1: Dark Fermentation (Clostridium butyricum)

• Converts glucose into hydrogen and volatile fatty acids (VFAs)
• Doubling time: ~3–6 hours
• Hydrogen production typically peaks at 18–24 hours post-inoculation.
• Full substrate conversion is generally achieved within 24 hours in batch mode.

10.1.2. Stage 2: Photofermentation (Rhodobacter sphaeroides)

• Converts VFAs into additional hydrogen via nitrogenase-driven metabolism.
• Process is light-dependent and slower.
• Photofermentation typically requires 2–5 days depending on light intensity and VFA concentration.
• Butyric acid breakdown extends total time requirements.

Based on this information, the full operating cycle is reasonably assumed to require approximately 72 hours (3 days).

10.2. Daily Breakdown of Core Outputs

Parameter	Per Cycle (3 Days)	Per Day (Averaged)
Hydrogen Output	12 mol	4 mol
Energy Output	3432 kJ	1144 kJ
Hydrogen Volume (STP)	~269 L	~89.7 L
CO₂ Emitted	6 mol	2 mol
CO₂ Mass	264 g	88 g
Glucose Consumed	180.16 g	60.05 g
Operational Cost	$0.637	~$0.212

Side Note: Hydrogen volume assumes STP = 22.4 L/mol at 0°C and 1 atm. Real values may vary depending on operational temperature and pressure.

10.3. Hydrogen Productivity Rate (Volume-Based)

Assuming a 2-litre total working reactor volume:

• Hydrogen output/day = 89.7 L / 2 L = 44.85 L H₂ per litre of reactor per day
• Energy output/day = 1144 kJ / 2 L = 572 kJ per litre per day

These values represent an idealized productivity scenario under fully optimized and stable microbial activity. Real-world rates may vary based on inoculum density, substrate purity, light intensity, and ambient conditions.

Note on Gas Volume Calculations:All hydrogen volume values in this paper are calculated assuming Standard Temperature and Pressure (STP) defined as:

• 0°C (273.15 K)
• 1 atm pressure (101.325 kPa)
• Molar volume of 22.4 L/mol

If operating conditions differ (e.g., 25°C or pressurized storage), gas volume may vary accordingly. For precision applications, corrected volumes using the Ideal Gas Law or real gas equations should be applied.

10.4. Microbial Biomass Yield Efficiency

Literature on Clostridium butyricum batch fermentation systems commonly report dry biomass concentrations ranging from 1.0 to 3.0 g/L, depending on substrate load, nutrient availability, and anaerobic conditions. For this model, a conservative estimate of 2.0 g/L was assumed, consistent with OD₆₀₀ values of ~1.5–2.0 reported in microbial hydrogen production studies.

Assuming a working bioreactor volume of 2 litres—representative of a lab-scale system operating at full aqueous capacity—the total microbial dry biomass is estimated at 4.0 grams. Under these conditions, the reactor yields:

• mol H₂ per gram of biomass per day
• 0.5 mol CO₂ per gram of biomass per day

This productivity aligns with the theoretical ceiling. Literature reports 0.15–0.75 mol/g/day depending on strain and light conditions. This model assumes optimized conditions and gene stack expression without metabolic burden.

Side Note: This system-wide production rate of ~4 mol H₂/day corresponds to a microbial biomass concentration of 2 g/L across 2 L of working volume. At this biomass level, approximately 5–15% of substrate carbon is redirected toward cellular material. This supports the application of a Biomass Diversion Factor in overall yield and energy calculations (see Section 8.2.2). Section 9.4 provides the biomass-normalized rates derived from this same baseline (1.0 mol H₂/g biomass/day and 0.5 mol CO₂/g biomass/day), ensuring direct alignment between overall gas output and microbial load.

10.5. Implications for Scale-Up

This daily throughput model allows for direct scalability calculations. For example:

• A 100-litre reactor system (assuming same conditions) would yield:

  ○

  ~8.97 m³ H₂/day

  ○

  ~114.4 MJ/day

  ○

  ~8.8 kg CO₂/day

Such normalized values are essential for evaluating deployment potential across decentralized grids, small-scale industrial applications, or off-grid energy systems with integrated carbon management strategies.

11. CO₂ Redirection to Passive Sugarcane Sink

To ensure net-zero emissions and maintain the system’s carbon-neutral profile, the carbon dioxide produced by microbial metabolism is redirected into a dedicated sugarcane greenhouse or plantation pod adjacent to the reactor setup.

11.1. Daily CO₂ Emissions

Based on stoichiometric fermentation pathways, the system emits:

• 6 mol CO₂ per cycle = 264 g CO₂ per 3 days
• Average per day = 2 mol CO₂ ≈ 88 grams/day

11.2. Sugarcane CO₂ Fixation Capacity

Sugarcane is a high-efficiency C₄ photosynthetic crop known for rapid carbon uptake.

• Literature indicates sugarcane can fix ~666 g CO₂ per kg of biomass produced.
• Required CO₂ offset: 88 g/day.

Required sugarcane growth per day:

$${\displaystyle \frac{88g}{666g/kg}}\approx 0.132kg=132gsugarcane/day$$

This requirement is modest and achievable in a small indoor greenhouse pod or planter-scale setup, especially under controlled light and nutrient conditions. Even minimal vegetative growth is sufficient to fully offset the reactor’s daily CO₂ emissions.

11.2.1. Why Sugarcane?

Sugarcane offers two-fold advantages:

• Exceptional CO₂ fixation due to its C₄ metabolic pathway
• Glucose-rich stalks make it a compatible future feedstock for microbial fermentation.

This symbiotic potential means that the same crop used for carbon offset could eventually be harvested to support future reactor cycles, contributing to a semi-closed bioresource loop.

11.3. Temporal Dynamics of CO₂ Fixation

While the daily CO₂ emission from microbial fermentation is estimated at ~88 g/day, the corresponding fixation by sugarcane biomass (~132 g/day growth for offset) assumes a smoothed average over an extended photoperiod and seasonal cycle. Sugarcane growth is not linear or instantaneous; it follows diurnal patterns and growth-phase-dependent rates influenced by light exposure, nutrient levels, and ambient conditions.

Therefore, the CO₂ fixation rate cited here represents an averaged sequestration benchmark, not a real-time sink. In practice, temporal mismatches between gas production and biological absorption may occur—particularly during early germination, cloudy periods, or seasonal dormancy.

To mitigate these mismatches, future designs could incorporate:

• Buffer tanks to store CO₂ until peak uptake hours.
• Multi-species sink systems (e.g., algae + sugarcane) for temporal balancing.
• Monitoring and throttle protocols based on plant health and CO₂ levels.

This correction ensures the system’s carbon neutrality is evaluated on a rolling timescale, not instantaneously. Accurate life-cycle sequestration tracking is essential for long-term deployment and ecological verification.

Author’s Note: The Industrial Concern

While this model proposes a sustainable and ecologically integrated CO₂ management strategy, it is important to acknowledge a concern: industrial-scale implementations may not uphold the passive offset system. For economic or spatial reasons, CO₂ redirection could be bypassed or ignored, undermining the environmental integrity of the design.

Therefore, responsible implementation must treat CO₂ fixation as a non-optional component, not merely a sustainability add-on. Without it, the system reverts to a net-positive carbon emitter, contradicting its intended purpose.

12. CO₂ Sink Implementation and Maintenance Considerations

The efficacy of the proposed microbial hydrogen production system relies not only on hydrogen output but also on the successful redirection and neutralization of byproduct CO₂. This section addresses the physical deployment, delivery method, and oversight mechanisms associated with the sugarcane-based carbon sink.

12.1. Physical Setup: Sugarcane Pod Location

The sugarcane plants are assumed to be grown within a dedicated greenhouse pod or adjacent indoor soil-based chamber. This setting allows for year-round growth under controlled temperature, humidity, and light conditions. Soil-based planting is prioritized over hydroponics to simplify nutrient cycling and reduce technological overhead. The pod is spatially integrated to minimize gas transfer distance.

12.2. CO₂ Transfer Mechanism

CO₂ generated during fermentation is captured via membrane separation and routed into the sugarcane pod through low-pressure, gas-tight tubing. Dispersion strategies include:

• Atmospheric diffusion: Gas is released into the plant chamber for natural absorption.
• Root-zone injection: CO₂ is directed into soil channels for partial dissolution and uptake.

For bench-scale systems, passive atmospheric mixing is likely sufficient.

12.3. Maintenance and Oversight

To prevent CO₂ buildup due to plant failure or neglect, future-scale systems may include:

• CO₂ sensors in the exhaust loop or pod
• Automatic throttling or shutdown if CO₂ thresholds are exceeded.
• Optional health monitoring of sugarcane growth using basic spectral sensors.

These features ensure carbon sequestration remains enforced, not optional.

12.4. Ethical and Industrial Caution

While the passive CO₂ sink model is designed for simplicity and reliability, there is a tangible risk that industrial-scale adopters may attempt to omit or disable the sequestration component to reduce costs or footprint. Such practices would fundamentally undermine the system’s sustainability claims and revert it to a net-emitter.

To avoid this outcome, system deployment should be guided by strict sustainability protocols that include the carbon sink as a non-optional element of the hydrogen production framework. Long-term success hinges not only on microbial efficiency, but also on upholding the environmental intent of the design.

13. Kinetic Considerations and Theoretical Modelling Framework

13.1. Objective of Kinetics Modelling

While stoichiometric analysis provides a snapshot of theoretical yield, it lacks the temporal resolution required to optimize real-world reactor performance. Kinetic modelling addresses this gap by characterizing the rate at which biochemical transformations occur under varying environmental and operational conditions.

In the context of this dual-stage microbial hydrogen system, kinetics modelling serves the following key purposes:

• Predicting hydrogen production rates over time for each microbial stage
• Identifying limiting factors, such as substrate saturation, product inhibition, or light flux variability
• Optimizing substrate feeding strategies to avoid bottlenecks or microbial shock
• Enabling scale-up modelling through time-dependent simulations of batch and continuous systems
• Supporting future automation by providing a framework for real-time feedback control

By introducing kinetic parameters into the system framework, the model transitions from a static yield estimate to a dynamic, time-resolved process capable of iterative refinement and control.

13.3. Dark Fermentation Kinetics (Clostridium butyricum)

The dark fermentation stage, facilitated by Clostridium butyricum, is responsible for the initial conversion of glucose into hydrogen and volatile fatty acids (VFAs). This process can be described using Monod kinetics, which relate microbial growth or substrate consumption rate to substrate concentration:

$$r=r_{max}{\displaystyle \frac{\left[S\right]}{K_s+\left[S\right]}}$$

Where:

• $r$: Rate of glucose consumption or hydrogen production
• $r_{max}$: Maximum rate under saturated conditions
• $\left[S\right]$: Substrate concentration (glucose)
• $K_s$: Half-Saturation constant

In specific conditions, substrate inhibition may occur at high glucose concentrations. This can be modelled by extended forms such as the Haldane equation to reflect non-linear inhibition kinetics.

13.3.1. Example Calculation using Empirical Parameters

Based on reported kinetic data for Clostridium butyricum, the following parameters were selected for modelling hydrogen production via Monod kinetics [Khanal et al., 2018]:

• $r_{max}=14.92{day}^{-1}$ (max growth rate)
• $K_s=1.02gCOD/L$ (substrate affinity)
• $\left[S\right]=1.5gCOD/L$ (Glucose feed concentration assumed

Substituting into Monod equation:

$$r=14.92\times {\displaystyle \frac{1.5}{1.02+1.5}}=14.92\times 0.595=8.87{day}^{-1}$$

Thus, under these conditions, the system achieves a glucose-specific growth rate of 8.87 day⁻¹, indicating strong substrate utilization.

[Source:] Khanal, S. K., et al. (2018). Kinetic modelling and microbial community analysis for high-rate hydrogen production in a dynamic membrane bioreactor. PubMed

13.4. Photofermentation Kinetics (Rhodobacter sphaeroides)

Due to the lack of real-time experimental data for photofermentation kinetics, light-saturation modelling is presented using standard literature values for similar strains under ideal lab conditions. These values are included solely for theoretical demonstration.

The second stage involves Rhodobacter sphaeroides, which uses VFAs produced during dark fermentation to generate additional hydrogen via photofermentation. This process is primarily limited by light intensity and substrate availability. The rate of hydrogen production can be approximated using light-dependent kinetics:

$$r=r_{max}{\displaystyle \frac{I}{K_I+I}}$$

Where:

• $I$: Incident light intensity (µmol photons/m²/s)
• $K_I$: Light saturation constant

In addition to light, nitrogenase enzyme activity and VFA accumulation influence the overall rate. Optimization may require tuning light exposure profiles and avoiding product buildup that inhibits microbial metabolism.

13.4.1. Example Calculation using Empirical Parameters

Light-dependent kinetics of Rhodobacter sphaeroides were modelled using data from Thaiwong et al. (2015), which report:

• Optimal light intensity: 7500 lux
• $r_{max}=~12.0ml{H}_2/L·h=0.288mol{H}_2/L/day$
• $K_I$ (Estimated light half-saturation constant): ~4500 lux.

Assuming light intensity $I$ = 6000 lux:

$$r=0.288\times {\displaystyle \frac{6000}{4500+6000}}=0.288\times 0.571=0.164mol{H}_2/L/day$$

This rate reflects mid-range illumination conditions achievable in bioreactor systems with LED panels.

[Source:] Thaiwong, N., et al. (2015). Photo-fermentational hydrogen production of Rhodobacter sp. KKU-PS1 under various illumination intensities. SciELO

Side Note: These empirical parameters are context-dependent and vary based on reactor design, inoculum age, and light delivery systems. They are presented here to demonstrate integration of real-world kinetics into the model framework.

13.5. Coupled Kinetic Modelling

To simulate the full dual-stage reactor dynamically, the two microbial stages can be coupled in a sequential or semi-continuous configuration. The output of the first stage (VFAs) becomes the input substrate for the second stage. An ODE-based model can be constructed to track glucose degradation, VFA production, and hydrogen evolution over time:

• Let $G\left(t\right)$: Glucose concentration over time.
• Let $V\left(t\right)$: Volatile fatty acid concentration.
• Let $H\left(t\right)$: Hydrogen yield curve.

Differential equations can describe rates of change in each metabolite pool and be fitted using experimental time-course data in future implementations.

13.6. Limitations and Future Directions

While the proposed kinetic framework outlines a theoretical basis for dynamic system modelling, it currently lacks empirical validation. Parameters such as $K_s$, $K_I$, and $r_{max}$ are not specified due to the absence of real-time growth or yield curve data. Future experimental studies are needed to:

• Calibrate model parameters using actual fermentation trials.
• Validate system response to variable substrate concentrations and light conditions.
• Develop predictive tools for system control and optimization.

These kinetic models are not intended as predictive tools in their current form. Rather, they provide a foundation for rate-dependent system refinement once empirical data becomes available.

Nonetheless, the inclusion of a kinetic modelling framework establishes the groundwork for advancing this hydrogen production system toward real-time performance tuning and intelligent reactor control.

14. Cost Efficiency Estimate and Benchmarking

To assess economic feasibility, an approximate cost breakdown was developed for a single hydrogen production cycle under bench-scale operation, assuming U.S.-based pricing standards [6,7].

14.1. U.S. Cost Assumptions:

• Electricity: $0.10 per kWh (industrial rate)
• Glucose (bulk, biotech-grade): ~$2 per kg
• Membrane unit: $100, amortized over 500 cycles.

14.2. Cost Breakdown per Cycle (1 mol glucose processed):

Item	Cost per Cycle (USD)	Notes
Electricity (0.267 kWh)	$0.027	960 kJ input total
Glucose (180.16 g)	$0.36	Based on $2/kg bulk rate
Membrane amortization	$0.20	$100 module over 500 cycles
Nutrients, buffers, misc.	$0.05	Estimated buffer cost
Total	$0.637

14.3. Cost per Energy Unit:

Based on the total operational cost of $0.637 per cycle, the system’s theoretical cost efficiency can be calculated relative to hydrogen energy output.

At the theoretical maximum (12 mol H₂), the system yields 3.43 MJ of energy:

Cost per MJ H₂ = $0.637 / 3.43 ≈ $0.186/MJ

Cost per mol H₂ = $0.637 / 12 mol ≈ $0.053/mol H₂

To incorporate metabolic inefficiencies, the Biological Realism Factor (BRF) is applied to simulate performance under non-ideal conditions:

BRF	Energy Output (MJ)	H2 Yield (mol)	Cost per MJ ($)	Cost per mol H2 ($)
1.00	3.43	12.0	0.186	0.0531
0.90	3.08	10.8	0.207	0.0590
0.80	2.75	9.6	0.232	0.0663
0.75	2.57	9.0	0.248	0.0708

These figures highlight how cost efficiency degrades as enzyme activity and metabolic performance fall short of ideal. The rise in cost per MJ reflects reduced hydrogen recovery, even though per-cycle expenses remain constant. This reinforces the importance of genetic stability, expression fidelity, and process control in maintaining economically viable operation.

15. Comparative Cost of Hydrogen Production Methods:

To contextualize the cost-effectiveness of the proposed system, the following table compares the theoretical model with conventional and emerging hydrogen production methods based on published estimates.

Method	Cost per MJ H2	Cost per Mole H2	CO2 Emissions	Energy Efficiency	Notes
This Model (Theoretical)	$0.186	$0.053	Medium (recycled)	~78% (ideal)	Bench-scale, bio-based
Steam Methane Reforming (SMR)	$0.80–$1.00	~$0.23–$0.29	High (direct emissions)	~65–75%	Industrial standard
Water Electrolysis (Renewables)	$0.50–$1.50	~$0.14–$0.43	Low (if powered by renewables)	~50–70%	Dependent on clean grid
Biomass Gasification	$0.40–$0.60	~$0.11–$0.17	Medium (feedstock-dependent)	~60–70%	Requires preprocessing
Alkaline Electrolysis (Grid)	~$1.20	~$0.34	High (if fossil-powered)	~55%	Mature but energy-intensive

This comparison highlights the economic promise of biologically driven hydrogen production when paired with low-cost substrates and efficient CO₂ handling. Although these values are theoretical and do not yet include full system capital expenditure or operational overhead, the model suggests a strong potential for competitive deployment in decentralized, sustainable energy frameworks.