Massive gravity

From Wikipedia, the free encyclopedia
Jump to navigation Jump to search

In theoretical physics, massive gravity is a theory of gravity that modifies general relativity by endowing the graviton with a nonzero mass. In the classical theory, this means that gravitational waves obey a massive wave equation and hence travel at speeds below the speed of light.

Massive gravity has a long and winding history, dating back to the 1930s when Wolfgang Pauli and Markus Fierz first developed a theory of a massive spin-2 field propagating on a flat spacetime background. It was later realized in the 1970s that theories of a massive graviton suffered from dangerous pathologies, including a ghost mode and a discontinuity with general relativity in the limit where the graviton mass goes to zero. While solutions to these problems had existed for some time in three spacetime dimensions,[1][2] they were not solved in four dimensions and higher until the work of Claudia de Rham, Gregory Gabadadze, and Andrew Tolley (dRGT model) in 2010.

One of the very early massive gravity theories was constructed in 1965 by Ogievetsky and Polubarinov (OP).[3] Despite OP model coincides with the ghost-free massive gravity models rediscovered in dRGT, OP model has been almost unknown among contemporary physicists who work on massive gravity, perhaps because the strategy followed in that model was quite different from what is generally adopted at present.[4] Massive dual gravity to OP model[5] can be obtained by coupling the dual graviton field to the curl of its own energy-momentum tensor.[6][7] The mathematical underpinning of the dual model in 4-D is very similar to that of Galileon theory.[8][9] This is manifested in the dual formulation of Galileon theory.[10][11]

The fact that general relativity is modified at large distances in massive gravity provides a possible explanation for the accelerated expansion of the Universe that does not require any dark energy. Massive gravity and its extensions, such as bimetric gravity,[12] can yield cosmological solutions which do in fact display late-time acceleration in agreement with observations.[13][14][15]

Observations of gravitational waves have constrained the Compton wavelength of the graviton to be λg > 1.6×1016 m, which can be interpreted as a bound on the graviton mass mg < 7.7×10−23 eV/c2.[16] The most recent observations pose a new constraint on λg > 1.66×1016 m for mg < 7.45×10−23 eV/c2 and λg > 1.83×1016 m for mg < 6.76×10−23 eV/c2.[17]

Linearized massive gravity[edit]

At the linear level, one can construct a theory of a massive spin-2 field propagating on Minkowski space. This can be seen as an extension of linearized gravity in the following way. Linearized gravity is obtained by linearizing general relativity around flat space, , where is the Planck mass with the gravitational constant. This leads to a kinetic term in the Lagrangian for which is consistent with diffeomorphism invariance, as well as a coupling to matter of the form


where is the stress–energy tensor. This kinetic term and matter coupling combined are nothing other than the Einstein-Hilbert action linearized about flat space.

Massive gravity is obtained by adding nonderivative interaction terms for . At the linear level (i.e., second order in ), there are only two possible mass terms:

Fierz and Pauli[18] showed in 1939 that this only propagates the expected five polarizations of a massive graviton (as compared to two for the massless case) if the coefficients are chosen so that . Any other choice will unlock a sixth, ghostly degree of freedom. A ghost is a mode with a negative kinetic energy. Its Hamiltonian is unbounded from below and it is therefore unstable to decay into particles of arbitrarily large positive and negative energies. The Fierz-Pauli mass term,

is therefore the unique consistent linear theory of a massive spin-2 field.

The vDVZ discontinuity[edit]

In the 1970s Hendrik van Dam and Martinus J. G. Veltman[19] and, independently, Valentin I. Zakharov[20] discovered a peculiar property of Fierz-Pauli massive gravity: its predictions do not uniformly reduce to those of general relativity in the limit . In particular, while at small scales (shorter than the Compton wavelength of the graviton mass), Newton's gravitational law is recovered, the bending of light is only three quarters of the result Albert Einstein obtained in general relativity. This is known as the vDVZ discontinuity.

We may understand the smaller light bending as follows. The Fierz-Pauli massive graviton, due to the broken diffeomorphism invariance, propagates three extra degrees of freedom compared to the massless graviton of linearized general relativity. These three degrees of freedom package themselves into a vector field, which is irrelevant for our purposes, and a scalar field. This scalar mode exerts an extra attraction in the massive case compared to the massless case. Hence, if one wants measurements of the force exerted between nonrelativistic masses to agree, the coupling constant of the massive theory should be smaller than that of the massless theory. But light bending is blind to the scalar sector, because the stress-energy tensor of light is traceless. Hence, provided the two theories agree on the force between nonrelativistic probes, the massive theory would predict a smaller light bending than the massless one.

