Polarization of vacuum fluctuations: source of the vacuum permittivity and speed of light

There are two types of fluctuations in the quantum vacuum: type 1 vacuum fluctuations are on shell and can interact with matter in specific, limited ways that have observable consequences; type 2 vacuum fluctuations are off shell and cannot interact with matter. A photon will polarize a type 1, bound, charged lepton-antilepton vacuum fluctuation in much the same manner that it would polarize a dielectric, suggesting the method used here for calculating the permittivity $\epsilon_0$ of the vacuum. In a model that retains only leading terms, $\epsilon_0 \cong (6\mu_0/\pi)(8e^2/\hbar)^2= 9.10\times 10^{-12}$ C/(Vm). The calculated value for $\epsilon_0$ is 2.7\% more than the accepted value. The permittivity of the vacuum, in turn, determines the speed $c$ of light in the vacuum. Since the vacuum is at rest with respect to every inertial frame of reference, $c$ is the same in every inertial reference frame.

dielectric constant of the vacuum: "Because of the polarizability of the nascent pairs, the dielectric constant of space into which no matter has been introduced differs from that of truly empty space." Also in 1934 the idea of treating the vacuum as a medium with electric and magnetic polarizability was discussed by Pauli and Weisskopf [2] and again two years later by Weisskopf [3]. In 1957 Dicke [4] wrote about the possibility that the vacuum could be considered as a dielectric medium.
Wilczek's 2008 book The Lightness of Being [5] expressed the fundamental characteristics of space and time as properties of the Grid, "the entity we perceive as empty space. Our deepest physical theories reveal it to be highly structured; indeed, it appears as the primary ingredient of reality." Recently the possibility that the properties of the vacuum determine, in the vacuum, the speed of light and the permittivity have been explored by a number of authors [6,7,8,9].
A basic postulate of physics is that the properties of a physical system are determined by the structure of that system. The only events that occur in the quantum vacuum are the spontaneous appearance and subsequent disappearance of vacuum fluctuations that occur for every type of particle-antiparticle pair, including pairs for which the particle is its own antiparticle, as is the case for photons.
The photon has spin 1; a vacuum fluctuation of the electromagnetic field must have spin 0. 2 It is important to point out that, other than the preceding comment, this article does not discuss fluctuations of the QED vacuum, on which there already exists considerable literature [11,12,13] and references therein. Accordingly, here it will not be possible to address questions as to whether the cosmological constant is affected by or results from QED fluctuations. In this regard, Leonhardt [14] has shown that the electromagnetic contribution to the cosmological constant need not include fluctuations of the QED quantum vacuum. It will also not be possible to address the Casimir force, long thought to be the result of zero-point energy of the QED vacuum. However, the energy of the QED vacuum is no longer considered to be responsible for the Casimir force [15,16,17]. A recent article by Meis [18] provides insights into understanding the electromagnetic quantum vacuum and the possibility of eliminating infinities associated with it.
Only fluctuations of the quantum fermion vacuum will be considered here, and a (real) test photon interacting with those fluctuations will be shown to determine the values of ǫ 0 and c. To minimize the violation of conservation of energy allowed by the Heisenberg uncertainty principle and to avoid violating conservation of angular momentum, in any inertial frame a vacuum fluctuation consisting of a charged fermion-antifermion pair must appear with its center of mass at rest and in the least energetic bound state that has zero angular momentum. In addition to the above theoretical reason, there is a second, physical reason for expecting charged fermion-antifermion vacuum fluctuations to appear as bound states: the value of the permittivity ǫ of a physical dielectric is determined by the electromagnetic properties of the physical dielectric and results from atoms or molecules in the dielectric being polarized by an electric field (or photons). Just as an effective spring constant is associated with the resultant relative motion of the constituents of the atoms or molecules, an effective spring constant exists for a bound, charged fermion-antifermion pair. Transient fermion-antifermion atoms -parapositronium for an electron-positron pair -are of primary interest in this article, and, as will be demonstrated, these transient atoms can be observed indirectly by the effect of their interaction with photons moving through the vacuum.
If the value of the permittivity of the vacuum ǫ 0 is determined by the structure of the vacuum, it should be possible to calculate ǫ 0 from that structure by examining the (polarizing) interaction of photons "introduced" into the vacuum as test particles. The model of the vacuum discussed here yields an approximate formula for the permittivity of the vacuum 3 ǫ 0 ∼ = 6µ 0 π 8e 2 2 = 9.10 × 10 −12 C Vm . (1) In the above equation e is the magnitude of the (renormalized) charge of an electron, is Planck's constant divided by 2π, and µ 0 4 is the permeability of the vacuum. The calculated value for ǫ 0 is 2.7% more than the accepted value. The calculation of ǫ 0 has been simplified, and the numerical accuracy has been reduced, by including (a) only contributions to lowest order in what turns out to be an expansion in powers of the fine-structure constant α and (b) only the interactions of photons with bound states of charged lepton-antilepton vacuum fluctuations.
Once ǫ 0 has been calculated, a formula for the speed c of light in the vacuum is obtained from c = 1/ √ µ 0 ǫ 0 and the formula for ǫ 0 in (1), The value of c calculated from this model is 1.4% less than the defined value c = 3.00 × 10 8 m/s. Thus the calculation in the following sections shows that a photon in the vacuum is slowed by interacting with and exerting a polarizing effect on a bound, charged lepton-antilepton vacuum fluctuation. An observer in any inertial frame is unable to detect the relative motion of the inertial frame and the quantum vacuum. Therefore the observer would conclude that the inertial frame and the vacuum are at rest with respect to each other. As Leonhardt et al. [19] state,"In free space the vacuum is Lorentz invariant, so a 3 SI units are used throughout this article except in the decay rate calculation in the "Appendix" where, following convention, and c are replaced by 1. After reinstating factors of and c, from (A.26) the electromagnetic decay rate for a photon-excited, parapositronium VF is Γ p−Ps = (α 5 mec 2 )/ . The mean lifetime of the state is 1/Γ p−Ps ∼ = 6.2 × 10 −11 s. 4 The value of µ 0 was originally arbitrarily chosen so that the rationalized meter-kilogramsecond unit of current was equal in size to the ampere in the "electromagnetic (emu)" system. As a consequence, µ 0 , which is a measurement-system constant, not a property of the vacuum, was defined to be exactly µ 0 ≡ 4π × 10 −7 H/m. Since physicists arbitrarily chose the value of µ 0 , µ 0 could neither be measured experimentally nor calculated theoretically. As of May 20, 2019, µ 0 is no longer the defined constant 4π × 10 −7 H/m. Instead µ 0 is defined by the equation µ 0 ≡ 2αh/(ce 2 ), which follows from the definition of the fine-structure constant and c = 1/ √ µ 0 ǫ 0 . Similarly, ǫ 0 ≡ e 2 /(2αhc), which follows from the definition of the fine-structure constant. These new definitions for µ 0 and ǫ 0 satisfy the condition µ 0 ǫ 0 = 1/c 2 , as they must. There is no effect so the value of µ 0 is unaffected by fluctuations of the charged fermion-antifermion vacuum.
uniformly moving observer would not see any effect due to motion, but an accelerated observer would. The latter is known as the Unruh effect [20]." The interaction of photons with fermion vacuum fluctuations determines the speed of light in the vacuum, and in the vacuum fermion vacuum fluctuations are at rest with respect to any inertial frame. It then follows that c is the same in any inertial frame, which is one of two postulates in the 1905 paper "On the Electrodynamics of Moving Bodies" [21] in which special relativity was introduced by Einstein. The derivation here of the value of c provides a theoretical explanation for why Einstein's postulate is true, obviating the need for the postulate as a separate concept. There is a second theoretical approach for calculating c that would satisfy Einstein's postulate that c is the same in any inertial frame: if c could be calculated from an internal consistency condition within QED, then c would be the same in every inertial frame since the equations of physics are the same in every inertial frame. Such a calculation has never been completed.

