Preprint
Article

This version is not peer-reviewed.

Path-Dependent Gravitation: Geodesic Levitation and Weight Anisotropy Induced by Spin-Curvature Coupling

Submitted:

29 January 2026

Posted:

30 January 2026

You are already at the latest version

Abstract
The Weak Equivalence Principle (WEP) postulates that the trajectory of a free-falling test body is independent of its internal structure. However, this universality formally applies only to structureless point masses. In this work, we re-examine the dynamics of extended spinning bodies within the Schwarzschild spacetime framework. Starting from the premise that gravity manifests as spacetime curvature---geometrically analogous to a slope---we derive that tangential motion modifies the effective geodesic path of an object.We demonstrate analytically that for a rotating body with angular momentum parallel to the local gravitational field, the geometric factors in the gravitational and inertial sectors of the radial geodesic equation exhibit an exact cancellation at the critical orbital velocity $v_c = \sqrt{GM/R}$. This implies that ``weight''---defined as the force required to deviate an object from its natural geodesic---is not an intrinsic invariant but a dynamic quantity dependent on the geometric alignment between the velocity vector and the spacetime curvature.We extend this finding to the microscopic regime, proposing that polarized atomic nuclei should exhibit weight anisotropy dependent on spin orientation. The thermodynamic consequence---a ``Geometric Buoyancy'' effect causing spontaneous stratification of spin-polarized gases---is derived and formally verified using the Lean~4 theorem prover.
Keywords: 
;  ;  ;  ;  ;  

1. Introduction

1.1. The Geometric Nature of Gravity

General Relativity fundamentally recast gravitation not as a force acting at a distance, but as the curvature of the spacetime manifold [1]. In this geometric paradigm, a stationary object resting on the surface of a massive sphere is not being “pulled” downward by an invisible force. Rather, it is being continuously accelerated upward—prevented from following its natural inertial path, which curves toward the center of the gravitating mass [2].
This natural path is called a geodesic: the trajectory of extremal proper time through curved spacetime. The profound insight of Einstein’s theory is that what we perceive as “gravitational force” is merely the sensation of being pushed away from geodesic motion. An apple falling from a tree is not being pulled down; it is finally free to follow its natural path. The ground, not the apple, is the accelerating reference frame.

1.2. The Slope Analogy

To build intuition for the arguments that follow, we introduce a geometric analogy that will prove central to our analysis. Consider spacetime curvature near a massive body as analogous to a curved slope or funnel. In this picture:
(1)
A stationary object on the ground is like a ball being held stationary on the side of a bowl. It requires continuous support (normal force) to prevent it from rolling toward the center.
(2)
A falling object is like a ball released on the slope—it follows the natural gradient toward the bottom.
(3)
An orbiting satellite is like a ball rolling around the inside of the bowl at precisely the right speed such that its tendency to “fall inward” is exactly balanced by its tangential motion carrying it “around” the curve.
The critical insight is this: tangential motion along a curved surface changes the effective “downward” tendency. A ball moving horizontally around the inside of a bowl does not fall toward the center as quickly as a stationary ball would. At sufficient speed, it maintains constant altitude indefinitely.

1.3. From Orbits to Spin

Classically, orbital motion (satellites) and local spin (gyroscopes) are treated as fundamentally distinct phenomena. A satellite traverses a macroscopic path through space, while a gyroscope’s rotation is confined to a local region. However, from the geometric perspective of General Relativity, both represent tangential transport of mass-energy along the curved manifold.
Consider a ring rotating about a vertical axis while resting on the Earth’s surface. Each mass element of the ring moves tangentially—horizontally, parallel to the ground. This is geometrically identical to orbital motion, merely confined to a smaller radius. If the tangential velocity reaches the orbital velocity at that radius, each mass element is instantaneously on a geodesic that does not intersect the ground.
This observation leads to our central hypothesis: the measurable “weight” of an object depends on the relationship between its internal motion and the local spacetime curvature. An object whose constituent mass elements move tangentially (parallel to equipotential surfaces) will exhibit reduced effective weight compared to a stationary object.

1.4. Scope and Structure

This paper develops the above intuition into a rigorous mathematical framework. Our analysis proceeds as follows:
  • Section 2 establishes the theoretical framework, deriving the geodesic equations in Schwarzschild spacetime and identifying the exact cancellation mechanism.
  • Section 3 derives the macroscopic weight reduction formula and establishes the conditions for geodesic levitation.
  • Section 4 extends the analysis to the microscopic regime, addressing the apparent paradox of nuclear spin velocities and introducing the geometric coupling factor.
  • Section 5 analyzes spin orientation effects and weight anisotropy.
  • Section 6 derives the thermodynamic consequences, specifically the “Geometric Buoyancy” effect.
  • Section 7 presents the formal verification of our mathematical derivations using the Lean 4 theorem prover.
  • Section 8 discusses experimental implications and proposes specific tests.
  • Section 9 addresses potential objections and reconciles our predictions with existing null results.
  • Section 10 provides theoretical implications and discussion.
  • Section 11 concludes with a summary of predictions and their broader implications.

2. Theoretical Framework

2.1. The Schwarzschild Metric

We consider the spacetime exterior to a spherically symmetric, non-rotating massive body of mass M. This geometry is described exactly by the Schwarzschild metric [3]. In coordinates
( t , r , θ , ϕ )
with geometric units ( G = c = 1 ), the line element is:
d s 2 = 1 2 M r d t 2 + 1 2 M r 1 d r 2 + r 2 d θ 2 + r 2 sin 2 θ d ϕ 2
For motion confined to the equatorial plane ( θ = π / 2 , d θ = 0 ), this simplifies to:
d s 2 = 1 2 M r d t 2 + 1 2 M r 1 d r 2 + r 2 d ϕ 2
The metric components are:
  • g t t = ( 1 2 M / r )
    — encodes gravitational time dilation
  • g r r = ( 1 2 M / r ) 1
    — encodes radial length contraction
  • g ϕ ϕ = r 2
    — encodes the geometry of angular motion

2.2. The Geodesic Equation

A freely falling particle follows a geodesic—a path that extremizes proper time. The geodesic equation is:
d 2 x μ d τ 2 + Γ α β μ d x α d τ d x β d τ = 0
where
τ
is proper time along the worldline, and
Γ α β μ
are the Christoffel symbols encoding the curvature of spacetime.
For the radial component ( μ = r ) in Schwarzschild spacetime, the relevant Christoffel symbols are:
Γ t t r = M r 2 1 2 M r
Γ r r r = M r 2 1 2 M r 1
Γ ϕ ϕ r = ( r 2 M ) = r 1 2 M r

