REVIEW 4 major objections 5 minor 48 references
Cyclical accretion regime change in the slow X-ray pulsar 4U 0114+65 observed with Chandra
T0 review · 4 major / 5 minor · reviewed 2026-08-10 · deepseek-v4-flash
Pith's one-line read The pulse-to-pulse luminosity swings of the slow X-ray pulsar 4U 0114+65 are explained by a cyclical accretion regime change, in which matter accumulates at the magnetosphere until Compton cooling becomes efficient.
desk verdict New Chandra constraints on 4U 0114+65 are solid, but the cyclical accretion-regime change rests on an inconsistent luminosity threshold and only three pulse trains. 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 central mechanism is the Compton-cooling-regulated accumulation cycle in the quasi-spherical settling accretion theory. In this model a hot convective shell sits above the magnetosphere, and plasma enters through Rayleigh-Taylor instability with a settling velocity about half the free-fall velocity; the entry rate is controlled by Compton and radiative cooling. When the accumulated density $\rho_1$ above the magnetosphere rises, the Atwood number $A \approx (\rho_1-\rho_2)/(\rho_1+\rho_2)$ approaches unity (the paper assumes $A\sim 1$, with interior density $\rho_2 \ll \rho_1$), and the RTI growth rate $\sigma = (Agk)^{1/2}$ increases, so the short-spike timescale $t_A \sim (R_A^3/GM)^{1/2} \sim 150$ s shortens. The long-spike timescale $t_B \sim (R_B^3/GM)^{1/2} \sim 2\times10^3$ s is set by the free-fall time from the Bondi radius. The cycle thus alternates between a matter-depleted, fainter state and a matter-rich, Compton-cooling-efficient brighter state.
What would settle it
A long continuous observation covering several complete 2.6-hour pulses should show a quasi-periodic bright-to-faint cycle whose recurrence time matches the magnetospheric accumulation timescale, with the absorption column building before each bright phase and short-spike durations shortening as brightness rises; observing bright pulses without the preceding NH increase, or finding no periodicity in pulse amplitudes, would falsify the regime-change scenario.
Extended reading notes
Core claim
The central claim is that the pulse-to-pulse variability of 4U 0114+65, including the observed bright-to-faint sequences, is the visible signature of an accumulation-discharge cycle in quasi-spherical settling accretion onto a slowly rotating magnetized neutron star. After a bright episode depletes the plasma above the magnetosphere, the inflow rate is set by a settling velocity roughly two times lower than free fall, so matter gradually builds up. When the density exceeds a critical value, Compton cooling becomes efficient, the Rayleigh-Taylor instability at the magnetospheric boundary grows faster, the mass accretion rate increases, and the source enters a brighter regime. The same cycle explains why the short spikes (timescale ~150 s) are shorter in brighter states: higher accumulated density lowers the magnetospheric radius and shortens the RTI timescale $t_A \sim R_A^{3/2}/(GM)^{1/2}$. Long spikes (~2e3 s) are attributed to density and velocity fluctuations in the captured stellar wind near the Bondi radius. Spectral support comes from the absorption column being about twice as high in the bright spectra than in the faint ones, and from the Comptonization index $\alpha$ changing from ~0.1 (bright) to ~4 (faint).
Load-bearing premise
The whole cycle rests on the assumption that this neutron star is in the quasi-spherical settling accretion regime, which requires its X-ray luminosity to sit below a Compton-cooling threshold; the paper quotes that threshold at $\simeq 4\times10^{36}$ erg/s in one section but ~$4\times10^{34}$ erg/s in another, while the source is measured at $0.1$–$1.3\times10^{36}$ erg/s.
Editorial extensions
If this is right
- If the cycle is real, the bright-to-faint sequences in the 4U 0114+65 light curve are deterministic outcomes of the accumulation–depletion cycle, not stochastic flaring, and the relative durations of bright and faint phases should reflect the ratio of the settling and free-fall timescales.
- The negative correlation between short-spike duration and source brightness follows directly from $t_A \propto R_A^{3/2}$ with $R_A$ decreasing at higher luminosity, so measuring spike durations in any additional luminosity state tests the RTI interpretation.
- The factor about two higher absorption column in the Bright spectra is a direct signature of matter accumulated in the shell before the bright phase; continued spectral monitoring across a full cycle would show the column rising before each bright pulse.
- Because 4U 0114+65 is near the boundary between settling and free-fall accretion regimes, the same cyclical behaviour should appear in other slow, wind-accreting X-ray pulsars whose luminosities approach the Compton-cooling threshold.
Reading between the lines
- One testable consequence the paper leaves implicit: the recurrence time of the bright phases should be quasi-periodic, equal to the time needed to re-fill the magnetospheric shell; a long continuous observation of this source should show a periodicity in pulse amplitudes at that accumulation timescale.
- The paper's conflicting threshold values ($\simeq 4\times10^{36}$ erg/s in the Introduction versus ~$4\times10^{34}$ erg/s in the Discussion) leave the source's regime membership uncertain at its measured $0.1$–$1.3\times10^{36}$ erg/s; a systematic study of how spike properties change as the source's average luminosity crosses these values would determine the true regime boundary.
- If the same Compton-cooling threshold controls flaring in supergiant fast X-ray transients, their short spikes should also anti-correlate with brightness; checking existing SFXT light curves for this pattern would tell whether a common mechanism is at work.
- The phase-resolved pulse profiles could be used to see a spectral drift within a single bright pulse from an efficient-Comptonization (small $\alpha$) start to a less efficient (large $\alpha$) end, which would confirm that the bright phase depletes the accumulated matter on timescales of one to two hours.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper analyzes 200 ks of Chandra HETG/ACIS observations of the slow X-ray pulsar 4U 0114+65, spread over nine uninterrupted exposures across four orbits. The authors measure a spin period of 9050±100 s, segment the light curves into peaks and valleys, classify the peaks into Bright, Intermediate, Faint, and Absorbed groups using k-means clustering, and detect short and long spikes using an autoencoder-based anomaly detector. From the orbital variation of the absorption column they infer an inclination of ~40° and a lower limit on the companion mass-loss rate of ~8.6e-7 solar masses per year. The central interpretive claim is that the pulse-to-pulse luminosity variations, in particular three decreasing bright-to-faint pulse trains, reflect a cyclical accretion regime change within the quasi-spherical settling accretion framework: matter accumulates at the magnetosphere until a critical density triggers more efficient Compton cooling, increasing the accretion rate and brightness until the reservoir is depleted.
Significance. If the regime-cycle interpretation holds, the paper would provide an observational diagnostic of the quasi-spherical settling accretion theory in a slow X-ray pulsar, connecting three independent datasets: pulse-to-pulse flux evolution, short-spike durations, and spectral parameters (NH and the bmc index alpha). The analysis is based on standard Chandra data reduction and spectral fitting, and the orbital NH modeling with a PSO fit is a reasonable approach that yields a plausible inclination and a lower limit on the wind mass-loss rate. The classification of long and short spikes into distinct clusters is a useful empirical result, and the short-spike timescale comparison with Eq. (6) is presented as an explicit, falsifiable check rather than a fitted parameter. However, the central regime-change claim currently rests on a small number of pulse trains and on a theoretical regime boundary that is stated inconsistently in the manuscript, so the significance of the paper as a test of settling-accretion theory remains conditional.
major comments (4)
- [§1 and §5, 'NS pulse to pulse variability'] The manuscript gives two mutually inconsistent luminosity boundaries for the quasi-spherical settling accretion regime: the Introduction states that this mode is present in sources with moderate X-ray luminosity 'below ≃ 4×10^36 erg s−1', while Section 5 states that 'in sources with moderate X-ray luminosity (below ∼ 4×10^34 erg s−1), a hot convective quasi-spherical shell forms above the NS's magnetosphere.' With the source luminosity in Table 1 of 0.1–1.3×10^36 erg s−1, the source lies inside the regime under the first boundary but 2.5–33 times above the second boundary. Because the entire cyclical Compton-cooling explanation presupposes that the source is in the settling regime with a hot shell above the magnetosphere, the lower threshold would invalidate the central mechanism. The authors must resolve this contradiction: if 4×10^34 is a typo, it must be corrected and the sentence reworded; if it is not, the proposed cycle cannot operate at this luminosity and a different theoretical framework is needed. This point is load-bearing, not stylistic.
- [§5, 'NS pulse to pulse variability'] The central claim of a cyclical regime change is not backed by a quantitative model. The text asserts that matter accumulates at the magnetospheric boundary until 'the density exceeds a critical threshold' and that this triggers more efficient Compton cooling, but no expression is given for this threshold, for the accumulation timescale, for the depletion phase, or for the resulting cycle period. The observational evidence consists of three decreasing pulse trains (Orb 1-1, Orb 14-3, Orb 31-1; Fig. 14), without a statistical test of whether such sequences are expected or unlikely under a null model of stochastic wind accretion. As written, the cycle scenario is a plausible narrative rather than a tested explanation, and competing mechanisms (e.g., wind clumps, temporary gating by the magnetosphere, or variable absorption) are not quantitatively excluded. A quantitative derivation, even order-of-magnitude, of the cycle threshold and timescale is required to support the regime-change claim.
- [§3.3 and §5, Eqs. (5)–(6)] The short-spike timescale estimate tA ∼ R_A^{3/2}/(GM)^{1/2} ∼ 150 s depends on two assumptions that are stated rather than derived: the Atwood number A is taken to be close to 1 without justification, and the magnetospheric radius R_A is taken from Table 1 without an uncertainty or a consistency check against the luminosity-dependent accretion rate. Since R_A depends on the mass accretion rate and the magnetic field of the neutron star, and since the paper itself notes that R_A is inversely proportional to luminosity, the agreement between tA and the observed spike durations should be verified across the full range of luminosities in Table 1, not just at a single representative value. This is not necessarily fatal—the comparison is an independent check—but it needs to be presented with the relevant ranges and caveats.
- [§5, spectral support, Table 6] The spectral support for the accumulation scenario is suggestive but not decisive. The paper notes that NH is about two times larger in the Bright spectrum than in the Faint spectrum (2.1±0.1 versus 1.1±0.2 ×10^22 cm−2) and that the bmc index alpha differs (∼0.01 versus ∼4), but these are four-point correlations within a single source, and the Absorbed spectrum has NH = 19×10^22 cm−2, which is dominated by the orbital phase around ϕ ∼ 0 rather than by the accumulation cycle. The claim that NH traces local matter accumulation above the magnetosphere would be strengthened by showing that the NH difference between Bright and Faint is not simply a phase effect, and by quantifying the expected column density contributed by the accumulating shell.
minor comments (5)
- [§3.1] The phrase 'To characterize the spin period is not straightforward thought' appears to contain a typo; 'thought' should be 'though'.
- [§1 and Table 1] The companion spectral type is printed as 'BI1a' in the Introduction but as 'B1Ia' in Table 1; the standard notation should be used consistently.
- [Appendix references] Several cross-references to tables and figures are ambiguous or inconsistent, for example 'Table Fig. C.1.', 'Table Fig. D.1.', and 'Table Fig. E.1.'; these should be cleaned up.
- [§2] The observation log lists nine exposures, but the prose says 'nine different uninterrupted light curves spread across four different orbits'; this is consistent, but the table would benefit from a column identifying the orbit number more explicitly to avoid confusion between observation ID and orbit number.
- [§3.2] The pulse fraction is quoted with very asymmetric or large uncertainties (e.g., 0.6 ± 0.6 for Orb 1-1), and the text states a correlation of r = 0.9 between pulse fraction and orbital evolution without giving a significance level or showing the correlation on the relevant figure; a quantitative statement of the correlation's uncertainty would be helpful.
Circularity Check
No significant circularity: the spike-timescale checks use independent radii, the accretion-regime framework is external theory, and no load-bearing step reduces by construction to its own input.
full rationale
The short- and long-spike timescale checks are genuine, non-circular comparisons. t_A is computed from the tabulated Alfvén radius R_A ≃ 1.7×10^10 cm and t_B from the Bondi radius R_B ≃ 10^11 cm using standard scalings (Eqs. 4–6); neither radius was fitted to the observed spike durations, and the predicted t_A ≃ 150 s is compared with independently measured spike durations. The anti-correlation between short-spike duration and brightness follows from the scaling R_A ∝ L^{-2/7} and is tested against the data, making it a prediction rather than a restatement of inputs. The quasi-spherical settling accretion framework is cited from Shakura et al. (2012, 2014, 2015) and Shakura & Postnov (2017); although one author is a coauthor, that framework is an external published theory with independent applications and is not invoked as a uniqueness theorem. The RTI mechanism is also anchored to Arons & Lea (1976) and Elsner & Lamb (1977). The spectral parameters (N_H, alpha) used to support the scenario are fitted from the same Chandra data, so that support is post-hoc rather than a predictive test, but the paper does not claim these parameters are derived from the cycle model, so this is not circularity by construction. The luminosity-regime inconsistency (4×10^36 vs 4×10^34 erg/s) is a correctness/consistency concern about whether the settling regime applies to 4U 0114+65, not a circularity in the derivation chain. No load-bearing step reduces to its own input.
Assumptions & free parameters
free parameters (6)
- Orbital inclination =
40 ± 3 degrees
- Companion mass-loss rate =
8.6 ± 1.7 x 10^-7 M_sun/yr (lower limit)
- BMC spectral parameters for Bright/Intermediate/Faint/Absorbed spectra =
NH ranges from about 1.1 to 23 x 10^22 cm^-2; alpha ranges from about 0.1 to 4
- Soft blackbody temperature and absorption =
kT = 0.06 keV, NH = 0.8 x 10^22 cm^-2
- Atwood number A =
~1 (assumed)
- Magnetospheric (Alfven) radius RA =
1.7 x 10^10 cm
assumptions (6)
- domain assumption CAK stellar wind model with beta = 0.8 and the wind is unionized and spherically symmetric.
- domain assumption Quasi-spherical settling accretion theory (Shakura et al. 2012, 2014) with settling velocity about 2 times smaller than free-fall and a Compton-cooling threshold.
- domain assumption Neutron star mass of 1.4 M_sun.
- standard math Rayleigh-Taylor instability dispersion relation sigma = sqrt(A g k) with A ~ 1.
- domain assumption Blackbody seed-photon radius formula R_W = 0.6 sqrt(L_34)(kT)^-2 km.
- ad hoc to paper The pulse-phase shift in Orb 32-1 is caused by pole switching.
Cite this review
Pith. "Pith review of Cyclical accretion regime change in the slow X-ray pulsar 4U 0114+65 observed with Chandra." pith.science (2026). https://pith.science/paper/ZBM5RE6I
@misc{pith2026250108702,
author = {Pith},
title = {Pith review of: Cyclical accretion regime change in the slow X-ray pulsar 4U 0114+65 observed with Chandra},
year = {2026},
howpublished = {\url{https://pith.science/paper/ZBM5RE6I}},
note = {Machine review of arXiv:2501.08702}
}
abstract
4U 0114+65 is a high-mass X-ray binary system formed by the luminous supergiant B1Ia, known as V{*} V662 Cas, and one of the slowest rotating neutron stars (NS) with a spin period of about 2.6 hours. This fact provides a rare opportunity to study interesting details of the accretion within each individual pulse of the compact object. In this paper, we analyze 200 ks of Chandra grating data, divided into 9 uninterrupted observations around the orbit. The changes in the circumstellar absorption column through the orbit suggest an orbital inclination of $\sim$ $40^{\circ}$ with respect to the observer and a companion mass-loss rate of $\sim$ 8.6 10$^{-7}$ solar masses yr$^{-1}$. The peaks of the NS pulse show a large pulse-to-pulse variability. Three of them show an evolution from a brighter regime to a weaker one. We propose that the efficiency of Compton cooling in this source fluctuates throughout an accumulation cycle. After significant depletion of matter within the magnetosphere, since the settling velocity is $\sim \times$ 2 times lower than the free-fall velocity, the source gradually accumulates matter until the density exceeds a critical threshold. This increase in density triggers a transition to a more efficient Compton cooling regime, leading to a higher mass accretion rate and consequently to an increased brightness.
Figures
Figures from the paper (12 more)
Reference graph
Works this paper leans on
-
[1]
Abdallah, M. H., Samir, R. M., Leahy, D. A., & Shaker, A. A. 2023, MNRAS, 522, 3271
work page 2023
- [2]
-
[3]
Bailer-Jones, C. A. L., Rybizki, J., Fouesneau, M., Demleitner, M., & Andrae, R. 2021, The Astronomical Journal, 161, 147
work page 2021
-
[4]
2020, arXiv e-prints, arXiv:2003.05991
Bank, D., Koenigstein, N., & Giryes, R. 2020, arXiv e-prints, arXiv:2003.05991
arXiv 2020
-
[5]
Bonning, E. W. & Falanga, M. 2005, A&A, 436, L31
work page 2005
-
[6]
Bonyadi, M. R. & Michalewicz, Z. 2017, Evolutionary Computation, 25, 1
work page 2017
- [7]
-
[8]
Bulik, T., Gondek-Rosi´nska, D., Santangelo, A., et al. 2003, A&A, 404, 1023
work page 2003
Show all 48 references
-
[9]
I., Abbott, D
Castor, J. I., Abbott, D. C., & Klein, R. I. 1975, ApJ, 195, 157
1975
-
[10]
B., & Cowley, A
Crampton, D., Hutchings, J. B., & Cowley, A. P. 1985, ApJ, 299, 839 D’Aì, A., Evans, P. A., Burrows, D. N., et al. 2016, ArXiv e-prints [arXiv:1607.04264]
1985 arXiv
-
[11]
1977, IAU Circ., 3144
Dower, R., Kelley, R., Margon, B., & Bradt, H. 1977, IAU Circ., 3144
1977
-
[12]
Duncan, R. C. & Thompson, C. 1992, ApJ, 392, L9
1992
-
[13]
Elsner, R. F. & Lamb, F. K. 1977, ApJ, 215, 897
1977
-
[14]
A., Sood, R
Farrell, S. A., Sood, R. K., & O’Neill, P. M. 2006, MNRAS, 367, 1457
2006
-
[15]
P., Belloni, T., & Cassinelli, J
Finley, J. P., Belloni, T., & Cassinelli, J. P. 1992, A&A, 262, L25
1992
-
[16]
C., Allen, G
Fruscione, A., McDowell, J. C., Allen, G. E., et al. 2006, in Society of Photo- Optical Instrumentation Engineers (SPIE) Conference Series, V ol. 6270, Ob- servatory Operations: Strategies, Processes, and Systems, ed. D. R. Silva & R. E. Doxsey, 62701V
2006
-
[17]
D., Blair, J
Grundstrom, E. D., Blair, J. L., Gies, D. R., et al. 2007, ApJ, 656, 431
2007
-
[18]
A., Finley, J
Hall, T. A., Finley, J. P., Corbet, R. H. D., & Thomas, R. C. 2000, ApJ, 536, 450
2000
-
[19]
J., Kahn, S
Hess, C. J., Kahn, S. M., & Paerels, F. B. 1993, in American Astronomical So- ciety Meeting Abstracts, V ol. 183, American Astronomical Society Meeting Abstracts, 65.02
1993
-
[20]
C., Narayan, R., & Kallman, T
Hickox, R. C., Narayan, R., & Kallman, T. R. 2004, ApJ, 614, 881
2004
-
[21]
Y ., Lin, L
Hu, C.-P., Chou, Y ., Ng, C. Y ., Lin, L. C.-C., & Yen, D. C.-C. 2017, ApJ, 844, 16
2017
-
[22]
C.-C., & Yen, D
Hu, C.-P., Chou, Y ., Ng, C.-Y ., Lin, L. C.-C., & Yen, D. C.-C. 2017, The Astro- physical Journal, 844, 16
2017
-
[23]
Ikhsanov, N. R. & Beskrovnaya, N. G. 2013, Astronomy Reports, 57, 287
2013
-
[24]
2006, A&A, 458, 513 Krtiˇcka, J., Kubát, J., & Krtiˇcková, I
Koenigsberger, G., Georgiev, L., Moreno, E., et al. 2006, A&A, 458, 513 Krtiˇcka, J., Kubát, J., & Krtiˇcková, I. 2015, A&A, 579, A111
2006
-
[25]
Lamers, H. J. G. L. M. & Cassinelli, J. P. 1999, Introduction to Stellar Winds, 452
1999
-
[26]
& van den Heuvel, E
Li, X.-D. & van den Heuvel, E. P. J. 1999, ApJ, 513, L45
1999
-
[27]
& Blomme, R
Lobel, A. & Blomme, R. 2008, ApJ, 678, 408
2008
-
[28]
1967, in Some methods for classification and analysis of multi- variate observations Martínez-Núñez, S., Kretschmar, P., Bozzo, E., et al
MacQueen, J. 1967, in Some methods for classification and analysis of multi- variate observations Martínez-Núñez, S., Kretschmar, P., Bozzo, E., et al. 2017, Space Sci. Rev.[arXiv:1701.08618]
1967 arXiv
-
[29]
W., et al
Massa, D., Oskinova, L., Fullerton, A. W., et al. 2014, MNRAS, 441, 2173
2014
-
[30]
M., Clarke, D., & Pollock, A
Oskinova, L. M., Clarke, D., & Pollock, A. M. T. 2001, A&A, 378, L21
2001
-
[31]
C., Bozzo, E., & Belloni, T
Pradhan, P., Paul, B., Paul, B. C., Bozzo, E., & Belloni, T. M. 2015, MNRAS, 454, 4467
2015
-
[32]
2016, ArXiv e-prints [arXiv:1607.04107]
Rea, N., Borghese, A., Esposito, P., et al. 2016, ArXiv e-prints [arXiv:1607.04107]
2016 arXiv
-
[33]
J., et al
Reig, P., Chakrabarty, D., Coe, M. J., et al. 1996, A&A, 311, 879
1996
-
[34]
Reig, P., Nersesian, A., Zezas, A., Gkouvelis, L., & Coe, M. J. 2016, A&A, 590, A122
2016
-
[35]
2009, Gaussian Mixture Models, ed
Reynolds, D. 2009, Gaussian Mixture Models, ed. S. Z. Li & A. Jain (Boston, MA: Springer US), 659–663
2009
-
[36]
Rousseeuw, P. J. 1987, Journal of Computational and Applied Mathematics, 20, 53 Sanjurjo-Ferrín, G., Torrejón, J. M., Postnov, K., et al. 2017, A&A, 606, A145 Sanjurjo-Ferrín, G., Villalva, J. P., Torrejón, J. M., & Rodes-Roca, J. J. 2024, Journal of Open Source Software, 9, 7220
1987
-
[37]
Scargle, J. D. 1982, ApJ, 263, 835
1982
- [38]
-
[39]
2012, MNRAS, 420, 216
Shakura, N., Postnov, K., Kochetkova, A., & Hjalmarsdotter, L. 2012, MNRAS, 420, 216
2012
-
[40]
I., Postnov, K
Shakura, N. I., Postnov, K. A., Kochetkova, A. Y ., & Hjalmarsdotter, L. 2014, in European Physical Journal Web of Conferences, V ol. 64, European Physical Journal Web of Conferences, 02001
2014
-
[41]
I., Postnov, K
Shakura, N. I., Postnov, K. A., Kochetkova, A. Y ., et al. 2015, Astronomy Re- ports, 59, 645
2015
-
[42]
A., Belfiore, A., et al
Sidoli, L., Postnov, K. A., Belfiore, A., et al. 2019, MNRAS, 487, 420
2019
-
[43]
Sood, R., Farrell, S., O’Neill, P., Manchanda, R., & Ashok, N. M. 2006, Ad- vances in Space Research, 38, 2779
2006
-
[44]
Titarchuk, L., Mastichiadis, A., & Kylafis, N. D. 1997, ApJ, 487, 834 Torrejón, J. M., Kreykenbohm, I., Orr, A., Titarchuk, L., & Negueruela, I. 2004, A&A, 423, 301
1997
-
[45]
S., de Koter, A., & Lamers, H
Vink, J. S., de Koter, A., & Lamers, H. J. G. L. M. 1999, A&A, 350, 181
1999
-
[46]
2011, MNRAS, 413, 1083
Wang, W. 2011, MNRAS, 413, 1083
2011
-
[47]
C., Tananbaum, H
Weisskopf, M. C., Tananbaum, H. D., Van Speybroeck, L. P., & O’Dell, S. L. 2000, in X-Ray Optics, Instruments, and Missions III, ed. J. E. Truemper & B. Aschenbach (SPIE)
2000
-
[48]
2000, ApJ, 542, 914 Article number, page 11 of 11
Wilms, J., Allen, A., & McCray, R. 2000, ApJ, 542, 914 Article number, page 11 of 11
2000
Reviewed August 10, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.