Vacuum fluctuations
Vacuum fluctuations in the quantum fermion vacuum refer to two distinct entities, labeled here as type 1 and type 2, with very different properties. The fieldtheoretical proofs that vacuum fluctuations exist allow the two types to be clearly distinguished. Based on those proofs, a type 1 vacuum fluctuation is defined to be a vacuum fluctuation (net charge zero) that appears as an external particle in a Feynman diagram, is on shell, and exists only for a time allowed by the Heisenberg uncertainty principle [22,23]. A type 1 vacuum fluctuation can consist of a charged or uncharged particle-antiparticle pair.
The appearance of a type 1 vacuum fluctuation is a stochastic process: as such, either a type 1 vacuum fluctuation does not appear at all or it appears on shell. Once charged lepton-antilepton type 1 vacuum fluctuations spontaneously appear, as will be discussed in Sec. 3 they can interact with physical particles, but only in very limited ways. And if they do not interact, the charged pair will simply annihilate, returning the borrowed energy to the vacuum. The calculation here relies on the "introduction" of a (real) photon into the quantum vacuum as a test particle. The interaction of that photon with a bound, charged lepton-antilepton type 1 vacuum fluctuation creates a polarized type 1 vacuum fluctuation that subsequently decays, emitting a real photon. The polarized, type 1 vacuum fluctuation is a quasi-stationary state, a state for which the lifetime of the state is large in comparison with the uncertainty in the lifetime [24]. In addition to interacting with a photon to form a quasi-stationary state, a type 1 vacuum fluctuation can also interact through the annihilation of a particle and antiparticle. For example, a real electron can annihilate with the positron of a type 1 vacuum fluctuation. The electron that was part of the type 1 vacuum fluctuation then becomes a real electron with a location different from that of the original, physical electron, giving rise to zitterbewegung ("trembling motion") [23].
Texts on field theory define type 2 vacuum fluctuations to be vacuum diagrams [25], which are the same as vacuum bubbles [22], a class of perturbation effects in quantum field theories with Feynman diagrams that look like bubbles that originate from and terminate in the vacuum. Vacuum diagrams do not make a contribution to physical processes [22,26,25]; consequently, type 2 vacuum fluctuations do not have observable effects of the type described here.
Vacuum fluctuations are sometimes confused with virtual particles, for which a widely accepted definition exists. For example Penrose [27] defines a virtual particle as follows: A virtual particle is "off-shell" and "can occur only in the interior of a Feynman graph." A virtual particle-antiparticle pair annihilates within a time that is "short enough that the energy required to produce the pair falls within the uncertainties of Heisenberg's time-energy relation. . . ." Both the definition of a virtual particle by Langacker [28] and the use of the term "virtual particle" by Berestetskii et al. [29] and by Gottfried and Weisskopf [30] are consistent with Penrose's definition.
Virtual particles, real particles, type 1 vacuum fluctuations and type 2 vacuum fluctuations all have different properties: virtual particles are off shell and transitory, real particles are on shell and their energy is permanently present, type 1 vacuum fluctuations are on shell and transitory and type 2 vacuum fluctuations are a collection of interacting, virtual particles that arise in and then disappear back into the vacuum so they are off shell and transitory. Perhaps because virtual particles and both types of vacuum fluctuations are transitory, the two terms "virtual particles" and "vacuum fluctuations" are sometimes used as synonyms. They are not. In this article the acronym VF refers exclusively to a type 1 vacuum fluctuation that consists of a charged particle bound with angular momentum zero to its antiparticle.
Field theory provides an explanation for the source of the energy available for the creation of a VF and a proof that VFs must exist. The structure of an atom that is a VF does not play a significant role in the discussion of the production of that VF. (Of course, the structure of the VF, which results from the electromagnetic interaction of the charged lepton-antilepton pair, is important for the calculation of the decay rate of the atom that is a VF.) Thus an atom consisting of a charged lepton and corresponding antilepton in its ground state with zero angular momentum can, as far as its creation is concerned, be approximately represented by a free, neutral, spin-0, Klein-Gordon field φ(x), as first suggested by Pauli and Weisskopf [2] and elaborated on by Wentzel [31]. Using a field to describe a particle with internal degrees of freedom is discussed in Ref. [23]; representing a particle by a quantum field is essential to understanding the mathematical structure of a VF of a quantum field. Indeed, it is the quanta of a free field that behave as free particles.
The vacuum expectation value of the free Hamiltonian H for a neutral, spin-0 field representing an atom with mass M and momentum p = k can be written in the form, For a lepton-antilepton VF the mass M in the above equation is twice the mass of the lepton minus the binding energy of the lepton-antilepton bound state. The energy in the vacuum, called the zero-point energy, appears on the right-hand side of (3). The sum over k is infinite because there are an infinite number of cells in k-space, implying that the energy in the vacuum is infinite. To ensure Lorentz invariance and translation invariance of the vacuum, VFs must be created at rest in any inertial frame, implying that k = 0. It then follows that the only zero-point energy available to create a VF at rest is the single, finite term ω k=0 = M c 2 . The factor of 1/2 in (3) occurs because VFs are, on average, present only for half of the time.
To show that VFs of charged lepton-antilepton pairs must exist, note that the expectation value of the product of the free field φ(x) at two different locations x and x ′ is [23,22] Eq. (4) has the feature that (0|φ 2 (x)|0) is infinite. However any vacuum fluctuation must have a finite size as will become apparent in the next section. As a consequence, the symbol (0|φ(x)φ(x ′ )|0) has no physical meaning unless x and x ′ are averaged over the size associated with the vacuum fluctuation. This point is particularly emphasized by [22]. The fact that the vacuum expectation value of the product φ(x)φ(x ′ ) is nonzero demonstrates that the free field φ(x) in the vacuum is nonzero. Because VFs occur for noninteracting fields, the quanta associated with the fields are on shell: note from (4) that the integral over d 3 k is subject to the constraint k 0 = ω k . It then immediately follows that VFs only affect properties of the vacuum and do not alter existing predictions of QED [32]. If a photon were emitted by some process, it could interact with a VF. But when the VF annihilated, a photon identical to the incident photon would be emitted. The photon emitted by the original process would just be traveling through the vacuum as usual. Also, a photon emitted when a VF annihilated could interact with matter. This would just be the process of an external photon being incident on and interacting with matter. Thus, the existence of fluctuations of charged fermion fields does not affect QED calculations. However, if a VF interacted in a manner that supplied energy equal to or greater than the energy originally borrowed from the vacuum for its creation, the VF could become a physical particle-antiparticle pair, and the pair would interact as would any other physical particles.

