Module 5

Hydrothermal Vents & Chemosynthesis

Life powered by Earth's internal heat: black smokers, giant tube worms, and ecosystems at the thermodynamic limit

5.1 Black Smokers: Submarine Volcanic Chimneys

Black smokers are hydrothermal vents that emit superheated fluid at 350–400°C, laden with dissolved minerals including H₂S, Fe²⁺, Mn²⁺, Cu²⁺, and Zn²⁺. When this scorching, acidic fluid meets near-freezing ($\sim 2$°C) ambient seawater, metal sulfides precipitate instantly, creating the dark billowing plumes that give black smokers their name.

The chimney structure grows from the inside out as minerals deposit: an inner core of chalcopyrite (CuFeS₂) and an outer layer of anhydrite (CaSO₄) and sphalerite (ZnS). Growth rates can exceed 30 cm per day, and chimneys can reach heights of 60 m.

Boiling Point Under Pressure

The reason vent fluids can exceed 100°C without boiling is the enormous hydrostatic pressure at depth. The boiling point of pure water as a function of pressure can be approximated as:

$ T_{\text{boil}}(P) \approx 100 + 28.6 \cdot \ln\!\left(\frac{P}{1\;\text{atm}}\right)\;\;°\text{C} $

At a typical vent depth of 2,500 m ($P \approx 250\;\text{atm}$):

$ T_{\text{boil}} \approx 100 + 28.6 \cdot \ln(250) \approx 100 + 28.6 \times 5.52 \approx 258\;°\text{C} $

The dissolved salts in vent fluid raise this further by $\sim 2$°C per wt% NaCl (vent fluids are typically 3–7 wt% NaCl). At extremely shallow vents or where pressures are lower, vent fluids can undergo phase separation, splitting into a low-salinity vapour phase and a high-salinity brine.

Vent Fluid Chemistry

Black smoker fluids are dramatically different from ambient seawater:

pH: 2–3 (vs. 8.1 for seawater) — strongly acidic due to dissolved H₂S and HCl
H₂S: 3–110 mmol/kg (seawater: $\sim 0$)
Fe²⁺: 0.01–24 mmol/kg (seawater: $\sim 0.001\;\mu\text{mol/kg}$)
Mn²⁺: 0.4–3.4 mmol/kg
O₂: 0 (completely anoxic)
Temperature: 350–407°C

Metal Speciation & Chimney Growth

The precipitation kinetics of metal sulfides depend on the mixing ratio of vent fluid with seawater. The rate of chimney growth can be estimated from the precipitation rate:

$ \frac{dh}{dt} = \frac{Q \cdot [\text{Fe}^{2+}] \cdot M_{\text{FeS}_2}}{\rho_{\text{chimney}} \cdot A_{\text{wall}}} $

where $Q$ is the fluid flow rate ($\sim 1\;\text{L/s}$),$[\text{Fe}^{2+}]$ is the iron concentration ($\sim 10\;\text{mmol/L}$),$M_{\text{FeS}_2}$ is the molar mass of pyrite (120 g/mol),$\rho_{\text{chimney}}$ is the chimney density ($\sim 3000\;\text{kg/m}^3$), and $A_{\text{wall}}$ is the chimney wall cross-sectional area.

5.2 Chemosynthesis: Primary Production Without Light

The discovery of hydrothermal vent ecosystems in 1977 by Corliss et al. on the Galapagos Rift fundamentally altered our understanding of life on Earth. Here was an entire ecosystem powered not by photosynthesis but by chemosynthesis —the fixation of CO₂ into organic matter using chemical energy from reduced inorganic compounds.

Sulfur Oxidation: The Primary Energy Source

The dominant chemosynthetic reaction at most hydrothermal vents is the oxidation of hydrogen sulfide:

$ \text{H}_2\text{S} + 2\text{O}_2 \rightarrow \text{SO}_4^{2-} + 2\text{H}^+ \qquad \Delta G^\circ = -798\;\text{kJ/mol} $

