ELMÉLETI KUTATÁSI TANULMÁNY
SKALÁRIS METRIKUS IDŐHULLÁMOK MINT A SÖTÉTANYAG-JELENSÉGEK LEHETSÉGES GEOMETRIAI MAGYARÁZATA
Absztrakt: A sötét anyag természetének meghatározása a modern kozmológia nyitott kérdése. A közvetlen detektálási kísérletek eddigi negatív eredményei korlátozzák egyes részecskefizikai modellek paraméterterét, de önmagukban nem cáfolják a sötét anyag létezését. Jelen tanulmány egy vizsgálandó hipotézist vázol fel: a galaktikus gravitációs anomáliák egy része a téridő metrikájának dinamikus skaláris szabadsági fokával is összefügghet. Kiindulópontként az \(f(R)\)-gravitáció metrikus formalizmusát és a gyenge téridő-perturbációk leírását használjuk. Bemutatjuk azokat a feltételeket, amelyeket egy ilyen modellnek teljesítenie kell a galaktikus dinamika, a gravitációs lencsézés, az órák frekvencia-összehasonlítása és a pulzár-időzítési hálózatok (PTA) adatainak magyarázatához. A kézirat nem állítja, hogy a hipotézis már levezetett vagy empirikusan igazolt sötétanyag-alternatíva: a konkrét \(f(R)\)-függvény, a hullámmegoldás, a megfigyelhető válaszfüggvények és az adatokon végzett modell-összehasonlítás további munkát igényel.
1. BEVEZETÉS ÉS KUTATÁSI KÉRDÉS
A standard kozmológiai modell (\(\Lambda\text{CDM}\)) szerint a kozmikus energiasűrűség jelentős részét hideg, nem barionos sötét anyag alkotja. A közvetlen detektálási kísérletek – köztük a LUX-ZEPLIN és a XENONnT – egyes jelölt részecskék esetében nem találtak meggyőző jelet a vizsgált érzékenységi tartományban. Ez korlátozza az adott modelleket, de nem jelenti a sötét anyag teljes körű cáfolatát. A kis léptékű kozmológiai problémák, például a cusp–core és a hiányzó műholdgalaxisok problémája, a standard modell és a numerikus szimulációk részleteinek további vizsgálatát igénylik; önmagukban nem bizonyítják a geometriai alternatívákat.
A tanulmány kutatási kérdése a következő: létrehozhat-e egy jól definiált módosított gravitációelméletben fellépő dinamikus skaláris metrikus szabadsági fok olyan megfigyelhető hatásokat, amelyek egyes körülmények között a sötét anyagnak tulajdonított gravitációs jelenségekhez hasonlítanak?
Ez jelenleg hipotézis. A geometriai modellnek ugyanazokat a megfigyeléseket kell reprodukálnia, mint a versengő modelleknek, és meg kell felelnie a már ellenőrzött relativisztikus, kozmológiai és laboratóriumi korlátoknak is.
2. METRIKUS FELBONTÁS ÉS AZ \(f(R)\)-GRAVITÁCIÓ
A metrikus \(f(R)\)-gravitáció általános hatása SI-konvencióban, a szokásos jelölésekkel:
$$S = \frac{c^3}{16\pi G}\int d^4x\,\sqrt{-g}\,f(R)+S_{\mathrm{matter}}.$$A metrika szerinti variáció mezőegyenlete:
$$f_R R_{\mu\nu}-\frac{1}{2}f(R)g_{\mu\nu} -\left(\nabla_\mu\nabla_\nu-g_{\mu\nu}\Box\right)f_R =\frac{8\pi G}{c^4}T_{\mu\nu},$$ahol \(f_R\equiv df/dR\), \(\Box\equiv\nabla^\alpha\nabla_\alpha\), \(T_{\mu\nu}\) pedig az energia-impulzus-tenzor. Az egyenlet nyoma:
$$3\Box f_R+f_RR-2f(R)=\frac{8\pi G}{c^4}T,$$vagy ekvivalensen
$$\Box f_R=\frac{1}{3}\left[2f(R)-f_RR+\frac{8\pi G}{c^4}T\right].$$Ez a nyomegyenlet a metrikus \(f(R)\)-elméletben dinamikus skaláris szabadsági fokot jelezhet. A pontos dinamika és a skaláris mód effektív tömege azonban az alkalmazott \(f(R)\)-függvénytől és a háttérmegoldástól függ. A fenti általános egyenletek önmagukban nem bizonyítják egy meghatározott frekvenciájú, koherens „időhullám” létezését.
A gyenge téridő-perturbációk egyik szokásos, Newtoni gauge-beli alakja:
$$ds^2=-\left(1+\frac{2\Psi(t,\mathbf{x})}{c^2}\right)c^2dt^2 +\left(1-\frac{2\Phi(t,\mathbf{x})}{c^2}\right)d\mathbf{x}^2.$$Itt \(\Psi\) és \(\Phi\) a két skaláris metrikus potenciál. A vizsgálandó ansatz:
$$\Psi(t,\mathbf{x})=\Psi_{\mathrm{barion}}(\mathbf{x})+\psi_{\mathrm{wave}}(t,\mathbf{x}).$$Ez felbontási javaslat, nem a mezőegyenletek megoldása. A további számítások előtt meg kell adni a konkrét \(f(R)\)-függvényt, a háttérgeometriát, a peremfeltételeket és a kezdeti feltételeket, majd ezekből kell meghatározni \(\psi_{\mathrm{wave}}\)-et.
Konkrét referencia-modell szükséges: Az általános \(f(R)\)-formalizmus önmagában nem ad számszerű előrejelzést. A további vizsgálatokhoz legalább egy referencia-modellt rögzíteni kell – például a Hu–Sawicki-féle \(f(R)=R-\frac{m^2c_1(R/m^2)^n}{1+c_2(R/m^2)^n}\) alakot \(n=1\) kitevővel –, valamint a háttér- és perturbációs paraméterek kezdőértékeit. A modellválasztást a Naprendszeri és kozmológiai korlátokhoz kell illeszteni.
2.1. SZŰRŐ (SCREENING) ÉS HELYI KORLÁTOK
A metrikus \(f(R)\)-elméletekben a skaláris módus nagy sűrűségű környezetben chameleon-szerű mechanizmus révén elnyomódhat. A modellnek számszerűen igazolnia kell, hogy a választott paraméterek mellett a hatás a következő korlátokon belül marad:
- Cassini-kísérlet (PPN \(\gamma\) paramétere, \(|\gamma-1|\lesssim 2\times10^{-5}\));
- Hold-lézeres távolságmérés (Nordtvedt-effektus);
- bináris pulzárok pericentrum-vándorlása (pl. PSR J0737−3039);
- gravitációs hullámok sebessége (GW170817: \(|c_{gw}/c-1|\lesssim 10^{-15}\));
- laboratóriumi és atominterferometriás tesztek.
Ha a konkrét \(f(R)\) paraméterei ezeket nem teljesítik, a hipotézis a jelen formájában kizárt.
2.2. STABILITÁS ÉS HÁTTÉRGEOMETRIA
A skaláris módus dinamikájának vizsgálatához legalább a következőket kell ellenőrizni:
- az effektív tömeg \(m_{\rm eff}^2\) előjelét és nagyságrendjét;
- a módus ghost- és gradiens-instabilitását;
- a háttérmegoldás stabilitását a releváns perturbációkra;
- a kozmológiai háttér és a lokális struktúrák összhangját.
3. GALAKTIKUS DINAMIKA: A LAPOS FORGÁSI GÖRBÉK FELTÉTELEI
A fenti metrikakonvenció mellett, lassú próbatestekre és gyenge, lassan változó mezőre a vezető rendű gyorsulás:
$$\frac{d^2\mathbf{x}}{dt^2}\approx-\nabla\Psi =-\nabla\Psi_{\mathrm{barion}}-\nabla\psi_{\mathrm{wave}}.$$Mivel \(g_{00}=-(1+2\Psi/c^2)\), ezzel ekvivalens alak:
$$\frac{d^2\mathbf{x}}{dt^2}\approx\frac{c^2}{2}\nabla g_{00}.$$Az előjel lényeges: a metrikából számított \(-c^2\nabla g_{00}/2\) nem a fenti gyorsulást adná. A körpályán mozgó, közel nemrelativisztikus próbatest esetén – tengelyszimmetrikus és közel stacionárius potenciál mellett –
$$v_c^2(r)=r\,\frac{\partial\Psi(r)}{\partial r}.$$A közel lapos forgási görbe tehát azt kívánja meg, hogy a teljes effektív potenciálra a vizsgált radiális tartományban közel teljesüljön:
$$\frac{\partial\Psi}{\partial r}\approx\frac{v_{\mathrm{flat}}^2}{r}, \qquad \Psi(r)\approx v_{\mathrm{flat}}^2\ln(r/r_0)+\mathrm{const}.$$Egy lehetséges időfüggő ansatz:
$$\psi_{\mathrm{wave}}(r,t)=A(r)\cos(\omega t+\delta).$$Ebből a radiális gradiens:
$$\frac{\partial\psi_{\mathrm{wave}}}{\partial r} =A'(r)\cos(\omega t+\delta).$$Az oszcilláló gradiens előjeles időátlaga ideális, teljes periódusú átlagolásban nulla. Az abszolút érték időátlaga nem nulla, de önmagában nem helyettesíthető be a körpálya-egyenletbe: a pályadinamika a tényleges, időfüggő erőtől és annak a pályára gyakorolt hatásától függ. Ezért a lapos forgási görbe nem következik pusztán a harmonikus időfüggés feltételezéséből.
Rezonancia lehetősége: Bár a harmonikus oszcilláció időátlaga nulla, rezonancia esetén – amikor a hullám frekvenciája közel esik a keringési frekvenciához vagy annak egész számú többszöröséhez – a pályaperturbáció felhalmozódhat, és átlagolás után nem tűnik el. Ezt a lehetőséget a modellnek külön meg kell vizsgálnia, mert a lapos forgási görbék szempontjából ez a releváns eset.
A modellnek a mezőegyenletek megoldásából kell levezetnie \(A(r)\)-t és \(\omega\)-t, majd numerikusan meg kell mutatnia, hogy a teljes potenciál reprodukálja-e a megfigyelt forgási görbéket. A Tully–Fisher-reláció magyarázata külön ellenőrzést igényel; azt a jelenlegi ansatz még nem vezeti le.
4. SAJÁTIDŐ, ÓRÁK ÉS MEGFIGYELHETŐ FREKVENCIAVÁLTOZÁS
Állandó térbeli koordinátán nyugvó megfigyelő világvonalán a sajátidő a metrikából:
$$d\tau=\sqrt{-g_{00}}\,dt =\sqrt{1+\frac{2\Psi(t,\mathbf{x})}{c^2}}\,dt.$$Ha a potenciál gyenge, \(|\Psi|/c^2\ll1\), akkor első rendben:
$$d\tau\approx\left(1+\frac{\Psi(t,\mathbf{x})}{c^2}\right)dt.$$Ha a barionos potenciál időben állandó, és a hullám amplitúdója az adott megfigyelő helyén \(A(\mathbf{x})\), akkor az előző ansatzból:
$$\tau(t)-\tau(0)\approx \left(1+\frac{\Psi_{\mathrm{barion}}}{c^2}\right)t +\frac{A(\mathbf{x})}{\omega c^2} \left[\sin(\omega t+\delta)-\sin\delta\right].$$Ez a kifejezés a feltételezett metrikából következő sajátidő-korrekciót mutatja; nem önálló bizonyíték arra, hogy a feltételezett hullám fizikailag létezik.
- Helyi mérés: egy ideális helyi óra a sajátidejét méri. A helyi fizika leírása ettől nem válik önellentmondóvá, de ebből nem következik automatikusan, hogy minden helyi vagy összehasonlító mérésben nulla az anomália.
- Távoli frekvencia-összehasonlítás: a kibocsátott és vett frekvencia aránya függ a kibocsátó és a vevő világvonalától, a metrikától, valamint a jel terjedésétől. Ezért a $$\frac{\Delta\nu(t)}{\nu_0}\sim\frac{A(\mathbf{x})}{c^2}\cos(\omega t+\delta)$$ alak legfeljebb szemléltető ansatz; általános megfigyelhető jóslatként csak a teljes jelterjedési és mérési modell levezetése után használható.
Az ellenőrizhető előrejelzéshez meg kell adni a forrás és a vevő helyét, a hullám fázisát mindkét végponton, a jel terjedési idejét, valamint a mérési protokollt. Az ismert gravitációs vöröseltolódást és a műszeres, illetve környezeti hatásokat külön kell modellezni.
4.1. NAGYSÁGRENDI BECSLÉSEK ÉS DETEKTÁLHATÓSÁG
A hipotézis tesztelhetőségéhez számszerűen meg kell adni a lehetséges \(A/c^2\) és \(\omega\) tartományokat. Például a PTA által mért maradékok szintjén (\(10^{-15}\) relatív frekvenciaváltozás, \(\omega\sim 10^{-8}\,\mathrm{Hz}\) nagyságrend) a szükséges amplitúdó \(A/c^2\sim 10^{-15}\). A Naprendszeri és laboratóriumi órák ennél jóval szigorúbb korlátokat adnak. A modellnek meg kell mondania, hogy a releváns paramétertartományban a jel detektálható-e a jelenlegi műszerekkel; ha nem, ez is fontos (negatív) eredmény.
5. PULZÁR-IDŐZÍTÉSI HÁLÓZATOK (PTA): TESZTELÉSI TERV
A PTA-k a pulzárok érkezési idejének maradékai közötti térszögfüggő korrelációkat mérik. A NANOGrav 15 éves adatsorában korrelált, sztochasztikus jelre találtak bizonyítékot; a mért korrelációs mintázat összhangban áll a gravitációshullám-háttér általános relativitáselméletből várt Hellings–Downs-görbéjével. Ez nem bizonyítja, hogy minden lehetséges alternatív hozzájárulás kizárt, de nem indokolja azt az állítást sem, hogy a megfigyelt jel már igazolja a jelen tanulmányban felvetett időhullámot.
Fontos pontosítás: nincs egyetlen, minden „skaláris időhullámra” érvényes PTA-szögkorrelációs függvény. A válaszfüggvény a konkrét gravitációelmélettől, a polarizációs módtól, a hullámterjedéstől, a forrás- és pulzártag kezelésétől, valamint a normalizálási konvenciótól függ. Ezért az alábbi kézirat nem közöl ellenőrizetlen zárt alakú skaláris korrelációs képletet, és nem ad meg belőle származtatott szögértékeket.
A GR tenzori háttér egyik szokásos, különböző pulzárokra használt Hellings–Downs-korrelációs alakja:
$$\Gamma_{\mathrm{HD}}(\zeta)=\frac{1}{2} +\frac{3}{2}x\ln x-\frac{1}{4}x,\qquad x=\frac{1-\cos\zeta}{2},$$ahol \(\zeta\) a két pulzár közötti szög. A pontos normalizálást és az azonos pulzárhoz tartozó sajátkorrelációs tagot az elemzésben következetesen kell kezelni. Más polarizációs módokra külön válaszfüggvényt kell levezetni vagy szakirodalomból, a használt konvenció pontos megjelölésével átvenni.
| Modell | Korrelációs minta | A kézirat jelenlegi státusza |
|---|---|---|
| GR tenzori háttér | Hellings–Downs típusú kvadrupólusos mintázat | Ismert referencia-modell; az adatokkal való illesztést az adott elemzésben kell elvégezni. |
| Skaláris transzverzális (breathing) | A polarizációs válaszfüggvényből származtatandó | Nem adható meg univerzális számsor a konkrét modell és konvenció nélkül. |
| Skaláris longitudinális | A polarizációs válaszfüggvényből származtatandó; egyes modellekben erős szögfüggés jelentkezhet | Különösen fontos a pulzártag, a normalizálás és a tényleges adatkorlátok kezelése. |
5.1. Javasolt verifikációs protokoll
A hipotézis teszteléséhez a nyilvános PTA-adatokon összehasonlítható, teljesen definiált modelleket kell illeszteni. Egy sematikus kovarianciamodell:
$$\mathbf{C}_{ab}(f)= A_{\mathrm{TT}}^2\,\Gamma_{\mathrm{HD}}(\zeta_{ab})\,S_{\mathrm{TT}}(f) +A_{\mathrm{scalar}}^2\,\Gamma_{\mathrm{scalar}}(\zeta_{ab})\,S_{\mathrm{scalar}}(f) +\mathbf{N}_{ab}(f).$$Ez csak vázlat: a gyakorlati PTA-elemzésben a kovariancia a pulzárok időzítési adataihoz, zajmodelljeihez, frekvenciafüggéséhez és a mérési eljáráshoz igazítandó; a fenti alak nem tekinthető kész illesztési implementációnak. A \(\Gamma_{\mathrm{scalar}}\) függvényt csak a konkrét elmélet és polarizációs válasz levezetése után szabad behelyettesíteni.
A tesztnek tartalmaznia kell a zaj- és efemerisz-modellek kezelését, a modellparaméterek priorjait, a posterior eloszlásokat, a Bayes-faktort vagy más előre rögzített modell-összehasonlító mérőszámot, valamint a robusztussági ellenőrzéseket. Az alacsony multipólusú komponensek órával vagy efemerisszel való összetévesztésének lehetőségét tesztelni kell, nem pedig előre a hipotézis bizonyítékaként értelmezni.
6. GRAVITÁCIÓS LENCSÉZÉS ÉS KOZMOLÓGIAI KORLÁTOK
A galaktikus forgási görbék és a gravitációs lencsézés külön megfigyelési próbák. A fényelhajlás a skaláris metrikus perturbációk megfelelő kombinációjától függ; ezért egy olyan modell, amely a nemrelativisztikus mozgást a \(\Psi\) potenciállal írja le, nem jósolja meg automatikusan a helyes lencsézést is. A modellből le kell vezetni a fényterjedést, majd a számított lencsézési jelet független adatokkal kell összevetni.
Hasonlóképpen vizsgálni kell a galaxishalmazok dinamikáját, a kozmikus mikrohullámú háttér megfigyeléseit, a szerkezetképződést és a gravitációelmélet helyi korlátait. A közvetlen sötétanyag-keresések negatív eredményei önmagukban nem bizonyítják a geometriai magyarázatot.
7. KÖVETKEZTETÉS ÉS KÖVETKEZŐ LÉPÉSEK
A skaláris metrikus időhullám ebben a dokumentumban kutatási hipotézis: azt vizsgálja, hogy a módosított gravitációelmélet dinamikus skaláris szabadsági foka adhat-e a sötét anyagnak tulajdonított jelenségekhez hasonló megfigyelhető hatásokat. A jelenlegi felírás nem bizonyítja, hogy a modell reprodukálja a lapos forgási görbéket, a Tully–Fisher-relációt, a gravitációs lencsézést vagy a NANOGrav által mért korrelációkat.
A következő, ellenőrizhető munkafázisok:
- egy konkrét \(f(R)\)-függvény, háttérmegoldás és stabilitási feltételek rögzítése;
- a skaláris perturbáció egyenletének megoldása és a hullám amplitúdójának, frekvenciájának levezetése;
- galaktikus forgásigörbék numerikus illesztése a barionos tömegeloszlás figyelembevételével;
- a lencsézési jel külön kiszámítása és független adatokkal való összehasonlítása;
- a PTA-válaszfüggvény levezetése, majd előre rögzített statisztikai protokoll szerinti adatillesztés.
Az elmélet értékét nem az adja, hogy a sötét anyagot már megoldott kérdésnek nyilvánítja, hanem az, hogy világos, számszerű és cáfolható előrejelzéseket tud-e adni. A hipotézis megerősítéséhez és elvetéséhez egyaránt reprodukálható számításokra és adatokra van szükség.
7.1. FALSZIFIKÁCIÓS KRITÉRIUMOK
Egy tudományos hipotézis értékét az adja, hogy milyen megfigyelés cáfolná. A jelen modell esetében legalább a következő kritériumokat kell előre rögzíteni:
- Ha a konkrét \(f(R)\) paraméterek mellett a Naprendszeri kísérletek (Cassini, LLI, Nordtvedt) korlátai sérülnek, a modell kizárt.
- Ha a galaktikus forgási görbék illesztése nem jobb, mint a barionos anyag + \(\Lambda\)CDM modell, a geometriai hozzájárulás nem szükséges.
- Ha a PTA szögkorreláció a konkrét skaláris válaszfüggvény alapján szignifikánsan eltér a mért mintázattól, a modell kizárt.
- Ha a gravitációs lencsézés mért jele nem egyezik a modell előrejelzésével, a modell kizárt.
8. HIVATKOZÁSOK ÉS TOVÁBBI OLVASMÁNY
- Agazie, G. et al. (NANOGrav Collaboration), “The NANOGrav 15-year Data Set: Evidence for a Gravitational-Wave Background”, The Astrophysical Journal Letters 951, L8 (2023). https://doi.org/10.3847/2041-8213/acdac6
- Chen, Z.-C. et al., “Search for Nontensorial Gravitational-Wave Backgrounds in the NANOGrav 15-year Dataset”, Physical Review D 109, 084045 (2024). https://doi.org/10.1103/PhysRevD.109.084045
- Sotiriou, T. P., & Faraoni, V., “\(f(R)\) Theories of Gravity”, Reviews of Modern Physics 82, 451 (2010). https://doi.org/10.1103/RevModPhys.82.451
- De Felice, A., & Tsujikawa, S., “\(f(R)\) Theories”, Living Reviews in Relativity 13, 3 (2010). https://doi.org/10.12942/lrr-2010-3
- Hu, W., & Sawicki, I., “Models of \(f(R)\) Cosmic Acceleration that Evade Solar-System Tests”, Physical Review D 76, 064004 (2007). https://doi.org/10.1103/PhysRevD.76.064004
- Will, C. M., “The Confrontation between General Relativity and Experiment”, Living Reviews in Relativity 17, 4 (2014). https://doi.org/10.12942/lrr-2014-4
THEORETICAL RESEARCH NOTE
SCALAR METRIC TIME WAVES AS A POSSIBLE GEOMETRIC EXPLANATION OF DARK-MATTER PHENOMENA
Abstract: The nature of dark matter remains an open question in modern cosmology. Null results from direct-detection experiments constrain portions of the parameter space of specific particle candidates, but do not by themselves disprove the existence of dark matter. This paper outlines a hypothesis to be tested: some galactic gravitational anomalies might be associated with a dynamical scalar degree of freedom in the spacetime metric. The framework starts from the metric formulation of \(f(R)\) gravity and weak-field metric perturbations. We identify the conditions that such a model would have to satisfy to account for galactic dynamics, gravitational lensing, clock-frequency comparisons, and pulsar timing array (PTA) observations. This manuscript does not claim that the hypothesis is already a derived or empirically validated alternative to dark matter. A specific \(f(R)\) function, a wave solution, observable response functions, and model comparisons against data remain necessary.
1. INTRODUCTION AND RESEARCH QUESTION
In the standard cosmological model (\(\Lambda\text{CDM}\)), a substantial fraction of the cosmic energy density is attributed to cold, non-baryonic dark matter. Direct-detection experiments, including LUX-ZEPLIN and XENONnT, have not found convincing signals for certain candidate particles within the parameter ranges probed. This constrains those models but does not disprove dark matter in general. Small-scale issues such as the cusp–core and missing-satellites problems motivate continued study of the details of the standard model and numerical simulations; they do not by themselves establish a geometric alternative.
The research question is: can a well-defined modified-gravity theory with a dynamical scalar metric degree of freedom produce observable effects that, in some regimes, resemble gravitational phenomena usually attributed to dark matter?
This remains a hypothesis. A geometric model must reproduce the relevant observations at least as well as competing models and must also satisfy existing relativistic, cosmological, and laboratory constraints.
2. METRIC DECOMPOSITION AND \(f(R)\) GRAVITY
In SI conventions, the action for metric \(f(R)\) gravity can be written as
$$S = \frac{c^3}{16\pi G}\int d^4x\,\sqrt{-g}\,f(R)+S_{\mathrm{matter}}.$$Variation with respect to the metric gives
$$f_R R_{\mu\nu}-\frac{1}{2}f(R)g_{\mu\nu} -\left(\nabla_\mu\nabla_\nu-g_{\mu\nu}\Box\right)f_R =\frac{8\pi G}{c^4}T_{\mu\nu},$$where \(f_R\equiv df/dR\), \(\Box\equiv\nabla^\alpha\nabla_\alpha\), and \(T_{\mu\nu}\) is the stress-energy tensor. Taking the trace yields
$$3\Box f_R+f_RR-2f(R)=\frac{8\pi G}{c^4}T,$$or equivalently
$$\Box f_R=\frac{1}{3}\left[2f(R)-f_RR+\frac{8\pi G}{c^4}T\right].$$This trace equation indicates a dynamical scalar degree of freedom in metric \(f(R)\) gravity. Its dynamics and effective mass depend on the chosen function \(f(R)\) and the background solution. These general equations alone do not establish the existence of a coherent “time wave” at any particular frequency.
A standard Newtonian-gauge form for weak metric perturbations is
$$ds^2=-\left(1+\frac{2\Psi(t,\mathbf{x})}{c^2}\right)c^2dt^2 +\left(1-\frac{2\Phi(t,\mathbf{x})}{c^2}\right)d\mathbf{x}^2.$$Here \(\Psi\) and \(\Phi\) are the two scalar metric potentials. The ansatz to be tested is
$$\Psi(t,\mathbf{x})=\Psi_{\mathrm{baryon}}(\mathbf{x})+\psi_{\mathrm{wave}}(t,\mathbf{x}).$$This is a proposed decomposition, not a solution of the field equations. Before drawing further conclusions, the specific \(f(R)\) function, background geometry, boundary conditions, and initial conditions must be specified and used to solve for \(\psi_{\mathrm{wave}}\).
A concrete reference model is required: The general \(f(R)\) formalism alone does not yield numerical predictions. At least one reference model must be fixed – for example the Hu–Sawicki form \(f(R)=R-\frac{m^2c_1(R/m^2)^n}{1+c_2(R/m^2)^n}\) with \(n=1\) – together with initial values for background and perturbation parameters. The model choice must be matched to Solar-System and cosmological constraints.
2.1. SCREENING AND LOCAL CONSTRAINTS
In metric \(f(R)\) theories, the scalar mode can be suppressed in dense environments through chameleon-like screening. The model must demonstrate numerically that, for the chosen parameters, the effect stays within the following bounds:
- Cassini experiment (PPN \(\gamma\) parameter, \(|\gamma-1|\lesssim 2\times10^{-5}\));
- lunar laser ranging (Nordtvedt effect);
- binary-pulsar periastron advance (e.g. PSR J0737−3039);
- gravitational-wave speed (GW170817: \(|c_{gw}/c-1|\lesssim 10^{-15}\));
- laboratory and atom-interferometry tests.
If the chosen \(f(R)\) parameters violate these bounds, the hypothesis is ruled out in its present form.
2.2. STABILITY AND BACKGROUND GEOMETRY
Investigating the scalar-mode dynamics requires checking at least:
- the sign and magnitude of the effective mass \(m_{\rm eff}^2\);
- ghost and gradient instabilities of the mode;
- stability of the background solution against relevant perturbations;
- consistency between the cosmological background and local structures.
3. GALACTIC DYNAMICS: CONDITIONS FOR FLAT ROTATION CURVES
With the metric convention above, the leading-order acceleration of slow test bodies in a weak, slowly varying field is
$$\frac{d^2\mathbf{x}}{dt^2}\approx-\nabla\Psi =-\nabla\Psi_{\mathrm{baryon}}-\nabla\psi_{\mathrm{wave}}.$$Since \(g_{00}=-(1+2\Psi/c^2)\), the equivalent expression is
$$\frac{d^2\mathbf{x}}{dt^2}\approx\frac{c^2}{2}\nabla g_{00}.$$The sign matters: \(-c^2\nabla g_{00}/2\) would not agree with the acceleration derived from the stated metric. For a nearly circular orbit in an axisymmetric, approximately stationary potential,
$$v_c^2(r)=r\,\frac{\partial\Psi(r)}{\partial r}.$$A flat rotation curve therefore requires the total effective potential to satisfy approximately, over the radial range of interest,
$$\frac{\partial\Psi}{\partial r}\approx\frac{v_{\mathrm{flat}}^2}{r}, \qquad \Psi(r)\approx v_{\mathrm{flat}}^2\ln(r/r_0)+\mathrm{const}.$$A possible time-dependent ansatz is
$$\psi_{\mathrm{wave}}(r,t)=A(r)\cos(\omega t+\delta).$$Its radial gradient is
$$\frac{\partial\psi_{\mathrm{wave}}}{\partial r} =A'(r)\cos(\omega t+\delta).$$The signed average of an ideal oscillating gradient over a full period is zero. The time average of its absolute value is nonzero, but it cannot simply be substituted into the circular-orbit equation: orbital dynamics depend on the actual time-dependent force and its effect on the trajectory. A flat rotation curve therefore does not follow from the assumed harmonic time dependence alone.
Possible resonance: Although the time average of a harmonic oscillation vanishes, resonant accumulation can occur when the wave frequency is close to the orbital frequency or an integer multiple of it; the orbital perturbation then does not average to zero. The model must examine this case separately, as it is the relevant one for flat rotation curves.
The model must derive \(A(r)\) and \(\omega\) from a solution of the field equations, then numerically test whether the total potential reproduces observed rotation curves. Explaining the Tully–Fisher relation requires a separate derivation; it is not established by the present ansatz.
4. PROPER TIME, CLOCKS, AND OBSERVABLE FREQUENCY SHIFTS
For an observer at fixed spatial coordinates, proper time follows from the metric:
$$d\tau=\sqrt{-g_{00}}\,dt =\sqrt{1+\frac{2\Psi(t,\mathbf{x})}{c^2}}\,dt.$$For a weak potential, \(|\Psi|/c^2\ll1\), the first-order expression is
$$d\tau\approx\left(1+\frac{\Psi(t,\mathbf{x})}{c^2}\right)dt.$$If the baryonic potential is time-independent and the wave amplitude at the observer is \(A(\mathbf{x})\), the ansatz gives
$$\tau(t)-\tau(0)\approx \left(1+\frac{\Psi_{\mathrm{baryon}}}{c^2}\right)t +\frac{A(\mathbf{x})}{\omega c^2} \left[\sin(\omega t+\delta)-\sin\delta\right].$$This expression gives the proper-time correction implied by the assumed metric; it is not independent evidence that the assumed wave physically exists.
- Local measurements: an ideal local clock measures its own proper time. This does not make local physics inconsistent, but it does not automatically imply that every local or comparative measurement will show zero anomaly.
- Remote frequency comparison: the emitted-to-received frequency ratio depends on the emitter and receiver worldlines, the metric, and signal propagation. Consequently, $$\frac{\Delta\nu(t)}{\nu_0}\sim\frac{A(\mathbf{x})}{c^2}\cos(\omega t+\delta)$$ is at most an illustrative ansatz. It can be used as an observable prediction only after deriving the full propagation and measurement model.
A testable prediction must specify the source and receiver locations, the wave phase at both endpoints, signal travel time, and measurement protocol. Known gravitational redshift and instrumental or environmental effects must be modelled separately.
4.1. ORDER-OF-MAGNITUDE ESTIMATES AND DETECTABILITY
Testing the hypothesis requires numerical ranges for \(A/c^2\) and \(\omega\). For example, at the level of PTA residuals (\(10^{-15}\) relative frequency shift, \(\omega\sim 10^{-8}\,\mathrm{Hz}\)), the required amplitude is \(A/c^2\sim 10^{-15}\). Solar-System and laboratory clocks give far tighter bounds. The model must state whether, in the relevant parameter range, the signal is detectable with current instruments; if not, this too is an important (negative) result.
5. PULSAR TIMING ARRAYS (PTAs): A TESTING PLAN
PTAs measure angular cross-correlations between pulsar timing residuals. The NANOGrav 15-year dataset provides evidence for a correlated stochastic signal; its correlation pattern is consistent with the Hellings–Downs curve expected for a gravitational-wave background in general relativity. This does not prove that every possible alternative contribution has been excluded, but it also does not justify claiming that the observed signal already verifies the time-wave hypothesis proposed here.
Important correction: there is no single angular PTA correlation function that applies to every “scalar time wave.” The response depends on the specific gravity theory, polarization mode, propagation, treatment of source and pulsar terms, and normalization convention. This manuscript therefore does not retain an unverified closed-form scalar correlation equation or numerical angular values derived from it.
A commonly used Hellings–Downs form for distinct pulsars in the GR tensor-background model is
$$\Gamma_{\mathrm{HD}}(\zeta)=\frac{1}{2} +\frac{3}{2}x\ln x-\frac{1}{4}x,\qquad x=\frac{1-\cos\zeta}{2},$$where \(\zeta\) is the angular separation between the pulsars. The normalization and the auto-correlation term for the same pulsar must be treated consistently in the analysis. Other polarization modes require their own response functions, derived for the model or adopted from the literature with the convention explicitly stated.
| Model | Correlation pattern | Status in this manuscript |
|---|---|---|
| GR tensor background | Hellings–Downs-type quadrupolar pattern | Established reference model; fit the data for the specific analysis. |
| Scalar transverse (breathing) | Must be derived from the polarization response function | No universal numerical values without a specific model and convention. |
| Scalar longitudinal | Must be derived from the polarization response function; strong angular dependence can arise in some models | Requires careful treatment of pulsar terms, normalization, and empirical constraints. |
5.1. Proposed verification protocol
Testing the hypothesis requires fitting fully specified, comparable models to public PTA data. A schematic covariance model is
$$\mathbf{C}_{ab}(f)= A_{\mathrm{TT}}^2\,\Gamma_{\mathrm{HD}}(\zeta_{ab})\,S_{\mathrm{TT}}(f) +A_{\mathrm{scalar}}^2\,\Gamma_{\mathrm{scalar}}(\zeta_{ab})\,S_{\mathrm{scalar}}(f) +\mathbf{N}_{ab}(f).$$This is only a schematic expression. In practical PTA analyses, the covariance must be tailored to the timing data, pulsar noise models, frequency dependence, and measurement procedure; the equation above is not a ready-to-run implementation. The function \(\Gamma_{\mathrm{scalar}}\) should be inserted only after deriving the response for a specified theory and polarization mode.
The test should specify treatment of noise and ephemeris models, parameter priors, posterior distributions, a Bayes factor or other pre-registered model-comparison statistic, and robustness checks. Possible degeneracy between low-multipole components, clock errors, and ephemeris errors must be tested rather than assumed to be evidence for the hypothesis.
6. GRAVITATIONAL LENSING AND COSMOLOGICAL CONSTRAINTS
Galactic rotation curves and gravitational lensing are separate observational tests. Light deflection depends on the appropriate combination of scalar metric perturbations; a model that describes non-relativistic motion using \(\Psi\) does not automatically predict the correct lensing signal. Light propagation must be derived from the same metric, and the predicted lensing must be compared with independent observations.
The model should also be tested against galaxy-cluster dynamics, cosmic microwave background observations, structure formation, and local constraints on gravity. Null results from direct dark-matter searches do not by themselves establish a geometric explanation.
7. CONCLUSION AND NEXT STEPS
In this document, scalar metric time waves are a research hypothesis: the question is whether a dynamical scalar degree of freedom in modified gravity can produce observable effects resembling phenomena attributed to dark matter. The present formulation does not demonstrate that the model reproduces flat rotation curves, the Tully–Fisher relation, gravitational lensing, or the correlations measured by NANOGrav.
The next verifiable steps are:
- specify a concrete \(f(R)\) function, background solution, and stability conditions;
- solve the scalar perturbation equation and derive the wave amplitude and frequency;
- numerically fit galactic rotation curves while accounting for the baryonic mass distribution;
- calculate the lensing signal separately and compare it with independent data;
- derive the PTA response function and then fit the data using a pre-specified statistical protocol.
The value of the theory will not come from declaring dark matter a solved problem, but from producing clear, quantitative, falsifiable predictions. Reproducible calculations and data are needed both to support and to reject the hypothesis.
7.1. FALSIFICATION CRITERIA
The value of a scientific hypothesis lies in the observations that would refute it. At least the following criteria must be pre-registered:
- If Solar-System constraints (Cassini, LLI, Nordtvedt) are violated for the chosen \(f(R)\) parameters, the model is ruled out.
- If the fit to galactic rotation curves is no better than baryons + \(\Lambda\)CDM, the geometric contribution is not required.
- If the PTA angular correlation predicted by the specific scalar response function departs significantly from the measured pattern, the model is ruled out.
- If the measured lensing signal disagrees with the model prediction, the model is ruled out.
8. REFERENCES AND FURTHER READING
- Agazie, G. et al. (NANOGrav Collaboration), “The NANOGrav 15-year Data Set: Evidence for a Gravitational-Wave Background”, The Astrophysical Journal Letters 951, L8 (2023). https://doi.org/10.3847/2041-8213/acdac6
- Chen, Z.-C. et al., “Search for Nontensorial Gravitational-Wave Backgrounds in the NANOGrav 15-year Dataset”, Physical Review D 109, 084045 (2024). https://doi.org/10.1103/PhysRevD.109.084045
- Sotiriou, T. P., & Faraoni, V., “\(f(R)\) Theories of Gravity”, Reviews of Modern Physics 82, 451 (2010). https://doi.org/10.1103/RevModPhys.82.451
- De Felice, A., & Tsujikawa, S., “\(f(R)\) Theories”, Living Reviews in Relativity 13, 3 (2010). https://doi.org/10.12942/lrr-2010-3
- Hu, W., & Sawicki, I., “Models of \(f(R)\) Cosmic Acceleration that Evade Solar-System Tests”, Physical Review D 76, 064004 (2007). https://doi.org/10.1103/PhysRevD.76.064004
- Will, C. M., “The Confrontation between General Relativity and Experiment”, Living Reviews in Relativity 17, 4 (2014). https://doi.org/10.12942/lrr-2014-4