Measurement of Earth’s Gravitational Acceleration on Antihydrogen with the ALPHA Experiment at CERN

Synopsis

Although the gravitational interaction between matter and antimatter has been the subject of theoretical speculation since antimatter was discovered in 1928, only recently has the ALPHA experiment at CERN been able to observe the effects of gravity on antimatter atoms, specifically antihydrogen.

This measurement is a test of the Weak Equivalence Principle, a foundational principle of Einstein’s general theory of relativity. The principle states that all masses respond identically to gravity, independent of internal structure.

To measure the gravitational acceleration of antihydrogen, anti-atoms are magnetically confined in the ALPHA apparatus and released under the combined influence of gravitational and magnetic forces.

The data analysis is based on a likelihood model of annihilation vertex positions. These vertices are produced when the magnetic confinement fields are lowered and antihydrogen atoms escape the trap.

The relevant likelihood parameter is the asymmetry between the number of anti-atoms released upward and the number released downward relative to the center of the electromagnetic trap.

The gravitational acceleration parameter is then obtained by regression against the data using a model derived from numerical simulations of antihydrogen motion inside the experiment.

The analysis includes statistical uncertainty estimation and systematic uncertainty treatment.

The results indicate that antihydrogen atoms behave consistently with gravitational attraction between antihydrogen and Earth. The experiment opens the path for more precise future tests of the Weak Equivalence Principle using anti-atoms.

Abstract

The gravitational interaction between matter and antimatter has been debated since antimatter’s theoretical prediction and experimental discovery. The ALPHA experiment at CERN, the Antihydrogen Laser Physics Apparatus, has performed a direct experimental observation of gravity acting on neutral antimatter atoms.

The experiment measures the response of antihydrogen to Earth’s gravitational field. Antihydrogen is the antimatter counterpart of hydrogen and consists of an antiproton bound to a positron.

The measurement tests whether antihydrogen responds to gravity in the same direction and approximate magnitude as ordinary matter, as required by the Weak Equivalence Principle.

Antihydrogen atoms are magnetically confined and later released by reducing the confining magnetic fields. The atoms escape the trap, annihilate on the apparatus walls, and produce detectable annihilation vertices.

The analysis constructs a likelihood based on these annihilation vertex positions. The core observable is the up-down asymmetry of released antihydrogen atoms.

The gravitational acceleration parameter is extracted through regression using a simulation-derived model of the antihydrogen release dynamics.

The result shows that antihydrogen behaves consistently with gravitational attraction toward Earth and does not support repulsive gravity at the tested scale.

Introduction

The observable universe is composed almost entirely of matter, while antimatter appears only in rare traces. Since antimatter’s prediction in 1928, its properties have been central to fundamental physics.

Two major antimatter research questions are directly relevant here:

• Whether antihydrogen energy levels match hydrogen energy levels, testing CPT symmetry.
• Whether antihydrogen responds to gravity like ordinary matter, testing the Weak Equivalence Principle.

The ALPHA experiment at CERN creates, traps, and studies antihydrogen. ALPHA has already performed several measurements of antihydrogen properties, including spectroscopy of antihydrogen energy levels for CPT symmetry testing and constraints on the net electric charge of antihydrogen.

In 2018, the ALPHA-g apparatus was constructed to measure Earth’s gravitational acceleration on antihydrogen. This thesis focuses on the analysis of data collected during the 2022 ALPHA-g data acquisition campaign.

The ALPHA-g measurement represents the first experimental test of the Weak Equivalence Principle using a neutral anti-atom.

Thesis Organization

Chapter 1 presents the matter-antimatter asymmetry problem and the theoretical background of the principles tested by ALPHA.

The chapter focuses on CPT symmetry, possible CPT violation, the Weak Equivalence Principle, and the concept of antigravity.

Chapter 2 describes the Antiproton Decelerator, the ELENA ring, and the ALPHA experimental apparatus.

The chapter focuses especially on the ALPHA-g apparatus, including the Penning-Malmberg trap embedded in an Ioffe-Pritchard trap and the physics of antihydrogen magnetic confinement.

Chapter 3 describes the technologies and techniques used to manipulate the charged plasmas required for antihydrogen synthesis. It also describes the 2022 gravity-measurement procedure.

Chapter 4 describes the data analysis performed on the 2022 dataset, including event cuts, likelihood construction, detector-efficiency estimation, simulation-derived regression, statistical uncertainties, and systematic uncertainties.

The results show that antihydrogen’s gravitational behavior is consistent with attraction between antihydrogen and Earth. Repulsive gravity between antihydrogen and Earth is ruled out at the tested level.

Chapter 1 — Matter-Antimatter Asymmetry, CPT Symmetry, and the Weak Equivalence Principle

One of the major unsolved problems in cosmology is the matter-antimatter asymmetry problem. It refers to the observed overabundance of matter relative to antimatter in the universe.

The standard cosmological model assumes that the universe was generated in the Big Bang. It is generally assumed that the early universe produced particles and antiparticles in nearly equal amounts.

However, cosmological observations show that the present universe is dominated by matter. This implies that some process created an asymmetry between matter and antimatter.

The observed matter asymmetry can in principle be explained by Sakharov’s three conditions:

• Violation of charge conjugation symmetry and charge-parity symmetry.
• Baryon number violation.
• Interactions occurring out of thermal equilibrium.

However, the amount of CP violation contained in the Standard Model appears insufficient to fully explain the observed asymmetry.

For this reason, ALPHA-2 studies antihydrogen energy levels to search for possible symmetry violations that could help explain why the universe is dominated by matter.

In parallel, some alternative cosmological theories propose that antimatter may experience gravitational acceleration differently from ordinary matter.

Such models attempt to explain matter dominance and may provide alternative perspectives on dark matter and dark energy. However, they are incompatible with the Weak Equivalence Principle if antimatter does not fall like matter.

The ALPHA-g experiment directly tests the Weak Equivalence Principle using antihydrogen.

1.1 Matter-Antimatter Asymmetry Problem

The standard cosmological model states that the universe began in the Big Bang and is presently composed approximately of:

Ordinary matter: 5%
Dark matter:     27%
Dark energy:     68%

It is widely assumed that the Big Bang produced particles and antiparticles in nearly equal quantities. This assumption leads to the matter-antimatter asymmetry problem because the observable universe is now almost entirely matter.

Within the Standard Model, elementary matter is composed of leptons and quarks, together with gauge bosons that mediate interactions.

Electrons, protons, and neutrons account for nearly all ordinary matter. These particles form the atoms observed throughout the universe.

The Big Bang theory arose from theoretical predictions and observational evidence including:

• Cosmic expansion.
• Cosmic microwave background radiation.
• Early-universe thermal plasma conditions.

In the early universe, when particle energies were high, pair production and annihilation reactions were in thermal equilibrium. Antimatter was therefore present in significant quantities.

As the universe expanded and cooled, particle energies eventually became too low to maintain pair production. Matter and antimatter annihilated, leaving behind a small residual excess of matter.

This residual matter excess forms the ordinary matter observed today.

Baryon Asymmetry

Within the Standard Model framework, the matter-antimatter asymmetry problem is usually expressed as the baryon asymmetry problem. This refers to the imbalance between baryonic matter and antibaryonic matter in the observable universe.

The baryon asymmetry can be defined using the difference between the number of baryons and antibaryons divided by their sum before antiprotons disappeared from the primordial plasma.

Since annihilation products are predominantly photons and no substantial antibaryon population exists today, the baryon-to-photon ratio is used to quantify this asymmetry.

η = NB / Nγ

A more explicit form relates the present baryon-to-photon ratio to the early baryon-antibaryon excess:

η = (NB − N̄B) / (NB + N̄B)

where:

• NB is the number of baryons.
• N̄B is the number of antibaryons.
• Nγ is the number of photons.
• η is the baryon-to-photon ratio.

The baryon-to-photon ratio is connected to the residual baryon density ΩB through:

ΩB ≃ η / (2.739 × 10⁻⁸ h²)

where h parameterizes the Hubble rate:

H0 = 100h km s⁻¹ Mpc⁻¹