2.3. The Radial Equation of Motion

For circular motion ( d r / d τ = 0 ,
d 2 r / d τ 2 = 0
for a geodesic), the radial geodesic equation becomes:
d 2 r d τ 2 + Γ t t r d t d τ 2 + Γ ϕ ϕ r d ϕ d τ 2 = 0
Introducing the four-velocity components
u t d t / d τ
and u ϕ d ϕ / d τ , we have:
d 2 r d τ 2 + Γ t t r ( u t ) 2 + Γ ϕ ϕ r ( u ϕ ) 2 = 0
Substituting the Christoffel symbols from Eqs. (4) and ():
d 2 r d τ 2 + M r 2 1 2 M r ( u t ) 2 r 1 2 M r ( u ϕ ) 2 = 0
This equation has a profound structure. The two terms represent:
1.
Gravitational Term:
Γ t t r ( u t ) 2 = M r 2 ( 1 2 M / r ) ( u t ) 2
— the “pull” toward the center, proportional to the time component of four-velocity squared.
2.
Centrifugal Term:
Γ ϕ ϕ r ( u ϕ ) 2 = r ( 1 2 M / r ) ( u ϕ ) 2
— the “push” away from the center, proportional to the angular component of four-velocity squared.

2.4. The Geometric Cancellation

The crucial observation is that the Schwarzschild factor
( 1 2 M / r )
appears linearly in both terms. This is not a coincidence but a deep geometric consequence of the metric structure.
For a stationary observer ( u ϕ = 0 ), there is no cancellation—the full gravitational term acts to accelerate the observer toward the center. But for an observer with tangential motion, the two terms can partially or fully cancel.
Factoring out the common term:
d 2 r d τ 2 = 1 2 M r M r 2 ( u t ) 2 r ( u ϕ ) 2
For an object outside the event horizon ( r > 2 M ), the factor ( 1 2 M / r ) > 0 . The condition for zero radial acceleration (geodesic levitation) is therefore:
M r 2 ( u t ) 2 = r ( u ϕ ) 2
Defining the coordinate angular velocity Ω u ϕ / u t :
Ω 2 = M r 3
This is precisely Kepler’s Third Law in relativistic form. The
( 1 2 M / r )
factors have canceled exactly, leaving a result that is formally identical to the Newtonian prediction.
Remark 1. 
The exact cancellation of the Schwarzschild factor is a key result. It demonstrates that the critical angular velocity for geodesic levitation isnotan approximation valid only in the weak-field limit, but an exact result of General Relativity.

3. Macroscopic Weight Reduction

3.1. Weight as Geodesic Deviation

We now establish a precise definition of “weight” within the geometric framework. Consider an object held stationary at radius
r
from the center of a massive body. The object is not following a geodesic—it requires a support force
F
to maintain its position.
Definition 1 
(Weight). Theweightof an object is the magnitude of the force required to prevent it from following its natural geodesic path.
For a stationary object, this force must counteract the entire gravitational acceleration. For an object with tangential motion, the required force is reduced because the natural geodesic is “less radial.”

3.2. The Weight Reduction Formula

From Equation (10), the radial acceleration for an object with angular velocity
Ω
is:
a r = 1 2 M r M r 2 ( u t ) 2 r ( Ω · u t ) 2
For a stationary object ( Ω = 0 ):
a r ( static ) = 1 2 M r M r 2 ( u t ) 2
The ratio of weight at angular velocity
Ω
to static weight is:
W ( Ω ) W static = a r ( Ω ) a r ( static ) = 1 r 3 Ω 2 M = 1 Ω 2 Ω c 2
where
Ω c 2 = M / r 3
is the critical angular velocity for geodesic levitation.
Converting to tangential velocity
v = r Ω
and restoring SI units:
W ( v ) = W static 1 v 2 v c 2
where
v c = G M / R
is the critical velocity (first cosmic velocity,
7.9
km/s for Earth).

3.3. Physical Interpretation

Equation (16) has several important implications:
1.
Continuous Reduction: Weight decreases parabolically with tangential velocity. There is no threshold effect—any tangential motion reduces effective weight.
2.
Orbital Equivalence: At v = v c , weight vanishes entirely. This is the condition for a circular orbit: the object is in free fall, following a geodesic that happens to maintain constant radius.
3.
Supercritical Motion: For v > v c , the formula predicts “negative weight”—a net outward force would be required to maintain circular motion. This corresponds to the requirement for a centripetal force beyond gravity’s provision (binding the object to prevent escape).
4.
Orientation Dependence: The derivation assumes tangential motion parallel to the equipotential surface. Radial motion does not contribute to this cancellation.

3.4. The Spinning Ring Model

Consider a thin ring of mass
m
and radius
ρ R
rotating with angular velocity
ω
about a vertical axis while resting on the Earth’s surface. Each mass element has tangential velocity v = ω ρ .
The effective weight of the ring is:
W ring = m g 1 ω 2 ρ 2 v c 2 = m g 1 ω 2 ρ 2 R G M
For a ring with
ρ = 1
m rotating at
ω = 7900
rad/s (giving v = v c ), the weight would vanish. However, such extreme angular velocities are impractical for macroscopic objects due to material stress limits.

3.5. Time Dilation at Levitation

For an object in geodesic levitation (satisfying Ω 2 = M / r 3 ), we can derive the time dilation factor. The four-velocity normalization condition is:
g μ ν u μ u ν = 1
In the equatorial plane:
1 2 M r ( u t ) 2 r 2 ( u ϕ ) 2 = 1
Substituting the levitation condition
u ϕ = Ω · u t
with Ω 2 = M / r 3 :
1 2 M r ( u t ) 2 r 2 · M r 3 ( u t ) 2 = 1
Simplifying:
1 2 M r M r ( u t ) 2 = 1 3 M r ( u t ) 2 = 1
Therefore:
( u t ) 2 = 1 1 3 M / r
This result requires
r > 3 M
for the solution to be physical (timelike worldline). The surface
r = 3 M
is the photon sphere, where circular null geodesics exist.
Remark 2. 
An object in geodesic levitation experiences greater time dilation than a stationary observer at the same radius. This additional dilation arises from the kinetic contribution to the metric.

