A gauge field treatment of a current oscillating at frequency ν of interacting neutral atoms leads to a set of matter-wave duals to Maxwell's equations for the electromagnetic field. In contrast to electromagnetics, the velocity of propagation has a lower limit rather than upper limit, and the wave impedance of otherwise free space is negative real-valued rather than 377 Ω. Quantization of the field leads to the matteron, the gauge boson dual to the photon. Unlike the photon, the matteron is bound to an atom and carries negative rather than positive energy, causing the source of the current to undergo cooling. Eigenstates of the combined matter and gauge field annihilation operator define the coherent state of the matter-wave field, which exhibits classical coherence in the limit of large excitation.

This work considers an oscillating current of interacting ultracold atoms through the lens of non-relativistic field theory. Over the past few decades, relativistic quantum field theory has led to the unification of three of the four fundamental forces and now provides our modern picture of the physics of fundamental particles as the Standard Model.1 Non-relativistic systems, in particular the connection between Hamilton–Jacobi theory and quantum mechanics, have been studied using a field-theoretic approach since the early days of quantum mechanics.2,3 In line with these early works and following the reasoning of Wheeler,4 the problem of interacting identical neutral particles can be treated using classical field theory; utilized in conjunction with a Hamilton–Jacobi formalism, the approach leads naturally to a wave description of particles and their interactions. Indeed, it was a significant revelation of the last century that Maxwell's equations themselves follow from the treatment of electrons as excitations of a field taken together with certain symmetry considerations.5,6 Our work draws on electromagnetism as a familiar application of classical field theory and offers a new as well as useful perspective on matter-waves.

The wave characteristics of quantum-mechanical particles, either individual or ensembles, are commonly referred to as “matter waves;” they are characterized by a de Broglie wavelength λ = h / p, where h is Planck's constant and p is the particle momentum. An excellent example of matter-waves is the flux of atoms extracted from a Bose–Einstein condensate (BEC)7,8 that was first achieved in 1995 in a diluted gas of rubidium and sodium atoms.9,10 Early experiments demonstrated the quantum coherence of these systems.11–13 The achievement of BEC prompted many-body physics to shift its focus toward ultracold gases, resulting in various models, analytical techniques, and numerical methods, such as the Gross–Pitaevskii equation (GPE), time-dependent Hartree–Fock–Bogoliubov (TDHFB), Bogoliubov–de Gennes equation, quantum kinetic theory, and stochastic methods, among others.14–17 

In parallel with advances in atom cooling and manipulation techniques, so have there been advances in atom interferometry18–29 that take advantage of the properties of ultracold atoms, particularly their coherence. Atom interferometry and its use of matter-waves has been very much inspired by its optical interferometry counterpart. At the same time, the nature of BEC systems also leads one naturally to think in terms of atom currents flowing from one part of a system to another, thus arose the field of atomtronics, which refers to an atom analog of electronics in which atom flux and chemical potential substitute for electric current and potential.8,30–37

As fundamentally non-thermal-equilibrium open quantum systems, even rather simple atomtronic circuits prove challenging for standard many-body methods. Yet, we can look more deeply at the analogy between atomtronic currents vs electronic currents and matter-waves vs electromagnetic waves. Maxwell's equations inform us that an oscillating electric current gives rise to a coherent electromagnetic wave. We can ask, therefore, whether an oscillating atomtronic current gives rise to a coherent matter-wave. The answer to this question sets the motivation for our work here. We will find that a field-theoretic description of matter-waves provides an additional set of tools in which to treat atomtronic circuits. The approach leverages the analytical, heuristic, and numerical tools that have been well-developed for electromagnetics. In particular is the appearance of “Maxwell matter waves”—the classical limit of coherent matter-waves governed by matter-wave duals to Maxwell equations. Like their electromagnetic counterparts, Maxwell matter-waves exhibit temporal coherence, and they have several other aspects that are familiar from electromagnetics.

Here, we consider an oscillating current of neutral atoms, such as 87Rb, which repel each other due to mutual van der Waals interactions. In field theory, atom interactions are accounted for by the introduction of a gauge field. In comparison to the Coulomb forces between electrons, van der Waals forces are weak, and they are short-ranged. One might conclude that, therefore, the gauge field can be dismissed if, say, particles are far apart. However, that conclusion proves to be misleading if not wrong. The gauge field that embodies the interaction between neutral atoms is associated with energy and, generally, the transmission of power. At low temperatures, atoms interact through s-wave collisions, and, thus, the gauge field that characterizes these interactions assumes a relatively simple form. We will sometimes refer to this field as the “matteron field” to provide it with an identity. The matteron field, then, accounts for the interactions among the cold neutral atoms.

With a motivation that parallels the electromagnetic case, we are particularly interested in the notion of a single-mode of the matteron field. Generally, a finite temperature gas, or even a pure BEC, will be comprised of a continuum of modes.

