REVIEW 4 major objections 4 minor 43 references
Black holes and their shadows in $F(R)$ gravity
T0 review · 4 major / 4 minor · reviewed 2026-08-11 · deepseek-v4-flash
Pith's one-line read F(R) gravity can shift black hole shadows into the observed EHT ranges.
desk verdict Novel ODE for F(R) black hole perturbations, but the EHT 'prediction' reduces to tuning a free metric seed; the technical core deserves review, the conclusion does not. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The load-bearing object is the third-order differential equation (15) that the F(R) field equations impose on $F_R(r)$ for a given static spherically symmetric metric function $\nu(r)$, together with its perturbative reduction (29) near the photon sphere. The reconstruction pipeline—solve (15) for $F_R(r)$, use (13) to obtain $\lambda(r)$, compute the scalar curvature $R(r)$, invert to $r(R)$, and integrate $F_R(R)$ to recover $F(R)$—converts an assumed geometry into a concrete F(R) theory. The observable side is the shadow formula $r_{\mathrm{sh}}=r e^{-\nu(r)}|_{r=r_{\mathrm{ph}}}$, which in the perturbative branch collapses to $r_{\mathrm{sh}}=3\sqrt{3}M(1-\varepsilon\nu_1(3M))$; the photon-sphere shift itself drops out of this expression.
What would settle it
For a specified F(R) model, solve the linearized field equation (29) (or the full equation (15)) to compute $\nu_1(r)$, then evaluate $-3\sqrt{3}\,\varepsilon\nu_1(3M)$ and check it against the inequalities (55) and (56); a value outside the relevant window shows that model is incompatible with the corresponding EHT shadow observation.
Extended reading notes
Core claim
The paper's central claim is that a perturbative reconstruction of F(R) from an assumed geometry near the photon sphere is enough to make black hole shadows compatible with observations. Writing $F_R=1+\frac12\varepsilon f_R$ and $\nu=\nu_0+\varepsilon\nu_1$, the vacuum field equations reduce to an inhomogeneous linear third-order equation (29) for $f_R(r)$. In the region $r\sim 3M$, a regular branch of solutions exists whenever $C_{\nu_1}\equiv(4M\nu_1'''-\frac{32}{3}\nu_1''+\frac{8}{3M}\nu_1')|_{r=3M}$ is finite and nonzero; along that branch the shadow radius is $r_{\mathrm{sh}}=3\sqrt{3}M(1-\varepsilon\nu_1(3M))$. Because $\nu_1(3M)$ is not fixed by the reconstruction, the M87* and Sgr A* bounds reduce to simple inequalities on the combination $-3\sqrt{3}\,\varepsilon\nu_1(3M)$, and parameter choices satisfying them exist. The paper therefore claims that F(R)-gravity black holes can easily pass the Event Horizon Telescope constraints while remaining perturbatively close to Schwarzschild.
Load-bearing premise
The load-bearing assumption is that the metric perturbation $\nu_1(r)$ may be chosen freely and then realized by some local F(R) theory; the paper never derives $\nu_1$ from a specified action, so the tuning freedom that fits the shadow data could disappear for a concrete model.
Editorial extensions
If this is right
- F(R) gravity remains viable at horizon scales: the current M87* and Sgr A* shadow measurements do not exclude the theory class.
- Any metric perturbation $\nu_1$ with $C_{\nu_1}$ finite and nonzero at $r=3M$ determines a locally reconstructed F(R) whose shadow radius is controlled by $\varepsilon\nu_1(3M)$ alone.
- The constraints become explicit numerical targets for model building: $-0.5<-3\sqrt{3}\,\varepsilon\nu_1(3M)<1.1$ for M87* and $-1.01\lesssim -3\sqrt{3}\,\varepsilon\nu_1(3M)\lesssim 0.36$ for Sgr A*.
- A shift in the photon-sphere radius at first order does not feed into the shadow radius, because $d(r/\sqrt{1-2M/r})/dr$ vanishes at $r=3M$; the resulting first-order shadow formula depends only on $\nu_1(3M)$.
- The non-perturbative $\nu=0$ power-law branch yields no finite photon sphere, so not every F(R) vacuum solution can describe observed shadows.
Reading between the lines
- The paper leaves open which explicit F(R) action realizes a given $\nu_1$; committing to a concrete model fixes $\nu_1$ and may remove the tuning freedom used to satisfy the inequalities.
- The same perturbative scheme should transfer to other metric-modified gravity theories: whenever the field equations linearize to a scalar equation for a metric perturbation near the photon sphere, EHT constraints reduce to a one-parameter test.
- The 'easily pass' conclusion is an existence statement about parameter choices, not a prediction for generic F(R) models; the method supplies a diagnostic for classifying individual models by their value of $\varepsilon\nu_1(3M)$.
- A single shadow measurement constrains only the combination $\varepsilon\nu_1(3M)$, not the full radial dependence of the deviation; distinguishing among F(R) models will require additional observations or a theoretical prior on $\nu_1$.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper investigates photon-sphere and black-hole-shadow radii in F(R) gravity. For a general static spherically symmetric metric (5) the authors derive, in Eq. (15), a third-order ODE for F_R(r); they then specialize to a perturbative setup in which the metric is a small deviation from Schwarzschild and F_R=1+(ε/2)f_R. Solving the linearized Eq. (29) locally near r=3M, they obtain f_R and λ_1 as expansions around the photon-sphere radius. Using the standard null-geodesic potential W(r), they find the shadow radius shift r_sh=3√3 M(1-ε ν_1(3M)), Eq. (54). The EHT constraints on M87* and Sgr A* are then translated into bounds on the combination δ=-3√3 ε ν_1(3M), Eqs. (55)-(56), and the paper concludes that F(R)-gravity black-hole shadows may easily pass EHT constraints.
Significance. If the central claim were established, the paper would show that a wide class of F(R) gravity models can evade EHT shadow-radius constraints, and the perturbative framework for deriving the third-order F_R equation could be a useful technical tool for modified-gravity shadow calculations. The leading-order shadow formula is a straightforward consequence of the null-geodesic calculation once the perturbative ansatz is accepted. However, the main conclusion is not supported by the presented derivation: the quantity constrained by EHT is an arbitrary seed function ν_1, not a parameter of a specified F(R) action, and no global asymptotically flat solution is constructed. The paper is therefore better viewed as a consistency relation for a free metric perturbation than as a prediction of F(R) gravity.
major comments (4)
- [II B 1 and III D, Eqs. (29)-(35), (55)-(56)] The EHT constraints are imposed on ε ν_1(3M), where ν_1 is a freely chosen seed function introduced in Eq. (27). The paper explicitly states that 'we did not specify ν_1 as long as Cν_1 is finite and does not vanish.' Thus Eqs. (55) and (56) restrict the choice of ν_1, not any parameter of a concrete F(R) action. The conclusion 'by tuning the parameters, the models can satisfy the constraints' is therefore circular with respect to the claim that F(R) gravity predicts shadow sizes consistent with EHT; a genuine test would require fixing an action F(R), deriving ν_1 from the resulting field equations, and then comparing the predicted ε ν_1(3M) with the bounds.
- [II, Eq. (15)] The derivation of Eq. (15) is not shown. The text states that substituting Eq. (13) into Eq. (14) yields Eq. (15), but the algebra is nontrivial and no intermediate steps or independent check are provided. Since Eq. (29) and all subsequent shadow results depend on Eq. (15), the derivation should be supplied in an appendix or at least sketched in sufficient detail for the reader to verify the reduction.
- [II B 1, Eqs. (32)-(35)] The solution for f_R and λ_1 is only a local Taylor expansion around r=3M. No argument establishes that these local expressions extend to a global, horizon-regular, asymptotically flat solution. The shadow formula Eq. (41) is evaluated at the photon sphere but refers the radius to an asymptotic observer; without asymptotic flatness the interpretation of r_sh as the observed shadow size is not justified. The paper neither constructs the global metric nor invokes any existence or asymptotic-flatness result that would warrant the local-to-global step.
- [IV, Summary and Conclusion] The final sentence, 'Hence, it is proved that BH shadows in F(R) gravity may easily pass the Event Horizon Telescope constraints,' overstates the result. What is proved, given the assumptions, is that for any sufficiently regular seed ν_1 one can locally solve for f_R near r=3M and obtain a shadow shift proportional to ν_1(3M). Since ν_1 is not derived from a specified F(R) model and no global solution is shown, the statement about F(R) gravity as a theory is not established by the manuscript.
minor comments (4)
- [II, Eq. (13)] The displayed formula for N has ambiguous fraction formatting; the denominator should be written explicitly, e.g. N = exp(-∫ dr1 [F_R''/(F_R/r1 + (1/2)F_R')]), to avoid confusion.
- [III A, Eqs. (43)-(44)] The symbol E is used for both the conserved energy and the Lagrangian, and the sentence 'E = L vanishes identically E = L = 0' is confusing; a distinct symbol for the Lagrangian would make the null-geodesic argument clearer.
- [III C, after Eq. (54)] The statement that 'the expression in (41) is general' is too strong; Eq. (54) is derived under the specific local perturbative ansatz of Section II B 1 and the assumption Cν_1 ≠ 0.
- [II A, Eqs. (21)-(22)] The notation α_±^3 = -10/27 ± 90/27 and the subsequent numerical values appear correct but are displayed in a compressed way; writing the cube roots explicitly would improve readability.
Circularity Check
The claimed EHT compatibility is a fit to a free seed function ε ν1(3M), not a prediction of any specified F(R) gravity theory.
-
fitted input called prediction
[Section III C-D, Eqs. (50)-(56)]
"rsh = re−ν(r)| r=rph = 3√3M (1 − ǫν1 (r = 3M)) . ... In the model of Subsubsection II B 1, we did not specify ν1 as long as Cν1 ≡ (4M ν'''1 − 32/3 ν''1 + 8/(3M) ν'1) | r=3M is finite and does not vanish. Therefore, the expression in (41) is general."
The shadow radius in Eq. (54) is expressed entirely in terms of the free first-order metric seed ν1, which Assumption (27) introduces and Section II B 1 explicitly leaves unspecified. The observational constraints (55) and (56) are inequalities on the same combination −3√3 ε ν1(3M); they do not test F(R) gravity but merely pick out allowed values of the free input. A concrete F(R) action would determine ν1 through the field equations, so the EHT agreement is achieved by construction, not derived.
-
fitted input called prediction
[Section III D, Eqs. (55)-(56) and conclusion]
"−0.5 < −3√3ǫν1 (r = 3M) < 1.1 , (55) for M87∗ and −1.01 ≲ −3√3ǫν1 (r = 3M) ≲ 0.36 , (56) for Sgr A∗. Therefore by tuning the parameters, the models can satisfy the constraints."
The paper itself states that the constraints are satisfied by tuning parameters. But ε and ν1(3M) are not parameters of a fixed F(R) model; ν1 was declared arbitrary in Subsubsection II B 1. Thus the conclusion 'Hence, it is proved that BH shadows in F(R) gravity may easily pass the Event Horizon Telescope constraints' reduces to the statement that a free function can be chosen to match the data.
full rationale
The central claim is circular in the fitted-input sense: the constrained quantity is a free input of the perturbative construction. The paper never constructs a concrete F(R) action for which ν1 is forced; instead, Eq. (29) solves locally for f_R given an arbitrary ν1, and Eq. (54) then gives the shadow shift as −3√3 M ε ν1(3M). The EHT inequalities (55)-(56) bound this same free combination, so satisfying them is a choice of ν1, not a prediction of F(R) gravity. There is no self-citation chain or uniqueness theorem forcing the result, so this is not pattern 3 or 4; the circularity is of the fitted-input-called-prediction type. The separate absence of a global asymptotically flat extension is a correctness gap but does not affect the circularity scoring. Overall score 7: one or more 'predictions' reduce by construction to choosing a free seed function, and the paper explicitly says the parameters are tuned to satisfy the constraints.
Assumptions & free parameters
free parameters (5)
- epsilon (perturbative expansion parameter) =
unspecified small value
- delta = -3*sqrt(3)*epsilon*nu1(3M) =
M87*: between -0.5 and 1.1; Sgr A*: between -1.01 and 0.36
- f_R^0 and f_R^1 =
integration constants
- lambda_1^0 =
integration constant
- r0 =
integration constant
assumptions (6)
- standard math The F(R) field equations in vacuum, Eq. (2), are the correct equations of motion.
- domain assumption The metric can be written in the static spherically symmetric form (5).
- ad hoc to paper The deviation from Schwarzschild and from Einstein gravity is small, with F_R = 1 + (epsilon/2) f_R and nu = nu0 + epsilon nu1.
- ad hoc to paper A local solution of the approximate equation near r=3M extends to a global, physically acceptable solution.
- standard math The photon sphere is determined by W(r)=W'(r)=0 and the shadow radius by r_sh = r e^{-nu(r)} at r=r_ph.
- domain assumption The observational constraints quoted from the literature apply to the asymptotic shadow radius.
Cite this review
Pith. "Pith review of Black holes and their shadows in $F(R)$ gravity." pith.science (2026). https://pith.science/paper/GMUSRC55
@misc{pith2026241213775,
author = {Pith},
title = {Pith review of: Black holes and their shadows in $F(R)$ gravity},
year = {2026},
howpublished = {\url{https://pith.science/paper/GMUSRC55}},
note = {Machine review of arXiv:2412.13775}
}
abstract
We investigate the radii of the photon sphere and the black hole shadow in the framework of $F(R)$ gravity. For this purpose, we derive the field equation for the corresponding theory when the general spherically symmetric and static configuration is considered. This equation is the third-order differential equation with respect to $F_R(r)\equiv \left. \frac{dF(R)}{dR}\right|_{R=R(r)}$, where $r$ is the radial coordinate. Solving the equation, we find $F(R)$ as a function of $r$, $F_R=F_R(r)$. By using the assumed and obtained geometry, one can calculate the scalar curvature $R$ as a function of $r$, $R=R(r)$, which could be solved with respect to $r$ as $r=r(R)$. Then one finds the functional form of $F_R$ as a function of the scalar curvature $R$, $F_R=F_R(R)=F_R\left( r=r\left(R\right)\right)$.
Reference graph
Works this paper leans on
-
[1]
Theory Center, High Energy Accelerator Research Organizat ion (KEK), Oho 1-1, Tsukuba, Ibaraki 305-0801, Japan
-
[2]
Kobayashi-Maskawa Institute for the Origin of Particles an d the Universe, Nagoya University, Nagoya 464-8602, Japan
-
[3]
Can Magrans s/n, 08193 Barcelona, Spain
Institute of Space Sciences (ICE, CSIC) C. Can Magrans s/n, 08193 Barcelona, Spain
-
[4]
ICREA, Passeig Lluis Companys, 23, 08010 Barcelona, Spain 1 Abstract We investigate the radii of the photon sphere and the black ho le shadow in the framework of F (R) gravity. For this purpose, we derive the field equation for t he corresponding theory when the general spherically symmetric and static configuration is c onsidered. This equation is the thir...
-
[5]
′′ ( A 2r + 3B 4r2 + B′ 4r + BC 2r ) + ( − A 4r − A′ 16 − CA 8 − 5B 8r2 + B′ 8r + BC 4r ) FR 2 ( ( FR
-
[6]
′) 2 + B 4r FR 4( FR
-
[7]
′′′ + ( A 8 − B 2r + B′ 16 + BC 8 ) ( FR
-
[8]
′ + ( − A 16 + B 8r ) ( ( FR
Show all 43 references
-
[9]
′) 3 + B 16 FR 2( FR
-
[10]
′′′ − 3B 16 FR 2 ( ( FR
-
[11]
(15) Here A ≡ ν′′ + 2ν′2 − 1 r2 , B ≡ ν′ − 1 r , C ≡ ν′ + 1 r , (16) Eq
′′ . (15) Here A ≡ ν′′ + 2ν′2 − 1 r2 , B ≡ ν′ − 1 r , C ≡ ν′ + 1 r , (16) Eq. (15) is the third-order differential equation with respect to FR(r) if ν = ν(r) is given. Therefore for the given ν(r), we can find FR = FR(r) by solving Eq. (15). Then using (13), one can find λ = λ(r)...
-
[12]
′′ + 3 4r3 FR 2 ( ( FR
-
[13]
′) 2 − 1 4r2 FR 4( FR
-
[14]
′) 3 − 1 16r FR 2( FR
-
[15]
′′′ + 3 16r FR 2 ( ( FR
-
[16]
′′) 2 − 1 16r ( ( FR
-
[17]
(18) By assuming FR 2 ∝ rα with a constant α, we find 0 = α4 − 2α3 + 8α2 − 4α , (19) whose solution is given by, α = 0 , 2 3 + β0 , 2 3 + β±
′′ . (18) By assuming FR 2 ∝ rα with a constant α, we find 0 = α4 − 2α3 + 8α2 − 4α , (19) whose solution is given by, α = 0 , 2 3 + β0 , 2 3 + β± . (20) Here β0 ≡ α+ + α− , β ± ≡ α±ζ + α∓ζ 2 , (β+ = (β−)∗) , ζ ≡ ei 2π 3 = −1 2 + i √ 3 2 , α ± 3 ≡ − 10 27 ± 90 27 . (21) The nume...
-
[18]
(33) Here f 0 R and f 1 R are constants
Example 1 When r = 3M if 4M ν′′′ 1 − 32 3 ν′′ 1 + 8 3M ν′ 1 is finite and does not vanish, by defining a constant Cν1 as follows, Cν1 ≡ ( 4M ν′′′ 1 − 32 3 ν′′ 1 + 8 3M ν′ 1 ) ⏐ ⏐ ⏐ ⏐ r=3M , (31) a solution of f ′ R is given by f ′ R ∼ f 0 R + ( Cν1 − 2f 0 R 9M ) (r − 3M) , (32) ...
-
[19]
dot” or “˙
Example 2 Next we consider the case that r ∼ 3M if 4 M ν′′′ 1 − 32 3 ν′′ 1 + 8 3M ν′ 1 vanishes. Because the solution for λ of the equation 0 = 4 M λ2 − 32 3 λ + 8 3M is given by λ = λ± ≡ 8 3 ± 1 3 √ 58 . (36) ν1 behaves as ν1 = ν0 1 + C + ν1eλ +(r−3M ) + C − ν1eλ −(r−3M ) . (...
-
[20]
Akiyama et al
K. Akiyama et al. [Event Horizon Telescope], Astrophys. J. Lett. 875 (2019), L1 doi:10.3847/2041-8213/ab0ec7 [arXiv:1906.11238 [astro -ph.GA]]. 13
2019 arXiv
-
[21]
Khodadi, S
M. Khodadi, S. Vagnozzi and J. T. Firouzjaee, Sci. Rep. 14 (2024) no.1, 26932 doi:10.1038/s41598-024-78264-y [arXiv:2408.03241 [gr- qc]]
2024 arXiv
-
[22]
Nojiri and S
S. Nojiri and S. D. Odintsov, Phys. Dark Univ. 46 (2024), 101669 doi:10.1016/j.dark.2024.101669 [arXiv:2408.05668 [gr- qc]]
2024
-
[23]
Nojiri and S
S. Nojiri and S. D. Odintsov, Phys. Dark Univ. 46 (2024), 101702 doi:10.1016/j.dark.2024.101702 [arXiv:2410.11530 [gr- qc]]
2024
-
[24]
Nojiri and G
S. Nojiri and G. G. L. Nashed, Phys. Lett. B 830 (2022), 137140 doi:10.1016/j.physletb.2022.137140 [arXiv:2202.03693 [gr-qc]]
2022
-
[25]
A. Held, R. Gold and A. Eichhorn, JCAP 06 (2019), 029 doi:10.1088/1475-7516/2019/06/029 [arXiv:1904.07133 [gr-qc]]
2019 arXiv
-
[26]
Perlick and O
V. Perlick and O. Y. Tsupko, Phys. Rept. 947 (2022), 1-39 doi:10.1016/j.physrep.2021.10.004 [arXiv:2105.07101 [gr-qc]]
2022 arXiv
-
[27]
S. Chen, J. Jing, W. L. Qian and B. Wang, Sci. China Phys. Me ch. Astron. 66 (2023) no.6, 260401 doi:10.1007/s11433-022-2059-5 [arXiv:2301.0011 3 [astro-ph.HE]]
2023
-
[28]
K. i. Maeda, Phys. Rev. D 39 (1989), 3159 doi:10.1103/PhysRevD.39.3159
1989 doi
-
[29]
Capozziello, Int
S. Capozziello, Int. J. Mod. Phys. D 11 (2002), 483-492 doi:10.1142/S0218271802002025 [arXiv:gr-qc/0201033 [gr-qc]]
2002 arXiv
-
[30]
Nojiri and S
S. Nojiri and S. D. Odintsov, Phys. Rev. D 68 (2003), 123512 doi:10.1103/PhysRevD.68.123512 [arXiv:hep-th/0307288 [hep-th]]
2003 arXiv
-
[31]
Capozziello and M
S. Capozziello and M. De Laurentis, Phys. Rept. 509 (2011), 167-321 doi:10.1016/j.physrep.2011.09.003 [arXiv:1108.6266 [g r-qc]]
2011 arXiv
-
[32]
Nojiri and S
S. Nojiri and S. D. Odintsov, Phys. Rept. 505 (2011), 59-144 doi:10.1016/j.physrep.2011.04.001 [arXiv:1011.0544 [g r-qc]]
2011 arXiv
-
[33]
Nojiri, S
S. Nojiri, S. D. Odintsov and V. K. Oikonomou, Phys. Rept . 692 (2017), 1-104 doi:10.1016/j.physrep.2017.06.001 [arXiv:1705.11098 [ gr-qc]]
2017 arXiv
-
[34]
Beyond Einstein Gravit y: A Survey of Gravitational Theories for Cosmology and Astrophysics,
V. Faraoni and S. Capozziello, “Beyond Einstein Gravit y: A Survey of Gravitational Theories for Cosmology and Astrophysics,” Springer, 2011, ISBN 978- 94-007-0164-9, 978-94-007-0165-6 doi:10.1007/978-94-007-0165-6
2011 doi
-
[35]
Multamaki and I
T. Multamaki and I. Vilja, Phys. Rev. D 74 (2006), 064022 doi:10.1103/PhysRevD.74.064022 [arXiv:astro-ph/0606373 [astro-ph]]. 14
2006 arXiv
-
[36]
Capozziello, A
S. Capozziello, A. Stabile and A. Troisi, Class. Quant. Grav. 24 (2007), 2153-2166 doi:10.1088/0264-9381/24/8/013 [arXiv:gr-qc/0703067 [ gr-qc]]
2007 arXiv
-
[37]
de la Cruz-Dombriz, A
A. de la Cruz-Dombriz, A. Dobado and A. L. Maroto, Phys. R ev. D 80 (2009), 124011 [erratum: Phys. Rev. D 83 (2011), 029903] doi:10.1103/PhysRevD.80.124011 [arXiv: 0907.3872 [gr-qc]]
2009 arXiv
-
[38]
Capozziello, M
S. Capozziello, M. De laurentis and A. Stabile, Class. Q uant. Grav. 27 (2010), 165008 doi:10.1088/0264-9381/27/16/165008 [arXiv:0912.5286 [ gr-qc]]
2010 arXiv
-
[39]
Sebastiani and S
L. Sebastiani and S. Zerbini, Eur. Phys. J. C 71 (2011), 1591 doi:10.1140/epjc/s10052-011- 1591-8 [arXiv:1012.5230 [gr-qc]]
2011 arXiv
-
[40]
Capozziello, N
S. Capozziello, N. Frusciante and D. Vernieri, Gen. Rel . Grav. 44 (2012), 1881-1891 doi:10.1007/s10714-012-1367-y [arXiv:1204.4650 [gr-qc ]]
2012 arXiv
-
[41]
G. G. L. Nashed and S. Nojiri, Phys. Rev. D 102 (2020), 124022 doi:10.1103/PhysRevD.102.124022 [arXiv:2012.05711 [gr -qc]]
2020 arXiv
-
[42]
Bambi, K
C. Bambi, K. Freese, S. Vagnozzi and L. Visinelli, Phys. Rev. D 100 (2019) no.4, 044057 doi:10.1103/PhysRevD.100.044057 [arXiv:1904.12983 [gr -qc]]
2019 arXiv
-
[43]
Vagnozzi, R
S. Vagnozzi, R. Roy, Y. D. Tsai, L. Visinelli, M. Afrin, A . Allahyari, P. Bambhaniya, D. Dey, S. G. Ghosh and P. S. Joshi, et al. Class. Quant. Grav. 40 (2023) no.16, 165007 doi:10.1088/1361-6382/acd97b [arXiv:2205.07787 [gr-qc ]]. 15
2023 arXiv
Reviewed August 11, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.