Calculation of dielectric properties of the vacuum: the characteristics of a VF
An electron-positron VF appears in the vacuum at rest as parapositronium [33], a singlet spin state with the lowest positronium bound-state energy and zero angular momentum. Here attention is restricted to parapositronium since the corresponding results for muon-antimuon and tau-antitau VFs immediately follow by replacing the electron mass with the muon or tau mass, respectively. The non-relativistic binding energy for parapositronium, E p−P s , is obtained from the n = 1 binding energy of hydrogen [34] by changing the reduced mass: µ hydrogen ≃ m e → µ p−P s = m e /2, where m e is the mass of an electron: In (5) α = e 2 /(4πǫ 0 c) is the fine-structure constant. From both the classical [35,36,37] and quantum [38,39] calculations of the permittivity of physical matter consisting of atoms (or molecules), it follows that for matter to possess permittivity, the atoms must be able to oscillate when interacting with an electric field (or photons). Taking the x-axis to point in the direction of the electric field, since the atom oscillates along the direction of the electric field, only the oscillatory properties of the atom along the x-axis are of significance in the interaction. Thus to first order in α, when calculating permittivity, the interaction of the electromagnetic field with the atom can be described by a one-dimensional oscillator. Although parapositronium VFs do not exhibit resonant behavior, the potential that binds the electron and positron into the parapositronium VF has a relative minimum and, for small oscillations, acts as if it were a harmonic oscillator with a resonant frequency.
As pointed out by Feynman [40], when interacting with an electric field, an atom in its ground state interacts with the electric field as if it were a harmonic oscillator with the first two energy levels separated by the binding energy of the atom. Since adjacent harmonic oscillator energy levels are separated by an energy ω 0 , from (5), The effective spring constant of parapositronium is K p−P s = µ p−P s (ω 0 ) 2 so the Hamiltonian H 0 describing an oscillating parapositronium VF is Using the formula for the energy [34] of a harmonic oscillator in one dimension, the energy eigenvalues E 0 n of the above Hamiltonian are For parapositronium that is a VF, the Heisenberg uncertainty principle is Denoting the mass of an electron (or positron) by m e , ∆E p−Ps is the energy 2m e c 2 for the production of parapositronium that is a VF 5 . The minimum time ∆t is the average lifetime ∆t p−Ps for the existence of parapositronium that is a VF. Thus (9) yields During the time ∆t p−Ps , a beam of light travels a distance L p−Ps given by Since a parapositronium VF appears from the vacuum at essentially a single location and since nothing can travel faster than the speed of light, while they exist the maximum distance between the electron and positron in parapositronium is L p−Ps . Having already borrowed energy 2m e c 2 from a volume V p−Ps ≡ (L p−Ps ) 3 of the vacuum, a second parapositronium VF is unlikely to form in the same volume while the first still exists. This suggests the ansatz, a result that can immediately be generalized to other charged lepton-antilepton VFs and quark-antiquark VFs. The center of mass of a VF will remain fixed, accompanied by an internal zitterbewegung of the VF, giving the VF its size [41,42]. The zitterbewegung internal to the VF has an amplitude [43] where the final equality follows from (11). Eq. (13) allows a length L Z p−Ps to be associated with the zitterbewegung. The volume (L Z p−Ps ) 3 = (L p−Ps ) 3 ≡ V p−Ps results from zitterbewegung and represents the volume of a VF. Requiring that there be only one VF in the volume of a VF as calculated from zitterbewegung yields the same formula previously obtained for the number density of VFs. As a result of zitterbewegung, the VF has a size V p−Ps over which it must be averaged, thus removing the infinity in (0|φ 2 (x)|0) discussed following (4).
The number of lepton-antilepton VFs per unit volume has values from 1.12 × 10 39 /m 3 for electron-positron VFs to 4.70 × 10 49 /m 3 for tau-antitau VFs. In contrast the number density of atoms or molecules of an ideal gas at STP is 2.68 × 10 25 /m 3 . It is possible for the number density of VFs to be many orders of magnitude greater than the number density of atoms or molecules of an ideal gas because a VF cannot exert a force.
First consider the electromagnetic force: if a VF has not already absorbed radiation, it cannot spontaneously emit radiation. If it did, the radiated photons would exist after the VF has disappeared back into the vacuum, permanently violating conservation of energy. If a VF has interacted with a photon, when the VF vanishes back into the vacuum, the VF must emit a photon identical to the incident photon in order to conserve energy, momentum, and angular momentum. Since a VF cannot "permanently" exchange a photon with either another VF or a physical quantum, it cannot exert a force on either. Similar arguments verify that VFs cannot exert a force of any type. There is also a simple, classical argument that a VF cannot exert a force on a physical particle: if it could, the energy associated with the work done by the force would remain after the VF vanished back into the vacuum, permanently violating conservation of energy.
The dipole moment (operator) p x of an atom in the presence an electromagnetic wave with its electric field along the x-axis is p x = e(x + − x − ) ≡ ex, where x − and x + are, respectively, the coordinates of the electron and positron. The Hamiltonian H 1 describing the interaction of the electric dipole of the atom with the electromagnetic wave E x = E 0 cos ωt is When examining (14), the following question immediately arises: when a VF interacts with a photon, does it typically interact with one or multiple photons? This is the essential question of photon-VF interactions; the answer is that if a lepton-antilepton VF interacts with a photon at all, it essentially always interacts with only one photon. As a result the electric field of each individual photon polarizes each individual VF with which it interacts, as suggested by the title of the paper.
To see this consider the photon number density N γ of a laser beam with crosssectional area A and wavelength λ. During a time ∆t the laser emits energy ∆E, which equals the energy per unit volume N γ hc/λ of the beam multiplied by the volume Ac∆t of the pulse, Solving (15) for N γ and noting that the power P = ∆E/∆t, In a 6,000 W, CO 2 cutting laser with a wavelength of 10µm and a beam radius of 0.16 mm, the number density N γ of photons is on the order of 10 22 photons/m 3 . Even in such an intense laser beam, the number density of lepton-antilepton VFs is is much greater than the number density of photons.