Our theoretical development is carried out in three sections. The first is a classical, field-theoretic derivation of the family of Maxwell equations characterized by a few formally introduced constants, where the classical theory does not constrain the dynamical parameters of wave propagation. While assuming the constants to be the speed of light and free-space impedance leads to the well-known Maxwell equations, the same constants with different meaningful values can represent matter-wave duals to Maxwell's equations. These values are provided by a quantum-mechanical treatment of the fields, which is carried out in Sec. III. This section also establishes the formal connection between the quantum-mechanical coherent state and the classical coherence of matter-waves implied by Maxwell equations duals. With the dynamical parameters in hand, Sec. III returns to the classical treatment, giving the fully constrained Maxwell equation duals from which a quantitative correspondence between experiment and theory can be established. We close with Sec. V to highlight the similarities and differences between coherent matter-wave and electromagnetic wave theories.

Our development of classical, non-relativistic field theory follows closely the excellent notes of Wheeler,4 which include substantially more discussion concerning the development than we provide. Typically, modern physics treatments of field theory are written to address the domain of high-energy physics, so the constraint of Lorentz covariance and its tie to the speed of light are imposed early on.38–41 In the realm of ultracold atomic physics, particle velocities of 10 m/s are already very high, placing us in a non-relativistic realm.

As a starting point, let us consider a Lagrangian density4 characterizing a flux of particles having mass m,
(1)
in which S is Hamilton's principal function, U ( x ) is an applied particle potential, R has units of ( length ) 3 and so describes the energy carried by the flux as a density, and i is a partial derivative over time t or spatial coordinates i { x , y , z }. (Our symbol choice R for density follows the notation of Wheeler;4 the perhaps more natural ρ is reserved for a later analog with electromagnetism.)
As an aside, we note that the Hamilton–Jacobi equation derived from the Lagrangian for S is
(2)
Considering for the moment a single particle, its momentum p = S. So one can recognize in the above its relationship to the single-particle Hamiltonian,
(3)
The Lagrangian Eq. (1) is invariant under the global gauge transformation,
(4)
where the constant has dimensions of action and w is dimensionless.
Establishing local gauge invariance begins with the introduction of vector and scalar gauge fields A and ϕ and incorporates the replacements,
(5)
Contriving to have the resulting equations as well as their units look familiar, here we have introduced in a formal way a “charge” q. In electromagnetics, the charge indicates the strength of the particle interactions; here, the choice to include the symbol is gratuitous, since the interaction strength per particle or per mass could equally well be incorporated into the units for the gauge field. The Lagrangian density is re-written as
(6)
The new Lagrangian is invariant with respect to a local gauge transformation characterized by W ( x , t ),
(7)
The Lagrangian Eq. (6) accounts for the energy of the particle flux as well as the energy of the interaction between the particles and the auxiliary fields A and ϕ. We need to also incorporate the energy associated with the auxiliary fields by themselves. To that end, we introduce the matter-wave duals F and G of the electric field E and magnetic field B,
(8)
(9)
The fact that these fields are gauge-invariant, as is the case for their electromagnetic duals, is easy to verify. With these fields so defined, we can write the combined Lagrangian density written in a manner contrived to be later familiar,
(10)
We will have more to say about the matter-wave duals ξ0 and η0 to the electric permittivity and permeability in Sec. IV.
The Hamilton–Jacobi equation for the field characterizing the density is
(11)
This equation for the density expresses a continuity relation for the Noetherean current.41 Define charge and current densities,
(12)
(13)
We then have a familiar continuity relation,
(14)
Since the divergence of a curl is identically zero and given the definitions of the fields F and G,
(15)
(16)
The Hamilton–Jacobi expression for the scalar potential is
(17)
which leads to
(18)
We have more to say about this Gauss' law dual in Sec. V. Continuing each component of the vector potential has a corresponding Hamilton–Jacobi equation, such as
(19)
Combining the spatial components and following the definitions of the fields and current lead to a final Maxwell equation,
(20)

We have now a complete set of Maxwell's equations: Eq. (18) as Gauss' law for the electric field, Eq. (15) as Gauss' law for the magnetic field, Eq. (16) as Faraday's law of induction, and Eq. (20) as Ampère's circuit law with the addition of Maxwell's displacement current, along with the continuity equation Eq. (14). In Wheeler's words, in fact, we have a family of Maxwell's equations.4 Substituting v 0 c , η 0 μ 0 , ξ 0 ε 0 and, of course, F E and G B returns us to electromagnetics in familiar nomenclature. The charge and current densities are sources for the electric and magnetic fields. Equivalently, we can view them as the source for the scalar and vector potentials.

Let us introduce the matter-wave dual of the Lorenz gauge,
(21)
In this case, a pair of wave equations govern the potentials,
(22)
(23)
The velocity v0 substitutes for the speed of light c, which is an experimentally determined fundamental physical quantity. The matter-wave dual will become evident as we blindly follow the course as we would for electromagnetics. In particular, we are now in a position to consider an oscillating matter current. Consider the positive frequency component of a one-dimensional oscillating current,
(24)
in which case, the continuity equation, Eq. (14), insists
(25)
along with
(26)
where we define the velocity
(27)
While v m is the phase velocity associated with the current density, thinking forward to the quantum aspects, we can know in advance that it is also the group velocity of the particles.36 Let us correspondingly write
(28)
(29)
and further define a refractive index n
(30)
Then, the amplitudes for the potentials are easily calculated,
(31)
(32)
Borrowing once again from electromagnetics, we write a matter-wave dual to Ohm's Law,
(33)
where the basic wave number k 0 k / n and the impedance is
(34)
The Theory of Relativity regards the speed of light as an upper speed limit. We see through the impedance that there is also a speed limit: as v m approaches v0, so that the refractive index approaches unity, the field amplitudes diverge. The minus sign associated with the impedance has been explicitly brought out because, as we shall see in Sec. III, the matter-wave speed limit is a lower bound rather than an upper bound, and the refractive index is typically less than unity rather than greater than unity as is the usual case for electromagnetics. The fact that the impedance is negative has important physical significance, namely, the power carried by the field,
(35)
is negative if the impedance is negative. Generally, the concept of negative real-valued impedance is familiar in electronics, and there are a few interpretations that are equivalent. We conclude that the power carried by the matter and field together is less than that carried by the matter alone. The relationship between power, particle flux, and current is further discussed in Sec. IV.