This is highly exergonic — almost as much energy as burning glucose ($\Delta G^\circ = -2870\;\text{kJ/mol}$). The energy is captured by chemolithoautotrophic bacteria (e.g., Thiomicrospira, Beggiatoa) and used to fix CO₂ via the Calvin cycle.

Comparison with Photosynthesis

To fix one molecule of glucose via the Calvin cycle requires approximately 18 ATP and 12 NADPH, equivalent to:

$ \Delta G_{\text{glucose}} = +2870\;\text{kJ/mol} $

Photosynthesis provides this energy through 8 photons of red light per CO₂ fixed:

$ E_{\text{photons}} = 8 \times \frac{hc}{\lambda} = 8 \times \frac{6.626 \times 10^{-34} \times 3 \times 10^8}{680 \times 10^{-9}} \times N_A \approx 1{,}410\;\text{kJ/mol} $

Per CO₂ fixed (6 per glucose), photosynthesis delivers $\sim 235\;\text{kJ/mol CO}_2$, while the oxidation of one H₂S provides 798 kJ. Thus, chemosynthesis can in principle fix more carbon per mole of substrate than photosynthesis, though the supply of reduced chemicals is far more limited than sunlight at the surface.

Other Chemosynthetic Reactions

Hydrogen oxidation: $\text{H}_2 + \frac{1}{2}\text{O}_2 \rightarrow \text{H}_2\text{O}$, $\Delta G^\circ = -237\;\text{kJ/mol}$
Iron oxidation: $4\text{Fe}^{2+} + \text{O}_2 + 10\text{H}_2\text{O} \rightarrow 4\text{Fe(OH)}_3 + 8\text{H}^+$, $\Delta G^\circ = -65\;\text{kJ/mol Fe}$
Methane oxidation: $\text{CH}_4 + 2\text{O}_2 \rightarrow \text{CO}_2 + 2\text{H}_2\text{O}$, $\Delta G^\circ = -818\;\text{kJ/mol}$
Ammonia oxidation: $\text{NH}_3 + \frac{3}{2}\text{O}_2 \rightarrow \text{NO}_2^- + \text{H}^+ + \text{H}_2\text{O}$, $\Delta G^\circ = -275\;\text{kJ/mol}$

5.3 Riftia pachyptila: The Giant Tube Worm

Riftia pachyptila is the iconic organism of hydrothermal vent ecosystems. Growing up to 2.4 m in length with growth rates of 85 cm per year (the fastest of any marine invertebrate), these remarkable annelids have no mouth, no gut, and no anus. They rely 100% on endosymbiotic sulfur-oxidizing bacteria of the genus Endoriftia, housed in a specialised organ called the trophosome.

Symbiont Carbon Fixation

The endosymbionts fix carbon using the Calvin–Benson–Bassham cycle, catalysed by RuBisCO. The rate of carbon fixation follows Michaelis–Menten kinetics:

$ A = \frac{V_{\max} \cdot [\text{CO}_2]}{K_m + [\text{CO}_2]} $

where $V_{\max}$ is the maximum fixation rate ($\sim 15\;\mu\text{mol CO}_2\;\text{g}^{-1}\;\text{h}^{-1}$),$[\text{CO}_2]$ is the dissolved CO₂ concentration, and$K_m \approx 50\;\mu\text{mol/L}$ is the Michaelis constant for the Form II RuBisCO used by Endoriftia.

Vent fluid CO₂ concentrations can reach 10–50 mmol/L, far exceeding$K_m$, so fixation operates near $V_{\max}$. The worm supplies H₂S and O₂ to the bacteria via its haemoglobin, which has the remarkable ability to bind both O₂ and H₂S simultaneously at different sites — a feat no mammalian haemoglobin can perform.

Riftia Haemoglobin

Riftia haemoglobin is a giant hexagonal bilayer molecule with molecular weight$\sim 3.5 \times 10^6\;\text{Da}$ (400$\times$ larger than human haemoglobin). It carries free zinc-bound sulfide ions on specific cysteine residues, preventing H₂S from poisoning the mitochondrial cytochrome c oxidase:

$ \text{Hb-Zn}^{2+} + \text{HS}^- \rightleftharpoons \text{Hb-Zn-S}^{-} + \text{H}^+ $

This allows the worm to transport both O₂ (on haem groups) and H₂S (on zinc sites) in the same blood, delivering them independently to the trophosome endosymbionts.

5.4 ANME Archaea & Anaerobic Oxidation of Methane

At cold seeps — locations where methane and other hydrocarbons seep from the seafloor —a remarkable metabolic process occurs: the anaerobic oxidation of methane (AOM). This process is carried out by consortia of anaerobic methanotrophic archaea (ANME) and sulfate-reducing bacteria (SRB).

The AOM Reaction: Life at the Thermodynamic Limit

The net reaction of AOM coupled to sulfate reduction is:

$ \text{CH}_4 + \text{SO}_4^{2-} \rightarrow \text{HCO}_3^{-} + \text{HS}^{-} + \text{H}_2\text{O} $

$ \Delta G^\circ = -16.6\;\text{kJ/mol} $

This is barely exergonic — just enough energy to synthesise approximately 0.3 mol of ATP. For comparison, aerobic methane oxidation yields$-818\;\text{kJ/mol}$. AOM organisms are living at the very edge of thermodynamic feasibility.

The free energy under in situ conditions (accounting for actual concentrations) is:

$ \Delta G = \Delta G^\circ + RT \ln \frac{[\text{HCO}_3^-][\text{HS}^-]}{[\text{CH}_4][\text{SO}_4^{2-}]} $

With typical cold seep concentrations ($[\text{CH}_4] \sim 1\;\text{mM}$,$[\text{SO}_4^{2-}] \sim 28\;\text{mM}$, $[\text{HCO}_3^-] \sim 2\;\text{mM}$,$[\text{HS}^-] \sim 5\;\text{mM}$), the in situ $\Delta G$ranges from $-20$ to $-35\;\text{kJ/mol}$, providing slightly more energy than the standard value.

ANME–SRB Consortia

ANME archaea and SRB grow in intimate physical association, forming structured aggregates visible under fluorescence in situ hybridisation (FISH) microscopy. The archaea oxidise methane and transfer electrons (possibly via direct interspecies electron transfer, DIET) to the sulfate-reducing bacterial partner. Three ANME clades have been identified:

ANME-1: rod-shaped, often without a bacterial partner, uses modified methanogenesis pathway in reverse
ANME-2: cocci, forms shell-type aggregates with Desulfosarcina/Desulfococcus SRB
ANME-3: closely related to Methanococcoides, partners with Desulfobulbus SRB

Doubling times for ANME–SRB consortia are extraordinarily slow: 2–7 months, reflecting the minimal energy available from AOM. This makes them among the slowest-growing organisms known.

5.5 Vent Ecosystem Food Web

Hydrothermal vent communities are among the most biomass-dense ecosystems on Earth, rivalling tropical rainforests in productivity per unit area. The key organisms include:

Primary Producers (Chemosynthetic)

Free-living bacteria: Thiomicrospira, Beggiatoa mats on chimney surfaces
Endosymbionts: sulfur-oxidizers in tube worms, mussels, and clams

Symbiont-Bearing Megafauna

Riftia pachyptila (giant tube worms): up to 2.4 m, growth rate 85 cm/yr
Bathymodiolus (vent mussels): dual symbionts — both sulfur-oxidizers and methanotrophs in same gill cells
Calyptogena magnifica (giant clams): up to 30 cm, sulfur-oxidizing symbionts in gills

Grazers & Predators

Rimicaris exoculata (vent shrimp): swarms of thousands, cultivate episymbiotic bacteria on modified gill covers
Bythograea thermydron (vent crabs): predator/scavenger, chemosensory-guided navigation to vents
Alvinella pompejana (Pompeii worm): the most heat-tolerant animal known, surviving at 80°C (possibly brief exposures to 105°C)