Calculation of dielectric properties of the vacuum: dipole moment of VFs
When a photon interacts with a parapositroniun VF, for example, and combines its energy, momentum and angular momentum with that of the parapositroniun VF, the electric field of the photon that interacts with the parapositronium VF is the electric field at the instant t i that the interaction takes place: E(t i ) = E 0 cos ωt i ≡ E 0 . Thus the electric field E 0x describes the interaction of a single photon with a VF and (14) must be replaced by The polarization of a bound, lepton-antilepton VF (a parapositronium vacuum fluctuation, for example) can be obtained from a quantum mechanics calculation that is similar to the calculation of the polarization of an ordinary dielectric. There are, however, two major differences: (1) At the instant a photon interacts with a vacuum fluctuation, the value of the electric field associated with the photon is the value at the moment of interaction. (2) Since an atom or molecule in a dielectric usually has an axis of orientation, a factor of one third is required when averaging over the three possible axes of orientation of the atom or molecule [45]. Because vacuum fluctuations are spherically symmetric and thus do not have an orientation axis, the factor of one third is not present.
The polarization that results when a parapositronium atom in its ground state interacts with a photon is readily calculated. Let ψ 0 n (x) be solutions to the timeindependent, one-dimensional, unperturbed Schrödinger equation, where H 0 and E 0 n are given by (7) and (8), respectively. From stationary perturbation theory, to lowest order in the perturbation, the ground-state solution ψ 1 0 (x) of the Schrödinger equation corresponding to the Hamiltonian H 0 + H 1 VF is [34] The only nonzero contribution to the sum occurs for n ′ = 1. Then in the above equation the integral equals −eE 0 /(2µω 0 ), and the energy difference The expectation value p VF of the electric dipole moment in the state characterized by ψ 1 0 (x) is VFs are not restricted to electron-positron pairs. To allow for more than one type of charged particle-antiparticle VF that binds into an atom, an index j is added to the charge, e → q j , the reduced mass, µ → µ j , and resonant frequency, ω 0 → ω 0 j . Thus (21) is replaced by Eq. (22) is the expectation value of an electric dipole that occurs in the vacuum as a result of a VF being present, interacting with a photon, and becoming polarized. The electric field of a photon that interacts with a parapositronium VF is the electric field of the photon at the instant t i that the photon becomes part of the polarized VF state. The photon released in the decay of that state is identical to the incident photon. (See the "Appendix".) The mathematical result is that the calculated expression for the permittivity ǫ 0 of the vacuum does not depend on the frequency ω of the incident photon, implying that the behavior is not resonant.
After a parapositronium VF interacts with a photon, the photon-excited parapositronum VF is much more stable than the original parapositronium VF. From (10) the average lifetime of a parapositronium VF is ∆t p−P s ∼ = 3.2 × 10 −22 s. In contrast, from footnote 3 or (A.26), the mean lifetime of photon-excited parapositronum VF is 1/Γ p−P s ∼ = 6.2 × 10 −11 s. The photon-excited parapositronium VF is much more stable than a parapositronium VF and is known as a quasistationary state [24].
The energy, momentum and angular momentum of the photon are all transferred to the parapositronium VF-photon combination during the interaction. Such a state with an added photon has been has been described by Dodonov [46]. When the parapositronium atom vanishes back into the vacuum, a photon is emitted that is identical to the incident photon as described by the QED calculation in the "Appendix".
As (7) indicates, in the center of mass system the parapositronium VF can be described as a one-dimensional harmonic oscillator. In that description, in addition to imposing the Heisenberg uncertainty principle of (9) as an equality, because the harmonic oscillator is one-dimensional there is a second Heisenberg uncertainly principle to which the equality sign can be applied, producing a state of minimum uncertainty [41] for the wave packet ψ(x) of the oscillator. In particular, the onedimensional wave packet will not change in character during the time for which the VF exists. When the photon interacts with the VF, a polarized, photon-excited VF is created, delaying the progress of the photon until the polarized, photon-excited VF decays. Note that the polarization of bound, charged lepton-antilepton VFs by the presence of a real test photon, which has just been calculated in (21) for parapositronium, is different from vacuum polarization, a subject originally discussed by Serber [47] and Uehling [48] and treated in many texts [30,49,50,51]. Vacuum polarization results from the presence of virtual, charged lepton-antilepton pairs created from either virtual or real photons and is calculated using perturbation theory in quantum electrodynamics. Feynman diagrams for vacuum polarization are, for example, given in [51]. Although vacuum polarization is a separate process from the polarization of a VF, vacuum polarization must be taken into account when considering charged lepton-antilepton VFs, in that it is the charge renormalized by vacuum polarization that constitutes the charge that is indirectly observed.