Vainshtein screening[edit]

It was argued by Vainshtein[21] two years later that the vDVZ discontinuity is an artifact of the linear theory, and that the predictions of general relativity are in fact recovered at small scales when one takes into account nonlinear effects, i.e., higher than quadratic terms in . Heuristically speaking, within a region known as the Vainshtein radius, fluctuations of the scalar mode become nonlinear, and its higher-order derivative terms become larger than the canonical kinetic term. Canonically normalizing the scalar around this background therefore leads to a heavily suppressed kinetic term, which damps fluctuations of the scalar within the Vainshtein radius. Because the extra force mediated by the scalar is proportional to (minus) its gradient, this leads to a much smaller extra force than we would have calculated just using the linear Fierz-Pauli theory.

This phenomenon, known as Vainshtein screening, is at play not just in massive gravity, but also in related theories of modified gravity such as DGP and certain scalar-tensor theories, where it is crucial for hiding the effects of modified gravity in the solar system. This allows these theories to match terrestrial and solar-system tests of gravity as well as general relativity does, while maintaining large deviations at larger distances. In this way these theories can lead to cosmic acceleration and have observable imprints on the large-scale structure of the Universe without running afoul of other, much more stringent constraints from observations closer to home.

The Boulware-Deser ghost[edit]

As a response to Freund-Maheshwari-Schonberg finite-range gravity model,[22] and around the same time as the vDVZ discontinuity and Vainshtein mechanism were discovered, David Boulware and Stanley Deser found in 1972 that generic nonlinear extensions of the Fierz-Pauli theory reintroduced the dangerous ghost mode;[23] the tuning which ensured this mode's absence at quadratic order was, they found, generally broken at cubic and higher orders, reintroducing the ghost at those orders. As a result, this Boulware-Deser ghost would be present around, for example, highly inhomogeneous backgrounds.

This is problematic because a linearized theory of gravity, like Fierz-Pauli, is well-defined on its own but cannot interact with matter, as the coupling breaks diffeomorphism invariance. This must be remedied by adding new terms at higher and higher orders, ad infinitum. For a massless graviton, this process converges and the end result is well-known: one simply arrives at general relativity. This is the meaning of the statement that general relativity is the unique theory (up to conditions on dimensionality, locality, etc.) of a massless spin-2 field.

In order for massive gravity to actually describe gravity, i.e., a massive spin-2 field coupling to matter and thereby mediating the gravitational force, a nonlinear completion must similarly be obtained. The Boulware-Deser ghost presents a serious obstacle to such an endeavor. The vast majority of theories of massive and interacting spin-2 fields will suffer from this ghost and therefore not be viable. In fact, until 2010 it was widely believed that all Lorentz-invariant massive gravity theories possessed the Boulware-Deser ghost.[24]

Ghost-free massive gravity[edit]

In 2010 a breakthrough was achieved when de Rham, Gabadadze, and Tolley constructed, order by order, a theory of massive gravity with coefficients tuned to avoid the Boulware-Deser ghost by packaging all ghostly (i.e., higher-derivative) operators into total derivatives which do not contribute to the equations of motion.[25][26] The complete absence of the Boulware-Deser ghost, to all orders and beyond the decoupling limit, was subsequently proven by Fawad Hassan and Rachel Rosen.[27][28]

The action for the ghost-free de Rham-Gabadadze-Tolley (dRGT) massive gravity is given by[29]

or, equivalently,

The ingredients require some explanation. As in standard general relativity, there is an Einstein-Hilbert kinetic term proportional to the Ricci scalar and a minimal coupling to the matter Lagrangian , with representing all of the matter fields, such as those of the Standard Model. The new piece is a mass term, or interaction potential, constructed carefully to avoid the Boulware-Deser ghost, with an interaction strength which is (if the nonzero are ) closely related to the mass of the graviton.

The interaction potential is built out of the elementary symmetric polynomials of the eigenvalues of the matrices or , parametrized by dimensionless coupling constants or , respectively. Here is the matrix square root of the matrix . Written in index notation, is defined by the relation