4. Microscopic Extension: Nuclear Spin and the Velocity Paradox

4.1. Intrinsic Angular Momentum of Elementary Particles

Having established the macroscopic framework for geodesic levitation, we now extend the analysis to the microscopic regime. Elementary particles and atomic nuclei possess intrinsic angular momentum—spin—which is a fundamental quantum property [4].
For a particle with spin quantum number s, the magnitude of the spin angular momentum is:
S = s ( s + 1 )
For spin-1/2 particles (protons, neutrons, electrons), this gives S = 3 2 .
While spin is fundamentally a quantum property without a classical analog, we can construct a semi-classical estimate of the “effective surface velocity” by treating the particle as a rotating sphere. For a proton with mass
m p
and charge radius
r p 0.87 × 10 15
m:
v s S m p r p = 2 m p r p 1.055 × 10 34 2 × 1.67 × 10 27 × 0.87 × 10 15
v s 3.6 × 10 7 m / s 0.12 c

4.2. The Velocity Paradox

This semi-classical surface velocity exceeds Earth’s critical levitation velocity by an enormous factor:
v s v c 3.6 × 10 7 7.9 × 10 3 4600
If the geometric coupling derived in Section 3 applied directly and without modification to nuclear spin, Equation (16) would predict:
W W static = 1 v s v c 2 1 2.1 × 10 7 2 × 10 7
This would imply an enormous “anti-gravity” effect—each proton would experience an outward force millions of times its normal weight. The manifest stability of atomic nuclei and the absence of spontaneous levitation of matter immediately falsifies this naive extrapolation.

4.3. Resolution: The Geometric Coupling Factor

The resolution of this paradox lies in recognizing that the geodesic analysis of Section 2 implicitly assumes that the spatial extent of the rotating system is comparable to the radius of curvature of spacetime. For a satellite orbiting Earth, this condition is satisfied: the orbital radius is comparable to Earth’s radius.
However, for a proton, the situation is radically different:
r p R Earth 10 15 10 7 = 10 22
The proton is
10 22
times smaller than the radius of curvature it is trying to “orbit.” At this scale, spacetime appears effectively flat over the extent of the particle. The curvature that would enable geodesic levitation is simply not resolved by the particle’s wave function.
We formalize this through a Geometric Coupling Factor  η , which quantifies the efficiency of spin-curvature coupling:
Δ W W = η v spin v c 2
The factor
η
must satisfy
η 1
for microscopic systems and
η 1
for macroscopic orbital motion. We propose that
η
scales with the ratio of the system size to the curvature radius:
η f r system R curvature
where
f
is a function satisfying
f ( x ) 0
as
x 0
and
f ( x ) 1
as x 1 .

4.4. Dimensional Analysis of the Coupling Factor

Several scaling relations are plausible on dimensional grounds:

4.4.0.1. Linear Scaling:

η r p R 10 22
This would give:
Δ W W 10 22 × 2 × 10 7 = 2 × 10 15

4.4.0.2. Quadratic Scaling:

η r p R 2 10 44
This would give:
Δ W W 10 44 × 2 × 10 7 = 2 × 10 37
The linear scaling produces effects potentially within reach of future precision gravimetry ( 10 15 ), while quadratic scaling would render the effect permanently unobservable.

4.5. Connection to Mathisson-Papapetrou Equations

Our phenomenological approach can be connected to the rigorous framework of spinning body dynamics in General Relativity. The Mathisson-Papapetrou-Dixon (MPD) equations [5,6] describe the motion of extended spinning bodies:
D p μ d τ = 1 2 R μ ν α β u ν S α β
D S μ ν d τ = p μ u ν p ν u μ
where
p μ
is the four-momentum,
S μ ν
is the spin tensor, and
R μ ν α β
is the Riemann curvature tensor.
Equation (35) shows explicitly that spin couples to curvature through the Riemann tensor. The magnitude of this coupling scales as:
D p μ d τ R μ ν α β · S M r 3 · S
For a proton in Earth’s field:
M Earth R Earth 3 6 × 10 24 ( 6 × 10 6 ) 3 3 × 10 4 kg / m 3
The spin-curvature force is thus:
F spin G M R 3 · S G M R 3 · 2
Comparing to gravitational force F g = m g = m · G M / R 2 :
F spin F g 2 m R 10 34 2 × 10 27 × 10 7 10 14 to 10 15
This independent estimate supports the linear scaling hypothesis and suggests effects at the
10 14
to
10 15
level may be physical.

5. Spin Orientation and Weight Anisotropy

5.1. The Geometric Alignment Condition

The geodesic cancellation mechanism derived in Section 2 has a crucial directional requirement: the tangential velocity must lie along the curvature of the manifold—that is, parallel to equipotential surfaces (horizontal).
For a macroscopic orbiting satellite, this condition is automatically satisfied: orbital motion is tangential to the gravitational potential. But for a spinning particle, the orientation of the spin vector determines whether the internal motion is horizontal” or vertical” relative to the gravitational field.

5.2. Spin Vector Orientations

Consider a nucleus with spin angular momentum vector S . The constituent quarks and gluons (in a semi-classical picture) circulate in planes perpendicular to S . We distinguish two limiting cases:

5.2.0.3. Case A: Vertical Spin ( S g )

When the spin vector is aligned with the gravitational field (pointing up or down), the internal mass-energy circulates in horizontal planes. This configuration maximizes the geometric alignment:
  • Circulation is tangent to equipotential surfaces
  • Internal motion mimics “micro-orbital” motion
  • Geodesic cancellation is maximally effective
Prediction: Reduced effective weight.

5.2.0.4. Case B: Horizontal Spin ( S g )

When the spin vector is perpendicular to the gravitational field (pointing horizontally), the internal mass-energy circulates in vertical planes. This configuration minimizes the geometric alignment:
  • Circulation crosses equipotential surfaces
  • Internal motion alternately falls” and rises”
  • No net geodesic cancellation
Prediction: Standard weight (no reduction).

5.3. The Weight Anisotropy Tensor

We can formalize this directional dependence through a weight anisotropy tensor. Let
g ^
be the unit vector in the gravitational direction and
s ^
be the unit vector along the spin axis. The effective weight is:
W eff = W 0 1 η v s v c 2 ( s ^ · g ^ ) 2
This formula interpolates between:
  • Full reduction when S g :
    ( s ^ · g ^ ) 2 = 1
  • No reduction when S g :
    ( s ^ · g ^ ) 2 = 0

