{"id":"df3b26bb-a8de-4689-8677-7b1abb9753c8","arxiv_id":"1908.06743","paper_version":1,"verdict":"REJECT","confidence":"HIGH","novelty_score":3.0,"correctness_risk":"high","formal_verification":"none","parameter_count":4,"one_line_summary":"In a toy anisotropic universe, gravitational waves acquire direction-dependent dispersion relations and anisotropy-dependent tidal acceleration, but the background equations contain a sign error that undermines the model.","lead":"The paper studies gravitational waves in a simple toy universe where space expands symmetrically in x and y but has an additional anisotropic stretch along z. It derives wave equations, dispersion relations, and tidal effects, and claims that this anisotropy changes the propagation of the waves.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Plane-wave ansatz with constant ω and k is inserted into wave equations with z- and t-dependent coefficients; the dispersion relations (38)–(39), (69), (74), (76) are therefore not valid solutions, leaving the direction-dependent speed claim unsupported.","rationale":"The reader's weakest-assumption identification, the invalid plane-wave ansatz, is the same point I would make first. For Eq. (35) with h=e^{i(ωt−kz)}, treating k as constant while b(z) varies gives an algebraic quadratic that the paper solves; but b(z) variation means the exponential with constant k is not a solution, and if k is z-dependent the phase must be an integral, not kz. The exact b=e^{mz} counterexample confirms that no plane wave of the stated form solves Eq. (35) unless m=0. The same structural problem affects the time dependence through a(t), which appears both in coefficients and in the explicit a^4 amplitudes of the x-direction modes. I also confirm the background sign error in Eq. (6): with the paper's own metric and sign convention, G00=3H^2=8πρ, not H^2=−8πρ/3, so the background energy density is inconsistent. The tracelessness argument in Sec. VI is additionally invalid because ∂0∂3(h/a^2)=0 implies h=a^2[f(z)+g(t)], not h=0. These are independent defects, but the ansatz failure alone blocks the central dispersion-relation and sub-luminal-speed claims, so the rejection is warranted.","tokens_in":13196,"tokens_out":14774,"duration_ms":143003,"concrete_test":"Specialize Eq. (35) to a(t)=1 and b(z)=e^{m z}; for h12=e^{iωt}Y(z) the equation becomes Y″−2mY′+e^{2mz}ω^2Y=0, whose exact solutions are Y=e^{m z [1±√(1−ω^2/m^2)]}. Substitute instead the paper's trial h12=ε exp{i(ωt−k12 z)} with k12 from Eq. (39): the residual is nonzero for m≠0, so (39) is not a solution of the stated wave equation. A WKB replacement with phase ∫^z k(ζ)dζ gives the corrected leading-order dispersion and differs from (39) at O(b′).","verdict_should_be":"REJECT","load_bearing_attack":"The central propagation results rest on substituting h_μν = ε_μν exp{i(ωt − kz)} with constant ω and k into equations whose coefficients depend on z through b(z) and on t through a(t). No WKB, slowly-varying, or locally-constant approximation is stated or justified. For an exponential with constant k to solve an equation with z-dependent coefficients, b(z) would have to be constant or nearly so on a wavelength; the paper never imposes this. If k is instead allowed to be the z-dependent quantity in (38)–(39), then the phase factor kz is not consistent: ∂_z exp{−i k(z) z} produces k′(z) and k″(z) terms, so the algebraic quadratic that led to (38)–(39) ceases to be the wave equation. The same defect appears in the x-direction: the amplitude in (68) carries a^4(t), so ̈h_23 contains ̈a^4 and 2iω ẏa^4 terms, and (69) is only an algebraic relation, not a solution of the time-dependent equation. Since the direction-dependent b(z)-dependence of the wave vector and the sub-luminal speed conclusion are exactly what (38)–(39) and (74)–(76) are used to assert, this invalidates the central claim. A separate independent problem is the sign error in the background G00 equation, Eq. (6), but the ansatz failure alone is sufficient to block the propagation claims.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies linearized gravitational waves in a simple anisotropic toy universe with metric ds^2 = dt^2 - a^2(t)(dx^2+dy^2) - a^2(t)b^2(z) dz^2. It imposes synchronous, transverse, traceless gauge conditions, computes first-order Ricci perturbations, writes the linearized Einstein equations for a perfect fluid, and proposes plane-wave solutions for GWs propagating along the anisotropy direction z and along the perpendicular direction x. The central claims are that the wave vector for propagation along z depends explicitly on the anisotropy function b(z) (Eqs. 38-39), that the wave vector for propagation along x is independent of b(z) (Eqs. 74-76), that tracelessness of the perturbations follows from the field equations (Eqs. 83 and 85), and that the GW speed is lower than c with the reduction controlled by a(t) and b(z) along z but only by a(t) along x. The paper ends with a computation of the tidal acceleration induced by the waves.","tokens_in":13553,"tokens_out":8450,"duration_ms":81873,"significance":"If the results were sound, the paper would provide a concrete analytical example of how a background anisotropy can modify the dispersion relation and speed of gravitational waves, complementing earlier studies of Bianchi type I universes. The paper is self-contained, does not fit parameters to data, and explicitly presents many intermediate Christoffel and Ricci perturbation expressions, which is a useful feature. However, the central propagation claims rest on a plane-wave ansatz that is not a valid solution of the variable-coefficient wave equations, and the background field equations contain a sign error. The significance is therefore currently overshadowed by technical issues that affect the main conclusions.","major_comments":[{"comment":"The background G00 equation has the wrong sign. Combining R00 = -3ä/a from Eq. (2) and R = -6ä/a - 6(ȧ/a)^2 from Eq. (3) with G00 = R00 - R/2 gives G00 = 3(ȧ/a)^2, so Einstein's equations imply 3(ȧ/a)^2 = 8πρ, not (ȧ/a)^2 = -8πρ/3. As written, Eq. (6) yields negative energy density for a standard expanding universe and is inconsistent with the perturbed 00 equation (29), which is the standard �/a = -4π/3(ρ+3p). This internal inconsistency affects the use of ρ in the dispersion relations such as Eq. (38).","section":"§II.A, Eq. (6)"},{"comment":"The claimed derivation of tracelessness is incorrect. From (1/2)∂0∂3(h/a^2) = 0 it follows that h/a^2 = f(z) + g(t), so h = a^2[f(z)+g(t)] is a general nonzero solution; the conclusion that a(t) being time-dependent forces h = 0 does not follow without additional boundary or initial conditions. Similarly, ∂0∂1(h/(a^2 b^2)) = 0 does not imply h = 0, since the most general solution is h = a^2 b^2[F(t,z) + G(x,z)] (or an equivalent separable form). Thus the statement that tracelessness follows from the perturbed Einstein equations in the manner claimed is not established.","section":"§VI, Eqs. (83) and (85)"}],"minor_comments":[{"comment":"There are numerous typographical errors, including 'Univers e' in the title, 'eqaution' in §V.B, and mislabelled references (e.g., 'arXiv:1608.01982 [gr-gc]' should be '[gr-qc]').","section":"Throughout"},{"comment":"The spatial components of the background Einstein equations appear to be missing factors of a^2 and b^2 relative to the stated index convention: with Tij = p gij and g11 = -a^2, Eq. (7) should contain an a^2 on the right-hand side, and Eq. (8) an a^2 b^2 factor. Please clarify the index convention used.","section":"§II.A, Eqs. (7)-(8)"},{"comment":"The notation is confusing: the same symbol k23 is used both as a constant wavenumber in the exponential and as a solution that depends on a(t) through Eq. (69). It would be helpful to state explicitly that Eq. (69) is a local algebraic relation at fixed t, not a global dispersion relation, if that is intended.","section":"§V, Eqs. (65)-(72)"}],"recommendation":"reject","confidential_remarks":"The manuscript has multiple load-bearing technical errors: the background 00 equation has a sign error, and the plane-wave ansatz used to derive the dispersion relations is not justified for the variable-coefficient wave equations. The tracelessness derivation is also mathematically incomplete. These issues go beyond local presentation problems and affect the main claims, so rejection seems appropriate. The paper may be resubmitted if the authors replace the inconsistent ansatz with a proper WKB treatment or solve for the t- and z-dependence of the perturbations consistently, and correct the background equations."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Dear Colleague,\n\nThe short version: this paper has a sign error in the background Einstein equations and its main dispersion relations are not actual solutions of the wave equations, so the headline claims about direction-dependent GW speed are unsupported. I'd desk-reject it.\n\nThat said, the authors are working on a legitimate topic. The toy metric—an LRS-like Bianchi I with a z-dependent anisotropy function b(z)—is a reasonable place to test how anisotropy affects GW propagation. They cite the relevant older literature (Hu, Miedema & van Leeuwen, Cho & Speliotopoulos, Adams et al.) and set up linear perturbation theory in a transparent way. The Ricci tensor computation for the background (Eqs. 2-3) is correct.\n\nThe problems start with Eq. (6). With their metric, G00 = R00 - R/2 = 3(ȧ/a)^2, so the dust density should be 8πρ = 3(ȧ/a)^2. The paper has the opposite sign, which forces ρ negative and contradicts their own perturbed equation (29). That's a clear internal inconsistency.\n\nThe bigger issue is the plane-wave ansatz. In Sec. IV, they substitute h ∝ e^{i(ωt - kz)} with constant k and ω into equations whose coefficients depend on t through a(t) and on z through b(z). The resulting dispersion relations, Eqs. (38)-(39), make k a function of z (and t), so the exponential is no longer a solution: derivatives of exp[-i k(z) z] produce k'(z) and k''(z) terms, and the t-derivatives of a time-dependent amplitude produce ȧ and ä terms. No WKB or slowly-varying justification is given. The same defect appears in the x-direction calculation, where the a^4(t) amplitude in Eq. (68) contributes extra terms to ̈h_{23}, so Eq. (69) is an algebraic relation that does not solve the PDE. Since the direction-dependent b(z)-dependence of the wave vector is the paper's central result, this is load-bearing.\n\nThere's also a logical error in the tracelessness validation. From ∂0∂3(h/a^2)=0 one cannot conclude h=0; h = a^2(t) B(z) is a nonzero solution. The same issue appears in Eq. (85).\n\nI don't see a hidden gem here. The paper is honest in presentation, and the errors are visible, but they are fundamental. The authors should fix the sign error and redo the propagation analysis with a proper mode decomposition or eikonal ansatz. As it stands, it doesn't deserve a referee's time.\n\nBest.","headline":"Sign error in background equations plus an invalid plane-wave ansatz sink the paper's central dispersion relations; the topic is fine but the execution is not sound.","tokens_in":14099,"tokens_out":13430,"would_cite":false,"duration_ms":110228,"reading_group":"maybe","serious_thinker":"no","would_accept_peer_review":false},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["04.30.-w","98.80.-k"],"model":"deepseek-v4-flash","headline":"In an anisotropic toy universe, gravitational wave speed depends on direction and falls below c.","keywords":["gravitational waves","anisotropic universe","toy model","linear perturbations","dispersion relation","transverse traceless gauge","tidal acceleration","cosmological perturbations"],"falsifier":"For a specific anisotropy $b(z) = 1 + \\varepsilon \\sin(z/\\lambda)$, substitute the plane-wave form $h_{11} = \\varepsilon_{11} \\exp\\{i(\\omega t - kz)\\}$ into Eq. (34) and check whether a single complex $k$ satisfies the equation at every $z$; if the equation demands $k(z)$, Eq. (38) is not a valid dispersion relation.","tokens_in":12918,"feed_emoji":"🌌","tokens_out":12423,"duration_ms":102257,"temperature":0.7,"pith_summary":"This paper studies how gravitational waves propagate in a simple anisotropic universe whose spatial metric is isotropic in x and y but stretched by a factor b(z) along z. The authors' central claim is that the wave number of a gravitational wave picks up an explicit dependence on b(z) when the wave travels along the anisotropy axis, while waves travelling perpendicular to that axis have wave numbers that depend only on the scale factor a(t). The paper also argues that the transverse-traceless gauge condition is not an extra assumption: the linearized Einstein equations themselves force the perturbation to be traceless, in both propagation directions, through equations (83) and (85). If correct, the model shows that cosmic anisotropy would leave a directional imprint on gravitational-wave dispersion, making wave speed direction-dependent and lower than c, which is qualitatively different from the FLRW case. The analysis is deliberately a toy model: pressure must vanish (p≃0) unless the anisotropy factor equals one, so the results are an illustration of mechanism rather than a realistic cosmology.","feed_headline":"Gravitational waves feel cosmic anisotropy on one axis only","feed_subtitle":"In a toy anisotropic universe, waves along the anisotropy axis pick up b(z), while perpendicular waves feel only the scale factor.","key_machinery":"The load-bearing device is the toy background metric (1), together with the plane-wave ansatz $h_{11}(t,z) = \\mathrm{Re}[\\varepsilon_{11} \\exp\\{i(\\omega t - k_{11} z)\\}]$ (and its $\\times$-polarization counterpart) inserted into the linearized Einstein equations. The background metric supplies $a(t)$, which sets the overall cosmological expansion, and $b(z)$, which breaks symmetry along one spatial direction; the plane-wave ansatz converts the partial differential equations (34) and (35) into algebraic dispersion relations for $k_{11}$ and $k_{12}$. The same procedure applied to $x$-propagation turns the perturbed equations (58), (63), and (64) into the wave-number formulas (69), (74), and (76). The tracelessness argument works through the extra Ricci perturbations, e.g. $\\delta R_{03}$ and $\\delta R_{01}$, which vanish only when the trace $h$ vanishes; equations (83) and (85) are the Einstein equations that enforce that vanishing.","core_discovery":"The central discovery is that the propagation of linear gravitational-wave perturbations in the metric $ds^2 = dt^2 - a^2(t)(dx^2 + dy^2) - a^2(t)b^2(z)\\,dz^2$ is direction-sensitive in a precise way. Along the anisotropy axis $z$, the two polarizations have distinct dispersion relations, $k_{11} = -ib'/b \\pm (1/b)\\sqrt{-b'^2 + a^2b^4(\\omega^2 + 8\\pi(\\rho-3p))}$ and $k_{12} = -ib'/b \\pm (1/b)\\sqrt{-b'^2 + a^2b^4\\omega^2}$, so both wave numbers depend on $b(z)$. Along a perpendicular direction $x$, the wave numbers $k_{22}$ and $k_{33}$ are equal in the pressureless case and independent of $b(z)$; they depend only on $a(t)$, through $\\omega$, $\\dot{a}$, $\\ddot{a}$, and $\\rho$. The same equations also yield the tracelessness of $h_{\\mu\\nu}$ as a consequence of particular components of Einstein's equations rather than as a gauge assumption, with the trace condition modified to $h_{33} + b^2 h_{22} = 0$ for $x$-propagation. The paper concludes that the speed of gravitational waves in this background is lower than $c$, that the lowering depends on both $a(t)$ and $b(z)$ along $z$ but only on $a(t)$ along $x$, and that in vacuum ($a \\to 1$, $b \\to 1$) the Minkowski results are recovered.","pith_inferences":["If $b(z)$ varies on scales comparable to the gravitational-wave wavelength, the constant-$k$ plane-wave ansatz is no longer an exact solution; a numerical or slowly-varying-wave treatment would be needed to see whether the $b$-dependence in equations (38)-(39) survives as a genuine dispersion effect or is an artifact of the ansatz.","The same calculation could be repeated for non-pressureless matter only if the constraint $p(b^2-1)=0$ is avoided, for instance by adding an anisotropic stress term; that would show whether the direction-dependent slowing is robust beyond dust.","A natural observational analogue is a gravitational wave crossing a region with mild anisotropic expansion: the model predicts a differential arrival time or phase between polarizations correlated with the direction to the source, though the effect would be tiny for realistic anisotropies.","The tidal-acceleration formulas imply that a detector's response depends on the anisotropy factor $b(z)$ when the wave travels perpendicular to the anisotropy axis; simulating the response for a prescribed $b(z)$ would give a concrete, testable waveform prediction."],"forward_implications":["Along the anisotropy axis the two polarizations travel with slightly different wave numbers, so an anisotropic background acts as a birefringent medium for gravitational waves.","Perpendicular to the anisotropy axis the wave number is independent of $b(z)$, so the directional dependence of dispersion could in principle be used to locate the anisotropy axis from gravitational-wave observations.","Because the wave speed is below $c$ in both directions, any observation of gravitational waves arriving with speed $c$ would constrain the allowed size of $a(t)$ and $b(z)$ in such a model.","The tracelessness of $h_{\\mu\\nu}$ follows from the field equations rather than being imposed, which means the synchronous transverse traceless gauge is consistent in this background.","In vacuum, $a(t)\\to 1$ and $b(z)\\to 1$ reproduce the standard Minkowski result, so the model's new effects are entirely due to the anisotropic expansion."],"supporting_citations":[{"why":"supplies the linearized perturbation and TT-gauge derivations that the paper follows for the connection and Ricci perturbations.","marker":"[13]"},{"why":"a companion derivation of perturbed connection and Ricci components on a similar background, used as a template for the calculation.","marker":"[18]"},{"why":"establishes that gravitational waves are non-material perturbations in Bianchi I and introduces the synchronous transverse traceless gauge used here.","marker":"[20]"},{"why":"provides the earlier result that polarizations couple and the gravitational wave acquires an effective mass in anisotropic backgrounds, which this paper compares against.","marker":"[21]"},{"why":"gives the FLRW baseline where polarizations decouple and behave as massless scalar fields.","marker":"[22]"},{"why":"the exact non-perturbative Bianchi formalism whose inhomogeneous wave equation the z-propagation result echoes.","marker":"[35]"},{"why":"supports the statement that the transverse traceless conditions remove scalar and vector pieces of the perturbation.","marker":"[41]"}],"fun_headline_variants":["Gravitational waves slow down on one axis only","Anisotropic universe yields direction-dependent GW speed","GW speed depends on direction in anisotropic cosmology","Cosmic anisotropy splits gravitational wave dispersion"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The plane-wave ansatz inserts a wave with constant frequency $\\omega$ and wave number $k$ into equations whose coefficients depend on $t$ and $z$; whenever $a(t)$ or $b(z)$ varies on scales comparable to the wavelength, the exponential is not a solution and the derived dispersion relations do not hold.","fun_headline_variants_meta":{"raw":{"variants":["Gravitational waves slow down on one axis only","Anisotropic universe yields direction-dependent GW speed","GW speed depends on direction in anisotropic cosmology","Cosmic anisotropy splits gravitational wave dispersion"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000696,"raw_usage":{"total_tokens":3198,"prompt_tokens":1051,"completion_tokens":2147,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":667,"completion_tokens_details":{"reasoning_tokens":2090}},"tokens_in":667,"tokens_out":2147,"duration_ms":17167,"temperature":1.0,"reasoning_tokens":2090,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T13:21:42.297866+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"For a specific anisotropy $b(z) = 1 + \\varepsilon \\sin(z/\\lambda)$, substitute the plane-wave form $h_{11} = \\varepsilon_{11} \\exp\\{i(\\omega t - kz)\\}$ into Eq. (34) and check whether a single complex $k$ satisfies the equation at every $z$; if the equation demands $k(z)$, Eq. (38) is not a valid dispersion relation.","supporting_citations":[{"cited_title":"(102) • For ‘+’ polarization, ε23 = 0, and we have ˆχ2(t) = 1 2a4ε2 2e{i(wt−kx)} ˆχ2 0 , ˆχ3(t) = − 1 2a4b2ε2 2e{i(wt−kx)} ˆχ3","cited_arxiv_id":null,"evidence_quote":"supplies the linearized perturbation and TT-gauge derivations that the paper follows for the connection and Ricci perturbations."},{"cited_title":"Gravitational Wave Solutions to Linearized Jordan-Brans-Dicke Theory on a Cosmological Background","cited_arxiv_id":"1707.01169","evidence_quote":"a companion derivation of perturbed connection and Ricci components on a similar background, used as a template for the calculation."},{"cited_title":"Abbott and R","cited_arxiv_id":null,"evidence_quote":"establishes that gravitational waves are non-material perturbations in Bianchi I and introduces the synchronous transverse traceless gauge used here."},{"cited_title":"Abbott and R","cited_arxiv_id":null,"evidence_quote":"provides the earlier result that polarizations couple and the gravitational wave acquires an effective mass in anisotropic backgrounds, which this paper compares against."},{"cited_title":"(33) Subtracting equation (","cited_arxiv_id":null,"evidence_quote":"gives the FLRW baseline where polarizations decouple and behave as massless scalar fields."}],"review_version":1}