The baryon asymmetry may be estimated from two independent cosmological measurements:

• Abundances of light elements in the intergalactic medium.
• Cosmic microwave background observations.

VillaUFO Administrative Scientific Finding

This thesis record is relevant because it connects antimatter gravity testing to three larger scientific questions:

• The origin of the matter-antimatter asymmetry.
• Possible CPT symmetry violation.
• Direct experimental testing of the Weak Equivalence Principle using neutral antimatter.

The ALPHA-g result is important for VillaUFO’s antigravity research archive because it provides direct experimental evidence against simple repulsive matter-antimatter gravity at the tested scale.

Measurement of Earth’s Gravitational Acceleration on Antihydrogen with the ALPHA Experiment at CERN

Supervisor: Prof. Germano Bonomi
Doctoral Candidate: Marta Urioni
Co-supervisor: Dr. Simone Stracka — INFN di Pisa

Synopsis

Although the gravitational interaction between matter and antimatter has been the subject of theoretical speculation since antimatter was discovered in 1928, only recently has the ALPHA experiment at CERN been able to observe the effects of gravity on antimatter atoms, specifically antihydrogen.

This measurement is a test of the Weak Equivalence Principle, a foundational principle of Einstein’s general theory of relativity. The principle states that all masses respond identically to gravity, independent of internal structure.

To measure the gravitational acceleration of antihydrogen, anti-atoms are magnetically confined in the ALPHA apparatus and released under the combined influence of gravitational and magnetic forces.

The data analysis is based on a likelihood model of annihilation vertex positions. These vertices are produced when the magnetic confinement fields are lowered and antihydrogen atoms escape the trap.

The relevant likelihood parameter is the asymmetry between the number of anti-atoms released upward and the number released downward relative to the center of the electromagnetic trap.

The gravitational acceleration parameter is then obtained by regression against the data using a model derived from numerical simulations of antihydrogen motion inside the experiment.

The analysis includes statistical uncertainty estimation and systematic uncertainty treatment.

The results indicate that antihydrogen atoms behave consistently with gravitational attraction between antihydrogen and Earth. The experiment opens the path for more precise future tests of the Weak Equivalence Principle using anti-atoms.

Abstract

The gravitational interaction between matter and antimatter has been debated since antimatter’s theoretical prediction and experimental discovery. The ALPHA experiment at CERN, the Antihydrogen Laser Physics Apparatus, has performed a direct experimental observation of gravity acting on neutral antimatter atoms.

The experiment measures the response of antihydrogen to Earth’s gravitational field. Antihydrogen is the antimatter counterpart of hydrogen and consists of an antiproton bound to a positron.

The measurement tests whether antihydrogen responds to gravity in the same direction and approximate magnitude as ordinary matter, as required by the Weak Equivalence Principle.

Antihydrogen atoms are magnetically confined and later released by reducing the confining magnetic fields. The atoms escape the trap, annihilate on the apparatus walls, and produce detectable annihilation vertices.

The analysis constructs a likelihood based on these annihilation vertex positions. The core observable is the up-down asymmetry of released antihydrogen atoms.

The gravitational acceleration parameter is extracted through regression using a simulation-derived model of the antihydrogen release dynamics.

The result shows that antihydrogen behaves consistently with gravitational attraction toward Earth and does not support repulsive gravity at the tested scale.

Introduction

The observable universe is composed almost entirely of matter, while antimatter appears only in rare traces. Since antimatter’s prediction in 1928, its properties have been central to fundamental physics.

Two major antimatter research questions are directly relevant here:

• Whether antihydrogen energy levels match hydrogen energy levels, testing CPT symmetry.
• Whether antihydrogen responds to gravity like ordinary matter, testing the Weak Equivalence Principle.

The ALPHA experiment at CERN creates, traps, and studies antihydrogen. ALPHA has already performed several measurements of antihydrogen properties, including spectroscopy of antihydrogen energy levels for CPT symmetry testing and constraints on the net electric charge of antihydrogen.

In 2018, the ALPHA-g apparatus was constructed to measure Earth’s gravitational acceleration on antihydrogen. This thesis focuses on the analysis of data collected during the 2022 ALPHA-g data acquisition campaign.

The ALPHA-g measurement represents the first experimental test of the Weak Equivalence Principle using a neutral anti-atom.

Thesis Organization

Chapter 1 presents the matter-antimatter asymmetry problem and the theoretical background of the principles tested by ALPHA.

The chapter focuses on CPT symmetry, possible CPT violation, the Weak Equivalence Principle, and the concept of antigravity.

Chapter 2 describes the Antiproton Decelerator, the ELENA ring, and the ALPHA experimental apparatus.

The chapter focuses especially on the ALPHA-g apparatus, including the Penning-Malmberg trap embedded in an Ioffe-Pritchard trap and the physics of antihydrogen magnetic confinement.

Chapter 3 describes the technologies and techniques used to manipulate the charged plasmas required for antihydrogen synthesis. It also describes the 2022 gravity-measurement procedure.

Chapter 4 describes the data analysis performed on the 2022 dataset, including event cuts, likelihood construction, detector-efficiency estimation, simulation-derived regression, statistical uncertainties, and systematic uncertainties.

The results show that antihydrogen’s gravitational behavior is consistent with attraction between antihydrogen and Earth. Repulsive gravity between antihydrogen and Earth is ruled out at the tested level.

Chapter 1 — Matter-Antimatter Asymmetry, CPT Symmetry, and the Weak Equivalence Principle

One of the major unsolved problems in cosmology is the matter-antimatter asymmetry problem. It refers to the observed overabundance of matter relative to antimatter in the universe.

The standard cosmological model assumes that the universe was generated in the Big Bang. It is generally assumed that the early universe produced particles and antiparticles in nearly equal amounts.

However, cosmological observations show that the present universe is dominated by matter. This implies that some process created an asymmetry between matter and antimatter.

The observed matter asymmetry can in principle be explained by Sakharov’s three conditions:

• Violation of charge conjugation symmetry and charge-parity symmetry.
• Baryon number violation.
• Interactions occurring out of thermal equilibrium.

However, the amount of CP violation contained in the Standard Model appears insufficient to fully explain the observed asymmetry.

For this reason, ALPHA-2 studies antihydrogen energy levels to search for possible symmetry violations that could help explain why the universe is dominated by matter.

In parallel, some alternative cosmological theories propose that antimatter may experience gravitational acceleration differently from ordinary matter.

Such models attempt to explain matter dominance and may provide alternative perspectives on dark matter and dark energy. However, they are incompatible with the Weak Equivalence Principle if antimatter does not fall like matter.

The ALPHA-g experiment directly tests the Weak Equivalence Principle using antihydrogen.

1.1 Matter-Antimatter Asymmetry Problem

The standard cosmological model states that the universe began in the Big Bang and is presently composed approximately of:

Ordinary matter: 5%
Dark matter:     27%
Dark energy:     68%

It is widely assumed that the Big Bang produced particles and antiparticles in nearly equal quantities. This assumption leads to the matter-antimatter asymmetry problem because the observable universe is now almost entirely matter.

Within the Standard Model, elementary matter is composed of leptons and quarks, together with gauge bosons that mediate interactions.

Electrons, protons, and neutrons account for nearly all ordinary matter. These particles form the atoms observed throughout the universe.

The Big Bang theory arose from theoretical predictions and observational evidence including:

• Cosmic expansion.
• Cosmic microwave background radiation.
• Early-universe thermal plasma conditions.

In the early universe, when particle energies were high, pair production and annihilation reactions were in thermal equilibrium. Antimatter was therefore present in significant quantities.

As the universe expanded and cooled, particle energies eventually became too low to maintain pair production. Matter and antimatter annihilated, leaving behind a small residual excess of matter.

This residual matter excess forms the ordinary matter observed today.

Baryon Asymmetry

Within the Standard Model framework, the matter-antimatter asymmetry problem is usually expressed as the baryon asymmetry problem. This refers to the imbalance between baryonic matter and antibaryonic matter in the observable universe.