5.4. Implications for Polarized Matter

For a sample of spin-polarized nuclei with polarization fraction
P
along direction n ^ :
( s ^ · g ^ ) 2 = P 2 ( n ^ · g ^ ) 2 + 1 P 2 3
The factor of
1 / 3
arises from averaging
cos 2 θ
over random orientations.
For fully polarized matter ( P = 1 ) with polarization along gravity ( n ^ = g ^ ):
( s ^ · g ^ ) 2 = 1
For fully polarized matter with polarization perpendicular to gravity:
( s ^ · g ^ ) 2 = 0
For unpolarized matter:
( s ^ · g ^ ) 2 = 1 3
The weight difference between vertically and horizontally polarized samples is:
Δ W = W W = W 0 · η v s v c 2

6. Thermodynamic Consequences: Geometric Buoyancy

6.1. The Boltzmann Distribution in a Gravitational Field

The equilibrium distribution of particles in a gravitational field follows the Boltzmann distribution [8]:
n ( h ) = n 0 exp m g h k B T
where
n ( h )
is the number density at height h,
n 0
is the density at the reference level,
m
is the particle mass,
g
is gravitational acceleration,
k B
is Boltzmann’s constant, and
T
is temperature.
The characteristic scale height
H
is defined as the height at which density falls to
1 / e
of its surface value:
H = k B T m g

6.2. Modified Scale Height for Spinning Particles

If the effective gravitational acceleration depends on spin orientation, the scale height becomes:
H eff = k B T m g eff = k B T m g ( 1 Δ )
where
Δ = η ( v s / v c ) 2 ( s ^ · g ^ ) 2
is the fractional weight reduction.
For small Δ :
H eff H 0 ( 1 + Δ )

6.3. The Geometric Buoyancy Effect

Consider a column of gas containing two species of spin-polarized particles:
  • Species A: Spin aligned vertically ( S g )
  • Species B: Spin aligned horizontally ( S g )
From Equation (49), Species A has a larger scale height than Species B:
H A H B = 1 Δ B 1 Δ A = 1 1 η ( v s / v c ) 2 > 1
This means that at any given altitude h > 0 , the density ratio is:
n A ( h ) n B ( h ) = n A , 0 n B , 0 exp m g h k B T 1 1 Δ A 1 1 Δ B
For equal surface densities and small Δ :
n A ( h ) n B ( h ) exp m g h Δ A k B T > 1
Physical Interpretation: Vertically-polarized particles are “lighter” and therefore preferentially found at higher altitudes. This is analogous to conventional buoyancy, but driven by geometric spin-curvature coupling rather than density differences. We term this Geometric Buoyancy.

6.4. Quantitative Estimate

For 129 Xe atoms (a common hyperpolarizable species) at room temperature:
  • m = 129 × 1.66 × 10 27
    kg
    = 2.14 × 10 25
    kg
  • T = 300
    K
  • g = 9.8
    m/s
    2
  • H 0 = k B T / m g 2.0 × 10 3
    m
    2
    km
The relative density enhancement at height
h
is:
n A n B n B h H 0 Δ = h H 0 × η × v s v c 2
For
h = 1
m and Δ 10 15 :
n A n B n B 1 2.0 × 10 3 × 10 15 5 × 10 19
This is an extremely small effect, likely undetectable with current technology. However, the effect scales with height and could potentially be enhanced in tall columns or centrifuge configurations.

6.5. Formal Statement of Geometric Buoyancy

We can state the geometric buoyancy effect as a theorem:
Theorem 1 
(Geometric Buoyancy). For a gas of spinning particles in thermal equilibrium within a gravitational field, if the effective gravitational acceleration is reduced by spin-curvature coupling, then the equilibrium density at any positive height is strictly greater than it would be for non-spinning particles of the same mass.
Mathematically:
g eff < g n spin ( h ) > n static ( h ) h > 0
This theorem will be formally verified in the next section.

7. Formal Verification via Lean 4

7.1. Motivation for Formal Verification

The mathematical derivations presented in this paper involve algebraic manipulations that, while straightforward, could potentially harbor subtle errors. To ensure absolute rigor, we have implemented formal proofs of our core results using the Lean 4 theorem prover [9].
Lean is an interactive theorem prover and programming language based on dependent type theory. Proofs in Lean are machine-checked: if the code compiles without errors, the theorem is guaranteed to be a valid logical consequence of the axioms and definitions.

7.2. Definitions

We begin by formalizing the key physical quantities:
Listing 1. Core Definitions in Lean 4
Preprints 196712 i001

7.3. Theorem 1: Geodesic Cancellation Implies Kepler’s Law

The first theorem verifies that the levitation condition (zero radial acceleration) implies the Keplerian angular velocity relation, with exact cancellation of the Schwarzschild factor.
Listing 2. Theorem 1: Geodesic Cancellation
Preprints 196712 i002

7.4. Theorem 2: Weight Reduction Formula

The second theorem verifies that weight decreases parabolically with angular velocity according to our derived formula.
Listing 3. Theorem 2: Weight Reduction Formula
Preprints 196712 i003

7.5. Theorem 3: Levitation Time Dilation

The third theorem derives the time dilation factor for an object in geodesic levitation.
Listing 4. Theorem 3: Levitation Time Dilation
Preprints 196712 i004a
Preprints 196712 i004b

7.6. Theorem 4: Geometric Buoyancy

The fourth theorem formally verifies that reduced effective gravity leads to increased density at altitude.
Listing 5. Theorem 4: Geometric Buoyancy
Preprints 196712 i005

7.7. Auxiliary Theorems: Levitation Equivalence

We also provide auxiliary theorems establishing the equivalence between the levitation condition and Kepler’s law.
Listing 6. Auxiliary Definitions and Theorems
Preprints 196712 i006

7.8. Verification Summary