Biomass Density: The Vent Oasis

The biomass contrast between vent sites and the surrounding abyssal plain is staggering:

$ \text{Biomass}_{\text{vent}} \approx 10\text{--}500\;\text{kg/m}^2 $

$ \text{Biomass}_{\text{abyssal}} \approx 0.001\;\text{kg/m}^2 $

This represents a $10{,}000\text{--}500{,}000\times$ enrichment. A single vent field of $\sim 1\;\text{km}^2$ may support more biomass than hundreds of square kilometres of surrounding abyssal plain.

5.6 Vent Fluid Mixing & Chemical Gradients

As hot, acidic, anoxic vent fluid mixes with cold, alkaline, oxygenated seawater, steep chemical gradients form. The mixing ratio $f$ (fraction of vent fluid) determines the local temperature and chemistry:

$ T(f) = f \cdot T_{\text{vent}} + (1 - f) \cdot T_{\text{sw}} $

$ [\text{H}_2\text{S}](f) = f \cdot [\text{H}_2\text{S}]_{\text{vent}} $

$ [\text{O}_2](f) = (1 - f) \cdot [\text{O}_2]_{\text{sw}} $

Life at vents clusters in the narrow zone where both H₂S (electron donor) and O₂(electron acceptor) coexist — typically at mixing ratios $f = 0.001\text{--}0.1$, corresponding to temperatures of 2–40°C. The organisms position themselves precisely in this chemical gradient:

Alvinella: closest to vent, $T \sim 60\text{--}80$°C, $f \sim 0.2$
Riftia: intermediate, $T \sim 15\text{--}35$°C, $f \sim 0.05\text{--}0.1$
Bathymodiolus mussels: periphery, $T \sim 5\text{--}15$°C, $f \sim 0.01\text{--}0.03$
Vent crabs & shrimp: mobile, sampling gradient at various $f$ values

5.7 Vent Biogeography & Larval Dispersal

Hydrothermal vents are ephemeral habitats: individual vents persist for only decades to centuries before the underlying magma source shifts or the chimney collapses. Yet vent species persist over geological timescales. This requires effective larval dispersal between vent fields separated by hundreds to thousands of kilometres.

The dispersal distance of a vent larva depends on the larval duration $\tau_L$and the ambient current speed $u$:

$ L_{\text{disp}} \sim u \cdot \tau_L + \sqrt{2 K_h \tau_L} $

where $K_h \sim 10^2\text{--}10^3\;\text{m}^2/\text{s}$ is the horizontal eddy diffusivity. For Riftia larvae with $\tau_L \sim 38\;\text{days}$and ridge-crest currents of $u \sim 5\;\text{cm/s}$:

$ L_{\text{disp}} \sim 0.05 \times 38 \times 86{,}400 + \sqrt{2 \times 500 \times 38 \times 86{,}400} \approx 164\;\text{km} + 57\;\text{km} \approx 220\;\text{km} $

This is sufficient to connect neighbouring vent fields along the East Pacific Rise (spaced 50–300 km apart) but insufficient to cross large gaps such as transform faults or the vast expanse between the Pacific and Atlantic ridge systems.

Six distinct vent biogeographic provinces have been identified (Rogers et al., 2012):

Northeast Pacific: dominated by Ridgeia tube worms, alvinellid polychaetes
East Pacific Rise: Riftia, Bathymodiolus thermophilus, Alvinella
Western Pacific back-arc basins: Alviniconcha snails, Ifremeria
Indian Ocean: mixture of Pacific and Atlantic taxa, recently discovered
Mid-Atlantic Ridge: Rimicaris shrimp-dominated, no tube worms
Southern Ocean (East Scotia Ridge): unique yeti crabs (Kiwa tyleri)

5.8 Hydrothermal Vents & the Origin of Life