The baryon asymmetry can be defined using the difference between the number of baryons and antibaryons divided by their sum before antiprotons disappeared from the primordial plasma.

Since annihilation products are predominantly photons and no substantial antibaryon population exists today, the baryon-to-photon ratio is used to quantify this asymmetry.

η = NB / Nγ

A more explicit form relates the present baryon-to-photon ratio to the early baryon-antibaryon excess:

η = (NB − N̄B) / (NB + N̄B)

where:

• NB is the number of baryons.
• N̄B is the number of antibaryons.
• Nγ is the number of photons.
• η is the baryon-to-photon ratio.

The baryon-to-photon ratio is connected to the residual baryon density ΩB through:

ΩB ≃ η / (2.739 × 10⁻⁸ h²)

where h parameterizes the Hubble rate:

H0 = 100h km s⁻¹ Mpc⁻¹

The baryon asymmetry may be estimated from two independent cosmological measurements:

• Abundances of light elements in the intergalactic medium.
• Cosmic microwave background observations.

VillaUFO Administrative Scientific Finding

This thesis record is relevant because it connects antimatter gravity testing to three larger scientific questions:

• The origin of the matter-antimatter asymmetry.
• Possible CPT symmetry violation.
• Direct experimental testing of the Weak Equivalence Principle using neutral antimatter.

The ALPHA-g result is important for VillaUFO’s antigravity research archive because it provides direct experimental evidence against simple repulsive matter-antimatter gravity at the tested scale.

1.1 Continued — Quark Mixing, CP Violation, and the Need for Beyond-Standard-Model Physics

The Standard Model contains three up-type quarks with electric charge +2/3 and three down-type quarks with electric charge −1/3. Up-type quarks include up, charm, and top. Down-type quarks include down, strange, and bottom.

Charged-current weak interactions allow up-type quarks to transition into down-type quarks, and down-type quarks to transition into up-type quarks.

These transition amplitudes are governed by the Cabibbo-Kobayashi-Maskawa matrix.

CKM matrix → transition amplitudes between quark generations

Cabibbo introduced early components of this mixing structure in 1963. Kobayashi and Maskawa later showed that charged-current weak interactions can violate CP symmetry when three quark generations are present.

The problem is that Kobayashi-Maskawa CP violation is too small to explain the observed matter-antimatter imbalance.

Even though the matter excess in the early universe may have been only about one part per billion, the CP violation supplied by the Standard Model is still many orders of magnitude too weak.

The Standard Model therefore contains mechanisms satisfying Sakharov’s three formal requirements, but those mechanisms do not appear strong enough to explain the observed baryon excess.

Beyond-Standard-Model Requirement

The prevailing interpretation is that the matter-antimatter asymmetry requires physics beyond the Standard Model.

Beyond-Standard-Model baryogenesis or leptogenesis scenarios may introduce:

• baryon number violation
• lepton number violation
• heavy gauge boson decay
• Higgs-sector extensions
• right-handed neutrino decay
• sterile neutrinos with heavy Majorana masses
• a first-order electroweak phase transition

These mechanisms can potentially generate the observed baryon asymmetry through baryogenesis or leptogenesis.

1.1.1 CPT Symmetry and CPT Violation

CPT invariance is the combined symmetry operation of charge conjugation, parity inversion, and time reversal.

CPT = C × P × T

Charge conjugation exchanges particles with antiparticles. Parity inversion reflects spatial coordinates.

P: r⃗ → −r⃗

Time reversal reverses the time coordinate.

T: t → −t

In known fundamental physics, CPT symmetry is closely tied to Lorentz symmetry. This connection makes CPT tests important probes of possible spacetime-symmetry violations.

CPT and the Sakharov Conditions

The standard Sakharov framework assumes CPT invariance. If CPT violation is permitted, a matter-antimatter asymmetry can in principle be generated even in thermal equilibrium.

In that case, CPT violation can replace the third Sakharov condition, which normally requires departure from thermal equilibrium.

However, even if CPT is violated, the model must still provide:

• C violation
• CP violation
• baryon-number or lepton-number violation

Standard Model Extension Framework

The Standard Model Extension is an effective field theory that incorporates:

• the Standard Model
• General Relativity
• possible Lorentz-violating spacetime operators
• possible CPT-violating spacetime operators

The SME is therefore a general framework for testing possible violations of CPT symmetry and the Weak Equivalence Principle while remaining compatible with established low-energy physics.

Minimal CPT-Violating Neutrino Term

A simplified leptogenesis model may include a minimal Lorentz- and CPT-violating neutrino term:

L = aμ(3) νL γμ νL

This term is:

• C violating
• CP odd
• CPT odd

The coefficient a₀⁽³⁾ acts like an effective chemical potential because it modifies the energy-momentum dispersion relation differently for neutrinos and antineutrinos.

At nonzero temperature, this produces different equilibrium number densities for neutrinos and antineutrinos.

μ = a0(3)

ηl ∼ μ TD²

where TD is the decoupling temperature.

When the universe cools and interaction rates fall below the Hubble expansion rate, the lepton number density freezes out.

ηl ∼ μ TD²

Since photon density scales approximately as:

ηγ ∼ T³

the lepton-to-photon ratio freezes at approximately:

ηl / ηγ ∼ a0(3) / TD

If the decoupling temperature is above the electroweak scale, sphaleron interactions can convert part of the lepton asymmetry into baryon asymmetry.

Three-Flavour Extension

For three neutrino flavours, the term generalizes to:

L = (aμ(3))ij νLi γμ νLj

This mechanism is conceptually similar to spontaneous leptogenesis, where a time-dependent vacuum expectation value of a scalar field plays the role of the effective CPT-violating coupling.

Experimental Constraint Problem

Current experimental limits on minimal SME coefficients make the simplest CPT-violating leptogenesis model too weak if the coefficients remain constant throughout cosmic history.

(a0(3))ij ≲ 10⁻²⁰ GeV

For decoupling temperatures above approximately:

TD > 100 GeV

the resulting asymmetry is many orders of magnitude too small.

Under the same assumptions, comparable constraints in the quark sector also rule out simple direct baryogenesis models of this type.

Non-Minimal SME Operators

Higher-dimension operators in the SME may avoid the minimal-model limitation because they can scale differently with decoupling temperature.

An example higher-order electron-sector coupling can be represented as:

L(5) = −aμρσ(5) ψ γμ ∂ρ ∂σ ψ

A term of this kind may generate a lepton asymmetry of approximate order:

|a0ρσ(5)| TD

Such non-minimal SME couplings are generally less constrained than minimal SME couplings.

Connection to Antihydrogen Spectroscopy

Nonrelativistic expansions of these operators can produce coefficients that affect the 1S-2S transition frequency of hydrogen and antihydrogen.

A nonzero measured coefficient in antihydrogen spectroscopy would be significant evidence for CPT-violation-based baryogenesis or leptogenesis models.

Hydrogen / antihydrogen 1S-2S comparison → CPT test channel

Administrative Scientific Finding

The thesis record establishes that Standard Model CP violation is insufficient to explain the observed matter-antimatter asymmetry. This motivates experimental searches for Beyond-Standard-Model effects through CPT tests and WEP tests.

ALPHA-2 contributes through antihydrogen spectroscopy and CPT symmetry testing. ALPHA-g contributes through direct gravitational testing of antihydrogen and the Weak Equivalence Principle.

For the VillaUFO antigravity research archive, the key finding is that speculative antimatter-gravity models must be judged against direct antihydrogen evidence. CPT-violation and WEP-violation models remain theoretically interesting, but the ALPHA-g result does not support simple repulsive gravity between antihydrogen and Earth.

Townsend Brown, the Biefeld–Brown Effect, and Electrogravitics

Thomas Townsend Brown was an American experimenter who investigated the interaction between high-voltage electrostatic fields and mechanical motion. Beginning in the 1920s, Brown reported that asymmetric capacitors subjected to high voltage appeared to generate a directional force. This phenomenon became known as the Biefeld–Brown effect, named after Brown and his mentor Paul Alfred Biefeld.

Brown proposed that strong electric fields interacting with dielectric materials could produce forces resembling gravitational behavior. The speculative engineering program based on this idea later became known as electrogravitics.