To provide some insight into the formalism, let us consider a classical picture of particle dynamics as illustrated in Fig. 1. Plotted along the x-axis is the particle potential energy due to an applied potential. In the case of alkali atoms, sculpting a particular potential is conceptually simple to implement, say, utilizing laser light tuned to the red or blue side of an atomic resonance to either raise or lower, respectively, the atomic potential energy. An ensemble of atoms is initially on the left, biased with high potential energy. They are shoved from the left toward the right, say, by the left-hand wall oscillating at the frequency ν. The atoms fall down the potential, propagate along a relatively flat region, and then climb back up the potential on the right. Of course, they have gained speed as they roll down the potential, and, therefore, they are far apart and then move closer together again as they climb the potential.

Fig. 1.

An illustration of the classical physics underlying gauge fields. Interactions among otherwise neutral atoms through van der Waals forces can be significant when atoms are close together, as in the left side of the figure. As atoms roll down the potential landscape, their speed increases and so does the inter-atom distance, so that the interaction forces become very small; nevertheless, they must “remember” the consequences of the earlier interactions. An oscillating current imposed on the left will be preserved even as the atom flux moves from left to right, so that it is coherent across the system.

Fig. 1.

An illustration of the classical physics underlying gauge fields. Interactions among otherwise neutral atoms through van der Waals forces can be significant when atoms are close together, as in the left side of the figure. As atoms roll down the potential landscape, their speed increases and so does the inter-atom distance, so that the interaction forces become very small; nevertheless, they must “remember” the consequences of the earlier interactions. An oscillating current imposed on the left will be preserved even as the atom flux moves from left to right, so that it is coherent across the system.

Close modal

Van der Waals forces are at play only when atoms are close together; therefore, at the bottom of the potential, their interactions are much weaker than at the top, where they are close together. Yet, both the applied and the interaction potential are conservative, so the atoms will “remember” their interaction energy such that they return to the same configuration on the right as they began on the left (more or less, since we have not given detail on how the walls work). If one imagines a continual stream of atoms from the left, the average flux will everywhere be the same, as is the current oscillating everywhere along the system at frequency ν. It is the gauge field that serves as the memory that invokes coherence of the current across the system.

The field theory we have applied treats particles in terms of a delocalized field. This field “knows” about the potential everywhere (say) in the + x direction, and, although the field is everywhere in the half-space, it is also propagating in the + x direction. In particular, in some sense, its current position and velocity are dependent on a distribution of the possible past positions at historical times. The view of particles as delocalized entities suits well a transition to a quantum-mechanical description. The notion of a short-range interaction is meaningful in a picture of localized particles undergoing a collision, yet when the particles are delocalized and, in fact, occupy an entire half-space, their interaction is manifest differently: if a steady oscillation is occurring in one location, one should expect it to occur in all locations, loosely in keeping with the depiction in Fig. 1.

As a practical matter, absent from the matter-wave duals of Maxwell's equations are values of the two key parameters, from which all dynamics follow. In electromagnetics, the vacuum speed of light and impedance (or the vacuum permeability and permittivity) are taken as fundamental constants of nature and are determined empirically (as is the charge of the electron). To determine the dynamical constants of our field theory, we must turn to a quantum-mechanical description, in which Hamiltonian steps in for the Lagrangian and operators step in for dynamical coordinates and their corresponding functions and functionals.

The classical theory anticipates that the potential energy associated with the gauge fields is negative. Negative energy associated with quantization of the matter-wave gauge field is in stark contrast with the electromagnetic case. While the physics itself is straightforward, conceptual navigation through negative energy territory is best done with some care. We, therefore, begin with the basics.

The transition from the classical to the quantum description of matter is typically accomplished with the introduction of massive particle creation and annihilation operators b ̂ m and b ̂ m, together referred to as “ladder” operators.42 They have the commutation property,
(36)
In contrast to the classical picture illustrated in Fig. 1, the creation operator b m creates a delocalized mode that exists everywhere in a half-space, while it propagates in the + x direction. For the purpose of this section, it can suffice to suppose that a mode propagates with a well-defined wavenumber, meaning that the particle potential is uniform. The operator treatment, however, equally well accommodates a spatially non-uniform potential, such as that depicted in Fig. 1. The Hamiltonian corresponding to a single-mode of the matter field is given by the number operator,
(37)
for which the Planck–Einstein frequency is
(38)
where km is the deBroglie wavenumber associated with the particle momentum and corresponding particle group velocity,
(39)
Eigenstates of the Hamiltonian are Fock (number) state. In particular, the ground state energy is zero,
(40)
Higher-lying Fock states are produced by repeated application of the creation operator,
(41)
while the energy associated with a single-mode is simply given by the number of particle excitations Nm of the matter field,
(42)