Hydrothermal vents are considered one of the most promising candidate environments for the origin of life on Earth. Martin & Russell (2003) proposed that alkaline hydrothermal vents (white smokers, like the Lost City field) provided the chemical gradients and mineral catalysts necessary for abiogenesis.

The key thermodynamic argument is that the mixing of alkaline vent fluid (pH $\sim 9\text{--}11$) with acidic Hadean ocean water (pH $\sim 5\text{--}6$) creates a proton gradient:

$ \Delta G_{\text{proton}} = -2.303 \cdot RT \cdot \Delta\text{pH} \approx -2.303 \times 8.314 \times 350 \times 5 \approx -33\;\text{kJ/mol} $

This is remarkably similar to the proton motive force used by modern cells to synthesise ATP ($\sim 30.5\;\text{kJ/mol}$). Iron–sulfur minerals (such as mackinawite, FeS) in vent chimneys could have acted as primitive catalysts, predating protein enzymes. The resemblance between Fe–S clusters in modern ferredoxin proteins and the crystal structure of greigite (Fe₃S₄) is often cited as evidence for this hypothesis.

Modern vent ecosystems may thus be living analogues of the earliest ecosystems on Earth — and by extension, potential models for life on icy moons such as Europa and Enceladus, where subsurface oceans may harbour hydrothermal activity.

5.9 Vent Temporal Dynamics & Succession

Hydrothermal vents are geologically ephemeral. On fast-spreading ridges (East Pacific Rise), individual vent sites persist for 10–100 years; on slow-spreading ridges (Mid-Atlantic Ridge), they may last 1,000–10,000 years. When a new vent opens, ecological succession follows a predictable pattern:

Microbial colonisation (days–weeks): free-living chemosynthetic bacteria (Thiomicrospira) form mats on chimney surfaces. These are the fastest colonisers of any deep-sea habitat.
Pioneer macrofauna (months): small polychaetes and gastropods settle. On the EPR, Tevnia jerichonana tube worms are typically first, colonising vigorous, high-temperature vents.
Climax community (1–5 years):Riftia replaces Tevnia as the dominant tube worm. Mussels (Bathymodiolus) colonise periphery. Species richness peaks.
Senescence (5–50 years): as vent flow decreases, mussels dominate over tube worms. Finally, suspension-feeders and background abyssal fauna replace vent-endemic species.

The 1991 eruption on the EPR at 9°50′N provided a natural experiment: an entire vent field was repaved by fresh lava, obliterating the previous community. Within 11 months, bacterial mats and Tevnia tube worms had colonised the site. Within 5 years, a full Riftia-dominated community had re-established, demonstrating remarkable resilience.

The colonisation rate $r$ can be modelled as a logistic function of time since eruption:

$ N(t) = \frac{K}{1 + \left(\frac{K - N_0}{N_0}\right) e^{-rt}} $

where $K$ is the carrying capacity, $N_0$ is the initial coloniser population, and $r \approx 2\text{--}5\;\text{yr}^{-1}$ for fast-growing vent species.

Hydrothermal Vent Ecosystem

A black smoker chimney with the major vent organisms, temperature gradient, and chemical zonation:

Simulation: Vent Chemistry & Biological Models

Four-panel simulation of chemosynthetic energy yields, vent fluid mixing gradients, Riftia carbon fixation kinetics, and ANME thermodynamic feasibility:

Hydrothermal Vents: Chemistry, Thermodynamics & Biology

Python

Energy yields, mixing models, carbon fixation, and ANME feasibility

script.py198 lines

import numpy as np
import matplotlib
matplotlib.use('Agg')
import matplotlib.pyplot as plt

fig, axes = plt.subplots(2, 2, figsize=(14, 11))
fig.patch.set_facecolor('#0a0a1a')
fig.suptitle('Hydrothermal Vents: Chemistry & Biology', color='white', fontsize=16, fontweight='bold', y=0.98)

# ===== Panel 1: Chemosynthetic Energy Yields =====
ax = axes[0, 0]
ax.set_facecolor('#0a0a2e')