All theorems compile successfully in Lean 4 with Mathlib, confirming:
Table 1. Summary of Formally Verified Theorems
Table 1. Summary of Formally Verified Theorems
Theorem Physical Content Status
1. Geodesic Cancellation ( 1 2 M / r ) cancels exactly; Ω 2 = M / r 3 Verified
2. Weight Reduction W ( Ω ) = W 0 ( 1 Ω 2 / Ω c 2 ) Verified
3. Time Dilation ( u t ) 2 = 1 / ( 1 3 M / r ) at levitation Verified
4. Geometric Buoyancy Spinning gas has higher density at altitude Verified
5. Levitation Equivalence Levitation ⇔ Keplerian velocity Verified
The formal verification provides mathematical certainty that our derivations contain no algebraic errors and that the conclusions follow necessarily from the premises.
Remark 3 
(Asymmetry in Hypotheses). The forward direction (levitation
Kepler) requires
r > 2 M
to ensure the Schwarzschild factor
( 1 2 M / r )
is non-zero and can be cancelled. The backward direction (Kepler
levitation) only requires r > 0 , as the algebraic identity holds regardless of whether we are inside or outside the horizon.

8. Experimental Proposals

8.1. Direct Weight Measurement: Polarized Gravimetry

The most straightforward test of spin-induced weight anisotropy is direct measurement of weight changes upon spin reorientation.

8.1.0.5. Experimental Configuration:

1.
Sample: A macroscopic quantity (1–100 cm 3 ) of hyperpolarized noble gas, preferably 129 Xe or 3 He, chosen for their long spin relaxation times and well-established polarization techniques [10].
2.
Polarization Method: Spin-exchange optical pumping (SEOP) or metastability-exchange optical pumping (MEOP) to achieve nuclear polarization P > 50 % .
3.
Containment: Non-magnetic, thermally stable container (e.g., aluminosilicate glass) suspended in a high-precision gravimeter.
4.
Gravimeter: Superconducting gravimeter or atom interferometer with sensitivity
δ g / g 10 12
or better [11].
5.
Modulation: Adiabatic rotation of the holding magnetic field
B 0
between vertical ( B g ) and horizontal ( B g ) orientations at periods of 10–100 seconds.

Expected Signal:

For
N
polarized nuclei with polarization P:
Δ W = N m n g · η v s v c 2 · P 2
For 1 mole of 129 Xe ( N 6 × 10 23 ,
m 2 × 10 25
kg) with P = 0.5 :
Δ W 6 × 10 23 × 2 × 10 25 × 10 × η × 2 × 10 7 × 0.25
Δ W 6 × 10 6 × η Newtons
For
η 10 22
(linear scaling):
Δ W 6 × 10 16 N 60 femtonewtons
This corresponds to a fractional weight change:
Δ W W η × 2 × 10 7 × 0.25 5 × 10 16
Current superconducting gravimeters achieve δ g / g 10 12 , which is four orders of magnitude short of the predicted signal. However, several enhancement strategies are possible.

8.2. Enhancement Strategies

A. Coherent Accumulation

If the weight modulation is periodic, lock-in detection over extended integration times can reduce noise. For white noise:
SNR t integration
To gain four orders of magnitude in sensitivity requires
10 8
seconds of integration (approximately 3 years), which is impractical for a single measurement but may be achievable with multiple parallel experiments.

B. Resonant Mechanical Amplification

Suspending the sample on a mechanical oscillator with quality factor
Q
and driving at the resonant frequency amplifies the force signal:
F measured = Q × F actual
For
Q 10 6
(achievable in cryogenic systems), this could bring the signal within reach of current technology.

C. Differential Measurement

Using two identical samples with opposite polarizations eliminates common-mode noise:
Δ W differential = W W = 2 × Δ W sin gle
More importantly, this configuration rejects environmental gravitational noise (from tides, atmospheric pressure changes, etc.) that would otherwise dominate the error budget.

8.3. Indirect Detection: Gravitational Fractionation

An alternative approach exploits the thermodynamic consequences derived in Section 6. Rather than measuring force directly, we measure the equilibrium distribution of spin-polarized particles in a gravitational potential.

Experimental Configuration:

1.
Sample: A tall column (1–10 m) of hyperpolarized 3 He gas at low pressure (to maximize mean free path and reduce collisional depolarization).
2.
Polarization: Continuous optical pumping to maintain polarization against relaxation.
3.
Detection: NMR imaging to map the spatial distribution of polarization as a function of height.
4.
Modulation: Periodic reversal of the magnetic field direction to switch between light” ( S g ) and heavy” ( S g ) configurations.

Expected Signal:

From Equation (54), the relative density difference between top and bottom of a column of height
h
is:
Δ n n h H 0 × η × v s v c 2
For 3 He at 300 K:
  • Scale height:
    H 0 = k B T / m g 8.5 × 10 4
    m
  • For
    h = 10
    m:
    h / H 0 1.2 × 10 4
Thus:
Δ n n 1.2 × 10 4 × 10 22 × 2 × 10 7 2.4 × 10 19
This is an extraordinarily small effect, well beyond current NMR sensitivity. However, the signal scales linearly with column height, suggesting that space-based experiments with very tall columns (kilometers) might eventually be feasible.

8.4. Centrifuge Enhancement

A more promising approach uses centrifugal acceleration to enhance the effective gravitational field. In a centrifuge with angular velocity
ω c
and radius R c :
g eff = ω c 2 R c
For
ω c = 10 4
rad/s and
R c = 0.1
m:
g eff = 10 8 × 0.1 = 10 7 m / s 2 10 6 g
This increases the signal by a factor of
10 6
while also decreasing the effective scale height by the same factor, leading to a net enhancement of
10 12
in the fractional density gradient:
Δ n n 2.4 × 10 19 × 10 12 2.4 × 10 7
This approaches the sensitivity of precision NMR techniques and may be detectable with careful experimental design.

8.5. Summary of Experimental Approaches

Table 2. Comparison of Experimental Approaches
Table 2. Comparison of Experimental Approaches
Method Predicted Signal Current Sensitivity Gap
Direct gravimetry r k δ g/g ∼ 5 × 10 16 r k r k 10 12 r k r k 10 4 r k
Resonant amplification r k δ g/g ∼ 5 × 10 10 r k r k 10 12 r k Achievable
Tall column fractionation r k δ n/n ∼ 2.4 × 10 19 r k r k 10 6 r k r k 10 14 r k
Centrifuge fractionation r k δ n/n ∼ 2.4 × 10 7 r k r k 10 6 r k r k 10 2 r k
The most promising near-term approach appears to be resonant mechanical amplification combined with differential measurement.

9. Reconciliation with Existing Experimental Constraints

9.1. The Hughes-Drever Experiments