The Hamiltonian Eq. (37) embodies the energy of a matter field, which is of the particles, alone. It does not account for the particle interactions. The classical development likewise began with the Lagrangian for the particles alone, then treated interactions through the introduction of a gauge field, and, finally, considered the energy of the gauge field alone. The quantization of the gauge field, i.e., the matteron field, is similar to that of the matter field. We introduce the ladder operators a ̂ f and a ̂ f, which obey commutations relations the same as those above for the matter, while they commute with each other, [ a ̂ f , b ̂ m ] = [ a ̂ f , b ̂ m ] = 0, etc.

The Hamiltonian associated with the matteron field, like the electromagnetic field, is itself analogous to that of the quantum-mechanical harmonic oscillator associated with harmonic frequency νf (for the time being, we take the oscillator frequency to be general, ν ν f). A notable difference between the matter field and the matteron field is that the ground state of the latter has nonzero energy. Moreover, the energy associated with the gauge field is negative; we comment further on this in Sec. V. Glauber has pointed out that an inverted harmonic oscillator can be described by ladder operators identical to the normal oscillator but is characterized by negative rather than positive energies,43 
(43)

The quantum of excitation of the matteron field is the dual of the photon of electromagnetics: a massless gauge boson. We have referred to the dual elsewhere as a “matteron,”36 hence the naming of this gauge field. In contrast to the photon, the matteron is associated with negative energy.

As is the case for the matter field, the energy eigenstates of the matteron field are Fock states. For a field comprised of exactly Nf matterons,
(44)

As we move to explicitly consider particle interactions, there is a subtle but important departure from the case of electrons interacting through the electromagnetic field. In the latter case, we can, philosophically at least, separate the photons from the electrons and then consider how they interact. We saw from the classical treatment that an oscillating particle current serves as a source for the gauge field-one cannot exist without the other. For the corresponding quantum case, an excitation of a matter field is necessarily accompanied by an excitation of the matteron field. We are, thus, interested in the combined matter and matteron field excitations—in particular, what is appropriate to call “a single-mode excitation of the matter-wave field.” Hereafter, we shall assume the excitation is one-for-one, i.e., Nm = Nf. In contrast to the familiar deBroglie matter-wave case, there are two distinct frequencies involved: the Planck–Einstein frequency ωm associated with the particle and the frequency νf associated with the oscillating field, which is set by some external agent (a circuit, say), applying a time-varying force and causing a propagating oscillating current.

We are now in a position to self-consistently determine that the limiting velocity v0 is a lower bound on the particle velocity: First, we can understand that the dynamics of the system must be such that the total energy is non-negative, meaning the particle energy must be greater than or equal to the matteron energy, since a negative energy would imply that the system lies in a bound state. This, in turn, means that the minimum velocity is dependent on the frequency νf associated with the gauge field,
(45)
Second, we see that physically, the minimum velocity is such that the particle is never moving in the negative direction as it oscillates. Doing so would contradict the assumptions, Eqs. (24) and (25), that the charge density and current density are described by plane waves propagating in the positive x direction. (On the other hand, particles moving in the negative x direction are associated with plane waves of current and charge density propagating in the negative x direction.) This is the rationale for clarifying the distinction between particle flux and particle current in the IV. Having a classical picture of an oscillating particle in mind, its motion, and hence particle flux, is always in the positive direction as long as its center of mass is moving sufficiently fast in the positive direction. Evidently
(46)
In contrast to electromagnetics (and optics in particular), our refractive index is less than unity, n ω ν < 1, and often much less. We note that at the minimum particle velocity, i.e., n = 1, the wavelength λ = 2 π / k = 2 π / ( n k 0 ) associated with the field becomes equal to the particle de Broglie wavelength λ = n 2 π / k 0.
Given the multiple contexts of the noun “field”, we clarify that we will utilize the operators b ̂ with the subscript m or another subscript to refer specifically to the matter field, and the operator a ̂ f specifically with the subscript f to refer to the matteron (gauge) field. The “matter-wave field” will refer to the two taken together. From here on the quantum-mechanical treatment, we will focus onto the single-mode matter-wave case. That means we shall fix the two frequencies, ν f ν and ω m ω, and drop subscripts when there is no ambiguity. Formally, matter and gauge fields occupy distinct Hilbert spaces. In the single-mode case, however, the two are so tightly entangled that for the purpose of this work we can treat them together. We, thus, use a ̂ ν , a ̂ ν as matter-wave field ladder operators,
(47)
In particular for a Fock state, N f N ν,
(48)
Next, we seek to formalize a means of exciting the matter-wave field from some ground state. In order to allay concerns of creating mass from vacuum, for simplicity, we have a reservoir that has no matteron energy ( ν f = 0) but provides a supply of massive particles having energy ω. We are thinking, in particular, of a Bose-condensate having a large and fixed number Nr of identical particles as such a supply (see Fig. 2). We set as the system single-mode ground state,
(49)
Fig. 2.