reactions = [
    'CH$_4$ + 2O$_2$
(aerobic)',
    'H$_2$S + 2O$_2$
(sulfide)',
    'NH$_3$ + 1.5O$_2$
(nitrification)',
    'H$_2$ + 0.5O$_2$
(hydrogen)',
    '4Fe$^{2+}$ + O$_2$
(iron)',
    'CH$_4$ + SO$_4^{2-}$
(AOM)',
]
dG_values = [-818, -798, -275, -237, -65, -16.6]
colors = ['#ff6644', '#44ccaa', '#ffaa44', '#4488ff', '#aa66ff', '#ff4488']

bars = ax.barh(range(len(reactions)), [-dG for dG in dG_values], color=colors, alpha=0.8, height=0.6)
ax.set_yticks(range(len(reactions)))
ax.set_yticklabels(reactions, fontsize=9)
ax.set_xlabel('|\u0394G\u00b0| (kJ/mol)', color='white', fontsize=11)
ax.set_title('Chemosynthetic Energy Yields', color='white', fontsize=13, fontweight='bold')
ax.tick_params(colors='white')
ax.grid(True, alpha=0.15, color='#4444aa', axis='x')

# Annotate values
for i, (val, bar) in enumerate(zip(dG_values, bars)):
    ax.text(bar.get_width() + 10, i, f'{abs(val):.1f}', va='center', color='white', fontsize=9)