Calculation of dielectric properties of the vacuum: permittivity of the vacuum
In a dielectric [52] the electric displacement D(t) satisfies where In (24) ǫ is the permittivity of the dielectric, and P (t) is the polarization density. In (25) N j is the number of oscillators per unit volume of the j th variety that are available to interact, and p j (t) is the average dipole moment of the j th variety of oscillator. The polarization density P (t) is responsible for the increase from ǫ 0 E(t) to ǫE(t) because of photons interacting with oscillators in the dielectric and results entirely from polarization of the atoms, molecules or both in the dielectric. It then follows that in the vacuum ǫ 0 E(t) must result entirely from the polarization density P VF (t) of atoms, molecules, or both that are VFs. Thus, Using (25) and then (22), (26) becomes The three types of VFs considered first are bound states of a charged lepton and antilepton, namely, parapositronium, muon-antimuon bound states, and tau-antitau bound states. Quark-antiquark states will be discussed later. Again, initially attention is restricted to parapositronium that is a VF. The progress of a photon traveling through the vacuum is slowed when it interacts with and has a polarizing effect on a VF consisting of a charged lepton and antilepton bound into an atom in its ground state. Isolated ordinary matter is kinematically forbidden from absorbing and then reemitting a photon. But an isolated bound state of a charged particle-antiparticle VF can be polarized by a photon and then, when the polarized VF annihilates, emit an identical photon, returning to the vacuum the energy originally borrowed for the creation of the VF. From the "Appendix", after reinstating factors of and c, the electromagnetic decay rate Γ p−Ps for a VF parapositronium atom after it has interacted with the incident photon to form a quasi-stationary state is The decay rate of ordinary parapositronium into two photons is half the above rate. [53,25].
The quantity e −Γ p−Ps t is the probability that a photon-excited parapositronium VF has not decayed electromagnetically during a time t, and 1 − e −Γ p−Ps t is the probability that it has decayed. The number density of parapositronium VFs with which a photon actually interacts is N V F j in (27). The rate for a polarized parapositronium VF to annihilate and emit a photon equals the rate for a parapositronium VF to interact with a photon and form a polarized VF state. As a result the quantity 1 − e −Γ p−Ps t is the probability that a parapositronium VF absorbs a photon during a time t.
The quantity N V F j for a parapositronium VF, denoted N p−Ps , equals the number density of parapositronium VFs times the probability that a parapositronium VF will interact with an incoming photon during the lifetime ∆t p−Ps : Since Γ p−Ps ∆t p−Ps ≪ 1, 1 − e −Γ p−Ps ∆t p−Ps can be replaced by Γ p−Ps ∆t p−Ps , with the result Substituting (6) and (30) into (27), Note that the mass of the electron has cancelled from the expression for ǫ 0 , implying that bound muon-antimuon and tau-antitau VFs each contribute the same amount to the value of ǫ 0 as parapositronium VFs. Including the contributions from the three types of bound, charged lepton-antilepton VFs yields Quark-antiquark VFs contribute little to the value of ǫ 0 : for the heavy quarks Q = c, b, or t it is appropriate to think in terms of static quark potentials for QQ bound states. An analysis similar to that leading to (31) shows that the contribution to ǫ 0 from η c (1S), the least massive cc bound state that has J = 0 and positive charge conjugation parity [54], is at least 10 −4 times smaller than the combined contribution from the three charged lepton-antilepton VFs. The corresponding bb and tt bound states are estimated to contribute even less. Consider the light quarks q = u, d, and s: the π 0 , η, and η ′ are the least massive qq bound states that have J = 0 and decay into two photons. Since the strong interactions are primarily responsible for the binding of these relativistic states, qq bound state VFs would have much higher natural frequencies than the electromagnetically bound, charged lepton-antilepton VFs. Accordingly qq bound state VFs contribute little to ǫ 0 and to lowest order need not be considered.
Ignoring the small contributions to ǫ 0 from quark-antiquark VFs, an approximate formula for ǫ 0 is immediately obtained from (32) using the defining formula α = e 2 /(4πǫ 0 c) to eliminate α and then using c = 1/ √ µ 0 ǫ 0 : Only the lowest-order terms in α have been retained in calculating ǫ 0 . The binding energy of parapositronium was neglected when calculating ∆t p−P s in (10), and only the leading term was retained when calculating ω 0 j . Moreover, only the leading term in the formula for Γ p−Ps has been included.