The potential energy landscape for a coherent matter-wave emitter. The source well contains a BEC that serves as a reservoir supply of particles. For the purposes of this work, the gate well can be treated as a tuned element that selects the frequency ν of the oscillating current. This configuration of potential energy and atoms is a self-oscillating system.36 The matter-wave emission is due to system dynamics rather than an imposed oscillating force as notionally depicted in Fig. 1. Reproduced with permission from Anderson, Phys. Rev. A 104, 033311 (2021). Copyright (2021) under a Creative Commons License.

Fig. 2.

The potential energy landscape for a coherent matter-wave emitter. The source well contains a BEC that serves as a reservoir supply of particles. For the purposes of this work, the gate well can be treated as a tuned element that selects the frequency ν of the oscillating current. This configuration of potential energy and atoms is a self-oscillating system.36 The matter-wave emission is due to system dynamics rather than an imposed oscillating force as notionally depicted in Fig. 1. Reproduced with permission from Anderson, Phys. Rev. A 104, 033311 (2021). Copyright (2021) under a Creative Commons License.

Close modal
As a condensate, we take all particles to having identical energy E = ω. Let us refer to the two sub-spaces comprising our system in terms of reservoir states (r) and matter-wave states (ν), as in practice, the former will involve trapped particles and the latter particles (massive and matterons) that propagate in free space. We introduce joint operators,
(50)
(51)
such that
(52)
and, for example,
(53)
such that the reservoir gives up a particle to the matter-wave field. Higher excitation of the field is accomplished with repeated action of c ̂ r ν . The joint operators entangle the reservoir and propagating fields. In the limit that the reservoir particle number is large, it can be treated as a constant, and we can drop the reservoir terms to recover the matter-wave Hamiltonian, Eq. (47).

Total system energy is, or course, not conserved, since for every particle contributed by the reservoir, ν of energy is removed by the gauge field. We keep in mind that the emitted massive particles themselves each carry the full ω worth of energy of the particle contributed by the reservoir—it is the negative energy carried by the matteron that is responsible for the energy shortage. This indicates that the emission of the matteron field itself induces cooling of the reservoir, although the emission of the particles proves to cause heating44 such that the next effect is heat added to the reservoir. We have omitted the detailed physics of these processes.

The familiar classically coherent wave corresponds to the quantum mechanically pure coherent state given as the eigenstate of the annihilation operator,
(54)
in which the eigenvalue is complex, α ν = | α ν | e i φ ν, and the phase is associated with the oscillation frequency ν. The expectation of the energy carried by the coherent state is
(55)

Our theoretical development has involved no reference to the origin of the particle interaction—only that it exists and can be characterized in terms of particle mass m and charge q. The pivotal requirement is that the interaction energy is sufficient to supply a matteron having energy ν. In Fig. 2, this corresponds to a BEC having an interaction mean-field contribution to its chemical potential μ B ν in order to supply matterons to the emission.

Finally, noting the role of impedance in Sec. II, we introduce the quantum impedance
(56)
In closing, we note that labeling of the charge q is convenient for the sake of analogy and for dimensional bookkeeping. In practical applications, however, one will find it convenient to set the charge either to the particle mass or to unity, corresponding to the impedance given in units of (energy-per-mass)-per-(mass-per-second) or in units of (energy-per-particle)-per-(particle-per-second).

The analogy with electromagnetics can provide insight into the matter-wave system through a simple radio frequency circuit. Figure 3 shows an oscillating microwave voltage source delivering power to a load RL through an electromagnetic transmission line.45 The voltage and current waves carry energy, and if the source, load resistance, and transmission line impedance are equal, all the energy is delivered to the load. Of note is the fact that the electric current oscillates between negative and positive values and, thus, flows in both directions. The voltage oscillates in phase with the current; as a consequence, the power flows continuously from left to right.

Fig. 3.

A microwave transmission line circuit is useful for highlighting the distinction between current and flux. The flux corresponds to the flow of power in the circuit. The maximum power P available from a microwave source VS is delivered to a load when the source resistance, load resistance, and characteristic impedance of the transmission line are all equal, R S = R L = Z c. In such a case, the power and the atom flux flowing from the source to the load are constant, while the voltage across and current through a given point along the line vary sinusoidally. In particular, the direction of the current through a point flows alternately toward the load and toward the source.

Fig. 3.

A microwave transmission line circuit is useful for highlighting the distinction between current and flux. The flux corresponds to the flow of power in the circuit. The maximum power P available from a microwave source VS is delivered to a load when the source resistance, load resistance, and characteristic impedance of the transmission line are all equal, R S = R L = Z c. In such a case, the power and the atom flux flowing from the source to the load are constant, while the voltage across and current through a given point along the line vary sinusoidally. In particular, the direction of the current through a point flows alternately toward the load and toward the source.

Close modal

The example highlights an important distinction between current and flux: We will use the term “flux” to refer to the flow of power, equivalently in the case of atomtronics the flow of particles, or quantum-mechanically to the flow of probability. As with the transmission line example, we consider a case in which ultracold atom flux is continuous, while atomtronic current and voltage are oscillating amplitudes. By contrast, a BEC whose center-of-mass is oscillating within a trap is described by an oscillating flux.