We have introduced a reference metric in order to construct the interaction term. There is a simple reason for this: it is impossible to construct a nontrivial interaction (i.e., nonderivative) term from alone. The only possibilities are and , both of which lead to a cosmological constant term rather than a bona fide interaction. Physically, corresponds to the background metric around which fluctuations take the Fierz-Pauli form. This means that, for instance, nonlinearly completing the Fierz-Pauli theory around Minkowski space given above will lead to dRGT massive gravity with , although the proof of absence of the Boulware-Deser ghost holds for general .[30]

In principle, the reference metric must be specified by hand, and therefore there is no single dRGT massive gravity theory, as the theory with a flat reference metric is different from one with a de Sitter reference metric, etc. Alternatively, one can think of as a constant of the theory, much like or . Instead of specifying a reference metric from the start, one can allow it to have its own dynamics. If the kinetic term for is also Einstein-Hilbert, then the theory remains ghost-free and we are left with a theory of massive bigravity,[12] propagating the two degrees of freedom of a massless graviton in addition to the five of a massive one.

In practice it is unnecessary to compute the eigenvalues of (or ) in order to obtain the . They can be written directly in terms of as

where brackets indicate a trace, . It is the particular antisymmetric combination of terms in each of the which is responsible for rendering the Boulware-Deser ghost nondynamical.

The choice to use or , with the identity matrix, is a convention, as in both cases the ghost-free mass term is a linear combination of the elementary symmetric polynomials of the chosen matrix. One can transform from one basis to the other, in which case the coefficients satisfy the relationship[29]

Massive gravity in the vielbein language[edit]

The presence of a square-root matrix is somewhat awkward and points to an alternative, simpler formulation in terms of vielbeins. Splitting the metrics into vielbeins as


and then defining one-forms

the ghost-free interaction terms above can be written simply as (up to numerical factors)[31]

In terms of vielbeins, rather than metrics, we can therefore see the physical significance of the ghost-free dRGT potential terms quite clearly: they are simply all the different possible combinations of wedge products of the vielbeins of the two metrics.

Note that massive gravity in the metric and vielbein formulations are only equivalent if the symmetry condition

is satisfied. While this is true for most physical situations, there may be cases, such as when matter couples to both metrics or in multimetric theories with interaction cycles, in which it is not. In these cases the metric and vielbein formulations are distinct physical theories, although each propagates a healthy massive graviton.


If the graviton mass is comparable to the Hubble rate , then at cosmological distances the mass term can produce a repulsive gravitational effect that leads to cosmic acceleration. Because, roughly speaking, the enhanced diffeomorphism symmetry in the limit protects a small graviton mass from large quantum corrections, the choice is in fact technically natural.[32] Massive gravity thus may provide a solution to the cosmological constant problem: why do quantum corrections not cause the Universe to accelerate at extremely early times?

However, it turns out that flat and closed Friedmann–Lemaître–Robertson–Walker cosmological solutions do not exist in dRGT massive gravity with a flat reference metric.[13] Open solutions and solutions with general reference metrics suffer from instabilities.[33] Therefore, viable cosmologies can only be found in massive gravity if one abandons the cosmological principle that the Universe is uniform on large scales, or otherwise generalizes dRGT. For instance, cosmological solutions are better behaved in bigravity,[14] the theory which extends dRGT by giving dynamics. While these tend to possess instabilities as well,[34][35] those instabilities might find a resolution in the nonlinear dynamics (through a Vainshtein-like mechanism) or by pushing the era of instability to the very early Universe.[15]

3D massive gravity[edit]

A special case exists in three dimensions, where a massless graviton does not propagate any degrees of freedom. Here several ghost-free theories of a massive graviton, propagating two degrees of freedom, can be defined. In the case of topologically massive gravity[1] one has the action

with the three-dimensional Planck mass. This is three-dimensional general relativity supplemented by a Chern-Simons-like term built out of the Christoffel symbols.

More recently, a theory referred to as new massive gravity has been developed,[2] which is described by the action

Relation to gravitational waves[edit]