Basic Asymmetric Capacitor Configuration

Brown’s typical device consisted of two electrodes separated by a dielectric medium. One electrode had a small radius of curvature while the other had a much larger surface area. A high-voltage potential difference was applied across the electrodes.

The electric field between the electrodes can be approximated by:

E = V / d

where:

• E is the electric field strength
• V is the applied potential difference
• d is the separation distance between electrodes

Electrostatic Energy Density

The energy density stored in an electric field within a dielectric medium is given by:

u = ½ ε E²

where:

• u is the electric field energy density
• ε is the permittivity of the dielectric
• E is the electric field strength

Substituting the electric field expression:

u = ½ ε (V² / d²)

This expression shows that the stored field energy increases quadratically with voltage.

Electrostatic Force on Capacitor Plates

The electrostatic force between capacitor plates can be derived from the gradient of stored energy with respect to separation distance.

F = (ε A V²) / (2 d²)

where:

• F is the force between electrodes
• A is the electrode area
• ε is the dielectric permittivity
• V is applied voltage
• d is electrode separation

This relationship shows that electrostatic force scales with the square of the applied voltage.

Brown’s Empirical Scaling Law

Brown experimentally reported that the thrust produced by asymmetric capacitors appeared to follow a relation proportional to the square of the applied voltage and inversely proportional to electrode separation.

F ≈ k (V² / d)

where:

• F is the observed thrust force
• k is an experimentally determined constant dependent on geometry and dielectric properties
• V is applied voltage
• d is electrode separation

Electrohydrodynamic Interpretation

Modern experimental work has shown that most observed thrust from such asymmetric capacitor systems results from ionized air acceleration, commonly referred to as ionic wind.

The momentum transfer to the surrounding gas can be expressed as:

F = ṁ v

where:

• ṁ is the mass flow rate of ionized air
• v is the drift velocity of the ions

Ion drift velocity is related to the electric field through ion mobility:

v = μ E

where:

• μ is the ion mobility coefficient
• E is the electric field strength

Combining these relations yields:

F = ṁ μ E

Electrogravitic Hypothesis

Brown interpreted the observed forces as evidence that electric fields could couple directly to gravitational fields through polarization of the vacuum or dielectric medium. He proposed that sufficiently strong electric fields could distort local spacetime curvature.

In speculative electrogravitic models, the effective gravitational acceleration produced by an electric field would scale with the field energy density.

g_eff ∝ ε E²

This relationship resembles the dependence of stress-energy density on electromagnetic fields in general relativity.

Electromagnetic Contribution to Stress-Energy

In general relativity the electromagnetic field contributes to spacetime curvature through the stress-energy tensor:

Tμν = (1 / μ₀) (Fμα Fνα − ¼ gμν Fαβ Fαβ)

where:

• Tμν is the stress-energy tensor
• Fμν is the electromagnetic field tensor
• μ₀ is the vacuum permeability
• gμν is the metric tensor

The Einstein field equations relate stress-energy to spacetime curvature:

Gμν = (8πG / c⁴) Tμν

This formal relation explains why extremely strong electromagnetic fields theoretically contribute to gravitational curvature, although the magnitude of the effect in laboratory systems is extremely small.

Electrogravitics Research Programs

During the 1950s several aerospace companies and military contractors investigated electrogravitic propulsion concepts. The research explored whether high-voltage dielectric systems could generate lift or propulsion without conventional reaction mass.

Typical electrogravitic device parameters reported in historical experiments included:

Voltage range: 30 kV – 300 kV
Electrode spacing: 1 cm – 10 cm
Field strength: 10⁶ – 10⁷ V/m

Most modern studies conclude that the observed thrust in atmospheric experiments arises from electrohydrodynamic ion flow rather than direct gravitational coupling.

Scientific Status

The Biefeld–Brown effect remains experimentally observable in atmospheric conditions due to ion wind propulsion. However, controlled vacuum experiments have generally failed to reproduce significant thrust once ionized air is eliminated.

As a result, mainstream physics interprets the effect primarily as an electrohydrodynamic phenomenon rather than evidence of controllable gravitational manipulation.

Nevertheless, the electrogravitics hypothesis continues to appear in speculative propulsion literature because the stress-energy relation between electromagnetic fields and spacetime curvature suggests that sufficiently intense fields could, in principle, influence gravity.

Administrative Scientific Note

The Townsend Brown experiments are historically significant in the development of alternative propulsion concepts. Although the dominant explanation is ionic wind, the research remains relevant in discussions of high-field electromagnetism, dielectric stress forces, and theoretical electromagnetic contributions to spacetime curvature.

Additional Equations, Inputs, and Outputs for the Biefeld–Brown Effect and Electrogravitics

The earlier section covered the principal electrostatic field relations and the Einstein stress–energy coupling. Several additional equations commonly used in electrogravitics analysis were not yet included. These relate to capacitor force derivation, dielectric polarization forces, ion current generation, and power–thrust relationships.

Capacitance of the Asymmetric Capacitor

The capacitance of a parallel-plate approximation used in many Biefeld–Brown devices is:

C = ε A / d

Inputs

• ε — dielectric permittivity
• A — electrode surface area
• d — electrode separation

Output

• C — capacitance

Stored Electrical Energy

Energy stored in the capacitor electric field is:

U = ½ C V²

Substituting the capacitance equation:

U = ½ (ε A / d) V²

Inputs

• C — capacitance
• V — applied voltage

Output

• U — stored field energy

Force Derived from Energy Gradient

Electrostatic force can be derived from the gradient of stored energy with respect to distance.

F = − dU / dd

Substituting the capacitor energy relation yields:

F = (ε A V²) / (2 d²)

Inputs

• ε — dielectric permittivity
• A — electrode area
• V — applied voltage
• d — separation distance

Output

• F — electrostatic force between electrodes

Electric Field Near a Curved Electrode

Because Biefeld–Brown devices use asymmetric electrodes, the field near a small-radius electrode increases due to curvature.

E ≈ V / (r ln(b/a))

Inputs

• V — applied voltage
• r — radial distance
• a — radius of inner electrode
• b — radius of outer electrode

Output

• E — local electric field strength

Ionization Threshold Field

Ionization begins when the electric field exceeds the breakdown field of air.

E_breakdown ≈ 3 × 10⁶ V/m

Inputs

• gas pressure
• electrode geometry

Output

• threshold electric field for ionization

Ion Current Density

Current density produced by ionized air flow can be approximated as:

J = ρ_q μ E

Inputs

• ρ_q — charge density
• μ — ion mobility
• E — electric field strength

Output

• J — electric current density

Total Ion Current

I = J A

Inputs

• J — current density
• A — effective flow area

Output

• I — ion current

Electrohydrodynamic Thrust Model

Momentum transfer from accelerated ions produces thrust:

F = I d / μ

Inputs

• I — ion current
• d — electrode spacing
• μ — ion mobility

Output

• F — thrust force

Power Input

Electrical power consumed by the device is:

P = V I

Inputs

• V — applied voltage
• I — current

Output

• P — electrical input power

Thrust Efficiency

Electrohydrodynamic propulsion efficiency can be expressed as:

η = F v / P

Inputs

• F — thrust force
• v — ion drift velocity
• P — input power

Output

• η — propulsion efficiency

Ion Drift Velocity

v = μ E

Inputs

• μ — ion mobility
• E — electric field

Output

• v — ion drift velocity

Maxwell Stress Force

The electromagnetic force acting on a dielectric medium can also be expressed through the Maxwell stress tensor.

F = ∮ T · dA

where the tensor is:

T_ij = ε (E_i E_j − ½ δ_ij E²)

Inputs

• E — electric field components
• ε — permittivity

Output

• F — net electromagnetic force on the system

Electrogravitic Field Hypothesis

Some speculative electrogravitic models propose that gravitational acceleration may scale with electric field energy density:

g_eff = k ε E²

Inputs

• ε — permittivity
• E — electric field strength
• k — coupling coefficient

Output

• g_eff — effective gravitational acceleration