Impedance reflects the energy cost of adding an additional charge to the current. The quantum impedance determines the duals to permittivity and permeability, since, in analogy with electromagnetics,
(57)
(58)
This leads us to
(59)
(60)
Finally
(61)
Having the impedance and velocity in hand allows one to calculate all field amplitudes, given the current density amplitude J0. In particular, the potential is given by Eq. (33), while the remaining field amplitudes are
(62)
The fact that the gauge and matter fields propagate together leads one naturally to think in high-frequency electronic circuit terms, in which currents and fields also coexist. In this context, it is convenient to work using the duals of electric voltage and current. To that end, we introduce
(63)
(64)
where A is an effective area, over which the fields and currents are transversely confined. The relationship between the two is set by Zf, which, in this context, is referred to as the characteristic impedance45,46 such that
(65)
While the current I 0 is an amplitude, generally, it is the flux I m, defined as the number of massive particles per second that traverse a given position, which is straightforward to measure. We are interested particularly in the energy stored in and carried by the gauge field. The power carried by the gauge field is simply
(66)
Hence, in the limit n 1,
(67)
(68)

In circuit applications, it will often be convenient to consider the energy of the particles ω carrying the currents as fixed, while the oscillation frequency ν is taken to be variable. In such cases, it will be convenient to make the substitution k 0 k / n in the various formulas for fields and impedance.

Through field theory, we have derived a set of matter-wave duals to Maxwell's equations describing electromagnetic waves. We have intentionally cast the theory to resemble electromagnetic theory. While valuable from an intuition standpoint, doing so also obscures the contrasts between the two.

Conspicuously absent in the treatment are the duals to the constitutive relations corresponding to D = ε E and H = B / μ. In particular, for example, the gauge field energy in the Lagrangian Eq. (10) uses the duals to the vacuum permeability and permittivity. That is because here there are no duals to polarization or magnetization of a medium. The velocity associated with the gauge field is tied to that of the atoms, which, in turn, is influenced by an external potential that acts directly on the atoms. By contrast, with electromagnetism, the gauge fields E and B themselves interact with a medium in which they propagate.

We have shown that an oscillating matter current gives rise to a coherent matter-wave that obeys the Maxwell equations duals. Of dramatic difference from the electromagnetic counterpart is that the oscillating matter current does not produce a radiated field that exists in a source-free region. On the contrary, through the negative sign in the Hamiltonian Eq. (43) for the matteron field, we see that the matter and its gauge field are bound together. The quantization of the field introduces the matteron, a massless gauge boson that is dual to the photon. In a particle picture, the matteron is bound to the massive particle, while the classical particle can be pictured as undergoing longitudinal oscillation. It is worth noting that if the matteron energy were taken to be positive, it would not be bound. As a massless boson, the matteron would then necessarily travel at the speed of light (see, for example, Schwichtenberg41) and be indicative of a long-range force, not of the short-ranged van der Waals forces at play in our particle interactions.

It is reasonable to consider whether our matter-waves are, in fact, sound waves, and whether the matteron should, in fact, be called a phonon. Sound waves, such as those that can propagate in a BEC,14 arise from the interaction among the condensate atoms, and the propagation velocity depends upon, among other things, the atomic density. We have seen in the discussion at the end of Sec. IV that the interaction energy associated with a matter-wave as it propagates in otherwise empty space can be arbitrarily low, and that the wave propagation speed is given by the particle group velocity. So there is nothing acoustic about the matter-waves under consideration.

The picture of longitudinal oscillation of the massive particle is perhaps at first objectionable because oscillation implies a restoring force. While matterons are massless, they carry momentum k 0, and they are bound to the massive particle with energy ν. Thus, we can think of a massive particle attached to a lighter one with a spring that supplies a restoring force as the two oscillate with respect to each other.

The contrast between matter-waves and electromagnetics is further highlighted, for example, by the dual to Gauss' law, Eq. (18), which might suggest field lines emanating from an isolated charge as it does for the electric case. Yet, the dispersion implicit in the dual to permittivity ξ Eq. (60) belies the complexity under the hood, so to speak, particularly at DC. In short, at a distance away from any current, no matteron field is detectable. Yet, within the flux of particles, their effects are measurable: the fact that the energy of the bound particle-plus-matteron state is less than the energy of the massive particle alone reveals the physics behind the flow of negative power carried by the gauge field.

The dynamical parameters for matter-waves, namely, the propagation speed and impedance, are imposed by quantum mechanics rather than by relativity as is the case for electromagnetism. The quantum treatment of the fields provides a clear understanding of the meaning of “classically coherent matter wave” as an eigenstate of the combined matter-field-plus-gauge-field annihilation operator, which has a well-defined phase relative to the oscillation frequency ν of the current. Keeping to the classical domain, it seems appropriate to refer to these waves as “Maxwell matter waves.” In this context, our duals to Maxwell's equations were contrived to appear the same as the familiar set, yet the fundamental conclusions are that the propagation velocity has a lower rather than upper limit, and that the impedance of free space has a negative real value for Maxwell matter-waves, compared with the positive 377 ohms for electromagnetic waves.

At the same time, Maxwell matter-waves do obey their dual Maxwell's equations, and thus the heuristic and analytical tools that already exist for electromagnetics can be put to effective use to address matter-wave calculations. Indeed, one can work with the fields F and G, but it is also natural to work in a circuit theory context involving scalar potential and current, connected to each other through impedance, as we have seen. Step changes in the refractive index, for example, lead to familiar coefficients of wave amplitude and power reflection and transmission that arise in the case of particle wavefunctions, light waves, microwaves, and so on.45–47 