The 2016 discovery of gravitational waves[36] and subsequent observations have yielded constraints on the maximum mass of gravitons, if they are massive at all. Following the GW170104 event, the graviton's Compton wavelength was found to be at least 1.6×1016 m, or about 1.6 light-years, corresponding to a graviton mass of no more than 7.7×10−23 eV/c2.[16] This relation between wavelength and energy is calculated with the same formula that relates electromagnetic wavelength to photon energy. However, photons, which have only energy and no mass, are fundamentally different from massive gravitons in this respect, since the Compton wavelength of the graviton is not equal to the gravitational wavelength. Instead, the lower-bound graviton Compton wavelength is about 9×109 times greater than the gravitational wavelength for the GW170104 event, which was ~ 1,700 km. The report[16] did not elaborate on the source of this ratio. It is possible that gravitons are not the quanta of gravitational waves, or that the two phenomena are related in a different way.

See also[edit]

Further reading[edit]

Review articles
  • de Rham, Claudia (2014), "Massive Gravity", Living Reviews in Relativity, 17 (1): 7, arXiv:1401.4173, Bibcode:2014LRR....17....7D, doi:10.12942/lrr-2014-7, PMC 5256007, PMID 28179850
  • Hinterbichler, Kurt (2012), "Theoretical Aspects of Massive Gravity", Reviews of Modern Physics, 84 (2): 671–710, arXiv:1105.3735, Bibcode:2012RvMP...84..671H, doi:10.1103/RevModPhys.84.671


  1. ^ a b Deser, Stanley; Jackiw, R.; Templeton, S. (1982). "Topologically Massive Gauge Theories". Annals of Physics. 140 (2): 372–411. Bibcode:1982AnPhy.140..372D. doi:10.1016/0003-4916(82)90164-6.
  2. ^ a b Bergshoeff, Eric A.; Hohm, Olaf; Townsend, Paul K. (2009). "Massive Gravity in Three Dimensions". Phys. Rev. Lett. 102 (20): 201301. arXiv:0901.1766. Bibcode:2009PhRvL.102t1301B. doi:10.1103/PhysRevLett.102.201301. PMID 19519014.
  3. ^ Ogievetsky, V.I; Polubarinov, I.V (November 1965). "Interacting field of spin 2 and the einstein equations". Annals of Physics. 35 (2): 167–208. Bibcode:1965AnPhy..35..167O. doi:10.1016/0003-4916(65)90077-1.
  4. ^ Mukohyama, Shinji; Volkov, Mikhail S. (2018-10-22). "The Ogievetsky-Polubarinov massive gravity and the benign Boulware-Deser mode". Journal of Cosmology and Astroparticle Physics. 2018 (10): 037. arXiv:1808.04292. Bibcode:2018JCAP...10..037M. doi:10.1088/1475-7516/2018/10/037. ISSN 1475-7516.
  5. ^ Ogievetsky, V. I; Polubarinov, I. V (1965-11-01). "Interacting field of spin 2 and the einstein equations". Annals of Physics. 35 (2): 167–208. Bibcode:1965AnPhy..35..167O. doi:10.1016/0003-4916(65)90077-1. ISSN 0003-4916.
  6. ^ Curtright, T. L.; Alshal, H. (2019-10-01). "Massive dual spin 2 revisited". Nuclear Physics B. 948: 114777. Bibcode:2019NuPhB.94814777C. doi:10.1016/j.nuclphysb.2019.114777. ISSN 0550-3213.
  7. ^ Alshal, H.; Curtright, T. L. (2019-09-10). "Massive dual gravity in N spacetime dimensions". Journal of High Energy Physics. 2019 (9): 63. arXiv:1907.11537. Bibcode:2019JHEP...09..063A. doi:10.1007/JHEP09(2019)063. ISSN 1029-8479.
  8. ^ Nicolis, Alberto; Rattazzi, Riccardo; Trincherini, Enrico (2009-03-31). "Galileon as a local modification of gravity". Physical Review D. 79 (6): 064036. arXiv:0811.2197. Bibcode:2009PhRvD..79f4036N. doi:10.1103/PhysRevD.79.064036.
  9. ^ Deffayet, C.; Esposito-Farèse, G.; Vikman, A. (2009-04-03). "Covariant Galileon". Physical Review D. 79 (8): 084003. arXiv:0901.1314. Bibcode:2009PhRvD..79h4003D. doi:10.1103/PhysRevD.79.084003.
  10. ^ Curtright, Thomas L.; Fairlie, David B. (2012). "A Galileon Primer". arXiv:1212.6972 [hep-th].
  11. ^ de Rham, Claudia; Keltner, Luke; Tolley, Andrew J. (2014-07-21). "Generalized Galileon duality". Physical Review D. 90 (2): 024050. arXiv:1403.3690. Bibcode:2014PhRvD..90b4050D. doi:10.1103/PhysRevD.90.024050.
  12. ^ a b Hassan, S.F.; Rosen, Rachel A. (2012). "Bimetric Gravity from Ghost-free Massive Gravity". JHEP. 1202 (2): 126. arXiv:1109.3515. Bibcode:2012JHEP...02..126H. doi:10.1007/JHEP02(2012)126.
  13. ^ a b D'Amico, G.; de Rham, C.; Dubovsky, S.; Gabadadze, G.; Pirtskhalava, D.; Tolley, A.J. (2011). "Massive Cosmologies". Phys. Rev. D84 (12): 124046. arXiv:1108.5231. Bibcode:2011PhRvD..84l4046D. doi:10.1103/PhysRevD.84.124046.
  14. ^ a b Akrami, Yashar; Koivisto, Tomi S.; Sandstad, Marit (2013). "Accelerated expansion from ghost-free bigravity: a statistical analysis with improved generality". JHEP. 1303 (3): 099. arXiv:1209.0457. Bibcode:2013JHEP...03..099A. doi:10.1007/JHEP03(2013)099.
  15. ^ a b Akrami, Yashar; Hassan, S.F.; Könnig, Frank; Schmidt-May, Angnis; Solomon, Adam R. (2015). "Bimetric gravity is cosmologically viable". Physics Letters B. 748: 37–44. arXiv:1503.07521. Bibcode:2015PhLB..748...37A. doi:10.1016/j.physletb.2015.06.062.
  16. ^ a b c B. P. Abbott; et al. (LIGO Scientific Collaboration and Virgo Collaboration) (1 June 2017). "GW170104: Observation of a 50-Solar-Mass Binary Black Hole Coalescence at Redshift 0.2". Physical Review Letters. 118 (22): 221101. arXiv:1706.01812. Bibcode:2017PhRvL.118v1101A. doi:10.1103/PhysRevLett.118.221101.
  17. ^ L. Bernus; et al. (18 October 2019). "Constraining the Mass of the Graviton with the Planetary Ephemeris INPOP". Physical Review Letters. 123 (16): 161103. arXiv:1901.04307. Bibcode:2019arXiv190104307B. doi:10.1103/PhysRevLett.123.161103. PMID 31702347.
  18. ^ Fierz, Markus; Pauli, Wolfgang (1939). "On relativistic wave equations for particles of arbitrary spin in an electromagnetic field". Proc. Roy. Soc. Lond. A. 173 (953): 211–232. Bibcode:1939RSPSA.173..211F. doi:10.1098/rspa.1939.0140.
  19. ^ van Dam, Hendrik; Veltman, Martinus J. G. (1970). "Massive and massless Yang-Mills and gravitational fields". Nucl. Phys. B. 22 (2): 397–411. Bibcode:1970NuPhB..22..397V. doi:10.1016/0550-3213(70)90416-5. hdl:1874/4816.
  20. ^ Zakharov, Valentin I. (1970). "Linearized gravitation theory and the graviton mass". JETP Lett. 12: 312. Bibcode:1970JETPL..12..312Z.
  21. ^ Vainshtein, A.I. (1972). "To the problem of nonvanishing gravitation mass". Phys. Lett. B. 39 (3): 393–394. Bibcode:1972PhLB...39..393V. doi:10.1016/0370-2693(72)90147-5.
  22. ^ Freund, Peter G. O.; Maheshwari, Amar; Schonberg, Edmond (August 1969). "Finite-Range Gravitation". The Astrophysical Journal. 157: 857. Bibcode:1969ApJ...157..857F. doi:10.1086/150118. ISSN 0004-637X.
  23. ^ Boulware, David G.; Deser, Stanley (1972). "Can gravitation have a finite range?" (PDF). Phys. Rev. D. 6 (12): 3368–3382. Bibcode:1972PhRvD...6.3368B. doi:10.1103/PhysRevD.6.3368.
  24. ^ Creminelli, Paolo; Nicolis, Alberto; Papucci, Michele; Trincherini, Enrico (2005). "Ghosts in massive gravity". JHEP. 0509 (9): 003. arXiv:hep-th/0505147. Bibcode:2005JHEP...09..003C. doi:10.1088/1126-6708/2005/09/003.
  25. ^ de Rham, Claudia; Gabadadze, Gregory (2010). "Generalization of the Fierz-Pauli Action". Phys. Rev. D. 82 (4): 044020. arXiv:1007.0443. Bibcode:2010PhRvD..82d4020D. doi:10.1103/PhysRevD.82.044020.
  26. ^ de Rham, Claudia; Gabadadze, Gregory; Tolley, Andrew J. (2011). "Resummation of Massive Gravity". Phys. Rev. Lett. 106 (23): 231101. arXiv:1011.1232. Bibcode:2011PhRvL.106w1101D. doi:10.1103/PhysRevLett.106.231101. PMID 21770493.
  27. ^ Hassan, S.F.; Rosen, Rachel A. (2012). "Resolving the Ghost Problem in non-Linear Massive Gravity". Phys. Rev. Lett. 108 (4): 041101. arXiv:1106.3344. Bibcode:2012PhRvL.108d1101H. doi:10.1103/PhysRevLett.108.041101. PMID 22400821.
  28. ^ Hassan, S.F.; Rosen, Rachel A. (2012). "Confirmation of the Secondary Constraint and Absence of Ghost in Massive Gravity and Bimetric Gravity". JHEP. 1204 (4): 123. arXiv:1111.2070. Bibcode:2012JHEP...04..123H. doi:10.1007/JHEP04(2012)123.
  29. ^ a b Hassan, S.F.; Rosen, Rachel A. (2011). "On Non-Linear Actions for Massive Gravity". JHEP. 1107 (7): 009. arXiv:1103.6055. Bibcode:2011JHEP...07..009H. doi:10.1007/JHEP07(2011)009.
  30. ^ Hassan, S.F.; Rosen, Rachel A.; Schmidt-May, Angnis (2012). "Ghost-free Massive Gravity with a General Reference Metric". JHEP. 1202 (2): 026. arXiv:1109.3230. Bibcode:2012JHEP...02..026H. doi:10.1007/JHEP02(2012)026.
  31. ^ Hinterbichler, Kurt; Rosen, Rachel A. (2012). "Interacting Spin-2 Fields". JHEP. 1207 (7): 047. arXiv:1203.5783. Bibcode:2012JHEP...07..047H. doi:10.1007/JHEP07(2012)047.
  32. ^ de Rham, Claudia; Heisenberg, Lavinia; Ribeiro, Raquel H. (2013). "Quantum Corrections in Massive Gravity". Phys. Rev. D. 88 (8): 084058. arXiv:1307.7169. Bibcode:2013PhRvD..88h4058D. doi:10.1103/PhysRevD.88.084058.
  33. ^ de Felice, Antonio; Gümrükçüoğlu, A. Emir; Lin, Chunshan; Mukohyama, Shinji (2013). "On the cosmology of massive gravity". Class. Quantum Grav. 30 (18): 184004. arXiv:1304.0484. Bibcode:2013CQGra..30r4004D. doi:10.1088/0264-9381/30/18/184004.
  34. ^ Comelli, Denis; Crisostomi, Marco; Pilo, Luigi (2012). "Perturbations in Massive Gravity Cosmology". JHEP. 1206 (6): 085. arXiv:1202.1986. Bibcode:2012JHEP...06..085C. doi:10.1007/JHEP06(2012)085.
  35. ^ Könnig, Frank; Akrami, Yashar; Amendola, Luca; Motta, Mariele; Solomon, Adam R. (2014). "Stable and unstable cosmological models in bimetric massive gravity". Phys. Rev. D. 90 (12): 124014. arXiv:1407.4331. Bibcode:2014PhRvD..90l4014K. doi:10.1103/PhysRevD.90.124014.
  36. ^ B. P. Abbott; et al. (LIGO Scientific Collaboration and Virgo Collaboration) (2016). "Observation of Gravitational Waves from a Binary Black Hole Merger". Phys. Rev. Lett. 116 (6): 061102. arXiv:1602.03837. Bibcode:2016PhRvL.116f1102A. doi:10.1103/PhysRevLett.116.061102. PMID 26918975.