Einstein Electromagnetic Stress–Energy Relation

Electromagnetic energy contributes to gravitational curvature through:

G_μν = (8πG / c⁴) T_μν

Inputs

• T_μν — electromagnetic stress–energy tensor
• G — gravitational constant
• c — speed of light

Output

• spacetime curvature tensor G_μν

Summary of Practical Inputs for Experimental Electrogravitic Devices

Voltage:            30 kV – 300 kV
Electrode spacing:  0.5 cm – 10 cm
Electric field:     10⁶ – 10⁷ V/m
Current:            microamp to milliamp range
Power input:        watts to kilowatts

Typical Observable Outputs

Ion drift velocity
Ion current
Thrust force
Momentum transfer to surrounding gas
Electrical power consumption

These equations complete the commonly used analytical framework for Biefeld–Brown experimental systems and electrogravitic propulsion analysis.

Equation Audit — Townsend Brown / Biefeld–Brown / Electrogravitics

Equations that are physically valid

E = V / d

Works as a first-order parallel-plate approximation. It is not accurate for highly asymmetric wire-to-plate geometry unless treated as a local approximation.

C = εA / d

Works for an ideal parallel-plate capacitor. It is not exact for Brown-style asymmetric capacitors.

U = ½CV²

Works for stored capacitor energy.

u = ½εE²

Works for electric-field energy density.

F = εAV² / (2d²)

Works for ideal electrostatic attraction between capacitor plates. It does not by itself produce net propulsion of a closed device.

P = VI

Works for electrical input power.

v = μE

Works as a simplified ion-drift relation in a gas under mobility-limited conditions.

J = ρq μE

Works as a simplified conduction-current density model for drifting charge.

Tij = ε(EiEj − ½δijE²)

Works as the electrostatic Maxwell stress tensor expression.

Gμν = (8πG / c⁴)Tμν

Works as the Einstein field equation form relating stress-energy to curvature.

Equations that need correction or limitation

F = I d / μ

This is not a universal thrust equation. It is a simplified electrohydrodynamic scaling relation sometimes used for ionic wind devices, but real thrust depends on electrode geometry, gas pressure, ionization region, recombination, space charge, leakage current, and boundary flow.

E ≈ V / (r ln(b/a))

This works for a coaxial cylindrical approximation, not for a generic asymmetric lifter geometry.

F ≈ k(V² / d)

This is empirical only. It can fit some observations, but k absorbs too many variables to be predictive.

Equations that are speculative

g_eff ∝ εE²

g_eff = kεE²

These are not established physics. They are speculative electrogravitic hypotheses. No accepted laboratory evidence shows that Townsend Brown-type devices create measurable gravitational acceleration through this mechanism.

Bottom Line

The equations work for electrostatics, ion flow, Maxwell stress, and ordinary electromagnetic field energy. They do not prove antigravity. In air, Biefeld–Brown thrust is best explained as electrohydrodynamic ionic wind. In vacuum, the claimed thrust generally disappears or becomes too small to support the electrogravitic claim.

Clarion–Clipperton Zone Rock Nodules (Polymetallic Nodules)

The Clarion–Clipperton Zone (CCZ) is a large abyssal region of the Pacific Ocean between Hawaii and Mexico. The seabed is covered with polymetallic nodules sometimes called manganese nodules. These are rounded mineral concretions that form extremely slowly on the ocean floor, typically at growth rates of millimeters per million years.

These nodules do not convert water into air. Their formation results from precipitation of dissolved metal ions from seawater and sediment pore water, followed by oxidation reactions that produce solid metal oxides around a nucleus such as a shark tooth, shell fragment, or basalt particle.

Typical Chemical Composition

Average bulk composition of CCZ nodules:

MnO2      ≈ 25 – 35 %
Fe2O3     ≈ 10 – 20 %
SiO2      ≈ 5 – 15 %
Al2O3     ≈ 3 – 8 %
Ni        ≈ 1 – 1.5 %
Cu        ≈ 0.5 – 1 %
Co        ≈ 0.2 – 0.4 %
H2O       ≈ 10 – 20 % (bound water)

Main mineral phases:

• Birnessite (δ-MnO₂)
• Vernadite
• Todorokite
• Iron oxyhydroxides

Basic Chemical Formation Reactions

Manganese Oxidation

Mn²⁺ + ½O₂ + H₂O → MnO₂(s) + 2H⁺

Inputs

• dissolved Mn²⁺ ions
• oxygen from seawater
• water

Output

• solid manganese dioxide (nodule growth layer)

Iron Oxidation

4Fe²⁺ + O₂ + 6H₂O → 4Fe(OH)₃

Iron hydroxide then dehydrates to form iron oxides.

2Fe(OH)₃ → Fe₂O₃ + 3H₂O

Nickel and Copper Adsorption

Transition metals are incorporated through surface adsorption onto manganese oxide layers:

MnO₂(surface) + Ni²⁺ → MnO₂–Ni(surface complex)

MnO₂(surface) + Cu²⁺ → MnO₂–Cu(surface complex)

Growth Rate Model

Nodule growth rate depends on metal deposition flux.

R = (dM / dt) / ρ

Inputs

• dM/dt — mass deposition rate
• ρ — density of nodule material

Output

• R — radial growth rate

Typical observed value:

R ≈ 1 – 10 mm per million years

Electrochemical Properties

Manganese oxide minerals behave as redox-active materials. Their electrochemical potential can be approximated using the Nernst equation.

E = E° − (RT / nF) ln(Q)

Inputs

• E° — standard redox potential
• R — gas constant
• T — temperature
• n — number of electrons transferred
• F — Faraday constant
• Q — reaction quotient

Output

• E — electrode potential

Gas Production Possibilities

If electrical current passes through seawater around manganese oxides, water electrolysis can occur.

Electrolysis Reaction

2H₂O → 2H₂ + O₂

Inputs

• electric current
• water

Outputs

• hydrogen gas
• oxygen gas

However, this requires an external electrical energy source. The nodules themselves do not spontaneously convert water to gas.

Magnetic and Electrical Properties

Some nodules contain iron oxides and trace magnetite.

Magnetization may be approximated as:

M = χH

Inputs

• χ — magnetic susceptibility
• H — magnetic field strength

Output

• M — magnetization

Typical susceptibility values are very small, indicating weak paramagnetic behavior.

Possible Electrochemical Energy Storage

Because manganese oxides are used in batteries, nodules can participate in electrochemical reactions similar to battery cathodes.

MnO₂ + H⁺ + e⁻ → MnOOH

Inputs

• MnO₂
• protons
• electrons

Output

• manganese oxyhydroxide

Antigravity Properties

No known physical theory or experimental observation indicates that polymetallic nodules possess antigravity properties.

The gravitational behavior of these minerals is governed by ordinary mass density and Newtonian gravity.

Gravitational Force

F = G (m₁ m₂) / r²

Inputs

• m₁ — mass of the nodule
• m₂ — mass of Earth
• r — distance between centers

Output

• gravitational force

Because nodules are dense metal oxides, their gravitational interaction is normal and slightly stronger per volume than typical silicate rocks due to higher density.

Summary

• Clarion–Clipperton nodules are manganese–iron oxide concretions formed by extremely slow oxidation and precipitation processes.
• Their chemistry is dominated by manganese dioxide, iron oxides, and adsorbed transition metals.
• They can participate in electrochemical reactions similar to battery materials.
• They do not naturally convert water into air.
• No known antigravity mechanism exists for these materials.

Clarion–Clipperton Zone Nodules — Detailed Geochemistry, Structure, and Physical Properties

Polymetallic nodules in the Clarion–Clipperton Zone (CCZ) are complex layered mineral concretions formed by extremely slow chemical precipitation and redox reactions at the seafloor. Their structure is typically concentric, with alternating layers of manganese oxide and iron oxyhydroxide phases accumulating around a central nucleus.

Internal Layer Structure

A typical nodule contains:

• nucleus (basalt fragment, shark tooth, shell, sediment grain)
• manganese oxide layers
• iron oxyhydroxide layers
• trace metal adsorption layers
• clay and silica inclusions

