<?xml version="1.0" encoding="UTF-8"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD JATS (Z39.96) Journal Publishing DTD v1.4 20241031//EN" "JATS-journalpublishing1-4.dtd">
<article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" article-type="research-article" dtd-version="1.4" xml:lang="en">
  <front>
    <journal-meta>
      <journal-id journal-id-type="publisher-id">jhepgc</journal-id>
      <journal-title-group>
        <journal-title>Journal of High Energy Physics, Gravitation and Cosmology</journal-title>
      </journal-title-group>
      <issn pub-type="epub">2380-4335</issn>
      <issn pub-type="ppub">2380-4327</issn>
      <publisher>
        <publisher-name>Scientific Research Publishing</publisher-name>
      </publisher>
    </journal-meta>
    <article-meta>
      <article-id pub-id-type="doi">10.4236/jhepgc.2026.123090</article-id>
      <article-id pub-id-type="publisher-id">jhepgc-152722</article-id>
      <article-categories>
        <subj-group>
          <subject>Article</subject>
        </subj-group>
        <subj-group>
          <subject>Physics</subject>
          <subject>Mathematics</subject>
        </subj-group>
      </article-categories>
      <title-group>
        <article-title>Quantum Vacuum Vortex Dynamics and Time Dilation: A Phenomenological Density-Perturbation Interpretation of the Hafele-Keating Experiment</article-title>
      </title-group>
      <contrib-group>
        <contrib contrib-type="author">
          <contrib-id contrib-id-type="orcid">0000-0001-9372-9113</contrib-id>
          <name name-style="western">
            <surname>Butto</surname>
            <given-names>Nader</given-names>
          </name>
          <xref ref-type="aff" rid="aff1">1</xref>
        </contrib>
      </contrib-group>
      <aff id="aff1"><label>1</label> Independent Researcher, Petah Tikva, Israel </aff>
      <author-notes>
        <fn fn-type="conflict" id="fn-conflict">
          <p>The author declares no conflicts of interest regarding the publication of this paper.</p>
        </fn>
      </author-notes>
      <pub-date pub-type="epub">
        <day>03</day>
        <month>06</month>
        <year>2026</year>
      </pub-date>
      <pub-date pub-type="collection">
        <month>06</month>
        <year>2026</year>
      </pub-date>
      <volume>12</volume>
      <issue>03</issue>
      <fpage>1759</fpage>
      <lpage>1784</lpage>
      <history>
        <date date-type="received">
          <day>29</day>
          <month>12</month>
          <year>2025</year>
        </date>
        <date date-type="accepted">
          <day>20</day>
          <month>07</month>
          <year>2026</year>
        </date>
        <date date-type="published">
          <day>23</day>
          <month>07</month>
          <year>2026</year>
        </date>
      </history>
      <permissions>
        <copyright-statement>© 2026 by the authors and Scientific Research Publishing Inc.</copyright-statement>
        <copyright-year>2026</copyright-year>
        <license license-type="open-access">
          <license-p> This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license ( <ext-link ext-link-type="uri" xlink:href="https://creativecommons.org/licenses/by/4.0/">https://creativecommons.org/licenses/by/4.0/</ext-link> ). </license-p>
        </license>
      </permissions>
      <self-uri content-type="doi" xlink:href="https://doi.org/10.4236/jhepgc.2026.123090">https://doi.org/10.4236/jhepgc.2026.123090</self-uri>
      <abstract>
        <p>Relativity accurately predicts the time shifts observed in the Hafele-Keating experiment through the combined effects of kinematic and gravitational time dilation. The present article proposes a complementary physical interpretation in which the same weak-field clock-shift structure is derived from small perturbations in an Earth-coupled structured quantum vacuum. In this framework, the Earth is assumed to be dynamically coupled to a large-scale vacuum vortex, and an airborne atomic clock experiences two independent perturbations: a directional kinetic perturbation caused by motion relative to the rotating vacuum flow, and a radial rarefaction caused by altitude. The vacuum-density perturbation is therefore derived from flight parameters, including Earth’s angular velocity, latitude, aircraft velocity, altitude, and flight duration, rather than being inferred from the observed clock shifts. The perturbation is then connected to atomic clock-rate variation through a wave-speed relation in which the effective electromagnetic propagation speed depends on the mechanically active vacuum-density scale. Using physically reasonable averaged flight parameters, the model reproduces the sign, order, and approximate magnitude of the eastward and westward Hafele-Keating clock shifts. The proposed framework does not challenge the empirical validity of special or general relativity, but offers a possible physical substrate interpretation of their weak-field clock-rate effects. The article concludes by identifying limitations of the present phenomenological formulation and proposing future tests involving latitude dependence, polar routes, altitude variation, counter-rotating satellites, and modern optical-clock experiments.</p>
      </abstract>
      <kwd-group kwd-group-type="author-generated" xml:lang="en">
        <kwd>Time Dilation</kwd>
        <kwd>Hafele-Keating Experiment</kwd>
        <kwd>Quantum Vacuum Vortex</kwd>
        <kwd>Vacuum Density</kwd>
        <kwd>Structured Vacuum</kwd>
        <kwd>Emergent Gravity</kwd>
        <kwd>Atomic Clocks</kwd>
        <kwd>Weak-Field Relativity</kwd>
      </kwd-group>
    </article-meta>
  </front>
  <body>
    <sec id="sec1">
      <title>1. Introduction</title>
      <p>Time dilation is one of the most precisely confirmed consequences of modern physics. Since Einstein’s formulation of special relativity in 1905 and general relativity in 1915, differences in clock rates have been understood as arising from relative motion and gravitational potential [<xref ref-type="bibr" rid="B1">1</xref>]. In special relativity, a moving clock ticks more slowly with respect to an inertial observer; in general relativity, a clock located at a higher gravitational potential ticks faster than one closer to a gravitating body. These effects are not merely theoretical. They are routinely confirmed in particle physics, satellite navigation, atomic-clock comparisons, and precision metrology. The Global Positioning System provides a particularly important practical example, since its operation requires relativistic corrections for both gravitational frequency shifts and motional time dilation of satellite and ground clocks.</p>
      <p>One of the historically important experimental confirmations of relativistic time dilation was the Hafele-Keating experiment, performed in 1971 and published in 1972 [<xref ref-type="bibr" rid="B2">2</xref>]. In this experiment, cesium-beam atomic clocks were flown eastward and westward around the Earth aboard commercial aircraft and then compared with reference clocks at the U.S. Naval Observatory. The eastward clocks lost time, whereas the westward clocks gained time relative to the ground reference clocks. The result was successfully explained by combining two effects: kinematic time dilation due to motion relative to the Earth-centered inertial frame, and gravitational time dilation due to the higher altitude of the aircraft.</p>
      <p>The standard relativistic calculation remains experimentally successful and is not challenged in the present work. However, relativity describes time dilation primarily through spacetime geometry and time-coordinate transformations. This description is mathematically powerful, yet it leaves open a deeper interpretive question: whether the geometric description may correspond to an underlying physical process occurring in the vacuum itself. In other words, can the same first-order clock-shift structure be interpreted as the manifestation of a structured physical substrate rather than only as a coordinate effect?</p>
      <p>This question has motivated several lines of theoretical research. Earlier approaches included Lorentz-type ether theories, pilot-wave interpretations, stochastic electrodynamics, and zero-point-field models of inertia [<xref ref-type="bibr" rid="B3">3</xref>]-[<xref ref-type="bibr" rid="B7">7</xref>]. More recent developments in analog gravity, emergent gravity, superfluid vacuum theory, and quantum fields in curved spacetime have suggested that spacetime behavior may, at least in some regimes, be understood as emerging from deeper vacuum or medium-like degrees of freedom [<xref ref-type="bibr" rid="B8">8</xref>]-[<xref ref-type="bibr" rid="B10">10</xref>]. These approaches do not necessarily contradict relativity; rather, they attempt to provide a physical interpretation for structures that relativity describes geometrically.</p>
      <p>The present article develops such a complementary interpretation using a quantum-vacuum vortex framework. In this model, the vacuum is treated as a structured, superfluid-like medium capable of sustaining large-scale vortex configurations. The Earth is assumed to be dynamically coupled to an Earth-centered vacuum vortex. A clock at the surface of the Earth co-rotates within this structured vacuum flow, while an aircraft clock moving eastward or westward experiences a different effective velocity relative to the vacuum vortex. This difference in effective velocity produces a small kinetic perturbation in the local vacuum-density state.</p>
      <p>A second contribution arises from altitude. In the proposed model, radial position within the Earth-centered vortex is associated with a small variation in effective vacuum density. At higher altitude, the aircraft clock experiences a slightly reduced vacuum-density state relative to the ground reference clock. Thus, the total perturbation affecting the airborne clock has two components: a directional kinetic component due to motion relative to the rotating vacuum vortex, and a radial component due to altitude.</p>
      <p>The central purpose of this article is to show that the vacuum-density perturbation can be derived independently from vortex kinematics and altitude-dependent radial rarefaction, without using the observed Hafele-Keating clock shifts as input. The derived perturbation is then connected to clock-rate variation through a wave-speed relation. The effective local propagation speed of electromagnetic interactions is assumed to depend on the mechanically active vacuum-density scale according to</p>
      <disp-formula id="FD1">
        <label>(1)</label>
        <mml:math>
          <mml:mrow>
            <mml:msub>
              <mml:mi>c</mml:mi>
              <mml:mrow>
                <mml:mtext>eff</mml:mtext>
              </mml:mrow>
            </mml:msub>
            <mml:mo>=</mml:mo>
            <mml:mi>c</mml:mi>
            <mml:msup>
              <mml:mrow>
                <mml:mrow>
                  <mml:mo>(</mml:mo>
                  <mml:mrow>
                    <mml:mn>1</mml:mn>
                    <mml:mo>+</mml:mo>
                    <mml:mfrac>
                      <mml:mrow>
                        <mml:mi>Δ</mml:mi>
                        <mml:mi>ρ</mml:mi>
                      </mml:mrow>
                      <mml:mrow>
                        <mml:msub>
                          <mml:mi>ρ</mml:mi>
                          <mml:mrow>
                            <mml:mtext>eff</mml:mtext>
                          </mml:mrow>
                        </mml:msub>
                      </mml:mrow>
                    </mml:mfrac>
                  </mml:mrow>
                  <mml:mo>)</mml:mo>
                </mml:mrow>
              </mml:mrow>
              <mml:mrow>
                <mml:mo>−</mml:mo>
                <mml:mrow>
                  <mml:mn>1</mml:mn>
                  <mml:mo>/</mml:mo>
                  <mml:mn>2</mml:mn>
                </mml:mrow>
              </mml:mrow>
            </mml:msup>
            <mml:mo>.</mml:mo>
          </mml:mrow>
        </mml:math>
      </disp-formula>
      <p>where <inline-formula><mml:math><mml:mrow><mml:msub><mml:mi> c </mml:mi><mml:mrow><mml:mtext> eff </mml:mtext></mml:mrow></mml:msub></mml:mrow></mml:math></inline-formula> is the effective local propagation speed of electromagnetic interactions, <inline-formula><mml:math><mml:mi> c </mml:mi></mml:math></inline-formula> is the reference speed of light, <inline-formula><mml:math><mml:mrow><mml:mi> Δ </mml:mi><mml:mi> ρ </mml:mi></mml:mrow></mml:math></inline-formula> is the local vacuum-density perturbation, and <inline-formula><mml:math><mml:mrow><mml:msub><mml:mi> ρ </mml:mi><mml:mrow><mml:mtext> eff </mml:mtext></mml:mrow></mml:msub></mml:mrow></mml:math></inline-formula> is the effective mechanically active vacuum-density scale.</p>
      <p>For small perturbations, where</p>
      <disp-formula id="FD2">
        <label>(2)</label>
        <mml:math>
          <mml:mrow>
            <mml:mrow>
              <mml:mo>|</mml:mo>
              <mml:mrow>
                <mml:mfrac>
                  <mml:mrow>
                    <mml:mi>Δ</mml:mi>
                    <mml:mi>ρ</mml:mi>
                  </mml:mrow>
                  <mml:mrow>
                    <mml:msub>
                      <mml:mi>ρ</mml:mi>
                      <mml:mrow>
                        <mml:mtext>eff</mml:mtext>
                      </mml:mrow>
                    </mml:msub>
                  </mml:mrow>
                </mml:mfrac>
              </mml:mrow>
              <mml:mo>|</mml:mo>
            </mml:mrow>
            <mml:mo>≪</mml:mo>
            <mml:mn>1.</mml:mn>
          </mml:mrow>
        </mml:math>
      </disp-formula>
      <p>the binomial approximation gives</p>
      <disp-formula id="FD3">
        <label>(3)</label>
        <mml:math>
          <mml:mrow>
            <mml:msub>
              <mml:mi>c</mml:mi>
              <mml:mrow>
                <mml:mtext>eff</mml:mtext>
              </mml:mrow>
            </mml:msub>
            <mml:mo>≃</mml:mo>
            <mml:mi>c</mml:mi>
            <mml:mrow>
              <mml:mo>(</mml:mo>
              <mml:mrow>
                <mml:mn>1</mml:mn>
                <mml:mo>−</mml:mo>
                <mml:mfrac>
                  <mml:mn>1</mml:mn>
                  <mml:mn>2</mml:mn>
                </mml:mfrac>
                <mml:mfrac>
                  <mml:mrow>
                    <mml:mi>Δ</mml:mi>
                    <mml:mi>ρ</mml:mi>
                  </mml:mrow>
                  <mml:mrow>
                    <mml:msub>
                      <mml:mi>ρ</mml:mi>
                      <mml:mrow>
                        <mml:mtext>eff</mml:mtext>
                      </mml:mrow>
                    </mml:msub>
                  </mml:mrow>
                </mml:mfrac>
              </mml:mrow>
              <mml:mo>)</mml:mo>
            </mml:mrow>
            <mml:mo>.</mml:mo>
          </mml:mrow>
        </mml:math>
      </disp-formula>
      <p>Since atomic transition rates are assumed to scale with the effective electromagnetic propagation speed, the fractional clock-rate variation becomes</p>
      <disp-formula id="FD4">
        <label>(4)</label>
        <mml:math>
          <mml:mrow>
            <mml:mfrac>
              <mml:mrow>
                <mml:mi>Δ</mml:mi>
                <mml:mi>τ</mml:mi>
              </mml:mrow>
              <mml:mi>T</mml:mi>
            </mml:mfrac>
            <mml:mo>≃</mml:mo>
            <mml:mo>−</mml:mo>
            <mml:mfrac>
              <mml:mn>1</mml:mn>
              <mml:mn>2</mml:mn>
            </mml:mfrac>
            <mml:mfrac>
              <mml:mrow>
                <mml:mi>Δ</mml:mi>
                <mml:mi>ρ</mml:mi>
              </mml:mrow>
              <mml:mrow>
                <mml:msub>
                  <mml:mi>ρ</mml:mi>
                  <mml:mrow>
                    <mml:mtext>eff</mml:mtext>
                  </mml:mrow>
                </mml:msub>
              </mml:mrow>
            </mml:mfrac>
            <mml:mo>.</mml:mo>
          </mml:mrow>
        </mml:math>
      </disp-formula>
      <p>This equation expresses the proposed link between vacuum-density variation and atomic clock-rate variation. The factor 1/2 arises from the square-root dependence of wave speed on density, as in ordinary wave propagation through an elastic medium.</p>
      <p>The vacuum-density perturbation itself is derived independently from vortex dynamics. The tangential velocity of the ground reference clock in an Earth-centered inertial frame is</p>
      <disp-formula id="FD5">
        <label>(5)</label>
        <mml:math>
          <mml:mrow>
            <mml:msub>
              <mml:mi>v</mml:mi>
              <mml:mo>⊕</mml:mo>
            </mml:msub>
            <mml:mrow>
              <mml:mo>(</mml:mo>
              <mml:mi>λ</mml:mi>
              <mml:mo>)</mml:mo>
            </mml:mrow>
            <mml:mo>=</mml:mo>
            <mml:msub>
              <mml:mi>Ω</mml:mi>
              <mml:mo>⊕</mml:mo>
            </mml:msub>
            <mml:msub>
              <mml:mi>R</mml:mi>
              <mml:mo>⊕</mml:mo>
            </mml:msub>
            <mml:mi>cos</mml:mi>
            <mml:mi>λ</mml:mi>
            <mml:mo>.</mml:mo>
          </mml:mrow>
        </mml:math>
      </disp-formula>
      <p>For an aircraft moving approximately eastward or westward with ground-relative speed <inline-formula><mml:math><mml:mi> u </mml:mi></mml:math></inline-formula> , its inertial-frame velocity is</p>
      <disp-formula id="FD6">
        <label>(6)</label>
        <mml:math>
          <mml:mrow>
            <mml:msub>
              <mml:mi>v</mml:mi>
              <mml:mrow>
                <mml:mtext>air</mml:mtext>
              </mml:mrow>
            </mml:msub>
            <mml:mo>=</mml:mo>
            <mml:msub>
              <mml:mi>v</mml:mi>
              <mml:mo>⊕</mml:mo>
            </mml:msub>
            <mml:mo>±</mml:mo>
            <mml:mi>u</mml:mi>
            <mml:mo>.</mml:mo>
          </mml:mrow>
        </mml:math>
      </disp-formula>
      <p>where the plus sign applies to eastward motion and the minus sign to westward motion.</p>
      <p>To reproduce the weak-field airborne-ground clock comparison, the structured-vacuum model postulates the following first-order phenomenological relation between the effective vacuum-density perturbation and the differences in inertial velocity and gravitational potential:</p>
      <disp-formula id="FD7">
        <label>(7)</label>
        <mml:math>
          <mml:mrow>
            <mml:mfrac>
              <mml:mrow>
                <mml:mtext>Δ</mml:mtext>
                <mml:mi>ρ</mml:mi>
              </mml:mrow>
              <mml:mrow>
                <mml:msub>
                  <mml:mi>ρ</mml:mi>
                  <mml:mrow>
                    <mml:mtext>eff</mml:mtext>
                  </mml:mrow>
                </mml:msub>
              </mml:mrow>
            </mml:mfrac>
            <mml:mo>=</mml:mo>
            <mml:mfrac>
              <mml:mrow>
                <mml:msubsup>
                  <mml:mi>v</mml:mi>
                  <mml:mrow>
                    <mml:mtext>air</mml:mtext>
                  </mml:mrow>
                  <mml:mn>2</mml:mn>
                </mml:msubsup>
                <mml:mo>−</mml:mo>
                <mml:msubsup>
                  <mml:mi>v</mml:mi>
                  <mml:mo>⊕</mml:mo>
                  <mml:mn>2</mml:mn>
                </mml:msubsup>
              </mml:mrow>
              <mml:mrow>
                <mml:msup>
                  <mml:mi>c</mml:mi>
                  <mml:mn>2</mml:mn>
                </mml:msup>
              </mml:mrow>
            </mml:mfrac>
            <mml:mo>−</mml:mo>
            <mml:mfrac>
              <mml:mrow>
                <mml:mn>2</mml:mn>
                <mml:mtext>ΔΦ</mml:mtext>
              </mml:mrow>
              <mml:mrow>
                <mml:msup>
                  <mml:mi>c</mml:mi>
                  <mml:mn>2</mml:mn>
                </mml:msup>
              </mml:mrow>
            </mml:mfrac>
          </mml:mrow>
        </mml:math>
      </disp-formula>
      <p>For a nearly constant gravitational acceleration and an altitude h ≪ R<sub>⊕</sub>,</p>
      <disp-formula id="FD8">
        <mml:math>
          <mml:mrow>
            <mml:mtext>ΔΦ</mml:mtext>
            <mml:mo>≃</mml:mo>
            <mml:mi>g</mml:mi>
            <mml:mi>h</mml:mi>
            <mml:mo>,</mml:mo>
          </mml:mrow>
        </mml:math>
      </disp-formula>
      <p>so that</p>
      <disp-formula id="FD9">
        <mml:math>
          <mml:mrow>
            <mml:mfrac>
              <mml:mrow>
                <mml:mtext>Δ</mml:mtext>
                <mml:mi>ρ</mml:mi>
              </mml:mrow>
              <mml:mrow>
                <mml:msub>
                  <mml:mi>ρ</mml:mi>
                  <mml:mrow>
                    <mml:mtext>eff</mml:mtext>
                  </mml:mrow>
                </mml:msub>
              </mml:mrow>
            </mml:mfrac>
            <mml:mo>=</mml:mo>
            <mml:mfrac>
              <mml:mrow>
                <mml:msubsup>
                  <mml:mi>v</mml:mi>
                  <mml:mrow>
                    <mml:mtext>air</mml:mtext>
                  </mml:mrow>
                  <mml:mn>2</mml:mn>
                </mml:msubsup>
                <mml:mo>−</mml:mo>
                <mml:msubsup>
                  <mml:mi>v</mml:mi>
                  <mml:mo>⊕</mml:mo>
                  <mml:mn>2</mml:mn>
                </mml:msubsup>
              </mml:mrow>
              <mml:mrow>
                <mml:msup>
                  <mml:mi>c</mml:mi>
                  <mml:mn>2</mml:mn>
                </mml:msup>
              </mml:mrow>
            </mml:mfrac>
            <mml:mo>−</mml:mo>
            <mml:mfrac>
              <mml:mrow>
                <mml:mn>2</mml:mn>
                <mml:mi>g</mml:mi>
                <mml:mi>h</mml:mi>
              </mml:mrow>
              <mml:mrow>
                <mml:msup>
                  <mml:mi>c</mml:mi>
                  <mml:mn>2</mml:mn>
                </mml:msup>
              </mml:mrow>
            </mml:mfrac>
            <mml:mo>.</mml:mo>
          </mml:mrow>
        </mml:math>
      </disp-formula>
      <p>where <inline-formula><mml:math><mml:mi> g </mml:mi></mml:math></inline-formula> is the gravitational acceleration near the Earth’s surface and <inline-formula><mml:math><mml:mi> h </mml:mi></mml:math></inline-formula> is the mean flight altitude. The first term represents the kinetic density perturbation caused by motion relative to the rotating vacuum vortex. The second term represents altitude-related radial rarefaction of the vacuum-density state.</p>
      <p>Combining Equation (4) with Equation (7) gives the first-order clock-shift equation</p>
      <disp-formula id="FD10">
        <label>(8)</label>
        <mml:math display="inline">
          <mml:mrow>
            <mml:mi>Δ</mml:mi>
            <mml:mi>τ</mml:mi>
            <mml:mo>=</mml:mo>
            <mml:mi>T</mml:mi>
            <mml:mrow>
              <mml:mo>[</mml:mo>
              <mml:mrow>
                <mml:mfrac>
                  <mml:mrow>
                    <mml:mi>g</mml:mi>
                    <mml:mi>h</mml:mi>
                  </mml:mrow>
                  <mml:mrow>
                    <mml:msup>
                      <mml:mi>c</mml:mi>
                      <mml:mn>2</mml:mn>
                    </mml:msup>
                  </mml:mrow>
                </mml:mfrac>
                <mml:mo>−</mml:mo>
                <mml:mfrac>
                  <mml:mrow>
                    <mml:msubsup>
                      <mml:mi>v</mml:mi>
                      <mml:mrow>
                        <mml:mtext>eff</mml:mtext>
                      </mml:mrow>
                      <mml:mn>2</mml:mn>
                    </mml:msubsup>
                    <mml:mo>−</mml:mo>
                    <mml:msubsup>
                      <mml:mi>v</mml:mi>
                      <mml:mo>⊕</mml:mo>
                      <mml:mn>2</mml:mn>
                    </mml:msubsup>
                  </mml:mrow>
                  <mml:mrow>
                    <mml:mn>2</mml:mn>
                    <mml:msup>
                      <mml:mi>c</mml:mi>
                      <mml:mn>2</mml:mn>
                    </mml:msup>
                  </mml:mrow>
                </mml:mfrac>
              </mml:mrow>
              <mml:mo>]</mml:mo>
            </mml:mrow>
            <mml:mo>.</mml:mo>
          </mml:mrow>
        </mml:math>
      </disp-formula>
      <p>Equation (8) has the same mathematical structure as the weak-field relativistic expression for combined gravitational and kinematic clock shifts. However, its physical interpretation is different. In the standard relativistic account, the two contributions are described as gravitational and kinematic time dilation. In the present framework, the same mathematical structure is interpreted as the combined effect of radial vacuum rarefaction and kinetic vacuum-density perturbation.</p>
      <p>The article therefore does not attempt to replace the standard relativistic prediction. Rather, it aims to reinterpret the weak-field relativistic clock-shift formula in terms of structured vacuum dynamics. This distinction is important. If the vacuum-density perturbation were obtained by first inserting the observed Hafele-Keating time shifts, the model would be circular. The present work avoids this by deriving the perturbation from physical flight parameters: Earth’s angular velocity, Earth’s radius, latitude, aircraft velocity, altitude, and flight duration. The Hafele-Keating results are then used only as an experimental comparison.</p>
      <p>The article is organized as follows. Section 2 reviews the Hafele-Keating experiment and the standard relativistic interpretation. Section 3 introduces the structured quantum-vacuum vortex hypothesis. Section 4 develops the relation between vacuum density, electromagnetic propagation speed, and clock rate. Section 5 derives the vacuum-density perturbation independently from vortex dynamics and altitude. Section 6 applies the model numerically to eastward and westward flight conditions. Section 7 compares the predicted values with the Hafele-Keating measurements. Section 8 discusses the relationship between the proposed interpretation and standard relativity. Section 9 presents limitations, objections, and future experimental tests.</p>
    </sec>
    <sec id="sec2">
      <title>2. The Hafele-Keating Experiment and the Standard Relativistic Interpretation</title>
      <p>The Hafele-Keating experiment remains one of the most important direct demonstrations of relativistic clock-rate differences under real transport conditions. In 1971, four cesium-beam atomic clocks were flown around the Earth aboard commercial aircraft, first eastward and then westward, and were compared after the flights with reference clocks maintained at the U.S. Naval Observatory. The experiment was designed to test the combined effects of special-relativistic kinematic time dilation and general-relativistic gravitational time dilation in a non-laboratory setting.</p>
      <p>The essential empirical result was directional. The eastward clocks, traveling in the direction of Earth rotation, lost time relative to the ground reference clocks. The westward clocks, traveling against Earth rotation, gained time. Hafele and Keating reported an observed eastward time change of approximately −59 ns and an observed westward time change of approximately +273 ns. The sign difference is crucial, because it reflects the combined influence of altitude and direction of motion relative to an Earth-centered inertial frame.</p>
      <p>In the standard relativistic account, the total clock shift is written as the sum of a gravitational contribution and a kinematic contribution:</p>
      <disp-formula id="FD11">
        <label>(9)</label>
        <mml:math>
          <mml:mrow>
            <mml:mi>Δ</mml:mi>
            <mml:msub>
              <mml:mi>τ</mml:mi>
              <mml:mrow>
                <mml:mtext>total</mml:mtext>
              </mml:mrow>
            </mml:msub>
            <mml:mo>=</mml:mo>
            <mml:mi>Δ</mml:mi>
            <mml:msub>
              <mml:mi>τ</mml:mi>
              <mml:mrow>
                <mml:mtext>grav</mml:mtext>
              </mml:mrow>
            </mml:msub>
            <mml:mo>+</mml:mo>
            <mml:mi>Δ</mml:mi>
            <mml:msub>
              <mml:mi>τ</mml:mi>
              <mml:mrow>
                <mml:mtext>kin</mml:mtext>
              </mml:mrow>
            </mml:msub>
            <mml:mo>.</mml:mo>
          </mml:mrow>
        </mml:math>
      </disp-formula>
      <p>The gravitational term arises because the aircraft clocks travel at a higher altitude than the ground reference clocks. In the weak-field approximation, for a mean altitude h much smaller than Earth’s radius, the fractional gravitational clock shift is</p>
      <disp-formula id="FD12">
        <label>(10)</label>
        <mml:math>
          <mml:mrow>
            <mml:mfrac>
              <mml:mrow>
                <mml:mi>Δ</mml:mi>
                <mml:msub>
                  <mml:mi>τ</mml:mi>
                  <mml:mrow>
                    <mml:mtext>grav</mml:mtext>
                  </mml:mrow>
                </mml:msub>
              </mml:mrow>
              <mml:mi>T</mml:mi>
            </mml:mfrac>
            <mml:mo>≃</mml:mo>
            <mml:mfrac>
              <mml:mrow>
                <mml:mi>g</mml:mi>
                <mml:mi>h</mml:mi>
              </mml:mrow>
              <mml:mrow>
                <mml:msup>
                  <mml:mi>c</mml:mi>
                  <mml:mn>2</mml:mn>
                </mml:msup>
              </mml:mrow>
            </mml:mfrac>
            <mml:mo>.</mml:mo>
          </mml:mrow>
        </mml:math>
      </disp-formula>
      <p>where <inline-formula><mml:math display="inline"><mml:mi> g </mml:mi></mml:math></inline-formula> is the gravitational acceleration near Earth’s surface, <inline-formula><mml:math display="inline"><mml:mi> h </mml:mi></mml:math></inline-formula> is the mean flight altitude, <inline-formula><mml:math display="inline"><mml:mi> c </mml:mi></mml:math></inline-formula> is the speed of light, and <inline-formula><mml:math display="inline"><mml:mi> T </mml:mi></mml:math></inline-formula> is the total elapsed flight time. This term is positive: at higher altitude, the clock experiences a weaker gravitational potential and therefore ticks faster than a ground clock.</p>
      <p>The kinematic term arises from motion. To first order in <inline-formula><mml:math display="inline"><mml:mrow><mml:mrow><mml:mi> v </mml:mi><mml:mo> / </mml:mo><mml:mi> c </mml:mi></mml:mrow></mml:mrow></mml:math></inline-formula> , the special-relativistic contribution is</p>
      <disp-formula id="FD13">
        <label>(11)</label>
        <mml:math>
          <mml:mrow>
            <mml:mfrac>
              <mml:mrow>
                <mml:mi>Δ</mml:mi>
                <mml:msub>
                  <mml:mi>τ</mml:mi>
                  <mml:mrow>
                    <mml:mtext>kin</mml:mtext>
                  </mml:mrow>
                </mml:msub>
              </mml:mrow>
              <mml:mi>T</mml:mi>
            </mml:mfrac>
            <mml:mo>≃</mml:mo>
            <mml:mo>−</mml:mo>
            <mml:mfrac>
              <mml:mrow>
                <mml:msup>
                  <mml:mi>v</mml:mi>
                  <mml:mn>2</mml:mn>
                </mml:msup>
              </mml:mrow>
              <mml:mrow>
                <mml:mn>2</mml:mn>
                <mml:msup>
                  <mml:mi>c</mml:mi>
                  <mml:mn>2</mml:mn>
                </mml:msup>
              </mml:mrow>
            </mml:mfrac>
            <mml:mo>.</mml:mo>
          </mml:mrow>
        </mml:math>
      </disp-formula>
      <p>where <inline-formula><mml:math display="inline"><mml:mi> v </mml:mi></mml:math></inline-formula> is the velocity of the moving clock relative to the relevant inertial frame. This term is negative: the greater the velocity, the slower the clock ticks relative to the reference frame.</p>
      <p>The first term represents the postulated velocity-dependent change in the effective vacuum state relative to the ground reference clock, while the second represents the postulated decrease associated with the higher gravitational potential of the aircraft. At this stage, both identifications are phenomenological assumptions rather than results derived from a microscopic vacuum-vortex field equation.</p>
      <p>For aircraft moving around the rotating Earth, the relevant velocity is not simply the speed of the aircraft relative to the surrounding air or to the Earth’s surface. Instead, the appropriate velocity is that measured in an <bold>Earth-Centered Inertial (ECI)</bold> reference frame. A clock fixed to the Earth’s surface already possesses a tangential velocity due to Earth’s rotation. At latitude <italic>λ</italic>, this velocity is</p>
      <disp-formula id="FD14">
        <label>(12)</label>
        <mml:math>
          <mml:mrow>
            <mml:msub>
              <mml:mi>v</mml:mi>
              <mml:mo>⊕</mml:mo>
            </mml:msub>
            <mml:mrow>
              <mml:mo>(</mml:mo>
              <mml:mi>λ</mml:mi>
              <mml:mo>)</mml:mo>
            </mml:mrow>
            <mml:mo>=</mml:mo>
            <mml:msub>
              <mml:mtext>Ω</mml:mtext>
              <mml:mo>⊕</mml:mo>
            </mml:msub>
            <mml:msub>
              <mml:mi>R</mml:mi>
              <mml:mo>⊕</mml:mo>
            </mml:msub>
            <mml:mtext>cos</mml:mtext>
            <mml:mi>λ</mml:mi>
            <mml:mo>,</mml:mo>
          </mml:mrow>
        </mml:math>
      </disp-formula>
      <p>where <inline-formula><mml:math><mml:mrow><mml:msub><mml:mtext> Ω </mml:mtext><mml:mo> ⊕ </mml:mo></mml:msub></mml:mrow></mml:math></inline-formula> is Earth’s angular velocity, <inline-formula><mml:math><mml:mrow><mml:msub><mml:mi> R </mml:mi><mml:mo> ⊕ </mml:mo></mml:msub></mml:mrow></mml:math></inline-formula> is Earth’s mean radius, and <inline-formula><mml:math><mml:mi> λ </mml:mi></mml:math></inline-formula> is the geographic latitude.</p>
      <p>If an aircraft moves with ground-relative speed <italic>u</italic> along an approximately east-west trajectory, its inertial-frame velocity is obtained by adding its velocity relative to the Earth's surface to the rotational velocity of the ground reference clock. For eastward flight,</p>
      <disp-formula id="FD15">
        <label>(13)</label>
        <mml:math>
          <mml:mrow>
            <mml:msub>
              <mml:mi>v</mml:mi>
              <mml:mrow>
                <mml:mtext>air</mml:mtext>
                <mml:mo>,</mml:mo>
                <mml:mtext>east</mml:mtext>
              </mml:mrow>
            </mml:msub>
            <mml:mo>=</mml:mo>
            <mml:msub>
              <mml:mi>v</mml:mi>
              <mml:mo>⊕</mml:mo>
            </mml:msub>
            <mml:mo>+</mml:mo>
            <mml:mi>u</mml:mi>
            <mml:mo>.</mml:mo>
          </mml:mrow>
        </mml:math>
      </disp-formula>
      <p>whereas for westward flight, </p>
      <disp-formula id="FD16">
        <label>(14)</label>
        <mml:math>
          <mml:mrow>
            <mml:msub>
              <mml:mi>v</mml:mi>
              <mml:mrow>
                <mml:mtext>air</mml:mtext>
                <mml:mo>,</mml:mo>
                <mml:mtext>west</mml:mtext>
              </mml:mrow>
            </mml:msub>
            <mml:mo>=</mml:mo>
            <mml:msub>
              <mml:mi>v</mml:mi>
              <mml:mo>⊕</mml:mo>
            </mml:msub>
            <mml:mo>−</mml:mo>
            <mml:mi>u</mml:mi>
            <mml:mo>.</mml:mo>
          </mml:mrow>
        </mml:math>
      </disp-formula>
      <p>These expressions assume idealized eastward and westward motion and neglect wind and small changes in latitude during the flight. More generally, the aircraft velocity is obtained by vector addition of the Earth’s rotational velocity and the aircraft’s velocity relative to the Earth’s surface.</p>
      <p>The first-order weak-field clock-shift expression therefore becomes</p>
      <disp-formula id="FD17">
        <label>(15a)</label>
        <mml:math>
          <mml:mrow>
            <mml:mtext>Δ</mml:mtext>
            <mml:mi>τ</mml:mi>
            <mml:mo>=</mml:mo>
            <mml:mi>T</mml:mi>
            <mml:mrow>
              <mml:mo>[</mml:mo>
              <mml:mrow>
                <mml:mfrac>
                  <mml:mrow>
                    <mml:mi>g</mml:mi>
                    <mml:mi>h</mml:mi>
                  </mml:mrow>
                  <mml:mrow>
                    <mml:msup>
                      <mml:mi>c</mml:mi>
                      <mml:mn>2</mml:mn>
                    </mml:msup>
                  </mml:mrow>
                </mml:mfrac>
                <mml:mfrac>
                  <mml:mrow>
                    <mml:msubsup>
                      <mml:mi>v</mml:mi>
                      <mml:mrow>
                        <mml:mtext>air</mml:mtext>
                      </mml:mrow>
                      <mml:mn>2</mml:mn>
                    </mml:msubsup>
                    <mml:mo>−</mml:mo>
                    <mml:msubsup>
                      <mml:mi>v</mml:mi>
                      <mml:mo>⊕</mml:mo>
                      <mml:mn>2</mml:mn>
                    </mml:msubsup>
                  </mml:mrow>
                  <mml:mrow>
                    <mml:mn>2</mml:mn>
                    <mml:msup>
                      <mml:mi>c</mml:mi>
                      <mml:mn>2</mml:mn>
                    </mml:msup>
                  </mml:mrow>
                </mml:mfrac>
              </mml:mrow>
              <mml:mo>]</mml:mo>
            </mml:mrow>
            <mml:mo>,</mml:mo>
          </mml:mrow>
        </mml:math>
      </disp-formula>
      <p>where <inline-formula><mml:math><mml:mrow><mml:msub><mml:mi> v </mml:mi><mml:mrow><mml:mtext> air </mml:mtext></mml:mrow></mml:msub></mml:mrow></mml:math></inline-formula> denotes the aircraft velocity in the Earth-Centered Inertial frame and <inline-formula><mml:math><mml:mrow><mml:msub><mml:mi> v </mml:mi><mml:mo> ⊕ </mml:mo></mml:msub></mml:mrow></mml:math></inline-formula> is the tangential velocity of the ground reference clock at the corresponding latitude. The first term represents the gravitational contribution associated with the higher altitude of the aircraft, while the second represents the kinematic contribution arising from the difference in inertial velocity between the airborne and ground clocks.</p>
      <p>For eastward motion,</p>
      <disp-formula id="FD18">
        <label>(15b)</label>
        <mml:math>
          <mml:mrow>
            <mml:msubsup>
              <mml:mi>v</mml:mi>
              <mml:mrow>
                <mml:mtext>air</mml:mtext>
                <mml:mo>,</mml:mo>
                <mml:mtext>east</mml:mtext>
              </mml:mrow>
              <mml:mn>2</mml:mn>
            </mml:msubsup>
            <mml:mo>−</mml:mo>
            <mml:msubsup>
              <mml:mi>v</mml:mi>
              <mml:mo>⊕</mml:mo>
              <mml:mn>2</mml:mn>
            </mml:msubsup>
            <mml:mo>=</mml:mo>
            <mml:mn>2</mml:mn>
            <mml:msub>
              <mml:mi>v</mml:mi>
              <mml:mo>⊕</mml:mo>
            </mml:msub>
            <mml:mi>u</mml:mi>
            <mml:mo>+</mml:mo>
            <mml:msup>
              <mml:mi>u</mml:mi>
              <mml:mn>2</mml:mn>
            </mml:msup>
            <mml:mo>,</mml:mo>
          </mml:mrow>
        </mml:math>
      </disp-formula>
      <p>which is always positive. Consequently, the kinematic contribution opposes the gravitational contribution and may become dominant, leading to a net loss of time for the airborne clock.</p>
      <p>For westward motion,</p>
      <disp-formula id="FD19">
        <label>(15c)</label>
        <mml:math>
          <mml:mrow>
            <mml:msubsup>
              <mml:mi>v</mml:mi>
              <mml:mrow>
                <mml:mtext>air</mml:mtext>
                <mml:mo>,</mml:mo>
                <mml:mtext>west</mml:mtext>
              </mml:mrow>
              <mml:mn>2</mml:mn>
            </mml:msubsup>
            <mml:mo>−</mml:mo>
            <mml:msubsup>
              <mml:mi>v</mml:mi>
              <mml:mo>⊕</mml:mo>
              <mml:mn>2</mml:mn>
            </mml:msubsup>
            <mml:mo>=</mml:mo>
            <mml:mo>−</mml:mo>
            <mml:mn>2</mml:mn>
            <mml:msub>
              <mml:mi>v</mml:mi>
              <mml:mo>⊕</mml:mo>
            </mml:msub>
            <mml:mi>u</mml:mi>
            <mml:mo>+</mml:mo>
            <mml:msup>
              <mml:mi>u</mml:mi>
              <mml:mn>2</mml:mn>
            </mml:msup>
            <mml:mo>,</mml:mo>
          </mml:mrow>
        </mml:math>
      </disp-formula>
      <p>which is negative for ordinary commercial aircraft because <inline-formula><mml:math display="inline"><mml:mrow><mml:mi> u </mml:mi><mml:mo> &lt; </mml:mo><mml:msub><mml:mi> v </mml:mi><mml:mo> ⊕ </mml:mo></mml:msub></mml:mrow></mml:math></inline-formula> . In this case the kinematic contribution acts in the same direction as the gravitational contribution, producing a net gain of time relative to the ground reference clock.</p>
      <p>Thus, the standard relativistic analysis naturally predicts opposite signs for the two flight directions: eastward flights experience a negative clock shift, whereas westward flights experience a positive clock shift. These experimentally confirmed results provide the benchmark that the structured-vacuum interpretation developed in the following sections seeks to reproduce while offering an alternative physical interpretation of the underlying mechanism.</p>
    </sec>
    <sec id="sec3">
      <title>3. Structured Quantum-Vacuum Vortex Hypothesis</title>
      <p>The preceding section reviewed the standard relativistic explanation of the Hafele-Keating experiment. The present section introduces the physical hypothesis used in this article to reinterpret the same weak-field clock-shift structure. The aim is not to modify the empirical prediction of relativity, but to explore whether the mathematical terms of the relativistic expression may correspond to physical changes in a structured vacuum medium.</p>
      <p>The central assumption is that the quantum vacuum is not an empty geometrical background, but a physically active substrate with superfluid-like properties. In this framework, the vacuum can sustain persistent vortex structures, density gradients, and pressure-like responses without ordinary viscous dissipation. Such a view is compatible in spirit with several approaches in analog gravity, emergent spacetime, and superfluid-vacuum theory, where relativistic phenomena are interpreted as macroscopic manifestations of deeper medium-like degrees of freedom.</p>
      <p>In the model developed here, the Earth is not assumed to create the vacuum vortex merely by mechanical rotation. Rather, the Earth is considered dynamically coupled to a pre-existing large-scale vortex configuration of the structured vacuum. The observed rotation of the Earth is therefore interpreted as an equilibrium motion within this vortex environment. A clock fixed to the Earth’s surface co-rotates with this local vacuum structure, while an aircraft moving eastward or westward experiences a different effective velocity relative to the vortex flow.</p>
      <p>The large-scale density scale used as the reference point is the observed cosmological vacuum-energy density associated with dark energy. This density may be written as</p>
      <disp-formula id="FD20">
        <label>(16)</label>
        <mml:math>
          <mml:mrow>
            <mml:msub>
              <mml:mi>ρ</mml:mi>
              <mml:mi>Λ</mml:mi>
            </mml:msub>
            <mml:mo>=</mml:mo>
            <mml:msub>
              <mml:mi>Ω</mml:mi>
              <mml:mi>Λ</mml:mi>
            </mml:msub>
            <mml:msub>
              <mml:mi>ρ</mml:mi>
              <mml:mi>c</mml:mi>
            </mml:msub>
            <mml:mo>≃</mml:mo>
            <mml:mn>6.3</mml:mn>
            <mml:mo>×</mml:mo>
            <mml:msup>
              <mml:mrow>
                <mml:mn>10</mml:mn>
              </mml:mrow>
              <mml:mrow>
                <mml:mo>−</mml:mo>
                <mml:mn>27</mml:mn>
              </mml:mrow>
            </mml:msup>
            <mml:mtext>
               
            </mml:mtext>
            <mml:mtext>kg</mml:mtext>
            <mml:mo>⋅</mml:mo>
            <mml:msup>
              <mml:mtext>m</mml:mtext>
              <mml:mrow>
                <mml:mo>−</mml:mo>
                <mml:mn>3</mml:mn>
              </mml:mrow>
            </mml:msup>
            <mml:mo>.</mml:mo>
          </mml:mrow>
        </mml:math>
      </disp-formula>
      <p>where <inline-formula><mml:math display="inline"><mml:mrow><mml:msub><mml:mi> ρ </mml:mi><mml:mi> Λ </mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> is the cosmological vacuum-density scale, <inline-formula><mml:math display="inline"><mml:mrow><mml:msub><mml:mi> Ω </mml:mi><mml:mi> Λ </mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> is the dark-energy fraction, and <inline-formula><mml:math display="inline"><mml:mrow><mml:msub><mml:mi> ρ </mml:mi><mml:mi> c </mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> is the critical density of the Universe. Numerically, this gives a value of order 10<sup>−27</sup> kg∙m<sup>−</sup><sup>3</sup>, which provides a physically measurable baseline scale for the proposed vacuum framework.</p>
      <p>However, the density that controls local electromagnetic propagation and atomic transition rates need not be exactly identical to the cosmological mean density. For this reason, the model distinguishes between the cosmological density scale and the local mechanically active vacuum-density scale. This is written as</p>
      <disp-formula id="FD21">
        <label>(17)</label>
        <mml:math>
          <mml:mrow>
            <mml:msub>
              <mml:mi>ρ</mml:mi>
              <mml:mrow>
                <mml:mtext>eff</mml:mtext>
              </mml:mrow>
            </mml:msub>
            <mml:mo>=</mml:mo>
            <mml:mi>κ</mml:mi>
            <mml:msub>
              <mml:mi>ρ</mml:mi>
              <mml:mi>Λ</mml:mi>
            </mml:msub>
            <mml:mo>.</mml:mo>
          </mml:mrow>
        </mml:math>
      </disp-formula>
      <p>where <inline-formula><mml:math display="inline"><mml:mrow><mml:msub><mml:mi> ρ </mml:mi><mml:mrow><mml:mtext> eff </mml:mtext></mml:mrow></mml:msub></mml:mrow></mml:math></inline-formula> is the effective local density scale relevant to clock-rate modulation, and <inline-formula><mml:math display="inline"><mml:mi> κ </mml:mi></mml:math></inline-formula> is a dimensionless factor that relates this local active density to the cosmological vacuum-density scale. In the simplest first-order model, <inline-formula><mml:math display="inline"><mml:mi> κ </mml:mi></mml:math></inline-formula> is expected to be close to unity, but a more complete theory should derive it from the vacuum compressibility, elasticity, and electromagnetic response.</p>
      <p>At the scale of the rotating Earth, the local vortex flow may be represented phenomenologically by an azimuthal velocity field. At radial distance <italic>r</italic> and latitude lambda, the co-rotating component of the Earth-coupled vacuum flow is written as</p>
      <disp-formula id="FD22">
        <label>(18)</label>
        <mml:math>
          <mml:mrow>
            <mml:msub>
              <mml:mstyle mathsize="normal" mathvariant="bold">
                <mml:mi>u</mml:mi>
              </mml:mstyle>
              <mml:mtext>v</mml:mtext>
            </mml:msub>
            <mml:mrow>
              <mml:mo>(</mml:mo>
              <mml:mrow>
                <mml:mi>r</mml:mi>
                <mml:mo>,</mml:mo>
                <mml:mi>λ</mml:mi>
              </mml:mrow>
              <mml:mo>)</mml:mo>
            </mml:mrow>
            <mml:mo>=</mml:mo>
            <mml:msub>
              <mml:mi>Ω</mml:mi>
              <mml:mo>⊕</mml:mo>
            </mml:msub>
            <mml:mi>r</mml:mi>
            <mml:mi>cos</mml:mi>
            <mml:mi>λ</mml:mi>
            <mml:mtext>
               
            </mml:mtext>
            <mml:mover accent="true">
              <mml:mi>ϕ</mml:mi>
              <mml:mo>^</mml:mo>
            </mml:mover>
            <mml:mo>.</mml:mo>
          </mml:mrow>
        </mml:math>
      </disp-formula>
      <p>where <inline-formula><mml:math><mml:mrow><mml:msub><mml:mi> Ω </mml:mi><mml:mo> ⊕ </mml:mo></mml:msub></mml:mrow></mml:math></inline-formula> is the angular velocity of the Earth and <italic>φ</italic>-hat is the azimuthal unit vector. This expression does not claim that the vacuum is a conventional material fluid. Rather, it provides an effective hydrodynamic description of the local rotational component of the structured vacuum state.</p>
      <p>A clock moving through this environment is affected not by its velocity relative to the aircraft cabin or the local air mass, but by its velocity relative to the structured vacuum flow. The relevant relative velocity is therefore</p>
      <disp-formula id="FD23">
        <label>(19)</label>
        <mml:math>
          <mml:mrow>
            <mml:msub>
              <mml:mstyle mathvariant="bold" mathsize="normal">
                <mml:mi>v</mml:mi>
              </mml:mstyle>
              <mml:mrow>
                <mml:mtext>rel</mml:mtext>
              </mml:mrow>
            </mml:msub>
            <mml:mo>=</mml:mo>
            <mml:msub>
              <mml:mstyle mathvariant="bold" mathsize="normal">
                <mml:mi>v</mml:mi>
              </mml:mstyle>
              <mml:mrow>
                <mml:mtext>clock</mml:mtext>
              </mml:mrow>
            </mml:msub>
            <mml:mo>−</mml:mo>
            <mml:msub>
              <mml:mstyle mathvariant="bold" mathsize="normal">
                <mml:mi>u</mml:mi>
              </mml:mstyle>
              <mml:mtext>v</mml:mtext>
            </mml:msub>
            <mml:mo>.</mml:mo>
          </mml:mrow>
        </mml:math>
      </disp-formula>
      <p>Here, the first velocity denotes the clock velocity in the chosen Earth-centered frame, and the second denotes the local vortex-flow velocity. This relative velocity determines the kinetic contribution to the vacuum-density perturbation experienced by the clock.</p>
      <p>The model assumes that small velocity-dependent and altitude-dependent perturbations of the vacuum density can be written as a dimensionless first-order functional of the relative kinetic term and the altitude-dependent gravitational term:</p>
      <disp-formula id="FD24">
        <label>(20)</label>
        <mml:math>
          <mml:mrow>
            <mml:mi>Δ</mml:mi>
            <mml:mi>ρ</mml:mi>
            <mml:mo>=</mml:mo>
            <mml:msub>
              <mml:mi>ρ</mml:mi>
              <mml:mrow>
                <mml:mtext>eff</mml:mtext>
              </mml:mrow>
            </mml:msub>
            <mml:mi>ℱ</mml:mi>
            <mml:mrow>
              <mml:mo>(</mml:mo>
              <mml:mrow>
                <mml:mfrac>
                  <mml:mrow>
                    <mml:msubsup>
                      <mml:mi>v</mml:mi>
                      <mml:mrow>
                        <mml:mtext>rel</mml:mtext>
                      </mml:mrow>
                      <mml:mn>2</mml:mn>
                    </mml:msubsup>
                  </mml:mrow>
                  <mml:mrow>
                    <mml:msup>
                      <mml:mi>c</mml:mi>
                      <mml:mn>2</mml:mn>
                    </mml:msup>
                  </mml:mrow>
                </mml:mfrac>
                <mml:mo>,</mml:mo>
                <mml:mfrac>
                  <mml:mrow>
                    <mml:mi>g</mml:mi>
                    <mml:mi>h</mml:mi>
                  </mml:mrow>
                  <mml:mrow>
                    <mml:msup>
                      <mml:mi>c</mml:mi>
                      <mml:mn>2</mml:mn>
                    </mml:msup>
                  </mml:mrow>
                </mml:mfrac>
              </mml:mrow>
              <mml:mo>)</mml:mo>
            </mml:mrow>
            <mml:mo>.</mml:mo>
          </mml:mrow>
        </mml:math>
      </disp-formula>
      <p>This expression summarizes the physical content of the hypothesis. The first argument describes the kinetic disturbance associated with motion relative to the vortex flow. The second describes the radial disturbance associated with altitude in the weak gravitational field of the Earth. Both quantities are extremely small under ordinary flight conditions, which justifies the use of first-order approximations.</p>
      <p>The distinction between the cosmological density scale and the effective local density scale is important for avoiding an overstatement of the model. The observed dark-energy density provides a cosmological density scale, but the local density that participates in electromagnetic propagation may represent only the mechanically active component of the structured vacuum. The Hafele-Keating comparison therefore serves not as a direct measurement of dark energy, but as a consistency test for whether the required effective density scale is of the same order as the observed vacuum-density scale.</p>
      <p>In the next section, this structured-vacuum hypothesis is connected to the propagation speed of electromagnetic interactions. The clock-rate shift will then be derived from the dependence of wave speed on effective vacuum density, leading to a direct relation between vacuum-density perturbation and atomic time dilation.</p>
    </sec>
    <sec id="sec4">
      <title>4. Vacuum Density, Electromagnetic Propagation Speed, and Clock Rate</title>
      <p>The previous section introduced the structured quantum-vacuum vortex hypothesis and distinguished the cosmological vacuum-density scale from the mechanically active local density scale. The present section develops the physical link between a small vacuum-density perturbation and the rate of an atomic clock. This link is essential because the Hafele-Keating experiment measured time differences, whereas the proposed mechanism begins with changes in the local vacuum-density state.</p>
      <p>The guiding assumption is that electromagnetic propagation through the structured vacuum behaves, to first order, like wave propagation through an elastic medium. In ordinary continuum physics, wave speed depends on the ratio between an elastic modulus and density. By analogy, the effective speed of electromagnetic propagation in the locally perturbed vacuum may be written as</p>
      <disp-formula id="FD25">
        <label>(21)</label>
        <mml:math>
          <mml:mrow>
            <mml:msub>
              <mml:mi>c</mml:mi>
              <mml:mrow>
                <mml:mtext>eff</mml:mtext>
              </mml:mrow>
            </mml:msub>
            <mml:mo>=</mml:mo>
            <mml:msqrt>
              <mml:mrow>
                <mml:mfrac>
                  <mml:mrow>
                    <mml:msub>
                      <mml:mi>K</mml:mi>
                      <mml:mrow>
                        <mml:mtext>vac</mml:mtext>
                      </mml:mrow>
                    </mml:msub>
                  </mml:mrow>
                  <mml:mrow>
                    <mml:msub>
                      <mml:mi>ρ</mml:mi>
                      <mml:mrow>
                        <mml:mtext>eff</mml:mtext>
                      </mml:mrow>
                    </mml:msub>
                    <mml:mo>+</mml:mo>
                    <mml:mi>Δ</mml:mi>
                    <mml:mi>ρ</mml:mi>
                  </mml:mrow>
                </mml:mfrac>
              </mml:mrow>
            </mml:msqrt>
            <mml:mo>.</mml:mo>
          </mml:mrow>
        </mml:math>
      </disp-formula>
      <p>where <inline-formula><mml:math><mml:mrow><mml:msub><mml:mi> c </mml:mi><mml:mrow><mml:mtext> eff </mml:mtext></mml:mrow></mml:msub></mml:mrow></mml:math></inline-formula> is the effective local propagation speed of electromagnetic interactions, <inline-formula><mml:math><mml:mrow><mml:msub><mml:mi> K </mml:mi><mml:mrow><mml:mtext> vac </mml:mtext></mml:mrow></mml:msub></mml:mrow></mml:math></inline-formula> is the effective elastic modulus of the structured vacuum, <inline-formula><mml:math><mml:mrow><mml:msub><mml:mi> ρ </mml:mi><mml:mrow><mml:mtext> eff </mml:mtext></mml:mrow></mml:msub></mml:mrow></mml:math></inline-formula> is the mechanically active reference density, and <inline-formula><mml:math><mml:mrow><mml:mi> Δ </mml:mi><mml:mi> ρ </mml:mi></mml:mrow></mml:math></inline-formula> is the local perturbation of that density.</p>
      <p>The unperturbed reference speed is obtained when the perturbation vanishes:</p>
      <disp-formula id="FD26">
        <label>(22)</label>
        <mml:math>
          <mml:mrow>
            <mml:mi>c</mml:mi>
            <mml:mo>=</mml:mo>
            <mml:msqrt>
              <mml:mrow>
                <mml:mfrac>
                  <mml:mrow>
                    <mml:msub>
                      <mml:mi>K</mml:mi>
                      <mml:mrow>
                        <mml:mtext>vac</mml:mtext>
                      </mml:mrow>
                    </mml:msub>
                  </mml:mrow>
                  <mml:mrow>
                    <mml:msub>
                      <mml:mi>ρ</mml:mi>
                      <mml:mrow>
                        <mml:mtext>eff</mml:mtext>
                      </mml:mrow>
                    </mml:msub>
                  </mml:mrow>
                </mml:mfrac>
              </mml:mrow>
            </mml:msqrt>
            <mml:mo>.</mml:mo>
          </mml:mrow>
        </mml:math>
      </disp-formula>
      <p>Dividing Equation (21) by Equation (22) gives</p>
      <disp-formula id="FD27">
        <label>(23)</label>
        <mml:math>
          <mml:mrow>
            <mml:mfrac>
              <mml:mrow>
                <mml:msub>
                  <mml:mi>c</mml:mi>
                  <mml:mrow>
                    <mml:mtext>eff</mml:mtext>
                  </mml:mrow>
                </mml:msub>
              </mml:mrow>
              <mml:mi>c</mml:mi>
            </mml:mfrac>
            <mml:mo>=</mml:mo>
            <mml:msup>
              <mml:mrow>
                <mml:mrow>
                  <mml:mo>(</mml:mo>
                  <mml:mrow>
                    <mml:mn>1</mml:mn>
                    <mml:mo>+</mml:mo>
                    <mml:mfrac>
                      <mml:mrow>
                        <mml:mi>Δ</mml:mi>
                        <mml:mi>ρ</mml:mi>
                      </mml:mrow>
                      <mml:mrow>
                        <mml:msub>
                          <mml:mi>ρ</mml:mi>
                          <mml:mrow>
                            <mml:mtext>eff</mml:mtext>
                          </mml:mrow>
                        </mml:msub>
                      </mml:mrow>
                    </mml:mfrac>
                  </mml:mrow>
                  <mml:mo>)</mml:mo>
                </mml:mrow>
              </mml:mrow>
              <mml:mrow>
                <mml:mo>−</mml:mo>
                <mml:mrow>
                  <mml:mn>1</mml:mn>
                  <mml:mo>/</mml:mo>
                  <mml:mn>2</mml:mn>
                </mml:mrow>
              </mml:mrow>
            </mml:msup>
            <mml:mo>.</mml:mo>
          </mml:mrow>
        </mml:math>
      </disp-formula>
      <p>Equation (23) is the central wave-speed relation used in this article. It states that a positive density perturbation reduces the effective propagation speed, whereas a negative density perturbation increases it. In physical terms, a slightly denser local vacuum state acts as a slightly more resistant propagation environment for electromagnetic processes, while a rarefied state allows a slightly faster effective propagation.</p>
      <p>The perturbations considered in the Hafele-Keating regime are extremely small. Therefore, the first-order approximation is valid whenever</p>
      <disp-formula id="FD28">
        <label>(24)</label>
        <mml:math>
          <mml:mrow>
            <mml:mrow>
              <mml:mo>|</mml:mo>
              <mml:mrow>
                <mml:mfrac>
                  <mml:mrow>
                    <mml:mi>Δ</mml:mi>
                    <mml:mi>ρ</mml:mi>
                  </mml:mrow>
                  <mml:mrow>
                    <mml:msub>
                      <mml:mi>ρ</mml:mi>
                      <mml:mrow>
                        <mml:mtext>eff</mml:mtext>
                      </mml:mrow>
                    </mml:msub>
                  </mml:mrow>
                </mml:mfrac>
              </mml:mrow>
              <mml:mo>|</mml:mo>
            </mml:mrow>
            <mml:mo>≪</mml:mo>
            <mml:mn>1.</mml:mn>
          </mml:mrow>
        </mml:math>
      </disp-formula>
      <p>Using the binomial approximation</p>
      <disp-formula id="FD29">
        <label>(25)</label>
        <mml:math>
          <mml:mrow>
            <mml:msup>
              <mml:mrow>
                <mml:mrow>
                  <mml:mo>(</mml:mo>
                  <mml:mrow>
                    <mml:mn>1</mml:mn>
                    <mml:mo>+</mml:mo>
                    <mml:mi>x</mml:mi>
                  </mml:mrow>
                  <mml:mo>)</mml:mo>
                </mml:mrow>
              </mml:mrow>
              <mml:mi>a</mml:mi>
            </mml:msup>
            <mml:mo>≃</mml:mo>
            <mml:mn>1</mml:mn>
            <mml:mo>+</mml:mo>
            <mml:mi>a</mml:mi>
            <mml:mi>x</mml:mi>
            <mml:mo>.</mml:mo>
          </mml:mrow>
        </mml:math>
      </disp-formula>
      <p>with <italic>x</italic> equal to Delta rho divided by <inline-formula><mml:math><mml:mrow><mml:msub><mml:mi> ρ </mml:mi><mml:mrow><mml:mtext> eff </mml:mtext></mml:mrow></mml:msub></mml:mrow></mml:math></inline-formula> and a equal to −1/2, Equation (23) becomes</p>
      <disp-formula id="FD30">
        <label>(26)</label>
        <mml:math>
          <mml:mrow>
            <mml:msub>
              <mml:mi>c</mml:mi>
              <mml:mrow>
                <mml:mtext>eff</mml:mtext>
              </mml:mrow>
            </mml:msub>
            <mml:mo>≃</mml:mo>
            <mml:mi>c</mml:mi>
            <mml:mrow>
              <mml:mo>(</mml:mo>
              <mml:mrow>
                <mml:mn>1</mml:mn>
                <mml:mo>−</mml:mo>
                <mml:mfrac>
                  <mml:mn>1</mml:mn>
                  <mml:mn>2</mml:mn>
                </mml:mfrac>
                <mml:mfrac>
                  <mml:mrow>
                    <mml:mi>Δ</mml:mi>
                    <mml:mi>ρ</mml:mi>
                  </mml:mrow>
                  <mml:mrow>
                    <mml:msub>
                      <mml:mi>ρ</mml:mi>
                      <mml:mrow>
                        <mml:mtext>eff</mml:mtext>
                      </mml:mrow>
                    </mml:msub>
                  </mml:mrow>
                </mml:mfrac>
              </mml:mrow>
              <mml:mo>)</mml:mo>
            </mml:mrow>
            <mml:mo>.</mml:mo>
          </mml:mrow>
        </mml:math>
      </disp-formula>
      <p>The corresponding fractional change in the effective propagation speed is therefore</p>
      <disp-formula id="FD31">
        <label>(27)</label>
        <mml:math>
          <mml:mrow>
            <mml:mfrac>
              <mml:mrow>
                <mml:mi>Δ</mml:mi>
                <mml:mi>c</mml:mi>
              </mml:mrow>
              <mml:mi>c</mml:mi>
            </mml:mfrac>
            <mml:mo>=</mml:mo>
            <mml:mfrac>
              <mml:mrow>
                <mml:msub>
                  <mml:mi>c</mml:mi>
                  <mml:mrow>
                    <mml:mtext>eff</mml:mtext>
                  </mml:mrow>
                </mml:msub>
                <mml:mo>−</mml:mo>
                <mml:mi>c</mml:mi>
              </mml:mrow>
              <mml:mi>c</mml:mi>
            </mml:mfrac>
            <mml:mo>≃</mml:mo>
            <mml:mo>−</mml:mo>
            <mml:mfrac>
              <mml:mn>1</mml:mn>
              <mml:mn>2</mml:mn>
            </mml:mfrac>
            <mml:mfrac>
              <mml:mrow>
                <mml:mi>Δ</mml:mi>
                <mml:mi>ρ</mml:mi>
              </mml:mrow>
              <mml:mrow>
                <mml:msub>
                  <mml:mi>ρ</mml:mi>
                  <mml:mrow>
                    <mml:mtext>eff</mml:mtext>
                  </mml:mrow>
                </mml:msub>
              </mml:mrow>
            </mml:mfrac>
            <mml:mo>.</mml:mo>
          </mml:mrow>
        </mml:math>
      </disp-formula>
      <p>This result provides the physical origin of the factor 1/2 in the clock-shift relation. It does not arise from an arbitrary fitting coefficient, but from the inverse square-root dependence of wave speed on density.</p>
      <p>Atomic clocks are based on the stability of atomic transition frequencies, with the SI second defined by the unperturbed hyperfine transition frequency of caesium-133 at 9,192,631,770 Hz [<xref ref-type="bibr" rid="B11">11</xref>]. The experimental realization of this definition has historically relied on molecular-beam and caesium-resonator methods, and later on caesium fountain primary frequency standards [<xref ref-type="bibr" rid="B12">12</xref>]-[<xref ref-type="bibr" rid="B14">14</xref>]. Modern optical clocks extend this metrological framework to even higher precision [<xref ref-type="bibr" rid="B15">15</xref>].</p>
      <p>In the present framework, these transition frequencies are assumed to scale, to first order, with the effective electromagnetic propagation speed of the local vacuum state. This gives</p>
      <disp-formula id="FD32">
        <label>(28)</label>
        <mml:math>
          <mml:mrow>
            <mml:mfrac>
              <mml:mrow>
                <mml:msub>
                  <mml:mi>ν</mml:mi>
                  <mml:mrow>
                    <mml:mtext>eff</mml:mtext>
                  </mml:mrow>
                </mml:msub>
              </mml:mrow>
              <mml:mrow>
                <mml:msub>
                  <mml:mi>ν</mml:mi>
                  <mml:mn>0</mml:mn>
                </mml:msub>
              </mml:mrow>
            </mml:mfrac>
            <mml:mo>≃</mml:mo>
            <mml:mfrac>
              <mml:mrow>
                <mml:msub>
                  <mml:mi>c</mml:mi>
                  <mml:mrow>
                    <mml:mtext>eff</mml:mtext>
                  </mml:mrow>
                </mml:msub>
              </mml:mrow>
              <mml:mi>c</mml:mi>
            </mml:mfrac>
            <mml:mo>.</mml:mo>
          </mml:mrow>
        </mml:math>
      </disp-formula>
      <p>where <inline-formula><mml:math><mml:mrow><mml:msub><mml:mi> ν </mml:mi><mml:mrow><mml:mtext> eff </mml:mtext></mml:mrow></mml:msub></mml:mrow></mml:math></inline-formula> is the effective atomic transition frequency in the perturbed vacuum state and <inline-formula><mml:math><mml:mrow><mml:msub><mml:mi> ν </mml:mi><mml:mn> 0 </mml:mn></mml:msub></mml:mrow></mml:math></inline-formula> is the corresponding unperturbed reference frequency. The proportionality between atomic transition frequency and effective electromagnetic propagation speed is a phenomenological assumption of the present model. A complete microscopic theory would need to derive this coupling from the dependence of atomic energy levels on the electromagnetic properties of the structured vacuum. Since the clock accumulates time by counting periodic atomic transitions, the fractional clock-rate variation follows the fractional change in transition frequency:</p>
      <disp-formula id="FD33">
        <label>(29)</label>
        <mml:math>
          <mml:mrow>
            <mml:mfrac>
              <mml:mrow>
                <mml:mi>Δ</mml:mi>
                <mml:mi>τ</mml:mi>
              </mml:mrow>
              <mml:mi>T</mml:mi>
            </mml:mfrac>
            <mml:mo>≃</mml:mo>
            <mml:mfrac>
              <mml:mrow>
                <mml:mi>Δ</mml:mi>
                <mml:mi>ν</mml:mi>
              </mml:mrow>
              <mml:mrow>
                <mml:msub>
                  <mml:mi>ν</mml:mi>
                  <mml:mn>0</mml:mn>
                </mml:msub>
              </mml:mrow>
            </mml:mfrac>
            <mml:mo>≃</mml:mo>
            <mml:mfrac>
              <mml:mrow>
                <mml:mi>Δ</mml:mi>
                <mml:mi>c</mml:mi>
              </mml:mrow>
              <mml:mi>c</mml:mi>
            </mml:mfrac>
            <mml:mo>.</mml:mo>
          </mml:mrow>
        </mml:math>
      </disp-formula>
      <p>Combining Equations (27) and (29) yields</p>
      <disp-formula id="FD34">
        <label>(30)</label>
        <mml:math>
          <mml:mrow>
            <mml:mfrac>
              <mml:mrow>
                <mml:mi>Δ</mml:mi>
                <mml:mi>τ</mml:mi>
              </mml:mrow>
              <mml:mi>T</mml:mi>
            </mml:mfrac>
            <mml:mo>≃</mml:mo>
            <mml:mo>−</mml:mo>
            <mml:mfrac>
              <mml:mn>1</mml:mn>
              <mml:mn>2</mml:mn>
            </mml:mfrac>
            <mml:mfrac>
              <mml:mrow>
                <mml:mi>Δ</mml:mi>
                <mml:mi>ρ</mml:mi>
              </mml:mrow>
              <mml:mrow>
                <mml:msub>
                  <mml:mi>ρ</mml:mi>
                  <mml:mrow>
                    <mml:mtext>eff</mml:mtext>
                  </mml:mrow>
                </mml:msub>
              </mml:mrow>
            </mml:mfrac>
            <mml:mo>.</mml:mo>
          </mml:mrow>
        </mml:math>
      </disp-formula>
      <p>Equation (30) is the required bridge between vacuum-density perturbation and measured clock-time variation. A positive density perturbation produces a negative time shift, meaning that the clock ticks more slowly relative to the reference clock. A negative density perturbation produces a positive time shift, meaning that the clock ticks faster relative to the reference clock.</p>
      <p>The same relation may also be expressed in terms of an effective refractive index. Since refractive index is defined as the ratio between the reference propagation speed and the effective propagation speed, one obtains</p>
      <disp-formula id="FD35">
        <label>(31)</label>
        <mml:math>
          <mml:mrow>
            <mml:msub>
              <mml:mi>n</mml:mi>
              <mml:mrow>
                <mml:mtext>eff</mml:mtext>
              </mml:mrow>
            </mml:msub>
            <mml:mo>=</mml:mo>
            <mml:mfrac>
              <mml:mi>c</mml:mi>
              <mml:mrow>
                <mml:msub>
                  <mml:mi>c</mml:mi>
                  <mml:mrow>
                    <mml:mtext>eff</mml:mtext>
                  </mml:mrow>
                </mml:msub>
              </mml:mrow>
            </mml:mfrac>
            <mml:mo>=</mml:mo>
            <mml:msup>
              <mml:mrow>
                <mml:mrow>
                  <mml:mo>(</mml:mo>
                  <mml:mrow>
                    <mml:mn>1</mml:mn>
                    <mml:mo>+</mml:mo>
                    <mml:mfrac>
                      <mml:mrow>
                        <mml:mi>Δ</mml:mi>
                        <mml:mi>ρ</mml:mi>
                      </mml:mrow>
                      <mml:mrow>
                        <mml:msub>
                          <mml:mi>ρ</mml:mi>
                          <mml:mrow>
                            <mml:mtext>eff</mml:mtext>
                          </mml:mrow>
                        </mml:msub>
                      </mml:mrow>
                    </mml:mfrac>
                  </mml:mrow>
                  <mml:mo>)</mml:mo>
                </mml:mrow>
              </mml:mrow>
              <mml:mrow>
                <mml:mrow>
                  <mml:mn>1</mml:mn>
                  <mml:mo>/</mml:mo>
                  <mml:mn>2</mml:mn>
                </mml:mrow>
              </mml:mrow>
            </mml:msup>
            <mml:mo>.</mml:mo>
          </mml:mrow>
        </mml:math>
      </disp-formula>
      <p>and, for small perturbations,</p>
      <disp-formula id="FD36">
        <label>(32)</label>
        <mml:math>
          <mml:mrow>
            <mml:msub>
              <mml:mi>n</mml:mi>
              <mml:mrow>
                <mml:mtext>eff</mml:mtext>
              </mml:mrow>
            </mml:msub>
            <mml:mo>≃</mml:mo>
            <mml:mn>1</mml:mn>
            <mml:mo>+</mml:mo>
            <mml:mfrac>
              <mml:mn>1</mml:mn>
              <mml:mn>2</mml:mn>
            </mml:mfrac>
            <mml:mfrac>
              <mml:mrow>
                <mml:mi>Δ</mml:mi>
                <mml:mi>ρ</mml:mi>
              </mml:mrow>
              <mml:mrow>
                <mml:msub>
                  <mml:mi>ρ</mml:mi>
                  <mml:mrow>
                    <mml:mtext>eff</mml:mtext>
                  </mml:mrow>
                </mml:msub>
              </mml:mrow>
            </mml:mfrac>
            <mml:mo>.</mml:mo>
          </mml:mrow>
        </mml:math>
      </disp-formula>
      <p>Thus, the density-based and refractive-index-based descriptions are mathematically equivalent at first order. An increase in effective vacuum density corresponds to a slight increase in the effective refractive index and a slight decrease in the effective speed of electromagnetic propagation.</p>
      <p>This formulation has two important consequences for the present article. First, the clock shift depends on the fractional perturbation Delta rho divided by rho<sub>eff</sub>, not only on the absolute perturbation Delta rho. Second, the effective density scale rho<sub>eff</sub> must be interpreted carefully. It is not necessarily identical to the cosmological vacuum-energy density, but it is expected to be of the same order if the cosmological vacuum state provides the underlying density scale of the structured medium.</p>
    </sec>
    <sec id="sec5">
      <title>5. Independent Derivation of Vacuum-Density Perturbation from Vortex Dynamics</title>
      <p>The previous section established the relationship between a small vacuum-density perturbation and the fractional change in atomic clock rate. The present section derives the perturbation itself from measurable flight parameters. This step is essential because the model must not infer the density perturbation from the observed Hafele-Keating clock shifts. Instead, the perturbation is obtained independently from the kinematics of motion in an Earth-Centered Inertial (ECI) frame together with the altitude-dependent radial rarefaction of the structured vacuum.</p>
      <p>In the structured-vacuum interpretation, a clock fixed to the Earth’s surface is not absolutely at rest. Owing to Earth’s rotation, it possesses a tangential velocity relative to an Earth-Centered Inertial frame. At latitude <inline-formula><mml:math><mml:mi> λ </mml:mi></mml:math></inline-formula> , this velocity is</p>
      <disp-formula id="FD37">
        <label>(33)</label>
        <mml:math>
          <mml:mrow>
            <mml:msub>
              <mml:mi>v</mml:mi>
              <mml:mo>⊕</mml:mo>
            </mml:msub>
            <mml:mrow>
              <mml:mo>(</mml:mo>
              <mml:mi>λ</mml:mi>
              <mml:mo>)</mml:mo>
            </mml:mrow>
            <mml:mo>=</mml:mo>
            <mml:msub>
              <mml:mi>Ω</mml:mi>
              <mml:mo>⊕</mml:mo>
            </mml:msub>
            <mml:msub>
              <mml:mi>R</mml:mi>
              <mml:mo>⊕</mml:mo>
            </mml:msub>
            <mml:mi>cos</mml:mi>
            <mml:mi>λ</mml:mi>
            <mml:mo>.</mml:mo>
          </mml:mrow>
        </mml:math>
      </disp-formula>
      <p>where <inline-formula><mml:math display="inline"><mml:mrow><mml:msub><mml:mi> Ω </mml:mi><mml:mo> ⊕ </mml:mo></mml:msub></mml:mrow></mml:math></inline-formula> is the angular velocity of the Earth, <inline-formula><mml:math display="inline"><mml:mrow><mml:msub><mml:mi> R </mml:mi><mml:mo> ⊕ </mml:mo></mml:msub></mml:mrow></mml:math></inline-formula> is the Earth’s radius, and <inline-formula><mml:math><mml:mi> λ </mml:mi></mml:math></inline-formula> is the geographic latitude. This velocity represents the inertial-frame velocity of the ground reference clock at the same latitude as the projected flight path..</p>
      <p>An aircraft moving with ground-relative speed <inline-formula><mml:math><mml:mi> u </mml:mi></mml:math></inline-formula> has an inertial-frame velocity obtained by combining its velocity relative to the Earth’s surface with the rotational velocity of the Earth. Assuming idealized east-west flight and neglecting wind, the aircraft velocity becomes</p>
      <disp-formula id="FD38">
        <label>(34)</label>
        <mml:math>
          <mml:mrow>
            <mml:msub>
              <mml:mi>v</mml:mi>
              <mml:mrow>
                <mml:mtext>air</mml:mtext>
              </mml:mrow>
            </mml:msub>
            <mml:mo>=</mml:mo>
            <mml:msub>
              <mml:mi>v</mml:mi>
              <mml:mo>⊕</mml:mo>
            </mml:msub>
            <mml:mo>±</mml:mo>
            <mml:mi>u</mml:mi>
            <mml:mo>.</mml:mo>
          </mml:mrow>
        </mml:math>
      </disp-formula>
      <p>where the positive plus sign applies to eastward motion and the negative sign to westward motion. Eastward flight therefore increases the aircraft’s inertial velocity relative to the ground reference clock, whereas westward flight decreases it.</p>
      <p>The kinetic contribution to the local vacuum-density perturbation is assumed to be proportional to the difference in squared inertial velocities between the airborne and ground clocks. To first order,</p>
      <disp-formula id="FD39">
        <label>(35)</label>
        <mml:math>
          <mml:mrow>
            <mml:msub>
              <mml:mrow>
                <mml:mrow>
                  <mml:mo>(</mml:mo>
                  <mml:mrow>
                    <mml:mfrac>
                      <mml:mrow>
                        <mml:mi>Δ</mml:mi>
                        <mml:mi>ρ</mml:mi>
                      </mml:mrow>
                      <mml:mrow>
                        <mml:msub>
                          <mml:mi>ρ</mml:mi>
                          <mml:mrow>
                            <mml:mtext>eff</mml:mtext>
                          </mml:mrow>
                        </mml:msub>
                      </mml:mrow>
                    </mml:mfrac>
                  </mml:mrow>
                  <mml:mo>)</mml:mo>
                </mml:mrow>
              </mml:mrow>
              <mml:mrow>
                <mml:mtext>kin</mml:mtext>
              </mml:mrow>
            </mml:msub>
            <mml:mo>=</mml:mo>
            <mml:mfrac>
              <mml:mrow>
                <mml:msubsup>
                  <mml:mi>v</mml:mi>
                  <mml:mrow>
                    <mml:mtext>air</mml:mtext>
                  </mml:mrow>
                  <mml:mn>2</mml:mn>
                </mml:msubsup>
                <mml:mo>−</mml:mo>
                <mml:msubsup>
                  <mml:mi>v</mml:mi>
                  <mml:mo>⊕</mml:mo>
                  <mml:mn>2</mml:mn>
                </mml:msubsup>
              </mml:mrow>
              <mml:mrow>
                <mml:msup>
                  <mml:mi>c</mml:mi>
                  <mml:mn>2</mml:mn>
                </mml:msup>
              </mml:mrow>
            </mml:mfrac>
            <mml:mo>.</mml:mo>
          </mml:mrow>
        </mml:math>
      </disp-formula>
      <p>This phenomenological relation represents the additional kinetic loading of the effective vacuum state experienced by the airborne clock relative to the ground reference clock. A positive value corresponds to an increase in the effective vacuum density and, according to the wave-speed relation developed in Section 4, leads to a slower clock rate. Conversely, a negative value corresponds to an effective rarefaction of the vacuum state and therefore to a faster clock rate.</p>
      <p>The second contribution arises from altitude. The aircraft travels at height <italic>h</italic> above the ground reference clock. Within the present framework, increasing altitude corresponds to a slight decrease in the effective vacuum-density state. In the weak-field and low-altitude limit, the altitude-dependent contribution is written as</p>
      <disp-formula id="FD40">
        <label>(36)</label>
        <mml:math>
          <mml:mrow>
            <mml:msub>
              <mml:mrow>
                <mml:mrow>
                  <mml:mo>(</mml:mo>
                  <mml:mrow>
                    <mml:mfrac>
                      <mml:mrow>
                        <mml:mi>Δ</mml:mi>
                        <mml:mi>ρ</mml:mi>
                      </mml:mrow>
                      <mml:mrow>
                        <mml:msub>
                          <mml:mi>ρ</mml:mi>
                          <mml:mrow>
                            <mml:mtext>eff</mml:mtext>
                          </mml:mrow>
                        </mml:msub>
                      </mml:mrow>
                    </mml:mfrac>
                  </mml:mrow>
                  <mml:mo>)</mml:mo>
                </mml:mrow>
              </mml:mrow>
              <mml:mrow>
                <mml:mtext>alt</mml:mtext>
              </mml:mrow>
            </mml:msub>
            <mml:mo>=</mml:mo>
            <mml:mo>−</mml:mo>
            <mml:mfrac>
              <mml:mrow>
                <mml:mn>2</mml:mn>
                <mml:mi>g</mml:mi>
                <mml:mi>h</mml:mi>
              </mml:mrow>
              <mml:mrow>
                <mml:msup>
                  <mml:mi>c</mml:mi>
                  <mml:mn>2</mml:mn>
                </mml:msup>
              </mml:mrow>
            </mml:mfrac>
            <mml:mo>.</mml:mo>
          </mml:mrow>
        </mml:math>
      </disp-formula>
      <p>where <inline-formula><mml:math display="inline"><mml:mi> g </mml:mi></mml:math></inline-formula> is the gravitational acceleration near the Earth’s surface and <inline-formula><mml:math display="inline"><mml:mi> c </mml:mi></mml:math></inline-formula> is the speed of light. The negative sign denotes radial rarefaction of the effective vacuum state with increasing altitude. This term is the density-based counterpart of the positive gravitational contribution appearing in the standard weak-field relativistic approximation.</p>
      <p>The total fractional vacuum-density perturbation experienced by the airborne clock is therefore:</p>
      <disp-formula id="FD41">
        <label>(37)</label>
        <mml:math>
          <mml:mrow>
            <mml:mfrac>
              <mml:mrow>
                <mml:mi>Δ</mml:mi>
                <mml:mi>ρ</mml:mi>
              </mml:mrow>
              <mml:mrow>
                <mml:msub>
                  <mml:mi>ρ</mml:mi>
                  <mml:mrow>
                    <mml:mtext>eff</mml:mtext>
                  </mml:mrow>
                </mml:msub>
              </mml:mrow>
            </mml:mfrac>
            <mml:mo>=</mml:mo>
            <mml:mfrac>
              <mml:mrow>
                <mml:msubsup>
                  <mml:mi>v</mml:mi>
                  <mml:mrow>
                    <mml:mtext>air</mml:mtext>
                  </mml:mrow>
                  <mml:mn>2</mml:mn>
                </mml:msubsup>
                <mml:mo>−</mml:mo>
                <mml:msubsup>
                  <mml:mi>v</mml:mi>
                  <mml:mo>⊕</mml:mo>
                  <mml:mn>2</mml:mn>
                </mml:msubsup>
              </mml:mrow>
              <mml:mrow>
                <mml:msup>
                  <mml:mi>c</mml:mi>
                  <mml:mn>2</mml:mn>
                </mml:msup>
              </mml:mrow>
            </mml:mfrac>
            <mml:mo>−</mml:mo>
            <mml:mfrac>
              <mml:mrow>
                <mml:mn>2</mml:mn>
                <mml:mi>g</mml:mi>
                <mml:mi>h</mml:mi>
              </mml:mrow>
              <mml:mrow>
                <mml:msup>
                  <mml:mi>c</mml:mi>
                  <mml:mn>2</mml:mn>
                </mml:msup>
              </mml:mrow>
            </mml:mfrac>
            <mml:mo>.</mml:mo>
          </mml:mrow>
        </mml:math>
      </disp-formula>
      <p>Equation (37) is the central result of the present phenomenological density derivation. It depends only on measurable flight parameters and does not use the observed Hafele-Keating clock shifts as input. Consequently, the model avoids the circular procedure of inferring the density perturbation from the measured time shifts and subsequently using that inferred perturbation to reproduce the same observations.</p>
      <p>The directional dependence follows immediately. For eastward motion, the effective velocity is</p>
      <disp-formula id="FD42">
        <label>(38)</label>
        <mml:math>
          <mml:mrow>
            <mml:msub>
              <mml:mi>v</mml:mi>
              <mml:mrow>
                <mml:mtext>air</mml:mtext>
                <mml:mo>,</mml:mo>
                <mml:mtext>east</mml:mtext>
              </mml:mrow>
            </mml:msub>
            <mml:mo>=</mml:mo>
            <mml:msub>
              <mml:mi>v</mml:mi>
              <mml:mo>⊕</mml:mo>
            </mml:msub>
            <mml:mo>+</mml:mo>
            <mml:mi>u</mml:mi>
            <mml:mo>.</mml:mo>
          </mml:mrow>
        </mml:math>
      </disp-formula>
      <p>which gives</p>
      <disp-formula id="FD43">
        <label>(39)</label>
        <mml:math>
          <mml:mrow>
            <mml:msubsup>
              <mml:mi>v</mml:mi>
              <mml:mrow>
                <mml:mtext>air</mml:mtext>
                <mml:mo>,</mml:mo>
                <mml:mtext>east</mml:mtext>
              </mml:mrow>
              <mml:mn>2</mml:mn>
            </mml:msubsup>
            <mml:mo>−</mml:mo>
            <mml:msubsup>
              <mml:mi>v</mml:mi>
              <mml:mo>⊕</mml:mo>
              <mml:mn>2</mml:mn>
            </mml:msubsup>
            <mml:mo>=</mml:mo>
            <mml:mn>2</mml:mn>
            <mml:msub>
              <mml:mi>v</mml:mi>
              <mml:mo>⊕</mml:mo>
            </mml:msub>
            <mml:mi>u</mml:mi>
            <mml:mo>+</mml:mo>
            <mml:msup>
              <mml:mi>u</mml:mi>
              <mml:mn>2</mml:mn>
            </mml:msup>
            <mml:mo>.</mml:mo>
          </mml:mrow>
        </mml:math>
      </disp-formula>
      <p>Since both terms on the right-hand side are positive, the kinetic contribution is always positive. Eastward motion therefore increases the effective vacuum-density loading. If this increase exceeds the altitude-induced rarefaction, the total density perturbation becomes positive and the airborne clock runs more slowly than the ground reference clock.</p>
      <p>For westward motion, </p>
      <disp-formula id="FD44">
        <label>(40)</label>
        <mml:math>
          <mml:mrow>
            <mml:msub>
              <mml:mi>v</mml:mi>
              <mml:mrow>
                <mml:mtext>air</mml:mtext>
                <mml:mo>,</mml:mo>
                <mml:mtext>west</mml:mtext>
              </mml:mrow>
            </mml:msub>
            <mml:mo>=</mml:mo>
            <mml:msub>
              <mml:mi>v</mml:mi>
              <mml:mo>⊕</mml:mo>
            </mml:msub>
            <mml:mo>−</mml:mo>
            <mml:mi>u</mml:mi>
            <mml:mo>.</mml:mo>
          </mml:mrow>
        </mml:math>
      </disp-formula>
      <p>and therefore</p>
      <disp-formula id="FD45">
        <label>(41)</label>
        <mml:math>
          <mml:mrow>
            <mml:msubsup>
              <mml:mi>v</mml:mi>
              <mml:mrow>
                <mml:mtext>air</mml:mtext>
                <mml:mo>,</mml:mo>
                <mml:mtext>west</mml:mtext>
              </mml:mrow>
              <mml:mn>2</mml:mn>
            </mml:msubsup>
            <mml:mo>−</mml:mo>
            <mml:msubsup>
              <mml:mi>v</mml:mi>
              <mml:mo>⊕</mml:mo>
              <mml:mn>2</mml:mn>
            </mml:msubsup>
            <mml:mo>=</mml:mo>
            <mml:mo>−</mml:mo>
            <mml:mn>2</mml:mn>
            <mml:msub>
              <mml:mi>v</mml:mi>
              <mml:mo>⊕</mml:mo>
            </mml:msub>
            <mml:mi>u</mml:mi>
            <mml:mo>+</mml:mo>
            <mml:msup>
              <mml:mi>u</mml:mi>
              <mml:mn>2</mml:mn>
            </mml:msup>
            <mml:mo>.</mml:mo>
          </mml:mrow>
        </mml:math>
      </disp-formula>
      <p>For ordinary commercial aircraft <inline-formula><mml:math><mml:mrow><mml:mi> u </mml:mi><mml:mo> &lt; </mml:mo><mml:msub><mml:mi> v </mml:mi><mml:mo> ⊕ </mml:mo></mml:msub></mml:mrow></mml:math></inline-formula> , so the term <inline-formula><mml:math><mml:mrow><mml:mo> − </mml:mo><mml:mn> 2 </mml:mn><mml:msub><mml:mi> v </mml:mi><mml:mo> ⊕ </mml:mo></mml:msub><mml:mi> u </mml:mi></mml:mrow></mml:math></inline-formula> dominates over <inline-formula><mml:math><mml:mrow><mml:msup><mml:mi> u </mml:mi><mml:mn> 2 </mml:mn></mml:msup></mml:mrow></mml:math></inline-formula> . Consequently, the kinetic contribution is negative, corresponding to a reduction in the effective vacuum-density loading. Combined with the altitude-induced rarefaction, this yields a negative total density perturbation and therefore a positive clock shift relative to the ground reference clock.</p>
      <p>Using the clock-rate relation derived in Section 4, the absolute time shift accumulated over a total flight duration <inline-formula><mml:math display="inline"><mml:mi> T </mml:mi></mml:math></inline-formula> is</p>
      <disp-formula id="FD46">
        <label>(42)</label>
        <mml:math>
          <mml:mrow>
            <mml:mi>Δ</mml:mi>
            <mml:mi>τ</mml:mi>
            <mml:mo>=</mml:mo>
            <mml:mi>T</mml:mi>
            <mml:mrow>
              <mml:mo>[</mml:mo>
              <mml:mrow>
                <mml:mfrac>
                  <mml:mrow>
                    <mml:mi>g</mml:mi>
                    <mml:mi>h</mml:mi>
                  </mml:mrow>
                  <mml:mrow>
                    <mml:msup>
                      <mml:mi>c</mml:mi>
                      <mml:mn>2</mml:mn>
                    </mml:msup>
                  </mml:mrow>
                </mml:mfrac>
                <mml:mo>−</mml:mo>
                <mml:mfrac>
                  <mml:mrow>
                    <mml:msubsup>
                      <mml:mi>v</mml:mi>
                      <mml:mrow>
                        <mml:mtext>eff</mml:mtext>
                      </mml:mrow>
                      <mml:mn>2</mml:mn>
                    </mml:msubsup>
                    <mml:mo>−</mml:mo>
                    <mml:msubsup>
                      <mml:mi>v</mml:mi>
                      <mml:mo>⊕</mml:mo>
                      <mml:mn>2</mml:mn>
                    </mml:msubsup>
                  </mml:mrow>
                  <mml:mrow>
                    <mml:mn>2</mml:mn>
                    <mml:msup>
                      <mml:mi>c</mml:mi>
                      <mml:mn>2</mml:mn>
                    </mml:msup>
                  </mml:mrow>
                </mml:mfrac>
              </mml:mrow>
              <mml:mo>]</mml:mo>
            </mml:mrow>
            <mml:mo>.</mml:mo>
          </mml:mrow>
        </mml:math>
      </disp-formula>
      <p>For eastward and westward motion, Equation (42) becomes respectively</p>
      <disp-formula id="FD47">
        <label>(43)</label>
        <mml:math>
          <mml:mrow>
            <mml:mi>Δ</mml:mi>
            <mml:msub>
              <mml:mi>τ</mml:mi>
              <mml:mrow>
                <mml:mtext>east</mml:mtext>
              </mml:mrow>
            </mml:msub>
            <mml:mo>=</mml:mo>
            <mml:mi>T</mml:mi>
            <mml:mrow>
              <mml:mo>[</mml:mo>
              <mml:mrow>
                <mml:mfrac>
                  <mml:mrow>
                    <mml:mi>g</mml:mi>
                    <mml:mi>h</mml:mi>
                  </mml:mrow>
                  <mml:mrow>
                    <mml:msup>
                      <mml:mi>c</mml:mi>
                      <mml:mn>2</mml:mn>
                    </mml:msup>
                  </mml:mrow>
                </mml:mfrac>
                <mml:mo>−</mml:mo>
                <mml:mfrac>
                  <mml:mrow>
                    <mml:mn>2</mml:mn>
                    <mml:msub>
                      <mml:mi>v</mml:mi>
                      <mml:mo>⊕</mml:mo>
                    </mml:msub>
                    <mml:mi>u</mml:mi>
                    <mml:mo>+</mml:mo>
                    <mml:msup>
                      <mml:mi>u</mml:mi>
                      <mml:mn>2</mml:mn>
                    </mml:msup>
                  </mml:mrow>
                  <mml:mrow>
                    <mml:mn>2</mml:mn>
                    <mml:msup>
                      <mml:mi>c</mml:mi>
                      <mml:mn>2</mml:mn>
                    </mml:msup>
                  </mml:mrow>
                </mml:mfrac>
              </mml:mrow>
              <mml:mo>]</mml:mo>
            </mml:mrow>
            <mml:mo>.</mml:mo>
          </mml:mrow>
        </mml:math>
      </disp-formula>
      <disp-formula id="FD48">
        <label>(44)</label>
        <mml:math>
          <mml:mrow>
            <mml:mi>Δ</mml:mi>
            <mml:msub>
              <mml:mi>τ</mml:mi>
              <mml:mrow>
                <mml:mtext>west</mml:mtext>
              </mml:mrow>
            </mml:msub>
            <mml:mo>=</mml:mo>
            <mml:mi>T</mml:mi>
            <mml:mrow>
              <mml:mo>[</mml:mo>
              <mml:mrow>
                <mml:mfrac>
                  <mml:mrow>
                    <mml:mi>g</mml:mi>
                    <mml:mi>h</mml:mi>
                  </mml:mrow>
                  <mml:mrow>
                    <mml:msup>
                      <mml:mi>c</mml:mi>
                      <mml:mn>2</mml:mn>
                    </mml:msup>
                  </mml:mrow>
                </mml:mfrac>
                <mml:mo>−</mml:mo>
                <mml:mfrac>
                  <mml:mrow>
                    <mml:mo>−</mml:mo>
                    <mml:mn>2</mml:mn>
                    <mml:msub>
                      <mml:mi>v</mml:mi>
                      <mml:mo>⊕</mml:mo>
                    </mml:msub>
                    <mml:mi>u</mml:mi>
                    <mml:mo>+</mml:mo>
                    <mml:msup>
                      <mml:mi>u</mml:mi>
                      <mml:mn>2</mml:mn>
                    </mml:msup>
                  </mml:mrow>
                  <mml:mrow>
                    <mml:mn>2</mml:mn>
                    <mml:msup>
                      <mml:mi>c</mml:mi>
                      <mml:mn>2</mml:mn>
                    </mml:msup>
                  </mml:mrow>
                </mml:mfrac>
              </mml:mrow>
              <mml:mo>]</mml:mo>
            </mml:mrow>
            <mml:mo>.</mml:mo>
          </mml:mrow>
        </mml:math>
      </disp-formula>
      <p>Equations (43) and (44) show why the two directions have opposite signs. The altitude term is positive in the clock-shift expression for both directions, because a clock at higher altitude runs faster relative to the ground clock. The velocity term, however, is direction-dependent. Eastward motion increases the effective velocity and can overcome the altitude term, producing a net negative shift. Westward motion decreases the effective velocity and adds to the altitude effect, producing a net positive shift.</p>
      <p>Thus, the east-west asymmetry is not inserted by hand and is not obtained by fitting the observed time shifts. It follows from the squared-velocity dependence of the kinetic density perturbation and from the sign difference between motion with and against the Earth-coupled vacuum vortex. The Hafele-Keating measurements can therefore be used in the next section as an experimental comparison rather than as an input to the derivation.</p>
    </sec>
    <sec id="sec6">
      <title>6. Numerical Estimate for Eastward and Westward Flights</title>
      <p>The previous section derived the density perturbation and the corresponding first-order clock-shift equation without using the observed Hafele-Keating results as input. The present section applies the model numerically to representative around-the-world flight conditions. The purpose is not to reconstruct every segment of the original flight trajectories, but to show that the independently derived equation gives the correct sign and magnitude using physically reasonable averaged parameters.</p>
      <p>For a first-order estimate, the following effective mean values are used: total flight duration <italic>T</italic>, surface tangential velocity <inline-formula><mml:math display="inline"><mml:mrow><mml:msub><mml:mi> v </mml:mi><mml:mo> ⊕ </mml:mo></mml:msub></mml:mrow></mml:math></inline-formula> near the effective flight latitude, aircraft ground-relative velocity <italic>u</italic>, and mean altitude <italic>h</italic>:</p>
      <disp-formula id="FD49">
        <label>(45)</label>
        <mml:math>
          <mml:mrow>
            <mml:mi>T</mml:mi>
            <mml:mo>=</mml:mo>
            <mml:mn>144000</mml:mn>
            <mml:mtext>
               
            </mml:mtext>
            <mml:mtext>s</mml:mtext>
            <mml:mo>,</mml:mo>
            <mml:mtext>
               
            </mml:mtext>
            <mml:mtext>
               
            </mml:mtext>
            <mml:mtext>
               
            </mml:mtext>
            <mml:msub>
              <mml:mi>v</mml:mi>
              <mml:mo>⊕</mml:mo>
            </mml:msub>
            <mml:mo>=</mml:mo>
            <mml:mn>465</mml:mn>
            <mml:mtext>
               
            </mml:mtext>
            <mml:mtext>m</mml:mtext>
            <mml:mo>⋅</mml:mo>
            <mml:msup>
              <mml:mtext>s</mml:mtext>
              <mml:mrow>
                <mml:mo>−</mml:mo>
                <mml:mn>1</mml:mn>
              </mml:mrow>
            </mml:msup>
            <mml:mo>,</mml:mo>
            <mml:mtext>
               
            </mml:mtext>
            <mml:mtext>
               
            </mml:mtext>
            <mml:mtext>
               
            </mml:mtext>
            <mml:mi>u</mml:mi>
            <mml:mo>=</mml:mo>
            <mml:mn>223</mml:mn>
            <mml:mtext>
               
            </mml:mtext>
            <mml:mtext>m</mml:mtext>
            <mml:mo>⋅</mml:mo>
            <mml:msup>
              <mml:mtext>s</mml:mtext>
              <mml:mrow>
                <mml:mo>−</mml:mo>
                <mml:mn>1</mml:mn>
              </mml:mrow>
            </mml:msup>
            <mml:mo>,</mml:mo>
            <mml:mtext>
               
            </mml:mtext>
            <mml:mtext>
               
            </mml:mtext>
            <mml:mtext>
               
            </mml:mtext>
            <mml:mi>h</mml:mi>
            <mml:mo>=</mml:mo>
            <mml:mn>9340</mml:mn>
            <mml:mtext>
               
            </mml:mtext>
            <mml:mtext>m</mml:mtext>
            <mml:mo>.</mml:mo>
          </mml:mrow>
        </mml:math>
      </disp-formula>
      <p>The value of u should be interpreted as an effective mean ground-relative speed averaged over the route, including climb, descent, non-equatorial latitude, non-uniform airspeed, and operational delays. The altitude h should likewise be interpreted as an effective mean cruising altitude over the clock-transport interval rather than a fixed geometric altitude.</p>
      <p>The duration <inline-formula><mml:math><mml:mrow><mml:mi> T </mml:mi><mml:mo> = </mml:mo><mml:mn> 1.44 </mml:mn><mml:mo> × </mml:mo><mml:msup><mml:mrow><mml:mn> 10 </mml:mn></mml:mrow><mml:mn> 5 </mml:mn></mml:msup></mml:mrow></mml:math></inline-formula> employed here is not intended to represent the full elapsed time of either Hafele-Keating journey. Rather, it serves as an effective clock-transport duration adopted for this first-order estimate. A more rigorous treatment would require a segment-by-segment integration over the actual flight logs, incorporating time-dependent variations in airborne intervals, ground stopovers, altitude, latitude, and aircraft velocity along the complete trajectories.</p>
      <p>For eastward motion, the aircraft moves in the same sense as the Earth-coupled vacuum vortex. The effective velocity is therefore</p>
      <disp-formula id="FD50">
        <label>(46)</label>
        <mml:math>
          <mml:mrow>
            <mml:msub>
              <mml:mi>v</mml:mi>
              <mml:mrow>
                <mml:mtext>eff</mml:mtext>
                <mml:mo>,</mml:mo>
                <mml:mtext>east</mml:mtext>
              </mml:mrow>
            </mml:msub>
            <mml:mo>=</mml:mo>
            <mml:msub>
              <mml:mi>v</mml:mi>
              <mml:mo>⊕</mml:mo>
            </mml:msub>
            <mml:mo>+</mml:mo>
            <mml:mi>u</mml:mi>
            <mml:mo>=</mml:mo>
            <mml:mn>688</mml:mn>
            <mml:mtext>
               
            </mml:mtext>
            <mml:mtext>m</mml:mtext>
            <mml:mo>⋅</mml:mo>
            <mml:msup>
              <mml:mtext>s</mml:mtext>
              <mml:mrow>
                <mml:mo>−</mml:mo>
                <mml:mn>1</mml:mn>
              </mml:mrow>
            </mml:msup>
            <mml:mo>.</mml:mo>
          </mml:mrow>
        </mml:math>
      </disp-formula>
      <p>For westward motion, the aircraft moves against the rotational flow of the Earth-coupled vortex. The effective velocity is therefore</p>
      <disp-formula id="FD51">
        <label>(47)</label>
        <mml:math>
          <mml:mrow>
            <mml:msub>
              <mml:mi>v</mml:mi>
              <mml:mrow>
                <mml:mtext>eff</mml:mtext>
                <mml:mo>,</mml:mo>
                <mml:mtext>west</mml:mtext>
              </mml:mrow>
            </mml:msub>
            <mml:mo>=</mml:mo>
            <mml:msub>
              <mml:mi>v</mml:mi>
              <mml:mo>⊕</mml:mo>
            </mml:msub>
            <mml:mo>−</mml:mo>
            <mml:mi>u</mml:mi>
            <mml:mo>=</mml:mo>
            <mml:mn>242</mml:mn>
            <mml:mtext>
               
            </mml:mtext>
            <mml:mtext>m</mml:mtext>
            <mml:mo>⋅</mml:mo>
            <mml:msup>
              <mml:mtext>s</mml:mtext>
              <mml:mrow>
                <mml:mo>−</mml:mo>
                <mml:mn>1</mml:mn>
              </mml:mrow>
            </mml:msup>
            <mml:mo>.</mml:mo>
          </mml:mrow>
        </mml:math>
      </disp-formula>
      <sec id="sec6dot1">
        <title>6.1. Density Perturbations</title>
        <p>Substituting the eastward effective velocity into the independently derived density-perturbation equation gives</p>
        <disp-formula id="FD52">
          <label>(48)</label>
          <mml:math>
            <mml:mrow>
              <mml:msub>
                <mml:mrow>
                  <mml:mrow>
                    <mml:mo>(</mml:mo>
                    <mml:mrow>
                      <mml:mfrac>
                        <mml:mrow>
                          <mml:mi>Δ</mml:mi>
                          <mml:mi>ρ</mml:mi>
                        </mml:mrow>
                        <mml:mrow>
                          <mml:msub>
                            <mml:mi>ρ</mml:mi>
                            <mml:mrow>
                              <mml:mtext>eff</mml:mtext>
                            </mml:mrow>
                          </mml:msub>
                        </mml:mrow>
                      </mml:mfrac>
                    </mml:mrow>
                    <mml:mo>)</mml:mo>
                  </mml:mrow>
                </mml:mrow>
                <mml:mrow>
                  <mml:mtext>east</mml:mtext>
                </mml:mrow>
              </mml:msub>
              <mml:mo>≃</mml:mo>
              <mml:mo>+</mml:mo>
              <mml:mn>8.22</mml:mn>
              <mml:mo>×</mml:mo>
              <mml:msup>
                <mml:mrow>
                  <mml:mn>10</mml:mn>
                </mml:mrow>
                <mml:mrow>
                  <mml:mo>−</mml:mo>
                  <mml:mn>13</mml:mn>
                </mml:mrow>
              </mml:msup>
              <mml:mo>.</mml:mo>
            </mml:mrow>
          </mml:math>
        </disp-formula>
        <p>The positive sign indicates a slight kinetic compression or loading of the local effective vacuum-density state. According to the wave-speed relation developed above, this positive density perturbation reduces the effective electromagnetic propagation rate and therefore slows the airborne clock relative to the ground reference clock.</p>
        <p>For westward motion, the corresponding calculation gives</p>
        <disp-formula id="FD53">
          <label>(49)</label>
          <mml:math>
            <mml:mrow>
              <mml:msub>
                <mml:mrow>
                  <mml:mrow>
                    <mml:mo>(</mml:mo>
                    <mml:mrow>
                      <mml:mfrac>
                        <mml:mrow>
                          <mml:mi>Δ</mml:mi>
                          <mml:mi>ρ</mml:mi>
                        </mml:mrow>
                        <mml:mrow>
                          <mml:msub>
                            <mml:mi>ρ</mml:mi>
                            <mml:mrow>
                              <mml:mtext>eff</mml:mtext>
                            </mml:mrow>
                          </mml:msub>
                        </mml:mrow>
                      </mml:mfrac>
                    </mml:mrow>
                    <mml:mo>)</mml:mo>
                  </mml:mrow>
                </mml:mrow>
                <mml:mrow>
                  <mml:mtext>west</mml:mtext>
                </mml:mrow>
              </mml:msub>
              <mml:mo>≃</mml:mo>
              <mml:mo>−</mml:mo>
              <mml:mn>3.79</mml:mn>
              <mml:mo>×</mml:mo>
              <mml:msup>
                <mml:mrow>
                  <mml:mn>10</mml:mn>
                </mml:mrow>
                <mml:mrow>
                  <mml:mo>−</mml:mo>
                  <mml:mn>12</mml:mn>
                </mml:mrow>
              </mml:msup>
              <mml:mo>.</mml:mo>
            </mml:mrow>
          </mml:math>
        </disp-formula>
        <p>The negative sign indicates an effective rarefaction of the local vacuum-density state relative to the ground reference clock. This occurs because westward motion partially cancels the rotational velocity associated with the Earth-coupled vacuum vortex. The result is a faster clock rate relative to the surface reference clock.</p>
      </sec>
      <sec id="sec6dot2">
        <title>6.2. Predicted Clock Shifts</title>
        <p>Using the first-order clock-shift equation derived in Section 5,</p>
        <disp-formula id="FD54">
          <label>(50)</label>
          <mml:math>
            <mml:mrow>
              <mml:mi>Δ</mml:mi>
              <mml:mi>τ</mml:mi>
              <mml:mo>=</mml:mo>
              <mml:mi>T</mml:mi>
              <mml:mrow>
                <mml:mo>[</mml:mo>
                <mml:mrow>
                  <mml:mfrac>
                    <mml:mrow>
                      <mml:mi>g</mml:mi>
                      <mml:mi>h</mml:mi>
                    </mml:mrow>
                    <mml:mrow>
                      <mml:msup>
                        <mml:mi>c</mml:mi>
                        <mml:mn>2</mml:mn>
                      </mml:msup>
                    </mml:mrow>
                  </mml:mfrac>
                  <mml:mo>−</mml:mo>
                  <mml:mfrac>
                    <mml:mrow>
                      <mml:msubsup>
                        <mml:mi>v</mml:mi>
                        <mml:mrow>
                          <mml:mtext>eff</mml:mtext>
                        </mml:mrow>
                        <mml:mn>2</mml:mn>
                      </mml:msubsup>
                      <mml:mo>−</mml:mo>
                      <mml:msubsup>
                        <mml:mi>v</mml:mi>
                        <mml:mo>⊕</mml:mo>
                        <mml:mn>2</mml:mn>
                      </mml:msubsup>
                    </mml:mrow>
                    <mml:mrow>
                      <mml:mn>2</mml:mn>
                      <mml:msup>
                        <mml:mi>c</mml:mi>
                        <mml:mn>2</mml:mn>
                      </mml:msup>
                    </mml:mrow>
                  </mml:mfrac>
                </mml:mrow>
                <mml:mo>]</mml:mo>
              </mml:mrow>
              <mml:mo>.</mml:mo>
            </mml:mrow>
          </mml:math>
        </disp-formula>
        <p>the eastward flight gives</p>
        <disp-formula id="FD55">
          <label>(51)</label>
          <mml:math>
            <mml:mrow>
              <mml:mi>Δ</mml:mi>
              <mml:msub>
                <mml:mi>τ</mml:mi>
                <mml:mrow>
                  <mml:mtext>east</mml:mtext>
                </mml:mrow>
              </mml:msub>
              <mml:mo>≃</mml:mo>
              <mml:mo>−</mml:mo>
              <mml:mn>5.92</mml:mn>
              <mml:mo>×</mml:mo>
              <mml:msup>
                <mml:mrow>
                  <mml:mn>10</mml:mn>
                </mml:mrow>
                <mml:mrow>
                  <mml:mo>−</mml:mo>
                  <mml:mn>8</mml:mn>
                </mml:mrow>
              </mml:msup>
              <mml:mtext>
                 
              </mml:mtext>
              <mml:mtext>s</mml:mtext>
              <mml:mo>=</mml:mo>
              <mml:mo>−</mml:mo>
              <mml:mn>59.2</mml:mn>
              <mml:mtext>
                 
              </mml:mtext>
              <mml:mtext>ns</mml:mtext>
              <mml:mo>.</mml:mo>
            </mml:mrow>
          </mml:math>
        </disp-formula>
        <p>The predicted eastward value is therefore essentially the same as the reported Hafele-Keating eastward result of approximately −59 ns.</p>
        <p>For the westward flight, the same equation gives</p>
        <disp-formula id="FD56">
          <label>(52)</label>
          <mml:math>
            <mml:mrow>
              <mml:mi>Δ</mml:mi>
              <mml:msub>
                <mml:mi>τ</mml:mi>
                <mml:mrow>
                  <mml:mtext>west</mml:mtext>
                </mml:mrow>
              </mml:msub>
              <mml:mo>≃</mml:mo>
              <mml:mo>+</mml:mo>
              <mml:mn>2.73</mml:mn>
              <mml:mo>×</mml:mo>
              <mml:msup>
                <mml:mrow>
                  <mml:mn>10</mml:mn>
                </mml:mrow>
                <mml:mrow>
                  <mml:mo>−</mml:mo>
                  <mml:mn>7</mml:mn>
                </mml:mrow>
              </mml:msup>
              <mml:mtext>
                 
              </mml:mtext>
              <mml:mtext>s</mml:mtext>
              <mml:mo>=</mml:mo>
              <mml:mo>+</mml:mo>
              <mml:mn>273.1</mml:mn>
              <mml:mtext>
                 
              </mml:mtext>
              <mml:mtext>ns</mml:mtext>
              <mml:mo>.</mml:mo>
            </mml:mrow>
          </mml:math>
        </disp-formula>
        <p>The predicted westward value is therefore essentially the same as the reported Hafele-Keating westward result of approximately +273 ns.</p>
        <p>The predicted difference between the westward and eastward clock shifts is</p>
        <disp-formula id="FD57">
          <label>(53)</label>
          <mml:math>
            <mml:mrow>
              <mml:mi>Δ</mml:mi>
              <mml:msub>
                <mml:mi>τ</mml:mi>
                <mml:mrow>
                  <mml:mtext>west</mml:mtext>
                </mml:mrow>
              </mml:msub>
              <mml:mo>−</mml:mo>
              <mml:mi>Δ</mml:mi>
              <mml:msub>
                <mml:mi>τ</mml:mi>
                <mml:mrow>
                  <mml:mtext>east</mml:mtext>
                </mml:mrow>
              </mml:msub>
              <mml:mo>≃</mml:mo>
              <mml:mn>332.3</mml:mn>
              <mml:mtext>
                 
              </mml:mtext>
              <mml:mtext>ns</mml:mtext>
              <mml:mo>.</mml:mo>
            </mml:mrow>
          </mml:math>
        </disp-formula>
        <p>This is the same east-west asymmetry measured in the Hafele-Keating experiment to first-order accuracy.</p>
        <table-wrap id="tbl1">
          <label>Table 1</label>
          <table>
            <tbody>
              <tr>
                <td>
                  <bold>Direction</bold>
                </td>
                <td>
                  <bold>Vortex-density</bold>
                  <bold>estimate</bold>
                </td>
                <td>
                  <bold>Hafele-Keating</bold>
                  <bold>measured</bold>
                  <bold>value</bold>
                </td>
                <td>
                  <bold>Agreement</bold>
                </td>
              </tr>
              <tr>
                <td>Eastward</td>
                <td>−59.2 ns</td>
                <td>−59 ns</td>
                <td>Correct sign and magnitude</td>
              </tr>
              <tr>
                <td>Westward</td>
                <td>+273.1 ns</td>
                <td>+273 ns</td>
                <td>Correct sign and magnitude</td>
              </tr>
              <tr>
                <td>West-East difference</td>
                <td>332.3 ns</td>
                <td>332 ns</td>
                <td>Correct east-west asymmetry</td>
              </tr>
            </tbody>
          </table>
        </table-wrap>
        <p>The close numerical agreement should not be interpreted as a precision trajectory-level prediction. The calculation uses effective averaged values for latitude, altitude, aircraft velocity, and clock-transport duration. Therefore, the result should be understood as a first-order consistency estimate showing that the model reproduces the correct sign, order of magnitude, and approximate east-west asymmetry of the Hafele-Keating measurements. A stronger test would require reconstructing the original flight paths segment by segment.</p>
      </sec>
      <sec id="sec6dot3">
        <title>6.3. Sensitivity to Averaged Flight Parameters</title>
        <p>The agreement depends on the use of effective mean parameters. This is expected because the original flights did not occur at a single velocity, altitude, or latitude. The model is sensitive to the altitude through the common positive clock-rate contribution</p>
        <disp-formula id="FD58">
          <label>(54)</label>
          <mml:math>
            <mml:mrow>
              <mml:mfrac>
                <mml:mrow>
                  <mml:mo>∂</mml:mo>
                  <mml:mi>Δ</mml:mi>
                  <mml:mi>τ</mml:mi>
                </mml:mrow>
                <mml:mrow>
                  <mml:mo>∂</mml:mo>
                  <mml:mi>h</mml:mi>
                </mml:mrow>
              </mml:mfrac>
              <mml:mo>=</mml:mo>
              <mml:mfrac>
                <mml:mrow>
                  <mml:mi>T</mml:mi>
                  <mml:mi>g</mml:mi>
                </mml:mrow>
                <mml:mrow>
                  <mml:msup>
                    <mml:mi>c</mml:mi>
                    <mml:mn>2</mml:mn>
                  </mml:msup>
                </mml:mrow>
              </mml:mfrac>
              <mml:mo>≃</mml:mo>
              <mml:mn>15.7</mml:mn>
              <mml:mtext>
                 
              </mml:mtext>
              <mml:mtext>ns</mml:mtext>
              <mml:mo>⋅</mml:mo>
              <mml:msup>
                <mml:mrow>
                  <mml:mtext>km</mml:mtext>
                </mml:mrow>
                <mml:mrow>
                  <mml:mo>−</mml:mo>
                  <mml:mn>1</mml:mn>
                </mml:mrow>
              </mml:msup>
              <mml:mo>.</mml:mo>
            </mml:mrow>
          </mml:math>
        </disp-formula>
        <p>Thus, an altitude change of only about 1 km changes the predicted clock shift by approximately 16 ns in both directions. This explains why using 10 km instead of 9.34 km produces a common offset of the order of tens of nanoseconds.</p>
        <p>The east-west separation is mainly controlled by the aircraft speed through</p>
        <disp-formula id="FD59">
          <label>(55)</label>
          <mml:math>
            <mml:mrow>
              <mml:mfrac>
                <mml:mrow>
                  <mml:mo>∂</mml:mo>
                  <mml:mrow>
                    <mml:mo>(</mml:mo>
                    <mml:mrow>
                      <mml:mi>Δ</mml:mi>
                      <mml:msub>
                        <mml:mi>τ</mml:mi>
                        <mml:mrow>
                          <mml:mtext>west</mml:mtext>
                        </mml:mrow>
                      </mml:msub>
                      <mml:mo>−</mml:mo>
                      <mml:mi>Δ</mml:mi>
                      <mml:msub>
                        <mml:mi>τ</mml:mi>
                        <mml:mrow>
                          <mml:mtext>east</mml:mtext>
                        </mml:mrow>
                      </mml:msub>
                    </mml:mrow>
                    <mml:mo>)</mml:mo>
                  </mml:mrow>
                </mml:mrow>
                <mml:mrow>
                  <mml:mo>∂</mml:mo>
                  <mml:mi>u</mml:mi>
                </mml:mrow>
              </mml:mfrac>
              <mml:mo>=</mml:mo>
              <mml:mfrac>
                <mml:mrow>
                  <mml:mn>2</mml:mn>
                  <mml:mi>T</mml:mi>
                  <mml:msub>
                    <mml:mi>v</mml:mi>
                    <mml:mo>⊕</mml:mo>
                  </mml:msub>
                </mml:mrow>
                <mml:mrow>
                  <mml:msup>
                    <mml:mi>c</mml:mi>
                    <mml:mn>2</mml:mn>
                  </mml:msup>
                </mml:mrow>
              </mml:mfrac>
              <mml:mo>≃</mml:mo>
              <mml:mn>1.49</mml:mn>
              <mml:mtext>
                 
              </mml:mtext>
              <mml:mtext>ns</mml:mtext>
              <mml:mo>⋅</mml:mo>
              <mml:msup>
                <mml:mrow>
                  <mml:mrow>
                    <mml:mo>(</mml:mo>
                    <mml:mrow>
                      <mml:mtext>m</mml:mtext>
                      <mml:mo>⋅</mml:mo>
                      <mml:msup>
                        <mml:mtext>s</mml:mtext>
                        <mml:mrow>
                          <mml:mo>−</mml:mo>
                          <mml:mn>1</mml:mn>
                        </mml:mrow>
                      </mml:msup>
                    </mml:mrow>
                    <mml:mo>)</mml:mo>
                  </mml:mrow>
                </mml:mrow>
                <mml:mrow>
                  <mml:mo>−</mml:mo>
                  <mml:mn>1</mml:mn>
                </mml:mrow>
              </mml:msup>
              <mml:mo>.</mml:mo>
            </mml:mrow>
          </mml:math>
        </disp-formula>
        <p>Thus, a change of only 10 to 15 m/s in the effective mean aircraft speed changes the east-west separation by approximately 15 to 22 ns. This explains why approximate values such as <italic>u</italic> = 250 m/s and <italic>h</italic> = 10 km produce a small overestimate, while more realistic effective means recover the measured values.</p>
        <p>These sensitivities show that the agreement is not obtained by introducing a new arbitrary density fitted to the observed clock shifts. Rather, the clock shifts follow from the independently derived first-order equation once realistic averaged flight parameters are used. A still stronger validation would require reconstructing the original flight paths segment by segment, including latitude, altitude, velocity, and stopover times. That trajectory-level calculation is beyond the scope of the present first-order analysis, but it represents a natural next step for testing the model more rigorously.</p>
      </sec>
    </sec>
    <sec id="sec7">
      <title>7. Comparison with the Hafele-Keating Measurements</title>
      <p>The previous section applied the independently derived vortex-density equation to representative mean flight parameters. The next step is to compare the resulting first-order clock shifts with the values reported in the Hafele-Keating experiment. This comparison is essential because the proposed model is not intended merely to reproduce the sign of the effect; it must also reproduce the observed order of magnitude and the east-west asymmetry.</p>
      <p>The experimentally reported clock shifts may be written as</p>
      <disp-formula id="FD60">
        <label>(56)</label>
        <mml:math>
          <mml:mrow>
            <mml:mi>Δ</mml:mi>
            <mml:msub>
              <mml:mi>τ</mml:mi>
              <mml:mrow>
                <mml:mtext>HK</mml:mtext>
                <mml:mo>,</mml:mo>
                <mml:mtext>east</mml:mtext>
              </mml:mrow>
            </mml:msub>
            <mml:mo>=</mml:mo>
            <mml:mo>−</mml:mo>
            <mml:mn>59</mml:mn>
            <mml:mtext>
               
            </mml:mtext>
            <mml:mtext>ns</mml:mtext>
            <mml:mo>,</mml:mo>
            <mml:mtext>
               
            </mml:mtext>
            <mml:mtext>
               
            </mml:mtext>
            <mml:mtext>
               
            </mml:mtext>
            <mml:mi>Δ</mml:mi>
            <mml:msub>
              <mml:mi>τ</mml:mi>
              <mml:mrow>
                <mml:mtext>HK</mml:mtext>
                <mml:mo>,</mml:mo>
                <mml:mtext>west</mml:mtext>
              </mml:mrow>
            </mml:msub>
            <mml:mo>=</mml:mo>
            <mml:mo>+</mml:mo>
            <mml:mn>273</mml:mn>
            <mml:mtext>
               
            </mml:mtext>
            <mml:mtext>ns</mml:mtext>
            <mml:mo>.</mml:mo>
          </mml:mrow>
        </mml:math>
      </disp-formula>
      <p>Using the effective mean parameters introduced in Section 6, the vortex-density model gives</p>
      <disp-formula id="FD61">
        <label>(57)</label>
        <mml:math>
          <mml:mrow>
            <mml:mi>Δ</mml:mi>
            <mml:msub>
              <mml:mi>τ</mml:mi>
              <mml:mrow>
                <mml:mtext>V</mml:mtext>
                <mml:mo>,</mml:mo>
                <mml:mtext>east</mml:mtext>
              </mml:mrow>
            </mml:msub>
            <mml:mo>≃</mml:mo>
            <mml:mo>−</mml:mo>
            <mml:mn>59.2</mml:mn>
            <mml:mtext>
               
            </mml:mtext>
            <mml:mtext>ns</mml:mtext>
            <mml:mo>,</mml:mo>
            <mml:mtext>
               
            </mml:mtext>
            <mml:mtext>
               
            </mml:mtext>
            <mml:mtext>
               
            </mml:mtext>
            <mml:mi>Δ</mml:mi>
            <mml:msub>
              <mml:mi>τ</mml:mi>
              <mml:mrow>
                <mml:mtext>V</mml:mtext>
                <mml:mo>,</mml:mo>
                <mml:mtext>west</mml:mtext>
              </mml:mrow>
            </mml:msub>
            <mml:mo>≃</mml:mo>
            <mml:mo>+</mml:mo>
            <mml:mn>273.1</mml:mn>
            <mml:mtext>
               
            </mml:mtext>
            <mml:mtext>ns</mml:mtext>
            <mml:mo>.</mml:mo>
          </mml:mrow>
        </mml:math>
      </disp-formula>
      <p>The residual difference between the model estimate and the reported measurement is therefore</p>
      <disp-formula id="FD62">
        <label>(58)</label>
        <mml:math>
          <mml:mrow>
            <mml:msub>
              <mml:mi>ϵ</mml:mi>
              <mml:mrow>
                <mml:mtext>east</mml:mtext>
              </mml:mrow>
            </mml:msub>
            <mml:mo>=</mml:mo>
            <mml:mi>Δ</mml:mi>
            <mml:msub>
              <mml:mi>τ</mml:mi>
              <mml:mrow>
                <mml:mtext>V</mml:mtext>
                <mml:mo>,</mml:mo>
                <mml:mtext>east</mml:mtext>
              </mml:mrow>
            </mml:msub>
            <mml:mo>−</mml:mo>
            <mml:mi>Δ</mml:mi>
            <mml:msub>
              <mml:mi>τ</mml:mi>
              <mml:mrow>
                <mml:mtext>HK</mml:mtext>
                <mml:mo>,</mml:mo>
                <mml:mtext>east</mml:mtext>
              </mml:mrow>
            </mml:msub>
            <mml:mo>≃</mml:mo>
            <mml:mo>−</mml:mo>
            <mml:mn>0.2</mml:mn>
            <mml:mtext>
               
            </mml:mtext>
            <mml:mtext>ns</mml:mtext>
            <mml:mo>,</mml:mo>
          </mml:mrow>
        </mml:math>
      </disp-formula>
      <p>And</p>
      <disp-formula id="FD63">
        <label>(59)</label>
        <mml:math>
          <mml:mrow>
            <mml:msub>
              <mml:mi>ϵ</mml:mi>
              <mml:mrow>
                <mml:mtext>west</mml:mtext>
              </mml:mrow>
            </mml:msub>
            <mml:mo>=</mml:mo>
            <mml:mi>Δ</mml:mi>
            <mml:msub>
              <mml:mi>τ</mml:mi>
              <mml:mrow>
                <mml:mtext>V</mml:mtext>
                <mml:mo>,</mml:mo>
                <mml:mtext>west</mml:mtext>
              </mml:mrow>
            </mml:msub>
            <mml:mo>−</mml:mo>
            <mml:mi>Δ</mml:mi>
            <mml:msub>
              <mml:mi>τ</mml:mi>
              <mml:mrow>
                <mml:mtext>HK</mml:mtext>
                <mml:mo>,</mml:mo>
                <mml:mtext>west</mml:mtext>
              </mml:mrow>
            </mml:msub>
            <mml:mo>≃</mml:mo>
            <mml:mo>+</mml:mo>
            <mml:mn>0.1</mml:mn>
            <mml:mtext>
               
            </mml:mtext>
            <mml:mtext>ns</mml:mtext>
            <mml:mo>.</mml:mo>
          </mml:mrow>
        </mml:math>
      </disp-formula>
      <p>Thus, the first-order comparison gives the following values: eastward, Hafele-Keating measurement = -59 ns, vortex-density estimate = −59.2 ns, residual = −0.2 ns; westward, Hafele-Keating measurement = +273 ns, vortex-density estimate = +273.1 ns, residual = +0.1 ns. The west-east asymmetry is 332 ns experimentally and approximately 332.3 ns in the vortex-density estimate.</p>
      <p>The agreement should be interpreted carefully. It does not mean that the model has replaced the full relativistic trajectory analysis of the Hafele-Keating experiment. The calculation uses effective mean values for altitude, aircraft velocity, latitude-dependent rotational velocity, and total flight duration. A rigorous reconstruction would require the complete flight trajectories, including changes in altitude, route latitude, aircraft speed, stopovers, and clock-comparison procedures.</p>
      <p>Nevertheless, the result is important for three reasons. First, the signs are obtained correctly: eastward motion leads to a clock loss, while westward motion leads to a clock gain. Second, the magnitude of the effect is reproduced at the nanosecond scale using physically measurable flight parameters rather than by inserting the observed clock shifts into the density calculation. Third, the east-west asymmetry follows naturally from the term involving the difference between the squared effective velocity of the airborne clock and that of the ground reference clock.</p>
      <p>This asymmetry can be written as</p>
      <disp-formula id="FD64">
        <label>(60)</label>
        <mml:math>
          <mml:mrow>
            <mml:mi>Δ</mml:mi>
            <mml:msub>
              <mml:mi>τ</mml:mi>
              <mml:mrow>
                <mml:mtext>west</mml:mtext>
              </mml:mrow>
            </mml:msub>
            <mml:mo>−</mml:mo>
            <mml:mi>Δ</mml:mi>
            <mml:msub>
              <mml:mi>τ</mml:mi>
              <mml:mrow>
                <mml:mtext>east</mml:mtext>
              </mml:mrow>
            </mml:msub>
            <mml:mo>≃</mml:mo>
            <mml:mfrac>
              <mml:mrow>
                <mml:mn>2</mml:mn>
                <mml:mi>T</mml:mi>
                <mml:msub>
                  <mml:mi>v</mml:mi>
                  <mml:mo>⊕</mml:mo>
                </mml:msub>
                <mml:mi>u</mml:mi>
              </mml:mrow>
              <mml:mrow>
                <mml:msup>
                  <mml:mi>c</mml:mi>
                  <mml:mn>2</mml:mn>
                </mml:msup>
              </mml:mrow>
            </mml:mfrac>
            <mml:mo>.</mml:mo>
          </mml:mrow>
        </mml:math>
      </disp-formula>
      <p>Equation (61) shows that the separation between westward and eastward clock shifts is controlled primarily by three measurable quantities: the flight duration, the rotational velocity of the Earth at the relevant latitude, and the aircraft velocity. The altitude contribution shifts both directions in the same positive direction, whereas the velocity contribution is direction-dependent. This is why altitude mainly controls the common offset of the two clock shifts, while aircraft speed mainly controls the east-west separation.</p>
      <p>In the present interpretation, the eastward aircraft has a larger effective velocity relative to the Earth-coupled vacuum vortex. This produces a positive kinetic density perturbation, slows the effective electromagnetic propagation rate, and leads to clock retardation. The westward aircraft partially cancels the rotational velocity of the ground reference frame, producing a negative kinetic density perturbation relative to the ground state. Combined with the altitude-related rarefaction, this leads to a net clock gain.</p>
      <p>The comparison therefore supports the internal consistency of the proposed mechanism: the same density-perturbation equation accounts for the sign, relative asymmetry, and approximate magnitude of the Hafele-Keating measurements. The result should be regarded as a first-order consistency test of the structured-vacuum interpretation, not as a final proof. The next step would be to perform a trajectory-resolved calculation using the original flight logs and clock-comparison data.</p>
    </sec>
    <sec id="sec8">
      <title>8. Relationship between the Vacuum-Density Interpretation and Standard Relativity</title>
      <p>The comparison with the Hafele-Keating measurements shows that the proposed structured-vacuum model reproduces the first-order sign and magnitude of the observed clock shifts. It is therefore necessary to clarify the relationship between this interpretation and standard relativity. The purpose of the present framework is not to reject the relativistic result, but to provide a possible physical-substrate interpretation for the same weak-field clock-shift structure.</p>
      <p>In standard relativity, the first-order difference between an airborne clock and a ground reference clock may be written as</p>
      <disp-formula id="FD65">
        <label>(61)</label>
        <mml:math>
          <mml:mrow>
            <mml:mfrac>
              <mml:mrow>
                <mml:mi>Δ</mml:mi>
                <mml:mi>τ</mml:mi>
              </mml:mrow>
              <mml:mi>T</mml:mi>
            </mml:mfrac>
            <mml:mo>≃</mml:mo>
            <mml:mfrac>
              <mml:mrow>
                <mml:mi>Δ</mml:mi>
                <mml:mi>Φ</mml:mi>
              </mml:mrow>
              <mml:mrow>
                <mml:msup>
                  <mml:mi>c</mml:mi>
                  <mml:mn>2</mml:mn>
                </mml:msup>
              </mml:mrow>
            </mml:mfrac>
            <mml:mo>−</mml:mo>
            <mml:mfrac>
              <mml:mrow>
                <mml:msubsup>
                  <mml:mi>v</mml:mi>
                  <mml:mrow>
                    <mml:mtext>air</mml:mtext>
                  </mml:mrow>
                  <mml:mn>2</mml:mn>
                </mml:msubsup>
                <mml:mo>−</mml:mo>
                <mml:msubsup>
                  <mml:mi>v</mml:mi>
                  <mml:mrow>
                    <mml:mtext>ground</mml:mtext>
                  </mml:mrow>
                  <mml:mn>2</mml:mn>
                </mml:msubsup>
              </mml:mrow>
              <mml:mrow>
                <mml:mn>2</mml:mn>
                <mml:msup>
                  <mml:mi>c</mml:mi>
                  <mml:mn>2</mml:mn>
                </mml:msup>
              </mml:mrow>
            </mml:mfrac>
            <mml:mo>.</mml:mo>
          </mml:mrow>
        </mml:math>
      </disp-formula>
      <p>where <inline-formula><mml:math display="inline"><mml:mrow><mml:mi> Δ </mml:mi><mml:mi> Φ </mml:mi></mml:mrow></mml:math></inline-formula> is the gravitational-potential difference between the airborne clock and the ground clock, <inline-formula><mml:math display="inline"><mml:mrow><mml:msub><mml:mi> v </mml:mi><mml:mrow><mml:mtext> air </mml:mtext></mml:mrow></mml:msub></mml:mrow></mml:math></inline-formula> is the velocity of the aircraft clock relative to the chosen Earth-centered inertial frame, and <inline-formula><mml:math display="inline"><mml:mrow><mml:msub><mml:mi> v </mml:mi><mml:mrow><mml:mtext> ground </mml:mtext></mml:mrow></mml:msub></mml:mrow></mml:math></inline-formula> is the corresponding velocity of the ground reference clock. For small altitudes above the Earth’s surface, the potential difference is approximated by</p>
      <disp-formula id="FD66">
        <label>(62)</label>
        <mml:math>
          <mml:mrow>
            <mml:mi>Δ</mml:mi>
            <mml:mi>Φ</mml:mi>
            <mml:mo>≃</mml:mo>
            <mml:mi>g</mml:mi>
            <mml:mi>h</mml:mi>
            <mml:mo>.</mml:mo>
          </mml:mrow>
        </mml:math>
      </disp-formula>
      <p>Substitution gives the familiar weak-field expression</p>
      <disp-formula id="FD67">
        <label>(63)</label>
        <mml:math>
          <mml:mrow>
            <mml:mfrac>
              <mml:mrow>
                <mml:mi>Δ</mml:mi>
                <mml:mi>τ</mml:mi>
              </mml:mrow>
              <mml:mi>T</mml:mi>
            </mml:mfrac>
            <mml:mo>≃</mml:mo>
            <mml:mfrac>
              <mml:mrow>
                <mml:mi>g</mml:mi>
                <mml:mi>h</mml:mi>
              </mml:mrow>
              <mml:mrow>
                <mml:msup>
                  <mml:mi>c</mml:mi>
                  <mml:mn>2</mml:mn>
                </mml:msup>
              </mml:mrow>
            </mml:mfrac>
            <mml:mo>−</mml:mo>
            <mml:mfrac>
              <mml:mrow>
                <mml:msubsup>
                  <mml:mi>v</mml:mi>
                  <mml:mrow>
                    <mml:mtext>eff</mml:mtext>
                  </mml:mrow>
                  <mml:mn>2</mml:mn>
                </mml:msubsup>
                <mml:mo>−</mml:mo>
                <mml:msubsup>
                  <mml:mi>v</mml:mi>
                  <mml:mo>⊕</mml:mo>
                  <mml:mn>2</mml:mn>
                </mml:msubsup>
              </mml:mrow>
              <mml:mrow>
                <mml:mn>2</mml:mn>
                <mml:msup>
                  <mml:mi>c</mml:mi>
                  <mml:mn>2</mml:mn>
                </mml:msup>
              </mml:mrow>
            </mml:mfrac>
            <mml:mo>.</mml:mo>
          </mml:mrow>
        </mml:math>
      </disp-formula>
      <p>where <inline-formula><mml:math display="inline"><mml:mrow><mml:msub><mml:mi> v </mml:mi><mml:mrow><mml:mtext> eff </mml:mtext></mml:mrow></mml:msub></mml:mrow></mml:math></inline-formula> corresponds to the effective velocity of the aircraft clock and <inline-formula><mml:math display="inline"><mml:mrow><mml:msub><mml:mi> v </mml:mi><mml:mo> ⊕ </mml:mo></mml:msub></mml:mrow></mml:math></inline-formula> to the tangential velocity of the ground clock at the relevant latitude. This is the expression used implicitly in the conventional interpretation of the Hafele-Keating experiment.</p>
      <p>In the vacuum-density interpretation, the clock-rate shift is written instead as</p>
      <disp-formula id="FD68">
        <label>(64)</label>
        <mml:math>
          <mml:mrow>
            <mml:mfrac>
              <mml:mrow>
                <mml:mi>Δ</mml:mi>
                <mml:mi>τ</mml:mi>
              </mml:mrow>
              <mml:mi>T</mml:mi>
            </mml:mfrac>
            <mml:mo>≃</mml:mo>
            <mml:mo>−</mml:mo>
            <mml:mfrac>
              <mml:mn>1</mml:mn>
              <mml:mn>2</mml:mn>
            </mml:mfrac>
            <mml:mfrac>
              <mml:mrow>
                <mml:mi>Δ</mml:mi>
                <mml:mi>ρ</mml:mi>
              </mml:mrow>
              <mml:mrow>
                <mml:msub>
                  <mml:mi>ρ</mml:mi>
                  <mml:mrow>
                    <mml:mtext>eff</mml:mtext>
                  </mml:mrow>
                </mml:msub>
              </mml:mrow>
            </mml:mfrac>
            <mml:mo>.</mml:mo>
          </mml:mrow>
        </mml:math>
      </disp-formula>
      <p>with the density perturbation derived from vortex kinematics and altitude as</p>
      <disp-formula id="FD69">
        <label>(65)</label>
        <mml:math>
          <mml:mrow>
            <mml:mfrac>
              <mml:mrow>
                <mml:mi>Δ</mml:mi>
                <mml:mi>ρ</mml:mi>
              </mml:mrow>
              <mml:mrow>
                <mml:msub>
                  <mml:mi>ρ</mml:mi>
                  <mml:mrow>
                    <mml:mtext>eff</mml:mtext>
                  </mml:mrow>
                </mml:msub>
              </mml:mrow>
            </mml:mfrac>
            <mml:mo>=</mml:mo>
            <mml:mfrac>
              <mml:mrow>
                <mml:msubsup>
                  <mml:mi>v</mml:mi>
                  <mml:mrow>
                    <mml:mtext>eff</mml:mtext>
                  </mml:mrow>
                  <mml:mn>2</mml:mn>
                </mml:msubsup>
                <mml:mo>−</mml:mo>
                <mml:msubsup>
                  <mml:mi>v</mml:mi>
                  <mml:mo>⊕</mml:mo>
                  <mml:mn>2</mml:mn>
                </mml:msubsup>
              </mml:mrow>
              <mml:mrow>
                <mml:msup>
                  <mml:mi>c</mml:mi>
                  <mml:mn>2</mml:mn>
                </mml:msup>
              </mml:mrow>
            </mml:mfrac>
            <mml:mo>−</mml:mo>
            <mml:mfrac>
              <mml:mrow>
                <mml:mn>2</mml:mn>
                <mml:mi>g</mml:mi>
                <mml:mi>h</mml:mi>
              </mml:mrow>
              <mml:mrow>
                <mml:msup>
                  <mml:mi>c</mml:mi>
                  <mml:mn>2</mml:mn>
                </mml:msup>
              </mml:mrow>
            </mml:mfrac>
            <mml:mo>.</mml:mo>
          </mml:mrow>
        </mml:math>
      </disp-formula>
      <p>Combining Equations (65) and (66) gives exactly the same first-order mathematical structure as Equation (64). This equivalence is not a weakness of the model; it is a necessary requirement. Any alternative physical interpretation of the Hafele-Keating experiment must reproduce the experimentally verified relativistic limit.</p>
      <sec id="sec8dot1">
        <title>8.1. Physical Meaning of the Equivalence</title>
        <p>The difference between the two accounts lies in interpretation. In the standard relativistic description, the first term in Equation (64) is gravitational time dilation and the second term is kinematic time dilation. In the present model, the same two terms are interpreted as consequences of changes in the local effective vacuum-density state.</p>
        <p>The altitude term corresponds to radial rarefaction of the structured vacuum. As the aircraft moves to a higher altitude, it occupies a slightly less dense effective vacuum state. According to the wave-speed relation, this rarefaction increases the effective propagation rate of electromagnetic processes and causes the airborne clock to gain time relative to a ground clock.</p>
        <p>The velocity term corresponds to kinetic loading or kinetic rarefaction relative to the rotating Earth-coupled vacuum vortex. Eastward motion increases the effective velocity relative to the vortex and therefore increases the local density loading. This slows the clock. Westward motion partially cancels the rotational velocity of the ground clock and therefore reduces the kinetic density loading. This permits a faster clock rate.</p>
        <p>The same identification may be written in compact form as</p>
        <disp-formula id="FD70">
          <label>(66)</label>
          <mml:math>
            <mml:mrow>
              <mml:mfrac>
                <mml:mrow>
                  <mml:mi>Δ</mml:mi>
                  <mml:mi>ρ</mml:mi>
                </mml:mrow>
                <mml:mrow>
                  <mml:msub>
                    <mml:mi>ρ</mml:mi>
                    <mml:mrow>
                      <mml:mtext>eff</mml:mtext>
                    </mml:mrow>
                  </mml:msub>
                </mml:mrow>
              </mml:mfrac>
              <mml:mo>≃</mml:mo>
              <mml:mfrac>
                <mml:mrow>
                  <mml:msubsup>
                    <mml:mi>v</mml:mi>
                    <mml:mrow>
                      <mml:mtext>air</mml:mtext>
                    </mml:mrow>
                    <mml:mn>2</mml:mn>
                  </mml:msubsup>
                  <mml:mo>−</mml:mo>
                  <mml:msubsup>
                    <mml:mi>v</mml:mi>
                    <mml:mrow>
                      <mml:mtext>ground</mml:mtext>
                    </mml:mrow>
                    <mml:mn>2</mml:mn>
                  </mml:msubsup>
                </mml:mrow>
                <mml:mrow>
                  <mml:msup>
                    <mml:mi>c</mml:mi>
                    <mml:mn>2</mml:mn>
                  </mml:msup>
                </mml:mrow>
              </mml:mfrac>
              <mml:mo>−</mml:mo>
              <mml:mfrac>
                <mml:mrow>
                  <mml:mn>2</mml:mn>
                  <mml:mi>Δ</mml:mi>
                  <mml:mi>Φ</mml:mi>
                </mml:mrow>
                <mml:mrow>
                  <mml:msup>
                    <mml:mi>c</mml:mi>
                    <mml:mn>2</mml:mn>
                  </mml:msup>
                </mml:mrow>
              </mml:mfrac>
              <mml:mo>.</mml:mo>
            </mml:mrow>
          </mml:math>
        </disp-formula>
        <p>This equation shows that the density perturbation is the vacuum-dynamic counterpart of the combined gravitational-potential and kinetic terms of the relativistic calculation.</p>
      </sec>
      <sec id="sec8dot2">
        <title>8.2. Why the Model Is Not Merely a Restatement of Relativity</title>
        <p>A possible objection is that the final time-shift formula resembles the weak-field relativistic expression. This resemblance is expected and should be emphasized rather than hidden. The proposed model is designed to reproduce the known relativistic limit. Its novelty is not a different numerical correction to Hafele-Keating, but a proposed physical interpretation of why that correction arises.</p>
        <p>Relativity describes the effect geometrically: moving clocks follow different spacetime intervals, and clocks at different gravitational potentials experience different proper times. The vacuum-density model describes the same first-order result dynamically: clock rates change because electromagnetic transition processes occur within a local vacuum state whose effective density is slightly modified by altitude and motion through the Earth-coupled vortex.</p>
        <p>The relationship may therefore be summarized conceptually: geometry and vacuum-density dynamics are complementary descriptions of the same weak-field clock-shift structure.</p>
        <p>In compact form, the relationship may be expressed as:</p>
        <p>geometry in relativity ↔ vacuum-density dynamics in the present model.</p>
        <p>This equivalence should be understood in the same spirit as analog-gravity models, where phenomena that are geometrically described in relativity can also emerge from the dynamics of a physical medium. In the present case, the structured vacuum is proposed as the medium-like substrate whose local density state gives rise to the observed clock-rate differences. This interpretive strategy is also consistent with broader emergent-gravity approaches, in which gravitational dynamics are viewed as arising from deeper vacuum, thermodynamic, or statistical degrees of freedom rather than being fundamental at the microscopic level [<xref ref-type="bibr" rid="B16">16</xref>][<xref ref-type="bibr" rid="B17">17</xref>].</p>
      </sec>
      <sec id="sec8dot3">
        <title>8.3. Interpretive Consequences of the Vacuum-Density Framework</title>
        <p>The model should therefore be presented cautiously. It does not claim that standard relativity is wrong. It does not claim that the Hafele-Keating experiment requires a new correction. It claims that the first-order relativistic clock-shift expression can be reinterpreted as an emergent consequence of structured vacuum dynamics.</p>
        <p>This distinction is essential for scientific credibility. If the manuscript is written as an attempt to replace relativity, reviewers will reject it because the standard theory already explains the experiment with high accuracy. If it is written as a complementary physical-substrate interpretation that reproduces the relativistic limit and proposes additional testable consequences, the argument becomes more defensible.</p>
        <p>The key criterion for future development is therefore not only reproduction of Hafele-Keating, but prediction beyond Hafele-Keating. The structured-vacuum interpretation should lead to trajectory-resolved, latitude-dependent, altitude-dependent, and orbital-direction-dependent predictions that can be compared with modern atomic-clock and satellite data. These future tests are discussed in the following section.</p>
      </sec>
    </sec>
    <sec id="sec9">
      <title>9. Effective Vacuum-Density Scale, Limitations, and Future Tests</title>
      <p>The preceding sections showed that the Hafele-Keating clock shifts can be reproduced, to first-order accuracy, from an independently derived vacuum-density perturbation. The density perturbation was obtained from vortex kinematics and altitude-dependent radial rarefaction, not from the observed clock shifts. This removes the main circularity that would otherwise weaken the argument. However, three important questions remain: the physical meaning of the effective density scale, the limitations of the present first-order formulation, and the predictions that could distinguish the model from a merely retrospective reinterpretation.</p>
      <sec id="sec9dot1">
        <title>9.1. Effective Vacuum-Density Scale</title>
        <p>In the present model, the cosmologically inferred dark-energy density is used as the natural order-of-magnitude reference for the vacuum-density scale. However, it is useful to distinguish between the large-scale cosmological density and the local mechanically active density that controls electromagnetic propagation and atomic transition rates. This distinction may be written as</p>
        <disp-formula id="FD71">
          <label>(67)</label>
          <mml:math>
            <mml:mrow>
              <mml:msub>
                <mml:mi>ρ</mml:mi>
                <mml:mrow>
                  <mml:mtext>eff</mml:mtext>
                </mml:mrow>
              </mml:msub>
              <mml:mo>=</mml:mo>
              <mml:mi>κ</mml:mi>
              <mml:msub>
                <mml:mi>ρ</mml:mi>
                <mml:mi>Λ</mml:mi>
              </mml:msub>
              <mml:mo>.</mml:mo>
            </mml:mrow>
          </mml:math>
        </disp-formula>
        <p>where <inline-formula><mml:math><mml:mrow><mml:msub><mml:mi> ρ </mml:mi><mml:mi> Λ </mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> denotes the cosmological vacuum-density scale and <inline-formula><mml:math display="inline"><mml:mi> κ </mml:mi></mml:math></inline-formula> is a dimensionless factor representing the coupling between the large-scale vacuum background and the local mechanically active density. In the simplest version of the model, <italic>κ</italic> is expected to be of order unity.</p>
        <p>If the absolute perturbation is expressed relative to the cosmological scale, then</p>
        <disp-formula id="FD72">
          <label>(68)</label>
          <mml:math>
            <mml:mrow>
              <mml:mi>Δ</mml:mi>
              <mml:mi>ρ</mml:mi>
              <mml:mo>=</mml:mo>
              <mml:msub>
                <mml:mi>ρ</mml:mi>
                <mml:mi>Λ</mml:mi>
              </mml:msub>
              <mml:mrow>
                <mml:mo>[</mml:mo>
                <mml:mrow>
                  <mml:mfrac>
                    <mml:mrow>
                      <mml:msubsup>
                        <mml:mi>v</mml:mi>
                        <mml:mrow>
                          <mml:mtext>eff</mml:mtext>
                        </mml:mrow>
                        <mml:mn>2</mml:mn>
                      </mml:msubsup>
                      <mml:mo>−</mml:mo>
                      <mml:msubsup>
                        <mml:mi>v</mml:mi>
                        <mml:mo>⊕</mml:mo>
                        <mml:mn>2</mml:mn>
                      </mml:msubsup>
                    </mml:mrow>
                    <mml:mrow>
                      <mml:msup>
                        <mml:mi>c</mml:mi>
                        <mml:mn>2</mml:mn>
                      </mml:msup>
                    </mml:mrow>
                  </mml:mfrac>
                  <mml:mo>−</mml:mo>
                  <mml:mfrac>
                    <mml:mrow>
                      <mml:mn>2</mml:mn>
                      <mml:mi>g</mml:mi>
                      <mml:mi>h</mml:mi>
                    </mml:mrow>
                    <mml:mrow>
                      <mml:msup>
                        <mml:mi>c</mml:mi>
                        <mml:mn>2</mml:mn>
                      </mml:msup>
                    </mml:mrow>
                  </mml:mfrac>
                </mml:mrow>
                <mml:mo>]</mml:mo>
              </mml:mrow>
              <mml:mo>.</mml:mo>
            </mml:mrow>
          </mml:math>
        </disp-formula>
        <p>and the effective density required by the measured clock shift can be written as</p>
        <disp-formula id="FD73">
          <label>(69)</label>
          <mml:math>
            <mml:mrow>
              <mml:msub>
                <mml:mi>ρ</mml:mi>
                <mml:mrow>
                  <mml:mtext>eff</mml:mtext>
                </mml:mrow>
              </mml:msub>
              <mml:mo>=</mml:mo>
              <mml:mo>−</mml:mo>
              <mml:mfrac>
                <mml:mrow>
                  <mml:mi>T</mml:mi>
                  <mml:mi>Δ</mml:mi>
                  <mml:mi>ρ</mml:mi>
                </mml:mrow>
                <mml:mrow>
                  <mml:mn>2</mml:mn>
                  <mml:mi>Δ</mml:mi>
                  <mml:mi>τ</mml:mi>
                </mml:mrow>
              </mml:mfrac>
              <mml:mo>.</mml:mo>
            </mml:mrow>
          </mml:math>
        </disp-formula>
        <p>This relation should not be used as the primary derivation of the model. It is only a consistency check after the density perturbation has already been derived from vortex dynamics. Using the first-order trajectory-adjusted parameters applied in Section 6, the values inferred separately from the eastward and westward clock shifts are</p>
        <disp-formula id="FD74">
          <label>(70)</label>
          <mml:math>
            <mml:mrow>
              <mml:msub>
                <mml:mi>ρ</mml:mi>
                <mml:mrow>
                  <mml:mtext>eff</mml:mtext>
                  <mml:mo>,</mml:mo>
                  <mml:mtext>east</mml:mtext>
                </mml:mrow>
              </mml:msub>
              <mml:mo>≃</mml:mo>
              <mml:mn>6.32</mml:mn>
              <mml:mo>×</mml:mo>
              <mml:msup>
                <mml:mrow>
                  <mml:mn>10</mml:mn>
                </mml:mrow>
                <mml:mrow>
                  <mml:mo>−</mml:mo>
                  <mml:mn>27</mml:mn>
                </mml:mrow>
              </mml:msup>
              <mml:mtext>
                 
              </mml:mtext>
              <mml:mtext>kg</mml:mtext>
              <mml:mo>⋅</mml:mo>
              <mml:msup>
                <mml:mtext>m</mml:mtext>
                <mml:mrow>
                  <mml:mo>−</mml:mo>
                  <mml:mn>3</mml:mn>
                </mml:mrow>
              </mml:msup>
              <mml:mo>.</mml:mo>
            </mml:mrow>
          </mml:math>
        </disp-formula>
        <disp-formula id="FD75">
          <label>(71)</label>
          <mml:math>
            <mml:mrow>
              <mml:msub>
                <mml:mi>ρ</mml:mi>
                <mml:mrow>
                  <mml:mtext>eff</mml:mtext>
                  <mml:mo>,</mml:mo>
                  <mml:mtext>west</mml:mtext>
                </mml:mrow>
              </mml:msub>
              <mml:mo>≃</mml:mo>
              <mml:mn>6.30</mml:mn>
              <mml:mo>×</mml:mo>
              <mml:msup>
                <mml:mrow>
                  <mml:mn>10</mml:mn>
                </mml:mrow>
                <mml:mrow>
                  <mml:mo>−</mml:mo>
                  <mml:mn>27</mml:mn>
                </mml:mrow>
              </mml:msup>
              <mml:mtext>
                 
              </mml:mtext>
              <mml:mtext>kg</mml:mtext>
              <mml:mo>⋅</mml:mo>
              <mml:msup>
                <mml:mtext>m</mml:mtext>
                <mml:mrow>
                  <mml:mo>−</mml:mo>
                  <mml:mn>3</mml:mn>
                </mml:mrow>
              </mml:msup>
              <mml:mo>.</mml:mo>
            </mml:mrow>
          </mml:math>
        </disp-formula>
        <p>The mean effective density is therefore</p>
        <disp-formula id="FD76">
          <label>(72)</label>
          <mml:math>
            <mml:mrow>
              <mml:msub>
                <mml:mover accent="true">
                  <mml:mi>ρ</mml:mi>
                  <mml:mo>¯</mml:mo>
                </mml:mover>
                <mml:mrow>
                  <mml:mtext>eff</mml:mtext>
                </mml:mrow>
              </mml:msub>
              <mml:mo>≃</mml:mo>
              <mml:mn>6.31</mml:mn>
              <mml:mo>×</mml:mo>
              <mml:msup>
                <mml:mrow>
                  <mml:mn>10</mml:mn>
                </mml:mrow>
                <mml:mrow>
                  <mml:mo>−</mml:mo>
                  <mml:mn>27</mml:mn>
                </mml:mrow>
              </mml:msup>
              <mml:mtext>
                 
              </mml:mtext>
              <mml:mtext>kg</mml:mtext>
              <mml:mo>⋅</mml:mo>
              <mml:msup>
                <mml:mtext>m</mml:mtext>
                <mml:mrow>
                  <mml:mo>−</mml:mo>
                  <mml:mn>3</mml:mn>
                </mml:mrow>
              </mml:msup>
              <mml:mo>.</mml:mo>
            </mml:mrow>
          </mml:math>
        </disp-formula>
        <p>Compared with the cosmological vacuum-density scale used in the model, this gives</p>
        <disp-formula id="FD77">
          <label>(73)</label>
          <mml:math>
            <mml:mrow>
              <mml:mfrac>
                <mml:mrow>
                  <mml:msub>
                    <mml:mover accent="true">
                      <mml:mi>ρ</mml:mi>
                      <mml:mo>¯</mml:mo>
                    </mml:mover>
                    <mml:mrow>
                      <mml:mtext>eff</mml:mtext>
                    </mml:mrow>
                  </mml:msub>
                </mml:mrow>
                <mml:mrow>
                  <mml:msub>
                    <mml:mi>ρ</mml:mi>
                    <mml:mi>Λ</mml:mi>
                  </mml:msub>
                </mml:mrow>
              </mml:mfrac>
              <mml:mo>≃</mml:mo>
              <mml:mn>1.00.</mml:mn>
            </mml:mrow>
          </mml:math>
        </disp-formula>
        <p>This result does not prove the model, but it supports its internal consistency: the mechanically active density scale required by the first-order clock-shift interpretation is of the same order of magnitude as the cosmologically inferred vacuum-density scale.</p>
      </sec>
      <sec id="sec9dot2">
        <title>9.2. Limitations of the Present Formulation</title>
        <p>Several limitations must be stated clearly. First, the relation between vacuum density and electromagnetic propagation speed has been modeled phenomenologically through a wave-speed analogy. A stronger theory would need to derive the effective vacuum elastic modulus, compressibility, and electromagnetic coupling from a deeper field model.</p>
        <p>Second, the present analysis uses effective mean values for velocity, altitude, latitude, and flight duration. The original Hafele-Keating flights followed complex trajectories, including changing altitude, latitude, ground speed, stopovers, and operational delays. A more rigorous calculation should therefore replace the averaged formula with a trajectory integral:</p>
        <disp-formula id="FD78">
          <label>(74)</label>
          <mml:math>
            <mml:mrow>
              <mml:mi>Δ</mml:mi>
              <mml:mi>τ</mml:mi>
              <mml:mo>=</mml:mo>
              <mml:mstyle displaystyle="true">
                <mml:mrow>
                  <mml:msubsup>
                    <mml:mo>∫</mml:mo>
                    <mml:mn>0</mml:mn>
                    <mml:mi>T</mml:mi>
                  </mml:msubsup>
                  <mml:mrow>
                    <mml:mrow>
                      <mml:mo>[</mml:mo>
                      <mml:mrow>
                        <mml:mfrac>
                          <mml:mrow>
                            <mml:mi>g</mml:mi>
                            <mml:mrow>
                              <mml:mo>(</mml:mo>
                              <mml:mi>t</mml:mi>
                              <mml:mo>)</mml:mo>
                            </mml:mrow>
                            <mml:mi>h</mml:mi>
                            <mml:mrow>
                              <mml:mo>(</mml:mo>
                              <mml:mi>t</mml:mi>
                              <mml:mo>)</mml:mo>
                            </mml:mrow>
                          </mml:mrow>
                          <mml:mrow>
                            <mml:msup>
                              <mml:mi>c</mml:mi>
                              <mml:mn>2</mml:mn>
                            </mml:msup>
                          </mml:mrow>
                        </mml:mfrac>
                        <mml:mo>−</mml:mo>
                        <mml:mfrac>
                          <mml:mrow>
                            <mml:msubsup>
                              <mml:mi>v</mml:mi>
                              <mml:mrow>
                                <mml:mtext>eff</mml:mtext>
                              </mml:mrow>
                              <mml:mn>2</mml:mn>
                            </mml:msubsup>
                            <mml:mrow>
                              <mml:mo>(</mml:mo>
                              <mml:mi>t</mml:mi>
                              <mml:mo>)</mml:mo>
                            </mml:mrow>
                            <mml:mo>−</mml:mo>
                            <mml:msubsup>
                              <mml:mi>v</mml:mi>
                              <mml:mo>⊕</mml:mo>
                              <mml:mn>2</mml:mn>
                            </mml:msubsup>
                            <mml:mrow>
                              <mml:mo>(</mml:mo>
                              <mml:mi>t</mml:mi>
                              <mml:mo>)</mml:mo>
                            </mml:mrow>
                          </mml:mrow>
                          <mml:mrow>
                            <mml:mn>2</mml:mn>
                            <mml:msup>
                              <mml:mi>c</mml:mi>
                              <mml:mn>2</mml:mn>
                            </mml:msup>
                          </mml:mrow>
                        </mml:mfrac>
                      </mml:mrow>
                      <mml:mo>]</mml:mo>
                    </mml:mrow>
                    <mml:mtext>d</mml:mtext>
                    <mml:mi>t</mml:mi>
                  </mml:mrow>
                </mml:mrow>
              </mml:mstyle>
              <mml:mo>.</mml:mo>
            </mml:mrow>
          </mml:math>
        </disp-formula>
        <p>where <inline-formula><mml:math display="inline"><mml:mrow><mml:mi> h </mml:mi><mml:mrow><mml:mo> ( </mml:mo><mml:mi> t </mml:mi><mml:mo> ) </mml:mo></mml:mrow></mml:mrow></mml:math></inline-formula> , <inline-formula><mml:math display="inline"><mml:mrow><mml:msub><mml:mi> v </mml:mi><mml:mrow><mml:mtext> eff </mml:mtext></mml:mrow></mml:msub><mml:mrow><mml:mo> ( </mml:mo><mml:mi> t </mml:mi><mml:mo> ) </mml:mo></mml:mrow></mml:mrow></mml:math></inline-formula> , and <inline-formula><mml:math display="inline"><mml:mrow><mml:msub><mml:mi> v </mml:mi><mml:mo> ⊕ </mml:mo></mml:msub><mml:mrow><mml:mo> ( </mml:mo><mml:mi> t </mml:mi><mml:mo> ) </mml:mo></mml:mrow></mml:mrow></mml:math></inline-formula> vary along the actual route. Such a trajectory-level calculation would be a stronger test of the model than the averaged first-order estimate presented here.</p>
        <p>Third, the model is constructed to reproduce the weak-field relativistic limit. This is a necessary condition, not a weakness. Any alternative physical interpretation of the Hafele-Keating experiment must recover the same first-order structure. The scientific value of the present framework therefore depends on whether it can generate clear physical meaning and additional testable predictions.</p>
        <p>Fourth, the effective density scale <inline-formula><mml:math display="inline"><mml:mrow><mml:msub><mml:mi> ρ </mml:mi><mml:mrow><mml:mtext> eff </mml:mtext></mml:mrow></mml:msub></mml:mrow></mml:math></inline-formula> must ultimately be derived from independent physical principles. In the present article it is constrained by consistency with the cosmological vacuum-density scale and by the Hafele-Keating comparison. Future work should attempt to derive ρ<sub>eff</sub> from vacuum elasticity, quantum-field vacuum structure, or vortex compressibility rather than treating it as an adjustable scale.</p>
      </sec>
      <sec id="sec9dot3">
        <title>9.3. Testable Predictions beyond Hafele-Keating</title>
        <p>The model becomes scientifically stronger only if it predicts effects beyond the original Hafele-Keating data. The most direct prediction concerns latitude. Since the tangential velocity of the Earth-coupled vacuum vortex depends on latitude,</p>
        <disp-formula id="FD79">
          <label>(75)</label>
          <mml:math>
            <mml:mrow>
              <mml:msub>
                <mml:mi>v</mml:mi>
                <mml:mo>⊕</mml:mo>
              </mml:msub>
              <mml:mrow>
                <mml:mo>(</mml:mo>
                <mml:mi>λ</mml:mi>
                <mml:mo>)</mml:mo>
              </mml:mrow>
              <mml:mo>=</mml:mo>
              <mml:msub>
                <mml:mi>Ω</mml:mi>
                <mml:mo>⊕</mml:mo>
              </mml:msub>
              <mml:msub>
                <mml:mi>R</mml:mi>
                <mml:mo>⊕</mml:mo>
              </mml:msub>
              <mml:mi>cos</mml:mi>
              <mml:mi>λ</mml:mi>
              <mml:mo>.</mml:mo>
            </mml:mrow>
          </mml:math>
        </disp-formula>
        <p>the east-west asymmetry should decrease as latitude increases. For two otherwise identical flights at the same altitude, duration, and aircraft speed, the predicted directional difference is</p>
        <disp-formula id="FD80">
          <label>(76)</label>
          <mml:math>
            <mml:mrow>
              <mml:mi>Δ</mml:mi>
              <mml:msub>
                <mml:mi>τ</mml:mi>
                <mml:mrow>
                  <mml:mtext>west</mml:mtext>
                </mml:mrow>
              </mml:msub>
              <mml:mo>−</mml:mo>
              <mml:mi>Δ</mml:mi>
              <mml:msub>
                <mml:mi>τ</mml:mi>
                <mml:mrow>
                  <mml:mtext>east</mml:mtext>
                </mml:mrow>
              </mml:msub>
              <mml:mo>=</mml:mo>
              <mml:mfrac>
                <mml:mrow>
                  <mml:mn>2</mml:mn>
                  <mml:mi>T</mml:mi>
                  <mml:msub>
                    <mml:mi>Ω</mml:mi>
                    <mml:mo>⊕</mml:mo>
                  </mml:msub>
                  <mml:msub>
                    <mml:mi>R</mml:mi>
                    <mml:mo>⊕</mml:mo>
                  </mml:msub>
                  <mml:mi>u</mml:mi>
                  <mml:mi>cos</mml:mi>
                  <mml:mi>λ</mml:mi>
                </mml:mrow>
                <mml:mrow>
                  <mml:msup>
                    <mml:mi>c</mml:mi>
                    <mml:mn>2</mml:mn>
                  </mml:msup>
                </mml:mrow>
              </mml:mfrac>
              <mml:mo>.</mml:mo>
            </mml:mrow>
          </mml:math>
        </disp-formula>
        <p>This equation gives a clear experimental signature: the east-west clock asymmetry should be maximal near the equator and progressively smaller at higher latitudes.</p>
        <p>In the limiting case of a polar route, where <inline-formula><mml:math display="inline"><mml:mrow><mml:mi> cos </mml:mi><mml:mi> λ </mml:mi></mml:mrow></mml:math></inline-formula> approaches zero, the directional asymmetry tends to vanish:</p>
        <disp-formula id="FD81">
          <label>(77)</label>
          <mml:math display="inline">
            <mml:mrow>
              <mml:munder>
                <mml:mrow>
                  <mml:mtext>lim</mml:mtext>
                </mml:mrow>
                <mml:mrow>
                  <mml:mi>λ</mml:mi>
                  <mml:mo>→</mml:mo>
                  <mml:msup>
                    <mml:mrow>
                      <mml:mn>90</mml:mn>
                    </mml:mrow>
                    <mml:mo>∘</mml:mo>
                  </mml:msup>
                </mml:mrow>
              </mml:munder>
              <mml:mrow>
                <mml:mo>(</mml:mo>
                <mml:mrow>
                  <mml:mi>Δ</mml:mi>
                  <mml:msub>
                    <mml:mi>τ</mml:mi>
                    <mml:mrow>
                      <mml:mtext>west</mml:mtext>
                    </mml:mrow>
                  </mml:msub>
                  <mml:mo>−</mml:mo>
                  <mml:mi>Δ</mml:mi>
                  <mml:msub>
                    <mml:mi>τ</mml:mi>
                    <mml:mrow>
                      <mml:mtext>east</mml:mtext>
                    </mml:mrow>
                  </mml:msub>
                </mml:mrow>
                <mml:mo>)</mml:mo>
              </mml:mrow>
              <mml:mo>=</mml:mo>
              <mml:mn>0.</mml:mn>
            </mml:mrow>
          </mml:math>
        </disp-formula>
        <p>The altitude contribution would remain, but the east-west rotational asymmetry would be strongly reduced. This prediction provides a simple way to distinguish the directional vortex component from the common altitude component.</p>
        <p>Additional tests may be designed using modern optical clocks, GPS-tracked aircraft trajectories, counter-rotating satellite orbits, and high-precision interferometric experiments. The most important tests are summarized below.</p>
        <table-wrap id="tbl2">
          <label>Table 2</label>
          <table>
            <tbody>
              <tr>
                <td>
                  <bold>Prediction</bold>
                </td>
                <td>
                  <bold>Expected</bold>
                  <bold>behavior</bold>
                </td>
                <td>
                  <bold>Experimental</bold>
                  <bold>test</bold>
                </td>
              </tr>
              <tr>
                <td>Latitude dependence</td>
                <td>
                  East-west asymmetry scales with cos
                  <italic>λ</italic>
                </td>
                <td>Repeat atomic-clock flights at different latitudes</td>
              </tr>
              <tr>
                <td>Polar-route reduction</td>
                <td>Directional asymmetry tends toward zero near the poles</td>
                <td>Compare polar and equatorial routes</td>
              </tr>
              <tr>
                <td>Altitude sensitivity</td>
                <td>Common positive shift changes by ≈16 ns per km for 40 h</td>
                <td>Repeat flights at different mean cruising altitudes</td>
              </tr>
              <tr>
                <td>Prograde-retrograde satellites</td>
                <td>Opposite orbital directions should show systematic asymmetry</td>
                <td>Compare clock rates in counter-rotating orbital configurations</td>
              </tr>
              <tr>
                <td>Trajectory-level reconstruction</td>
                <td>Segment-by-segment calculation should improve agreement</td>
                <td>Use historical or modern GPS-tracked flight data</td>
              </tr>
            </tbody>
          </table>
        </table-wrap>
      </sec>
    </sec>
    <sec id="sec10">
      <title>10. Conclusions</title>
      <p>The Hafele-Keating experiment remains a landmark demonstration of relativistic clock behavior. Its eastward and westward atomic-clock shifts are accurately described by the standard combination of kinematic and gravitational time dilation. The purpose of the present article has not been to challenge that empirical success, but to explore whether the same weak-field structure can be given a complementary physical interpretation in terms of structured quantum-vacuum dynamics.</p>
      <p>The central result of the present work is that the clock-shift structure can be obtained from an independently derived vacuum-density perturbation. The perturbation is not inferred from the observed time shifts. Instead, it is derived from the effective velocity of the airborne clock relative to an Earth-coupled vacuum vortex, together with the altitude-dependent radial rarefaction of the local vacuum-density state. In this way, the model avoids the circular reasoning that would arise if the Hafele-Keating data were used first to define the density perturbation and then again to reproduce the same data.</p>
      <p>Within this interpretation, eastward and westward flights differ because they interact differently with the rotating structured vacuum. Eastward motion increases the effective velocity relative to the Earth-coupled vortex and therefore increases the kinetic density loading of the local vacuum state. Westward motion partially cancels the rotational component and therefore reduces the effective kinetic loading. The altitude term acts in the same direction for both flights, producing a common contribution associated with reduced vacuum density at higher radial position. The observed signs of the Hafele-Keating shifts therefore arise naturally from the combined directional and radial components of the model.</p>
      <p>The numerical estimates show that, using physically reasonable effective flight parameters, the model reproduces the order, sign, and magnitude of the measured time shifts. The required effective vacuum-density scale lies close to the cosmologically inferred vacuum-density scale, suggesting that the proposed mechanism is at least internally consistent at the level of first-order approximation. This agreement should not be interpreted as definitive proof of the model, but it does show that the framework is quantitatively plausible and worthy of further development.</p>
      <p>A key feature of the proposed interpretation is its relationship to standard relativity. The final clock-shift equation has the same first-order mathematical structure as the weak-field relativistic expression. This similarity is not accidental and should not be regarded as a weakness. Any physically viable reinterpretation must reproduce the established relativistic limit. The difference lies in the meaning assigned to the terms: in the standard account, they are interpreted geometrically as kinematic and gravitational time dilation; in the present account, they are interpreted dynamically as kinetic compression and radial rarefaction of an effective vacuum-density field.</p>
      <p>Several limitations remain. The relation between effective vacuum density and electromagnetic propagation speed is presently phenomenological, although it is motivated by the general wave-speed dependence on medium density. A complete theory would need to derive the effective vacuum elastic modulus, compressibility, and electromagnetic coupling from a deeper vacuum-field model. In addition, the numerical comparison presented here uses averaged flight parameters rather than a full segment-by-segment reconstruction of the original trajectories. Such a trajectory-level calculation should be a priority for future work.</p>
      <p>The value of the framework will ultimately depend on whether it can generate predictions beyond the historical Hafele-Keating experiment. The most direct tests involve latitude dependence, polar-route suppression of east-west asymmetry, altitude-dependent clock shifts, counter-rotating satellite comparisons, and modern optical-clock experiments using precisely tracked trajectories. These tests could determine whether the proposed structured-vacuum interpretation contains physical content beyond a reformulation of the relativistic weak-field limit.</p>
      <p>In conclusion, the quantum-vacuum vortex model offers a possible physical substrate interpretation of relativistic time dilation. It reproduces the Hafele-Keating clock-shift structure from independently derived density perturbations, preserves agreement with the weak-field relativistic formula, and suggests new experimental directions. The model should therefore be regarded not as a replacement for relativity, but as a proposed deeper physical reading of its clock-rate effects in terms of structured vacuum dynamics.</p>
    </sec>
  </body>
  <back>
    <ref-list>
      <title>References</title>
      <ref id="B1">
        <label>1.</label>
        <citation-alternatives>
          <mixed-citation publication-type="other">Einstein, A. (1916) The Foundation of the General Theory of Relativity. <italic>Annalen der Physik</italic>, 49, 769-822. https://doi.org/10.1002/andp.19163540702 <pub-id pub-id-type="doi">10.1002/andp.19163540702</pub-id><ext-link ext-link-type="uri" xlink:href="https://doi.org/10.1002/andp.19163540702">https://doi.org/10.1002/andp.19163540702</ext-link></mixed-citation>
          <element-citation publication-type="other">
            <person-group person-group-type="author">
              <string-name>Einstein, A.</string-name>
            </person-group>
            <year>1916</year>
            <article-title>The Foundation of the General Theory of Relativity</article-title>
            <source>Annalen der Physik</source>
            <volume>49</volume>
            <pub-id pub-id-type="doi">10.1002/andp.19163540702</pub-id>
          </element-citation>
        </citation-alternatives>
      </ref>
      <ref id="B2">
        <label>2.</label>
        <citation-alternatives>
          <mixed-citation publication-type="other">Hafele, J.C. and Keating, R.E. (1972) Around-the-World Atomic Clocks: Predicted Relativistic Time Gains. <italic>Science</italic>, 177, 166-168. <underline> https://doi.org/10.1126/science.177.4044.166 </underline><pub-id pub-id-type="doi">10.1126/science.177.4044.166</pub-id><pub-id pub-id-type="pmid">17779917</pub-id><ext-link ext-link-type="uri" xlink:href="https://doi.org/10.1126/science.177.4044.166">https://doi.org/10.1126/science.177.4044.166</ext-link></mixed-citation>
          <element-citation publication-type="other">
            <person-group person-group-type="author">
              <string-name>Hafele, J.C.</string-name>
              <string-name>Keating, R.E.</string-name>
            </person-group>
            <year>1972</year>
            <article-title>Around-the-World Atomic Clocks: Predicted Relativistic Time Gains</article-title>
            <source>Science</source>
            <volume>177</volume>
            <pub-id pub-id-type="doi">10.1126/science.177.4044.166</pub-id>
            <pub-id pub-id-type="pmid">17779917</pub-id>
          </element-citation>
        </citation-alternatives>
      </ref>
      <ref id="B3">
        <label>3.</label>
        <citation-alternatives>
          <mixed-citation publication-type="confproc">Lorentz, H.A. (1904) Electromagnetic Phenomena in a System Moving with Any Velocity Less than That of Light. <italic>Proceedings of the Royal Netherlands Academy of Arts and Sciences</italic>, 6, 809-831.</mixed-citation>
          <element-citation publication-type="confproc">
            <person-group person-group-type="author">
              <string-name>Lorentz, H.A.</string-name>
            </person-group>
            <year>1904</year>
            <article-title>Electromagnetic Phenomena in a System Moving with Any Velocity Less than That of Light</article-title>
            <source>Proceedings of the Royal Netherlands Academy of Arts and Sciences</source>
            <volume>6</volume>
          </element-citation>
        </citation-alternatives>
      </ref>
      <ref id="B4">
        <label>4.</label>
        <citation-alternatives>
          <mixed-citation publication-type="other">Bohm, D. (1952) A Suggested Interpretation of the Quantum Theory in Terms of “Hidden” Variables. II. <italic>Physical Review</italic>, 85, 180-193. <underline> https://doi.org/10.1103/physrev.85.180 </underline><pub-id pub-id-type="doi">10.1103/physrev.85.180</pub-id><ext-link ext-link-type="uri" xlink:href="https://doi.org/10.1103/physrev.85.180">https://doi.org/10.1103/physrev.85.180</ext-link></mixed-citation>
          <element-citation publication-type="other">
            <person-group person-group-type="author">
              <string-name>Bohm, D.</string-name>
            </person-group>
            <year>1952</year>
            <article-title>A Suggested Interpretation of the Quantum Theory in Terms of “Hidden” Variables</article-title>
            <source>II. Physical Review</source>
            <volume>85</volume>
            <pub-id pub-id-type="doi">10.1103/physrev.85.180</pub-id>
          </element-citation>
        </citation-alternatives>
      </ref>
      <ref id="B5">
        <label>5.</label>
        <citation-alternatives>
          <mixed-citation publication-type="other">Boyer, T.H. (1975) Random Electrodynamics: The Theory of Classical Electrodynamics with Classical Electromagnetic Zero-Point Radiation. <italic>Physical Review D</italic>, 11, 790-808. <underline> https://doi.org/10.1103/physrevd.11.790 </underline><pub-id pub-id-type="doi">10.1103/physrevd.11.790</pub-id><ext-link ext-link-type="uri" xlink:href="https://doi.org/10.1103/physrevd.11.790">https://doi.org/10.1103/physrevd.11.790</ext-link></mixed-citation>
          <element-citation publication-type="other">
            <person-group person-group-type="author">
              <string-name>Boyer, T.H.</string-name>
            </person-group>
            <year>1975</year>
            <article-title>Random Electrodynamics: The Theory of Classical Electrodynamics with Classical Electromagnetic Zero-Point Radiation</article-title>
            <source>Physical Review D</source>
            <volume>11</volume>
            <pub-id pub-id-type="doi">10.1103/physrevd.11.790</pub-id>
          </element-citation>
        </citation-alternatives>
      </ref>
      <ref id="B6">
        <label>6.</label>
        <citation-alternatives>
          <mixed-citation publication-type="other">Haisch, B., Rueda, A. and Puthoff, H.E. (1994) Inertia as a Zero-Point-Field Lorentz Force. <italic>Physical Review A</italic>, 49, 678-694. <underline> https://doi.org/10.1103/physreva.49.678 </underline><pub-id pub-id-type="doi">10.1103/physreva.49.678</pub-id><pub-id pub-id-type="pmid">9910287</pub-id><ext-link ext-link-type="uri" xlink:href="https://doi.org/10.1103/physreva.49.678">https://doi.org/10.1103/physreva.49.678</ext-link></mixed-citation>
          <element-citation publication-type="other">
            <person-group person-group-type="author">
              <string-name>Haisch, B.</string-name>
              <string-name>Rueda, A.</string-name>
              <string-name>Puthoff, H.E.</string-name>
            </person-group>
            <year>1994</year>
            <article-title>Inertia as a Zero-Point-Field Lorentz Force</article-title>
            <source>Physical Review A</source>
            <volume>49</volume>
            <pub-id pub-id-type="doi">10.1103/physreva.49.678</pub-id>
            <pub-id pub-id-type="pmid">9910287</pub-id>
          </element-citation>
        </citation-alternatives>
      </ref>
      <ref id="B7">
        <label>7.</label>
        <citation-alternatives>
          <mixed-citation publication-type="other">Rueda, A. and Haisch, B. (1998) Gravity and the Quantum Vacuum Inertia Hypothesis. <italic>Annalen der Physik</italic>, 10, 393-414.</mixed-citation>
          <element-citation publication-type="other">
            <person-group person-group-type="author">
              <string-name>Rueda, A.</string-name>
              <string-name>Haisch, B.</string-name>
            </person-group>
            <year>1998</year>
            <article-title>Gravity and the Quantum Vacuum Inertia Hypothesis</article-title>
            <source>Annalen der Physik</source>
            <volume>10</volume>
          </element-citation>
        </citation-alternatives>
      </ref>
      <ref id="B8">
        <label>8.</label>
        <citation-alternatives>
          <mixed-citation publication-type="other">Barceló, C., Liberati, S. and Visser, M. (2011) Analogue Gravity. <italic>Living Reviews in Relativity</italic>, 14, 3. <underline> https://doi.org/10.12942/lrr-2011-3 </underline><pub-id pub-id-type="doi">10.12942/lrr-2011-3</pub-id><pub-id pub-id-type="pmid">28179830</pub-id><ext-link ext-link-type="uri" xlink:href="https://doi.org/10.12942/lrr-2011-3">https://doi.org/10.12942/lrr-2011-3</ext-link></mixed-citation>
          <element-citation publication-type="other">
            <person-group person-group-type="author">
              <string-name>Liberati, S.</string-name>
              <string-name>Visser, M.</string-name>
            </person-group>
            <year>2011</year>
            <article-title>Analogue Gravity</article-title>
            <source>Living Reviews in Relativity</source>
            <volume>14</volume>
            <pub-id pub-id-type="doi">10.12942/lrr-2011-3</pub-id>
            <pub-id pub-id-type="pmid">28179830</pub-id>
          </element-citation>
        </citation-alternatives>
      </ref>
      <ref id="B9">
        <label>9.</label>
        <citation-alternatives>
          <mixed-citation publication-type="book">Volovik, G.E. (2003) The Universe in a Helium Droplet. Oxford University Press.</mixed-citation>
          <element-citation publication-type="book">
            <person-group person-group-type="author">
              <string-name>Volovik, G.E.</string-name>
            </person-group>
            <year>2003</year>
            <article-title>The Universe in a Helium Droplet</article-title>
          </element-citation>
        </citation-alternatives>
      </ref>
      <ref id="B10">
        <label>10.</label>
        <citation-alternatives>
          <mixed-citation publication-type="book">Birrell, N.D. and Davies, P.C.W. (1982) Quantum Fields in Curved Space. Cambridge University Press. <underline> https://doi.org/10.1017/cbo9780511622632 </underline><pub-id pub-id-type="doi">10.1017/cbo9780511622632</pub-id><ext-link ext-link-type="uri" xlink:href="https://doi.org/10.1017/cbo9780511622632">https://doi.org/10.1017/cbo9780511622632</ext-link></mixed-citation>
          <element-citation publication-type="book">
            <person-group person-group-type="author">
              <string-name>Birrell, N.D.</string-name>
              <string-name>Davies, P.C.W.</string-name>
            </person-group>
            <year>1982</year>
            <article-title>Quantum Fields in Curved Space</article-title>
            <pub-id pub-id-type="doi">10.1017/cbo9780511622632</pub-id>
          </element-citation>
        </citation-alternatives>
      </ref>
      <ref id="B11">
        <label>11.</label>
        <citation-alternatives>
          <mixed-citation publication-type="book">Bureau International des Poids et Mesures (BIPM) (2019) The International System of Units (SI). 9th Edition, Bureau International des Poids et Mesures.</mixed-citation>
          <element-citation publication-type="book">
            <person-group person-group-type="author">
              <string-name>Edition, B</string-name>
            </person-group>
            <year>2019</year>
            <article-title>The International System of Units (SI)</article-title>
            <source>9th Edition</source>
          </element-citation>
        </citation-alternatives>
      </ref>
      <ref id="B12">
        <label>12.</label>
        <citation-alternatives>
          <mixed-citation publication-type="book">Ramsey, N.F. (1956) Molecular Beams. Oxford University Press.</mixed-citation>
          <element-citation publication-type="book">
            <person-group person-group-type="author">
              <string-name>Ramsey, N.F.</string-name>
            </person-group>
            <year>1956</year>
            <article-title>Molecular Beams</article-title>
          </element-citation>
        </citation-alternatives>
      </ref>
      <ref id="B13">
        <label>13.</label>
        <citation-alternatives>
          <mixed-citation publication-type="other">Essen, L. and Parry, J.V.L. (1955) An Atomic Standard of Frequency and Time Interval: A Cæsium Resonator. <italic>Nature</italic>, 176, 280-282. <underline> https://doi.org/10.1038/176280a0 </underline><pub-id pub-id-type="doi">10.1038/176280a0</pub-id><ext-link ext-link-type="uri" xlink:href="https://doi.org/10.1038/176280a0">https://doi.org/10.1038/176280a0</ext-link></mixed-citation>
          <element-citation publication-type="other">
            <person-group person-group-type="author">
              <string-name>Essen, L.</string-name>
              <string-name>Parry, J.V.L.</string-name>
            </person-group>
            <year>1955</year>
            <article-title>An Atomic Standard of Frequency and Time Interval: A Cæsium Resonator</article-title>
            <source>Nature</source>
            <volume>176</volume>
            <pub-id pub-id-type="doi">10.1038/176280a0</pub-id>
          </element-citation>
        </citation-alternatives>
      </ref>
      <ref id="B14">
        <label>14.</label>
        <citation-alternatives>
          <mixed-citation publication-type="other">Heavner, T.P., Jefferts, S.R., Donley, E.A., Shirley, J.H. and Parker, T.E. (2005) NIST-F1: Recent Improvements and Accuracy Evaluations. <italic>Metrologia</italic>, 42, 411-422. <underline> https://doi.org/10.1088/0026-1394/42/5/012 </underline><pub-id pub-id-type="doi">10.1088/0026-1394/42/5/012</pub-id><ext-link ext-link-type="uri" xlink:href="https://doi.org/10.1088/0026-1394/42/5/012">https://doi.org/10.1088/0026-1394/42/5/012</ext-link></mixed-citation>
          <element-citation publication-type="other">
            <person-group person-group-type="author">
              <string-name>Heavner, T.P.</string-name>
              <string-name>Jefferts, S.R.</string-name>
              <string-name>Donley, E.A.</string-name>
              <string-name>Shirley, J.H.</string-name>
              <string-name>Parker, T.E.</string-name>
            </person-group>
            <year>2005</year>
            <article-title>NIST-F1: Recent Improvements and Accuracy Evaluations</article-title>
            <source>Metrologia</source>
            <volume>42</volume>
            <pub-id pub-id-type="doi">10.1088/0026-1394/42/5/012</pub-id>
          </element-citation>
        </citation-alternatives>
      </ref>
      <ref id="B15">
        <label>15.</label>
        <citation-alternatives>
          <mixed-citation publication-type="other">Ludlow, A.D., Boyd, M.M., Ye, J., Peik, E. and Schmidt, P.O. (2015) Optical Atomic Clocks. <italic>Reviews of Modern Physics</italic>, 87, 637-701. <underline> https://doi.org/10.1103/revmodphys.87.637 </underline><pub-id pub-id-type="doi">10.1103/revmodphys.87.637</pub-id><ext-link ext-link-type="uri" xlink:href="https://doi.org/10.1103/revmodphys.87.637">https://doi.org/10.1103/revmodphys.87.637</ext-link></mixed-citation>
          <element-citation publication-type="other">
            <person-group person-group-type="author">
              <string-name>Ludlow, A.D.</string-name>
              <string-name>Boyd, M.M.</string-name>
              <string-name>Ye, J.</string-name>
              <string-name>Peik, E.</string-name>
              <string-name>Schmidt, P.O.</string-name>
            </person-group>
            <year>2015</year>
            <article-title>Optical Atomic Clocks</article-title>
            <source>Reviews of Modern Physics</source>
            <volume>87</volume>
            <pub-id pub-id-type="doi">10.1103/revmodphys.87.637</pub-id>
          </element-citation>
        </citation-alternatives>
      </ref>
      <ref id="B16">
        <label>16.</label>
        <citation-alternatives>
          <mixed-citation publication-type="other">Sakharov, A.D. (1968) Vacuum Quantum Fluctuations in Curved Space and the Theory of Gravitation. <italic>Soviet Physics Doklady</italic>, 12, 1040-1041.</mixed-citation>
          <element-citation publication-type="other">
            <person-group person-group-type="author">
              <string-name>Sakharov, A.D.</string-name>
            </person-group>
            <year>1968</year>
            <article-title>Vacuum Quantum Fluctuations in Curved Space and the Theory of Gravitation</article-title>
            <source>Soviet Physics Doklady</source>
            <volume>12</volume>
          </element-citation>
        </citation-alternatives>
      </ref>
      <ref id="B17">
        <label>17.</label>
        <citation-alternatives>
          <mixed-citation publication-type="other">Jacobson, T. (1995) Thermodynamics of Spacetime: The Einstein Equation of State. <italic>Physical Review Letters</italic>, 75, 1260-1263. <underline> https://doi.org/10.1103/physrevlett.75.1260 </underline><pub-id pub-id-type="doi">10.1103/physrevlett.75.1260</pub-id><pub-id pub-id-type="pmid">10060248</pub-id><ext-link ext-link-type="uri" xlink:href="https://doi.org/10.1103/physrevlett.75.1260">https://doi.org/10.1103/physrevlett.75.1260</ext-link></mixed-citation>
          <element-citation publication-type="other">
            <person-group person-group-type="author">
              <string-name>Jacobson, T.</string-name>
            </person-group>
            <year>1995</year>
            <article-title>Thermodynamics of Spacetime: The Einstein Equation of State</article-title>
            <source>Physical Review Letters</source>
            <volume>75</volume>
            <pub-id pub-id-type="doi">10.1103/physrevlett.75.1260</pub-id>
            <pub-id pub-id-type="pmid">10060248</pub-id>
          </element-citation>
        </citation-alternatives>
      </ref>
    </ref-list>
  </back>
</article>