The Hughes-Drever experiments [12,13] are often cited as stringent tests of Local Lorentz Invariance (LLI) and, by extension, potential spin-gravity couplings. These experiments search for anisotropies in atomic energy levels as the Earth rotates, changing the orientation of nuclear spins relative to fixed stars.

Original Results:

Hughes et al. (1960) and Drever (1961) found no energy level splitting in 7 Li and 23 Na nuclei at the level of:
Δ E E < 10 20
This null result is sometimes interpreted as ruling out any spin-direction-dependent physics at this level.

Resolution:

We emphasize a crucial distinction between inertial mass and gravitational weight:
1.
Hughes-Drever experiments test inertial mass anisotropy. They measure whether the binding energy (and hence rest mass) of a nucleus depends on the orientation of its spin relative to some preferred frame.
2.
Our prediction concerns gravitational weight anisotropy. We predict that the coupling between mass and the gravitational field depends on spin orientation, not that the mass itself is anisotropic.
In the language of the parametrized post-Newtonian (PPN) formalism, Hughes-Drever experiments constrain violations of Local Lorentz Invariance (parameters α 1 , α 2 , etc.), while our prediction concerns the structure of the gravitational coupling itself.
Mathematically:
  • Hughes-Drever: Tests whether
    m i ( s ^ )
    depends on spin direction
  • Our prediction: Tests whether
    W = m i · g eff ( s ^ )
    depends on spin direction
These are logically independent. A theory in which
m i
is isotropic but
g eff
is anisotropic would satisfy Hughes-Drever constraints while still predicting weight anisotropy.

9.2. Torsion Balance Experiments

Torsion balance experiments, such as those performed by the Eot-Wash group [14], have tested the universality of free fall to extraordinary precision:
| a 1 a 2 | | a | < 10 13
for various material pairs.

9.2.0.14. Relevance to Our Prediction:

These experiments compare the gravitational acceleration of macroscopic test masses with random spin orientations. In unpolarized matter:
( s ^ · g ^ ) 2 = 1 3
for all materials. Thus, even if spin-gravity coupling exists, it would affect all unpolarized test masses equally and would not be detected in differential free-fall measurements.
Our prediction specifically requires polarized samples to observe anisotropy. Standard Eot-Wash experiments using unpolarized materials are not sensitive to this effect.

9.3. Neutron Gravity Experiments

Experiments measuring the gravitational quantum states of ultracold neutrons [15] have verified the gravitational acceleration of neutrons to:
δ g g < 10 3
These neutrons are unpolarized, so the same argument applies: any spin-gravity coupling is averaged to its isotropic mean value and cannot be distinguished from standard gravity.
Polarized neutron gravity experiments have been proposed [16] but not yet performed at sufficient precision to test our predictions.

9.4. Gyroscope Experiments

Gravity Probe B [17] measured the geodetic and frame-dragging precession of orbiting gyroscopes, confirming General Relativity at the level of:
δ Ω Ω < 10 3
However, these experiments measured precession rates, not weight changes. The gyroscope’s center of mass followed a geodesic regardless of spin orientation—which is consistent with our prediction that weight (deviation from geodesic) depends on spin, not the geodesic itself.

9.5. Summary of Existing Constraints

Table 3. Existing Experimental Constraints and Their Relevance
Table 3. Existing Experimental Constraints and Their Relevance
Experiment Type What It Tests Constraint Conflict?
Hughes-Drever Inertial mass isotropy < 10 20 No
Eot-Wash Universality of free fall < 10 13 No
Neutron gravity Neutron gravitational mass < 10 3 No
Gravity Probe B Gyroscope precession < 10 3 No
Conclusion: No existing experiment directly constrains the spin-orientation-dependent weight effect predicted in this paper. The effect specifically requires polarized samples and differential weight measurements—a combination not yet explored experimentally.

10. Theoretical Implications

10.1. Reinterpretation of the Weak Equivalence Principle

The Weak Equivalence Principle (WEP) states that the trajectory of a freely falling test body is independent of its internal structure and composition. Our analysis suggests a refinement:
  • Standard WEP:
    a grav
    is independent of internal structure.
  • Refined WEP:
    a grav
    is independent of internal scalar properties (mass, binding energy, composition) but may depend on internal vector properties (spin orientation) through geometric coupling.
This refinement is not a violation of General Relativity but rather a consequence of the Mathisson-Papapetrou equations, which have been part of GR since the 1930s. What we have done is to extract a potentially measurable macroscopic consequence of these equations.

10.2. The Spin-Curvature Coupling Hierarchy

Our analysis reveals a hierarchy of spin-curvature effects:
1.
Orbital Motion: Full geodesic coupling; weight can be reduced to zero (satellites).
2.
Macroscopic Rotation: Reduced coupling due to material stress limits; maximal weight reduction limited by mechanical failure.
3.
Microscopic Spin: Strongly suppressed coupling due to scale mismatch; effect survives but is reduced by factor η 1 .
This hierarchy explains why the effect has not been noticed historically: macroscopic rotation cannot reach orbital velocity, and microscopic spin, while fast enough, couples too weakly to produce obvious effects.

10.3. Connection to Quantum Gravity

The spin-curvature coupling explored here exists entirely within classical General Relativity. However, it may have implications for quantum gravity:
1.
Spin is fundamentally quantum. The “internal velocity” of a spin-1/2 particle is a semi-classical approximation. A full quantum treatment would involve the coupling between the spin operator
S ^
and the curvature tensor.
2.
Measurement implications. If weight depends on spin orientation, and spin orientation can be in superposition, then weight itself may exhibit quantum superposition. This connects to ongoing research on gravitationally-induced decoherence [18].
3.
Geometric phase effects. The coupling between spin and curvature may induce geometric (Berry) phases in gravitational contexts, analogous to the Aharonov-Bohm effect in electromagnetism.

10.4. Cosmological Implications

If spin-gravity coupling affects the motion of polarized matter, there may be cosmological consequences:
1.
Galactic magnetic fields may induce partial polarization of interstellar gas, affecting its gravitational dynamics.
2.
Neutron stars are highly magnetized and may have significant spin polarization in their crusts, potentially affecting their structure.
3.
Early universe phase transitions may have created domains of polarized matter with anomalous gravitational properties.
These speculations require further theoretical development but suggest that spin-curvature coupling may have broader astrophysical relevance.

10.5. Summary of Key Results