Layer formation occurs by two primary mechanisms:

• hydrogenetic precipitation directly from seawater
• diagenetic precipitation from sediment pore fluids

Hydrogenetic Growth Reaction

Hydrogenetic growth occurs when dissolved metal ions oxidize in oxygenated seawater.

Mn²⁺ + ½O₂ + H₂O → MnO₂ + 2H⁺

Inputs

• Mn²⁺ ions in seawater
• dissolved oxygen
• water

Outputs

• manganese dioxide solid layer
• protons

Diagenetic Growth Reaction

Diagenetic growth occurs when metals diffuse upward from sediments and precipitate on the underside of nodules.

Mn²⁺(sediment pore water) + O₂ → MnO₂(s)

Typical Mineral Phases

Primary manganese minerals:

Birnessite      (δ-MnO₂)
Vernadite       (amorphous Mn oxide)
Todorokite      (MnO₂ tunnel structure)

Iron minerals:

FeOOH (goethite)
Fe₂O₃ (hematite)

Trace Metal Incorporation

Transition metals are captured by adsorption onto manganese oxide surfaces.

MnO₂(surface) + Ni²⁺ → MnO₂–Ni
MnO₂(surface) + Co²⁺ → MnO₂–Co
MnO₂(surface) + Cu²⁺ → MnO₂–Cu

These surface complexes give nodules their economic value as polymetallic ore deposits.

Elemental Distribution Model

Metal concentration profiles in nodules can be described by diffusion-limited deposition:

J = −D (dC/dx)

Inputs

• D — diffusion coefficient
• dC/dx — concentration gradient

Output

• J — diffusive metal flux

Growth Kinetics

Radial growth rate can be expressed as:

dr/dt = J / ρ

Inputs

• J — metal deposition flux
• ρ — mineral density

Output

• radial growth rate

Typical observed value:

dr/dt ≈ 1–5 mm per million years

Density and Mass

Average nodule density:

ρ ≈ 2.1 – 3.0 g/cm³

Mass can be calculated as:

m = ρV

For spherical approximation:

V = (4/3)πr³

Electrical Conductivity

Manganese oxides are semiconducting materials. Electrical conductivity follows:

σ = nqμ

Inputs

• n — charge carrier density
• q — electron charge
• μ — carrier mobility

Output

• σ — electrical conductivity

Electrochemical Redox Behavior

Manganese oxides can participate in redox reactions similar to battery cathodes.

MnO₂ + H⁺ + e⁻ → MnOOH

This reaction underlies alkaline battery chemistry.

Potential Difference Across Redox Interfaces

Electrochemical potential difference can be modeled using the Nernst equation.

E = E° − (RT/nF) ln(Q)

Inputs

• E° — standard potential
• R — gas constant
• T — temperature
• n — electrons transferred
• F — Faraday constant
• Q — reaction quotient

Magnetic Susceptibility

Iron oxides and trace magnetite create weak magnetic behavior.

M = χH

Inputs

• χ — magnetic susceptibility
• H — applied magnetic field

Output

• M — magnetization

Thermal Conductivity

Heat transfer through nodules follows Fourier’s law.

q = −k ∇T

Inputs

• k — thermal conductivity
• ∇T — temperature gradient

Output

• q — heat flux

Mechanical Strength

Stress response follows Hooke’s law for elastic deformation.

σ = E ε

Inputs

• E — Young’s modulus
• ε — strain

Output

• σ — mechanical stress

Electromagnetic Energy Density

If nodules are exposed to strong electromagnetic fields:

u = ½ (εE² + B²/μ)

Inputs

• E — electric field
• B — magnetic field
• ε — permittivity
• μ — permeability

Output

• electromagnetic energy density

Gravitational Interaction

The gravitational force on a nodule near Earth is:

F = mg

Inputs

• m — mass
• g ≈ 9.81 m/s²

Output

• gravitational force

Assessment of Exotic Properties

No experimental evidence shows that CCZ nodules possess:

• antigravity behavior
• gravity shielding
• spontaneous gas production
• propulsion properties

Their physical behavior is fully consistent with known geochemistry, electrochemistry, and mineral physics.

Clarion–Clipperton Zone Nodules — Advanced Geochemistry, Thermodynamics, and Physical Modeling

Polymetallic nodules represent a complex coupled system involving geochemistry, thermodynamics, electrochemistry, diffusion transport, and mineral crystallization. The following relations describe deeper physical mechanisms governing their formation and properties.

Thermodynamic Stability of Oxide Formation

Mineral formation is governed by Gibbs free energy minimization.

ΔG = ΔH − TΔS

Inputs

• ΔH — enthalpy change
• T — temperature
• ΔS — entropy change

Output

• ΔG — Gibbs free energy

A reaction proceeds spontaneously when:

ΔG < 0

Equilibrium Constant for Mineral Precipitation

ΔG° = −RT ln(K)

Inputs

• R — gas constant
• T — temperature
• K — equilibrium constant

Output

• standard Gibbs free energy change

Manganese Oxide Redox Potential

Redox potential for manganese oxide formation:

MnO₂ + 4H⁺ + 2e⁻ → Mn²⁺ + 2H₂O

Electrode potential follows the Nernst equation:

E = E° − (RT / nF) ln([Mn²⁺] / [H⁺]⁴)

Inputs

• E° — standard redox potential
• [Mn²⁺] — manganese ion concentration
• [H⁺] — proton concentration

Output

• E — redox potential

Diffusion Transport in Sediments

Metal ions migrate through sediment pores by diffusion.

Fick’s First Law

J = −D (dC/dx)

Inputs

• D — diffusion coefficient
• dC/dx — concentration gradient

Output

• J — diffusive flux

Fick’s Second Law

∂C/∂t = D ∂²C/∂x²

Inputs

• D — diffusion coefficient
• C — concentration field

Output

• time evolution of concentration

Adsorption of Trace Metals

Trace metals attach to manganese oxide surfaces following adsorption isotherms.

Langmuir Isotherm

θ = (KC) / (1 + KC)

Inputs

• K — adsorption constant
• C — metal ion concentration

Output

• θ — fraction of occupied adsorption sites

Surface Concentration

Γ = Γmax θ

Inputs

• Γmax — maximum surface coverage

Output

• adsorbed metal concentration

Crystal Growth Kinetics

Mineral layer growth can follow a reaction-rate controlled model:

dr/dt = k (C − Ceq)

Inputs

• k — kinetic growth constant
• C — dissolved metal concentration
• Ceq — equilibrium concentration

Output

• radial growth velocity

Porosity and Permeability of Sediment Matrix

Darcy’s Law

Q = −kA (ΔP / μL)

Inputs

• k — permeability
• A — cross-sectional area
• ΔP — pressure difference
• μ — fluid viscosity
• L — flow length

Output

• fluid flow rate

Electrical Conductivity of Nodules

Electron conduction occurs through hopping mechanisms in manganese oxides.

σ = σ₀ exp(−Ea / kBT)

Inputs

• σ₀ — conductivity constant
• Ea — activation energy
• kB — Boltzmann constant
• T — temperature

Output

• electrical conductivity

Electrochemical Capacitance of Oxide Surfaces

Surface capacitance can be approximated as:

C = εA / d

Inputs

• ε — permittivity
• A — surface area
• d — effective charge separation distance

Output

• capacitance

Magnetic Susceptibility of Iron Oxides

Iron-bearing phases contribute weak magnetic response.

χ = M / H

Inputs

• M — magnetization
• H — magnetic field

Output

• magnetic susceptibility

Thermal Diffusion in Nodules

Heat Diffusion Equation

∂T/∂t = α ∇²T

Inputs

• α — thermal diffusivity

Output

• temperature evolution

Radiogenic Heat Production

Trace uranium and thorium may generate small internal heat through radioactive decay.

Q = λN E

Inputs

• λ — decay constant
• N — number of radioactive atoms
• E — decay energy

Output

• heat generation rate

Gravitational Potential Energy

U = mgh

Inputs

• m — mass of nodule
• g — gravitational acceleration
• h — height relative to reference

Output

