REVIEW 3 major objections 5 minor 28 references
About rectified sigmoid function for enhancing the accuracy of Physics-Informed Neural Networks
T0 review · 3 major / 5 minor · reviewed 2026-08-10 · deepseek-v4-flash
Pith's one-line read A rectified sigmoid makes single-hidden-layer PINNs an order of magnitude more accurate on ODEs.
desk verdict A correct and neat derivation of a piecewise-linear activation trick for small PINNs, but the headline accuracy claim rests on interpolation with undisclosed reference data and is broader than the Lorenz numbers support. 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 rectified sigmoid activation, defined piecewise as $\mathrm{Re}\text{-}\sigma(x)=0$ for $x<-1$, $\tfrac12(x+1)$ for $-1\le x\le 1$, and $1$ for $x>1$, which is a sum of two ReLU functions: $\mathrm{Re}\text{-}\sigma(x)=\tfrac12[\mathrm{ReLU}(x+1)-\mathrm{ReLU}(x-1)]$. Because its derivative is the box function with $\mathrm{Re}\text{-}\sigma'(0)=1/2$ and with the chosen spacing $\Delta\zeta=1$ the derivative of one neuron does not overlap the neighbouring neuron, the initialization formula collapses to $W^{(2)}_k = -\Delta x\, N_l[u_k,x_k]/\Delta\zeta$, a direct closed-form expression for the output weights from the ODE residual at grid points. This converts the physics equations into an explicit, gradient-free construction of the network.
What would settle it
Take a new ODE with a known exact solution, initialize the Re-$\sigma$ network without supplying any grid values $\{u_k\}$ — setting $W^{(2)}_k$ only from the ODE residual and the initial condition — and compare against the sigmoid network. If the Re-$\sigma$ network no longer beats sigmoid by an order of magnitude, the reported gain depends on the precomputed solution points.
Extended reading notes
Core claim
The central claim is that substituting the sigmoid activation with the rectified sigmoid $\mathrm{Re}\text{-}\sigma(x) = 0$ for $x<-1$, $\tfrac12(x+1)$ for $-1\le x\le 1$, and $1$ for $x>1$ reduces the relative $L^2$ error of the predicted ODE solution by at least an order of magnitude for any given value, after the same initialization and training procedures. The paper reports this for PIDD initialization and for NbN training (three epochs), on the harmonic oscillator, the relativistic slingshot, and the Lorenz system. In the harmonic-oscillator case the Re-$\sigma$ network's error is close to the error of the standard high-precision reference solver itself.
Load-bearing premise
The initialization algorithm needs the solution values $\{u_k\}$ at the grid points before training, and the paper never states where they come from; if they come from the same reference solver used to measure error, the reported accuracy may reflect interpolation of that solver's output rather than physics-informed prediction.
Editorial extensions
If this is right
- Single-hidden-layer networks with the rectified sigmoid reach errors near $10^{-6}$ on the harmonic oscillator, close to the reference solver's own error in the reported setup.
- The reported order-of-magnitude gap over the sigmoid appears both after initialization alone (PIDD) and after the gradient-free NbN training, and across linear, nonlinear, and chaotic systems.
- Training stays cheap: the PIDD runs in milliseconds and the NbN fitting in a few minutes in the reported experiments, with three epochs sufficient.
- Because the same algorithms transfer from sigmoid to Re-$\sigma$ with only the spacing constant $\Delta\zeta$ changed, the method applies directly to the earlier single-hidden-layer toolbox.
Reading between the lines
- The paper compares Re-$\sigma$ only against the sigmoid; a direct comparison with plain ReLU under the same PIDD initialization would isolate whether the benefit comes from the rectified shape or from the smooth central slope, an experiment the paper does not run.
- The 'any given value' strength of the claim rests on three systems at two widths; a width and parameter sweep would show whether the gap persists for very small networks and for PDEs.
- Because the output weights are set by the residual formula, the scheme could be used for fast parameter sweeps or surrogate modeling of a family of ODEs by re-initializing with new residuals, an application left implicit.
- If the solution values consumed by the PIDD algorithm come from a reference solver, the reported errors measure interpolation of that solver's output; a natural test is to run the same initialization with the grid values withheld and watch whether the gap survives.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper proposes a rectified sigmoid activation function Re-sigma (Eq. 4) for single-hidden-layer physics-informed neural networks solving ODEs, together with a physics-informed data-driven (PIDD) initialization (Algorithm 1) and a neuron-by-neuron (NbN) training procedure (Algorithm 2) adapted from previous work. The main technical claim is that, because Re-sigma has compact support and Delta-zeta = 1, each neuron affects the derivative only at its own collocation point, so the output weights are determined exactly by the ODE residual, as shown in Eqs. (6)-(8). The numerical experiments in Section 3 compare relative L2 errors of networks with sigmoid and Re-sigma activations on a harmonic oscillator, a relativistic slingshot, and a Lorenz system, and Section 3.2 claims an order-of-magnitude error reduction for any given value. The central unresolved issue is the provenance of the solution values {u_k} consumed by Algorithm 1, which the paper never states.
Significance. If the central claim were fully established, the paper would offer a practically useful, almost parameter-free way to obtain accurate ODE solutions from tiny single-hidden-layer networks, with a simple and elegant derivative construction. The Section 2 derivation is correct as far as it goes: the compact-support argument makes Eq. (8) exact at the grid points, and the explicit algorithms make the method easy to reproduce. However, the claimed superiority over the sigmoid activation is only as strong as the source of the input data in Algorithm 1; if {u_k} are the reference solver values, the reported errors are interpolation errors, not physics-informed solution errors. The paper also does not quantify uncertainty, but because the algorithms are deterministic this is less important than the data-source issue.
major comments (3)
- [Algorithm 1; Section 3.2; Eq. (8)] The PIDD initialization in Algorithm 1 requires the solution values {u_k} at every collocation point, and Eq. (8) sets W_m^(2) = -Delta x N_l[u_m, x_m]/Delta zeta from those values. The paper never states where these values come from; in Section 3 the only numerical reference mentioned is scipy's odeint solver, used for the 'reference solutions'. If {u_k} are the odeint outputs, then the networks in Table 1 are constructed by interpolation of the reference trajectory, not by solving the ODE from the initial condition alone, and the comparison with the sigmoid activation is a comparison of two interpolation bases. Because Section 3.2's order-of-magnitude claim rests entirely on Table 1, this missing provenance is load-bearing. Please state the source of {u_k} explicitly, or re-run the PIDD experiments using only the initial condition and the ODE residual.
- [Section 3.2 vs. Section 3.3, Tables 1 and 2] The claim in Section 3.2 that replacing sigmoid with rectified sigmoid 'results in a reduction of the relative error by at least an order of magnitude for any given value' is not supported by the NbN results: in Table 2, for the Lorenz system the Re-sigma relative errors are only about 1.3 times smaller than the sigma errors (e.g., 2.63e-3 vs 3.47e-3 for x). If the claim is meant to cover only the PIDD initialization, the text should say so; if it is meant to cover the proposed method generally, the statement needs to be revised to match Table 2.
- [Section 3.2, harmonic oscillator paragraph] The observation that Re-sigma's relative L2 errors for the harmonic oscillator (3.88e-6 and 2.63e-6) are close to the odeint-vs-exact errors (1.03e-6 and 1.04e-6) is consistent with the interpretation that the Re-sigma network is reproducing the odeint reference solution to the accuracy of that solver. The paper should state whether the reported errors therefore measure the network's fit to the reference data, and should discuss how this relates to the claim that the network 'solves' the ODE.
minor comments (5)
- [Throughout] The system in Section 3.1.3 and in Tables 1 and 2 is called 'Lorentz system'; the standard name is 'Lorenz system' (also used in Eq. (13)).
- [Section 2, Eq. (6)] There is a typo: 'derivate' should be 'derivative', and the phrase 'Pay attention Re-sigma'(0) = 1/2' needs a comma or semicolon.
- [Algorithm 1 vs Eq. (8)] The subscript l is used in Eq. (8) (N_l[u_m, x_m]) but omitted in Algorithm 1's W_k^(2) update; make the component index consistent.
- [Algorithm 2] The initialization in Algorithm 2 uses N_l[u(0), x_k] but the notation u(0) is ambiguous for a system: it should be stated whether this means the constant initial vector evaluated at every collocation point, and how this differs from the correction term b0 used in Algorithm 1.
- [Section 3.3, Figure 1] The text says the absolute errors are 'significantly lower' than those of vanilla PINN, but no definition, settings, or results for vanilla PINN are given; either cite the specific comparison in [21] or add the missing details.
Circularity Check
PIDD initialization consumes the reference solution as input, so the Section 3.2 'order-of-magnitude' comparison measures in-sample quadrature accuracy, not independent PINN prediction.
-
fitted input called prediction
[Algorithm 1 (Data line, Eq. 8) and Section 3.2 (Table 1, reference-solution sentence)]
"Data: {u_k}_{k=0}^{N-1} for {x_k}_{k=0}^{N-1} ... W_k^(2) ← −Δx/Δζ N [u_k, x_k]; ... The reference solutions were obtained using the odeint solver of scipy.integrate library. ... Column 3 (relative L2 errors) shows the relative errors between solutions, provided by neural networks and solutions obtained using odeint solver of scipy.integrate library."
Algorithm 1 sets every second-layer weight from N[u_k, x_k] at every collocation point, so the solution values u_k must be known before the network is constructed. The only solution values named in Section 3 are the odeint reference values. With the Re-σ tent basis (W_k^(1)=2/Δx, b_k^(1)=−2k), the initialized network output at x_m is a trapezoidal quadrature of the derivative field built from those same u_k values. Table 1 then reports error against the odeint values that supplied Eq. 8. The network is thus constructed from the reference trajectory; the error measures quadrature accuracy of a data-driven fit, not independent PINN solving. The σ-vs-Re-σ comparison compares two interpolation bases, so the 'order of magnitude' sentence is an in-sample interpolation claim.
full rationale
The main circularity is concentrated in the PIDD experiments. The algorithm's input is the solution grid {u_k}; the experiments' only named source of such values is the odeint reference solution, and the accuracy metric in Table 1 is the deviation from that same reference. Because the Re-sigmoid basis has disjoint local support on the grid (each neuron activates only on its own cell and the neighboring cell), Eq. 8 turns the initialization into a quadrature formula for the reference derivative, so the 'prediction' is largely determined by the reference values themselves. This does not make the comparison between σ and Re-σ meaningless as a study of basis functions, but it means the strong claim in Section 3.2 ('reduction of the relative error by at least an order of magnitude for any given value') is not a demonstration that Re-sigmoid PINNs solve ODEs better from first principles. The NbN results (Table 2, Algorithm 2) are not circular in the same way: Algorithm 2 states 'Data: —' and updates weights from the network's own current solution, so those experiments provide partially independent support, although for the chaotic Lorenz case the Re-σ advantage is only about 1.3×, not an order of magnitude. No load-bearing self-citation or uniqueness-import pattern was found; Algorithms 3 and 4 are reproduced in the appendices, so the comparison rests on the paper's own equations rather than on [21] as an unverified authority. Overall, because the central quantitative claim rests on a data-driven initialization evaluated on the same data, the circularity score is 6.
Assumptions & free parameters
free parameters (3)
- Number of hidden neurons / collocation points N =
20000
- Number of training epochs E (NbN) =
3
- Number of time windows for NbN =
20
assumptions (4)
- domain assumption The ODE initial value problems are well-posed and have unique solutions on [0,T].
- domain assumption The reference solutions from scipy.integrate.odeint are accurate enough to serve as ground truth.
- ad hoc to paper The input data {u_k} for the PIDD initialization are available and consistent with the ODE.
- ad hoc to paper Three epochs of neuron-by-neuron updates suffice for convergence.
Cite this review
Pith. "Pith review of About rectified sigmoid function for enhancing the accuracy of Physics-Informed Neural Networks." pith.science (2026). https://pith.science/paper/LKKLNPAV
@misc{pith2026241220851,
author = {Pith},
title = {Pith review of: About rectified sigmoid function for enhancing the accuracy of Physics-Informed Neural Networks},
year = {2026},
howpublished = {\url{https://pith.science/paper/LKKLNPAV}},
note = {Machine review of arXiv:2412.20851}
}
read the original abstract
The article is devoted to the study of neural networks with one hidden layer and a modified activation function for solving physical problems. A rectified sigmoid activation function has been proposed to solve physical problems described by the ODE with neural networks. Algorithms for physics-informed data-driven initialization of a neural network and a neuron-by-neuron gradient-free fitting method have been presented for the neural network with this activation function. Numerical experiments demonstrate the superiority of neural networks with a rectified sigmoid function over neural networks with a sigmoid function in the accuracy of solving physical problems (harmonic oscillator, relativistic slingshot, and Lorentz system).
Figures
Reference graph
Works this paper leans on
-
[1]
Highly accurate protein structure prediction with AlphaFold,
J. Jumper, R. Evans, A. Pritzel, T. Green, M. Figurnov, O. Ronneberger, K. Tunyasuvunakool, R. Bates, A. Žídek, A. Potapenko, A. Bridgland, C. Meyer, S. A. A. Kohl, A. J. Ballard, A. Cowie, B. Romera-Paredes, S. Nikolov, R. Jain, J. Adler, T. Back, S. Petersen, D. Reiman, E. Clancy, M. Zielinski, M. Steinegger, M. Pacholska, T. Berghammer, S. Bodenstein, ...
work page 2021
-
[2]
A foundation model for the earth system,
C. Bodnar, W. P. Bruinsma, A. Lucic, M. Stanley, A. Vaughan, J. Brandstetter, P. Garvan, M. Riechert, J. A. Weyn, H. Dong, J. K. Gupta, K. Thambiratnam, A. T. Archibald, C.-C. Wu, E. Heider, M. Welling, R. E. Turner, and P. Perdikaris, “A foundation model for the earth system,” 2024. [Online]. Available: https://arxiv.org/abs/2405.13063
arXiv 2024
-
[3]
M. Raissi, P. Perdikaris, and G. E. Karniadakis, “Physics-informed neural networks: A deep learning framework for solving forward and inverse problems involving nonlinear partial differential equations,” Journal of Computational Physics , vol. 378, pp. 686–707, 2019. [Online]. Available: https://linkinghub.elsevier.com/retrieve/pii/S0021999118307125
work page 2019
-
[4]
Respecting causality is all you need for training physics-informed neural networks,
S. Wang, S. Sankaran, and P. Perdikaris, “Respecting causality is all you need for training physics-informed neural networks,” 2022. [Online]. Available: http://arxiv.org/abs/2203.07404
arXiv 2022
-
[5]
H. Wang, X. Qian, Y . Sun, and S. Song, “A Modified Physics Informed Neural Networks for Solving the Partial Differential Equation with Conservation Laws,” —. [Online]. Available: https://ssrn.com/abstract=4274376
-
[6]
Learning nonlinear operators via DeepONet based on the universal approximation theorem of operators,
L. Lu, P. Jin, G. Pang, Z. Zhang, and G. E. Karniadakis, “Learning nonlinear operators via DeepONet based on the universal approximation theorem of operators,” Nature Machine Intelligence, vol. 3, no. 3, pp. 218–229, mar
-
[7]
M. A. Krinitskiy, V . M. Stepanenko, A. O. Malkhanov, and M. E. Smorkalov, “A General Neural- Networks-Based Method for Identification of Partial Differential Equations, Implemented on a Novel AI Accelerator,” Supercomputing Frontiers and Innovations , vol. 9, no. 3, sep 2022. [Online]. Available: https://superfri.org/index.php/superfri/article/view/439
work page 2022
-
[8]
V . Fanaskov and I. Oseledets, “Spectral Neural Operators,” 2022. [Online]. Available: https: //arxiv.org/abs/2205.10573
arXiv 2022
Show all 28 references
-
[9]
Physics Informed Deep Learning (Part I): Data- driven Solutions of Nonlinear Partial Differential Equations,
M. Raissi, P. Perdikaris, and G. E. Karniadakis, “Physics Informed Deep Learning (Part I): Data- driven Solutions of Nonlinear Partial Differential Equations,” no. Part I, pp. 1–22. [Online]. Available: https://arxiv.org/abs/1711.10561 8 A preprint
-
[10]
About optimal loss function for training physics-informed neural networks under respecting causality,
V . A. Es’kin, D. V . Davydov, E. D. Egorova, A. O. Malkhanov, M. A. Akhukov, and M. E. Smorkalov, “About optimal loss function for training physics-informed neural networks under respecting causality,” 2023, arXiv:2304.02282
2023 arXiv
-
[11]
hp-VPINNs: Variational physics-informed neural networks with domain decomposition,
E. Kharazmi, Z. Zhang, and G. E. Karniadakis, “hp-VPINNs: Variational physics-informed neural networks with domain decomposition,” Computer Methods in Applied Mechanics and Engineering , vol. 374, p. 113547, 2021. [Online]. Available: https://doi.org/10.1016/j.cma.2020.113547
2021
-
[12]
B-PINNs: Bayesian physics-informed neural networks for forward and inverse PDE problems with noisy data,
L. Yang, X. Meng, and G. E. Karniadakis, “B-PINNs: Bayesian physics-informed neural networks for forward and inverse PDE problems with noisy data,” Journal of Computational Physics , vol. 425, p. 109913, 2021. [Online]. Available: https://doi.org/10.1016/j.jcp.2020.109913
2021
-
[13]
Scientific Machine Learning through Physics-Informed Neural Networks: Where we are and What’s next,
S. Cuomo, V . S. di Cola, F. Giampaolo, G. Rozza, M. Raissi, and F. Piccialli, “Scientific Machine Learning through Physics-Informed Neural Networks: Where we are and What’s next,” jan 2022. [Online]. Available: http://arxiv.org/abs/2201.05624
2022 arXiv
-
[14]
Thermodynamically consistent physics-informed neural networks for hyperbolic systems,
R. G. Patel, I. Manickam, N. A. Trask, M. A. Wood, M. Lee, I. Tomas, and E. C. Cyr, “Thermodynamically consistent physics-informed neural networks for hyperbolic systems,” Journal of Computational Physics , vol. 449, p. 110754, jan 2022. [Online]. Available: https://doi.org/10...
2022
-
[15]
Physics-informed neural networks (PINNs) for fluid mechanics: a review,
S. Cai, Z. Mao, Z. Wang, M. Yin, and G. E. Karniadakis, “Physics-informed neural networks (PINNs) for fluid mechanics: a review,” Acta Mechanica Sinica , vol. 37, no. 12, pp. 1727–1738, dec 2021. [Online]. Available: https://link.springer.com/10.1007/s10409-021-01148-1
2021 doi
-
[16]
Solving the wave equation with physics-informed deep learning,
B. Moseley, A. Markham, and T. Nissen-Meyer, “Solving the wave equation with physics-informed deep learning,” jun 2020. [Online]. Available: http://arxiv.org/abs/2006.11894
2020 arXiv
-
[17]
A-PINN: Auxiliary physics informed neural networks for forward and inverse problems of nonlinear integro-differential equations,
L. Yuan, Y .-Q. Ni, X.-Y . Deng, and S. Hao, “A-PINN: Auxiliary physics informed neural networks for forward and inverse problems of nonlinear integro-differential equations,” Journal of Computational Physics, vol. 462, p. 111260, aug 2022. [Online]. Available: https://doi.org...
2022
-
[18]
Ai-aristotle: A physics-informed framework for systems biology gray-box identification,
N. A. Daryakenari, M. D. Florio, K. Shukla, and G. E. Karniadakis, “Ai-aristotle: A physics-informed framework for systems biology gray-box identification,” 2023. [Online]. Available: https://arxiv.org/abs/2310.01433
2023 arXiv
-
[19]
Piratenets: Physics-informed deep learning with residual adaptive networks,
S. Wang, B. Li, Y . Chen, and P. Perdikaris, “Piratenets: Physics-informed deep learning with residual adaptive networks,” 2024. [Online]. Available: https://arxiv.org/abs/2402.00326
2024 arXiv
-
[20]
Element-wise multiplication based deeper physics-informed neural networks,
F. Jiang, X. Hou, and M. Xia, “Element-wise multiplication based deeper physics-informed neural networks,”
-
[21]
Are two hidden layers still enough for the physics-informed neural networks?
V . A. Es’kin, A. O. Malkhanov, and M. E. Smorkalov, “Are two hidden layers still enough for the physics-informed neural networks?” 2024. [Online]. Available: https://arxiv.org/abs/2412.19235
2024 arXiv
-
[22]
Three decades of activations: A comprehensive survey of 400 activation functions for neural networks,
V . Kunc and J. Kléma, “Three decades of activations: A comprehensive survey of 400 activation functions for neural networks,” 2024. [Online]. Available: https://arxiv.org/abs/2402.09092
2024 arXiv
-
[23]
Separable physics-informed neural networks for the solution of elasticity problems,
V . A. Es’kin, D. V . Davydov, J. V . Gur’eva, A. O. Malkhanov, and M. E. Smorkalov, “Separable physics-informed neural networks for the solution of elasticity problems,” 2024. [Online]. Available: https://arxiv.org/abs/2401.13486
2024 arXiv
-
[24]
Paszke, S
A. Paszke, S. Gross, F. Massa, A. Lerer, J. Bradbury, G. Chanan, T. Killeen, Z. Lin, N. Gimelshein, L. Antiga, A. Desmaison, A. Köpf, E. Yang, Z. DeVito, M. Raison, A. Tejani, S. Chilamkurthy, B. Steiner, L. Fang, J. Bai, and S. Chintala, PyTorch: An Imperative Style, High-Per...
2019
-
[25]
Relativistic slingshot: A source for single circularly polarized attosecond x-ray pulses,
J. Wang, S. V . Bulanov, M. Chen, B. Lei, Y . Zhang, R. Zagidullin, V . Zorina, W. Yu, Y . Leng, R. Li, M. Zepf, and S. G. Rykovanov, “Relativistic slingshot: A source for single circularly polarized attosecond x-ray pulses,”Phys. Rev. E, vol. 102, p. 061201, Dec 2020. [Online...
2020 doi
-
[26]
Deterministic nonperiodic flow,
E. N. Lorenz, “Deterministic nonperiodic flow,” Journal of the Atmospheric Sciences , vol. 20, pp. 130–141, 1963. [Online]. Available: https://api.semanticscholar.org/CorpusID:15359559 9
1963
-
[2021]
Available: https://doi.org/10.1038%2Fs42256-021-00302-5
[Online]. Available: https://doi.org/10.1038%2Fs42256-021-00302-5
-
[2024]
Available: https://arxiv.org/abs/2406.04170
[Online]. Available: https://arxiv.org/abs/2406.04170
Reviewed August 10, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.