It is relevant to note that the interference of Maxwell matter-waves reveals substantially different character than the interference of de Broglie matter-waves. For example, the fringe spacing of interfering de Broglie waves will decrease with increasing particle velocity, while that of Maxwell matter-waves will increase.

Traditional is the de Broglie view in which thinks of matter classically as a particle phenomenon whose wave character is revealed by quantum mechanics. Here, we see that Maxwell matter-waves provide the complementary view in which matter can be viewed classically as a wave phenomenon whose particle aspects are revealed by quantum mechanics.

The work of K. Krzyzanowska was supported by the Quantum Science Center (QSC), a National Quantum Information Science Research Center of the U.S. Department of Energy (DOE). The work of D. Z. Anderson was supported by Infleqtion. We are grateful to S. Du and V. Colussi for valuable discussions.

D.Z.A. owns stock in ColdQuanta Inc. and serves on the Board of Directors of ColdQuanta Inc. (doing business as Infleqtion).

Dana Z. Anderson: Conceptualization (equal); Formal analysis (equal); Funding acquisition (equal); Investigation (equal); Writing – original draft (equal); Writing – review & editing (equal). Katarzyna Anna Krzyzanowska: Conceptualization (equal); Formal analysis (equal); Funding acquisition (equal); Investigation (equal); Writing – original draft (equal); Writing – review & editing (equal).

Data sharing is not applicable to this article as no new data were created or analyzed in this study.