• gravitational potential energy

Stress from Overlying Sediment

σ = ρ g h

Inputs

• ρ — sediment density
• g — gravitational acceleration
• h — burial depth

Output

• compressive stress

Electromagnetic Energy Density in Ocean Environment

u = ½ (εE² + B²/μ)

Inputs

• E — electric field
• B — magnetic field
• ε — permittivity
• μ — permeability

Output

• electromagnetic energy density

Conclusion

The Clarion–Clipperton nodules are governed by well-understood chemical and physical processes including redox reactions, diffusion transport, adsorption chemistry, and mineral growth kinetics. These processes explain their composition and layered structure without invoking exotic mechanisms.

No known equation, experimental measurement, or theoretical framework indicates antigravity behavior in these materials.

Stubblefield–Tesla Earth Battery System

The earth battery concept uses the electrical potential differences naturally present in soil and groundwater. Two dissimilar conductors placed in moist ground form a galvanic system that converts chemical and electrochemical energy into electrical current.

Nathan Stubblefield demonstrated a coil-based earth battery in the late 1800s using buried coils composed of dissimilar metals. Nikola Tesla later proposed similar ground-coupled electrical systems in which the Earth itself functions as part of the electrical circuit.

Basic Galvanic Earth Battery Reaction

A buried galvanic pair functions similarly to a simple electrochemical cell.

Zn → Zn²⁺ + 2e⁻

(anode oxidation)

Cu²⁺ + 2e⁻ → Cu

(cathode reduction)

Inputs

• soil moisture
• dissolved electrolytes
• dissimilar metal electrodes
• temperature

Outputs

• electric current
• electric potential difference

Cell Voltage

The theoretical voltage follows the electrochemical potential difference between the electrodes.

Ecell = Ecathode − Eanode

Inputs

• standard electrode potentials
• electrolyte concentration

Output

• generated voltage

Nernst Equation (Real Soil Conditions)

E = E° − (RT / nF) ln(Q)

Inputs

• E° standard electrode potential
• R gas constant
• T temperature
• n electron transfer number
• F Faraday constant
• Q reaction quotient

Output

• actual cell voltage

Current Production

Electrical current depends on internal soil resistance.

I = V / R

Inputs

• V generated voltage
• R soil resistance

Output

• current

Soil Resistivity

Electrical conductivity of soil controls battery performance.

R = ρL / A

Inputs

• ρ soil resistivity
• L electrode separation
• A cross-section area

Output

• resistance of ground path

Power Output

P = IV

Inputs

• I current
• V voltage

Output

• electrical power

Electromagnetic Induction in Stubblefield Coil

Stubblefield’s design used bifilar coils of iron and copper wire buried underground. The coil can produce additional electrical activity through electromagnetic induction.

V = −N (dΦ / dt)

Inputs

• N number of coil turns
• Φ magnetic flux

Output

• induced voltage

Magnetic Flux

Φ = BA cos(θ)

Inputs

• B magnetic field strength
• A area of coil
• θ angle between field and surface

Output

• magnetic flux

Earth–Ionosphere Resonance (Tesla Concept)

Tesla believed the Earth could be used as a giant resonant electrical conductor.

f = c / (2πR)

Inputs

• c speed of light
• R Earth radius

Output

• resonant frequency

This relation approximates the Earth-ionosphere cavity resonance frequency.

Energy Stored in the Earth–Capacitor System

U = ½ CV²

Inputs

• C capacitance of Earth-ionosphere system
• V voltage difference

Output

• stored electric energy

Capacitance of Earth-Ionosphere System

C ≈ 4π ε₀ R

Inputs

• ε₀ vacuum permittivity
• R Earth radius

Output

• global capacitance

Electric Field Strength

E = V / d

Inputs

• V voltage
• d electrode spacing

Output

• electric field intensity

Earth Battery Power Density

Pd = σE²

Inputs

• σ soil conductivity
• E electric field

Output

• power density in soil

Summary

The Stubblefield–Tesla earth battery operates through a combination of electrochemical galvanic reactions, soil conductivity, and electromagnetic induction in buried coils. The system converts chemical energy from soil electrolytes and metal corrosion into low-voltage electrical power.

These systems typically generate millivolt to volt level potentials and milliamp currents depending on soil chemistry, moisture content, electrode spacing, and metal composition.

Stubblefield–Tesla Earth Battery — Full Math Process

This section models the earth battery as an electrochemical cell placed into conductive soil, with optional coil induction and Earth-ground coupling.

Core Inputs

Electrode 1: zinc or iron
Electrode 2: copper or carbon
Soil moisture: electrolyte medium
Electrode spacing: d
Electrode buried area: A
Soil resistivity: ρ
Cell voltage: Vcell
Load resistance: Rload
Internal resistance: Rint

Step A — Electrode Voltage

The open-circuit voltage is the difference between cathode and anode potentials.

Ecell = Ecathode − Eanode

For a zinc-copper approximation:

Ecathode = +0.34 V
Eanode = −0.76 V

Ecell = 0.34 − (−0.76)

Ecell = 1.10 V

Output:

Open-circuit voltage ≈ 1.10 V

Step B — Soil Internal Resistance

The ground path resistance is:

Rint = ρL / A

Example inputs:

ρ = 50 Ω·m
L = 1.0 m
A = 0.10 m²

Calculation:

Rint = (50 × 1.0) / 0.10

Rint = 500 Ω

Output:

Internal soil resistance = 500 Ω

Step C — Current Through a Load

Total circuit resistance:

Rtotal = Rint + Rload

Example:

Rint = 500 Ω
Rload = 1000 Ω

Rtotal = 500 + 1000

Rtotal = 1500 Ω

Current:

I = Vcell / Rtotal

I = 1.10 / 1500

I = 0.000733 A

I = 0.733 mA

Output:

Current ≈ 0.733 mA

Step D — Voltage Across the Load

Vload = I Rload

Vload = 0.000733 × 1000

Vload = 0.733 V

Output:

Load voltage ≈ 0.733 V

Step E — Power Output

P = I Vload

P = 0.000733 × 0.733

P = 0.000537 W

P = 0.537 mW

Output:

Power output ≈ 0.537 milliwatts

Step F — Maximum Power Transfer

Maximum power occurs when:

Rload = Rint

So:

Rload = 500 Ω

Maximum power:

Pmax = Vcell² / (4Rint)

Pmax = 1.10² / (4 × 500)

Pmax = 1.21 / 2000

Pmax = 0.000605 W

Pmax = 0.605 mW

Output:

Maximum available power ≈ 0.605 milliwatts

Step G — Chemical Reaction Rate

Faraday’s law connects current to material consumption.

m = (I t M) / (nF)

Inputs:

I = current
t = time
M = molar mass
n = electrons transferred
F = Faraday constant = 96485 C/mol

For zinc:

Zn → Zn²⁺ + 2e⁻

MZn = 65.38 g/mol
n = 2

Example: one day of operation:

I = 0.000733 A
t = 86400 s

m = (0.000733 × 86400 × 65.38) / (2 × 96485)

m = 4138.7 / 192970

m = 0.0214 g

Output:

Zinc consumed per day ≈ 0.0214 grams

Step H — Coil Induction Add-On

If the Stubblefield coil experiences a changing magnetic flux:

Vinduced = −N(dΦ/dt)

Magnetic flux:

Φ = B A cos(θ)

Example inputs:

N = 500 turns
ΔB = 20 μT = 20 × 10⁻⁶ T
A = 0.02 m²
Δt = 1 s
θ = 0°

Flux change:

ΔΦ = ΔB × A × cos(0)

ΔΦ = 20 × 10⁻⁶ × 0.02 × 1

ΔΦ = 4.0 × 10⁻⁷ Wb

Induced voltage:

Vinduced = −500 × (4.0 × 10⁻⁷ / 1)

Vinduced = −2.0 × 10⁻⁴ V

Vinduced = −0.2 mV

Output:

Induced voltage ≈ 0.2 millivolts

Step I — Total Voltage Estimate

If galvanic voltage and induction are in series:

Vtotal = Vcell + Vinduced