Permeability of the vacuum
In (2) a formula for the speed of light c was calculated using (1) and c = 1/ √ µ 0 ǫ 0 , which is obtained from Maxwell's equations provided that the electric displacement D = ǫ 0 E, where E is the electric field; the magnetic field B = µ 0 H, where H is the magnetic field intensity. However, in a medium [52] where M is the magnetic moment density of the medium. In the fermion quantum vacuum the medium consists of VFs. Parapositronium has no magnetic moment resulting from the orbital motion of the electron and positron because parapositronium is a singlet state, the electron and positron have equal masses and the charges on the electron and positron are equal in magnitude and opposite in sign. Also, for parapositronium the expectation value of the combined magnetic moments of the electron and positron is zero [55]. Thus for parapositronium M = 0. (Similarly M = 0 for muon-antimuon, tau-antitau and quark-antiquark VFs.) Thus, in spite of the presence of VFs, the vacuum satisfies the standard vacuum relation B = µ 0 H, justifying the statement in footnote 3 that for lepton-antilepton VFs, µ = µ 0 , and justifying the calculations of ǫ 0 and c in (1) and (2), respectively.

Results and Discussion
In the model of the the vacuum discussed here, VFs of bound, charged leptonantilepton pairs are polarized by photons much the way that ordinary matter is polarized. As a consequence, the permittivity of the vacuum (33) and then c (2) are calculated. The calculation of c presented here automatically satisfies the condition that c is the same in any direction in every inertial frame. This eliminates the need for one of the two postulates on which special relativity is based.
In the early universe when the temperature was sufficiently high that it was difficult for charged lepton-antilepton VFs to either bind into atoms or remain bound once they had formed atoms, the number density of charged lepton-antilepton atoms that are VFs would have been less than today. From (27) it then follows that ǫ 0 would also have been smaller. Since c = 1/ √ µ 0 ǫ 0 , a decrease in the value of ǫ 0 would increase the speed of light. Variable-speed-of-light cosmology, for example, has been discussed, originally by Moffat [56] and later by Magueijo [57]. The finite temperature and density corrections to vacuum polarization must be taken into account when describing the early universe [58]. The results obtained here should provide additional insight into such investigations.
Various books have summarized physicists' significantly expanded understanding of processes involving photons. These include the books by Milonni [11], Meis [13], Leonhardt[12], Quantum Optics by Garrison and Chiao [59], and Quantum Optics by Grynberg, Aspect and Fabre [60]. There is, however, much yet to be learned. In particular, a further analysis is needed of fluctuations of the photon vacuum, carried out along the lines presented here for fermion-antifermion VFs.
A Calculation of the decay rate of a photon-excited, parapositronium vacuum fluctuation Here the decay rate is calculated for a photon-excited, VF of parapositronium. Single-photon decay is kinematically forbidden for a photon-excited, quasi-stationary state of ordinary parapositronium but is allowed for a photon-excited, quasi-stationary state of parapositronium that is a VF since its energy and momentum are transitory and thus do not enter into overall energy-momentum conservation. The formula for the decay rate immediately generalizes to yield decay rates for VFs of muon-antimuon and tau-antitau that are bound into atoms in their singlet, ground states.
Labeling the initial (incident) and final (emitted) photons, respectively, by γ i and γ f , to lowest order the two Feynman diagrams that contribute to the process γ i +positronium that Fig. 1 (a) Photon γ i interacts with a positron that then annihilates with an electron, emitting photon γ f . (b) Photon γ i interacts with an electron that then annihilates with a positron, emitting photon γ f . four-momenta of the electron, positron, initial photon, and final photon.
As already mentioned, for ordinary positronium the process is kinematically forbidden. In the center-of-mass rest frame of positronium, p − + p + = 0. Therefore, in this frame, where me is the mass of an electron. Conservation of energy and momentum requires p + + p − + k i = k f . Squaring both sides of the above equation yields Eq. (A.2) cannot be satisfied for ordinary positronium since both terms on the left-hand side are positive. However, after a photon excites positronuim that is a VF, a photon can be emitted, but only when the positronium vanishes into the vacuum because only then does E − → 0, allowing (A.2) to be satisfied. When performing an electrodynamics calculation, if a factor exp±(ip · x) is associated with a particle that is part of a VF when it appears in its initial state, a factor exp∓(ip · x) is associated with a particle that is part of a VF when it vanishes into the vacuum. This just eliminates the contribution of the particle that is part of a VF to overall energy-momentum conservation. When progressing along an energy-momentum line in a Feynman diagram, the energy-momentum associated with a particle that is part of a VF is not further used after the particle vanishes. Since the parapositronium atom is a VF, it is on shell. That is, as is the case for physical parapositronium, the electron and positron are on shell and they are bound by 6.8 eV, the binding energy of parapositronium. Using the notation of Ref. [49], the S-matrix for the transition photon + positronium that is a VF → photon is 8 In (A.3) the fermion wave functions are normalized to unit probability in a box of volume V . The equality 1/( / p ± me) = ( / p ∓ me)/(p 2 − m 2 e ) is used to rewrite the above two propagators. From the identity {γ µ , γ ν } = 2g µν I, the equation immediately follows where a and b are four-vectors and I is the identity matrix. Using (A.4) To obtain the decay rate to lowest order, in S f i the respective velocities v − and v + of the electron and positron can be neglected. Thus The respective polarization vectors ǫ i and ǫ f of the initial and final photons are chosen to be space-like: Using (A.6), (A.7), and k i · k i = 0, With the aid of (A.8), (A.5) becomes Since k i = k f , it follows from (A.4) and (A.7) that / k i anti-commutes with both / ǫ i and / ǫ f , allowing (A.9) to be rewritten as Then, In (A.11) the interaction is assumed to occur in the time interval −T /2 < t < T /2. Summing over the electron and positron spins, which converts |S fi | 2 into a trace denoted by T r, Using (A.4) to reverse the order of / k i / p − and using / k i / k i = k i · k i I = 0, the following term that appears in the second line above can be simplified: As mentioned previously, / k i anti-commutes with both / ǫ i and / ǫ f , with the result that (A.12) becomes (A.14) Eq. (A.4) is used as necessary to reverse the order of / ǫ i and / ǫ f so as to obtain terms of the form / ǫ i / ǫ i = ǫ i · ǫ i I = −ǫ i · ǫ i I = −I and / ǫ f / ǫ f = −I. All nonzero traces in (A.14) are then either of the form T r(/ a/ b) or T r(/ a/ b/ c/ d) that are easily simplified, yielding the result The average cross section is now calculated for photon-excited positronium that is a VF to annihilate and emit a photon: |S f i | 2 is divided by V T to form a rate per unit volume, divided by the electron-positron flux |v + − v − |/V , and divided by the number of target particles per unit volume 1/V . Averaging over the spins of the electron and positron ( 1 4 s ± ), summing over the polarizations of the final photon ( ǫ f ), averaging over the polarizations of the initial photon ( 1 2 ǫ i ), summing over the number of states of the final photon in the momentum interval d 3 k f ( V d 3 k f /(2π) 3 ), and averaging over the number of states of the initial photon in the momentum interval d 3 k i ( V d 3 k i /(2π) 3 ), an expression for the cross section σ is obtained: Using (A.15), (A.16) becomes Choosing the z-axis to point in the direction of k i , the unit polarization vectors for the initial photon ǫ a i and ǫ b i are chosen in the x-and y-direction, respectively. Because the delta function in (A.17) imposes the condition k f = k i , the unit polarization vectors for the final photon ǫ a f and ǫ b f can also be chosen in the x-and y-direction, respectively. The sum over polarizations in (A.16) is now easily performed: it is straightforward to show that The theta function θ(k i0 ) = 0 if k i0 0 and θ(k i0 ) = 1 if k i0 0. With the aid of (A.20), the second line in (A.17) can be rewritten as Factoring k 2 f = k 2 f 0 −|k f | 2 in the above δ-function, using (A.19), rewriting d 3 k f as dΩ f |k f | 2 d|k f |, performing the angular integration over dΩ f , which yields a factor of 4π, and then integrating over |k f |, From the formula for the cross section, a formula for the decay rate is readily obtained. The logic is the same as that used to calculate the decay rate for parapositronium decaying into two photons [53,25]: parapostronium, orthopositronium, and a photon have respective charge conjugation parities of +1, -1, and -1. Thus photon-excited parapositronium has charge conjugation parity of -1 while photon-excited orthopositronium has charge conjugation parity of +1. Since electromagnetic interactions are invariant under charge conjugation, photon-excited parapositronium, but not photon-excited orthopositronium, can decay into a single photon.
In obtaining (A.23) the electron and positron spins were averaged over all four spins, resulting in the sum being divided by four. But the annihilating state is parapositronium, the singlet state. Orthopositronium, the triplet state, does not contribute. Since only one of the four spin states contributes to the cross section, the formula for the cross section should not have been divided by four, it should have been divided by the number one. Thus the formula for σ in (A.23) should be multiplied by a factor of four to obtain the cross section, abbreviated σ p−Ps , for the annihilation into a photon of photon-excited parapositronium that is a VF, For the annihilation of photon-excited parapositronium that is a VF into a photon, the electromagnetic decay rate Γ p−P s is calculated using the mechanism for the annihilation of ordinary parapositronium [53,25]. The Schrödinger wave function ψ(x) for parapositronium is just the ground-state hydrogen atom wave function with the reduced mass of hydrogen, which is approximately me, replaced by me/2, the reduced mass of parapositronium: In the above formula x is the magnitude of the vector x = xe − xp where xe and xp are, respectively, the positions of the electron and the positron. The decay rate Γ p−P s is the product of σ p−Ps and the flux of a parapositronium atom, which is the relative velocity of approach of the electron and positron in parapositronium multiplied by |ψ(0)| 2 , the probability density that the electron and positron collide and annihilate. It immediately follows that the corresponding decay rates for photon-excited, muon-antimuon or tau-antitau VFs bound into a singlet, ground state are obtained by replacing the electron mass in (A.26) with the muon or tau mass, respectively. As photons travel through the vacuum, these three decay rates characterize how photons interact with charged lepton-antilepton VFs.