Nonminimal Lorentz Violation in Macroscopic Matter

The effects of Lorentz and CPT violations on macroscopic objects are explored. Effective composite coefficients for Lorentz violation are derived in terms of coefficients for electrons, protons, and neutrons in the Standard-Model Extension, including all minimal and nonminimal violations. The hamiltonian and modified Newton's second law for a test body are derived. The framework is applied to free-fall and torsion-balance tests of the weak equivalence principle and to orbital motion. The effects on continuous media are studied, and the frequency shifts in acoustic resonators are calculated.


Introduction
Lorentz invariance is one of the few principles at the heart of both General Relativity (GR) and the Standard Model (SM) of particle physics. However, attempts to reconcile gravity with quantum mechanics suggest this fundamental symmetry of Nature may be broken slightly at low energies [1,2]. While Lorentz violations are expected to be minuscule, simple estimates imply they may be within the reach of high-precision experiments. Spurred by this observation and the development of the Standard-Model Extension (SME) [3][4][5], hundreds of searches for Lorentz violations in a wide variety of systems have been performed in recent decades [6][7][8][9].
The SME is a framework that is designed to characterize all realistic violations of Lorentz and CPT invariance in an effective field theory. It contains both the SM and GR as a Lorentz-invariant limit, which is augmented by all possible terms involving conventional fields. The SME includes terms that violate Lorentz invariance and CPT invariance as well as other fundamental principles, such as diffeomorphism invariance [10,11] and the equivalence principle [12]. A term in the SME action consists of combinations of SM fields, the spacetime metric g µν , and their derivatives contracted with a tensor coefficient for Lorentz violation to form an observer-independent coordinate scalar. The coefficients for Lorentz violation may vary in space and time and could be dynamical in nature. This is especially important when considering Lorentz violations in GR [5,[10][11][12]. However, empirical studies generally assume the coefficients for Lorentz violation are constant in inertial frames, in which case the coefficients impart a nontrivial tensor structure to the vacuum. The dynamics of particles and fields are altered by interactions with this Lorentz-violating background.
A term in the action of the SME is classified, in part, by the mass dimension d of its conventional piece. The restriction to the lowest dimensions d = 3 and d = 4 is call the minimal SME [3][4][5]. The full theory contains an infinite series of terms with d ≥ 3 [13][14][15][16], which when taken together should encompass the low-energy effective limit of any fundamental theory unifying gravity and particle physics. The effects of Lorentz violation typically scale by d-dependent powers of the energy and momentum. As a result, higher-energy particles generally give better sensitivity to nonminimal d > 4 violations. Most tests involving ordinary matter are highly nonrelativistic, reducing their sensitivity to nonminimal violations. However, since the energy is bounded below by the mass of the particle, a subset of effects remain finite in the limit of zero velocity. Nonrelativistic experiments like those discussed below are particularly sensitive to these forms of Lorentz violation.
In this work, we examine the effects of Lorentz violation on macroscopic matter due to microscopic violations in free Dirac fermions. The primary goal is to connect Lorentz violations of arbitrary d in electrons, protons, and neutrons to signals in large objects comprised of these particles.
Ignoring violations that lead to spin-dependent effects, the modified Dirac lagrangian for a fermion of species w is given by [15] The Lorentz violation is controlled by the operatorsâ w µ eff andĉ w µ eff which depend on the four-momentum p µ = i∂ µ . Expanding in p µ , they take the form where a w eff coefficients are nonzero for d = even ≥ 4. They are CPT even and generally produce the same effects in particles and antiparticles.
The above theory yields a modified Dirac equation for electrons, protons, and neutrons, which affects the dynamics of any object made of these particles. The result for ordinary matter is a modified Newton's second law, which depends on macroscopic coefficients c T for the observed test body T. The c T coefficients are linear combinations of the a w The d = 3 and d = 4 violations in (1) are part of the minimal SME [3,4] and have received intense scrutiny in the intervening decades since its construction [9]. It has been shown that the d = 3 violations associated with the a w (3)µ eff coefficients can be removed from the theory through a field redefinition and have no physical effects [3]. We will therefore restrict attention to violations with d ≥ 4. Note, however, that a w (3)µ eff violations are observable through Lorentz-violating matter-gravity couplings [17]. The d = 4 coefficients c w (4)µν eff are observable. They do, however, mimic a species-specific defect in the spacetime metric η µν , which can be removed from one particle through a coordinate transformation [13,18,19]. We use this freedom to eliminate analogous coefficients from the photon sector. Other minimal violations in photons produce birefringence and are strictly limited by astrophysical tests [20][21][22][23][24][25][26][27][28][29][30]. We therefore can safely neglect the effects of minimal Lorentz violations in the pure-photons sector of the SME. We will also neglect violations in electromagnetic interactions [31,32], matter-gravity couplings [17,33,34], and nonminimal violations in photons [13]. Including these would be of interest but would complicate the analysis. They are expected to produce similar effects to those found here.
To date, constraints on minimal c w (4) eff coefficients have been placed in studies involving astrophysics [15,[34][35][36][37][38][39][40][41], tests of the equivalence principle [34,42,43], gravimeters [44], accelerators [45][46][47], electromagnetic cavities [48][49][50], atomic systems [51][52][53][54][55][56][57][58][59][60][61][62], and acoustic resonators [63]. The sensitivities in electrons have reached levels of parts in 10 20 in experiments involving atomic clocks [58] and trapped ions [59,60]. Constraints on protons and neutrons have reached the 10 −29 level in tests using comagnetometers [61,62]. Many of the atomic constraints have been translated into similarly stringent bounds on nonminimal d > 4 violations [64,65]. Tight constraints on nonminimal violations have also been inferred from laboratory [66] and astrophysical [13] tests of relativistic kinematics and from Penning-trap experiments [67]. See [9] for an extensive list of constraints on Lorentz violation in other sectors. While experiments based on microscopic physics give sensitivities that are orders of magnitude beyond what has been demonstrated with macroscopic matter, each experiment is based on different assumptions. So macroscopic tests of Lorentz invariance play an important complementary role. This paper is organized as follows. The basic theory is discussed in Section 2. An effective hamiltonian for ordinary macroscopic matter is constructed in Section 2.1, and a modified Newton's second law is derived in Section 2.2. Section 2.3 provides a brief review of observer rotations of spherical SME coefficients and derives, for the first time, the boosts of the spherical coefficients to first order in boost velocity. Several applications are discussed in Section 3, including tests of the weak equivalence principle in Section 3.1 and tests involving planetary orbits in Section 3.2. A lagrangian for continuous media is given in Section 3.3 and used to derive the frequency shift in piezoelectric acoustic resonators. Section 4 summarizes the results of the work. A useful product identity for spherical-harmonic tensors is derived in the Appendix.

Hamiltonian
Ignoring spin-dependent violations, the leading-order effects of Lorentz violation on a free Dirac fermion of species w are described by the effective hamiltonian h w = E w + δh w , where [15] Here, E w = p 2 w + M 2 w is the conventional free-particle energy for species mass M w . The hamiltonian for antiparticles is given by (3) with the opposite sign onâ w ν eff . Other forms of the hamiltonian (3) may be convenient in practice. A common signal in searches for Lorentz violation is unexpected direction dependence, indicating a violation of rotational symmetry. The prominent role played by rotations in the field motivates the spherical-harmonic expansion where ( m n ) are binomial coefficients, and Y k jm are the orthonormal spherical-harmonic tensors recently derived in Ref. [68]. We use Latin indices a, b, . . . to indicate the restriction to spatial dimensions. In many cases, a nonrelativistic approximation is warranted, leading to a third version, where the nonrelativistic spherical coefficients for Lorentz violation are with k − j = even ≥ 0. All spherical coefficients, including the composite coefficients derived below, obey the complex conjugation relation K * jm = (−1) m K j(−m) . Next, we envision a small but macroscopic volume of matter of mass M containing a large number of electrons, protons, and neutrons. We write the Lorentz-violating change in the hamiltonian for the volume as δH = δH e + δH p + δH n , where δH e , δH p , and δH n represent the total free-particle hamiltonians for electrons, protons, and neutrons, respectively. Each can be written as where N w is the number of particles of species w in the volume, and brackets w indicate the average over all the particles of that species. We then split the bulk motion from the internal motion of the particles. Let p be the total momentum, which is conjugate to the center-of-mass position x. Then p ′ w = p w − M w M p is the conventional momentum of a particle in the center of mass frame. Normally p/M is the velocity of the center of mass and v ′ w = p ′ w /M w is the velocity of a particle relative to the center of mass, but this may no longer be true in the Lorentz-violating case. However, in leading-order calculations, we can assume the usual relations in the Lorentz-violating contributions to the hamiltonian since corrections to the velocity would produce higher-order effects.
The average in (8) can be written as where ⊙ represents the symmetric tensor product, and Y k jm are the rank-k spherical-harmonic tensors [68]. The product p ′⊙q w can be expanded in spherical-harmonic tensors, giving We then make the simplifying assumption that the internal-momentum distribution is approximately isotropic. This implies the j ′ = m ′ = 0 term in the sum dominates when averaged over the particles, yielding the approximation p ′⊙q w w ≈ |p ′ w | q w Y q * 00 / q + 1, which vanishes for odd values of q. Replacing q with 2q, calculation then gives where the sum is limited to integers q for which the arguments of all factorials are nonnegative. The Lorentz-violating hamiltonian for the small macroscopic volume of matter then reads where we define composite coefficients for Lorentz violation The q index sums over all nonnegative integers, the index k is restricted to nonnegative values, and The macroscopic coefficients for Lorentz violation in (13) are the dimensionless combinations of the particle coefficients that affect macroscopic matter at leading order. They are defined so that they depend on the material's particle content and not the size or shape of the object. The macroscopic coefficients are controlled by the mass fractions ρ w /ρ, where ρ w is the mass density for each species and ρ = ∑ w ρ w is the total density. Materials with different particle content have different c T kjm coefficients, so constraints on Lorentz violation in multiple different forms of matter could be combined to separately constrain violations in electrons, protons, and neutrons.
Writing the c T kjm in terms of the relativistic coefficients, we note that both a w kjm coefficients for all dimensions d contribute at each k, so the effects of CPT-even and CPT-odd Lorentz violation cannot be disentangled using only normal matter. We therefore define a single set of effective c T kjm coefficients for matter. These will, however, differ from the macroscopic coefficients for antimatter due to CPT violation. Also note that for fixed d, the c T kjm combinations depend on the internal velocity v ′ Under normal circumstances, the electron internal energy is less than a keV. The internal energies of protons and neutrons is typically on the order of 10 MeV. So while the nonrelativistic internal motions contribute to macroscopic hamiltonian, their effects for fixed d are highly suppressed relative the p ′ w -independent terms with q = 0. The suppressed terms could, however, be used to access different combinations of coefficients. Ignoring the contributions for the internal velocities gives the simplifying approximation Combined with (12), this connects the underlying coefficients for Lorentz violation to the dynamics of macroscopic matter. In ordinary neutral matter made of atoms with atomic number Z and atomic mass M a , The composite coefficients simplify to where we define coefficient combinations for even d, with a similar expression for odd-d a npe kjm coefficients. These coefficient combinations lead to different effects in different types of matter and can be tested in experiments comparing test bodies made of different elements. These include the equivalence-principle experiments discussed in Section 3.1. In contrast, the remaining parts of c T kjm involving the neutron coefficients for Lorentz violation produce identical affects in all matter and are not testable through matter-comparison experiments. Note that the above combinations mirror ones arising in studies of matter-gravity coupling in the SME [34]. A partial match to this work is given in the next section.
Lorentz violation introduces signatures other than composition-dependent dynamics. Searches for these signatures in ordinary matter are less dependent on precise makeup of the mass. In this case, it may suffice to assume roughly equal numbers of neutrons, protons, and electron, which leads to the approximation where with a similar expression for a npe kjm . Note that this approximation also neglects the difference in the proton and neutron masses.
A cartesian version of the macroscopic coefficients for Lorentz violation may be convenient in some applications. Defining the Lorentz-violating hamiltonian becomes Inverting the relation gives The c T k a 1 ...a k tensors are real and totally symmetric. The k index on composite coefficients is restricted to nonnegative integers, and the angular-momentum indices obey j = k, k − 2, k − 4, . . . ≥ 0 and |m| ≤ j.

Equations of motion
The Lorentz-violating hamiltonian δH depends on the center-of-mass momentum p, but is independent of the center-of-mass position x. This implies that the net force, defined as the rate of change in the canonical momentum, is unchanged by Lorentz violation: F = ∂ t p = −∇ x δH. The Lorentz violation enters through the altered relationship between the momentum p and the center-of-mass velocity:ẋ = ∇ p H = p/M + ∇ p δH. Combining the two Hamilton's equations, we arrive at a modified Newton's second law, where Lorentz violation is governed by the symmetric dimensionless tensor We write this in terms of the conventional velocity v = p/M for convenience. Note that v can be taken asẋ in leading-order calculations. We then find a Lorentz-violating force δF ≈ −C · F that depends on the velocity v and the conventional force F. Alternatively, we can write the equations of motion as where the effects of Lorentz violation can be viewed as a velocity-dependent anisotropic mass matrix M(1 + C). This generalizes the modifications from d = 4 violations found in Refs. [69] and [70]. The velocity-dependent C tensor can be expanded in spin-weighted spherical harmonics. This is done by expanding the tensor in the helicity basis [13]: e r =ê r =v = sin θ cos φê x + sin θ sin φê y + cos θê z , The velocity directionv = v/|v| defines the "radial" direction, andê θ andê φ are the usual unit vectors associated with spherical-coordinate angles θ and φ. The components C ab =ê a · C ·ê b in the helicity basis are spin-weighted functions and can be expanded in spin-weighted spherical harmonics s Y jm . The result is All of the objects appearing in these expressions are dimensionless. We note that only composite coefficients with k ≥ 2 affect the macroscopic dynamics at leading-order. For these coefficients, the effects are proportional to |v| k−2 . Since v is the velocity relative to the speed of light, the Lorentz violation from k > 2 will be highly suppressed in most applications. Consequently, the dominant effects are likely those from the k = 2 macroscopic coefficients c T 2 . In the k = 2 restriction, the 8 of 27 tensor C ab = 2 c T 2 ab is constant, and its cartesian components are linear combinations of the spherical coefficients c T 2jm : The resulting velocity-independent effects are then limited to j = 0 isotropic violations and j = 2 quadrupole anisotropies. The k = 2 case can be partially mapped onto previous analyses of composite matter in the SME. In particular, Ref. [34] derives the effects of d = 3 and d = 4 violations in matter and gravity, including the matter-gravity coupling. Dropping the violations involving gravity and those that are cubic in velocity, the coefficients in that work correspond to C ab = −2c Tab −c Ttt δ ab . This provides a map for d = 4 coefficients in the two approaches, which can be extended to higher-d violations through their contributions to c T 2 ab coefficients. This connection could, in principle, be used to convert bounds onc T to bounds on d > 4 violations. Note, however, that this may be problematic in analyses involving boosts of the apparatus since the contributions from different d transform differently. Rotations and boosts of the coefficients are described in the next section.

Lorentz Transformations
Many tests of Lorentz invariance search for changes in a signal with changes in the orientation or velocity of the apparatus. The changes in orientation are typically due to the daily rotation of the Earth, but can be achieved through the use of turntables. The changes in velocity are usually those resulting from the orbital motion of the Earth. Assuming constant coefficients for Lorentz violation in inertial frames, the above motions produce variations in the coefficients in noninertial apparatus-fixed frames. These variations lead to variations in experimental observables, producing potential signals of Lorentz violation. In this section, we review rotations of spherical coefficients and discuss common frames used in tests of Lorentz invariance. We also derive the boosts of spherical coefficients for Lorentz violation to first order in boost velocity.
Rotations of spherical-harmonic expansion coefficients are given by Wigner matrices. Consider the expansion of a spin-weighted function f (n) = ∑ jm f jm s Y jm (n), wheren is a direction unit vector We then consider two frames whose coordinates {x, y, z} and {x ′ , The connection between spherical-harmonic components of f in these two frames is given where D (j) mm ′ (0, β, 0) are Wigner matrices. By convention, tests involving the SME report results in a Sun-centered celestial equatorial inertial reference frame with spacetime coordinates {T, X, Y, Z}. The Z axis points along the Earth's rotation axis, X points towards the vernal equinox, and Y completes the system. The standard time T is define so that T = 0 at vernal equinox in the year 2000. A standard noninertial laboratory frame {t, x, y, z} is defined with x and y horizontal and z pointing vertically up. The x axis points at an angle ϕ east of south. In order to incorporate boosts, we define an intermediate Earth-centered frame with coordinates {T ′ , X ′ , Y ′ , Z ′ } that is boosted but not rotated relative to the Sun frame. The rotation relating the lab frame and the Earth frame is where χ is the colatitude of the laboratory, and α is the right ascension of the laboratory zenith. The transformation of spherical coefficients leads to The right ascension increases at Earth's sidereal rateα ≈ ω ⊕ = 2π/23 hr 56 min, producing a sidereal variation in laboratory-frame coefficients. Horizontal turntables would give variations in ϕ as well.
Ignoring boosts, the Earth-frame and Sun-frame coefficients are the same up to a time-dependent translation. We next determine the relationship between these frames to first order in the boost velocity β. The transformation depends on the tensor structure of the underlying coefficients for Lorentz violation.
Here, we focus on the relativistic spherical a w (d) kjm and c w (d) kjm coefficients in (5). The derivation is the same for both of these sets of coefficients, so we show the calculation for a w (d) kjm only. Inverting relationship (5) between the spherical and cartesian representations of the coefficients for Lorentz violation, we can write Consider a "primed" frame that is moving with small boost velocity β relative to a second "unprimed" frame. The difference between the cartesian coefficients δa w eff in the two frames is then given by Using (32), we can expand the right-hand side of the above expression in terms of spherical coefficients.
Expanding the boost velocity in spherical-harmonic tensors, β a = 4π 3 β ∑ m Y 1m (β)(Y * 1m ) a , and using the product identity derived in the Appendix, we can expand (33) in spherical-harmonic tensors. Combining the result with (5) gives the change in the spherical coefficients for Lorentz violation. The result is the first-order boost of the spherical coefficients where β = |β|,β = β/β, and in terms of the A

The velocity of the Earth in the Sun frame is approximately given by
where the orbital speed is β ≈ 9.9 × 10 −5 , η ≈ 23.4 • is the inclination of the orbit, and Ω ⊕ = 2π/year is the orbital frequency. With this, we can write Combining the boost between the Sun and Earth frame with the rotation between the Earth and laboratory frame, we find the coefficients transform according to The indices m r , m s , and m a are the harmonic numbers for variations at the turntable rotation rateφ, the sidereal rateα, and the annual rate Ω ⊕ , respectively. The c w kjm coefficients in the laboratory and Sun frames obey the same relationship.

Tests of the Equivalence Principle
The modified equations of motion (23) imply that the acceleration of a mass under the influence of gravity depends on its particle content, giving an apparent violation of the weak equivalence principle (WEP). Consequently, tests of the WEP may be adapted for searches for nonminimal Lorentz violations in particle sectors of the SME. Most WEP tests involving macroscopic masses fall into two main categories, free-fall experiments and torsion-pendulum experiments. We consider each of these in turn.
In first class of experiment, the vertical free-fall acceleration of two test masses made of different material are compared. In the experiments performed to date, each mass is released from rest in a vacuum chamber. Accelerations are measured using lasers and corner-cube reflectors attached to the masses [71][72][73][74][75][76].
We work in the standard laboratory frame described in Section 2.3 and assume the gravitational field g = −gê z is approximately uniform. The leading-order change in the acceleration is given by δẍ = 200 + 2a npe (−6 ± 15) × 10 −9 GeV −4 (4 ± 34 ± 34) × 10 −14 GeV −4 (4 ± 10) × 10 −12 GeV −4 Table 1. Limits on isotropic SME coefficients from tests of the equivalence principle. The first column gives the coefficient combinations. The second column contains the combined constraints from ground-based free-fall experiments. The third column lists constraints from the space-based MICROSCOPE experiment. The last column gives the combined constraints from torsion-pendulum experiments.
g C ·ê z , where the C tensor is a function of the unperturbed velocity v = −gtê z . For fixed k ≥ 2, we can integrate twice to get the change in position after time t, δx l = g −1 (−gt) k c T,lab k lz...z , in terms of the laboratory-frame cartesian composite coefficients for Lorentz violation. In principle, free-fall experiments could search for ∼ t k displacements from the velocity-dependent k > 2 violations. These effects, however, are highly suppressed by the small velocities. For the velocity-independent k = 2 cases, the vertical acceleration is constant and given by This includes the isotropic c T,lab 200 and the quadrupole c T,lab 220 . The quadrupole violations lead to sidereal variations, which could be sought in future analyses. However, in order to understand the reach of these experiments, we will ignore boosts and focus on isotropic violations c T,lab 200 ≈ c T 200 . In this particular limit, the acceleration is proportional to the gravitational field,ẍ = (1 − 1 π c T 200 )g. So isotropic Lorentz violation mimics a difference in the inertial mass M I and gravitational mass M g . This type of behavior is traditionally characterized using the Eötvös parameter, defined as the difference in free-fall acceleration ∆a divided by the averageā for the two test masses. In terms of the masses, this gives η = ∆a/ā = ∆r/r, where r = M g /M I . In the present context, we find get an effective Lorentz-violating Eötvös parameter for fixed dimension d. The differences in the test bodies enter through the difference ∆(Z/M a ) in the ratio of the atomic number and atom mass. Note that the simple correspondence (39)  . Translating these to constraints on the above combinations of SME coefficients, we find (3 ± 13) × 10 −9 using copper and uranium [72], (21 ± 60) × 10 −9 using aluminum and beryllium [74], (−10 ± 58) × 10 −9 using aluminum and copper [74], (17 ± 131) × 10 −9 using aluminum and carbon [74], (21 ± 53 ± 74) × 10 −9 using aluminum and copper [76], and (23 ± 30 ± 33) × 10 −9 copper and tungsten [76]. We combine these to produce a "best-fit" ground-based measurement and then translate this to constraints on the isotropic-coefficient combinations. The results up to d = 8 are given the second column Table 1. [77][78][79] could also be used to search for Lorentz violation. For example, the T-SAGE instrument aboard the MICROSCOPE satellite has placed a constraint of η = (−1 ± 9 ± 9) × 10 −15 on the difference between the free-fall accelerations of titanium and platinum [79]. The resulting constraint on isotropic Lorentz violations are included in the third column of Table 1, demonstrating that sensitivities to Lorentz violation on the order of 10 −13 GeV 4−d are possible.

Similar tests of WEP in space
The classic Eötvös experiment [80] and its descendants represent another class of WEP tests based on torsion pendulums. In the prototypical experiment, two tests masses of different composition are attached to the ends of a rod hanging from a horizontal fiber. A difference in the gravitational and inertial masses would lead to an imbalance in the horizontal components of the gravitational and centrifugal forces, leading to a net torque about the fiber. Modern versions achieve high sensitive by seeking modulated signals due to the changing field from the Sun over the day [81,82] or by rotating the apparatus in the laboratory [83][84][85].
Anisotropic violations could also be tested in these experiments. The daily rotation of the laboratory in ground-based experiments will lead to variations at multiples of the sidereal frequency. Signals in space-based experiments will arise at harmonics of the satellite rotation rate. Boosts will introduce additional frequencies in the variations of the signal, including the annual frequency due to the motion of the Earth around the Sun. A search for these types of variations in MICROSCOPE data was recently carried out [86], where d = 3 and d = 4 violations in the matter-gravity couplings of the SME are constrained down to the expected 10 −13 GeV 4−d range. A similar analysis could be used to constrain higher-order Lorentz violation.
Future space-based tests include STEP [77] and GG [78], which could reach sensitivities two or three orders of magnitude better than MICROSCOPE. Other promising opportunities for future studies include experiments utilizing drop towers [87], balloons [88], bouncing masses [89], and sounding rockets [90].

Orbits
In this section, we consider the effects of Lorentz violation on a satellite in a gravitational orbit around a larger body. For simplicity, we'll restrict attention to approximately circular orbits and neglect effects depending on both eccentricity and coefficients for Lorentz violation. We work in a fixed orbit-centered frame with the z axis along the orbital axis and z = 0 in the orbital plane of the Lorentz-invariant limit. We also use an orbit time t defined so that the satellite velocity is along the x direction at t = 0. The position of the satellite can be described using cylindrical coordinates {ρ, ϕ, z} in the orbit frame. We denote the corresponding basis unit vectors as {ê ρ ,ê ϕ ,ê z }.
At leading order, the acceleration of the satellite is given byẍ = g − C · g, where g = −R 3 ω 2 |x| −3 x is the gravitational field in terms of the usual semimajor axis R and the orbit angular frequency ω. This leads to the change in position δx of the satellite due to Lorentz violation that satisfies δẍ = −ω 2 δx + 3ω 2 (ê ρ · δx)ê ρ + Rω 2 C ·ê ρ , neglecting terms involving the small eccentricity. The above implies the cylindrical-basis components obey the coupled differential equations The components C ab =ê a · C ·ê b can be taken as functions of the conventional velocity v ≈ Rωê ϕ . Since the velocity v is periodic with period 2π/ω, the C ab components in (40) drive changes in the motion at harmonics of ω. We characterize this driving force using the form ab driving amplitudes arise naturally out of the spherical-harmonic expansion of the C ab components. Recall that the spherical-harmonic expansion of C uses the helicity vectors (25) defined with respect to the velocity v. Matching to the cylindrical basis givesê r =ê ϕ andê in the orbital plane. The velocity vector points at polar angle θ = π/2 and azimuthal angle φ = ωt. This leads to the driving amplitudes in terms of the orbit-frame coefficients for Lorentz violation c T,orb kjm . These depend on spherical harmonics for theê x direction, which lies at θ = π/2 and φ = 0.
The orbit-specific C

[m]
ab coefficient combinations determine the leading-order effects of Lorentz violation for a particular satellite. However, they must be connected to the standard Sun-frame coefficients to be useful. Ignoring boosts, the Sun-frame coefficients c T kjm and the orbit-frame coefficients c T,orb kjm are related through the rotation where α, γ, and η are a convenient set of Euler angles and are illustrated in Figure 1. The angle α is between the Sun-frame X axis and the line of nodes, η is the inclination of the orbit relative to the X-Y plane and γ is the angle between the orbit-frame x axis and the line of nodes. For our analysis, γ is somewhat arbitrary. It can be chosen, for example, so that the pericenter lies on theê x axis. The α and η angles for the eight planets are given in Table 2.  To solve (40), we first seek force-free homogeneous solutions. Several homogeneous solutions exist that can be connected to conventional perturbations to circular orbits. The general homogeneous solution is The constant δϕ represents a small translation along the orbit. The constant δR gives a transition to a circular orbit that is larger or smaller by δR. To see this note that while R and ω may change, the combination R 3 ω 2 is constant for Kepler orbits, so δ(R 3 ω 2 ) = 3R 2 ω 2 δR + 2R 3 ωδω = 0. Using this, we find that the change in position due to a change in radius is δx = δRê ρ + Rδω tê ϕ = δRê ρ − 3 2 ωtδRê ϕ , matching the above result. This variation is only valid for sufficiently small times. The constant ε adds small eccentricity with a pericenter at time t 1 . It can be characterized using an eccentricity vector ε = Re iεe −iωt 1 (ê x + iê y ) , which has magnitude ε and points to the pericenter. The perturbations in the position can be written as δx ρ = −Rε ·ê ρ , δx ϕ = −2Rε ·ê ϕ . Setting t 1 to one quarter of the orbit period places the pericenter on the x axis. The angle η ′ corresponds to a small inclination from the orbit-frame x-y plane with the ascending node at time t 2 . Using these results, we can distinguish between conventional perturbations of the orbit and those caused by Lorentz violation.
For the Lorentz-violating inhomogeneous problem, we first consider the m = 0 and m = 1 special cases separately since these harmonics also appear in the homogeneous solutions. For m = 0, we find that C ρρ . This gives a change in the size of the orbit without the corresponding change in the frequency required by Kepler's third law. The effect mimics a small change in either Newton's constant G or the source mass. While it may be difficult to detect, this form of Lorentz violation could be sought, in principle, by comparing the third-law constant R 3 ω 2 of different satellites orbiting the same source.
For the m = ±1 case, two distinct effects arise. We first note that C [1] ρρ = iC [1] ϕρ and neither of C [1] ρρ or C [1] ϕρ contain velocity-independent k = 2 contributions. So only C [1] zρ gives effects that are unsuppressed by small orbital velocities. Nonetheless, the C [1] ρρ = iC [1] ϕρ matrix element drives a speed-dependent change in the eccentricity. The inhomogeneous solution can be taken as where we parameterize the phase as C [1] ρρ = i C [1] ρρ e −iωt 3 . The above is only valid for small t, but the result is an eccentricity that increases a rate ofε = ω C [1] ρρ with pericenter at time t 3 . The rate of change in the eccentricity vector isε = ω Re C [1] ρρ (ê x + iê y ) . (45) Note that this will add to the conventional eccentricity vector that may point in a different direction. The small eccentricities of the planets lead to crude limits on the above effect. As an example, consider the Earth. The eccentricity added per orbit is 2π C [1] ρρ , and the Earth has made roughly N ≈ 4.5 × 10 9 orbits over the age of the Solar System. Taking the Earth's eccentricity ε ≈ 0.017 as an upper bound on the eccentricity due to Lorentz violation, we find the constraint C Since k = 2 Lorentz violations do not contribute, the dominant effects would likely be from k = 3, which are linear in speed. Earth's speed is about v ≈ 10 −4 , implying potential sensitivity on the order of 10 −9 to c T 3jm coefficients. Venus, with its smaller eccentricity and shorter year, is the only planet yielding a slightly better sensitivity.
The C [1] zρ coefficient combination gives the modification where we parameterize C [1] zρ = − C [1] zρ e −iωt 4 . This implies a rotation of the orbital plane at an instantaneous rate of Ω = ω C [1] zρ about an in-plane axis pointing towards the satellite at time t 4 . We can account for both the rate and the direction by defining a rotation vector The orbit axes will gradually rotate according toė a = Ω ×ê a . The resulting secular variations in the Euler angles areα = ω sin γ Im C [1] zρ + cos γ Re C Note that the C [1] zρ coefficient combination also depends on the Euler angles. A demonstration of the above rotation is provided in Ref. [70], where the effects of d = 4 violations on binary systems are simulated.
Assuming a small net effect, we can approximate the total rotation of the orbital plane after N orbits as 2πN C [1] zρ . Given that the orbital planes of the eight planets and the Sun's equator all differ by no more than about 10 • , we take this as an upper bound on the change in inclination for each planet over the age of the Solar System. Mercury, with its short year, then produces the tightest constraint, C Planetary ephemerides could be used to place more rigorous bounds at levels comparable to the simple estimate given above. For example, limits on the evolution of planetary orbits have been shown to constrain the dimensionlesss µν coefficients in the pure-gravity sector of the SME down to parts in 10 12 [92,93]. Lunar laser ranging has also been used to test Lorentz symmetry in gravity [94,95] and in matter-gravity couplings [96] down to parts in 10 12 . Thes µν coefficients produce effects similar to those found above [12], and we expect similar constraints on c T 22m coefficients. This implies sensitivities at the level of ∼ 10 −12 GeV 4−d to a ′ npe (d) kjm SME coefficient combinations. Binary pulsars provide another test of Lorentz invariance in orbital dynamics [12,97]. These systems have been used to test Lorentz invariance to parts in 10 11 in gravity [98,99] and matter-gravity couplings [40]. They've also been used to search for velocity-dependent effects from dimension d = 5 terms [100] and from d = 8 cubic terms in the gravity sector of the SME [101]. We also note that Lorentz violation can be constrained with non-binary pulsars [41,102]. Binary pulsar are unique among orbital tests in that they provide clean access to neutron coefficients for Lorentz violation and are therefore complementary to tests involving ordinary matter.
The perturbations driven at higher frequencies with m > 1 are ϕρ e imωt , Unlike the m = 1 case, which gave a secular evolution of the orbit, violations with m ≥ 2 produce periodic deviations from the conventional orbit. For example, δx z produces periodic oscillations about of the average orbital plane. Among the effects from δx ρ and δx ϕ displacements is a time-dependent change in the areal velocity δȦ = Rωδx ρ + 1 ϕρ e imωt , violating Kepler's second law. Again, these periodic variations are similar to ones arising in the gravity sector of the SME [12] and could be sought in planetary motion or in lunar laser ranging. Note that since k ≥ j ≥ |m|, the effects of Lorentz violation at frequencies greater than 2ω necessarily involve the orbital speed of the satellite. Consequently, unsuppressed speed-independent periodic variations only arise at twice the orbital frequency.

Acoustic Resonators
This section considers the effects of Lorentz violation in continuous media with particular focus on acoustic resonances in piezoelectric materials. For continuous media, the Lorentz-violating hamiltonian density can be taken as where ρ is the mass density of the material, P is the momentum density, andP = P /|P |. Using the above, we can find modifications to the Hamilton equations of motion. Alternatively, we could instead employ a lagrangian approach, where the leading-order change to the Lagrange density can be taken as where v = v(x, t) is the local velocity of the medium, andv = v/|v|. Denoting the mechanical displacement of the medium at equilibrium position x in a body-fixed frame as u = u(x, t), the velocity is then given by v =u. Note that v ≈ P /ρ, but differs slightly from the usual result due to Lorentz violation. The conventional Lagrange density for a piezoelectric material is given by where u a ,b = ∂u a /∂x b are spatial derivatives of the displacements u a , φ ,a = ∂φ/∂x a is the gradient of the electric potential φ, S abcd is the stiffness tensor, ǫ ab is the permittivity tensor, and e abc is the piezoelectric tensor. The stiffness tensor S abcd is taken to be symmetric in the first pair of indices and the last pair of indices and symmetric under interchange of the pairs, giving twenty-one independent components. The permittivity tensor ǫ ab is symmetric, and the piezoelectric tensor e abc is symmetric in the last two indices. The equations of motion for the system including Lorentz violation are given by where T ab = S abcd u c ,d + e cab φ ,c is the stress tensor, D a = −ǫ ab φ ,b + e abc u b ,c is the electric displacement field, and C ab is the Lorentz-violating tensor from (24) evaluated at velocity v =u.
Periodic solutions to (53) can be found using methods similar to those used for orbits. Solutions with period 2π/ω will in general includes various harmonics of the fundamental frequency ω. To find them, first expand each variable in Fourier modes: The equations of motion (53) lead to a set of coupled equations relating the various Fourier components, which can be solved perturbatively. However, we are primarily interested in changes to the frequency ω, which can be found using a simpler method.
We begin by assuming the solution with Lorentz violation u has frequency and amplitudes that are close to those for a conventional solution u 0 . Manipulating the equations of motion for u and u 0 , one can show the relation where the left-hand side is integrated over the volume V of the resonator, and the right-hand side is integrated over the surface ∂V. The conventional stress tensor T 0 depends on u 0 , the conventional potential φ 0 , and the conventional displacement field D 0 . We then assume the surface terms vanish, giving Assuming simple harmonic conventional solutions, this expression oscillates at frequencies mω ± ω 0 = (m ± 1)ω 0 + mδω, where δω = ω − ω 0 is the shift in the fundamental frequency from the usual frequency ω 0 . Writing the amplitudes as u [m] = u [m] 0 + δu [m] , where δu [m] is the change due to Lorentz violation, we can expand the frequency components of (55) in small parameters depending on coefficients for Lorentz violation. The zeroth-order equations are identically satisfied. The first-order equations give Note that u [m] 0 = 0 for m = ±1 since we assume u 0 is simple harmonic. The shift in frequency δω can be isolated by taking m = ∓1, which gives where brackets t indicate the time average. The C tensor in this expression is calculated using the conventional velocity v 0 =u 0 . The leading-order frequency shift is then completely determined by the coefficients for Lorentz violation and the usual solution u 0 . The time averages in (57) may be difficult to calculate in general but are relatively simple in the case of standing waves with local linear polarization, where we can take u 0 (x, t) → u 0 (x) sin(ω 0 t). The velocity is replaced with v 0 (x, t) → ω 0 u(x) cos(ω 0 t) and is parallel to the displacement u 0 (x). The time average in the denominator of (57) becomes u 0 · u 0 t → 1 2 u 2 0 . The numerator can be shown to vanish for odd values of k. For fixed even values of k, the time average in the numerator becomes The frequency shift is then given by where k is restricted to even values k ≥ 2, and c T,lab kjm are laboratory-frame coefficients. The dimensionless factors multiplying the c T kjm coefficients determine the sensitivity of an acoustic-resonator experiment. Assuming oscillation amplitudes on the order of 100 angstrom and frequencies on the order of a MHz, these factors scale as ∼ 10 −10(k−2) . This drastically reduces the sensitivity to violations with k = 2. We therefore focus on the k = 2 case. The problem simplifies even further for cases in which the vibration directionû 0 (x) is relatively constant over the volume of the resonator: assuming in the last line that the medium is comprised of roughly equal numbers of neutrons, protons, and electrons. The frequency shift is then limited to quadrupole j = 2 and isotropic j = 0 violations. The rotational and orbital motion of the Earth implies that the laboratory frame is noninertial. As a result, the laboratory-frame coefficients will change as the orientation and velocity of the laboratory changes, producing periodic variations in the frequency shift. We account for these changes using the transformation between the laboratory frame and Sun-centered frame discussed in Section 2.3. The rotations introduce sidereal variations in the laboratory-frame coefficients. The coefficients also vary with the angle of the resonator in the laboratory. In experiments involving rotating turntables, this produces variations at the turn rate. Annual changes in the velocity of Earth lead to annual variations in the signal. These, however, enter through boosts and are suppressed by the boost velocity β ≈ 10 −4 relative to the other variations.
The fluctuations in the frequency shift take the form where ω ⊕ = 2π/23hr 56min and Ω ⊕ = 2π/1yr are respectively the sidereal and annual frequencies.
The time T ⊕ is defined so that the laboratory zenith points at right ascension α = 0 when T ⊕ = 0, and time T = 0 at the vernal equinox. The angle ϕ is between the laboratory-frame x axis and south. The laboratory frame is fixed to the resonator, which may be affixed to a turntable. In this case, ϕ changes at the turntable rotation rate ω r . The indices m r , m s , and m a are the harmonic numbers for variations at respectively the turntable rotation frequency ω r , sidereal frequency ω ⊕ , and annual frequency Ω ⊕ . The amplitudes obey the relation A * m r m s m a = A (−m r )(−m s )(−m a ) , ensuring that the frequency shift is real. Applying the Lorentz transformations outlined in Section 2.3 to the laboratory-frame coefficients in (60), we find that the modulation amplitudes due to rotations only are in terms of the Sun-frame a ′ npe (d) kjm coefficients. The isotropic j = 0 violations produce a constant shift. The quadrupole j = 2 violations give variations at frequencies m r ω r + m s ω ⊕ up to the second harmonic in both the turntable and sidereal rates: |m r |, |m s | ≤ 2. Note, however, that |m r | = 1 variations will be absent in oscillators with horizontal or vertical vibrations.
Including leading-order boost effects due to the orbital motion of the Earth gives variations at frequencies m r ω r + m s ω ⊕ + m a Ω ⊕ with m a = ±1. The amplitudes for these are given by where the numerical Γ (d)k ′ j ′ m ′ m ′′ (2−2l)jm s constants are given in (35), and the B m ′′ m a boost factors are in (36). This gives sensitivity to other coefficients for Lorentz violation but at levels suppressed by the small boost velocity β ≈ 10 −4 of the Earth.
Searches for Lorentz violation in quartz resonators have demonstrated sensitivities on the order of parts in 10 14 to d = 4 violations [63] and are expected to improve by two orders of magnitude [103]. We therefore expect sensitivities near 10 −16 GeV 4−d to the dimension-d combinations a ′

Summary
A violation of Lorentz invariance would necessarily indicate new physics with potential origins in quantum gravity. High-precision experiments have limited violations in a large variety of systems [9]. In this paper, we derive the effects of Lorentz violation on dynamics of ordinary matter. We include all linear dimension-d violations in the electrons, protons, and neutrons, excluding violations involving electromagnetic and gravitational interactions.
The effective hamiltonian for a macroscopic test body is derived in Section 2.1. The Lorentz-violating contributions are given in (12) in terms of macroscopic coefficients for Lorentz violation c T kjm . Equation (14) relates these coefficients to underlying SME coefficients for electrons, protons, and neutrons. Ignoring internal kinetic energy, the result reduces to (15). Equation (16) gives c T kjm for matter with equal numbers electrons and protons, and (18) is for matter with equal numbers of electrons, protons, and neutrons. The equations of motion are discussed in Section 2.2. A modified Newton's second law is given in (23). Section 2.3 discusses observer Lorentz transformations of the coefficients, relating coefficients in the Sun-centered equatorial frame to a standard laboratory frame. The boosts are calculated to first order in velocity, resulting in (37). Section 3 contains several applications. Tests of the the weak equivalence principle are discussed in Section 3.1, including tests involving free-fall experiments [72][73][74][75][76], the space-based MICROSCOPE experiment [79] and torsion-balance experiments [81][82][83][84][85]. Implied bounds on isotropic Lorentz violation from these experiments are given in Table 1, demonstrating sensitivities down to ∼ 10 −13 GeV 4−d to dimension d violations.
Planetary orbits are discussed in Section 3.2. The effects of Lorentz violation include a drift in eccentricity, a rotation of the orbital plane, and periodic variations about conventional orbits. The small eccentricities of Earth and Venus limit Lorentz violation at the ∼ 10 −9 GeV 4−d level. The approximate alignment of the planets' orbital planes leads to bounds of ∼ 10 −12 GeV 4−d . Improvements on these rough constraints are expected in detailed studies of planetary ephemerides [92,93] and through lunar laser ranging [94][95][96]. Binary pulsars offer another promising area of study that is particularly sensitive to Lorentz violations in neutrons [12,40,[97][98][99][100][101]. Section 3.3 gives the Lorentz-violating Lagrange density for continuous media (51). The shift in resonant frequency in piezoelectric acoustic resonators is calculated, including boost effects. The shifts vary periodically at frequencies involving the the turntable rotation rate, the Earth's sidereal rotation rate, and the Earth's orbital frequency. Experiments have demonstrated sensitivities at parts in 10 14 to dimension d = 4 Lorentz violations [63] and are expected to reach 10 −16 GeV 4−d to arbitrary dimension d violations [103].
These results show that extreme precision can be achieve in studies of spacetime symmetries in macroscopic matter. While not as sensitive as the best of the microscopic tests [58][59][60][61][62], experiments involving ordinary matter rely on different assumptions and may provide access to different combinations of SME coefficients and therefore represent a powerful tool in our search for new physics.

Conflicts of Interest:
The author declares no conflicts of interest.

Appendix A
This Appendix derives the symmetric product identity for spherical-harmonic tensors. See Ref. [68] for a detailed discussion of the Y ̺ jm tensors and the notation used here. Expanding the symmetric product of two spherical-harmonic tensors in the basis of spherical-harmonic tensors, we can write (A1) The A ̺ 1 ̺ 2 j 1 m 1 j 2 m 2 J M coefficients are nonzero for the usual angular-momentum-addition relations M = m 1 + m 2 and j 1 + j 2 ≥ J ≥ |j 1 − j 2 | and for j 1 + j 2 − J = even. The nonzero values are real and given by for convenience. Our derivation starts by considering the vector for arbitrary complex number ζ. The unit vectorsê ↑ = (ê x + iê y )/ √ 2,ê ↓ = (ê x − iê y )/ √ 2, andê z form a spin-eigenbasis for quantization along the z axis. A short calculation reveals that the j-fold symmetric product of ξ is where Y jm = Y for |m| ≤ j. The symmetric product in (A5) serves as a generating function for the traceless spherical-harmonic tensors. Note that ξ · ξ = 0, which confirming that it is traceless.
Next consider the inner product ξ ⊙J · (Y * j 1 m 1 ⊙ Y * j 2 m 2 ) = (ξ ⊙j 1 · Y * j 1 m 1 )(ξ ⊙j 2 · Y * j 2 m 2 ) , where J = j 1 + j 2 . The two sides of this equation evaluate to which implies The complex conjugate of this gives the A j 1 j 2 j 1 m 1 j 2 m 2 J M coefficients for traceless tensors. To find the inner product for tensors of nonzero trace, we consider traces of the product Using ξ 1 · ξ 2 = −(ζ 1 − ζ 2 ) 2 /2ζ 1 ζ 2 , one can show that taking N traces gives where g is the euclidean metric, and x n indicates the falling factorial. Matching terms by their powers in ξ 1 and ξ 2 , we find where can be used to show that