Vtotal ≈ 1.10 + 0.0002

Vtotal ≈ 1.1002 V

Output:

Total voltage remains dominated by the galvanic earth battery.

Step J — Earth-Ionosphere Capacitance Approximation

Cearth ≈ 4π ε0 Rearth

Inputs:

ε0 = 8.854 × 10⁻¹² F/m
Rearth = 6.371 × 10⁶ m

Calculation:

Cearth = 4π × 8.854 × 10⁻¹² × 6.371 × 10⁶

Cearth ≈ 7.09 × 10⁻⁴ F

Cearth ≈ 709 μF

Output:

Earth capacitance approximation ≈ 709 microfarads

Step K — Stored Earth-Capacitor Energy

U = ½ C V²

If the Earth-ionosphere potential is approximated as:

V = 300,000 V
C = 7.09 × 10⁻⁴ F

Then:

U = ½ × 7.09 × 10⁻⁴ × (300000)²

U = 0.0003545 × 90,000,000,000

U ≈ 31,905,000 J

Output:

Stored atmospheric electrical energy estimate ≈ 31.9 MJ

Administrative Finding

The Stubblefield earth battery works as a low-voltage electrochemical generator. The core output is produced by galvanic corrosion and soil-electrolyte conduction. Coil induction can add small transient voltage if the magnetic flux changes. Tesla-style Earth resonance concepts involve much larger ground-coupled electrical systems, but they do not make a simple buried earth battery a high-power source by themselves.

The practical earth-battery result is:

High voltage: no
High current: no
Continuous low-voltage current: yes
Chemical consumption: yes
Antigravity property: no verified mechanism

VillaEquation Registry — Electric Propulsion Physics

This registry organizes equations used in electric propulsion systems such as ion thrusters, Hall-effect thrusters, magnetoplasmadynamic engines, and plasma accelerators. The equations originate from classical electromagnetism, plasma physics, and rocket propulsion theory.

1. Lorentz Force

F = q(E + v × B)

2. Electric Field from Potential

E = −∇φ

3. Magnetic Force Component

F = q(v × B)

4. Particle Kinetic Energy from Electric Potential

½mv² = qV

5. Ion Exhaust Velocity

v = √(2qV/m)

6. Ion Beam Current

I = q n A v

7. Mass Flow Rate

ṁ = ρ A v

8. Thrust Equation

T = ṁ v_e

9. Rocket Momentum Equation

F = d(mv)/dt

10. Specific Impulse

Isp = v_e / g0

11. Propulsion Power

P = ½ ṁ v_e²

12. Beam Power

Pb = I V

13. Propulsive Efficiency

η = T v_e / (2P)

14. Rocket Equation

Δv = v_e ln(m0 / mf)

15. Plasma Current Density

J = nqv

16. Charge Continuity

∇·J + ∂ρ/∂t = 0

17. Maxwell Equation (Gauss Law)

∇·E = ρ / ε0

18. Maxwell Equation (Faraday)

∇ × E = −∂B/∂t

19. Maxwell Equation (Ampere-Maxwell)

∇ × B = μ0J + μ0ε0 ∂E/∂t

20. Plasma Frequency

ωp = √(ne e² / ε0 me)

21. Debye Length

λD = √(ε0 kTe / ne e²)

22. Thermal Velocity

v_th = √(2kT / m)

23. Plasma Pressure

p = nkT

24. Ion Sound Speed

Cs = √(kTe / mi)

25. Bohm Sheath Criterion

vi ≥ √(kTe / mi)

26. Plasma Sheath Potential

eφ ≈ kTe ln(√(mi / 2πme))

27. Child–Langmuir Space-Charge Limit

J = (4/9) ε0 √(2q/m) V^(3/2) / d²

28. Ion Optics Perveance

K = I / V^(3/2)

29. Beam Divergence

θ ≈ √(T_perp / T_parallel)

30. Beam Current Density

J = I / A

31. Hall Thruster Electron Drift

v_d = E × B / B²

32. Cyclotron Frequency

ωc = qB / m

33. Larmor Radius

rL = mv_perp / qB

34. Magnetized Plasma Momentum Equation

ρ dv/dt = −∇p + J × B

35. Plasma Conductivity

σ = ne² / (me ν)

36. Ohm’s Law for Plasma

J = σ(E + v × B)

37. Magnetic Diffusion Equation

∂B/∂t = η ∇²B

38. Alfven Velocity

VA = B / √(μ0ρ)

39. Alfven Mach Number

MA = v / VA

40. Magnetic Pressure

Pm = B² / (2μ0)

41. Plasma Beta

β = p / (B² / 2μ0)

42. Ionization Rate

R = n_e n_n ⟨σv⟩

43. Mean Free Path

λ = 1 / (n σ)

44. Collision Frequency

ν = v / λ

45. Magnetic Reynolds Number

Rm = μ0 σ v L

46. Power-to-Thrust Relation

T = 2P / v_e

47. Exhaust Power Flux

Φ = ½ ρ v_e³

48. Spacecraft Acceleration

a = T / m

49. Ion Beam Momentum Flux

Π = ρ v²

50. Electric Field from Charge Density

E = ρ / ε0

Summary

These equations form the mathematical framework of electric propulsion research, describing ion acceleration, plasma confinement, thrust generation, and beam dynamics used in modern spacecraft propulsion systems.

VillaEquation Engine — Electric Propulsion Simulator

VillaEquation Propulsion Laboratory

VillaPropulsionLab — VillaEquation Propulsion Laboratory

Administrative Record Framing Statement: This module preserves the existing VillaEquation propulsion calculator backbone and expands it into a deeper electric-propulsion laboratory. Each equation is rendered with inputs, substitution logic, solved output, and a browser-executed JavaScript function.

ε₀ = 8.854187817 × 10⁻¹² F/m
μ₀ = 1.25663706212 × 10⁻⁶ N/A²
e  = 1.602176634 × 10⁻¹⁹ C
k  = 1.380649 × 10⁻²³ J/K
mₑ = 9.1093837 × 10⁻³¹ kg
g₀ = 9.80665 m/s²

vₑ = √(2qV / mᵢ)
T = ṁvₑ
Isp = vₑ / g₀
P = ½ṁvₑ²
Ibeam = qṁ / mᵢ
ωc = qB / mᵢ

λD = √(ε₀kTe / ne²)
φs ≈ (kTe/e) ln √(mᵢ / 2πmₑ)

J = (4/9)ε₀√(2q/mᵢ)V^(3/2)/d²
K = I / V^(3/2)
θ ≈ √(T⊥ / T∥)

vd = E / B
ωc = qB / mᵢ
rL = mᵢv⊥ / qB
βH = ωc / ν

Pm = B² / 2μ₀
VA = B / √(μ₀ρ)
β = p / Pm
aMHD ≈ (J×B − ∇p) / ρ

λ = 1 / nσ
ν = v / λ
σp = ne² / mₑν
Rm = μ₀σp v L
Rion = ne nn <σv>

Δv = vₑ ln(m₀/mf)
tburn = mp / ṁ
a = T / m
distance ≈ ½at²

T/P = 2 / vₑ
η = T v / 2P
ηideal = useful jet power / input power

VillaEquation Antigravity Test Layer

Administrative Record Framing Statement

This section does not claim verified antigravity.

It tests whether a system produces acceleration not explained by known forces.

Core test:

a_anomalous = a_measured − a_known

where:

a_known = (T + FE + FB + Fbuoyancy − mg) / m

If a_anomalous ≈ 0, no antigravity is detected.
If a_anomalous > 0 after all known forces are removed, the result is an anomaly requiring independent verification.

Fg = mg
T/W = T / mg
a_net = (T − mg) / m

FE = qE
aE = FE / m

FB = qvB
aB = FB / m

a_known = (T + FE + FB + Fbuoyancy − mg) / m

a_anomalous = a_measured − a_known

VillaGravityControlLab

VillaSpacetimeLab

VillaUnifiedFieldLab

VillaQuantumVacuumLab

VillaQuantumFieldSpectrumLab

VillaCosmologyLab

VillaGravityAnomalyLab