# Add ATP yield annotation
ax.axvline(x=30.5, color='#ffcc44', linewidth=1.5, linestyle='--', alpha=0.6)
ax.text(35, 5.5, '\u0394G for 1 ATP
(30.5 kJ/mol)', color='#ffcc44', fontsize=8)

# ===== Panel 2: Vent Fluid Mixing Model =====
ax = axes[0, 1]
ax.set_facecolor('#0a0a2e')

f = np.linspace(0, 0.3, 500)  # mixing fraction (vent fluid)
T_vent = 350  # C
T_sw = 2  # C
H2S_vent = 50  # mM
O2_sw = 0.15  # mM

T_mix = f * T_vent + (1 - f) * T_sw
H2S_mix = f * H2S_vent
O2_mix = (1 - f) * O2_sw

# Normalize for dual y-axis
ax_T = ax
ax_chem = ax.twinx()

l1, = ax_T.plot(f, T_mix, color='#ff4444', linewidth=2.5, label='Temperature (\u00b0C)')
ax_T.set_ylabel('Temperature (\u00b0C)', color='#ff4444', fontsize=11)
ax_T.tick_params(axis='y', labelcolor='#ff6644')

l2, = ax_chem.plot(f, H2S_mix, color='#44ccaa', linewidth=2, label='H\u2082S (mM)')
l3, = ax_chem.plot(f, O2_mix * 1000, color='#4488ff', linewidth=2, linestyle='--', label='O\u2082 (\u03bcM)')
ax_chem.set_ylabel('Concentration', color='white', fontsize=11)
ax_chem.tick_params(axis='y', labelcolor='white')

# Life zone
ax_T.axvspan(0.001, 0.1, alpha=0.15, color='#44cc88')
ax_T.text(0.05, 320, 'Life zone', color='#44cc88', fontsize=11, fontweight='bold', ha='center')

# Organism positions
for name, f_val, color in [('Mussels', 0.02, '#6688aa'), ('Riftia', 0.07, '#cc4444'), ('Alvinella', 0.2, '#cc8844')]:
    ax_T.axvline(x=f_val, color=color, linewidth=1, linestyle=':', alpha=0.6)
    ax_T.text(f_val, 10, name, color=color, fontsize=8, rotation=90, va='bottom')

ax.set_xlabel('Mixing Fraction f (vent fluid)', color='white', fontsize=11)
ax_T.set_title('Vent Fluid Mixing Model', color='white', fontsize=13, fontweight='bold')
ax.tick_params(axis='x', colors='white')
lines = [l1, l2, l3]
labs = [l.get_label() for l in lines]
ax.legend(lines, labs, facecolor='#1a1a3e', edgecolor='#4444aa', labelcolor='white', fontsize=9, loc='center right')

# ===== Panel 3: Riftia Carbon Fixation (Michaelis-Menten) =====
ax = axes[1, 0]
ax.set_facecolor('#0a0a2e')

CO2 = np.linspace(0, 500, 500)  # uM
Vmax = 15  # umol CO2 / g / h
Km = 50  # uM

A = Vmax * CO2 / (Km + CO2)

ax.plot(CO2, A, color='#cc4444', linewidth=2.5, label='$A = V_{max} [CO_2] / (K_m + [CO_2])$')
ax.axhline(y=Vmax, color='#ff8888', linewidth=1, linestyle='--', alpha=0.5, label=f'V$_{{max}}$ = {Vmax} \u03bcmol/g/h')
ax.axvline(x=Km, color='#88aaff', linewidth=1, linestyle='--', alpha=0.5, label=f'K$_m$ = {Km} \u03bcM')

# Mark half-Vmax
ax.plot(Km, Vmax/2, 'o', color='#88aaff', markersize=8, zorder=5)
ax.annotate(f'V$_{{max}}$/2 at K$_m$', xy=(Km, Vmax/2), xytext=(150, 5), color='#88aaff', fontsize=10,
            arrowprops=dict(arrowstyle='->', color='#88aaff', lw=1.5))

# Vent CO2 range
ax.axvspan(100, 500, alpha=0.1, color='#44cc88')
ax.text(300, 3, 'Typical vent
CO\u2082 range', color='#44cc88', fontsize=10, ha='center')

ax.set_xlabel('[CO\u2082] (\u03bcM)', color='white', fontsize=11)
ax.set_ylabel('Fixation Rate (\u03bcmol CO\u2082 g\u207b\u00b9 h\u207b\u00b9)', color='white', fontsize=11)
ax.set_title('Riftia Carbon Fixation Kinetics', color='white', fontsize=13, fontweight='bold')
ax.tick_params(colors='white')
ax.grid(True, alpha=0.15, color='#4444aa')
ax.legend(facecolor='#1a1a3e', edgecolor='#4444aa', labelcolor='white', fontsize=9)

# ===== Panel 4: ANME Thermodynamic Feasibility =====
ax = axes[1, 1]
ax.set_facecolor('#0a0a2e')

R = 8.314e-3  # kJ/(mol K)
T = 275  # K (deep sea ~2C)
dG0 = -16.6  # kJ/mol

# Vary CH4 and SO4 concentrations
CH4 = np.logspace(-4, -1, 200)  # mol/L (0.1 uM to 100 mM)
SO4_fixed = 0.028  # mol/L (28 mM seawater)
HCO3_fixed = 0.002  # mol/L
HS_fixed = 0.005  # mol/L

# deltaG = deltaG0 + RT ln Q
Q = (HCO3_fixed * HS_fixed) / (CH4 * SO4_fixed)
dG = dG0 + R * T * np.log(Q)

ax.plot(CH4 * 1000, dG, color='#ff4488', linewidth=2.5, label='\u0394G vs [CH\u2084]')
ax.fill_between(CH4 * 1000, dG, 0, where=dG < 0, alpha=0.15, color='#44cc88')
ax.fill_between(CH4 * 1000, dG, 0, where=dG >= 0, alpha=0.15, color='#ff4444')

ax.axhline(y=0, color='white', linewidth=1, linestyle='-', alpha=0.3)
ax.axhline(y=-10, color='#ffcc44', linewidth=1, linestyle='--', alpha=0.5)
ax.text(50, -8, 'Min. energy for
ATP synthesis
(~-10 kJ/mol)', color='#ffcc44', fontsize=8)

# Typical cold seep range
ax.axvspan(0.5, 10, alpha=0.1, color='#4488ff')
ax.text(3, 5, 'Cold seep
[CH\u2084] range', color='#4488ff', fontsize=9, ha='center')

ax.text(0.15, 10, 'Unfavourable
(\u0394G > 0)', color='#ff6644', fontsize=10, ha='center')
ax.text(30, -28, 'Thermodynamically
favourable', color='#44cc88', fontsize=10, ha='center')

ax.set_xlabel('[CH\u2084] (mM)', color='white', fontsize=11)
ax.set_ylabel('\u0394G (kJ/mol)', color='white', fontsize=11)
ax.set_title('AOM Thermodynamic Feasibility', color='white', fontsize=13, fontweight='bold')
ax.set_xscale('log')
ax.tick_params(colors='white')
ax.grid(True, alpha=0.15, color='#4444aa')
ax.legend(facecolor='#1a1a3e', edgecolor='#4444aa', labelcolor='white', fontsize=9)

plt.tight_layout(rect=[0, 0, 1, 0.96])
plt.savefig('output.png', dpi=150, bbox_inches='tight', facecolor='#0a0a1a')

# Print key results
print("=== Chemosynthetic Energy Yields ===")
for rxn, dg in zip(['CH4+O2', 'H2S+O2', 'NH3+O2', 'H2+O2', 'Fe2++O2', 'AOM'], dG_values):
    atp_yield = abs(dg) / 30.5
    print(f"  {rxn}: dG = {dg} kJ/mol, ~{atp_yield:.1f} ATP equivalent")

print("
=== Boiling Point Under Pressure ===")
for P in [1, 100, 250, 500]:
    T_boil = 100 + 28.6 * np.log(P) if P > 1 else 100
    print(f"  P = {P} atm: T_boil = {T_boil:.0f} C")

print("
=== Vent Mixing ===")
for frac in [0.001, 0.01, 0.05, 0.1, 0.2]:
    T = frac * 350 + (1-frac) * 2
    h2s = frac * 50
    o2 = (1-frac) * 150  # uM
    print(f"  f={frac}: T={T:.1f}C, H2S={h2s:.2f}mM, O2={o2:.0f}uM")

print("
=== AOM at Cold Seep ===")
CH4_seep = 0.001  # 1 mM
Q_seep = (0.002 * 0.005) / (CH4_seep * 0.028)
dG_seep = -16.6 + R * T * np.log(Q_seep)
print(f"  [CH4]=1mM: dG = {dG_seep:.1f} kJ/mol")
print(f"  This is {'feasible' if dG_seep < 0 else 'not feasible'} for life")
print(f"  ATP yield: ~{abs(dG_seep)/30.5:.2f} mol ATP per mol CH4")

Click Run to execute the Python code

Code will be executed with Python 3 on the server

References

Boetius, A., et al. (2000). A marine microbial consortium apparently mediating anaerobic oxidation of methane. Nature, 407(6804), 623–626.

Bright, M., & Lallier, F. H. (2010). The biology of vestimentiferan tubeworms. Oceanography and Marine Biology: An Annual Review, 48, 213–266.

Childress, J. J., & Fisher, C. R. (1992). The biology of hydrothermal vent animals: physiology, biochemistry, and autotrophic symbioses. Oceanography and Marine Biology: An Annual Review, 30, 337–441.

Corliss, J. B., et al. (1979). Submarine thermal springs on the Galapagos Rift. Science, 203(4385), 1073–1083.

Knittel, K., & Boetius, A. (2009). Anaerobic oxidation of methane: progress with an unknown process. Annual Review of Microbiology, 63, 311–334.

Lutz, R. A., & Kennish, M. J. (1993). Ecology of deep-sea hydrothermal vent communities: a review. Reviews of Geophysics, 31(3), 211–242.

Tivey, M. K. (2007). Generation of seafloor hydrothermal vent fluids and associated mineral deposits. Oceanography, 20(1), 50–67.

Van Dover, C. L. (2000). The Ecology of Deep-Sea Hydrothermal Vents. Princeton University Press.