This paper has developed a geometric framework for understanding spin-gravity coupling based entirely on the standard apparatus of General Relativity. Our principal results are:
1.
Geodesic Levitation: A body rotating with tangential velocity
v c = G M / R
experiences zero effective weight because its constituent mass elements follow geodesics that do not intersect the ground. This is exact in Schwarzschild spacetime, with the metric factor
( 1 2 M / r )
canceling completely.
2.
Weight Reduction Formula: For velocities below critical, weight decreases parabolically:
W ( v ) = W 0 1 v 2 v c 2
3.
Microscopic Extension: Nuclear spin represents internal tangential motion at velocities far exceeding v c , but the coupling is suppressed by a geometric factor η ( r nucleus / R Earth ) n 1 .
4.
Weight Anisotropy: The suppressed coupling is maximized when spin is aligned with gravity (horizontal rotation plane) and minimized when spin is perpendicular to gravity (vertical rotation plane).
5.
Geometric Buoyancy: Spin-polarized gases should exhibit height-dependent density anomalies, with “light” (horizontally-spinning) particles enriched at higher altitudes.
6.
Formal Verification: All mathematical derivations have been machine-verified using the Lean 4 theorem prover.

10.6. Strengths of the Analysis

1.
Derivation from first principles. No new physics is postulated; all results follow from the Schwarzschild metric and geodesic equation.
2.
Exact results. The geodesic cancellation is not an approximation but an exact algebraic identity, as confirmed by formal verification.
3.
Clear experimental predictions. The theory makes specific, falsifiable predictions about the weight of polarized matter.
4.
Consistency with existing experiments. No current experimental null result directly contradicts our predictions.

10.7. Limitations and Open Questions

1.
The coupling factor
η
is not derived from first principles. We have argued on dimensional grounds that
η 1
for microscopic systems, but the precise functional form remains unknown.
2.
Quantum treatment is incomplete. Our semi-classical treatment of spin may miss important quantum effects. A full quantum field theory treatment of spin-curvature coupling is needed.
3.
Experimental detection is challenging. The predicted effects are small and require significant advances in experimental technique to detect.
4.
Astrophysical consequences are speculative. While we have suggested possible cosmological implications, these require detailed modeling to assess.

10.8. Relation to Other Work

The coupling between spin and gravitational curvature has been studied extensively in the context of the Mathisson-Papapetrou-Dixon equations [5,6,7]. Our contribution is to extract a potentially measurable consequence—weight anisotropy—from this established formalism.
Recent work on “gravitomagnetism” [19] and spin-gravity coupling in atom interferometry [20] explores related territory but focuses on different observables (precession rather than weight).
The formal verification aspect of our work connects to growing efforts to apply proof assistants to physics [21], ensuring mathematical rigor in theoretical derivations.

11. Conclusions

11.1. Principal Conclusions

We have demonstrated that General Relativity predicts a path-dependent coupling between spin angular momentum and gravitational weight. The key insights are:
1.
Gravity as slope: The geometric picture of gravity as spacetime curvature implies that tangential motion (parallel to equipotential surfaces) reduces the effective “downward” tendency.
2.
Exact cancellation: In the Schwarzschild metric, the factor
( 1 2 M / r )
appears identically in both the gravitational and centrifugal terms of the radial geodesic equation, leading to exact cancellation at the critical (Keplerian) velocity.
3.
Spin as micro-orbit: Nuclear spin represents internal tangential motion that, despite being suppressed by geometric factors, may produce measurable weight anisotropy in polarized samples.
4.
Weight is not intrinsic: The “weight” of an object—defined as the force required to deviate it from geodesic motion—depends on the geometric relationship between internal motion and external curvature.

11.2. Experimental Outlook

The most promising experimental approaches are:
1.
Resonantly-amplified differential gravimetry of hyperpolarized gases with modulated spin orientation.
2.
Centrifuge-enhanced gravitational fractionation measuring density gradients in polarized samples under high effective gravity.
3.
Polarized neutron interferometry exploiting the long coherence lengths available in neutron optics.
Detection will require pushing current experimental precision by 2–4 orders of magnitude, which appears challenging but not fundamentally impossible.

11.3. Theoretical Outlook

If spin-gravity coupling is confirmed experimentally, it would:
1.
Validate an overlooked prediction of GR present in the Mathisson-Papapetrou equations since 1937.
2.
Refine the Weak Equivalence Principle to distinguish scalar and vector internal properties.
3.
Open new avenues for quantum gravity research by connecting spin (a quantum property) to spacetime curvature (a classical geometric property).
4.
Provide new tools for precision metrology based on spin-gravity effects.

11.4. Closing Remarks

The falling apple and the orbiting moon, as Newton recognized, are manifestations of the same underlying phenomenon. What we have shown is that the rotating nucleus may also partake in this unity—not through macroscopic translation through space, but through the geometric equivalence of local spin and orbital motion in curved spacetime.
If gravity is truly geometry, then the way a particle moves through that geometry—including its internal spin—must affect its gravitational behavior. The “weight” we measure on a scale is not a fixed property of mass, but a dynamic geometric quantity reflecting the relationship between matter’s motion and spacetime’s curvature.
We conclude with a prediction and a challenge:
Prediction: Polarized atomic nuclei exhibit weight anisotropy at the level of
Δ W / W 10 15
to 10 14 , with vertically-polarized spins (horizontal rotation) being lighter than horizontally-polarized spins (vertical rotation).
Challenge: Design and execute an experiment to test this prediction at the required precision.
The answer will illuminate either a hidden consequence of General Relativity or the limitations of semiclassical reasoning about quantum spin in curved spacetime. Either outcome would advance our understanding of the deep connection between gravity and the quantum world.

Acknowledgments

The author thanks [Collaborator Names] for insightful discussions on the geometric interpretation of spin-curvature coupling and the feasibility of experimental tests. The formal verification component benefited from the Mathlib community’s development of analysis libraries for Lean 4. [Funding acknowledgment if applicable.]

Appendix A. Complete Lean 4 Source Code

The following is the complete, self-contained Lean 4 source code that formally verifies the mathematical results presented in this paper. This code compiles successfully with Lean 4 and Mathlib4.
Listing 7. Complete Lean 4 Verification Code
Preprints 196712 i007a
Preprints 196712 i007b
Preprints 196712 i007c
Preprints 196712 i007d
Preprints 196712 i007e
Preprints 196712 i007f
Preprints 196712 i007g