1.
T.-P.
Cheng
and
L.-F.
Li
,
Gauge Theory of Elementary Particle Physics
(
Oxford University Press
,
New York
,
1984
).
2.
E. T.
Whittaker
,
Proc. R. Soc. Edinburgh, Sect. A
61
,
1
(
1941
).
4.
N.
Wheeler
, see https://www.reed.edu/physics/faculty/wheeler/documents/index.html for “
Classical field theory class notes
” (
1995
).
5.
D. E.
Soper
,
Classical Field Theory
(
Wiley
,
New York
,
1976
).
6.
J.
Franklin
,
Classical Field Theory
(
Cambridge University Press
,
New York
,
2017
).
7.
M. O.
Mewes
,
M. R.
Andrews
,
D. M.
Kurn
,
D. S.
Durfee
,
C. G.
Townsend
, and
W.
Ketterle
,
Phys. Rev. Lett.
78
,
582
(
1997
).
8.
C.
Ryu
and
M. G.
Boshier
,
New J. Phys.
17
,
092002
(
2015
).
9.
M. H.
Anderson
,
J. R.
Ensher
,
M. R.
Matthews
,
C. E.
Wieman
, and
E. A.
Cornell
,
Science
269
,
198
(
1995
).
10.
K.
Davis
,
M.
Mewes
,
M.
Andrews
,
N. v
Druten
,
D.
Durfee
,
D.
Kurn
, and
W.
Ketterle
,
Phys. Rev. Lett.
75
,
3969
(
1995
).
11.
M. R.
Andrews
,
C. G.
Townsend
,
H.-J.
Miesner
,
D. S.
Durfee
,
D. M.
Kurn
, and
W.
Ketterle
,
Science
275
,
637
(
1997
).
12.
E. A.
Burt
,
R. W.
Ghrist
,
C. J.
Myatt
,
M. J.
Holland
,
E. A.
Cornell
, and
C. E.
Wieman
,
Phys. Rev. Lett.
79
,
337
(
1997
).
13.
W.
Ketterle
and
H.-J.
Miesner
,
Phys. Rev. A
56
,
3291
(
1997
).
14.
C.
Pethick
and
H.
Smith
,
Bose-Einstein Condensation in Dilute Gases
(
Cambridge University Press
,
2002
).
15.
L.
Pitaevskii
and
S.
Stringari
,
Bose-Einstein Condensation
, The International Series of Monographs on Physics (
Oxford University Press
,
New York
,
2003
).
16.
A.
Griffin
,
T.
Nikuni
, and
E.
Zaremba
,
Bose-Condensed Gases at Finite Temperatures
(
Cambridge University Press
,
2009
).
17.
P.
Blakie
,
A.
Bradley
,
M.
Davis
,
R.
Ballagh
, and
C.
Gardiner
,
Adv. Phys.
57
,
363
(
2008
).
18.
P. R.
Berman
,
Atom Interferometry
, edited by
P. R.
Berman
(
Academic Press
,
San Diego, CA
,
1997
).
19.
M.
Kasevich
, “
Precision atom interferometry
,” in
Conference on Precision Electromagnetic Measurements (CPEM 2016)
,
2016
.
20.
G. W.
Biedermann
,
H. J.
McGuinness
,
A. V.
Rakholia
,
Y. Y.
Jau
,
D. R.
Wheeler
,
J. D.
Sterk
, and
G. R.
Burns
, “
Atom interferometry in a warm vapor
,”
Phys. Rev. Lett.
118
,
163601
(
2017
).
21.
S.-B.
Lee
,
T.
YongKwon
,
S.
EonPark
,
M.-S.
Heo
, and
H.-G.
Hong
, “
Atomic gravimeter based on atom interferometry being developed at KRISS
,” in
IEEE International Frequency Control Symposium (IFCS)
(
IEEE
,
2018
).
22.
E. R.
Elliott
,
M. C.
Krutzik
,
J. R.
Williams
,
R. J.
Thompson
, and
D. C.
Aveline
,
npj Microgravity
4
,
16
(
2018
).
23.
M. J.
Mazon
,
G. H.
Iyanu
, and
H.
Wang
, “
A portable, compact cold atom physics package for atom interferometry
,” in
Joint Conference of the IEEE International Frequency Control Symposium and European Frequency and Time Forum (EFTF/IFC
) (
IEEE
,
2019
).
24.
V. S.
Malinovsky
,
M. H.
Goerz
,
P. D.
Kunz
, and
M. A.
Kasevich
, “
High-precision atom interferometry using optimal quantum control
,” in
Quantum Information and Measurement (QIM) V: Quantum Technologies
,
2019
.
25.
R.
Kumar
and
A.
Rakonjac
,
Adv. Opt. Technol.
9
,
221
(
2020
).
26.
S.
Abend
,
M.
Gersemann
,
C.
Schubert
,
D.
Schlippert
,
E. M.
Rasel
,
M.
Zimmermann
,
M. A.
Efremov
,
A.
Roura
,
F. A.
Narducci
et al, “
Atom interferometry and its applications
,” arXiv:2001.10976 (
2020
).
27.
J.
Thom
,
C.
Picken
,
J.
Malcolm
,
A.
Kelly
,
X.
Cheng
,
A.
Hinton
,
A.
Bunting
,
G.
Maker
,
N.
Hempler
et al,
Opt. Quantum Sens. Precis. Metrol.
11700
,
1170007
(
2021
).
28.
M. D.
Lachmann
,
H.
Ahlers
,
D.
Becker
,
A. N.
Dinkelaker
,
J.
Grosse
,
O.
Hellmig
,
H.
Müntinga
,
V.
Schkolnik
,
S. T.
Seidel
et al,
Nat. Commun.
12
,
1317
(
2021
).
29.
J.
Li
,
G. R. M. d
Silva
,
W. C.
Huang
,
M.
Fouda
,
J.
Bonacum
,
T.
Kovachy
, and
S. M.
Shahriar
,
Atoms
9
,
51
(
2021
).
30.
B. T.
Seaman
,
M.
Kramer
,
D. Z.
Anderson
, and
M. J.
Holland
,
Phys. Rev. A
75
,
023615
(
2007
).
31.
R. A.
Pepino
,
J.
Cooper
,
D. Z.
Anderson
, and
M. J.
Holland
, “
Atomtronic circuits of diodes and transistors
,” arXiv:0705.3268 (
2007
).
32.
R. A.
Pepino
,
J.
Cooper
,
D.
Meiser
,
D. Z.
Anderson
, and
M. J.
Holland
,
Phys. Rev. A
82
,
013640
(
2010
).
33.
L.
Amico
,
G.
Birkl
,
M.
Boshier
, and
L.-C.
Kwek
,
New J. Phys.
19
,
020201
(
2017
).
34.
L.
Amico
,
M.
Boshier
,
G.
Birkl
,
A.
Minguzzi
,
C.
Miniatura
,
L.-C.
Kwek
,
D.
Aghamalyan
,
V.
Ahufinger
,
D.
Anderson
et al,
AVS Quantum Sci.
3
,
039201
(
2021
).
36.
D. Z.
Anderson
,
Phys. Rev. A
104
,
033311
(
2021
).
37.
L.
Amico
,
D.
Anderson
,
M.
Boshier
,
J.-P.
Brantut
,
L.-C.
Kwek
,
A.
Minguzzi
, and
W.
von Klitzing
,
Rev. Mod. Phys.
94
,
041001
(
2022
).
38.
D.
Bailin
and
A.
Love
,
Introduction to Gauge Field Theory
,
revised ed.
(
Taylor & Francis Group
,
New York
,
1993
).
39.
M. E.
Peskin
and
D. V.
Schroeder
,
An Introduction to Quantum Field Theory
(
CRC Press/Taylor & Francis Group
,
Boca Raton, FL
,
2018
).
40.
L.
Susskind
and
A.
Friedman
,
Special Relatively and Classical Field Theory: The Theoretical Minimum
(
Basic Books
,
New York
,
2017
).
41.
J.
Schwichtenberg
,
No Nonsense Quantum Field Theory
(
No-Nonsense Books
,
2020
).
42.
M. O.
Scully
and
M. S.
Zubairy
,
Quantum Optics
(
Cambridge University Press
,
New York
,
1997
).
43.
R. J.
Glauber
,
Ann. N. Y. Acad. Sci.
480
,
336
(
1986
).
44.
A. A.
Zozulya
and
D. Z.
Anderson
,
Phys. Rev. A
88
,
043641
(
2013
).
45.
G.
Gonzalez
,
Microwave Transistor Amplifiers: Analysis and Design
,
2nd ed
. (
Prentice Hall
,
Upper Saddle River, NJ
,
1997
).
46.
P. M.
Morse
and
H.
Feshback
,
Methods of Theoretical Physics
(
McGraw-Hill Book Company
,
New York
,
1953
), Vol.
1
.
47.
E.
Hecht
,
Optics
(
Pearson
,
Essex England
,
2016
).