Appendix B. Notation and Conventions

Table A1. Symbol Definitions.
Table A1. Symbol Definitions.
Symbol Definition Units (SI)
M Gravitating mass kg
R, r Radial coordinate m
G Gravitational constant m 3 /(kg·s 2 )
c Speed of light m/s
Ω Coordinate angular velocity rad/s
v Tangential velocity m/s
v c Critical (orbital) velocity m/s
u μ Four-velocity components dimensionless
Γ α β μ Christoffel symbols m 1
S Spin angular momentum J·s
η Geometric coupling factor dimensionless
H Atmospheric scale height m
k B Boltzmann constant J/K
Throughout this paper, we use the metric signature
( , + , + , + )
and, where noted, geometric units with G = c = 1 .

Appendix C. Detailed Calculation of Christoffel Symbols

For completeness, we provide the explicit calculation of the Christoffel symbols used in this paper.
The Christoffel symbols are defined by:
Γ α β μ = 1 2 g μ ν α g ν β + β g ν α ν g α β
For the Schwarzschild metric in the equatorial plane:
d s 2 = 1 2 M r d t 2 + 1 2 M r 1 d r 2 + r 2 d ϕ 2
The non-zero metric components and their inverses are:
g t t = 1 2 M r , g t t = 1 2 M r 1
g r r = 1 2 M r 1 , g r r = 1 2 M r
g ϕ ϕ = r 2 , g ϕ ϕ = r 2
For Γ t t r :
Γ t t r = 1 2 g r r r g t t
= 1 2 1 2 M r 2 M r 2 1
= M r 2 1 2 M r
For Γ ϕ ϕ r :
Γ ϕ ϕ r = 1 2 g r r r g ϕ ϕ
= 1 2 1 2 M r 2 r
= r 1 2 M r
These expressions confirm Equations (4) and () in the main text.

References

  1. Einstein, A. Die Grundlage der allgemeinen Relativitatstheorie. Annalen der Physik 1916, 354(7), 769–822. [Google Scholar] [CrossRef]
  2. Misner, C. W.; Thorne, K. S.; Wheeler, J. A. Gravitation; W. H. Freeman, 1973. [Google Scholar]
  3. Schwarzschild, K. Uber das Gravitationsfeld eines Massenpunktes nach der Einsteinschen Theorie. Sitzungsberichte der Koniglich Preußischen Akademie der Wissenschaften 1916, 189–196. [Google Scholar]
  4. Griffiths, D. J. Introduction to Quantum Mechanics, 3rd ed.; Cambridge University Press, 2018. [Google Scholar]
  5. Mathisson, M. Neue mechanik materieller systeme. Acta Physica Polonica 1937, 6, 163–200. [Google Scholar]
  6. Papapetrou, A. Spinning test-particles in general relativity. I. Proceedings of the Royal Society of London. Series A 1951, 209(1097), 248–258. [Google Scholar] [CrossRef]
  7. Dixon, W. G. Dynamics of extended bodies in general relativity. I. Momentum and angular momentum. Proceedings of the Royal Society of London. Series A 1970, 314(1519), 499–527. [Google Scholar] [CrossRef]
  8. Reif, F. Fundamentals of Statistical and Thermal Physics; McGraw-Hill, 1965. [Google Scholar]
  9. de Moura, L.; et al. “The Lean 4 Theorem Prover and Programming Language.” Automated Deduction – CADE 28. In Lecture Notes in Computer Science; 2021; vol. 12699. [Google Scholar]
  10. Walker, T. G.; Happer, W. Spin-exchange optical pumping of noble-gas nuclei. Reviews of Modern Physics 1997, 69(2), 629–642. [Google Scholar] [CrossRef]
  11. Crossley, D.; Hinderer, J.; Riccardi, U. The measurement of surface gravity. Reports on Progress in Physics 2013, 76(4), 046101. [Google Scholar] [CrossRef] [PubMed]
  12. Hughes, V. W.; Robinson, H. G.; Beltran-Lopez, V. Upper Limit for the Anisotropy of Inertial Mass from Nuclear Resonance Experiments. Physical Review Letters 1960, 4(7), 342–344. [Google Scholar] [CrossRef]
  13. Drever, R. W. P. A search for anisotropy of inertial mass using a free precession technique. Philosophical Magazine 1961, 6(65), 683–687. [Google Scholar] [CrossRef]
  14. Schlamminger, S.; et al. Test of the equivalence principle using a rotating torsion balance. Physical Review Letters 2008, 100(4), 041101. [Google Scholar] [CrossRef] [PubMed]
  15. Nesvizhevsky, V. V.; et al. Quantum states of neutrons in the Earth’s gravitational field. Nature 2002, 415(6869), 297–299. [Google Scholar] [CrossRef] [PubMed]
  16. Abele, H.; Leeb, H. Gravitation and quantum interference experiments with neutrons. New Journal of Physics 2012, 14(5), 055010. [Google Scholar] [CrossRef]
  17. Everitt, C. W. F.; et al. Gravity Probe B: Final results of a space experiment to test general relativity. Physical Review Letters 2011, 106(22), 221101. [Google Scholar] [CrossRef] [PubMed]
  18. Penrose, R. On gravity’s role in quantum state reduction. General Relativity and Gravitation 1996, 28(5), 581–600. [Google Scholar] [CrossRef]
  19. Mashhoon, B. Gravitoelectromagnetism: A brief review. The Measurement of Gravitomagnetism: A Challenging Enterprise 2003, 29–39. [Google Scholar]
  20. Tarallo, M. G.; et al. Test of Einstein equivalence principle for 0-spin and half-integer-spin atoms: Search for spin-gravity coupling effects. Physical Review Letters 2014, 113(2), 023005. [Google Scholar] [CrossRef] [PubMed]
  21. Buzzard, K. Formalising mathematics in simple type theory. In Proceedings of the 13th International Conference on Intelligent Computer Mathematics, 2020. [Google Scholar]
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.
Copyright: This open access article is published under a Creative Commons CC BY 4.0 license, which permit the free download, distribution, and reuse, provided that the author and preprint are cited in any reuse.
Prerpints.org logo

Preprints.org is a free preprint server supported by MDPI in Basel, Switzerland.

Subscribe

Disclaimer

Terms of Use

Privacy Policy

Privacy Settings

© 2026 MDPI (Basel, Switzerland) unless otherwise stated